Changeset 1425

Timestamp:

07/22/02 17:32:22 (22 years ago)

Author:

bigongia

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• manalysis/MCerPhotCalc2.cc (modified) (8 diffs)
• manalysis/MCerPhotCalc2.h (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -1 +1 @@
                                                                   -*-*- END -*-*-
+
+ 2002/07/22: Abelardo Moralejo
+
+   * manalysis/MCerPhotCalc2.[h,cc]:
+     -added procedure SetWeights.
+
+   * macros/MagicHillas.C:
+     -added example on how to use MCerPhotCalc2
 
  2002/07/22: Thomas Bretz

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc2.cc

@@ -28 +28 @@
 //                                                                          //
 //   This is a task which calculates the number of photons from the FADC    //
-//   time slices. It weights the each slice according to the fraction       //
-//   of signal which, in average, falls within it. This average has been    //
-//   determined over a run (about 1000 triggers, 0 deg z.a. gammas).        // 
+//   time slices. It weights the each slice according to the numbers in     //
+//   the array fWeight. The default values (see end of this file) are       //
+//   the fraction of signal which, in average, falls within each slice.     //
+//   This averages have been determined over a run (about 1000 triggers,    //
+//   0 deg z.a. gammas).                                                    // 
 //                                                                          //
 //  Input Containers:                                                       //
-//   MRawEvtData, MPedestalCam                                             //
+//   MRawRunHeader, MRawEvtData, MPedestalCam                               //
 //                                                                          //
 //  Output Containers:                                                      //
@@ -65 +67 @@
 {
     fName  = name  ? name  : "MCerPhotCalc2";
-    fTitle = title ? title : "Task to calculate Cerenkov photons from raw data";
+    fTitle = title ? title : "Task to calculate pixel signal from raw data";
 
     AddToBranchList("MRawEvtData.fHiGainPixId");
@@ -71 +73 @@
     AddToBranchList("MRawEvtData.fHiGainFadcSamples");
     AddToBranchList("MRawEvtData.fLoGainFadcSamples");
-
-    fSumQuadWeights = fSumWeights = 0.;
-    for (Int_t i = 0; i < 15; i++)
-      {
-        fSumQuadWeights += fWeight[i]*fWeight[i];
-        fSumWeights += fWeight[i];
-      }
-    fSumQuadWeights = sqrt(fSumQuadWeights);
-
 }
 
@@ -120 +113 @@
         return kFALSE;
 
+    // Calculate quadratic sum of weights:
+    fSumQuadWeights = 0.;
+    for (Int_t i = 0; i < 15; i++)
+      fSumQuadWeights += fWeight[i]*fWeight[i];
+
+    fSumQuadWeights = sqrt(fSumQuadWeights);
+
     return kTRUE;
 }
@@ -138 +138 @@
     }
 
+    if (fRunHeader->GetNumSamplesHiGain() != 15)
+    {
+        *fLog << dbginf << "Number of FADC slices (" << fRunHeader->GetNumSamplesHiGain() <<") is different from expected (15)... aborting." << endl;
+        return kFALSE;
+    }
+
     if (runheader->GetRunType() != kRTMonteCarlo)
         return kTRUE;
@@ -165 +171 @@
     MRawEvtPixelIter pixel(fRawEvt);
 
-    Float_t ADCBins[15];
+    Float_t BinSignal[15];
 
     while (pixel.Next())
@@ -187 +193 @@
 
         Float_t nphot = 0.;
-        for(Int_t i = 0; i<15;i++)
+
+        for(Int_t i = 0; i<15; i++)
           {
-            ADCBins[i] = (Float_t) *(ptr+i);
-            nphot += ADCBins[i] * fWeight[i];
+            BinSignal[i] =  (Float_t) ptr[i] - mean;
+            nphot       +=  BinSignal[i] * fWeight[i];
           }
-
-        //
-        // We check that the pixel is not empty, if it is empty
-        // we won't substract the pedestal. Empty means that it has
-        // 0 signal in all the slices.
-        //
-        if (nphot!=0)
-          nphot -= mean*fSumWeights;
 
         Float_t nphoterr = ped.GetSigma()* fSumQuadWeights;
@@ -213 +212 @@
 }
 
-Float_t MCerPhotCalc2::fWeight[15] = {0, 0.0809835, 0.289593, 0.366926, 0.211665, 0.0508328, 0., 0., 0., 0., 0., 0., 0., 0., 0.};
+// Default weights: 
+Float_t MCerPhotCalc2::fWeight[] = {0, 0.0809835, 0.289593, 0.366926, 0.211665, 0.0508328, 0., 0., 0., 0., 0., 0., 0., 0., 0.};

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc2.h

@@ -25 +25 @@
     MRawEvtData    *fRawEvt;     // raw event data (time slices)
     MCerPhotEvt    *fCerPhotEvt; // Cerenkov Photon Event used for calculation
-    MRawRunHeader  *fRunHeader;  //  RunHeader information
+    MRawRunHeader  *fRunHeader;  // RunHeader information
  
     Bool_t          fEnableFix;  // fix for a bug in files from older camera versions (<=40)
 
-    Float_t         fSumQuadWeights, fSumWeights;
-    static Float_t  fWeight[15];
+    static Float_t  fWeight[15];  // Weights for adding up the ADC slices
+    Float_t         fSumQuadWeights;
 
 public:
@@ -37 +37 @@
     Bool_t PreProcess(MParList *pList);
     Bool_t Process();
+    Bool_t ReInit(MParList *pList);
 
-    Bool_t ReInit(MParList *pList);
+    void   SetWeights(Float_t *w) {memcpy(fWeight,w,15*sizeof(Float_t));}
 
     ClassDef(MCerPhotCalc2, 0)   // Task to calculate cerenkov photons from raw data
 };
-
+ 
 
 #endif