Changes made to support PreDigits (SDigits) plus other helpful changes.
[u/mrichter/AliRoot.git] / ITS / AliITSresponseSDD.h
1 #ifndef ALIITSRESPONSESDD_H
2 #define ALIITSRESPONSESDD_H
3
4 #include "TArrayF.h"
5 #include <TString.h>
6 #include <iostream.h>
7 #include "AliITSresponse.h"
8
9 // response for SDD
10
11 class AliITSresponseSDD :
12   public AliITSresponse {
13 public:
14   //
15   // Configuration methods
16   //
17   
18   AliITSresponseSDD();
19   virtual ~AliITSresponseSDD();
20
21   void SetElectronics(Int_t p1=1) {
22     // Electronics: Pascal (1) or OLA (2)
23     fElectronics=p1;
24   }
25   
26   Int_t Electronics() {
27     // Electronics: 1 = Pascal; 2 = OLA
28     return fElectronics;
29   }
30   
31   void    SetMaxAdc(Float_t p1=1024.) {
32     // Adc-count saturation value
33     fMaxAdc=p1;
34   }
35   Float_t MaxAdc()  {
36     // Get maximum Adc-count value
37     return fMaxAdc;
38   }                       
39   
40   void    SetChargeLoss(Float_t p1=0.0) {
41     // Set Linear Charge Loss Steepness  // 0.01 for 20%
42     fChargeLoss=p1;
43   }
44   Float_t ChargeLoss()  {
45     // Get Charge Loss Coefficient
46     return fChargeLoss;
47   }                       
48   
49   void    SetDynamicRange(Float_t p1=132.) {
50     // Set Dynamic Range
51     fDynamicRange=p1;
52   }
53   Float_t DynamicRange()  {
54     // Get Dynamic Range
55     return fDynamicRange;
56   }                       
57   
58   void    SetDiffCoeff(Float_t p1=3.23,Float_t p2=30.) {
59     // Diffusion coefficients
60     fDiffCoeff=p1;
61     fDiffCoeff1=p2;
62   }
63   void DiffCoeff(Float_t&diff,Float_t&diff1) {
64     // Get diffusion coefficients
65     diff = fDiffCoeff;
66     diff1 = fDiffCoeff1;
67   } 
68   
69   void    SetDriftSpeed(Float_t p1=7.3) {
70     // Drift velocity
71     fDriftSpeed=p1;
72   }
73   Float_t DriftSpeed() {
74     // drift speed
75     return fDriftSpeed;
76   } 
77   
78   void    SetTemperature(Float_t p1=23.) {
79     // Temperature
80     fTemperature=p1;
81   }
82   Float_t Temperature() {
83     // Get temperature
84     return fTemperature;
85   } 
86   
87   void    SetDataType(const char *data="simulated") {
88     // Type of data - real or simulated
89     fDataType=data;
90   }
91   const char  *DataType() const {
92     // Get data type
93     return fDataType.Data();
94   } 
95   
96   void SetParamOptions(const char *opt1="same",const char *opt2="same"){
97     // Parameters: "same" or read from "file" 
98     fParam1=opt1; fParam2=opt2;
99   }
100   void   ParamOptions(char *opt1,char *opt2) {
101     // options
102     strcpy(opt1,fParam1.Data()); strcpy(opt2,fParam2.Data());
103   }
104   
105   void  SetNoiseParam(Float_t n=0., Float_t b=20.){
106     // Noise and baseline  // 8.3 for ALICE with beam test measurements
107     fNoise=n; fBaseline=b;
108   }   
109   void  SetNoiseAfterElectronics(Float_t n=0.){
110     // Noise after electronics (ADC units) // 1.6 for ALICE from beam test measurements
111     fNoiseAfterEl=n;
112   }   
113   void  GetNoiseParam(Float_t &n, Float_t &b) {
114     // get noise param
115     n=fNoise; b=fBaseline;
116   }  
117   Float_t  GetNoiseAfterElectronics(){
118     // Noise after electronics (ADC units)
119     return fNoiseAfterEl;
120   }   
121
122   void  SetDo10to8(Bool_t bitcomp=kTRUE) {
123     // set the option for 10 to 8 bit compression
124     fBitComp = bitcomp;
125   }
126
127   Bool_t Do10to8() {
128     // get 10 to 8 compression option
129     return fBitComp;
130   }   
131   
132   void    SetZeroSupp (const char *opt="1D") {
133     // Zero-suppression option - could be 1D, 2D or non-ZS 
134     fOption=opt;
135   }
136   const char *ZeroSuppOption() const {
137     // Get zero-suppression option
138     return fOption.Data();
139   }
140   void  SetMinVal(Int_t mv=4) {
141     // Min value used in 2D - could be used as a threshold setting
142     fMinVal = mv;
143   }
144   Int_t  MinVal() {
145     // min val
146     return fMinVal;
147   }
148   
149   void   SetFilenames(const char *f1="",const char *f2="",const char *f3="") {
150     // Set filenames - input, output, parameters ....
151     fFileName1=f1; fFileName2=f2; fFileName3=f3;
152   }
153   void   Filenames(char *input,char *baseline,char *param) {
154     // Filenames
155    strcpy(input,fFileName1.Data());  strcpy(baseline,fFileName2.Data());  
156    strcpy(param,fFileName3.Data());
157   }     
158   
159   
160   void  SetOutputOption(Bool_t write=kFALSE) {
161     // set output option
162     fWrite = write;
163   }
164   Bool_t OutputOption()  {
165     // output option
166     return fWrite;
167   }
168   // 
169   // Compression parameters
170   void  SetCompressParam(Int_t cp[8]); 
171   void  GiveCompressParam(Int_t *x);
172   
173   //  
174   // Detector type response methods
175   void    SetNSigmaIntegration(Float_t p1=3.) {
176     // Set number of sigmas over which cluster disintegration is performed
177     fNsigmas=p1;
178   }
179   Float_t NSigmaIntegration() {
180     // Get number of sigmas over which cluster disintegration is performed
181     return fNsigmas;
182   }
183   void SetNLookUp(Int_t p1=121) {
184     // Set number of sigmas over which cluster disintegration is performed
185     fNcomps=p1;
186     fGaus = new TArrayF(fNcomps+1);
187     for(Int_t i=0; i<=fNcomps; i++) {
188       Float_t x = -fNsigmas + (2.*i*fNsigmas)/(fNcomps-1);
189       (*fGaus)[i] = exp(-((x*x)/2));
190  //     cout << "fGaus[" << i << "]: " << fGaus->At(i) << endl;
191     }
192   }
193     // Get number of intervals in which the gaussian lookup table is divided
194   Int_t GausNLookUp() {return fNcomps;}
195   
196   Float_t IntPH(Float_t eloss) {
197     // Pulse height from scored quantity (eloss)
198     return 0.;
199   }
200   Float_t IntXZ(AliITSsegmentation *) {
201     // Charge disintegration 
202     return 0.;
203   }
204   Float_t GausLookUp(Int_t i) {
205     if(i<0 || i>=fNcomps) return 0.;
206     return fGaus->At(i);
207   }
208   void SetDeadChannels(Int_t nmodules=0, Int_t nchips=0, Int_t nchannels=0);
209   void    PrintGains();
210   void    Print();
211
212
213 private:
214
215   AliITSresponseSDD(const AliITSresponseSDD &source); // copy constructor
216   AliITSresponseSDD& operator=(const AliITSresponseSDD &source); // ass. op.
217     
218 protected:
219   
220   static const Int_t fModules = 520;     // Total number of SDD modules
221   static const Int_t fChips = 4;        // Number of chips/module
222   static const Int_t fChannels = 64;    // Number of channels/chip
223   Float_t   fGain[fModules][fChips][fChannels];   // Array for channel gains
224
225   Int_t     fCPar[8];        // Hardware compression parameters
226   Float_t   fNoise;          // Noise
227   Float_t   fBaseline;       // Baseline
228   Float_t   fNoiseAfterEl;   // Noise after electronics
229   Float_t   fDynamicRange;   // Set Dynamic Range 
230   Float_t   fChargeLoss;     // Set Linear Coefficient for Charge Loss 
231   Float_t   fTemperature;    // Temperature 
232   Float_t   fDriftSpeed;     // Drift velocity
233   Int_t     fElectronics;    // Electronics
234   
235   Float_t    fMaxAdc;        // Adc saturation value
236   Float_t    fDiffCoeff;     // Diffusion Coefficient (scaling the time)
237   Float_t    fDiffCoeff1;    // Diffusion Coefficient (constant term)
238   Float_t    fNsigmas;       // Number of sigmas over which charge disintegration 
239                              // is performed 
240   TArrayF   *fGaus;          // Gaussian lookup table for signal generation
241   Int_t      fNcomps;        // Number of samplings along the gaussian
242   
243   Int_t      fMinVal;        // Min value used in 2D zero-suppression algo
244   
245   Bool_t     fWrite;         // Write option for the compression algorithms
246   Bool_t     fBitComp;       // 10 to 8 bit compression option
247
248   TString    fOption;        // Zero-suppresion option (1D, 2D or none)
249   TString    fParam1;        // Read baselines from file option
250   TString    fParam2;        // Read compression algo thresholds from file 
251   
252   TString         fDataType;         // data type - real or simulated
253   TString         fFileName1;        // input keys : run, module #
254   TString         fFileName2;        // baseline & noise val or output coded                                                 // signal or monitored bgr.
255   TString         fFileName3;        // param values or output coded signal 
256   
257   ClassDef(AliITSresponseSDD,2) // SDD response 
258     
259     };
260 #endif
261
262
263
264
265
266
267