X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSCalibrationSDD.h;h=d71f5d7238b88207329202654aba0c8fc40724f9;hb=f03c632cb03c3386432de0a20aa7e29ba2a928ba;hp=4144da2b4b7bedeed1114c13c456d8a30ead146a;hpb=a66a0eb6d2f5f7920015e3924e9caf067a558149;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSCalibrationSDD.h b/ITS/AliITSCalibrationSDD.h index 4144da2b4b7..d71f5d7238b 100644 --- a/ITS/AliITSCalibrationSDD.h +++ b/ITS/AliITSCalibrationSDD.h @@ -8,9 +8,11 @@ #include "AliITSCalibration.h" #include "AliITSresponseSDD.h" +#include "AliITSsegmentationSDD.h" #include "TArrayI.h" class AliITSMapSDD; +class AliITSDriftSpeedArraySDD; class AliITSresponseSDD; /////////////////////////////////////////////////////// // Response for SDD // @@ -35,12 +37,16 @@ class AliITSCalibrationSDD : public AliITSCalibration { virtual Float_t GetNoise(Int_t anode) const {return fNoise[anode];} virtual void SetNoise(Int_t anode, Double_t noise) {fNoise[anode]=noise;} - virtual void SetThresholds(Double_t mv, Double_t /* b */){ - // Min value used in 2D - could be used as a threshold setting - fMinVal = mv;} - virtual void Thresholds(Double_t & mv, Double_t & /* b */) const - {mv = fMinVal;} - virtual void GiveCompressParam(Int_t *x,Int_t ian) const; + virtual void SetThresholds(Double_t /* mv */, Double_t /* b */) { + NotImplemented("SetThresholds");} + virtual void Thresholds(Double_t & /* mv */, Double_t & /* b */) const { + NotImplemented("Thresholds");} + virtual void GiveCompressParam(Int_t *x) const; + + void SetZSLowThreshold(Int_t iWing, Int_t thr=25){fZSTL[iWing]=thr;} + void SetZSHighThreshold(Int_t iWing, Int_t thr=29){fZSTH[iWing]=thr;} + Int_t GetZSLowThreshold(Int_t iWing) const {return fZSTL[iWing];} + Int_t GetZSHighThreshold(Int_t iWing) const {return fZSTH[iWing];} void SetNoiseAfterElectronics(Int_t anode,Double_t n=2.38){ // Noise after electronics (ADC units) @@ -53,10 +59,24 @@ class AliITSCalibrationSDD : public AliITSCalibration { void SetDeadChannels(Int_t ndead=0){fDeadChannels=ndead; fBadChannels.Set(ndead);} Int_t GetDeadChips() const { return fDeadChips; } Int_t GetDeadChannels() const { return fDeadChannels; } - Float_t Gain(Int_t wing,Int_t chip,Int_t ch)const - {return fGain[wing][chip][ch]; } - virtual void SetGain(Double_t g,Int_t wing,Int_t chip, Int_t ch) - {fGain[wing][chip][ch]=g;} + Float_t GetChannelGain(Int_t anode) const {return fGain[anode];} + virtual void SetGain(Int_t anode,Double_t g){fGain[anode]=g;} + + + Int_t GetWing(Int_t anode) const{ + if(anode>=fgkChips*fgkChannels) return 1; + else return 0; + } + Int_t GetChipChannel(Int_t anode) const {return anode%fgkChannels;} + Int_t GetChip(Int_t anode) const {return anode/fgkChannels;} + Int_t GetAnodeNumber(Int_t iwing, Int_t ichip03, Int_t ichan) const { + if(iwing>=2 || ichip03>=4 || ichan>=64) return -1; + else return iwing*fgkChips*fgkChannels+ichip03*fgkChannels+ichan; + } + Int_t GetAnodeNumber(Int_t ichip07, Int_t ichan) const { + if(ichip07>=8 || ichan>=64) return -1; + else return ichip07*fgkChannels+ichan; + } void PrintGains() const; void Print(); @@ -76,19 +96,34 @@ class AliITSCalibrationSDD : public AliITSCalibration { virtual void SigmaSpread(Double_t & /* p1 */,Double_t & /* p2 */) const {NotImplemented("SigmaSpread");} - void SetDead() { fIsDead = kTRUE; }; - Bool_t IsDead() const { return fIsDead; }; + void SetBad() { + fIsBad = kTRUE; + for(Int_t i=0;iGetCellContent(i,j); - else return fMapAW1->GetCellContent(i,j); + else return fMapAW1->GetCellContent(i-256,j); } virtual void SetMapA(Int_t wing,AliITSMapSDD* mapA) { if(wing==0) fMapAW0=mapA; @@ -96,13 +131,17 @@ class AliITSCalibrationSDD : public AliITSCalibration { } Float_t GetMapTCell(Int_t i,Int_t j) const { if(i<256) return fMapTW0->GetCellContent(i,j); - else return fMapTW1->GetCellContent(i,j); + else return fMapTW1->GetCellContent(i-256,j); } virtual void SetMapT(Int_t wing,AliITSMapSDD* mapT) { if(wing==0) fMapTW0=mapT; else fMapTW1=mapT; } - static Int_t GetMapTimeNBin() {return fgkMapTimeNBin;} + + virtual void SetDriftSpeed(Int_t wing, AliITSDriftSpeedArraySDD* arr){ + if(wing==0) fDrSpeed0=arr; + else fDrSpeed1=arr; + } virtual void SetElectronics(Int_t p1=1) {((AliITSresponseSDD*)fResponse)->SetElectronics(p1);} virtual Int_t GetElectronics() const {return ((AliITSresponseSDD*)fResponse)->Electronics();} @@ -112,22 +151,17 @@ class AliITSCalibrationSDD : public AliITSCalibration { virtual Float_t GetChargeLoss() const {return ((AliITSresponseSDD*)fResponse)->ChargeLoss();} virtual void SetDynamicRange(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetDynamicRange(p1);} virtual Float_t GetDynamicRange() const {return ((AliITSresponseSDD*)fResponse)->DynamicRange();} - virtual void SetDriftSpeed(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetDriftSpeed(p1);} - virtual void SetDriftSpeedParam(Int_t iWing, Float_t* p); + virtual Float_t GetTimeOffset() const {return ((AliITSresponseSDD*)fResponse)->TimeOffset();} virtual Float_t GetADC2keV() const {return ((AliITSresponseSDD*)fResponse)->ADC2keV();} - virtual Float_t GetDriftSpeed() const {return ((AliITSresponseSDD*)fResponse)->DriftSpeed();} - virtual Float_t GetDriftSpeedAtAnode(Float_t nAnode) const { - if(nAnode<256){ - return fDriftVelParW0[0]+fDriftVelParW0[1]*nAnode+fDriftVelParW0[2]*nAnode*nAnode+fDriftVelParW0[3]*nAnode*nAnode*nAnode; - }else{ - nAnode-=256; - return fDriftVelParW1[0]+fDriftVelParW1[1]*nAnode+fDriftVelParW1[2]*nAnode*nAnode+fDriftVelParW1[3]*nAnode*nAnode*nAnode; - } + virtual Float_t GetDriftSpeedAtAnode(Float_t nAnode) const{ + if(fDrSpeed0==0 || fDrSpeed1==0) AliFatal("Drift speed not set\n"); + if(nAnode<256) return fDrSpeed0->GetDriftSpeed(0,nAnode); + else return fDrSpeed1->GetDriftSpeed(0,nAnode-256); } virtual void SetParamOptions(const char *opt1,const char *opt2) {((AliITSresponseSDD*)fResponse)->SetParamOptions(opt1,opt2);} virtual void GetParamOptions(char *opt1,char *opt2) const {((AliITSresponseSDD*)fResponse)->ParamOptions(opt1,opt2);} - virtual Bool_t Do10to8() const {return ((AliITSresponseSDD*)fResponse)->Do10to8();} + virtual void SetZeroSupp (const char *opt) {((AliITSresponseSDD*)fResponse)->SetZeroSupp(opt);} virtual const char *GetZeroSuppOption() const {return ((AliITSresponseSDD*)fResponse)->ZeroSuppOption();} virtual void SetNSigmaIntegration(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetNSigmaIntegration(p1);} @@ -135,14 +169,15 @@ class AliITSCalibrationSDD : public AliITSCalibration { virtual void SetNLookUp(Int_t p1) {((AliITSresponseSDD*)fResponse)->SetNLookUp(p1);} virtual Int_t GetGausNLookUp() const {return ((AliITSresponseSDD*)fResponse)->GausNLookUp();} virtual Float_t GetGausLookUp(Int_t i) const {return ((AliITSresponseSDD*)fResponse)->GausLookUp(i);} - virtual Int_t Convert8to10(Int_t signal) const {return ((AliITSresponseSDD*)fResponse)->Convert8to10(signal);} + virtual void SetJitterError(Double_t jitter=20) {((AliITSresponseSDD*)fResponse)->SetJitterError(jitter);} virtual Float_t GetJitterError() const {return ((AliITSresponseSDD*)fResponse)->JitterError();} - virtual Float_t GetDriftPath(Float_t time,Float_t /*anodecoord*/) const {return time*GetDriftSpeed();} + virtual Float_t GetDriftPath(Float_t time, Float_t xAnode) const {return time*GetDriftSpeedAtAnode(xAnode);} + void GetCorrections(Float_t z, Float_t x, Float_t &devz, Float_t &devx, AliITSsegmentationSDD* seg); virtual Float_t GetThresholdAnode(Int_t anode,Int_t nsigma=3) const { return nsigma*fNoiseAfterEl[anode];} - virtual void SetDo10to8(Bool_t bitcomp=kTRUE) {((AliITSresponseSDD*)fResponse)->SetDo10to8(bitcomp);} + protected: @@ -153,35 +188,39 @@ class AliITSCalibrationSDD : public AliITSCalibration { static const Float_t fgkTemperatureDefault; // default for fT (Kelvin) static const Float_t fgkNoiseDefault; // default for fNoise static const Float_t fgkBaselineDefault; // default for fBaseline - static const Float_t fgkMinValDefault; // default for fMinVal static const Float_t fgkGainDefault; //default for gain - static const Int_t fgkMapTimeNBin = 72; //map granularity along drift direction + Int_t fDeadChips; // Number of dead chips Int_t fDeadChannels; // Number of dead channels - Float_t fGain[fgkWings][fgkChips][fgkChannels];//Array for channel gains + Float_t fGain[fgkWings*fgkChips*fgkChannels]; //Array for channel gains Float_t fNoise[fgkWings*fgkChips*fgkChannels]; // Noise array Float_t fBaseline[fgkWings*fgkChips*fgkChannels]; // Baseline array Float_t fNoiseAfterEl[fgkWings*fgkChips*fgkChannels]; // Noise after electronics - Float_t fMinVal; // Min value used in 2D zero-suppression algo - Bool_t fIsDead; // module is dead or alive ? - TArrayI fBadChannels; //Array with bad anodes number (0-512) + Int_t fZSTL[2]; // Low threshold in 2D zero-suppression (2 hybrids) + Int_t fZSTH[2]; // High threshold in 2D zero-suppression (2 hybrids) + + Bool_t fIsBad; // module is dead or alive ? + Bool_t fIsChipBad[fgkWings*fgkChips]; // chip is dead or alive ? + TArrayI fBadChannels; //Array with bad anodes number (0-512) - Float_t fDriftVelParW0[4]; // Coeff. of pol3 fit to drift speed vs. anode (wing0) - Float_t fDriftVelParW1[4]; // Coeff. of pol3 fit to drift speed vs. anode (wing1) + Bool_t fUseACorrMap; // flag for the use of correction maps (anode) + Bool_t fUseTCorrMap; // flag for the use of correction maps (drift) + AliITSMapSDD* fMapAW0; //! map of residuals on anode coord. wing 0 AliITSMapSDD* fMapAW1; //! map of residuals on anode coord. wing 1 AliITSMapSDD* fMapTW0; //! map of residuals on time coord. wing 0 AliITSMapSDD* fMapTW1; //! map of residuals on time coord. wing 1 - + AliITSDriftSpeedArraySDD* fDrSpeed0; //! drift speed for wing 0 + AliITSDriftSpeedArraySDD* fDrSpeed1; //! drift speed for wing 1 private: AliITSCalibrationSDD(const AliITSCalibrationSDD &ob); // copy constructor AliITSCalibrationSDD& operator=(const AliITSCalibrationSDD & /* source */); // ass. op. - ClassDef(AliITSCalibrationSDD,7) // SDD response + ClassDef(AliITSCalibrationSDD,11) // SDD response }; #endif