Changeset 6232 for trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MFindSupercutsONOFFThetaLoop.cc

Timestamp:

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

Author:

mazin

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MFindSupercutsONOFFThetaLoop.cc (modified) (25 diffs)

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();
+
         
     }