Changeset 3274 for trunk/MagicSoft/Mars/macros

Timestamp:

02/24/04 15:49:17 (21 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/ONOFFAnalysis.C (modified) (16 diffs)

trunk/MagicSoft/Mars/macros/ONOFFAnalysis.C

@@ -172 +172 @@
       //-----------------------------------------------
       //const char *offfile = "~magican/ct1test/wittek/offdata.preproc"; 
-      const char *offfile = "20040126_OffCrab_"; 
+      //const char *offfile = "20040126_OffCrab_"; 
+      //const char *offfile  = "*.OFF";
+      const char *offfile  = "12*.OFF";
 
       //const char *onfile  = "~magican/ct1test/wittek/mkn421_on.preproc"; 
       //      const char *onfile  = "~magican/ct1test/wittek/mkn421_00-01"; 
-      const char *onfile  = "20040126_Crab_";
+      //const char *onfile  = "20040126_Crab_";
+      //const char *onfile  = "MCerPhot_output";
+      //const char *onfile  = "*.ON";
+      const char *onfile  = "12*.ON";
 
      const char *mcfile  = "/data/MAGIC/mc_eth/magLQE_3/gh/0/0/G_M0_00_0_550*.root";
@@ -184 +189 @@
 
       // path for input for Mars
-      TString inPath = "/.magic/magicserv01/scratch/";
+      //TString inPath = "/.magic/magicserv01/scratch/";
+     TString inPath = "/mnt/data17a/hbartko/";
+     //TString inPath = "~wittek/datacrab_feb04/";
 
       // path for output from Mars
-      TString outPath = "~wittek/datacrab/";
+      TString outPath = "~wittek/datacrab_feb04/";
 
       //-----------------------------------------------
@@ -235 +242 @@
     Bool_t JobB_SC_UP  = kFALSE;
     Bool_t CMatrix     = kFALSE;  // create training and test matrices 
-    Bool_t RMatrix     = kTRUE;  // read training and test matrices from file
-    Bool_t WOptimize   = kTRUE;  // do optimization using the training sample
+    Bool_t RMatrix     = kFALSE;  // read training and test matrices from file
+    Bool_t WOptimize   = kFALSE;  // do optimization using the training sample
                                   // and write supercuts parameter values 
                                   // onto the file parSCfile
     Bool_t RTest       = kFALSE;  // test the supercuts using the test matrix
-    Bool_t WSC         = kFALSE;  // update input root file ?
+    Bool_t WSC         = kTRUE;  // update input root file ?
 
 
@@ -316 +323 @@
     TString fileON  = inPath;
     fileON += onfile;
-    fileON += "CalibratedEvts";
     fileON += ".root";
 
     TString fileOFF = inPath;
     fileOFF += offfile;
-    fileOFF += "CalibratedEvts";
     fileOFF += ".root";
 
     TString fileMC  = mcfile;
     fileMC += mcfile;
-    fileMC += "CalibratedEvts";
     fileMC += ".root";
 
@@ -339 +343 @@
     //--------------------------------------------------
     // type of data to be padded 
-    //TString typeInput = "ON";
-    TString typeInput = "OFF";
+    TString typeInput = "ON";
+    //TString typeInput = "OFF";
     //TString typeInput = "MC";
     gLog << "typeInput = " << typeInput << endl;
@@ -366 +370 @@
     outNameImage += "Hillas";
     outNameImage += typeInput;
-    outNameImage += "1.root";
+    outNameImage += "1a.root";
     gLog << "padded data to be written onto : " << outNameImage << endl;
 
@@ -680 +684 @@
     MGeomApply        apply;
 
-    MPedestalWorkaround waround;
+    //MPedestalWorkaround waround;
 
     // a way to find out whether one is dealing with MC :
@@ -693 +697 @@
     //                                             "MCT1PointingCorrCalc");
     //}
+    MSourcePosfromStarPos sourcefromstar;
+    sourcefromstar.SetSourceAndStarPosition("Crab", 22, 0, 52, 5, 34, 32,
+                                        "Zeta-Tau", 21, 8, 33, 5, 37, 38.7);
+    sourcefromstar.AddFile("~wittek/datacrab_feb04/positions.4.txt", 0);
 
     MBlindPixelCalc blindbeforepad;
@@ -706 +714 @@
     contbasic.SetName("SelBasic");
 
-    //MFillH fillblind("BlindPixels[MHBlindPixels]", "MBlindPixels");
-    //fillblind.SetName("HBlind");
+    MFillH fillblind("BlindPixels[MHBlindPixels]", "MBlindPixels");
+    fillblind.SetName("HBlind");
 
     MSigmabarCalc sigbarcalc;
@@ -750 +758 @@
     selstandard.SetHillasName(fHilName);
     selstandard.SetImgParName(fImgParName);
-    selstandard.SetCuts(100, 4, 20, 0.0, 1.0, 0.0, 0.0);
+    selstandard.SetCuts(200, 6, 600, 0.4, 1.1, 0.0, 0.0);
     MContinue contstandard(&selstandard);
     contstandard.SetName("SelStandard");
@@ -786 +794 @@
     //tliston.AddToList(&f2);
     tliston.AddToList(&apply);
-    tliston.AddToList(&waround);
+    //tliston.AddToList(&waround);
+    tliston.AddToList(&sourcefromstar);
 
     //tliston.AddToList(&blindbeforepad);
@@ -793 +802 @@
     //  tliston.AddToList(&pointcorr);
 
-    //tliston.AddToList(&blind);
-    tliston.AddToList(&contbasic);
-
-    //tliston.AddToList(&fillblind);
+    tliston.AddToList(&blind);
+    //tliston.AddToList(&contbasic);
+
+    tliston.AddToList(&fillblind);
     tliston.AddToList(&sigbarcalc);
     tliston.AddToList(&fillsigtheta);
@@ -1783 +1792 @@
     TString parSCinit = outPath;
     //parSCinit += "parSC_060204";
+    //parSCinit = "parSC_240204a";
     parSCinit = "";
 
@@ -1802 +1812 @@
 
     TString parSCfile = outPath;
-    parSCfile += "parSC_130204a";
+    parSCfile += "parSC_240204a";
 
     gLog << "parSCfile = " << parSCfile << endl;
@@ -1843 +1853 @@
     filenameTrain += typeInput;
     filenameTrain += "1.root";
-    Int_t howManyTrain = 7000;
+    Int_t howManyTrain = 20000;
     gLog << "filenameTrain = " << filenameTrain << ",   howManyTrain = "
          << howManyTrain  << endl; 
@@ -1853 +1863 @@
     filenameTest += typeInput;
     filenameTest += "1.root";
-    Int_t howManyTest = 7000;
+    Int_t howManyTest = 20000;
 
     gLog << "filenameTest = " << filenameTest << ",   howManyTest = "