Changeset 6232 for trunk/MagicSoft/Mars

Timestamp:

02/03/05 12:42:21 (20 years ago)

Author:

mazin

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mmpi

Files:

• SupercutsONOFFClasses/MFindSupercutsONOFF.cc (modified) (18 diffs)
• SupercutsONOFFClasses/MFindSupercutsONOFF.h (modified) (5 diffs)
• SupercutsONOFFClasses/MFindSupercutsONOFFThetaLoop.cc (modified) (25 diffs)
• SupercutsONOFFClasses/MFindSupercutsONOFFThetaLoop.h (modified) (2 diffs)
• SupercutsONOFFClasses/MHFindSignificanceONOFF.cc (modified) (16 diffs)
• SupercutsONOFFClasses/MSupercutsCalcONOFF.cc (modified) (14 diffs)
• SupercutsONOFFClasses/MSupercutsCalcONOFF.h (modified) (2 diffs)
• SupercutsONOFFClasses/MTSupercutsApplied.cc (modified) (2 diffs)
• macros/SuperCutsONOFFMacro.C (modified) (12 diffs)

trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MFindSupercutsONOFF.cc

@@ -645 +645 @@
     // The default is not cut, i.e. all values (0-1) are taken
     
-    fThetaMin = 0; // in miliradians  // FIXME: change name
-    fThetaMax = 90; // new: in deg; old: (1570) in miliradians // FIXME: change name
-    
+    //fThetaMin = 0; // in miliradians  // FIXME: change name
+    //fThetaMax = 1570; // in miliradians // FIXME: change name
+    
+    fThetaMin = 0; // degrees (new in magic)
+    fThetaMax = 90; // degrees (new in magic)
+
+
 
     fAlphaSig = 20; // By default, signal is expected in alpha<20 for CT1
@@ -770 +774 @@
 
     // Angular cuts not yet implemented ...
-    Double_t ConvMMToDeg = 0.00337034;
-    
-    fDistCutLow = 0.4/ConvMMToDeg;
-    fDistCutUp = 1.5/ConvMMToDeg;
-    
-    fLengthCutLow = 0.1/ConvMMToDeg;
-    fLengthCutUp = 1/ConvMMToDeg;
-
-    fWidthCutLow = 0.07/ConvMMToDeg;
-    fWidthCutUp = 1/ConvMMToDeg;
+    fConvMMToDeg = 0.00337034;
+
+    fDistCutLow = 0.4/fConvMMToDeg;
+    fDistCutUp = 1.5/fConvMMToDeg;
+    
+    fLengthCutLow = 0.1/fConvMMToDeg;
+    fLengthCutUp = 1/fConvMMToDeg;
+
+    fWidthCutLow = 0.07/fConvMMToDeg;
+    fWidthCutUp = 1/fConvMMToDeg;
 
     // ENDTMP
@@ -1141 +1145 @@
 // in previous functions...
 
-    fThetaMin = int(ThetaMin*1000.0);
-    fThetaMax = int(ThetaMax*1000.0);
+  //    fThetaMin = int(ThetaMin*1000.0);
+  //  fThetaMax = int(ThetaMax*1000.0);
+
+
+  // in new magic theta is given in deg
+  // 0.5 is added so that the rounding to integer is correct
+  fThetaMin = int(ThetaMin + 0.5); 
+  fThetaMax = int(ThetaMax + 0.5); 
 
     fThetaRangeString = ("ThetaRange");
@@ -1148 +1158 @@
     fThetaRangeString += ("_");
     fThetaRangeString += (fThetaMax);
-    fThetaRangeString += ("mRad");
+    // fThetaRangeString += ("mRad");
+    fThetaRangeString += ("_Degrees"); // new in magic
+
 
     return kTRUE;
