Changeset 7094 for trunk/MagicSoft/Mars/mhflux

Timestamp:

05/27/05 09:46:22 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MHCollectionArea.cc (modified) (8 diffs)
• MHCollectionArea.h (modified) (1 diff)
• MMcSpectrumWeight.cc (added)
• MMcSpectrumWeight.h (added)

trunk/MagicSoft/Mars/mhflux/MHCollectionArea.cc

@@ -31 +31 @@
 #include "MHCollectionArea.h" 
 
+#include <TLatex.h>
 #include <TCanvas.h>
 #include <TPaveStats.h>
@@ -81 +82 @@
     fHistAll.SetTitle("Number of events produced");
     fHistAll.SetXTitle("\\Theta [deg]");
-    fHistAll.SetYTitle("E [GeV]");
+    fHistAll.SetYTitle("E_{mc} [GeV]");
     fHistAll.SetDirectory(NULL);
     fHistAll.UseCurrentStyle();
 
     fHEnergy.SetName("CollEnergy");
-    fHEnergy.SetTitle("Collection Area vs Energy");
+    fHEnergy.SetTitle("Collection Area");
     fHEnergy.SetXTitle("E [GeV]");
-    fHEnergy.SetYTitle("A [m^{2}]");
+    fHEnergy.SetYTitle("A_{eff} [m^{2}]");
     fHEnergy.SetDirectory(NULL);
     fHEnergy.UseCurrentStyle();
@@ -119 +120 @@
     // read from run header, but it is not yet in!!
     //
-    const Float_t r1 = 0;
-    const Float_t r2 = fMcAreaRadius;
-
-    //*fLog << warn << endl << dbginf << "WARNING! I will assume a maximum impact parameter of " << fMCAreaRadius << " meters for the MC events. Check that this is the true one!" <<endl<<endl;
-    const Float_t total_area = TMath::Pi() * (r2*r2 - r1*r1);
+    const Float_t totalarea = GetCollectionAreaAbs();//TMath::Pi() * (r2*r2 - r1*r1);
 
     for (Int_t ix=1; ix<=nbinx; ix++)
@@ -141 +138 @@
             continue;
 
-        const Double_t eff = Ns/Na;
-
-        const Double_t efferr = TMath::Sqrt((1.-eff)*Ns)/Na;
+        const Double_t eff         = Ns/Na;
+        const Double_t efferr      = TMath::Sqrt((1.-eff)*Ns)/Na;
         
-        const Float_t col_area  =  eff * total_area;
-        const Float_t col_area_error =  efferr * total_area;
-
-        fHEnergy.SetBinContent(ix, col_area);
-        fHEnergy.SetBinError(ix, col_area_error);
+        const Float_t colarea      =  eff    * totalarea;
+        const Float_t colareaerror =  efferr * totalarea;
+
+        fHEnergy.SetBinContent(ix, colarea);
+        fHEnergy.SetBinError(ix,   colareaerror);
     }
 
@@ -275 +271 @@
 void MHCollectionArea::Paint(Option_t *option)
 {
+    if (TString(option)=="paint3")
+    {
+        /*
+        TH1 *h = dynamic_cast<TH1*>(gPad->FindObject("Efficiency"));
+        if (h)
+        {
+            const TString txt = Form("N/N_{0}=%.2f",
+                                 GetCollectionAreaEff(),
+                                 GetCollectionAreaAbs(), fMcAreaRadius);
+
+        TLatex text(0.31, 0.95, txt);
+        text.SetBit(TLatex::kTextNDC);
+        text.SetTextSize(0.04);
+        text.Paint();*/
+        return;
+    }
+    if (TString(option)=="paint4")
+    {
+        const TString txt = Form("A_{eff}=%.0fm^{2}  A_{abs}=%.0fm^{2}  r=%.0fm",
+                                 GetCollectionAreaEff(),
+                                 GetCollectionAreaAbs(), fMcAreaRadius);
+
+        TLatex text(0.31, 0.95, txt);
+        text.SetBit(TLatex::kTextNDC);
+        text.SetTextSize(0.04);
+        text.Paint();
+        return;
+    }
+
     TVirtualPad *pad;
 
@@ -405 +430 @@
         h->SetNameTitle("Efficiency", "Combined cut and trigger efficiency");
         h->SetDirectory(NULL);
+        AppendPad("paint3");
     }
     else
@@ -416 +442 @@
         gPad->SetGridy();
         fHEnergy.Draw();
+        AppendPad("paint4");
     }
     else
@@ -428 +455 @@
     //*fLog << energy << " " << theta << endl;
 
-    fHistSel.Fill(theta, energy);
+    fHistSel.Fill(theta, energy, weight);
 
     return kTRUE;

trunk/MagicSoft/Mars/mhflux/MHCollectionArea.h

@@ -56 +56 @@
     const TH1D &GetHEnergy() const { return fHEnergy; }
 
+    Double_t GetEntries() const { return fHistAll.Integral(); }
+    Double_t GetCollectionAreaEff() const { return fHEnergy.Integral(); }
+    Double_t GetCollectionAreaAbs() const { return TMath::Pi()*fMcAreaRadius*fMcAreaRadius; }
+
 /*
     void DrawAll(Option_t *option="");