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