@@ -1168 +1180 @@
 
 
+
+    return kTRUE;
+}
+
+
+Bool_t MFindSupercutsONOFF::SetFilters(Double_t LeakageMax, Double_t DistMax, Double_t DistMin)
+{
+
+    fDistCutLow = DistMin/fConvMMToDeg;
+    fDistCutUp  = DistMax/fConvMMToDeg;
+    fLeakageMax = LeakageMax;
+    
+    *fLog << "MFindSupercutsONOFF::SetSizeRange" << endl
+          << "Data matrices will be filled with events whose MHillasSrc.fDist " << endl
+          << "is in the range " 
+          << fDistCutLow <<"-"<< fDistCutUp << " degrees" << endl
+          << "and fLeakage " << "< " << fLeakageMax << endl;
 
     return kTRUE;
@@ -1682 +1711 @@
 
 
-    // TMP
+    
     // Cuts in Size,
 
@@ -1696 +1725 @@
     SizeCutMax.SetInverted();
     
-
-
-
-    // ENDTMP
+    // Cuts in Dist
+    TString DistCutMinString ("MHillasSrc.fDist>"); 
+    DistCutMinString += fDistCutLow;
+    MContinue DistCutMin(DistCutMinString);
+    DistCutMin.SetInverted();
+
+    TString DistCutMaxString ("MHillasSrc.fDist<"); 
+    DistCutMaxString += fDistCutUp; 
+    MContinue DistCutMax(DistCutMaxString);
+    DistCutMax.SetInverted();
+
+
+    // Cuts in leakage
+    TString LeakCutMaxString ("MNewImagePar.fLeakage1<"); 
+    LeakCutMaxString += fLeakageMax; 
+    MContinue LeakCutMax(LeakCutMaxString);
+    LeakCutMax.SetInverted();
 
 
@@ -1706 +1748 @@
     Double_t prob = whichfractiontrain /(whichfractiontrain+whichfractiontest);
                    
+    
+    
+
     MFRandomSplit split(prob);
 
@@ -1735 +1780 @@
 
     tlist.AddToList(&read);
-//    tlist.AddToList(&ThetaCutMin);
-//    tlist.AddToList(&ThetaCutMax);
-    
-    //TMP
+    tlist.AddToList(&ThetaCutMin);
+    tlist.AddToList(&ThetaCutMax);
+    
+    
     tlist.AddToList(&SizeCutMin);
     tlist.AddToList(&SizeCutMax);
 
-    //ENDTMP
+    tlist.AddToList(&DistCutMin);
+    tlist.AddToList(&DistCutMax);
+    tlist.AddToList(&LeakCutMax);
+
+    
     tlist.AddToList(&split);
     tlist.AddToList(&filltrain);
@@ -1848 +1897 @@
 
 
-    // TMP
+    
     // Cuts in Size,
 
@@ -1863 +1912 @@
     
 
-
-
-    // ENDTMP
+    // Cuts in Dist
+    TString DistCutMinString ("MHillasSrc.fDist>");
+    DistCutMinString += fDistCutLow;
+    MContinue DistCutMin(DistCutMinString);
+    DistCutMin.SetInverted();
+
+    TString DistCutMaxString ("MHillasSrc.fDist<");
+    DistCutMaxString += fDistCutUp;
+    MContinue DistCutMax(DistCutMaxString);
+    DistCutMax.SetInverted();
+
+
+    // Cuts in leakage
+    TString LeakCutMaxString ("MNewImagePar.fLeakage1<");
+    LeakCutMaxString += fLeakageMax;
+    MContinue LeakCutMax(LeakCutMaxString);
+    LeakCutMax.SetInverted();
+
+
+
+
+
+    
 
 
@@ -1871 +1940 @@
     
     Double_t prob = whichfractiontrain/(whichfractiontrain + whichfractiontest);
-                   
+
+    
+
     MFRandomSplit split(prob);
 
@@ -1903 +1974 @@
 
     tlist.AddToList(&read);
-//    tlist.AddToList(&ThetaCutMin);
-//    tlist.AddToList(&ThetaCutMax);
+    tlist.AddToList(&ThetaCutMin);
+    tlist.AddToList(&ThetaCutMax);
     
        
-    //TMP
+    
     tlist.AddToList(&SizeCutMin);
     tlist.AddToList(&SizeCutMax);
 
-    //ENDTMP
+    tlist.AddToList(&DistCutMin);
+    tlist.AddToList(&DistCutMax);
+    tlist.AddToList(&LeakCutMax);
+
+    
 
     tlist.AddToList(&split);
@@ -1975 +2050 @@
 
 
+
+/// **********///
+
+// --------------------------------------------------------------------------
+//
+// Read only training matrices ON and OFF
+//
+//
+
+Bool_t MFindSupercutsONOFF::ReadMatrixTrain(const TString &filetrainON, const TString &filetrainOFF)
+{
+  //--------------------------
+  // read in training matrix ON sample
+
+  TFile filetr(filetrainON);
+  fMatrixTrain->Read("MatrixTrain");
+  fMatrixTrain->Print("SizeCols");
+
+  *fLog << "MFindSupercuts::ReadMatrixTrain; Training matrix for ON sample was read in from file '"
+        << filetrainON << "'" << endl;
+  filetr.Close();
+
+
+  // read in training matrix OFF sample
+
+  TFile filetrOFF(filetrainOFF);
+  fMatrixTrainOFF->Read("MatrixTrainOFF");
+  fMatrixTrainOFF->Print("SizeCols");
+
+  *fLog << "MFindSupercutsONOFF::ReadMatrixTrain; Training matrix for OFF sample was read in from file '"
+        << filetrainOFF << "'" << endl;
+  filetrOFF.Close();
+
+  return kTRUE;  
+
+}
+
+// --------------------------------------------------------------------------
+//
+// Read only test matrices ON and OFF
+//
+//
+
+Bool_t MFindSupercutsONOFF::ReadMatrixTest(const TString &filetestON, const TString &filetestOFF)
+{
+  //--------------------------
+  // read in testing matrix ON sample
+
+//--------------------------
+  // read in test matrix for ON sample
+  
+  TFile filete(filetestON);
+  fMatrixTest->Read("MatrixTest");
+  fMatrixTest->Print("SizeCols");
+
+  *fLog << "MFindSupercuts::ReadMatrixTest; Test matrix for ON sample was read in from file '"
+        << filetestON << "'" << endl;
+  filete.Close();
+  
+
+   //--------------------------
+  // read in test matrix for OFF sample
+  
+  TFile fileteOFF(filetestOFF);
+  fMatrixTestOFF->Read("MatrixTestOFF");
+  fMatrixTestOFF->Print("SizeCols");
+
+  *fLog << "MFindSupercutsONOFF::ReadMatrixTest; Test matrix for OFF sample was read in from file '"
+        << filetestOFF << "'" << endl;
+  filete.Close();
+
+
+
+  return kTRUE;  
+
+}
+
+
 /// **********///
 
@@ -1999 +2152 @@
   //--------------------------
   // read in test matrix
-
+  
   TFile filete(filetest);
   fMatrixTest->Read("MatrixTest");
@@ -2007 +2160 @@
         << filetest << "'" << endl;
   filete.Close();
-
+  
   return kTRUE;  
 }
@@ -2032 +2185 @@
   //--------------------------
   // read in test matrix
-
+  
   TFile filete(filetest);
   fMatrixTestOFF->Read("MatrixTestOFF");
@@ -2040 +2193 @@
         << filetest << "'" << endl;
   filete.Close();
-
+  
   return kTRUE;  
 }

trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MFindSupercutsONOFF.h

@@ -59 +59 @@
     Double_t  fWidthCutUp;
 
+    Double_t  fLeakageMax;
+
+    Double_t fConvMMToDeg;
 
     // ENDTMP
@@ -85 +88 @@
   Double_t fAlphaBkgMax; 
   
-  Int_t fThetaMin; // Cuts in ThetaOrig.fVal (in mili radians!!!)
-  Int_t fThetaMax; // Cuts in ThetaOrig.fVal (in mili radians !!!)
+  Int_t fThetaMin; // Cut in Theta value (in degrees)
+  Int_t fThetaMax; // Cuts in Theta value (in degrees)
   TString fThetaRangeString;
 
@@ -314 +317 @@
   Bool_t CheckAlphaSigBkg();
 
+
+
+
+
   void SetUseFittedQuantities (Bool_t b)
       {fUseFittedQuantities = b;}
@@ -371 +378 @@
   Bool_t ReadMatrixOFF( const TString &filetrainOFF, const TString &filetestOFF);
 
+
+  Bool_t ReadMatrixTrain(const TString &filetrainON, const TString &filetrainOFF);
+  Bool_t ReadMatrixTest(const TString &filetestON, const TString &filetestOFF);
+
+
+
   Bool_t ComputeNormFactorTrain ();
   Bool_t ComputeNormFactorTest ();
@@ -395 +408 @@
 
   Bool_t SetSizeRange(Double_t SizeMin, Double_t SizeMax);
+
+  Bool_t SetFilters(Double_t LeakageMax, Double_t DistMax, Double_t DistMin);
+
 
 

trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MFindSupercutsONOFFThetaLoop.cc

@@ -110 +110 @@
     //---------------------------
 
-
+    // *********************************
+    // ** NOT TRUE ANY LONGER... ** //
     // Cuts (low and up) in variable ThetaOrig.fVal
     // The default is not cut, i.e. all values (0-1) are taken
     // fThetaOrig.fVal is measured in Radians; thus 1 = 57 degrees.
 
+    //    fThetaMin = 0.0;
+    //    fThetaMax = 1.0;
+    // *********************************
+
+
+    // For MAGIC new (July 2004) Theta is measured in degrees
+
     fThetaMin = 0.0;
-    fThetaMax = 1.0;
+    fThetaMax = 90.0;
+
+
 
     fAlphaSig = 20; // By default, signal is expected in alpha<20 for CT1
@@ -545 +555 @@
 
 
-    // Give a "soft warning" if Vector has increasing values 
-    // of CosThetas; which means decreasing values of Thetas
-    /*
-    if (fCosThetaRangeVector[0] <= fCosThetaRangeVector[1])
-    {
-        *fLog << "MFindSupercutsONOFFThetaLoop::SetCosThetaRangeVector; " 
-              << endl
-              << "Values in vector fCosThetaRangeVector are in increasing "
-              << "order, i.e., values in theta will be in decreasing order"
-              << endl
-              << "DO NOT forget it !!" << endl;
-    }
-
-    */
+    
 
     // fCosThetaBinCenterVector is defined
