Changeset 1443 for trunk/MagicSoft/Mars/manalysis/MCerPhotCalc2.cc

Timestamp:

07/25/02 17:14:55 (22 years ago)

Author:

bigongia

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MCerPhotCalc2.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc2.cc

@@ -73 +73 @@
     AddToBranchList("MRawEvtData.fHiGainFadcSamples");
     AddToBranchList("MRawEvtData.fLoGainFadcSamples");
+
+    SetDefaultWeights();
 }
 
@@ -115 +117 @@
     // Calculate quadratic sum of weights:
     fSumQuadWeights = 0.;
-    for (Int_t i = 0; i < 15; i++)
+    for (Int_t i = 0; i < fWeight.GetSize(); i++)
       fSumQuadWeights += fWeight[i]*fWeight[i];
 
@@ -138 +140 @@
     }
 
-    if (fRunHeader->GetNumSamplesHiGain() != 15)
-    {
-        *fLog << dbginf << "Number of FADC slices (" << fRunHeader->GetNumSamplesHiGain() <<") is different from expected (15)... aborting." << endl;
+    if (fRunHeader->GetNumSamplesHiGain() != fWeight.GetSize())
+    {
+        *fLog << dbginf << "Number of FADC slices (" << fRunHeader->GetNumSamplesHiGain() <<") is different from assumed one (" << fWeight.GetSize() << ")... aborting." << endl;
         return kFALSE;
     }
@@ -171 +173 @@
     MRawEvtPixelIter pixel(fRawEvt);
 
-    Float_t BinSignal[15];
+    TArrayF BinSignal(fWeight.GetSize());
 
     while (pixel.Next())
@@ -194 +196 @@
         Float_t nphot = 0.;
 
-        for(Int_t i = 0; i<15; i++)
+        for(Int_t i = 0; i<fWeight.GetSize(); i++)
           {
             BinSignal[i] =  (Float_t) ptr[i] - mean;
@@ -212 +214 @@
 }
 
-// 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.};
+//
+// Set default values for the number of slices and weights:
+//
+
+void MCerPhotCalc2::SetDefaultWeights()
+{
+  const Float_t dummy[15] = {0, 0.0809835, 0.289593, 0.366926, 0.211665, 0.0508328, 0., 0., 0., 0., 0., 0., 0., 0., 0.};
+
+  fWeight.Set(15,dummy);
+  return;
+}
+
+
+
+