Changeset 3848 for trunk/MagicSoft/Mars/mhistmc

Timestamp:

04/27/04 14:03:12 (20 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhistmc

Files:

• MHMcCollectionArea.cc (modified) (2 diffs)
• MHMcCollectionArea.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mhistmc/MHMcCollectionArea.cc

@@ -55 +55 @@
 //
 MHMcCollectionArea::MHMcCollectionArea(const char *name, const char *title)
+  : fMCAreaRadius(300.)
 { 
     //   initialize the histogram for the distribution r vs E
@@ -256 +257 @@
     //
     const Float_t r1 = 0;
-    const Float_t r2 = 300;
-
-    *fLog << warn << endl << dbginf << "WARNING! I will assume a maximum impact parameter of 300 meters for the MC events. Check that this is the true one!" <<endl<<endl;
+    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);
 

trunk/MagicSoft/Mars/mhistmc/MHMcCollectionArea.h

@@ -19 +19 @@
 
     TH1D *fHistCol; //  the collection area
+
+    Float_t fMCAreaRadius; // [m] Radius of circle within which the cores 
+                           // of Corsika events have been uniformly 
+                           // distributed.
 
     void Calc(TH2D &hsel, TH2D &hall);
@@ -45 +49 @@
     void Calc(const MHMcEfficiency &heff);
 
+    void SetMCAreaRadius(Float_t x) { fMCAreaRadius = x; } 
+
     ClassDef(MHMcCollectionArea, 1)  // Data Container to calculate Collection Area
 };