@@ -622 +619 @@
     // Now costheta values have to be converted to Radians
     // Function TMath::ACos(x) gives angle in Radians
-
+    
     
     fThetaMin = TMath::ACos(fThetaMin);
     fThetaMax = TMath::ACos(fThetaMax);
+
+   
+    // For new MAGIC (July 2004), angles are given in degrees.
+    // So I convert angles into degrees by multiplying by a the 
+    // conversion rad/degree
+
+    fThetaMin = fThetaMin * TMath::RadToDeg();
+    fThetaMax = fThetaMax * TMath::RadToDeg();
+
+
 
     // If fThetaMin > fThetaMax means that fCosThetaRangeVector
@@ -644 +651 @@
           << "MFindSupercutsONOFFThetaLoop::SetThetaRange; " << endl
           << "Theta range for this iteration is set to "
-          << fThetaMin << "-" << fThetaMax << " Radians." << endl 
+          << fThetaMin << "-" << fThetaMax << " degrees." << endl 
           << "-------------------------------------------------"<< endl;
     
@@ -932 +939 @@
     
 
-    // Histograms references are removed from the current directory
-
-    /* IT DOES NOT WORK !!! THESE STATEMENTS PRODUCE SEGMENTATION FAULT.
-    fNormFactorTrainHist -> SetDirectory(NULL);
-    
-    fNormFactorTestHist -> SetDirectory(NULL);
-    
-    fSigmaLiMaTrainHist -> SetDirectory(NULL);
-    
-    fSigmaLiMaTestHist -> SetDirectory(NULL);
-    
-    fSigmaLiMaTestHist -> SetDirectory(NULL);
-    
-    fNexTestHist -> SetDirectory(NULL);
-    */
+    
 
 
@@ -1059 +1052 @@
     Double_t fofftest)
 {
+
+  // Check that fractions set by user are within the expected limits
+  
+  if (fontrain < 0.00 || fontrain > 1.00 ||
+      fontest < 0.00 || fontest > 1.00 ||
+      fofftrain < 0.00 || fofftrain > 1.00 ||
+      fofftest < 0.00 || fontest > 1.00)
+    {
+      *fLog << err 
+            << "MFindSupercutsONOFFThetaLoop::SetFractionTrainTestOnOffEvents: " 
+            << endl
+            << "Fractions of ON/OFF data for TRAIN/TEST sample which were specified by user " << endl 
+            << "are outside the natural limits 0.0-1.0. " << endl
+            << "Set reasonable fraction numbers for TRAIN/TEST please. The program will be aborted."
+            << endl;
+      
+      exit(1);
+
+
+    }
+  
+
+  if((fontrain+fontest) > 1.00 ||
+     (fofftrain+fofftest) >1.00)
+    
+    {
+      *fLog << err 
+            << "MFindSupercutsONOFFThetaLoop::SetFractionTrainTestOnOffEvents: " 
+            << endl
+            << "The sum of the fractions of TRAIN/TEST for ON or OFF data which were specified by user " << endl 
+            << "are larger than 1.0. " << endl
+            << "Set reasonable fraction numbers for TRAIN/TEST please. The program will be aborted."
+            << endl;
+      
+      exit(1);
+
+
+    }
+
 
     fWhichFractionTrain = fontrain;
@@ -1154 +1186 @@
 
     // Path + "OptSCParametersONOFF" + "ThetaRange" 
-    // + int {(fThetaMin_fThetaMax)*1000} + ".root"
+    // + int {(fThetaMin_fThetaMax)} + ".root"
 
 
@@ -1160 +1192 @@
 
     // Path + "ON(OFF)Data" + "Train(Test)" + "Matrix" + Fraction 
-    // + "ThetaRange" + int {(fThetaMin_fThetaMax)*1000} + ".root"
+    // + "ThetaRange" + int {(fThetaMin_fThetaMax)} + ".root"
 
 
@@ -1176 +1208 @@
 
 
-    Int_t ThetaMinTimes1000 = 0;
-    Int_t ThetaMaxTimes1000 = 0;
+    //    Int_t ThetaMinTimes1000 = 0;
+    // Int_t ThetaMaxTimes1000 = 0;
+
+    // For MAGIC new (July 2004) Theta is given in deg
+
+    Int_t ThetaMinInt = 0;
+    Int_t ThetaMaxInt = 0;
+
+
     Int_t FractionONTrainTimes100 = 0;
     Int_t FractionONTestTimes100 = 0;
@@ -1200 +1239 @@
         
         
-        ThetaMinTimes1000 = int(fThetaMin*1000);
-        ThetaMaxTimes1000 = int(fThetaMax*1000);
-        sprintf(ThetaRangeString, "ThetaRange%d_%dmRad", 
-                ThetaMinTimes1000, ThetaMaxTimes1000);
+        //ThetaMinTimes1000 = int(fThetaMin*1000);
+        //ThetaMaxTimes1000 = int(fThetaMax*1000);
+        //sprintf(ThetaRangeString, "ThetaRange%d_%dmRad", 
+        //      ThetaMinTimes1000, ThetaMaxTimes1000);
+
+        // For New MAGIC (July 2004)
+        // Theta is given in deg, no need to multiply by 1000
+        
+        // 0.5 is added so that the rounding to integer is correct
+        ThetaMinInt = int(fThetaMin + 0.5);
+        ThetaMaxInt = int(fThetaMax + 0.5);
+        sprintf(ThetaRangeString, "ThetaRange%d_%d_Degrees", 
+                ThetaMinInt, ThetaMaxInt);
+
+        
 
         fThetaRangeStringVector[i] = (ThetaRangeString);
@@ -1211 +1261 @@
         // endtmp
 
-        Double_t tmpdouble = 0.0;
-        tmpdouble = fWhichFractionTrain*100 - int(fWhichFractionTrain*100); 
-        
-        FractionONTrainTimes100 = tmpdouble < 0.5 ? int(fWhichFractionTrain*100):
-            int(fWhichFractionTrain*100) + 1;
-        
-
-        tmpdouble = fWhichFractionTest*100 - int(fWhichFractionTest*100);
-
-        FractionONTestTimes100 = tmpdouble < 0.5 ? int(fWhichFractionTest*100):
-            int(fWhichFractionTest*100) + 1;
-        
-        tmpdouble = fWhichFractionTrainOFF*100 - int(fWhichFractionTrainOFF*100);
-
-        FractionOFFTrainTimes100 = tmpdouble < 0.5 ? int(fWhichFractionTrainOFF*100):
-            int(fWhichFractionTrainOFF*100) + 1;
-
-        tmpdouble = fWhichFractionTestOFF*100 - int(fWhichFractionTestOFF*100);
-        FractionOFFTestTimes100 = tmpdouble < 0.5 ? int(fWhichFractionTestOFF*100):
-            int(fWhichFractionTestOFF*100) + 1;
+        // 0.5 is added so that the rounding to integer is correct
+
+        FractionONTrainTimes100 = int(fWhichFractionTrain*100 + 0.5);
+        FractionONTestTimes100 = int(fWhichFractionTest*100 + 0.5);
+        FractionOFFTrainTimes100 = int(fWhichFractionTrainOFF*100 + 0.5);
+        FractionOFFTestTimes100 = int(fWhichFractionTestOFF*100 + 0.5);
+          
+
         
         
@@ -1525 +1563 @@
         FindSupercuts.SetSizeRange(fSizeCutLow, fSizeCutUp);
 
+        FindSupercuts.SetFilters(fLeakMax, fDistMax, fDistMin);
 
 
@@ -1627 +1666 @@
                   << "Reading Data matrices from root files... " << endl;
 
+            /*
             FindSupercuts.ReadMatrix(fTrainMatrixONFilenameVector[i], 
                                      fTestMatrixONFilenameVector[i]);
             FindSupercuts.ReadMatrixOFF(fTrainMatrixOFFFilenameVector[i], 
                                         fTestMatrixOFFFilenameVector[i]);
+            */
+
+            FindSupercuts.ReadMatrixTrain(fTrainMatrixONFilenameVector[i], 
+                                          fTrainMatrixOFFFilenameVector[i]);
+
+            if (fTestParameters)
+              {
+                FindSupercuts.ReadMatrixTest(fTestMatrixONFilenameVector[i], 
+                                             fTestMatrixOFFFilenameVector[i]);
+              }
 
         }
