Changeset 1789 for trunk/MagicSoft/Mars

Timestamp:

02/24/03 12:26:32 (22 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mhist/MHMcRate.cc (modified) (2 diffs)
• mhist/MHMcRate.h (modified) (3 diffs)
• mmontecarlo/MMcTriggerRateCalc.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -6 +6 @@
       - Fixed mistake when analysing camera files containing all events 
         (including non-triggers): fShowers was always zero.
+      - Added reading of MMcTrigHeaders in ReInit
 
     * macros/trigrate.C
       - Added some explanations. Style improvements in function GetNSBEvents.
+
+    * mhist/MHMcRate.[h,cc]
+      - Added GetTriggerRate() and GetTriggerRateError()
+      - Added members fMultiplicity and fMeanThreshold, with their
+        corresponding Set and Get functions
+      - Added info on discriminator threshold and L1 trigger multiplicity
+        on the printout.
+
+    * mmc/MMcTrigHeader.hxx
+      - Added GetMultiplicity() and GetMeanThreshold()
 
  2003/02/22: Wolfgang Wittek

trunk/MagicSoft/Mars/mhist/MHMcRate.cc

@@ -61 +61 @@
     fTriggerRate= -1.0;      // Trigger rate in Hz
     fTriggerRateError= -1.0; // Estimated error for the trigger rate in Hz
+
+    fMultiplicity = -1;
+    fMeanThreshold = -1.;
+
 }
 
@@ -233 +237 @@
 {
     *fLog << all << "Incident rate " << fShowerRate << " Hz " << endl;
+    *fLog << "Multiplicity: " << fMultiplicity << ",  Discr. threshold: " << fMeanThreshold << endl;
     *fLog << "Trigger Rate " << fTriggerRate << " +- " << fTriggerRateError << " Hz" << endl;
 }

trunk/MagicSoft/Mars/mhist/MHMcRate.h

@@ -38 +38 @@
     Float_t fTriggerRateError;  // Estimated error for the trigger rate in Hz
 
+    Float_t fMeanThreshold;     // Mean discriminator threshold (mV) of trigger
+                                // pixels.
+
+    Short_t fMultiplicity;      // L1 trigger multiplicity.
+
+
     void Init(const char *name, const char *title);
 
@@ -65 +71 @@
     void SetEnergyMin(Float_t Energy) {fEnergyMin=Energy;}
 
+    void SetMultiplicity(Short_t nMul) {fMultiplicity = nMul;}
+    void SetMeanThreshold(Float_t thresh) {fMeanThreshold = thresh;}
+
     void UpdateBoundaries(Float_t energy, Float_t theta, Float_t phi, Float_t impact);
+
+    Float_t GetTriggerRate() {return fTriggerRate;}
+    Float_t GetTriggerRateError() {return fTriggerRateError;}
+
+    Short_t GetMultiplicity() {return fMultiplicity;}
+    Float_t  GetMeanThreshold() {return fMeanThreshold;}
 
     void CalcRate(Float_t trig, Float_t anal, Float_t simu);
@@ -78 +93 @@
 
 #endif 
+
+
+

trunk/MagicSoft/Mars/mmontecarlo/MMcTriggerRateCalc.cc

@@ -39 +39 @@
 #include "MMcTrig.hxx"
 #include "MMcRunHeader.hxx"
+#include "MMcTrigHeader.hxx"
 #include "MMcCorsikaRunHeader.h"
 
@@ -122 +123 @@
       "(5 degrees) for the original showers" << endl;
 
+
     for (UInt_t i=0; i<fNum; i++)
       {
@@ -136 +138 @@
         GetRate(i)->SetEnergyMin(fMcCorRunHeader->GetELowLim()/1000.);
         GetRate(i)->SetEnergyMax(fMcCorRunHeader->GetEUppLim()/1000.);
+
+        TString th("MMcTrigHeader;");
+        th += i+1;
+
+        MMcTrigHeader* trighead = (MMcTrigHeader*) pList->FindObject(th);
+        GetRate(i)->SetMeanThreshold(trighead->GetMeanThreshold());
+        GetRate(i)->SetMultiplicity(trighead->GetMultiplicity());
 
       }