Changeset 1886 for trunk/MagicSoft/Mars/manalysis

Timestamp:

04/01/03 17:49:19 (22 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MEnergyEstParam.cc (modified) (3 diffs)
• MEnergyEstParam.h (modified) (2 diffs)

trunk/MagicSoft/Mars/manalysis/MEnergyEstParam.cc

@@ -58 +58 @@
     fA[3] =  0.006;
     fA[4] =   58.3;
+
     /*
-     fA[0] =   39.667402; // [cm]
-     fA[1] =  143.081912; // [cm]
-     fA[2] = -395.501677; // [cm]
-     fA[3] =    0.006193;
-     fA[4] =    0;
-     */
+      fA[0] =   39.667402; // [cm]
+      fA[1] =  143.081912; // [cm]
+      fA[2] = -395.501677; // [cm]
+      fA[3] =    0.006193;
+      fA[4] =    0;
+    */
 
     fB.Set(7);
@@ -74 +75 @@
     fB[5] =  4.13e-8;   // [GeV]
     fB[6] =     0.05;
+
     /*
-     fB[0] =   45.037867; // [GeV]
-     fB[1] =    0.087222; // [GeV]
-     fB[2] =   -0.208065; // [GeV/cm]
-     fB[3] =   78.164942; // [GeV]
-     fB[4] = -159.361283; // [GeV]
-     fB[5] =   -0.130455; // [GeV/cm]
-     fB[6] =    0.051139;
-     */
+      fB[0] =   45.037867; // [GeV]
+      fB[1] =    0.087222; // [GeV]
+      fB[2] =   -0.208065; // [GeV/cm]
+      fB[3] =   78.164942; // [GeV]
+      fB[4] = -159.361283; // [GeV]
+      fB[5] =   -0.130455; // [GeV/cm]
+      fB[6] =    0.051139;
+    */
 }
 
@@ -340 +342 @@
         const Double_t dist = fMatrix ? GetVal(col++) : ((MHillasSrc*)NextH())->GetDist();
 
-
         /* DANIEL - MARCOS */
         const Double_t ir = i0 * (i1 + fA[1]*dist); // [cm]
         Double_t er = e0 * (e1 + e2*ir);       // [GeV]
 
-
-        /* MY PARAM */
+        /* THOMAS BRETZ */
         // if (width==0) return kCONTINUE;
         // const Double_t ir = i0 * (i1 + dist*(fA[1]/width + fA[4]/log10(size))); // [cm]

trunk/MagicSoft/Mars/manalysis/MEnergyEstParam.h

@@ -40 +40 @@
     Double_t GetVal(Int_t i) const;
 
+
 public:
     MEnergyEstParam(const char *hil="MHillas", const char *name=NULL, const char *title=NULL);
@@ -58 +59 @@
     void SetCoeffB(const TArrayD &arr);
 
+    Double_t GetCoeff(Int_t i) { return i<fA.GetSize()? fA[i] : fB[i-fA.GetSize()]; }
+
     void Print(Option_t *o=NULL);