@@ -2325 +2375 @@
         
             // TEMP    
-            // cout << NexVSAlphaNameTmp << endl;
-            // cout << SignificanceVSAlphaNameTmp << endl;
+            //cout << NexVSAlphaNameTmp << endl;
+            //cout << SignificanceVSAlphaNameTmp << endl;
             // ENDTEMP
               
@@ -2388 +2438 @@
 
             AlphaOFFAfterCuts = (TH1F*) rootfile.Get(TmpAlphaOFFHistoName);
+
+           
+
             
 
@@ -2885 +2938 @@
 {
 
+
+   *fLog << endl << endl 
+        << "****************************************************************************" << endl
+        << "Combining number of excess events and significance for the  " << endl
+        << "selected bins in theta angle using the function " << endl
+        << " MFindSupercutsONOFFThetaLoop::ComputeOverallSignificance() " << endl
+        << "****************************************************************************" << endl
+        << endl;
+  
+
+
+ 
+
+
     // First of all let's check that we have the proper the "drinks" for the party
 
@@ -2917 +2984 @@
     // Remember that histogram index START AT 1 !!!
 
+
+
+    // Histograms that will contain the alpha distributions 
+    // for the combined theta bins
+
+
     TH1F* CombinedAlphaTrainON = NULL;
-    TString CombinedAlphaTrainONName = ("CombinedAlphaTrainON");
+    TString CombinedAlphaTrainONName = ("TrainSampleAlphaONCombinedThetaBins");
     TH1F* CombinedAlphaTrainOFF = NULL;
-    TString CombinedAlphaTrainOFFName = ("CombinedAlphaTrainOFF");
+    TString CombinedAlphaTrainOFFName = ("TrainSampleAlphaOFFCombinedThetaBins");
    
     TH1F* CombinedAlphaTestON = NULL;
-    TString CombinedAlphaTestONName = ("CombinedAlphaTestON");
+    TString CombinedAlphaTestONName = ("TestSampleAlphaONCombinedThetaBins");
    
     TH1F* CombinedAlphaTestOFF = NULL;
-    TString CombinedAlphaTestOFFName = ("CombinedAlphaTestOFF");
+    TString CombinedAlphaTestOFFName = ("TestSampleAlphaONCombinedThetaBins");
    
-
-    
+    
+
 
     
@@ -3089 +3162 @@
     {
 
+
+      
+      *fLog << endl 
+        << "****************************************************************************" << endl
+        << "Combining data for the TRAIN sample " << endl
+        << "****************************************************************************" << endl
+        << endl;
+
         // Retrieving NormFactorTrainHisto from root file. 
         
@@ -3489 +3570 @@
         CombinedAlphaTrainOFF->SetEntries(tmpint);
 
-
+        // Check that the number of events in the combined normalized histogram is 
+        // that of the initial sample
+
+        /*
         cout << "SILLY INFO FOR TRAIN SAMPLE: Total number of events Before nomralization and " 
              << " re-re normalized" << endl
              << TotalNumberOFFTrain << " - " << EventCounter << endl;
 
+        */
     }
 
@@ -3506 +3591 @@
     if (CombineTestData)
     {
+
+      *fLog << endl 
+        << "****************************************************************************" << endl
+        << "Combining data for the TEST sample " << endl
+        << "****************************************************************************" << endl
+        << endl;
+
 
         // Retrieving NormFactorTestHisto from root file. 
@@ -3930 +4022 @@
 
 
+        // Check that the number of events in the combined normalized histogram is 
+        // that of the initial sample
+
+        /*
         cout << "SILLY INFO FOR TEST SAMPLE: Total number of events Before nomralization and " 
              << " re-re normalized" << endl
              << TotalNumberOFFTest << " - " << EventCounter << endl;
+
+        */
+
 
         /*
@@ -3982 +4081 @@
     if (SuccessfulThetaBinsVector.GetSize() >= 1)
     {
+
+
+      *fLog << endl 
+        << "**************************************************************************************" << endl
+        << "Computation of excess events and singnificance for the combined alpha histograms " << endl
+        << "**************************************************************************************" << endl
+        << endl;
+
         // There is, at least, ONE theta bin for which optimization of 
         // supercuts was SUCCESSFUL, and thus, it is worth to 
@@ -4026 +4133 @@
 
             // Name of psfile where Alpha plot will be stored
-            TString psfilename = ("TRAINSampleCombined"); 
+          
+            TString psfilename = (fPathForFiles);
+            psfilename += ("TRAINSampleCombinedThetaBins"); 
             
             findsig.FindSigmaONOFF(CombinedAlphaTrainON, CombinedAlphaTrainOFF, 
@@ -4053 +4162 @@
 
             // Name of psfile where Alpha plot will be stored
-            TString psfilename = ("TESTSampleCombined"); 
+          TString psfilename = (fPathForFiles);
+          psfilename += ("TESTSampleCombinedThetaBins"); 
             
             findsig.FindSigmaONOFF(CombinedAlphaTestON, CombinedAlphaTestOFF, 
@@ -4075 +4185 @@
         
 
+        *fLog << endl 
+              << "**************************************************************************************" << endl
+              << "Alpha distributions (Train/Test and ON/OFF) for the combined  " << endl
+              << " theta bins are stored in the root file " << endl 
+              << fAlphaDistributionsRootFilename << endl
+              << "**************************************************************************************" << endl
+        << endl;
+
+
+
+        TFile rootfiletowrite(fAlphaDistributionsRootFilename, "UPDATE");
+
+
+        
+        // Histograms that will contain the Norm factors for the 
+        // combined TRAIN/TEST theta bins are created
+
+        
+        TString CombinedNormFactorTrainName = ("TrainSampleNormFactorCombinedThetaBins");
+        TH1F CombinedNormFactorTrain (CombinedNormFactorTrainName, 
+                                      CombinedNormFactorTrainName, 
+                                      1, 0.0, 1.0);
+
+        TString CombinedNormFactorTestName = ("TestSampleNormFactorCombinedThetaBins");
+        TH1F CombinedNormFactorTest (CombinedNormFactorTestName, 
+                                     CombinedNormFactorTestName, 
+                                     1, 0.0, 1.0);
+
+
+        if (CombineTrainData)
+        {
+
+        // Histogram that will contain the Norm factors for the 
+        // combined TRAIN theta bins are filled
+        
+        
+        CombinedNormFactorTrain.Fill(0.5, MeanNormFactorTrain);
+
+
+        // Train histos are written to root file
+        
+        CombinedAlphaTrainON -> Write();
+        CombinedAlphaTrainOFF -> Write();
+        CombinedNormFactorTrain.Write();
+
+
+
+        }
+
+
+        if (CombineTestData)
+        {
+
+          // Histogram that will contain the Norm factors for the 
+        // combined TEST theta bins are filled
+
+        
+          CombinedNormFactorTest.Fill(0.5, MeanNormFactorTest);
+
+          // Test histos are written to root file
+
+        CombinedAlphaTestON -> Write();
+        CombinedAlphaTestOFF -> Write();
+        CombinedNormFactorTest.Write();
+        
+
+        }
+
+
+        // Root file is closed
+
+        rootfiletowrite.Close();
+
         
     }

trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MFindSupercutsONOFFThetaLoop.h

@@ -73 +73 @@
   Double_t fSizeCutUp;
 
+  Double_t fLeakMax;
+  Double_t fDistMax;
+  Double_t fDistMin;
   
  // Variables for binning of alpha plots
@@ -373 +376 @@
    void SetSizeRange(Double_t SizeMin, Double_t SizeMax)
         {fSizeCutLow = SizeMin; fSizeCutUp = SizeMax; }
+
+//DM:
+   void SetFilters(Double_t LeakageMax, Double_t DistMax, Double_t DistMin)
+        {fLeakMax = LeakageMax; fDistMax = DistMax; fDistMin = DistMin; }
 
    // Function that loops over the theta ranges defined by 

trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MHFindSignificanceONOFF.cc

@@ -132 +132 @@
 
       //-----------------------------
-      // ignore unwanted points
-      if (error > 1.e19)
+      // ignore unwanted points      if (error > 1.e19)
         continue;
 
@@ -373 +372 @@
 
     // allow rebinning of the alpha plot
-    fRebin = kTRUE;
+//    fRebin = kTRUE;
+
+    fRebin = kFALSE;
 
     // allow reducing the degree of the polynomial
@@ -396 +397 @@
 
     fPrintResultsOntoAlphaPlot = kTRUE;
+    //fPrintResultsOntoAlphaPlot = kFALSE;
 
 
@@ -872 +874 @@
   if (fDraw)
   {
-
-
        // TEMPORALLY I will plot fHistOFF and fHistOFFNormalized (with the fits)
 
@@ -1046 +1046 @@
 
       // count bins with zero entry
-      if (content <= 0.0)
-        {
-          fNzeroOFF++;
-          // The error of the bin is set to a huge number, 
-          // so that it does not have any weight in the fit
-          fHistOFF->SetBinError(i, dummy); 
-        }
+      if (content < 0.0001)
+        fNzeroOFF++;
+     
       // set minimum error
-      if (content < 9.0)
+      // modified by K. Mase
+      if (content < 0.0)
       {
         fMlowOFF += 1;
@@ -1081 +1078 @@
     mean /= ((Double_t) fMbinsOFF);
     rms  /= ((Double_t) fMbinsOFF);
-  }
-
-  rms = sqrt( rms - mean*mean );
+    rms = sqrt( rms - mean*mean );
+  }
+
+  
 
   // if there are no events in the background region
@@ -1155 +1153 @@
 
   fConstantBackg = kFALSE;
-
-  // Condition for disabling the fitting procedure and 
-  // assuming a constant background (before Nov 2004)
-
-  // if ( fNzeroOFF > 0  ||  (Double_t)fMlowOFF>0.05*(Double_t)fMbinsOFF )
-
-
-  // Condition for disabling the fitting procedure and 
-  // assuming a constant background (After Nov 01 2004)
-  // I softened the condition to allow the fit also in situations 
-  // where the reduction of the background is such that very 
-  // few events survived; which is 
-  // Specially frequent with Random Forest at high Sizes)
-
-   if ( (Double_t)fNzeroOFF > 0.1*(Double_t)fMbinsOFF ||  
-       (Double_t)fMlowOFF > 0.2*(Double_t)fMbinsOFF )
+  if ( fNzeroOFF > 0  ||  (Double_t)fMlowOFF>0.05*(Double_t)fMbinsOFF ) 
   {
     *fLog << "MHFindSignificanceONOFF::FitPolynomialOFF; polynomial fit not possible,  fNzeroOFF, fMlowOFF, fMbinsOFF = "
@@ -1201 +1184 @@
     Double_t val, err;
     val = mean;
-    err = rms; // sqrt( mean / (Double_t)fMbinsOFF );
+    err = rms;
+    
 
     fPolyOFF->SetParameter(0, val);
@@ -1571 +1555 @@
 
       // count bins with zero entry
-      if (content <= 0.0)
-        {
-          fNzero++;
-          // The error of the bin is set to a huge number, 
-          // so that it does not have any weight in the fit
-          fHistOFF->SetBinError(i, dummy); 
-        }
-
+      /*if (content <= 0.0)
+        fNzero++;*/
+     
       // set minimum error
-      if (content < 9.0)
+      /*if (content < 9.0)
       {
         fMlow += 1;
         fHist->SetBinError(i, 3.0);
-      }
+        }*/
 
       //*fLog << "Take : i, content, error = " << i << ",  " 
@@ -1681 +1660 @@
 
   fConstantBackg = kFALSE;
-
-  // Condition for disabling the fitting procedure and 
-  // assuming a constant background (before Nov 2004)
-
-  // if ( fNzero > 0  ||  (Double_t)fMlow>0.05*(Double_t)fMbins )
-
-
-  // Condition for disabling the fitting procedure and 
-  // assuming a constant background (After Nov 01 2004)
-  // I softened the condition to allow the fit also in situations 
-  // where the reduction of the background is such that very 
-  // few events survived; which is 
-  // Specially frequent with Random Forest at high Sizes)
-
-  if ( (Double_t)fNzero > 0.1*(Double_t)fMbins ||  
-       (Double_t)fMlow > 0.2*(Double_t)fMbins )
-
+  if ( fNzero > 0  ||  (Double_t)fMlow>0.05*(Double_t)fMbins ) 
   {
     *fLog << "MHFindSignificanceONOFF::FitPolynomial; polynomial fit not possible,  fNzero, fMlow, fMbins = "
@@ -1728 +1691 @@
     Double_t val, err;
     val = mean;
-    err = rms; // sqrt( mean / (Double_t)fMbins );
+    err = sqrt( mean / (Double_t)fMbins );
 
     fPoly->SetParameter(0, val);
@@ -2776 +2739 @@
      }
      sum  *= fac*fac;
-     
+
      if (sum < 0.0)
      {
@@ -3231 +3194 @@
 
     CanvasHistOFFNormalized -> Update();
-
-
-    */
+*/
 
     return kTRUE;
@@ -3292 +3253 @@
             fpoly->DrawCopy("same");
         }
+
+        fHistOFFNormalized->SetLineColor(kRed);
+        fHistOFFNormalized->DrawCopy("same");
+
     }
 
@@ -3316 +3281 @@
       {
         // 10, 1 is white and solid
-        fbackg->SetLineColor(10);
+        fbackg->SetLineColor(kRed);
         fbackg->SetLineStyle(1);
         fbackg->SetLineWidth(4);
-        // fbackg->DrawCopy("same"); I do not want to draw it... already too many things.
+        //fbackg->DrawCopy("same"); // I do not want to draw it... already too many things.
       }
     }
