| 1 |
#ifndef DIGITIZER_H |
#ifndef DIGITIZER_H |
| 2 |
#define DIGITIZER_H |
#define DIGITIZER_H |
|
#define nspe 200 |
|
|
#define ntof 200 |
|
|
#define ncat 50 |
|
|
#define ncas 50 |
|
|
#define ncar 100 |
|
|
#define ncal 1000 |
|
|
#define nnd 200 |
|
|
#define nstr 1000 |
|
| 3 |
|
|
| 4 |
#include <iostream> |
#include <iostream> |
| 5 |
|
#include <istream> |
| 6 |
#include <fstream> |
#include <fstream> |
| 7 |
|
|
| 8 |
|
|
| 17 |
#include "TLeafI.h" |
#include "TLeafI.h" |
| 18 |
#include "TH1.h" |
#include "TH1.h" |
| 19 |
#include "TH2.h" |
#include "TH2.h" |
| 20 |
|
#include "TF1.h" |
| 21 |
#include "TMath.h" |
#include "TMath.h" |
| 22 |
#include "TSQLServer.h" |
#include "TSQLServer.h" |
| 23 |
// |
#include "TGraph.h" |
| 24 |
|
#include "GLTables.h" |
| 25 |
|
#include "TSystem.h" |
| 26 |
|
#include "TRandom.h" |
| 27 |
|
|
| 28 |
#include "CRC.h" |
#include "CRC.h" |
| 29 |
|
#include <CalibCalPedEvent.h> |
| 30 |
|
#include "CalibTrk1Event.h" |
| 31 |
|
#include "CalibTrk2Event.h" |
| 32 |
|
|
| 33 |
class Digitizer{ |
class Digitizer{ |
| 34 |
|
|
|
|
|
| 35 |
protected: |
protected: |
| 36 |
|
|
| 37 |
// in and output data members |
// in and output data members |
| 38 |
|
|
| 39 |
TTree* fhBookTree; |
TTree* fhBookTree; |
| 40 |
ofstream fOutputfile; |
ofstream fOutputfile; |
| 41 |
|
int *nspe; |
| 42 |
|
int *ntof; |
| 43 |
|
int *ncat; |
| 44 |
|
int *ncas; |
| 45 |
|
int *ncar; |
| 46 |
|
int *ncal; |
| 47 |
|
int *nnd; |
| 48 |
|
int *nstr; |
| 49 |
|
|
| 50 |
static const int fPSCUbuffer = 8; |
static const int fPSCUbuffer = 8; |
| 51 |
static const int fTRIGGERbuffer = 152; //(bytes) |
static const int fTRIGGERbuffer = 152; //(bytes) |
| 52 |
static const int fTOFbuffer = 276; //(bytes) |
static const int fTOFbuffer = 276; //(bytes) |
|
static const int fRunHeaderbuffer = 27; // ?? |
|
|
static const int fRunTrailerbuffer = 17; |
|
| 53 |
static const int fACbuffer = 128; |
static const int fACbuffer = 128; |
| 54 |
static const int fNDbuffer = 6; |
static const int fNDbuffer = 6; |
| 55 |
static const int fS4buffer = 3; |
static const int fS4buffer = 3; |
| 91 |
UShort_t fDataPadding[fPADbuffer]; |
UShort_t fDataPadding[fPADbuffer]; |
| 92 |
UInt_t fPadding; |
UInt_t fPadding; |
| 93 |
|
|
|
UShort_t fDataRunHeader[fRunHeaderbuffer]; |
|
|
UShort_t fDataRunTrailer[fRunTrailerbuffer]; |
|
|
|
|
| 94 |
UChar_t fDataTrigger[fTRIGGERbuffer]; |
UChar_t fDataTrigger[fTRIGGERbuffer]; |
| 95 |
UChar_t fDataTof[fTOFbuffer]; |
UChar_t fDataTof[fTOFbuffer]; |
| 96 |
|
|
| 146 |
Float_t Phi; |
Float_t Phi; |
| 147 |
Float_t P0; |
Float_t P0; |
| 148 |
Int_t Nthtof; |
Int_t Nthtof; |
| 149 |
UChar_t Ipltof[ntof]; |
UChar_t *Ipltof; |
| 150 |
UChar_t Ipaddle[ntof]; |
UChar_t *Ipaddle; |
| 151 |
UChar_t Ipartof[ntof]; |
// UChar_t *Ipartof;//DPMJET |
| 152 |
Float_t Xintof[ntof]; |
UShort_t *Ipartof; |
| 153 |
Float_t Yintof[ntof]; |
Float_t *Xintof; |
| 154 |
Float_t Zintof[ntof]; |
Float_t *Yintof; |
| 155 |
Float_t Xouttof[ntof]; |
Float_t *Zintof; |
| 156 |
Float_t Youttof[ntof]; |
Float_t *Xouttof; |
| 157 |
Float_t Zouttof[ntof]; |
Float_t *Youttof; |
| 158 |
Float_t Ereltof[ntof]; |
Float_t *Zouttof; |
| 159 |
Float_t Timetof[ntof]; |
Float_t *Ereltof; |
| 160 |
Float_t Pathtof[ntof]; |
Float_t *Timetof; |
| 161 |
Float_t P0tof[ntof]; |
Float_t *Pathtof; |
| 162 |
|
Float_t *P0tof; |
| 163 |
Int_t Nthcat; |
Int_t Nthcat; |
| 164 |
UChar_t Iparcat[ncat]; |
UChar_t *Iparcat; |
| 165 |
UChar_t Icat[ncat]; |
UChar_t *Icat; |
| 166 |
Float_t Xincat[ncat]; |
Float_t *Xincat; |
| 167 |
Float_t Yincat[ncat]; |
Float_t *Yincat; |
| 168 |
Float_t Zincat[ncat]; |
Float_t *Zincat; |
| 169 |
Float_t Xoutcat[ncat]; |
Float_t *Xoutcat; |
| 170 |
Float_t Youtcat[ncat]; |
Float_t *Youtcat; |
| 171 |
Float_t Zoutcat[ncat]; |
Float_t *Zoutcat; |
| 172 |
Float_t Erelcat[ncat]; |
Float_t *Erelcat; |
| 173 |
Float_t Timecat[ncat]; |
Float_t *Timecat; |
| 174 |
Float_t Pathcat[ncat]; |
Float_t *Pathcat; |
| 175 |
Float_t P0cat[ncat]; |
Float_t *P0cat; |
| 176 |
Int_t Nthcas; |
Int_t Nthcas; |
| 177 |
UChar_t Iparcas[ncas]; |
UChar_t *Iparcas; |
| 178 |
UChar_t Icas[ncas]; |
UChar_t *Icas; |
| 179 |
Float_t Xincas[ncas]; |
Float_t *Xincas; |
| 180 |
Float_t Yincas[ncas]; |
Float_t *Yincas; |
| 181 |
Float_t Zincas[ncas]; |
Float_t *Zincas; |
| 182 |
Float_t Xoutcas[ncas]; |
Float_t *Xoutcas; |
| 183 |
Float_t Youtcas[ncas]; |
Float_t *Youtcas; |
| 184 |
Float_t Zoutcas[ncas]; |
Float_t *Zoutcas; |
| 185 |
Float_t Erelcas[ncas]; |
Float_t *Erelcas; |
| 186 |
Float_t Timecas[ncas]; |
Float_t *Timecas; |
| 187 |
Float_t Pathcas[ncas]; |
Float_t *Pathcas; |
| 188 |
Float_t P0cas[ncas]; |
Float_t *P0cas; |
| 189 |
Int_t Nthspe; |
Int_t Nthspe; |
| 190 |
UChar_t Iparspe[nspe]; |
// UChar_t *Iparspe;//DPMJET |
| 191 |
UChar_t Itrpb[nspe]; |
// UShort_t *Iparspe; |
| 192 |
UChar_t Itrsl[nspe]; |
UChar_t *Itrpb; |
| 193 |
UChar_t Itspa[nspe]; |
UChar_t *Itrsl; |
| 194 |
Float_t Xinspe[nspe]; |
UChar_t *Itspa; |
| 195 |
Float_t Yinspe[nspe]; |
Float_t *Xinspe; |
| 196 |
Float_t Zinspe[nspe]; |
Float_t *Yinspe; |
| 197 |
Float_t Xoutspe[nspe]; |
Float_t *Zinspe; |
| 198 |
Float_t Youtspe[nspe]; |
Float_t *Xoutspe; |
| 199 |
Float_t Zoutspe[nspe]; |
Float_t *Youtspe; |
| 200 |
Float_t Xavspe[nspe]; |
Float_t *Zoutspe; |
| 201 |
Float_t Yavspe[nspe]; |
Float_t *Xavspe; |
| 202 |
Float_t Zavspe[nspe]; |
Float_t *Yavspe; |
| 203 |
Float_t Erelspe[nspe]; |
Float_t *Zavspe; |
| 204 |
Float_t Pathspe[nspe]; |
Float_t *Erelspe; |
| 205 |
Float_t P0spe[nspe]; |
Float_t *Pathspe; |
| 206 |
UChar_t Nxmult[nspe]; |
Float_t *P0spe; |
| 207 |
UChar_t Nymult[nspe]; |
UChar_t *Nxmult; |
| 208 |
|
UChar_t *Nymult; |
| 209 |
Int_t Nstrpx; |
Int_t Nstrpx; |
| 210 |
UChar_t Npstripx[nstr]; |
UChar_t *Npstripx; |
| 211 |
UChar_t Ntstripx[nstr]; |
UChar_t *Ntstripx; |
| 212 |
UShort_t Istripx[nstr]; |
UShort_t *Istripx; |
| 213 |
//UInt_t Istripx[nstr]; |
Float_t *Qstripx; |
| 214 |
Float_t Qstripx[nstr]; |
Float_t *Xstripx; |
|
Float_t Xstripx[nstr]; |
|
| 215 |
Int_t Nstrpy; |
Int_t Nstrpy; |
| 216 |
UChar_t Npstripy[nstr]; |
UChar_t *Npstripy; |
| 217 |
UChar_t Ntstripy[nstr]; |
UChar_t *Ntstripy; |
| 218 |
UShort_t Istripy[nstr]; |
UShort_t *Istripy; |
| 219 |
//UInt_t Istripy[nstr]; |
Float_t *Qstripy; |
| 220 |
Float_t Qstripy[nstr]; |
Float_t *Ystripy; |
|
Float_t Ystripy[nstr]; |
|
| 221 |
Int_t Nthcali; |
Int_t Nthcali; |
| 222 |
UChar_t Icaplane[4224]; |
UChar_t Icaplane[4224]; |
| 223 |
UChar_t Icastrip[4224]; |
UChar_t Icastrip[4224]; |
| 224 |
Int_t Icamod[4224]; |
Int_t Icamod[4224]; |
| 225 |
Float_t Enestrip[4224]; |
Float_t Enestrip[4224]; |
| 226 |
Int_t Nthcal; |
Int_t Nthcal; |
| 227 |
UChar_t Icapl[ncal]; |
UChar_t *Icapl; |
| 228 |
UChar_t Icasi[ncal]; |
UChar_t *Icasi; |
| 229 |
UChar_t Icast[ncal]; |
UChar_t *Icast; |
| 230 |
Float_t Xincal[ncal]; |
Float_t *Xincal; |
| 231 |
Float_t Yincal[ncal]; |
Float_t *Yincal; |
| 232 |
Float_t Zincal[ncal]; |
Float_t *Zincal; |
| 233 |
Float_t Erelcal[ncal]; |
Float_t *Erelcal; |
| 234 |
Int_t Nthnd; |
Int_t Nthnd; |
| 235 |
UChar_t Itubend[nnd]; |
UChar_t *Itubend; |
| 236 |
UChar_t Iparnd[nnd]; |
UChar_t *Iparnd; |
| 237 |
Float_t Xinnd[nnd]; |
Float_t *Xinnd; |
| 238 |
Float_t Yinnd[nnd]; |
Float_t *Yinnd; |
| 239 |
Float_t Zinnd[nnd]; |
Float_t *Zinnd; |
| 240 |
Float_t Xoutnd[nnd]; |
Float_t *Xoutnd; |
| 241 |
Float_t Youtnd[nnd]; |
Float_t *Youtnd; |
| 242 |
Float_t Zoutnd[nnd]; |
Float_t *Zoutnd; |
| 243 |
Float_t Erelnd[nnd]; |
Float_t *Erelnd; |
| 244 |
Float_t Timend[nnd]; |
Float_t *Timend; |
| 245 |
Float_t Pathnd[nnd]; |
Float_t *Pathnd; |
| 246 |
Float_t P0nd[nnd]; |
Float_t *P0nd; |
| 247 |
Int_t Nthcard; |
Int_t Nthcard; |
| 248 |
UChar_t Iparcard[ncar]; |
UChar_t *Iparcard; |
| 249 |
UChar_t Icard[ncar]; |
UChar_t *Icard; |
| 250 |
Float_t Xincard[ncar]; |
Float_t *Xincard; |
| 251 |
Float_t Yincard[ncar]; |
Float_t *Yincard; |
| 252 |
Float_t Zincard[ncar]; |
Float_t *Zincard; |
| 253 |
Float_t Xoutcard[ncar]; |
Float_t *Xoutcard; |
| 254 |
Float_t Youtcard[ncar]; |
Float_t *Youtcard; |
| 255 |
Float_t Zoutcard[ncar]; |
Float_t *Zoutcard; |
| 256 |
Float_t Erelcard[ncar]; |
Float_t *Erelcard; |
| 257 |
Float_t Timecard[ncar]; |
Float_t *Timecard; |
| 258 |
Float_t Pathcard[ncar]; |
Float_t *Pathcard; |
| 259 |
Float_t P0card[ncar]; |
Float_t *P0card; |
| 260 |
|
|
| 261 |
// List of branches |
// List of branches |
| 262 |
TBranch *b_Irun; //! |
TBranch *b_Irun; //! |
| 379 |
TBranch *b_Pathcard; //! |
TBranch *b_Pathcard; //! |
| 380 |
TBranch *b_P0card; //! |
TBranch *b_P0card; //! |
| 381 |
// |
// |
| 382 |
|
TF1 *attenAC; |
| 383 |
public: |
public: |
| 384 |
|
|
| 385 |
Digitizer(TTree *tree, char* &file_raw); // constructor |
Digitizer(TTree*,char*&,int,int,int,int,int,int,int,int); // constructor |
| 386 |
|
|
| 387 |
void Close(); |
void Close(); |
| 388 |
void SetTree(TTree *tree) {fhBookTree = tree;}; |
void SetTree(TTree *tree) {fhBookTree = tree;}; |
| 390 |
|
|
| 391 |
void DigitizePSCU(UInt_t, UChar_t, UShort_t*); // digitizes the CPU header |
void DigitizePSCU(UInt_t, UChar_t, UShort_t*); // digitizes the CPU header |
| 392 |
void AddPadding(); |
void AddPadding(); |
|
void DigitizeRunHeader(); |
|
|
void DigitizeRunTrailer(); |
|
|
// void DigitizeTRIGGER(); // SO nov'07: included in DigitizeTOF |
|
| 393 |
void DigitizeTOF(int np,float *atte1,float *atte2,float *lambda1,float *lambda2); // digitizes the tof |
void DigitizeTOF(int np,float *atte1,float *atte2,float *lambda1,float *lambda2); // digitizes the tof |
| 394 |
void LoadTOFCalib(int np,float *atte1,float *atte2,float *lambda1,float *lambda2); //load TOF calibration |
void LoadTOFCalib(int np,float *atte1,float *atte2,float *lambda1,float *lambda2); //load TOF calibration |
| 395 |
void DigitizeAC(); // digitizes the anti-counter |
void DigitizeAC(); // digitizes the anti-counter |
| 433 |
Float_t SaturationTracky(Float_t ADC); // Saturation simulation |
Float_t SaturationTracky(Float_t ADC); // Saturation simulation |
| 434 |
// end methods for Tracker |
// end methods for Tracker |
| 435 |
|
|
|
void WriteRunHeader(); |
|
|
void WriteRunTrailer(); |
|
| 436 |
void WriteData(); // writes data to file |
void WriteData(); // writes data to file |
| 437 |
void ReadData(); |
void ReadData(); |
| 438 |
|
|