| 1 | #ifndef MARS_MExtractor | 
|---|
| 2 | #define MARS_MExtractor | 
|---|
| 3 |  | 
|---|
| 4 | ///////////////////////////////////////////////////////////////////////////// | 
|---|
| 5 | //                                                                         // | 
|---|
| 6 | // MExtractor                                                              // | 
|---|
| 7 | //                                                                         // | 
|---|
| 8 | // Base class for the signal extractors                                    // | 
|---|
| 9 | //                                                                         // | 
|---|
| 10 | ///////////////////////////////////////////////////////////////////////////// | 
|---|
| 11 |  | 
|---|
| 12 | #ifndef MARS_MTask | 
|---|
| 13 | #include "MTask.h" | 
|---|
| 14 | #endif | 
|---|
| 15 |  | 
|---|
| 16 | class MRawEvtData; | 
|---|
| 17 | class MRawRunHeader; | 
|---|
| 18 |  | 
|---|
| 19 | class MPedestalSubtractedEvt; | 
|---|
| 20 | class MExtractedSignalCam; | 
|---|
| 21 |  | 
|---|
| 22 | class MExtractor : public MTask | 
|---|
| 23 | { | 
|---|
| 24 | private: | 
|---|
| 25 | static const Float_t fgOffsetLoGain;     //! Default for fOffsetLoGain (now set to 1.51 (= 5ns) | 
|---|
| 26 |  | 
|---|
| 27 | Bool_t  fNoiseCalculation;               //! Flag if extractor determines noise contribution from pedestal file. | 
|---|
| 28 |  | 
|---|
| 29 | Float_t fResolutionPerPheHiGain;         // Extractor-dependent charge resolution per phe for high-gain (see TDAS-0502). | 
|---|
| 30 | Float_t fResolutionPerPheLoGain;         // Extractor-dependent charge resolution per phe for low-gain  (see TDAS-0502). | 
|---|
| 31 |  | 
|---|
| 32 | protected: | 
|---|
| 33 | static const UInt_t  fgSaturationLimit;  //! Default for fSaturationLimit (now set to: 254) | 
|---|
| 34 | static const char   *fgNameSignalCam;    //! "MExtractedSignalCam" | 
|---|
| 35 |  | 
|---|
| 36 | Float_t fOffsetLoGain;                   // Offset of the low-gain signal w.r.t. the High-Gain slices | 
|---|
| 37 |  | 
|---|
| 38 | MExtractedSignalCam *fSignals;           //! Extracted signal of all pixels in the camera | 
|---|
| 39 |  | 
|---|
| 40 | MRawEvtData         *fRawEvt;            //! Raw event data (time slices) | 
|---|
| 41 | MRawRunHeader       *fRunHeader;         //! RunHeader information | 
|---|
| 42 |  | 
|---|
| 43 | MPedestalSubtractedEvt *fSignal;         //! | 
|---|
| 44 |  | 
|---|
| 45 | UShort_t fHiGainFirst;                   // First FADC slice nr. to extract the High Gain signal | 
|---|
| 46 | UShort_t fHiGainLast;                    // Last FADC slice nr. to extract the High Gain signal | 
|---|
| 47 | Int_t    fLoGainFirst;                   // First FADC slice nr. to extract the Low Gain signal | 
|---|
| 48 | Byte_t   fLoGainLast;                    // Last FADC slice nr. to extract the Low Gain signal | 
|---|
| 49 |  | 
|---|
| 50 | //  Byte_t   fHiLoLast;                      // Number of slices in fLoGainSamples counted for the High-Gain signal | 
|---|
| 51 |  | 
|---|
| 52 | Float_t  fNumHiGainSamples;              // Number High Gain FADC slices used to extract the signal | 
|---|
| 53 | Float_t  fNumLoGainSamples;              // Number Low  Gain FADC slices used to extract the signal | 
|---|
| 54 |  | 
|---|
| 55 | Float_t  fSqrtHiGainSamples;             // Sqrt. nr. High Gain FADC slices used to extract the signal | 
|---|
| 56 | Float_t  fSqrtLoGainSamples;             // Sqrt. nr. Low  Gain FADC slices used to extract the signal | 
|---|
| 57 |  | 
|---|
| 58 | UInt_t   fSaturationLimit;               // Highest FADC slice value until being declared saturated | 
|---|
| 59 |  | 
|---|
| 60 | TString  fNameSignalCam;                 // Name of the 'MExtractedSignalCam' container | 
|---|
| 61 |  | 
|---|
| 62 | // MExtractor | 
|---|
| 63 | virtual void FindSignalHiGain(Byte_t *firstused, Byte_t *lowgain, Float_t &sum, Byte_t &sat) const { } | 
|---|
| 64 | virtual void FindSignalLoGain(Byte_t *firstused, Float_t &sum, Byte_t &sat) const { } | 
|---|
| 65 |  | 
|---|
| 66 | void SetResolutionPerPheHiGain(Float_t f) { fResolutionPerPheHiGain=f; } | 
|---|
| 67 | void SetResolutionPerPheLoGain(Float_t f) { fResolutionPerPheLoGain=f; } | 
|---|
| 68 |  | 
|---|
| 69 | Int_t   PreProcessStd(MParList *pList); | 
|---|
| 70 |  | 
|---|
| 71 | // MTask | 
|---|
| 72 | Int_t   PreProcess( MParList *pList ); | 
|---|
| 73 | Bool_t  ReInit    ( MParList *pList ); | 
|---|
| 74 | Int_t   Process   (); | 
|---|
| 75 | void    StreamPrimitive(ostream &out) const; | 
|---|
| 76 | Int_t   ReadEnv(const TEnv &env, TString prefix, Bool_t print); | 
|---|
| 77 |  | 
|---|
| 78 |  | 
|---|
| 79 | public: | 
|---|
| 80 | MExtractor(const char *name=NULL, const char *title=NULL); | 
|---|
| 81 |  | 
|---|
| 82 | // getter | 
|---|
| 83 | UShort_t GetHiGainFirst()      const { return fHiGainFirst;      } | 
|---|
| 84 | UShort_t GetHiGainLast ()      const { return fHiGainLast ;      } | 
|---|
| 85 | Int_t    GetLoGainFirst()      const { return fLoGainFirst;      } | 
|---|
| 86 | Byte_t   GetLoGainLast ()      const { return fLoGainLast ;      } | 
|---|
| 87 | Float_t  GetNumHiGainSamples() const { return fNumHiGainSamples; } | 
|---|
| 88 | Float_t  GetNumLoGainSamples() const { return fNumLoGainSamples; } | 
|---|
| 89 | Float_t  GetOffsetLoGain()     const { return fOffsetLoGain;     } | 
|---|
| 90 | Float_t  GetResolutionPerPheHiGain() const { return fResolutionPerPheHiGain; } | 
|---|
| 91 | Float_t  GetResolutionPerPheLoGain() const { return fResolutionPerPheLoGain; } | 
|---|
| 92 | UInt_t   GetSaturationLimit() const { return fSaturationLimit; } | 
|---|
| 93 |  | 
|---|
| 94 | Bool_t  HasLoGain() const { return fLoGainFirst>0 || fLoGainLast>0; } | 
|---|
| 95 | Bool_t  IsNoiseCalculation () const { return fNoiseCalculation; } | 
|---|
| 96 |  | 
|---|
| 97 | // Setter | 
|---|
| 98 | Float_t SetResolutionPerPheHiGain() const { return fResolutionPerPheHiGain; } | 
|---|
| 99 | Float_t SetResolutionPerPheLoGain() const { return fResolutionPerPheLoGain; } | 
|---|
| 100 |  | 
|---|
| 101 | virtual void SetRange(UShort_t hifirst=0, UShort_t hilast=0, Int_t lofirst=0, Byte_t lolast=0); | 
|---|
| 102 |  | 
|---|
| 103 | void SetOffsetLoGain    ( const Float_t  f=fgOffsetLoGain          ) { fOffsetLoGain     = f;    } | 
|---|
| 104 | void SetSaturationLimit ( const UInt_t lim=fgSaturationLimit       ) { fSaturationLimit  = lim;  } | 
|---|
| 105 | void SetNameSignalCam   ( const char *name=fgNameSignalCam ) { fNameSignalCam    = name; } | 
|---|
| 106 | void SetNoiseCalculation( const Bool_t   b=kTRUE                   ) { fNoiseCalculation = b;    } | 
|---|
| 107 |  | 
|---|
| 108 | // TObject | 
|---|
| 109 | void Clear(Option_t *o="") | 
|---|
| 110 | { | 
|---|
| 111 | fHiGainFirst = fHiGainLast = fLoGainFirst = fLoGainLast = /*fHiLoLast =*/ 0; | 
|---|
| 112 | } | 
|---|
| 113 |  | 
|---|
| 114 | void Print(Option_t *o="") const; //*MENU* | 
|---|
| 115 |  | 
|---|
| 116 | ClassDef(MExtractor, 9) // Signal Extractor Base Class | 
|---|
| 117 | }; | 
|---|
| 118 |  | 
|---|
| 119 | #endif | 
|---|