@@ -3495 +3460 @@
             // fPsFilenameString
             
-            filename += ("AlphaPlotAfterSupercuts.ps");
+            filename += ("AlphaPlot.ps");
             cout << filename << endl;   
             fCanvas -> SaveAs(filename);

trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MSupercutsCalcONOFF.cc

@@ -133 +133 @@
     fWidthUpperLimit = 0.4;
 
+    // Temp upper limit to get rid of triangular events
+    fLog10ConcUpperLimit = -0.35; 
+
     fLeakage1UpperLimit = 0.25; 
 
@@ -138 +141 @@
 
     fDistLowerLimit = 0.1;
-    fLengthLowerLimit = 0.09;
-    fWidthLowerLimit = 0.05;
+    fLengthLowerLimit = 0.05; // values at October 21th 2004
+    fWidthLowerLimit = 0.03;  // values at October 21th 2004
 
     
@@ -237 +240 @@
 
 } 
+
+
+// Function implementing a filter for muon induced images 
+// (Values from Keichi cuts October 18th 2004)
+
+// Function returns kTRUE if the event does NOT survive the 
+// muon filter. So kTRUE means event MUST be rejected
+
+
+Bool_t MSupercutsCalcONOFF::MuonFilter (Double_t Size, 
+                                        Double_t Length)
+{
+  // Parametrization for cut in length. 
+  // Event is rejected if 
+  // Log10(Length) < k1*(log10(Size)) + k2   (Size<10000photons)
+  // k1 = 0.286; k2 = -1.91 (October 18th 2004)
+
+  // Length < 0.2 deg (Size > 10000)
+
+
+  // Parametrization for cut in Length/Size 
+  // Event is rejected if
+  // log10(LenOverSize) < k3*(log10(Size)) + k4 (Size<4000photons)
+  // k3 = -0.780; k4 = -1.71
+
+
+  Double_t LengthUpperLimit = 0.0;
+  Double_t LengthOverSizeUpperLimit = 0.0;
+  
+
+  Double_t SizeUpperLimitForCutInLength = 10000; // in photons
+  Double_t SizeUpperLimitForCutInLengthOverSize = 4000; // in photons
+  
+  Double_t k1 =  0.286;
+  Double_t k2 =  -1.91;
+  Double_t k3 =  -0.78;
+  Double_t k4 =  -1.71;
+  
+
+  
+
+
+  if (Size <= SizeUpperLimitForCutInLength)
+    {
+      LengthUpperLimit = k1 * TMath::Log10(Size) + k2;
+      LengthUpperLimit = TMath::Power(10, LengthUpperLimit);
+    }
+  else
+    {
+      LengthUpperLimit = 0.2;
+    }
+
+  // Apply cut in Length
+  
+  if (Length <  LengthUpperLimit)
+        return kTRUE;
+
+
+  
+  if (Size < SizeUpperLimitForCutInLengthOverSize)
+    {
+      LengthOverSizeUpperLimit = k3 * TMath::Log10(Size) + k4;
+      LengthOverSizeUpperLimit = TMath::Power(10,  LengthOverSizeUpperLimit);
+  
+      if (Length/Size <  LengthOverSizeUpperLimit)
+        return kTRUE;
+    }
+
+
+   return kFALSE;
+
+}
+
+
+
+
+
 
 
@@ -340 +420 @@
     //    ls2: ls^2
     //     ct: Cos(ZA.) - Cos(fThetaOffset)
