/[PAMELA software]/PamelaDigitizer/Digitizer.h
ViewVC logotype

Annotation of /PamelaDigitizer/Digitizer.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.5 - (hide annotations) (download)
Wed May 21 09:50:44 2008 UTC (16 years, 6 months ago) by pamelats
Branch: MAIN
Changes since 1.4: +449 -442 lines
File MIME type: text/plain
*** empty log message ***

1 pamelats 1.5 #ifndef DIGITIZER_H
2     #define DIGITIZER_H
3     #define nspe 200
4     #define ntof 200
5     #define ncat 50
6     #define ncas 50
7     #define ncar 100
8     #define ncal 1000
9     #define nnd 200
10     #define nstr 1000
11    
12     #include <iostream>
13     #include <fstream>
14    
15    
16     #include <stdlib.h>
17     #include <string.h>
18     #include <ctype.h>
19    
20     #include "Riostream.h"
21     #include "TFile.h"
22     #include "TDirectory.h"
23     #include "TTree.h"
24     #include "TLeafI.h"
25     #include "TH1.h"
26     #include "TH2.h"
27     #include "TMath.h"
28     #include "TSQLServer.h"
29     //
30     #include "CRC.h"
31    
32     class Digitizer{
33    
34    
35     protected:
36    
37     // in and output data members
38    
39     TTree* fhBookTree;
40     ofstream fOutputfile;
41    
42    
43     static const int fPSCUbuffer = 8;
44     static const int fTRIGGERbuffer = 152; //(bytes)
45     static const int fTOFbuffer = 276; //(bytes)
46     static const int fRunHeaderbuffer = 27; // ??
47     static const int fRunTrailerbuffer = 17;
48     static const int fACbuffer = 128;
49     static const int fNDbuffer = 6;
50     static const int fS4buffer = 3;
51     static const int fPADbuffer = 64; // changed from 32!
52     static const int fCALObuffer = 10000;
53     static const float fCALOGeV2MIPratio = 0.0001059994;
54    
55     static const int fDummybuffer = 64;
56    
57     // dimensions of vectors and buffer, parameter for tracker
58     static const int fTRACKbuffer = 50000; // to be optimized
59     static const int fNviews = 12;
60     static const int fNstrips_view = 3072;
61     static const int fNladder=3;
62     static const int fNstrips_ladder = 1024;
63     static const int fCutzop= 4;
64     static const int fCutclu= 7;
65     static const int fNclst = 2;
66    
67    
68     static const float fAvePedex = 1200.; // to be adjusted
69     static const float fSigmaPedex = 1.; // to be adjusted
70     static const float fAveSigmax = 4.3; // to be adjusted
71     static const float fSigmaSigmax =0.1; // to be adjusted
72     static const float fAvePedey = 2800.; // to be adjusted
73     static const float fSigmaPedey = 1.; // to be adjusted
74     static const float fAveSigmay = 8.; // to be adjusted
75     static const float fSigmaSigmay =0.1; // to be adjusted
76     static const float fSigmaCommon = 25.; // to be adjusted
77    
78    
79     // end dimensions of vectors and buffer, parameter for tracker
80    
81     UShort_t fDataPSCU[fPSCUbuffer];
82     UInt_t fCounter;
83     UInt_t fCounterPhys;
84     UInt_t fOBT;
85    
86     UShort_t fDataPadding[fPADbuffer];
87     UInt_t fPadding;
88    
89     UShort_t fDataRunHeader[fRunHeaderbuffer];
90     UShort_t fDataRunTrailer[fRunTrailerbuffer];
91    
92     UChar_t fDataTrigger[fTRIGGERbuffer];
93     UChar_t fDataTof[fTOFbuffer];
94    
95     UShort_t fDataAC[fACbuffer];
96     UShort_t fDataND[fNDbuffer];
97     UShort_t fDataS4[fS4buffer];
98     UShort_t fDataDummy[fDummybuffer];
99    
100     UShort_t fDataCALO[fCALObuffer];
101     UShort_t fSecCalo[4];
102     UShort_t fSecCALOLength[4];
103     UShort_t fCALOlength;
104     UShort_t fModCalo;
105     Float_t fCalomip[2][22][96]; // ADC to MIP conversion values
106     Float_t fcalped[4][11][96]; // pedestals
107     Float_t fcalgood[4][11][96]; // used or not in the common noise calculation
108     Float_t fcalthr[4][11][6]; // thresholds
109     Float_t fcalrms[4][11][96]; // rms
110     Float_t fcalbase[4][11][6]; // baselines (common noise) during calibration
111     Float_t fcalvar[4][11][6]; // variance
112     Float_t fcstwerr[4]; // status word
113     Float_t fcperror[4]; // unpacking error flag
114     UInt_t fGivenCaloCalib;
115    
116     char* fFilename;
117     UInt_t fLen;
118    
119    
120    
121     // vectors, buffer, and variables for tracker
122     UShort_t fDataTrack[fTRACKbuffer];
123     UInt_t fTracklength; // Actual length of buffer
124     Float_t fPedeTrack[fNviews][fNstrips_view];
125     Float_t fSigmaTrack[fNviews][fNstrips_view];
126     // UShort_t fBadTrack[fNviews][fNstrips_view];
127     Float_t fMipCor[fNladder][fNviews] ;
128     // end vectors, buffer, and variables for tracker
129    
130    
131    
132     //
133     // db variables
134     //
135     TSQLServer *fDbc;
136     // Declaration of leave types
137     Int_t Irun;
138     Int_t Ievnt;
139     Int_t Ipa;
140     Float_t X0;
141     Float_t Y0;
142     Float_t Z0;
143     Float_t Theta;
144     Float_t Phi;
145     Float_t P0;
146     Int_t Nthtof;
147     UChar_t Ipltof[ntof];
148     UChar_t Ipaddle[ntof];
149     UChar_t Ipartof[ntof];
150     Float_t Xintof[ntof];
151     Float_t Yintof[ntof];
152     Float_t Zintof[ntof];
153     Float_t Xouttof[ntof];
154     Float_t Youttof[ntof];
155     Float_t Zouttof[ntof];
156     Float_t Ereltof[ntof];
157     Float_t Timetof[ntof];
158     Float_t Pathtof[ntof];
159     Float_t P0tof[ntof];
160     Int_t Nthcat;
161     UChar_t Iparcat[ncat];
162     UChar_t Icat[ncat];
163     Float_t Xincat[ncat];
164     Float_t Yincat[ncat];
165     Float_t Zincat[ncat];
166     Float_t Xoutcat[ncat];
167     Float_t Youtcat[ncat];
168     Float_t Zoutcat[ncat];
169     Float_t Erelcat[ncat];
170     Float_t Timecat[ncat];
171     Float_t Pathcat[ncat];
172     Float_t P0cat[ncat];
173     Int_t Nthcas;
174     UChar_t Iparcas[ncas];
175     UChar_t Icas[ncas];
176     Float_t Xincas[ncas];
177     Float_t Yincas[ncas];
178     Float_t Zincas[ncas];
179     Float_t Xoutcas[ncas];
180     Float_t Youtcas[ncas];
181     Float_t Zoutcas[ncas];
182     Float_t Erelcas[ncas];
183     Float_t Timecas[ncas];
184     Float_t Pathcas[ncas];
185     Float_t P0cas[ncas];
186     Int_t Nthspe;
187     UChar_t Iparspe[nspe];
188     UChar_t Itrpb[nspe];
189     UChar_t Itrsl[nspe];
190     UChar_t Itspa[nspe];
191     Float_t Xinspe[nspe];
192     Float_t Yinspe[nspe];
193     Float_t Zinspe[nspe];
194     Float_t Xoutspe[nspe];
195     Float_t Youtspe[nspe];
196     Float_t Zoutspe[nspe];
197     Float_t Xavspe[nspe];
198     Float_t Yavspe[nspe];
199     Float_t Zavspe[nspe];
200     Float_t Erelspe[nspe];
201     Float_t Pathspe[nspe];
202     Float_t P0spe[nspe];
203     UChar_t Nxmult[nspe];
204     UChar_t Nymult[nspe];
205     Int_t Nstrpx;
206     UChar_t Npstripx[nstr];
207     UChar_t Ntstripx[nstr];
208     UShort_t Istripx[nstr];
209     //UInt_t Istripx[nstr];
210     Float_t Qstripx[nstr];
211     Float_t Xstripx[nstr];
212     Int_t Nstrpy;
213     UChar_t Npstripy[nstr];
214     UChar_t Ntstripy[nstr];
215     UShort_t Istripy[nstr];
216     //UInt_t Istripy[nstr];
217     Float_t Qstripy[nstr];
218     Float_t Ystripy[nstr];
219     Int_t Nthcali;
220     UChar_t Icaplane[4224];
221     UChar_t Icastrip[4224];
222     Int_t Icamod[4224];
223     Float_t Enestrip[4224];
224     Int_t Nthcal;
225     UChar_t Icapl[ncal];
226     UChar_t Icasi[ncal];
227     UChar_t Icast[ncal];
228     Float_t Xincal[ncal];
229     Float_t Yincal[ncal];
230     Float_t Zincal[ncal];
231     Float_t Erelcal[ncal];
232     Int_t Nthnd;
233     UChar_t Itubend[nnd];
234     UChar_t Iparnd[nnd];
235     Float_t Xinnd[nnd];
236     Float_t Yinnd[nnd];
237     Float_t Zinnd[nnd];
238     Float_t Xoutnd[nnd];
239     Float_t Youtnd[nnd];
240     Float_t Zoutnd[nnd];
241     Float_t Erelnd[nnd];
242     Float_t Timend[nnd];
243     Float_t Pathnd[nnd];
244     Float_t P0nd[nnd];
245     Int_t Nthcard;
246     UChar_t Iparcard[ncar];
247     UChar_t Icard[ncar];
248     Float_t Xincard[ncar];
249     Float_t Yincard[ncar];
250     Float_t Zincard[ncar];
251     Float_t Xoutcard[ncar];
252     Float_t Youtcard[ncar];
253     Float_t Zoutcard[ncar];
254     Float_t Erelcard[ncar];
255     Float_t Timecard[ncar];
256     Float_t Pathcard[ncar];
257     Float_t P0card[ncar];
258    
259     // List of branches
260     TBranch *b_Irun; //!
261     TBranch *b_Ievnt; //!
262     TBranch *b_Ipa; //!
263     TBranch *b_X0; //!
264     TBranch *b_Y0; //!
265     TBranch *b_Z0; //!
266     TBranch *b_Theta; //!
267     TBranch *b_Phi; //!
268     TBranch *b_P0; //!
269     TBranch *b_Nthtof; //!
270     TBranch *b_Ipltof; //!
271     TBranch *b_Ipaddle; //!
272     TBranch *b_Ipartof; //!
273     TBranch *b_Xintof; //!
274     TBranch *b_Yintof; //!
275     TBranch *b_Zintof; //!
276     TBranch *b_Xouttof; //!
277     TBranch *b_Youttof; //!
278     TBranch *b_Zouttof; //!
279     TBranch *b_Ereltof; //!
280     TBranch *b_Timetof; //!
281     TBranch *b_Pathtof; //!
282     TBranch *b_P0tof; //!
283     TBranch *b_Nthcat; //!
284     TBranch *b_Iparcat; //!
285     TBranch *b_Icat; //!
286     TBranch *b_Xincat; //!
287     TBranch *b_Yincat; //!
288     TBranch *b_Zincat; //!
289     TBranch *b_Xoutcat; //!
290     TBranch *b_Youtcat; //!
291     TBranch *b_Zoutcat; //!
292     TBranch *b_Erelcat; //!
293     TBranch *b_Timecat; //!
294     TBranch *b_Pathcat; //!
295     TBranch *b_P0cat; //!
296     TBranch *b_Nthcas; //!
297     TBranch *b_Iparcas; //!
298     TBranch *b_Icas; //!
299     TBranch *b_Xincas; //!
300     TBranch *b_Yincas; //!
301     TBranch *b_Zincas; //!
302     TBranch *b_Xoutcas; //!
303     TBranch *b_Youtcas; //!
304     TBranch *b_Zoutcas; //!
305     TBranch *b_Erelcas; //!
306     TBranch *b_Timecas; //!
307     TBranch *b_Pathcas; //!
308     TBranch *b_P0cas; //!
309     TBranch *b_Nthspe; //!
310     TBranch *b_Iparspe; //!
311     TBranch *b_Itrpb; //!
312     TBranch *b_Itrsl; //!
313     TBranch *b_Itspa; //!
314     TBranch *b_Xinspe; //!
315     TBranch *b_Yinspe; //!
316     TBranch *b_Zinspe; //!
317     TBranch *b_Xoutspe; //!
318     TBranch *b_Youtspe; //!
319     TBranch *b_Zoutspe; //!
320     TBranch *b_Xavspe; //!
321     TBranch *b_Yavspe; //!
322     TBranch *b_Zavspe; //!
323     TBranch *b_Erelspe; //!
324     TBranch *b_Pathspe; //!
325     TBranch *b_P0spe; //!
326     TBranch *b_Nxmult; //!
327     TBranch *b_Nymult; //!
328     TBranch *b_Nstrpx; //!
329     TBranch *b_Npstripx; //!
330     TBranch *b_Ntstripx; //!
331     TBranch *b_Istripx; //!
332     TBranch *b_Qstripx; //!
333     TBranch *b_Xstripx; //!
334     TBranch *b_Nstrpy; //!
335     TBranch *b_Npstripy; //!
336     TBranch *b_Ntstripy; //!
337     TBranch *b_Istripy; //!
338     TBranch *b_Qstripy; //!
339     TBranch *b_Ystripy; //!
340     TBranch *b_Nthcali; //!
341     TBranch *b_Icaplane; //!
342     TBranch *b_Icastrip; //!
343     TBranch *b_Icamod; //!
344     TBranch *b_Enestrip; //!
345     TBranch *b_Nthcal; //!
346     TBranch *b_Icapl; //!
347     TBranch *b_Icasi; //!
348     TBranch *b_Icast; //!
349     TBranch *b_Xincal; //!
350     TBranch *b_Yincal; //!
351     TBranch *b_Zincal; //!
352     TBranch *b_Erelcal; //!
353     TBranch *b_Nthnd; //!
354     TBranch *b_Itubend; //!
355     TBranch *b_Iparnd; //!
356     TBranch *b_Xinnd; //!
357     TBranch *b_Yinnd; //!
358     TBranch *b_Zinnd; //!
359     TBranch *b_Xoutnd; //!
360     TBranch *b_Youtnd; //!
361     TBranch *b_Zoutnd; //!
362     TBranch *b_Erelnd; //!
363     TBranch *b_Timend; //!
364     TBranch *b_Pathnd; //!
365     TBranch *b_P0nd; //!
366     TBranch *b_Nthcard; //!
367     TBranch *b_Iparcard; //!
368     TBranch *b_Icard; //!
369     TBranch *b_Xincard; //!
370     TBranch *b_Yincard; //!
371     TBranch *b_Zincard; //!
372     TBranch *b_Xoutcard; //!
373     TBranch *b_Youtcard; //!
374     TBranch *b_Zoutcard; //!
375     TBranch *b_Erelcard; //!
376     TBranch *b_Timecard; //!
377     TBranch *b_Pathcard; //!
378     TBranch *b_P0card; //!
379     //
380    
381     public:
382    
383     Digitizer(TTree *tree, char* &file_raw); // constructor
384    
385     void Close();
386     void SetTree(TTree *tree) {fhBookTree = tree;};
387     void Loop(); // does the looping over events
388    
389     void DigitizePSCU(UInt_t, UChar_t, UShort_t*); // digitizes the CPU header
390     void AddPadding();
391     void DigitizeRunHeader();
392     void DigitizeRunTrailer();
393     // void DigitizeTRIGGER(); // SO nov'07: included in DigitizeTOF
394     void DigitizeTOF(int np,float *atte1,float *atte2,float *lambda1,float *lambda2); // digitizes the tof
395     void LoadTOFCalib(int np,float *atte1,float *atte2,float *lambda1,float *lambda2); //load TOF calibration
396     void DigitizeAC(); // digitizes the anti-counter
397     void DigitizeND(); // digitizes the ND
398     void DigitizeS4(); // digitizes the S4
399     void DigitizeCALO(); // digitizes the calorimeter
400     void DigitizeCALOCALIB(); // digitizes the calorimeter calibration
401     void DigitizeCALORAW(); // digitizes the calorimeter
402     void DigitizeCALOCOMPRESS(); // digitizes the calorimeter
403     void DigitizeCALOFULL(); // digitizes the calorimeter
404    
405     Float_t GetCALOen(Int_t sec, Int_t plane, Int_t strip); // retrieve calorimeter strip energy from the simulation
406     void ClearCaloCalib(Int_t s); // clear calorimeter calibration variables for section s
407     Int_t CaloLoadCalib(Int_t s,TString fcalname, UInt_t calibno); // load calorimeter calibration variables for section s from file fcalname entry calibno
408     void CaloLoadCalib();
409    
410     void DigitizeDummy(); // digitizes the dummy detector
411    
412     // ToF methods:
413     UChar_t EvaluateCrcTof(UChar_t*);
414     void Paddle2Pmt(Int_t, Int_t, Int_t*, Int_t*);
415     //void Paddle2Pmt(Int_t, Int_t, Int_t &pmtleft, Int_t &pmtright);
416     UChar_t Bin2GrayTof(UChar_t, UChar_t);
417     void Crc8Tof(UChar_t*, UChar_t*);
418    
419     UChar_t EvaluateCrcTrigger(UChar_t*, Int_t);
420    
421     UShort_t* GetDataAC(){return(fDataAC);}; // retrieve AC data (not used right now)
422    
423     UShort_t* GetDataCALO(){return(fDataCALO);}; // retrieve AC data (not used right now)
424    
425     // methods for Tracker
426     void DigitizeTrack(); // digitizes the Tracker data
427     void DigitizeTrackCalib(Int_t ii); // digitizes the Tracker calibration (x view or y view)
428     void WriteTrackCalib(); // writes Tracker calibration to file
429     void ClearTrackCalib(); // clear tracker calibration
430     void LoadTrackCalib(); // load in memmory tracker calibration
431     void CompressTrackData(Float_t AdcTrack[fNviews][fNstrips_view]); // Compress and Digitize data of one Ladder in turn
432     void LoadMipCor(); // Load the strip calibration in memory
433     Float_t SaturationTrackx(Float_t ADC); // Saturation simulation
434     Float_t SaturationTracky(Float_t ADC); // Saturation simulation
435     // end methods for Tracker
436    
437     void WriteRunHeader();
438     void WriteRunTrailer();
439     void WriteData(); // writes data to file
440     void ReadData();
441    
442     //
443    
444     };
445    
446     #endif /* DIGITIZER_H */
447    
448    
449    

  ViewVC Help
Powered by ViewVC 1.1.23