Changeset 2178 for trunk/MagicSoft/Mars/manalysis

Timestamp:

06/16/03 13:45:01 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MCerPhotCalc.cc (modified) (4 diffs)
• MCerPhotCalc.h (modified) (1 diff)
• MCerPhotEvt.cc (modified) (1 diff)
• MCerPhotEvt.h (modified) (3 diffs)
• MCurrents.h (modified) (3 diffs)
• MMcPedestalNSBAdd.cc (modified) (1 diff)
• MPedestalCam.cc (modified) (1 diff)
• MPedestalCam.h (modified) (3 diffs)

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc

@@ -32 +32 @@
 //   the array fWeight (default: all slices added up with weight 1).
 //
+//   The weights are rescaled, such that sum(weigths)=num slices
+//
 //  Input Containers:
 //   MRawRunHeader, MRawEvtData, MPedestalCam
@@ -112 +114 @@
 
     // Calculate quadratic sum of weights:
-    fSumQuadWeights = 0.;
-    for (Int_t i = 0; i < fWeight.GetSize(); i++)
+    fSumWeights = 0;
+    fSumQuadWeights = 0;
+    for (Int_t i=0; i<fWeight.GetSize(); i++)
+    {
+        fSumWeights += fWeight[i];
         fSumQuadWeights += fWeight[i]*fWeight[i];
+    }
 
     fSumQuadWeights = sqrt(fSumQuadWeights);
@@ -169 +175 @@
 
     MRawEvtPixelIter pixel(fRawEvt);
-
-    TArrayF binsignal(fWeight.GetSize());
 
     while (pixel.Next())
@@ -194 +198 @@
         // Calculate pixel signal unless it has all FADC slices empty:
         //
-        const Byte_t *ptr = pixel.GetHiGainSamples();
-
-        Float_t nphot = 0;
-        Float_t nphoterr = 0;
-
-        if (pixel.GetSumHiGainSamples()>0)
-        {
-            for (Int_t i=0; i<fWeight.GetSize(); i++)
-            {
-                binsignal[i] =  ptr[i] - mean;
-                nphot       +=  binsignal[i] * fWeight[i];
-            }
-            nphoterr = ped.GetSigma() * fSumQuadWeights;
-        }
+        Byte_t *ptr = pixel.GetHiGainSamples();
+
+        Float_t nphot = 0;
+        for(Int_t i=0; i<fWeight.GetSize(); i++)
+            nphot += ptr[i]*fWeight[i];
+
+        nphot -= mean*fSumWeights;
+
+        const Float_t nphoterr = ped.GetSigma()* fSumQuadWeights;
 
         fCerPhotEvt->AddPixel(pixid, nphot, nphoterr);

trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.h

@@ -33 +33 @@
 
     TArrayF         fWeight;  // Weights for adding up the ADC slices
+    Float_t         fSumWeights;
     Float_t         fSumQuadWeights;
 

trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.cc

@@ -384 +384 @@
     fNumPixels=fPixels->GetEntriesFast();
 }
+
+Bool_t MCerPhotEvt::GetPixelContent(Float_t &val, Int_t idx, Float_t ratio, Int_t type) const
+{
+    MCerPhotPix *pix = GetPixById(idx);
+    if (!pix || !pix->IsPixelUsed())
+        return kFALSE;
+
+    switch (type)
+    {
+    case 1:
+        val = pix->GetErrorPhot()*TMath::Sqrt(ratio);
+        return kTRUE;
+    case 2:
+        if (pix->GetErrorPhot()<=0)
+            return kFALSE;
+        val = pix->GetNumPhotons()*TMath::Sqrt(ratio)/pix->GetErrorPhot();
+        return kTRUE;
+    default:
+        val = pix->GetNumPhotons()*ratio;
+        return kTRUE;
+    }
+    return kTRUE;
+}

trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.h

@@ -4 +4 @@
 #ifndef ROOT_TClonesArray
 #include <TClonesArray.h>
+#endif
+#ifndef MARS_MCamEvent
+#include "MCamEvent.h"
 #endif
 #ifndef MARS_MCerPhotPix
@@ -12 +15 @@
 class MCerPhotPix;
 
-class MCerPhotEvt : public MParContainer
+class MCerPhotEvt : public MCamEvent
 {
 private:
@@ -61 +64 @@
     void Clear(Option_t *opt=NULL) { Reset(); }
 
+    Bool_t GetPixelContent(Float_t &val, Int_t idx, Float_t ratio=1, Int_t type=0) const;
+    void DrawPixelContent(Int_t num) const
+    {
+    }
+
     ClassDef(MCerPhotEvt, 1)    // class for an event containing cerenkov photons
 };

trunk/MagicSoft/Mars/manalysis/MCurrents.h

@@ -2 +2 @@
 #define MARS_MCurrents
 
-#ifndef MARS_MParContainer
-#include "MParContainer.h"
+#ifndef MARS_MCamEvent
+#include "MCamEvent.h"
 #endif
 
@@ -10 +10 @@
 #endif
 
-class MCurrents : public MParContainer
+#include <iostream.h>
+class MCurrents : public MCamEvent
 {
 private:
@@ -37 +38 @@
     void Print(Option_t *opt=NULL) const;
 
+    Bool_t GetPixelContent(Float_t &val, Int_t idx, Float_t ratio=1, Int_t type=0) const
+    {
+        val = (*this)[idx];
+        return val>0;
+    }
+    void DrawPixelContent(Int_t num) const
+    {
+    }
+
     ClassDef(MCurrents, 1) // Storage Container for the Currents (PRELIMINARY)
 };

trunk/MagicSoft/Mars/manalysis/MMcPedestalNSBAdd.cc

@@ -16 +16 @@
 !
 !
-!   Author(s): Oscar Blanch  11/2001 < mailto:blanch@ifae.es>
+!   Author(s): Oscar Blanch, 11/2001 <mailto:blanch@ifae.es>
 !
 !   Copyright: MAGIC Software Development, 2000-2001

trunk/MagicSoft/Mars/manalysis/MPedestalCam.cc

@@ -181 +181 @@
     return maxval;
 }
+
+Bool_t MPedestalCam::GetPixelContent(Float_t &val, Int_t idx, Float_t ratio=1, Int_t type=0) const
+{
+    val = (*this)[idx].GetMean()*ratio;
+    return kTRUE;
+}

trunk/MagicSoft/Mars/manalysis/MPedestalCam.h

@@ -2 +2 @@
 #define MARS_MPedestalCam
 
-#ifndef MARS_MParContainer
-#include "MParContainer.h"
+#ifndef MARS_MCamEvent
+#include "MCamEvent.h"
 #endif
 
@@ -11 +11 @@
 class MPedestalPix;
 
-class MPedestalCam : public MParContainer
+class MPedestalCam : public MCamEvent
 {
 private:
@@ -35 +35 @@
     void Print(Option_t *o="") const;
 
+    Bool_t GetPixelContent(Float_t &val, Int_t idx, Float_t ratio=1, Int_t type=0) const;
+    void DrawPixelContent(Int_t num) const
+    {
+    }
+
     ClassDef(MPedestalCam, 1)   // Storage Container for all pedestal information of the camera
 };