| | 1 | #ifndef ALIITSONLINESDDBASE_H |
| | 2 | #define ALIITSONLINESDDBASE_H |
| | 3 | |
| | 4 | /////////////////////////////////////////////////////////////////// |
| | 5 | // // |
| | 6 | // Class used for SDD baseline and noise analysis // |
| | 7 | // Origin: F.Prino, Torino, prino@to.infn.it // |
| | 8 | // // |
| | 9 | /////////////////////////////////////////////////////////////////// |
| | 10 | #include "AliITSOnlineSDD.h" |
| | 11 | #include <TMath.h> |
| | 12 | |
| | 13 | class TH2F; |
| | 14 | class TGraph; |
| | 15 | class AliITSOnlineSDDBase : public AliITSOnlineSDD { |
| | 16 | |
| | 17 | public: |
| | 18 | AliITSOnlineSDDBase(); |
| | 19 | AliITSOnlineSDDBase(Int_t nddl, Int_t ncarlos, Int_t sid); |
| | 20 | virtual ~AliITSOnlineSDDBase(); |
| | 21 | void Reset(); |
| | 22 | void AddEvent(TH2F* hrawd); |
| | 23 | void ValidateAnodes(); |
| | 24 | |
| | 25 | void SetMinBaseline(Float_t bas=10.){fMinBaseline=bas;} |
| | 26 | void SetMaxBaseline(Float_t bas=150.){fMaxBaseline=bas;} |
| | 27 | void SetMinRawNoise(Float_t ns=0.5){fMinRawNoise=ns;} |
| | 28 | void SetMaxRawNoise(Float_t ns=9.){fMaxRawNoise=ns;} |
| | 29 | void SetNSigmaNoise(Float_t ns=4.){fNSigmaNoise=ns;} |
| | 30 | void SetGoldenBaselineValue(Float_t val=20.){fGoldenBaseline=val;} |
| | 31 | void SetZeroSuppThresholds(Float_t vall=2.2,Float_t valh=4.){ |
| | 32 | fLowThrFact=vall; |
| | 33 | fHighThrFact=valh; |
| | 34 | } |
| | 35 | |
| | 36 | Bool_t IsAnodeGood(Int_t iAnode)const{ return fGoodAnode[iAnode];} |
| | 37 | Float_t GetAnodeBaseline(Int_t iAnode) const{ |
| | 38 | if(fNEvents>0) return fSumBaseline[iAnode]/fNEvents; |
| | 39 | else return 0; |
| | 40 | } |
| | 41 | void GetMinAndMaxBaseline(Float_t &basMin, Float_t &basMax) const; |
| | 42 | Float_t GetMinimumBaseline() const; |
| | 43 | Float_t GetAnodeRawNoise(Int_t iAnode) const{ |
| | 44 | if(fNEvents>0) return TMath::Sqrt(fSumRawNoise[iAnode]/fNEvents-TMath::Power(GetAnodeBaseline(iAnode),2)); |
| | 45 | |
| | 46 | else return 0; |
| | 47 | } |
| | 48 | |
| | 49 | Int_t CountGoodAnodes() const{ |
| | 50 | Int_t nGdAn=0; |
| | 51 | for(Int_t ian=0;ian<fgkNAnodes;ian++) if(fGoodAnode[ian]) nGdAn++; |
| | 52 | return nGdAn; |
| | 53 | } |
| | 54 | Float_t CalcMeanRawNoise() const; |
| | 55 | Float_t GetAnodeCommonMode(Int_t iAnode) const{ |
| | 56 | if(fNEvents>0) return fSumCMN[iAnode]/fNEvents; |
| | 57 | else return 0; |
| | 58 | } |
| | 59 | Int_t GetNEvents() const {return fNEvents;} |
| | 60 | void WriteToASCII(); |
| | 61 | |
| | 62 | protected: |
| | 63 | |
| | 64 | private: |
| | 65 | static const Int_t fgkMaxCorr; // maximum baseline correction in AMBRA (=63) |
| | 66 | |
| | 67 | Int_t fNEvents; // number of events |
| | 68 | Bool_t fGoodAnode[fgkNAnodes]; // anode quality: good(1) - bad (0) |
| | 69 | Float_t fSumBaseline[fgkNAnodes]; // baseline summed over events |
| | 70 | Float_t fSumRawNoise[fgkNAnodes]; // noise summed over events |
| | 71 | Float_t fSumCMN[fgkNAnodes]; // common mode noise coeff. |
| | 72 | Float_t fMinBaseline; // Cut value for minimum baseline |
| | 73 | Float_t fMaxBaseline; // Cut value for maximum baseline |
| | 74 | Float_t fMinRawNoise; // Cut value for minimum noise |
| | 75 | Float_t fMaxRawNoise; // Cut value for maximum noise |
| | 76 | Float_t fNSigmaNoise; // Cut value for noise (n*sigma) |
| | 77 | Float_t fGoldenBaseline; // golden value for equalizing baselines |
| | 78 | Float_t fLowThrFact; // factor for low threshold |
| | 79 | Float_t fHighThrFact; // factor for high threshold |
| | 80 | |
| | 81 | ClassDef(AliITSOnlineSDDBase,2); |
| | 82 | }; |
| | 83 | #endif |
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental