#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 AliITSMapSDD;
+class AliITSDriftSpeedArraySDD;
class AliITSresponseSDD;
///////////////////////////////////////////////////////
// Response for SDD //
//
AliITSCalibrationSDD();
AliITSCalibrationSDD(const char *dataType);
- virtual ~AliITSCalibrationSDD() {;}
- virtual void SetNoiseParam(Double_t n, Double_t b){
- fNoise=n; fBaseline=b;}
+ virtual ~AliITSCalibrationSDD();
+ virtual void SetNoiseParam(Double_t /*n*/, Double_t /*b*/){
+ NotImplemented("SetNoiseParam");}
- virtual void GetNoiseParam(Double_t &n, Double_t &b) const {
- n=fNoise; b=fBaseline;}
+ virtual void GetNoiseParam(Double_t &/*n*/, Double_t &/*b*/) const {
+ NotImplemented("GetNoiseParam");}
+
+ 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 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) const;
+ virtual void GiveCompressParam(Int_t *x,Int_t ian) const;
- void SetNoiseAfterElectronics(Double_t n=2.38){
+ 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);}
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) 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-256,j);
+ }
+ virtual void SetMapA(Int_t wing,AliITSMapSDD* 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,AliITSMapSDD* 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 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 Float_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 Float_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 Float_t GetDynamicRange() const {return ((AliITSresponseSDD*)fResponse)->DynamicRange();}
+
+ 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{
+ 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);}
- virtual Double_t GetNSigmaIntegration() const {return ((AliITSresponseSDD*)fResponse)->NSigmaIntegration();}
+ virtual Float_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 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 Double_t GetJitterError() const {return ((AliITSresponseSDD*)fResponse)->JitterError();}
+ 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];}
+
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 fgkMinValDefault; // default for fMinVal
+ static const Float_t fgkGainDefault; //default for gain
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
+ Float_t fMinVal; // Min value used in 2D zero-suppression algo
+
+ 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)
+
+
+ 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,1) // SDD response
+ ClassDef(AliITSCalibrationSDD,10) // SDD response
};
#endif