X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSCalibrationSDD.h;h=315476201735e693dc5213bf824e2295f45ee45f;hb=bc07eb4bbd90189149b21b698e3adbe290603d35;hp=19cf9d28c5e102b922197d88e50b945d3364586a;hpb=fcf95fc7a1a40644de272f60ebe2a9059a8ba9a7;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSCalibrationSDD.h b/ITS/AliITSCalibrationSDD.h index 19cf9d28c5e..31547620173 100644 --- a/ITS/AliITSCalibrationSDD.h +++ b/ITS/AliITSCalibrationSDD.h @@ -1,14 +1,18 @@ #ifndef ALIITSCALIBRATIONSDD_H #define ALIITSCALIBRATIONSDD_H -/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * +/* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. * * See cxx source for full Copyright notice */ +/* $Id$ */ #include "AliITSCalibration.h" -#include "AliITSresponseSDD.h" +#include "AliITSsegmentationSDD.h" +#include "TArrayI.h" + +class AliITSCorrMapSDD; +class AliITSDriftSpeedArraySDD; -class AliITSresponseSDD; /////////////////////////////////////////////////////// // Response for SDD // /////////////////////////////////////////////////////// @@ -20,33 +24,50 @@ class AliITSCalibrationSDD : public AliITSCalibration { // AliITSCalibrationSDD(); AliITSCalibrationSDD(const char *dataType); - virtual ~AliITSCalibrationSDD() {;} - virtual void SetNoiseParam(Double_t n, Double_t b){ - fNoise=n; fBaseline=b;} - - virtual void GetNoiseParam(Double_t &n, Double_t &b) const { - n=fNoise; b=fBaseline;} - - 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 ~AliITSCalibrationSDD(); + + virtual Float_t GetBaseline(Int_t anode) const {return fBaseline[anode];} + virtual void SetBaseline(Int_t anode,Double_t bas) {fBaseline[anode]=bas;} + 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 GiveCompressParam(Int_t *x) const; - void SetNoiseAfterElectronics(Double_t n=2.38){ + 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) // 2.36 for ALICE from beam test measurements 2001 - fNoiseAfterEl=n;} - Double_t GetNoiseAfterElectronics() const { + fNoiseAfterEl[anode]=n;} + Float_t GetNoiseAfterElectronics(Int_t anode) const { // Noise after electronics (ADC units) - return fNoiseAfterEl;} - void SetCompressParam(Int_t cp[8]); - void SetDeadChannels(Int_t nchips=0, Int_t nchannels=0); + return fNoiseAfterEl[anode];} + //void SetDeadChannels(Int_t nchips=0, Int_t nchannels=0); + void SetDeadChannels(Int_t ndead=0){fDeadChannels=ndead; fBadChannels.Set(ndead);} Int_t GetDeadChips() const { return fDeadChips; } Int_t GetDeadChannels() const { return fDeadChannels; } - Double_t Gain(Int_t wing,Int_t chip,Int_t ch)const - {return fGain[wing][chip][ch]; } + 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(); virtual void Print(ostream *os) const {AliITSCalibrationSDD::Print(os);} @@ -65,64 +86,119 @@ 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;i1) return kFALSE; + if(IsChipBad(nwing*4) && IsChipBad(nwing*4+1) && IsChipBad(nwing*4+2) && IsChipBad(nwing*4+3)) return kTRUE; + return kFALSE; + } + Int_t Wings()const{return fgkWings;}//Total number of SDD wings Int_t Chips() const{return fgkChips;} // Number of chips/module Int_t Channels() const{ return fgkChannels;}//Number of channels/chip + Int_t NOfAnodes() const {return fgkChannels*fgkChips*fgkWings;} + + virtual void SetBadChannel(Int_t i,Int_t anode); + Int_t GetBadChannel(Int_t i) const {return fBadChannels[i];} + Bool_t IsBadChannel(Int_t anode) const{ + if(GetChannelGain(anode)==0) return kTRUE; + else return kFALSE; + } + Float_t GetMapACell(Int_t i,Int_t j) const { + if(i<256) return fMapAW0->GetCellContent(i,j); + else return fMapAW1->GetCellContent(i-256,j); + } + virtual void SetMapA(Int_t wing,AliITSCorrMapSDD* mapA) { + if(wing==0) fMapAW0=mapA; + else fMapAW1=mapA; + } + Float_t GetMapTCell(Int_t i,Int_t j) const { + if(i<256) return fMapTW0->GetCellContent(i,j); + else return fMapTW1->GetCellContent(i-256,j); + } + virtual void SetMapT(Int_t wing,AliITSCorrMapSDD* mapT) { + if(wing==0) fMapTW0=mapT; + else fMapTW1=mapT; + } + + virtual void SetDriftSpeed(Int_t wing, AliITSDriftSpeedArraySDD* arr){ + if(wing==0) fDrSpeed0=arr; + else fDrSpeed1=arr; + } + + + 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 SetZeroSupp(Bool_t opt=kTRUE) {fZeroSupp=opt;} + virtual Bool_t GetZeroSupp() const {return fZeroSupp;} + + virtual void SetAMAt40MHz() {fAMAt20MHz=kFALSE;} + virtual void SetAMAt20MHz() {fAMAt20MHz=kTRUE;} + virtual Bool_t IsAMAt20MHz() const {return fAMAt20MHz;} + + void GetCorrections(Float_t z, Float_t x, Float_t &devz, Float_t &devx, AliITSsegmentationSDD* seg); + void GetShiftsForSimulation(Float_t z, Float_t x, Float_t &devz, Float_t &devx, AliITSsegmentationSDD* seg); + virtual Float_t GetThresholdAnode(Int_t anode, Double_t nsigma=2.2) const { + return nsigma*fNoiseAfterEl[anode];} - virtual void SetElectronics(Int_t p1=1) {((AliITSresponseSDD*)fResponse)->SetElectronics(p1);} - virtual Int_t GetElectronics() const {return ((AliITSresponseSDD*)fResponse)->Electronics();} - virtual void SetMaxAdc(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetMaxAdc(p1);} - virtual Double_t GetMaxAdc() const {return ((AliITSresponseSDD*)fResponse)->MaxAdc();} - virtual void SetChargeLoss(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetChargeLoss(p1);} - virtual Double_t GetChargeLoss() const {return ((AliITSresponseSDD*)fResponse)->ChargeLoss();} - virtual void SetDynamicRange(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetDynamicRange(p1);} - virtual Double_t GetDynamicRange() const {return ((AliITSresponseSDD*)fResponse)->DynamicRange();} - virtual void SetDriftSpeed(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetDriftSpeed(p1);} - virtual Double_t GetDriftSpeed() const {return ((AliITSresponseSDD*)fResponse)->DriftSpeed();} - 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);} - virtual Double_t GetNSigmaIntegration() const {return ((AliITSresponseSDD*)fResponse)->NSigmaIntegration();} - virtual void SetNLookUp(Int_t p1) {((AliITSresponseSDD*)fResponse)->SetNLookUp(p1);} - virtual Int_t GetGausNLookUp() const {return ((AliITSresponseSDD*)fResponse)->GausNLookUp();} - virtual Double_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 Double_t GetJitterError() const {return ((AliITSresponseSDD*)fResponse)->JitterError();} - virtual void SetDo10to8(Bool_t bitcomp=kTRUE) {((AliITSresponseSDD*)fResponse)->SetDo10to8(bitcomp);} protected: + + // these statis const should be move to AliITSsegmentationSDD static const Int_t fgkWings = 2; // Number of wings per module static const Int_t fgkChips = 4; // Number of chips/module static const Int_t fgkChannels = 64; // Number of channels/chip - static const Double_t fgkTemperatureDefault; // default for fT (Kelvin) - static const Double_t fgkNoiseDefault; // default for fNoise - static const Double_t fgkBaselineDefault; // default for fBaseline - static const Double_t fgkMinValDefault; // default for fMinVal + 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 fgkGainDefault; //default for gain + Bool_t fZeroSupp; // zero suppression + Bool_t fAMAt20MHz; // flag for Analog memory of Pascal at 20 MHz Int_t fDeadChips; // Number of dead chips Int_t fDeadChannels; // Number of dead channels - Double_t fGain[fgkWings][fgkChips][fgkChannels];//Array for channel gains - Int_t fCPar[8]; // Hardware compression parameters - Double_t fNoise; // Noise - Double_t fBaseline; // Baseline - Double_t fNoiseAfterEl; // Noise after electronics - Double_t fMinVal; // Min value used in 2D zero-suppression algo - - Bool_t fIsDead; // module is dead or alive ? - + 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 + + 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) + + + AliITSCorrMapSDD* fMapAW0; //! map of residuals on anode coord. wing 0 + AliITSCorrMapSDD* fMapAW1; //! map of residuals on anode coord. wing 1 + AliITSCorrMapSDD* fMapTW0; //! map of residuals on time coord. wing 0 + AliITSCorrMapSDD* 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,1) // SDD response + ClassDef(AliITSCalibrationSDD,17) }; #endif