-    //    dd2: (DIST- fDistOffset)^2
+    //    dd2: DIST^2 - fDistOffset^2
+
+  //     lconc: log10CONC   
+
 
     Double_t limit;
@@ -488 +571 @@
 
 
+
+    // The parametrization of upper cut in CONC is just provisional. 
+    // The aim is to get rid of triangular events. 
+    // The parametrization is 
+    // Log10(Conc)Uppercut = m*(log(Size) - log(SizeOffset)) + b
+
+
+    Double_t log10conc = TMath::Log10(conc);
+
+
     // in the parameterization of the cuts 
     // the dist parameter used is the one computed from the source position, 
@@ -508 +601 @@
     const Double_t asym    = asym0   * fMm2Deg;
 
-
+    
+    // MuuonLikeEvent is a variable which is set 
+    // to kFALSE/kTRUE by the filter cuts implemented 
+    // in MuonFilter function. If kTRUE, the event 
+    // will be rejected.
+
+// DM:
+    //Bool_t MuonLikeEvent = MuonFilter (size, length); 
+    
     
     // computation of the cut limits
@@ -524 +625 @@
     Double_t asymlow  = CtsMCut (fSuper->GetAsymLo(),   dmls, dmcza, dmls2, dd2);
 
-    Double_t concup   = CtsMCut (fSuper->GetConcUp(),   dmls, dmcza, dmls2, dd2);
-    Double_t conclow  = CtsMCut (fSuper->GetConcLo(),   dmls, dmcza, dmls2, dd2);
+    Double_t log10concup   = CtsMCut (fSuper->GetConcUp(),   dmls, dmcza, dmls2, dd2);
+    Double_t log10conclow  = CtsMCut (fSuper->GetConcLo(),   dmls, dmcza, dmls2, dd2);
 
     Double_t  leakageup = CtsMCut (fSuper->GetLeakage1Up(),dmls, dmcza, dmls2, dd2);
@@ -533 +634 @@
     Double_t lengthoverwidth = length/width;
 
+
+    // Cut in CONC is parametrized
+
     Double_t hadronness = 0.0;
 
 
-    // If the cut values computed before for Dist, Length and Width exceed the upper limits set 
-    // by the user through function MSupercutsCalcONOFF::SetHillasDistLengthWidthUpperLowerLimits,
-    // (or the default values defined in constructor); such cut values are replaced by the 
+    // If the cut values computed before for Dist, 
+    //  Length and Width exceed the upper limits set 
+    // by the user through function 
+    // MSupercutsCalcONOFF::SetHillasDistLengthWidthUpperLowerLimits,
+    // (or the default values defined in constructor); such cut values 
+    // are replaced by the 
     // the limit values
 
@@ -576 +683 @@
       }
 
