Index: trunk/MagicSoft/Mars/Changelog
===================================================================
--- trunk/MagicSoft/Mars/Changelog	(revision 3847)
+++ trunk/MagicSoft/Mars/Changelog	(revision 3848)
@@ -18,4 +18,12 @@
 
                                                  -*-*- END OF LINE -*-*-
+
+ 2004/04/26: Abelardo Moralejo
+
+   * mhistmc/MHMcCollectionArea.[h,cc]
+     - Added fMCAreaRadius and setter to allow changing the radius of
+       the circular area in which MC events are produced. Useful for
+       instance for stereo configuration studies. By default it is 300 m
+       as was assumed before.
 
  2004/04/26: Markus Gaug
Index: trunk/MagicSoft/Mars/mhistmc/MHMcCollectionArea.cc
===================================================================
--- trunk/MagicSoft/Mars/mhistmc/MHMcCollectionArea.cc	(revision 3847)
+++ trunk/MagicSoft/Mars/mhistmc/MHMcCollectionArea.cc	(revision 3848)
@@ -55,4 +55,5 @@
 //
 MHMcCollectionArea::MHMcCollectionArea(const char *name, const char *title)
+  : fMCAreaRadius(300.)
 { 
     //   initialize the histogram for the distribution r vs E
@@ -256,7 +257,7 @@
     //
     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);
 
Index: trunk/MagicSoft/Mars/mhistmc/MHMcCollectionArea.h
===================================================================
--- trunk/MagicSoft/Mars/mhistmc/MHMcCollectionArea.h	(revision 3847)
+++ trunk/MagicSoft/Mars/mhistmc/MHMcCollectionArea.h	(revision 3848)
@@ -19,4 +19,8 @@
 
     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,4 +49,6 @@
     void Calc(const MHMcEfficiency &heff);
 
+    void SetMCAreaRadius(Float_t x) { fMCAreaRadius = x; } 
+
     ClassDef(MHMcCollectionArea, 1)  // Data Container to calculate Collection Area
 };