]> git.uio.no Git - u/mrichter/AliRoot.git/blob - ITS/AliITSresponseSDD.h
Step size region by region.
[u/mrichter/AliRoot.git] / ITS / AliITSresponseSDD.h
1 #ifndef ALIITSRESPONSESDD_H
2 #define ALIITSRESPONSESDD_H
3  
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6  
7 /*
8   $Id$
9 */
10
11 #include "TArrayF.h"
12 #include "AliITSresponse.h"
13
14 // response for SDD
15
16 class AliITSresponseSDD : public AliITSresponse {
17   public:
18     //
19     // Configuration methods
20     //
21     AliITSresponseSDD();
22     AliITSresponseSDD(const char *dataType);
23     virtual ~AliITSresponseSDD();
24
25     // Implementation of virtual member functions declared in AliITSresponse 
26     virtual void SetElectronics(Int_t p1=1) 
27       {fElectronics=p1;  /* Electronics: Pascal (1) or OLA (2) */ }
28     virtual Int_t Electronics() const {// Electronics: 1 = Pascal; 2 = OLA
29         return fElectronics;}
30     virtual void    SetMaxAdc(Double_t p1) {// Adc-count saturation value
31         fMaxAdc=p1;}
32     virtual Double_t MaxAdc() const {// Get maximum Adc-count value
33         return fMaxAdc;}
34     virtual void    SetChargeLoss(Double_t p1) {
35         // Set Linear Charge Loss Steepness  // 0.01 for 20%
36         fChargeLoss=p1;}
37     Double_t ChargeLoss() const {// Get Charge Loss Coefficient
38         return fChargeLoss;}
39     virtual void    SetDynamicRange(Double_t p1) {// Set Dynamic Range
40         fDynamicRange = p1;}
41     virtual Double_t DynamicRange() const {// Get Dynamic Range
42         return fDynamicRange;}
43     // Diffusion coefficients
44     virtual void    SetDiffCoeff(Double_t p1, Double_t p2) 
45         {fDiffCoeff=p1;fDiffCoeff1=p2;}
46     // Get diffusion coefficients
47     virtual void DiffCoeff(Double_t&diff,Double_t&diff1) const 
48         {diff = fDiffCoeff;diff1 = fDiffCoeff1;}
49     virtual void    SetDriftSpeed(Double_t p1) {// Drift velocity
50         fDriftSpeed=p1;}
51     virtual Double_t DriftSpeed() const {// drift speed
52         return fDriftSpeed;}
53     virtual void SetParamOptions(const char *opt1,const char *opt2){
54         // Parameters: "same" or read from "file" 
55         fParam1=opt1; fParam2=opt2;}
56     virtual void   ParamOptions(char *opt1,char *opt2) const {// options
57         strcpy(opt1,fParam1.Data()); strcpy(opt2,fParam2.Data());}
58     virtual void  SetNoiseParam(Double_t n, Double_t b){
59         // Noise and baseline  // 10 for ALICE with beam test measurements 2001
60         fNoise=n; fBaseline=b;}
61         // get noise param
62     virtual void  GetNoiseParam(Double_t &n, Double_t &b) const
63         {n=fNoise; b=fBaseline;}
64     virtual Bool_t Do10to8() const {// get 10 to 8 compression option
65         return fBitComp;}
66     void    SetZeroSupp (const char *opt) {
67         // Zero-suppression option - could be 1D, 2D or non-ZS 
68         fOption=opt;}
69     const char *ZeroSuppOption() const {// Get zero-suppression option
70         return fOption.Data();}
71     virtual void   SetThresholds(Double_t  mv, Double_t /* b */){
72        // Min value used in 2D - could be used as a threshold setting
73         fMinVal = mv;}
74     virtual void   Thresholds(Double_t &  mv, Double_t & /* b */) const 
75       {mv = fMinVal;}
76     virtual void  GiveCompressParam(Int_t *x) const;
77     // Detector type response methods
78     virtual void    SetNSigmaIntegration(Double_t p1) {
79         // Set number of sigmas over which cluster disintegration is performed
80         fNsigmas=p1;}
81     virtual Double_t NSigmaIntegration() const {
82         // Get number of sigmas over which cluster disintegration is performed
83         return fNsigmas;}
84     virtual void SetNLookUp(Int_t p1);
85     // Get number of intervals in which the gaussian lookup table is divided
86     virtual Int_t GausNLookUp() const {return fNcomps;}
87     virtual Double_t GausLookUp(Int_t i) const  {
88         if(i<0 || i>=fNcomps) return 0.;return fGaus->At(i);}
89    
90 //Declaration of member functions peculiar to this class
91     Int_t Convert8to10(Int_t signal) const; //undo 10 to 8 bit SDD compresion
92     void  SetNoiseAfterElectronics(Double_t n=2.38){
93         // Noise after electronics (ADC units)
94         // 2.36 for ALICE from beam test measurements 2001
95         fNoiseAfterEl=n;}
96     Double_t  GetNoiseAfterElectronics() const {
97         // Noise after electronics (ADC units)
98         return fNoiseAfterEl;}
99     void  SetJitterError(Double_t jitter=20) {
100         // set Jitter error (20 um for ALICE from beam test measurements 2001)
101         fJitterError=jitter;}
102     Double_t  JitterError() const {// set Jitter error
103         return fJitterError;}
104     void  SetDo10to8(Bool_t bitcomp=kTRUE) {
105         // set the option for 10 to 8 bit compression
106         fBitComp = bitcomp;}
107     // Compression parameters
108     void  SetCompressParam(Int_t cp[8]); 
109     void SetDeadChannels(Int_t nmodules=0, Int_t nchips=0, Int_t nchannels=0);
110     Int_t GetDeadModules() const { return fDeadModules; }
111     Int_t GetDeadChips() const { return fDeadChips; }
112     Int_t GetDeadChannels() const { return fDeadChannels; }
113     Double_t Gain(Int_t mod,Int_t chip,Int_t ch)const 
114         {return fGain[mod][chip][ch]; }
115     // these functions should be move to AliITSsegmentationSDD
116     const Int_t Modules()const{return fgkModules;}//Total number of SDD modules
117     const Int_t Chips() const{return fgkChips;} // Number of chips/module
118     const Int_t Channels() const{ return fgkChannels;}//Number of channels/chip
119     //********
120     void    PrintGains();
121     void    Print();
122
123     // not implemented virtual methods (devlared in the mother class
124     virtual  void   SetDetParam(Double_t *)
125       {NotImplemented("SetDetParam");}
126     virtual void   GetDetParam(Double_t *) const 
127       {NotImplemented("GetDetParam");}
128     virtual  void   SetNDetParam(Int_t /* n */)
129       {NotImplemented("SetNDetParam");}
130     virtual Int_t  NDetParam() const
131       {NotImplemented("NDetParam"); return 0;}
132     virtual void    SetSigmaSpread(Double_t, Double_t) 
133       {NotImplemented("SetSigmaSpread");}
134     virtual void    SigmaSpread(Double_t & /* p1 */,Double_t & /* p2 */) const 
135       {NotImplemented("SigmaSpread");}
136
137  protected:
138     // these statis const should be move to AliITSsegmentationSDD
139     static const Int_t fgkModules = 520;     // Total number of SDD modules
140     static const Int_t fgkChips = 4;        // Number of chips/module
141     static const Int_t fgkChannels = 64;    // Number of channels/chip
142     //*******
143     static const Int_t fgkMaxAdcDefault; // default for fMaxAdc
144     static const Double_t fgkDynamicRangeDefault; // default for fDynamicRange
145     static const Double_t fgkfChargeLossDefault; // default for fChargeLoss
146     static const Double_t fgkDiffCoeffDefault; // default for fDiffCoeff
147     static const Double_t fgkDiffCoeff1Default; // default for fDiffCoeff1 
148     static const Double_t fgkTemperatureDefault; // default for fT (Kelvin)
149     static const TString fgkParam1Default; // default for fParam1
150     static const TString fgkParam2Default; // default for fParam2
151     static const Double_t fgkNoiseDefault; // default for fNoise
152     static const Double_t fgkBaselineDefault; // default for fBaseline
153     static const TString fgkOptionDefault; // default for fOption
154     static const Double_t fgkMinValDefault; // default for fMinVal
155     static const Double_t fgkDriftSpeedDefault; // default for fDriftSpeed
156     static const Double_t fgkNsigmasDefault; //default for fNsigmas
157     static const Int_t fgkNcompsDefault; //default for fNcomps
158     //
159
160     Int_t fDeadModules;                   // Total number of dead SDD modules
161     Int_t fDeadChips;                     // Number of dead chips
162     Int_t fDeadChannels;                  // Number of dead channels
163     Double_t fGain[fgkModules][fgkChips][fgkChannels];//Array for channel gains
164     Int_t     fCPar[8];        // Hardware compression parameters
165     Double_t   fNoise;          // Noise
166     Double_t   fBaseline;       // Baseline
167     Double_t   fNoiseAfterEl;   // Noise after electronics
168     Double_t   fJitterError;    // jitter error
169     Double_t   fDynamicRange;   // Set Dynamic Range 
170     Double_t   fChargeLoss;     // Set Linear Coefficient for Charge Loss 
171     Double_t   fDriftSpeed;     // Drift velocity
172     Int_t     fElectronics;    // Electronics
173     Double_t   fMaxAdc;         // Adc saturation value
174     Double_t   fDiffCoeff;      // Diffusion Coefficient (scaling the time)
175     Double_t   fDiffCoeff1;     // Diffusion Coefficient (constant term)
176     Double_t   fNsigmas;   // Number of sigmas over which charge disintegration
177                                // is performed
178     TArrayF   *fGaus;          // Gaussian lookup table for signal generation
179     Int_t      fNcomps;        // Number of samplings along the gaussian
180     Double_t      fMinVal;        // Min value used in 2D zero-suppression algo
181     Bool_t     fBitComp;       // 10 to 8 bit compression option
182     TString    fOption;        // Zero-suppresion option (1D, 2D or none)
183     TString    fParam1;        // Read baselines from file option
184     TString    fParam2;        // Read compression algo thresholds from file
185
186  private:
187     AliITSresponseSDD(const AliITSresponseSDD &ob); // copy constructor
188     AliITSresponseSDD& operator=(const AliITSresponseSDD & /* source */); // ass. op.
189
190
191     ClassDef(AliITSresponseSDD,3) // SDD response 
192     
193     };
194 #endif