Ignore:
Timestamp:
06/29/04 11:10:51 (21 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/mpedestal
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mpedestal/MPedCalcFromData.cc

    r4356 r4362  
    6262    AddToBranchList("fHiGainFadcSamples");
    6363
    64     SetRange(fgLoGainFirst, fgLoGainLast);
     64    SetRange(fLoGainFirst, fLoGainLast);
    6565    Clear();
    6666}
     
    7676void MPedCalcFromData::SetLoRange(Byte_t lofirst, Byte_t lolast)
    7777{
    78     fgLoGainFirst = lofirst;
    79     fgLoGainLast = lolast;
     78    fLoGainFirst = lofirst;
     79    fLoGainLast = lolast;
    8080   
    8181    fWindowSizeLoGain = lolast - lofirst;
     
    141141      const UInt_t idx    = pixel.GetPixelId();
    142142     
    143       if ( (UInt_t)pixel.GetMaxHiGainSample() < HighGainThreshold ) {
     143      if ( (UInt_t)pixel.GetMaxHiGainSample() < fHiGainThreshold ) {
    144144         
    145145          fEvtCounter[idx]++;
     
    167167          fSumx2[idx]          += sqrsum;
    168168
    169           if (fEvtCounter[idx] == fgDump){
     169          if ((UInt_t)fEvtCounter[idx] == fDump){
    170170             
    171171              // Compute pedestals and rms from the sample
    172               const ULong_t n     = fWindowSizeLoGain*fgDump;
     172              const ULong_t n     = fWindowSizeLoGain*fDump;
    173173
    174174              const Float_t sum  = fSumx.At(idx);
     
    177177
    178178              // 1. Calculate the Variance of the sums:
    179               Float_t higainVar = (sum2-sum*sum/fgDump)/(fgDump-1.);
     179              Float_t higainVar = (sum2-sum*sum/fDump)/(fDump-1.);
    180180              // 2. Scale the variance to the number of slices:
    181181              higainVar /= (Float_t)(fWindowSizeLoGain);
     
    196196}
    197197
    198 
    199 Int_t MPedCalcFromData::PostProcess()
    200 {
    201   return kTRUE;
    202 }
  • trunk/MagicSoft/Mars/mpedestal/MPedCalcFromData.h

    r4356 r4362  
    2222#endif
    2323
    24 #ifndef MARS_MTask
    25 #include "MTask.h"
    26 #endif
    27 
    2824#ifndef MARS_MExtractor
    2925#include "MExtractor.h"
     
    3329{
    3430
    35     UInt_t fgDump; // Number for dumping.
    36     Byte_t fgLoGainFirst;      // First FADC slice Lo-Gain (currently set to: 3)
    37     Byte_t fgLoGainLast;       // Last FADC slice Lo-Gain (currently set to: 14)
     31    UInt_t fDump; // Number for dumping.
     32    Byte_t fLoGainFirst;      // First FADC slice Lo-Gain (currently set to: 3)
     33    Byte_t fLoGainLast;       // Last FADC slice Lo-Gain (currently set to: 14)
    3834    Byte_t fWindowSizeLoGain;             // Number of Lo Gain slices in window
    39     Byte_t HighGainThreshold;
     35    Byte_t fHiGainThreshold;
    4036
    4137    TArrayD fSumx;         // sum of values
     
    4844    Bool_t ReInit     ( MParList *pList );
    4945    Int_t  Process    ();
    50     Int_t  PostProcess();
    5146
    5247public:
     
    5550 
    5651    void Clear(const Option_t *o="");
    57     void SetDumpEvents(UInt_t dumpevents = 0){fgDump = dumpevents;}
    58     void SetHighGainThreshold(Byte_t Threshold = 0){HighGainThreshold = Threshold;}
     52    void SetDumpEvents(UInt_t dumpevents = 0){fDump = dumpevents;}
     53    void SetfHiGainThreshold(Byte_t Threshold = 0){fHiGainThreshold = Threshold;}
    5954    void SetLoRange(Byte_t lofirst=0, Byte_t lolast=0);
    6055
Note: See TracChangeset for help on using the changeset viewer.