Fixes for Coverity warnings (M. van Leeuwen)
[u/mrichter/AliRoot.git] / ITS / AliITSCalibrationSDD.h
index 19cf9d2..971c45b 100644 (file)
@@ -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,112 @@ 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;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
+    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