| 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 MPedestalCam;
|
|---|
| 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 Byte_t fgSaturationLimit; //! Default for fSaturationLimit (now set to: 254)
|
|---|
| 34 | static const TString fgNamePedestalCam; //! "MPedestalCam"
|
|---|
| 35 | static const TString fgNameSignalCam; //! "MExtractedSignalCam"
|
|---|
| 36 |
|
|---|
| 37 | Float_t fOffsetLoGain; // Offset of the low-gain signal w.r.t. the High-Gain slices
|
|---|
| 38 |
|
|---|
| 39 | MPedestalCam *fPedestals; //! Pedestals of all pixels in the camera
|
|---|
| 40 | MExtractedSignalCam *fSignals; //! Extracted signal of all pixels in the camera
|
|---|
| 41 |
|
|---|
| 42 | MRawEvtData *fRawEvt; //! Raw event data (time slices)
|
|---|
| 43 | MRawRunHeader *fRunHeader; //! RunHeader information
|
|---|
| 44 |
|
|---|
| 45 | Byte_t fHiGainFirst; // First FADC slice nr. to extract the High Gain signal
|
|---|
| 46 | Byte_t fHiGainLast; // Last FADC slice nr. to extract the High Gain signal
|
|---|
| 47 | Byte_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 | Byte_t fSaturationLimit; // Highest FADC slice value until being declared saturated
|
|---|
| 59 |
|
|---|
| 60 | TString fNamePedestalCam; // Name of the 'MPedestalCam' container
|
|---|
| 61 | TString fNameSignalCam; // Name of the 'MExtractedSignalCam' container
|
|---|
| 62 |
|
|---|
| 63 | // MExtractor
|
|---|
| 64 | virtual void FindSignalHiGain(Byte_t *firstused, Byte_t *lowgain, Float_t &sum, Byte_t &sat) const { }
|
|---|
| 65 | virtual void FindSignalLoGain(Byte_t *firstused, Float_t &sum, Byte_t &sat) const { }
|
|---|
| 66 |
|
|---|
| 67 | void SetResolutionPerPheHiGain(Float_t f) { fResolutionPerPheHiGain=f; }
|
|---|
| 68 | void SetResolutionPerPheLoGain(Float_t f) { fResolutionPerPheLoGain=f; }
|
|---|
| 69 |
|
|---|
| 70 | Int_t PreProcessStd(MParList *pList);
|
|---|
| 71 |
|
|---|
| 72 | // MTask
|
|---|
| 73 | Int_t PreProcess( MParList *pList );
|
|---|
| 74 | Bool_t ReInit ( MParList *pList );
|
|---|
| 75 | Int_t Process ();
|
|---|
| 76 | void StreamPrimitive(ostream &out) const;
|
|---|
| 77 | Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print);
|
|---|
| 78 |
|
|---|
| 79 |
|
|---|
| 80 | public:
|
|---|
| 81 | MExtractor(const char *name=NULL, const char *title=NULL);
|
|---|
| 82 |
|
|---|
| 83 | // getter
|
|---|
| 84 | Byte_t GetHiGainFirst() const { return fHiGainFirst; }
|
|---|
| 85 | Byte_t GetHiGainLast () const { return fHiGainLast ; }
|
|---|
| 86 | Byte_t GetLoGainFirst() const { return fLoGainFirst; }
|
|---|
| 87 | Byte_t GetLoGainLast () const { return fLoGainLast ; }
|
|---|
| 88 | Float_t GetNumHiGainSamples() const { return fNumHiGainSamples; }
|
|---|
| 89 | Float_t GetNumLoGainSamples() const { return fNumLoGainSamples; }
|
|---|
| 90 | Float_t GetOffsetLoGain() const { return fOffsetLoGain; }
|
|---|
| 91 | Float_t GetResolutionPerPheHiGain() const { return fResolutionPerPheHiGain; }
|
|---|
| 92 | Float_t GetResolutionPerPheLoGain() const { return fResolutionPerPheLoGain; }
|
|---|
| 93 |
|
|---|
| 94 | Bool_t IsNoiseCalculation () const { return fNoiseCalculation; }
|
|---|
| 95 |
|
|---|
| 96 | // Setter
|
|---|
| 97 | Float_t SetResolutionPerPheHiGain() const { return fResolutionPerPheHiGain; }
|
|---|
| 98 | Float_t SetResolutionPerPheLoGain() const { return fResolutionPerPheLoGain; }
|
|---|
| 99 |
|
|---|
| 100 | virtual void SetRange(Byte_t hifirst=0, Byte_t hilast=0, Byte_t lofirst=0, Byte_t lolast=0);
|
|---|
| 101 |
|
|---|
| 102 | void SetOffsetLoGain ( const Float_t f=fgOffsetLoGain ) { fOffsetLoGain = f; }
|
|---|
| 103 | void SetSaturationLimit ( const Byte_t lim=fgSaturationLimit ) { fSaturationLimit = lim; }
|
|---|
| 104 | void SetNamePedestalCam ( const char *name=fgNamePedestalCam.Data()) { fNamePedestalCam = name; }
|
|---|
| 105 | void SetNameSignalCam ( const char *name=fgNameSignalCam.Data() ) { fNameSignalCam = name; }
|
|---|
| 106 | void SetNoiseCalculation( const Bool_t b=kTRUE ) { fNoiseCalculation = b; }
|
|---|
| 107 |
|
|---|
| 108 | void SetPedestals (MPedestalCam *pedcam) { fPedestals = pedcam; }
|
|---|
| 109 |
|
|---|
| 110 | // TObject
|
|---|
| 111 | void Clear(Option_t *o="")
|
|---|
| 112 | {
|
|---|
| 113 | fHiGainFirst = fHiGainLast = fLoGainFirst = fLoGainLast = fHiLoLast = 0;
|
|---|
| 114 | }
|
|---|
| 115 |
|
|---|
| 116 | void Print(Option_t *o="") const; //*MENU*
|
|---|
| 117 |
|
|---|
| 118 | ClassDef(MExtractor, 6) // Signal Extractor Base Class
|
|---|
| 119 | };
|
|---|
| 120 |
|
|---|
| 121 | #endif
|
|---|
