Changeset 8709 for trunk/MagicSoft/Mars/mhflux

Timestamp:

08/25/07 16:30:26 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MHCollectionArea.cc (modified) (3 diffs)
• MHEnergyEst.cc (modified) (4 diffs)
• MMcSpectrumWeight.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/mhflux/MHCollectionArea.cc

@@ -99 +99 @@
 
     MBinning binsa, binse, binst;
-    binse.SetEdgesLog(15, 10, 1000000);
+    binse.SetEdgesLog(21, 6.3, 100000);
     binst.SetEdgesASin(67, -0.005, 0.665);
 
@@ -136 +136 @@
     MH::SetBinomialErrors(fHEnergy, *hsel, *hall);
 #endif
+
     fHEnergy.Scale(totalarea);
 
@@ -198 +199 @@
     }
 
-    // FIXME: Does this need some weighting with the number of produced events?
     if (runheader->GetImpactMax()>fMcAreaRadius*100)
     {

trunk/MagicSoft/Mars/mhflux/MHEnergyEst.cc

@@ -95 +95 @@
 
     MBinning binsi, binse, binst, binsr;
-    binse.SetEdgesLog(25, 10, 1000000);
+    binse.SetEdgesLog(21, 6.3, 100000);
     binst.SetEdgesASin(51, -0.005, 0.505);
-    binsi.SetEdges(10, 0, 400);
     binsr.SetEdges(75, -1.75, 1.75);
+
+    // Use the binning in impact to do efficiency studies
+    binsi.SetEdges(1, 0, 1000);
 
     SetBinning(&fHEnergy,     &binse, &binse, &binst);
@@ -513 +515 @@
     line.DrawLine(0,0,1,1);
 
+    line.SetLineColor(kBlue);
+    line.SetLineWidth(2);
+    line.SetLineStyle(kDashed);
+
     pad2->cd(2);
     h = MakePlot(fHResolution, "zy");
@@ -521 +527 @@
     h->SetMaximum(1.3);
 
+    line.DrawLine(h->GetXaxis()->GetXmin(), 0, h->GetXaxis()->GetXmax(), 0);
+
     pad2->cd(3);
     h = MakePlot(fHResolution, "zx");
@@ -528 +536 @@
     h->SetMinimum(-1.3);
     h->SetMaximum(1.3);
+
+    line.DrawLine(h->GetXaxis()->GetXmin(), 0, h->GetXaxis()->GetXmax(), 0);
 }
 

trunk/MagicSoft/Mars/mhflux/MMcSpectrumWeight.cc

@@ -493 +493 @@
     if (hasold)
     {
-        *fLog << " Old Spectrum: " << GetFormulaSpecOldX();
+        *fLog << " Old Spectrum:  " << GetFormulaSpecOldX();
         if (fEnergyMin>=0 && fEnergyMax>0)
             *fLog << "   (I=" << GetSpecOldIntegral() << ")";
@@ -500 +500 @@
     if (hasnew)
     {
-        *fLog << " New Spectrum: " << GetFormulaSpecNewX();
+        *fLog << " New Spectrum:  " << GetFormulaSpecNewX();
         if (fEnergyMin>=0 && fEnergyMax>0)
             *fLog << "   (I=" << GetSpecNewIntegral() << ")";