+    if (log10concup > fLog10ConcUpperLimit)
+      {
+        log10concup = fLog10ConcUpperLimit;
+      }
+
     
     if (leakageup > fLeakage1UpperLimit)
@@ -582 +694 @@
       }
     
-
-   
+    
 
 
     // Upper cut in leakage 1 also implemented
 
-    
-    
-   
     if (// 
         newdist > distup ||
@@ -600 +708 @@
         leakage > leakageup ||
         leakage < leakagelow ||
-        lengthoverwidth < fLengthOverWidthUpperLimit
+        lengthoverwidth < fLengthOverWidthUpperLimit ||
         //
         //asym    < asymup &&
@@ -608 +716 @@
         //dist    > distlow &&
         
-        //conc    < concup &&
-        //conc    > conclow &&  
+        log10conc    > log10concup 
+        // log10conc    < log10conclow  
         //
+        //log10conc    > log10concup ||
+        // DM MuonLikeEvent // if KTRUE, event is rejected      
         ) 
           
@@ -641 +751 @@
         ShowerParams[3] = newdist;
         ShowerParams[4] = asym;
-        ShowerParams[5] = conc;
+        ShowerParams[5] = log10conc;
         ShowerParams[6] = leakage;
         ShowerParams[7] = theta;
@@ -666 +776 @@
         SuperCutParams[6] = asymup;
         SuperCutParams[7] = asymlow;
-        SuperCutParams[8] = concup;
-        SuperCutParams[9] = conclow;
+        SuperCutParams[8] = log10concup;
+        SuperCutParams[9] = log10conclow;
         SuperCutParams[10] = leakageup;
         SuperCutParams[11] = leakagelow;

trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MSupercutsCalcONOFF.h

@@ -94 +94 @@
     Double_t fWidthUpperLimit;
 
+    Double_t fLog10ConcUpperLimit;
+
      Double_t fDistLowerLimit;
     Double_t fLengthLowerLimit;
@@ -149 +151 @@
                                                   Double_t widthlow);
 
+      // Function implementing a filter for muon induced images 
+// (Keichi cuts October 18th 2004)
+
+  Bool_t MuonFilter (Double_t Size, 
+                     Double_t Length);
+
+
+
     void InitMapping(MHMatrix *mat); // use quantity ThetaOrig.fVal at theta
     

trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MTSupercutsApplied.cc

@@ -130 +130 @@
     StringVector[6] = ("AsymUp/D");
     StringVector[7] = ("AsymLow/D");
-    StringVector[8] = ("ConcUp/D");
-    StringVector[9] = ("ConcLow/D");
+    StringVector[8] = ("Log10ConcUp/D");
+    StringVector[9] = ("Log10ConcLow/D");
     StringVector[10] = ("LeakageUp/D");
     StringVector[11] = ("LeakageLow/D");
@@ -155 +155 @@
     StringVector2[3] = ("NewDist/D");
     StringVector2[4] = ("Asym/D");
-    StringVector2[5] = ("Conc/D");
+    StringVector2[5] = ("Log10Conc/D");
     StringVector2[6] = ("Leakage/D");
     StringVector2[7] = ("Theta/D");

trunk/MagicSoft/Mars/mtemp/mmpi/macros/SuperCutsONOFFMacro.C

@@ -8 +8 @@
 gROOT -> Reset();
 
