Changeset 2268 for trunk/MagicSoft/Mars/mhistmc/MHMcCT1CollectionArea.cc

Timestamp:

07/07/03 16:22:51 (21 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhistmc/MHMcCT1CollectionArea.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mhistmc/MHMcCT1CollectionArea.cc

@@ -66 +66 @@
   fTitle = title ? title : "Collection Area vs. log10 Energy";
 
+  fEaxis = kLog10Energy;
+
   fHistAll = new TH2D;
   fHistSel = new TH2D;
@@ -93 +95 @@
   fHistCol->SetYTitle("theta [deg]");
   fHistCol->SetZTitle("A [m^{2}]");
+
 }
 
@@ -130 +133 @@
     fHistCol->Sumw2();
 
+    if (fEaxis == kEnergy)
+      {
+        fTitle = "Collection Area vs. Energy";
+        fHistCol->SetTitle(fTitle);
+        fHistAll->SetTitle("All showers - Theta vs Energy distribution");
+        fHistSel->SetTitle("Selected showers - Theta vs Energy distribution");
+        fHistCol->SetXTitle("E [GeV]");
+        fHistAll->SetXTitle("E [GeV]");
+        fHistSel->SetXTitle("E [GeV]");
+      }
+
     return kTRUE;
 }
@@ -142 +156 @@
   MMcEvt &mcevt = *(MMcEvt*)par;
 
-  fHistSel->Fill(log10(mcevt.GetEnergy()), kRad2Deg*mcevt.GetTelescopeTheta(), w);
+  if (fEaxis == kLog10Energy)
+    fHistSel->Fill(log10(mcevt.GetEnergy()), kRad2Deg*mcevt.GetTelescopeTheta(), w);
+  else
+    fHistSel->Fill(mcevt.GetEnergy(), kRad2Deg*mcevt.GetTelescopeTheta(), w);
+
   return kTRUE;
 }
@@ -193 +211 @@
 
   c.cd(1);
+  if (fEaxis == kEnergy)
+    gPad->SetLogx();
   fHistCol->SetDirectory(NULL);
   fHistCol->DrawCopy(option);
 
   c.cd(2);
+  if (fEaxis == kEnergy)
+    gPad->SetLogx();
   fHistSel->SetDirectory(NULL);
   fHistSel->DrawCopy(option);
 
   c.cd(3);
+  if (fEaxis == kEnergy)
+    gPad->SetLogx();
   fHistAll->SetDirectory(NULL);
   fHistAll->DrawCopy(option);
@@ -435 +459 @@
       for (Int_t i=1; i <= fHistAll->GetNbinsX(); i++)
         {
-          Float_t e1 = pow(10.,fHistAll->GetXaxis()->GetBinLowEdge(i));
-          Float_t e2 = pow(10.,fHistAll->GetXaxis()->GetBinLowEdge(i+1));
+          Float_t e1;
+          Float_t e2;
+
+          if (fEaxis == kLog10Energy)
+            {
+              e1 = pow(10.,fHistAll->GetXaxis()->GetBinLowEdge(i));
+              e2 = pow(10.,fHistAll->GetXaxis()->GetBinLowEdge(i+1));
+            }
+          else
+            {
+              e1 = fHistAll->GetXaxis()->GetBinLowEdge(i);
+              e2 = fHistAll->GetXaxis()->GetBinLowEdge(i+1);
+            }
 
           Float_t events = 0.;