#include "AliITSCalibration.h"
#include "AliITSresponseSDD.h"
+#include "AliITSsegmentationSDD.h"
#include "TArrayI.h"
class AliITSMapSDD;
+class AliITSDriftSpeedArraySDD;
class AliITSresponseSDD;
///////////////////////////////////////////////////////
// Response for SDD //
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]; }
- Float_t GetChannelGain(Int_t anode) const;
- 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 {
- Int_t chip=anode/fgkChannels;
- if(GetWing(anode)==1)chip-=fgkChips;
- return chip;
+ 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;
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;i<fgkChips*fgkWings;i++) fIsChipBad[i]=kTRUE;
+ }
+ virtual Bool_t IsBad() const { return fIsBad; }
+ void SetChipBad(Int_t nChip) {
+ fIsChipBad[nChip] = kTRUE;
+ }
+ virtual Bool_t IsChipBad(Int_t nChip) const {
+ return fIsChipBad[nChip];
+ }
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
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){
+ Bool_t IsBadChannel(Int_t anode) const{
if(GetChannelGain(anode)==0) return kTRUE;
else return kFALSE;
}
+ void SetUseCorrectionMaps(Bool_t useAnodeMap, Bool_t useDriftMap){
+ fUseACorrMap=useAnodeMap;
+ fUseTCorrMap=useDriftMap;
+ }
Float_t GetMapACell(Int_t i,Int_t j) const {
if(i<256) return fMapAW0->GetCellContent(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;
}
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();}
virtual void SetDynamicRange(Double_t p1) {((AliITSresponseSDD*)fResponse)->SetDynamicRange(p1);}
virtual Float_t GetDynamicRange() const {return ((AliITSresponseSDD*)fResponse)->DynamicRange();}
- 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 GetDriftSpeedAtAnode(Float_t nAnode) const;
-
+ 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 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 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];}
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)
+ 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,10) // SDD response
};
#endif