-void SuperCutsONOFFMacro()
+void SuperCutsONOFFMacroNew()
 {
   gLog.SetNoColors();
@@ -15 +15 @@
 
     // From magicserv01
-    TString  ONDataFilename("/.magic/data16a/mazin/data/Mrk421/2004_04_22/4slices/Hillas_20040422_4sl_time_clean/Mrk421_*_HillasON.root");
+//    TString  ONDataFilename("/.magic/data16a/mazin/data/mcdata/2004_06_30/*HillasON.root");
                                        
-    TString  OFFDataFilename("/.magic/data16a/mazin/data/Mrk421/2004_04_22/4slices/Hillas_20040422_4sl_time_clean/Mrk421_*_HillasOFF.root");
-
-    
-    TString PathForFiles ("/mnt/magicserv01/scratch/David/SillyTestForCommiting_July20_2004/");
-
-    
-    // **********************************************
+//    TString  OFFDataFilename("/.magic/data16a/mazin/data/mcdata/2004_06_30/*HillasOFF.root");
+
+//   TString  ONDataFilename("/.magic/magicserv01/scratch/Period21/HillasFiles/CrabNebula_20040914.root");
+//   TString  ONDataFilename("/.magic/magicserv01/scratch/Period21/HillasFiles/CrabNebula_20040921.root");
+//   TString  ONDataFilename("/.magic/magicserv01/scratch/Period21/HillasFiles/CrabNebula_20040922.root");
+
+//   TString  ONDataFilename("/.magic/magicserv01/scratch/Period21/HillasFiles/CrabNebula_2004*.root");
+  
+   //TString  ONDataFilename("/.magic/data22a/mazin/HillasParam/Period24/2004_12_18/Crab_40_25_1_AbsHillas.root");
+   TString  ONDataFilename("/.magic/data04a/nadia/Analysis/starfiles/DF/40_25_1/CrabNebula_*_Hillas_Abs.root");
+   TString  OFFDataFilename("/.magic/data04a/nadia/Analysis/starfiles/DF/40_25_1/OffCrab1*Hillas_Abs.root");
+   //TString ONDataFilename("/.magic/data03a/mazin/results/2004_09_21/CrabNebulaNadiaHillas.root");  
+   //TString OFFDataFilename("/.magic/data03a/mazin/results/2004_09_21/CrabNebulaNadiaHillas.root");  
+
+
+    TString PathForFiles ("/.magic/data03a/mazin/results/Crab/DF/CrabNadia/SuperCuts/Size2000/");  
+
+//    TString PathForFiles ("/.magic/magicserv01/scratch/hbartko/SC/"); // was: Daniel/SuperCuts/Mrk421/2004_04_22/4slices_3520_nc/E800_1200_Opt_MC_Test/
+
+    
+
     // Boolean variables defining the job of the 
     // macro
-    // **********************************************
-
-    
-
+ 
     // Boolean variable that decides wether data is read from files specified above
     // (ON/OFF) or read from already existing Matrices (which are obviously stored
@@ -36 +47 @@
     // values must be specified by user.
 
-    // kTRUE reads alredy existing matrices, and kFALSE read data and produce matrices.
-
-    Bool_t ReadMatrixFromRootFiles = kTRUE; 
+    Bool_t ReadMatrixFromRootFiles = kTRUE;  
   
 
-    // Boolean variable that controls wether to use the 
-    // TRAIN sample or not.
-    
-
-    Bool_t TrainParams   = kTRUE;  
-
-
-
-    // Variable that allows the user to skip the optimization on the 
-    // train sample. If optimization is skipped (value kTRUE), the
-    // previously optimized supercuts (stored in root file
-    // which is called OptSCParametersONOFFThetaRangeXXXXXmRad.root, and located 
-    // in the directory specified by variable PathForFiles) are used 
-    // on the train and/or the test sample.
-
-    // If value kFALSE, the cuts are optimized.
+    // Boolean variable that controls the supercuts 
+    // optimization using the training sample
     // The optimized cuts will be  written in root file
     // located in directory specified before. Name of 
     // the root files is created automatically.
-    
-    Bool_t SkipOptimization = kTRUE;
-  
-
-
-
+
+    Bool_t TrainParams   = kTRUE;  
+
+    // Variable that allows the user to skip the optimization on the 
+    // train sample. If optimization is skipped (value kTRUE), the
+    // previously optimized supercuts (stored in root file) are used 
+    // on the train sample.
+    
+    Bool_t SkipOptimization = kFALSE;
 
     // Boolean variable that allows the user to write the initial parameters 
     // into the root file that will be used to store the optimum cuts.
-    // If ApplyInitialParams = kTRUE , the initial 
-    // parameters are written into this root file, and they 
-    // will be applied to the data (TRAIN/TEST ) 
-    // IF NO OPTIMIZATION PROCEDURE IS PERFORMED.
-    
-    // If cuts are optimized (ie, variable SkipOptimization = kFALSE), 
-    // the cuts applied to the data are the optimized cuts.
+    // If fUseInitialSCParams = kTRUE , parameters are written.
+    // In this way, the initial SC parameters can be applied on the data (train/test) 
+    
+    // The initial parameters are ONLY written to the root file if 
+    // there is NO SC params optimization, i.e., if variable 
+    // fSkipOptimization = kTRUE;
   
-    // NOTE: be aware that, if ApplyInitialSCParams = kTRUE and 
-    // there was a root file with the optimized cuts 
+    // NOTE: be aware that, if there was a file with the optimized cuts 
     // (previously computed), it will be overwritten with the initial 
-    // SC parameters. 
-
-    Bool_t ApplyInitialSCParams = kTRUE;
-    
-
-
-    // Boolean variable that controls wether to use the 
-    // TEST sample or not.
-    
-
-    Bool_t TestParams = kFALSE;  
+    // SC parameters. This is something that I WILL HAVE TO CHANGE IN 
+    // future. Yet for the time being...
+
+    Bool_t UseInitialSCParams = kTRUE;  // was kFALSE 
+    
+
+
+    // Variable that decides whether the optimized cuts are used 
+    // in the test sample.
+
+    Bool_t TestParams = kTRUE;  
 
 
@@ -107 +103 @@
 
     // Fraction of ON events used for the training/testing
-    Double_t whichfractiontrain = 0.999;
-    Double_t whichfractiontest = 0.001;
+    Double_t whichfractiontrain = 0.5;
+    Double_t whichfractiontest = 0.5;
 
     // Fraction of OFF events used for the training/testing 
-    Double_t whichfractiontrainOFF = 0.999;
-    Double_t whichfractiontestOFF = 0.001;
+    Double_t whichfractiontrainOFF = 0.5;
+    Double_t whichfractiontestOFF = 0.5;
 
 
@@ -126 +122 @@
     // Alpha value (degrees) below which signal is expected
     
-    Double_t alphasig = 12; 
+    Double_t alphasig = 9; 
 
     // Definition of alpha bkg region (where no signal is expected)
@@ -135 +131 @@
     // Definition of the Size range
 
-    Double_t SizeLow = 800;
-    Double_t SizeUp = 1200;
-//    Double_t SizeUp = 1000000;
-    
+    Double_t SizeLow = 2000;
+//    Double_t SizeUp =  2000;
+    Double_t SizeUp = 1000000;
+    
+    Double_t LeakageMax = 0.04;
+    Double_t DistMax    = 1.5;
+    Double_t DistMin    = 0.1;
 
      // Definition of binning of alpha plots
@@ -167 +166 @@
     // in the computation of teh dynamical cuts that take place within
     // class MCT1SupercutsCalc
-    Bool_t NotUseTheta = kTRUE; // kTRUE renoves theta from the parameterization of cuts
+    //Bool_t NotUseTheta = kFALSE; // kTRUE removes theta from the parameterization of cuts
+    Bool_t NotUseTheta = kTRUE; // kTRUE removes theta from the parameterization of cuts
 
     // Boolean variable used to decide wether to use dynamical cuts or static cuts
@@ -192 +192 @@
     // Boolean variable used to decide wether initial parameters are 
     // read from ascii file or not. If kTRUE, parameters are retrieved 
-    // from ascii file. Otherwise, default parameters from MSupercuts 
-    // class are used.
+    // from ascii file. 
 
     Bool_t ReadInitParamsFromAsciiFile = kTRUE;
@@ -219 +218 @@
       // {"../InitialSCParametersSteps/StartingValuesForOptimizationMkn421DynCutsOnSize.txt"};
       // {"../InitialSCParametersSteps/StartingValuesForOptimizationMkn421DynCutsOnSizeAndDist.txt"};
-      {"mtemp/mmpi/asciifiles/OptimizedMkn421DynCutsGridWithSelected22pointsMay19.txt"};
+      //{"mtemp/mmpi/asciifiles/OptimizedMkn421DynCutsGridWithSelected22pointsMay19.txt"};
+      //{"/home/pcmagic16/mazin/mars/Mars260804/mtemp/mmpi/asciifiles/StartingValuesForSmallSizes1.txt"};
+//      {"/home/pcmagic16/mazin/mars/Mars260804/mtemp/mmpi/asciifiles/SmallSizeRFStartValues.txt"};
+      {"/home/pcmagic16/mazin/mars/Mars260804/mtemp/mmpi/asciifiles/OptimizedCrabFrom2000_a.txt"};
+      //{"mtemp/mmpi/asciifiles/VeryGoodDummyValuesWithConcCut.txt"};
 
 
@@ -244 +247 @@
 
     TArrayD CosThetaRangeVector(2);
-    CosThetaRangeVector[0] = 0.0;
+    CosThetaRangeVector[1] = 0.866;  //30
+    //CosThetaRangeVector[1] = 0.74;  //42
+    //CosThetaRangeVector[1] = 0.82;  //35
+    //CosThetaRangeVector[1] = 0;
     //CosThetaRangeVector[1] = 0.825;
     //CosThetaRangeVector[2] = 0.921; 
     //CosThetaRangeVector[3] = 0.961; 
-    CosThetaRangeVector[1] = 1.0;
+    CosThetaRangeVector[0] = 0.5;  // 60
+    //CosThetaRangeVector[0] = 1;
 
 
@@ -303 +310 @@
 
     FindSupercuts.SetSizeRange(SizeLow, SizeUp);
+    FindSupercuts.SetFilters(LeakageMax, DistMax, DistMin);
+//    FindSupercuts.SetFilters(0.001, 1.5, 0.2);
 
 
@@ -318 +327 @@
     FindSupercuts.SetTrainParameters(TrainParams);
     FindSupercuts.SetSkipOptimization(SkipOptimization);
-    FindSupercuts.SetUseInitialSCParams(ApplyInitialSCParams);
+    FindSupercuts.SetUseInitialSCParams(UseInitialSCParams);
 
     FindSupercuts.SetTestParameters(TestParams);