]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/AliITSCalibrationSSD.h
SSD calibration
[u/mrichter/AliRoot.git] / ITS / AliITSCalibrationSSD.h
index b83584e8d879f81e8aae00906bf60b80fc83aaa1..e4ddf2899c38bc80ddf9a0ccc71d3e65d263ca8d 100644 (file)
@@ -3,26 +3,67 @@
  
 #include "AliITSCalibration.h"
 #include "AliITSresponseSSD.h"
+#include "TArrayD.h"
+#include "TArrayI.h"
 
 //////////////////////////////////////////////
 // Response class for SSD                   //
 //                                          //
 //////////////////////////////////////////////
+
+
 class AliITSCalibrationSSD : public AliITSCalibration {
 
  public:
     AliITSCalibrationSSD();
     AliITSCalibrationSSD(const char *dataType);
     virtual ~AliITSCalibrationSSD();
+
     virtual  void   SetNoiseParam(Double_t np, Double_t nn) {
-       // set noise par
-       fNoiseP=np; fNoiseN=nn;
+       // set noise par
+       fNoiseP=np; fNoiseN=nn;
     }
+
     virtual void    GetNoiseParam(Double_t &np, Double_t &nn) const {
-       // get noise par
-       np=fNoiseP; nn=fNoiseN;
+       // get noise par
+       np=fNoiseP; nn=fNoiseN;
     }
 
+    // EF
+    void SetNNoiseP(Int_t n) { fNoisP.Set(n); }
+    void AddNoiseP(Int_t c, Double_t n) { fNoisP.AddAt(n,c);}       
+    TArrayD GetNoiseP() const {return fNoisP; }
+    Double_t GetNoiseP(Int_t n) {return fNoisP.At(n); }
+    void SetNNoiseN(Int_t n) { fNoisN.Set(n); }
+    void AddNoiseN(Int_t c, Double_t n) { fNoisN.AddAt(n,c);}
+    TArrayD GetNoiseN() const {return fNoisN; }
+    Double_t GetNoiseN(Int_t n) {return fNoisN.At(n); }
+
+    void SetNGainP(Int_t n) { fGainP.Set(n); }
+    void AddGainP(Int_t c, Double_t n) { fGainP.AddAt(n,c);}       
+    TArrayD GetGainP() const {return fGainP; }
+    Double_t GetGainP(Int_t n) {return fGainP.At(n); }
+    void SetNGainN(Int_t n) { fGainN.Set(n); }
+    void AddGainN(Int_t c, Double_t n) { fGainN.AddAt(n,c);}
+    TArrayD GetGainN() const {return fGainN; }
+    Double_t GetGainN(Int_t n) {return fGainN.At(n); }
+
+    void SetNoisePThreshold(Int_t threshold) { fNoisePThreshold = threshold;}
+    void AddNoisyPChannel(Int_t c, Int_t n) { fNoisyPChannelsList.AddAt(n,c);}
+    TArrayI GetNoisyPChannelsList() const {return fNoisyPChannelsList; }
+    void SetNoiseNThreshold(Int_t threshold) { fNoiseNThreshold = threshold;}
+    void AddNoisyNChannel(Int_t c, Int_t n) { fNoisyNChannelsList.AddAt(n,c);}
+    TArrayI GetNoisyNChannelsList() const {return fNoisyNChannelsList; }
+
+    void SetNDeadPChannelsList(Int_t n) { fDeadPChannelsList.Set(n); }
+    void AddDeadPChannel(Int_t c, Int_t n) { fDeadPChannelsList.AddAt(n,c);}
+    TArrayI GetDeadPChannelsList() const {return fDeadPChannelsList; }
+    void SetNDeadNChannelsList(Int_t n) { fDeadNChannelsList.Set(n); }
+    void AddDeadNChannel(Int_t c, Int_t n) { fDeadNChannelsList.AddAt(n,c);}
+    TArrayI GetDeadNChannelsList() const {return fDeadNChannelsList; }
+    //
+
+
     virtual  void   SetNDetParam(Int_t npar) {
        // set number of param
        fNPar=npar;
@@ -37,13 +78,14 @@ class AliITSCalibrationSSD : public AliITSCalibration {
     virtual void    GetDetParam(Double_t *dpar) const; 
 
     virtual void    SetSigmaSpread(Double_t p1, Double_t p2) {
-       // Set sigmas of the charge spread function: Pside-Nside
-       // square of (microns)
-       fSigmaP=p1; fSigmaN=p2;
+       // Set sigmas of the charge spread function: Pside-Nside
+       // square of (microns)
+       fSigmaP=p1; fSigmaN=p2;
     }
+
     virtual void    SigmaSpread(Double_t &sP, Double_t &sN) const {
-       // Get sigmas for the charge spread 
-       sP=fSigmaP; sN=fSigmaN;
+      // Get sigmas for the charge spread 
+      sP=fSigmaP; sN=fSigmaN;
     }
 
     virtual void   SetThresholds(Double_t /* a */, Double_t /* b */)
@@ -53,25 +95,50 @@ class AliITSCalibrationSSD : public AliITSCalibration {
   
     virtual void SetParamOptions(const char *opt1, const char *opt2) {((AliITSresponseSSD*)fResponse)->SetParamOptions(opt1,opt2);}
     virtual void GetParamOptions(char *opt1,char *opt2) const {((AliITSresponseSSD*)fResponse)->ParamOptions(opt1,opt2);}
-    virtual void SetADCpereV(Double_t a=50./30000.0) {((AliITSresponseSSD*)fResponse)->SetADCpereV(a);}
+    virtual void SetADCpereV(Double_t a=200./30000.0) {((AliITSresponseSSD*)fResponse)->SetADCpereV(a);}
     virtual Double_t GetDEvToADC(Double_t eV) const {return ((AliITSresponseSSD*)fResponse)->DEvToADC(eV);}
     virtual Int_t IEvToADC(Double_t eV) const {return ((AliITSresponseSSD*)fResponse)->IEvToADC(eV);} 
 
- protected:
+  virtual Double_t GetCouplingPR() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingPR();}
+  virtual Double_t GetCouplingPL() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingPL();}
+  virtual Double_t GetCouplingNR() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingNR();}
+  virtual Double_t GetCouplingNL() const {return ((AliITSresponseSSD*)fResponse)->GetCouplingNL();}
+  virtual void SetCouplings(Double_t pr, Double_t pl, Double_t nr, Double_t nl) 
+    { ((AliITSresponseSSD*)fResponse)->SetCouplings(pr,pl,nr,nl);}
+
+  virtual Int_t GetZSThreshold() const {return ((AliITSresponseSSD*)fResponse)->GetZSThreshold();}
+  virtual void SetZSThreshold(Int_t zsth) 
+    { ((AliITSresponseSSD*)fResponse)->SetZSThreshold(zsth);}
+
+protected:
     static const Double_t fgkNoiseNDefault; // default for fNoiseN
     static const Double_t fgkNoisePDefault; // default for fNoiseP
     static const Int_t fgkNParDefault; // default for fNPar
     static const Double_t fgkSigmaPDefault; //default for fSigmaP
     static const Double_t fgkSigmaNDefault; //default for fSigmaP
+
     Int_t   fNPar;            // Number of detector param 
     Double_t *fDetPar;         //[fNPar] Array of parameters
 
     Double_t fNoiseP;          // Noise on Pside
     Double_t fNoiseN;          // Noise on Nside
-
     Double_t fSigmaP;          // Sigma charge spread on Pside
     Double_t fSigmaN;          // Sigma charge spread on Nside
     
+    TArrayD fGainP;           // Gain for P side channels
+    TArrayD fGainN;           // Gain for N side channels
+
+    TArrayD fNoisP;           // Noise for P side channels
+    TArrayD fNoisN;           // Noise for N side channels
+
+    Double_t fNoisePThreshold;     // need to decide if channel is noisy  
+    TArrayI  fNoisyPChannelsList; // list of P side noisy channels
+    Double_t fNoiseNThreshold;     // need to decide if channel is noisy  
+    TArrayI  fNoisyNChannelsList; // list of N side noisy channels
+
+    TArrayI  fDeadNChannelsList;  // list of P side dead channels
+    TArrayI  fDeadPChannelsList;  // list of N side dead channels
+
  private:
     AliITSCalibrationSSD(const AliITSCalibrationSSD &source); // copy constructor
     AliITSCalibrationSSD& operator=(const AliITSCalibrationSSD &source); // ass. op.