Changeset 8173 for trunk/MagicSoft

Timestamp:

10/27/06 20:12:39 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mjobs/MJSpectrum.cc (modified) (3 diffs)
• mjobs/MJSpectrum.h (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -34 +34 @@
      - added reading and usage of "CalcHadronness" when processing 
        the MCs
+     - added a plot for the Cut efficiency
 
 

trunk/MagicSoft/Mars/mjobs/MJSpectrum.cc

@@ -36 +36 @@
 #include <TH1.h>
 #include <TH2.h>
+#include <TLine.h>
 #include <TFile.h>
 #include <TChain.h>
@@ -972 +973 @@
 
     return kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
+void MJSpectrum::DisplayCutEfficiency(const MHCollectionArea &area0, const MHCollectionArea &area1) const
+{
+    if (!fDisplay)
+        return;
+
+    const TH1D &trig = area0.GetHEnergy();
+    TH1D &cut = (TH1D&)*area1.GetHEnergy().Clone();
+
+    fDisplay->AddTab("CutEff");
+
+    gPad->SetBorderMode(0);
+    gPad->SetFrameBorderMode(0);
+    gPad->SetLogx();
+    gPad->SetGridx();
+    gPad->SetGridy();
+
+    cut.Divide(&trig);
+    cut.Scale(100);
+    cut.SetNameTitle("CutEff", "Background Supression: Cut efficiency (after star)");
+    cut.SetYTitle("\\eta [%]");
+    cut.SetDirectory(0);
+    cut.SetMinimum(0);
+    cut.SetMaximum(100);
+    cut.SetBit(kCanDelete);
+    cut.Draw();
+
+    TLine line;
+    line.SetLineColor(kBlue);
+    line.SetLineWidth(2);
+    line.SetLineStyle(kDashed);
+    line.DrawLine(cut.GetBinLowEdge(1), 50, cut.GetBinLowEdge(cut.GetNbinsX()+1), 50);
 }
 
@@ -1274 +1310 @@
     cout << "MC normalization factor:  " << scale << endl;
 
+    // Display cut efficiency
+    DisplayCutEfficiency(area0, area1);
+
     // Overlay normalized plots
     DisplaySize(plist, scale);

trunk/MagicSoft/Mars/mjobs/MJSpectrum.h

@@ -51 +51 @@
     TArrayD DisplaySpectrum(MHCollectionArea &area, TH1D &excess, MHEnergyEst &hest, Double_t ontime) const;
     Bool_t  DisplaySize(MParList &plist, Double_t scale) const;
+    void    DisplayCutEfficiency(const MHCollectionArea &area0, const MHCollectionArea &area1) const;
     Bool_t  PlotSame(MStatusArray &arr, MParList &plist, const char *name, const char *tab, const char *plot, Double_t scale) const;