]>
Commit | Line | Data |
---|---|---|
f3ca8792 | 1 | /* |
2 | ||
039a08c3 | 3 | This program reads the DAQ data files passed as argument using the monitoring library. |
4 | ||
039a08c3 | 5 | The program reports about its processing progress. |
f3ca8792 | 6 | |
7 | Messages on stdout are exported to DAQ log system. | |
8 | ||
b3b2b9dd | 9 | DA for ZDC standalone pedestal runs |
ea628de7 | 10 | |
b3b2b9dd | 11 | Contact: Chiara.Oppedisano@to.infn.it |
442e1b18 | 12 | Link: |
b3b2b9dd | 13 | Run Type: STANDALONE_PEDESTAL_RUN |
442e1b18 | 14 | DA Type: LDC |
b3b2b9dd | 15 | Number of events needed: no constraint (tipically ~10^3) |
27afc0c8 | 16 | Input Files: none |
442e1b18 | 17 | Output Files: ZDCPedestal.dat, ZDCChMapping.dat |
b3b2b9dd | 18 | Trigger Types Used: Standalone Trigger |
f3ca8792 | 19 | |
20 | */ | |
218f916a | 21 | #define PEDDATA_FILE "ZDCPedestal.dat" |
65448375 | 22 | #define PEDHISTO_FILE "ZDCPedHisto.root" |
218f916a | 23 | #define MAPDATA_FILE "ZDCChMapping.dat" |
f3ca8792 | 24 | |
25 | #include <stdio.h> | |
039a08c3 | 26 | #include <stdlib.h> |
f3ca8792 | 27 | #include <Riostream.h> |
28 | ||
29 | // DATE | |
f3ca8792 | 30 | #include <event.h> |
31 | #include <monitor.h> | |
039a08c3 | 32 | #include <daqDA.h> |
f3ca8792 | 33 | |
34 | //ROOT | |
65448375 | 35 | #include <TROOT.h> |
36 | #include <TPluginManager.h> | |
f3ca8792 | 37 | #include <TH1F.h> |
38 | #include <TH2F.h> | |
39 | #include <TProfile.h> | |
40 | #include <TF1.h> | |
41 | #include <TFile.h> | |
442e1b18 | 42 | #include <TFitter.h> |
65448375 | 43 | #include "TMinuitMinimizer.h" |
f3ca8792 | 44 | |
45 | //AliRoot | |
46 | #include <AliRawReaderDate.h> | |
442e1b18 | 47 | #include <AliRawEventHeaderBase.h> |
f3ca8792 | 48 | #include <AliZDCRawStream.h> |
49 | ||
50 | ||
51 | /* Main routine | |
039a08c3 | 52 | Arguments: list of DATE raw data files |
f3ca8792 | 53 | */ |
54 | int main(int argc, char **argv) { | |
442e1b18 | 55 | |
65448375 | 56 | |
57 | gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", | |
58 | "*", | |
59 | "TStreamerInfo", | |
60 | "RIO", | |
61 | "TStreamerInfo()"); | |
62 | ||
63 | TMinuitMinimizer m; | |
64 | gROOT->GetPluginManager()->AddHandler("ROOT::Math::Minimizer", "Minuit","TMinuitMinimizer", | |
65 | "Minuit", "TMinuitMinimizer(const char *)"); | |
66 | TVirtualFitter::SetDefaultFitter("Minuit"); | |
f3ca8792 | 67 | |
039a08c3 | 68 | int status = 0; |
9d0b658a | 69 | int const kNChannels = 24; |
c3986754 | 70 | int const kNScChannels = 32; |
8700404e | 71 | Int_t kFirstADCGeo=0, kLastADCGeo=3; |
c3986754 | 72 | |
73 | Int_t ich=0; | |
74 | Int_t adcMod[2*kNChannels], adcCh[2*kNChannels], sigCode[2*kNChannels]; | |
75 | Int_t det[2*kNChannels], sec[2*kNChannels]; | |
76 | for(Int_t y=0; y<2*kNChannels; y++){ | |
77 | adcMod[y]=adcCh[y]=sigCode[y]=det[y]=sec[y]=0; | |
78 | } | |
039a08c3 | 79 | |
80 | /* log start of process */ | |
63604324 | 81 | printf("\n ZDC PEDESTAL program started\n"); |
039a08c3 | 82 | |
83 | /* check that we got some arguments = list of files */ | |
65448375 | 84 | if (argc<2){ |
039a08c3 | 85 | printf("Wrong number of arguments\n"); |
86 | return -1; | |
87 | } | |
88 | ||
f3ca8792 | 89 | // --- Histograms for ADC pedestals |
b3b2b9dd | 90 | // [22 signal channels + 2 reference PTMs] x 2 gain chains |
f3ca8792 | 91 | // |
039a08c3 | 92 | TH1F::AddDirectory(0); |
9d0b658a | 93 | // |
b3b2b9dd | 94 | TH1F *hPedhg[kNChannels], *hPedOutOfTimehg[kNChannels]; |
95 | TH2F *hPedCorrhg[kNChannels]; | |
96 | TH1F *hPedlg[kNChannels], *hPedOutOfTimelg[kNChannels]; | |
97 | TH2F *hPedCorrlg[kNChannels]; | |
f3ca8792 | 98 | // |
b3b2b9dd | 99 | char namhist1hg[50], namhist2hg[50], namhist3hg[50]; |
100 | char namhist1lg[50], namhist2lg[50], namhist3lg[50]; | |
101 | for(Int_t j=0; j<kNChannels; j++){ | |
7f4bde92 | 102 | if(j<=4){ // ZNC |
103 | sprintf(namhist1hg,"PedZNChg_%d",j); | |
104 | sprintf(namhist2hg,"PedZNChgOutOfTime_%d",j); | |
105 | sprintf(namhist3hg,"PedCorrZNChg_%d",j); | |
b3b2b9dd | 106 | // |
7f4bde92 | 107 | sprintf(namhist1lg,"PedZNClg_%d",j); |
108 | sprintf(namhist2lg,"PedZNClgOutOfTime_%d",j); | |
109 | sprintf(namhist3lg,"PedCorrZNClg_%d",j); | |
b3b2b9dd | 110 | } |
7f4bde92 | 111 | else if(j>=5 && j<=9){ // ZPC |
112 | sprintf(namhist1hg,"PedZPChg_%d",j-5); | |
113 | sprintf(namhist2hg,"PedZPChgOutOfTime_%d",j-5); | |
114 | sprintf(namhist3hg,"PedCorrZPChg_%d",j-5); | |
b3b2b9dd | 115 | // |
7f4bde92 | 116 | sprintf(namhist1lg,"PedZPClg_%d",j-5); |
117 | sprintf(namhist2lg,"PedZPClgOutOfTime_%d",j-5); | |
118 | sprintf(namhist3lg,"PedCorrZPClg_%d",j-5); | |
f3ca8792 | 119 | } |
7f4bde92 | 120 | else if(j==10 || j==11){ // ZEM |
3018cd7f | 121 | sprintf(namhist1hg,"PedZEMhg_%d",j-9); |
122 | sprintf(namhist2hg,"PedZEMhgOutOfTime_%d",j-9); | |
123 | sprintf(namhist3hg,"PedCorrZEMhg_%d",j-9); | |
b3b2b9dd | 124 | // |
3018cd7f | 125 | sprintf(namhist1lg,"PedZEMlg_%d",j-9); |
126 | sprintf(namhist2lg,"PedZEMlgOutOfTime_%d",j-9); | |
127 | sprintf(namhist3lg,"PedCorrZEMlg_%d",j-9); | |
f3ca8792 | 128 | } |
7f4bde92 | 129 | else if(j>=12 && j<=16){ // ZNA |
130 | sprintf(namhist1hg,"PedZNAhg_%d",j-12); | |
131 | sprintf(namhist2hg,"PedZNAhgOutOfTime_%d",j-12); | |
132 | sprintf(namhist3hg,"PedCorrZNAhg_%d",j-12); | |
b3b2b9dd | 133 | // |
7f4bde92 | 134 | sprintf(namhist1lg,"PedZNAlg_%d",j-12); |
135 | sprintf(namhist2lg,"PedZNAlgOutOfTime_%d",j-12); | |
136 | sprintf(namhist3lg,"PedCorrZNAlg_%d",j-12); | |
f3ca8792 | 137 | } |
7f4bde92 | 138 | else if(j>=17 && j<=21){ // ZPA |
139 | sprintf(namhist1hg,"PedZPAhg_%d",j-17); | |
140 | sprintf(namhist2hg,"PedZPAhgOutOfTime_%d",j-17); | |
141 | sprintf(namhist3hg,"PedCorrZPAhg_%d",j-17); | |
b3b2b9dd | 142 | // |
7f4bde92 | 143 | sprintf(namhist1lg,"PedZPAlg_%d",j-17); |
144 | sprintf(namhist2lg,"PedZPAlgOutOfTime_%d",j-17); | |
145 | sprintf(namhist3lg,"PedCorrZPAlg_%d",j-17); | |
f3ca8792 | 146 | } |
07b2ecec | 147 | else if(j==22 || j==23){ //Reference PMs |
b3b2b9dd | 148 | sprintf(namhist1hg,"PedRefhg_%d",j-22); |
149 | sprintf(namhist2hg,"PedRefhgOutOfTime_%d",j-22); | |
150 | sprintf(namhist3hg,"PedCorrRefhg_%d",j-22); | |
151 | // | |
152 | sprintf(namhist1lg,"PedReflg_%d",j-22); | |
153 | sprintf(namhist2lg,"PedReflgOutOfTime_%d",j-22); | |
154 | sprintf(namhist3lg,"PedCorrReflg_%d",j-22); | |
f3ca8792 | 155 | } |
b3b2b9dd | 156 | // --- High gain chain histos |
057ed567 | 157 | hPedhg[j] = new TH1F(namhist1hg, namhist1hg, 200, 0., 200.); |
158 | hPedOutOfTimehg[j] = new TH1F(namhist2hg, namhist2hg, 200, 0., 200.); | |
b3b2b9dd | 159 | hPedCorrhg[j] = new TH2F(namhist3hg,namhist3hg,100,0.,200.,100,0.,200.); |
160 | // --- Low gain chain histos | |
057ed567 | 161 | hPedlg[j] = new TH1F(namhist1lg, namhist1lg, 100, 0., 1000.); |
162 | hPedOutOfTimelg[j] = new TH1F(namhist2lg, namhist2lg, 100, 0., 1000.); | |
163 | hPedCorrlg[j] = new TH2F(namhist3lg,namhist3lg,100,0.,1000.,100,0.,1000.); | |
f3ca8792 | 164 | } |
165 | ||
f3ca8792 | 166 | |
167 | /* open result file */ | |
168 | FILE *fp=NULL; | |
9d0b658a | 169 | fp=fopen("./result.txt","w"); |
f3ca8792 | 170 | if (fp==NULL) { |
171 | printf("Failed to open file\n"); | |
172 | return -1; | |
173 | } | |
174 | ||
442e1b18 | 175 | FILE *mapFile4Shuttle; |
f3ca8792 | 176 | |
039a08c3 | 177 | /* report progress */ |
178 | daqDA_progressReport(10); | |
f3ca8792 | 179 | |
f3ca8792 | 180 | |
f3ca8792 | 181 | /* init some counters */ |
182 | int nevents_physics=0; | |
183 | int nevents_total=0; | |
184 | ||
27afc0c8 | 185 | struct eventHeaderStruct *event; |
186 | eventTypeType eventT; | |
187 | ||
039a08c3 | 188 | /* read the data files */ |
189 | int n; | |
442e1b18 | 190 | for(n=1;n<argc;n++){ |
039a08c3 | 191 | |
192 | status=monitorSetDataSource( argv[n] ); | |
f3ca8792 | 193 | if (status!=0) { |
039a08c3 | 194 | printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status)); |
195 | return -1; | |
f3ca8792 | 196 | } |
197 | ||
039a08c3 | 198 | /* report progress */ |
199 | /* in this example, indexed on the number of files */ | |
200 | daqDA_progressReport(10+80*n/argc); | |
201 | ||
202 | /* read the file */ | |
203 | for(;;) { | |
039a08c3 | 204 | |
205 | /* get next event */ | |
206 | status=monitorGetEventDynamic((void **)&event); | |
442e1b18 | 207 | if(status==MON_ERR_EOF) break; /* end of monitoring file has been reached */ |
208 | if(status!=0) { | |
039a08c3 | 209 | printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status)); |
210 | return -1; | |
211 | } | |
212 | ||
213 | /* retry if got no event */ | |
442e1b18 | 214 | if(event==NULL) { |
039a08c3 | 215 | break; |
216 | } | |
442e1b18 | 217 | |
218 | // Initalize raw-data reading and decoding | |
219 | AliRawReader *reader = new AliRawReaderDate((void*)event); | |
220 | reader->Select("ZDC"); | |
221 | // --- Reading event header | |
222 | //UInt_t evtype = reader->GetType(); | |
223 | //printf("\n\t ZDCPEDESTALda -> ev. type %d\n",evtype); | |
224 | //printf("\t ZDCPEDESTALda -> run # %d\n",reader->GetRunNumber()); | |
225 | // | |
226 | AliZDCRawStream *rawStreamZDC = new AliZDCRawStream(reader); | |
227 | ||
f3ca8792 | 228 | |
039a08c3 | 229 | /* use event - here, just write event id to result file */ |
230 | eventT=event->eventType; | |
442e1b18 | 231 | |
c3986754 | 232 | Int_t iScCh=0; |
27afc0c8 | 233 | Int_t scMod[kNScChannels], scCh[kNScChannels], scSigCode[kNScChannels]; |
234 | Int_t scDet[kNScChannels], scSec[kNScChannels]; | |
235 | for(Int_t y=0; y<kNScChannels; y++){ | |
c3986754 | 236 | scMod[y]=scCh[y]=scSigCode[y]=scDet[y]=scSec[y]=0; |
237 | } | |
238 | // | |
239 | Int_t modNum=-1, modType=-1; | |
7f4bde92 | 240 | |
442e1b18 | 241 | if(eventT==START_OF_DATA){ |
242 | ||
7f4bde92 | 243 | rawStreamZDC->SetSODReading(kTRUE); |
244 | ||
9d0b658a | 245 | // -------------------------------------------------------- |
246 | // --- Writing ascii data file for the Shuttle preprocessor | |
247 | mapFile4Shuttle = fopen(MAPDATA_FILE,"w"); | |
442e1b18 | 248 | if(!rawStreamZDC->Next()) printf(" \t No raw data found!! \n"); |
249 | else{ | |
c3986754 | 250 | while((rawStreamZDC->Next())){ |
251 | if(rawStreamZDC->IsHeaderMapping()){ // mapping header | |
252 | modNum = rawStreamZDC->GetADCModule(); | |
253 | modType = rawStreamZDC->GetModType(); | |
254 | } | |
255 | if(rawStreamZDC->IsChMapping()){ | |
256 | if(modType==1){ // ADC mapping ---------------------- | |
257 | adcMod[ich] = rawStreamZDC->GetADCModFromMap(ich); | |
258 | adcCh[ich] = rawStreamZDC->GetADCChFromMap(ich); | |
259 | sigCode[ich] = rawStreamZDC->GetADCSignFromMap(ich); | |
260 | det[ich] = rawStreamZDC->GetDetectorFromMap(ich); | |
261 | sec[ich] = rawStreamZDC->GetTowerFromMap(ich); | |
262 | // | |
263 | fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\t%d\t%d\n", | |
264 | ich,adcMod[ich],adcCh[ich],sigCode[ich],det[ich],sec[ich]); | |
265 | // | |
266 | //printf(" Mapping in DA -> %d ADC: mod %d ch %d, code %d det %d, sec %d\n", | |
267 | // ich,adcMod[ich],adcCh[ich],sigCode[ich],det[ich],sec[ich]); | |
268 | // | |
269 | ich++; | |
270 | } | |
271 | else if(modType==2){ //VME scaler mapping -------------------- | |
272 | scMod[iScCh] = rawStreamZDC->GetScalerModFromMap(iScCh); | |
273 | scCh[iScCh] = rawStreamZDC->GetScalerChFromMap(iScCh); | |
274 | scSigCode[iScCh] = rawStreamZDC->GetScalerSignFromMap(iScCh); | |
275 | scDet[iScCh] = rawStreamZDC->GetScDetectorFromMap(iScCh); | |
276 | scSec[iScCh] = rawStreamZDC->GetScTowerFromMap(iScCh); | |
277 | // | |
278 | fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\t%d\t%d\n", | |
279 | iScCh,scMod[iScCh],scCh[iScCh],scSigCode[iScCh],scDet[iScCh],scSec[iScCh]); | |
280 | // | |
281 | //printf(" Mapping in DA -> %d Scaler: mod %d ch %d, code %d det %d, sec %d\n", | |
282 | // iScCh,scMod[iScCh],scCh[iScCh],scSigCode[iScCh],scDet[iScCh],scSec[iScCh]); | |
283 | // | |
284 | iScCh++; | |
285 | } | |
442e1b18 | 286 | } |
287 | } | |
288 | } | |
442e1b18 | 289 | fclose(mapFile4Shuttle); |
9d0b658a | 290 | }// SOD event |
442e1b18 | 291 | |
d3f256ca | 292 | else if(eventT==PHYSICS_EVENT){ |
442e1b18 | 293 | // --- Reading data header |
294 | reader->ReadHeader(); | |
039a08c3 | 295 | const AliRawDataHeader* header = reader->GetDataHeader(); |
442e1b18 | 296 | if(header){ |
b3b2b9dd | 297 | UChar_t message = header->GetAttributes(); |
a5747d8f | 298 | if((message & 0x20) == 0x20){ // PEDESTAL RUN |
442e1b18 | 299 | //printf("\t STANDALONE_PEDESTAL RUN raw data found\n"); |
300 | } | |
301 | else{ | |
27afc0c8 | 302 | printf("ZDCPEDESTALda.cxx -> NO STANDALONE_PEDESTAL RUN raw data found\n"); |
8700404e | 303 | printf(" CDH attributes: %x\n", message); |
27afc0c8 | 304 | return -1; |
442e1b18 | 305 | } |
306 | } | |
307 | else{ | |
308 | printf("\t ATTENTION! No Raw Data Header found!!!\n"); | |
309 | return -1; | |
310 | } | |
7f4bde92 | 311 | |
312 | rawStreamZDC->SetSODReading(kTRUE); | |
442e1b18 | 313 | |
314 | if(!rawStreamZDC->Next()) printf(" \t No raw data found!! \n"); | |
315 | // | |
316 | // ----- Setting ch. mapping ----- | |
9d0b658a | 317 | for(Int_t jk=0; jk<2*kNChannels; jk++){ |
7f4bde92 | 318 | //printf("ZDCPEDESTALDA.cxx -> ch.%d mod %d, ch %d, code %d det %d, sec %d\n", |
319 | // jk,adcMod[jk],adcCh[jk],sigCode[jk],det[jk],sec[jk]); | |
442e1b18 | 320 | rawStreamZDC->SetMapADCMod(jk, adcMod[jk]); |
321 | rawStreamZDC->SetMapADCCh(jk, adcCh[jk]); | |
322 | rawStreamZDC->SetMapADCSig(jk, sigCode[jk]); | |
323 | rawStreamZDC->SetMapDet(jk, det[jk]); | |
324 | rawStreamZDC->SetMapTow(jk, sec[jk]); | |
325 | } | |
326 | // | |
327 | Int_t iraw=0; | |
039a08c3 | 328 | Int_t RawADChg[kNChannels], RawADCoothg[kNChannels]; |
329 | Int_t RawADClg[kNChannels], RawADCootlg[kNChannels]; | |
330 | for(Int_t j=0; j<kNChannels; j++){ | |
331 | RawADChg[j]=0; RawADCoothg[j]=0; | |
332 | RawADClg[j]=0; RawADCootlg[j]=0; | |
333 | } | |
442e1b18 | 334 | // |
039a08c3 | 335 | while(rawStreamZDC->Next()){ |
442e1b18 | 336 | Int_t index=-1; |
7f4bde92 | 337 | Int_t detector = rawStreamZDC->GetSector(0); |
9d0b658a | 338 | Int_t sector = rawStreamZDC->GetSector(1); |
a5747d8f | 339 | // |
340 | //printf(" rawData: det %d sec %d value %d\n", | |
341 | // detector, sector,rawStreamZDC->GetADCGain(),rawStreamZDC->GetADCValue() ); | |
9d0b658a | 342 | |
8700404e | 343 | if((rawStreamZDC->IsADCDataWord()) && (detector!=-1) && |
344 | (rawStreamZDC->GetADCModule()>=kFirstADCGeo && rawStreamZDC->GetADCModule()<=kLastADCGeo)){ | |
9d0b658a | 345 | if(sector!=5){ // Physics signals |
346 | if(detector==1) index = sector; // *** ZNC | |
347 | else if(detector==2) index = sector+5; // *** ZPC | |
348 | else if(detector==3) index = sector+9; // *** ZEM | |
349 | else if(detector==4) index = sector+12; // *** ZNA | |
350 | else if(detector==5) index = sector+17; // *** ZPA | |
f3ca8792 | 351 | } |
b3b2b9dd | 352 | else{ // Reference PMs |
7f4bde92 | 353 | index = (detector-1)/3+22; |
f3ca8792 | 354 | } |
b3b2b9dd | 355 | // |
c3986754 | 356 | if(index==-1) printf("ERROR in ZDCPEDESTALda.cxx -> det %d quad %d index %d \n", |
357 | detector,sector,index); | |
7f4bde92 | 358 | |
442e1b18 | 359 | // |
360 | if(iraw<2*kNChannels){ // --- In-time pedestals (1st 48 raw data) | |
b3b2b9dd | 361 | if(rawStreamZDC->GetADCGain()==0){ |
362 | hPedhg[index]->Fill(rawStreamZDC->GetADCValue()); | |
363 | RawADChg[index] = rawStreamZDC->GetADCValue(); | |
442e1b18 | 364 | // |
365 | //printf("\t filling histo hPedhg[%d]\n",index); | |
b3b2b9dd | 366 | } |
367 | else{ | |
368 | hPedlg[index]->Fill(rawStreamZDC->GetADCValue()); | |
369 | RawADClg[index] = rawStreamZDC->GetADCValue(); | |
442e1b18 | 370 | // |
371 | //printf("\t filling histo hPedlg[%d]\n",index); | |
b3b2b9dd | 372 | } |
442e1b18 | 373 | } |
374 | else{ // --- Out-of-time pedestals | |
375 | if(rawStreamZDC->GetADCGain()==0){ | |
376 | hPedOutOfTimehg[index]->Fill(rawStreamZDC->GetADCValue()); | |
377 | RawADCoothg[index] = rawStreamZDC->GetADCValue(); | |
378 | // | |
379 | //printf("\t filling histo hPedOutOfTimehg[%d]\n",index); | |
039a08c3 | 380 | } |
442e1b18 | 381 | else{ |
382 | hPedOutOfTimelg[index]->Fill(rawStreamZDC->GetADCValue()); | |
383 | RawADCootlg[index] = rawStreamZDC->GetADCValue(); | |
384 | // | |
385 | //printf("\t filling histo hPedOutOfTimelg[%d]\n",index); | |
039a08c3 | 386 | } |
442e1b18 | 387 | } |
388 | iraw++; | |
389 | }//IsADCDataWord() | |
d3f256ca | 390 | } |
5aa49948 | 391 | // |
392 | if(iraw==4*kNChannels ){ // Last ADC channel -> Filling correlation histos | |
393 | for(Int_t k=0; k<kNChannels; k++){ | |
394 | hPedCorrhg[k]->Fill(RawADCoothg[k], RawADChg[k]); | |
395 | hPedCorrlg[k]->Fill(RawADCootlg[k], RawADClg[k]); | |
396 | } | |
397 | //printf(" ev. %d -> Filling correlation histo %d\n",nevents_physics, kNChannels); | |
398 | } | |
d3f256ca | 399 | nevents_physics++; |
400 | // | |
401 | delete reader; | |
402 | delete rawStreamZDC; | |
039a08c3 | 403 | |
404 | }//(if PHYSICS_EVENT) | |
d3f256ca | 405 | |
406 | /* exit when last event received, no need to wait for TERM signal */ | |
407 | else if(eventT==END_OF_RUN) { | |
408 | printf(" -> EOR event detected\n"); | |
409 | break; | |
410 | } | |
411 | ||
039a08c3 | 412 | nevents_total++; |
f3ca8792 | 413 | |
f3ca8792 | 414 | } |
d3f256ca | 415 | |
416 | /* free resources */ | |
417 | free(event); | |
039a08c3 | 418 | } |
f3ca8792 | 419 | |
420 | /* Analysis of the histograms */ | |
421 | // | |
422 | FILE *fileShuttle; | |
218f916a | 423 | fileShuttle = fopen(PEDDATA_FILE,"w"); |
f3ca8792 | 424 | // |
b3b2b9dd | 425 | Float_t MeanPed[2*kNChannels], MeanPedWidth[2*kNChannels], |
63604324 | 426 | MeanPedOOT[2*kNChannels], MeanPedWidthOOT[2*kNChannels]; |
b3b2b9dd | 427 | // --- In-time pedestals |
428 | TF1 *ADCfunchg[kNChannels]; | |
429 | for(Int_t i=0; i<kNChannels; i++){ | |
f9641c3b | 430 | if(hPedhg[i]->GetEntries() == 0){ |
431 | printf("\n WARNING! Empty histos -> ending DA WITHOUT writing output\n\n"); | |
432 | return -1; | |
433 | } | |
b3b2b9dd | 434 | hPedhg[i]->Fit("gaus","Q"); |
435 | ADCfunchg[i] = hPedhg[i]->GetFunction("gaus"); | |
442e1b18 | 436 | MeanPed[i] = (Double_t) ADCfunchg[i]->GetParameter(1); |
437 | MeanPedWidth[i] = (Double_t) ADCfunchg[i]->GetParameter(2); | |
218f916a | 438 | fprintf(fileShuttle,"\t%f\t%f\n",MeanPed[i],MeanPedWidth[i]); |
057ed567 | 439 | //printf("\t MeanPedhg[%d] = %f\n",i, MeanPed[i]); |
b3b2b9dd | 440 | } |
441 | TF1 *ADCfunclg[kNChannels]; | |
442 | for(Int_t i=0; i<kNChannels; i++){ | |
443 | hPedlg[i]->Fit("gaus","Q"); | |
444 | ADCfunclg[i] = hPedlg[i]->GetFunction("gaus"); | |
442e1b18 | 445 | MeanPed[i+kNChannels] = (Double_t) ADCfunclg[i]->GetParameter(1); |
446 | MeanPedWidth[i+kNChannels] = (Double_t) ADCfunclg[i]->GetParameter(2); | |
218f916a | 447 | fprintf(fileShuttle,"\t%f\t%f\n",MeanPed[i+kNChannels],MeanPedWidth[i+kNChannels]); |
057ed567 | 448 | //printf("\t MeanPedlg[%d] = %f\n",i+kNChannels, MeanPed[i+kNChannels]); |
f3ca8792 | 449 | } |
450 | // --- Out-of-time pedestals | |
b3b2b9dd | 451 | TF1 *ADCootfunchg[kNChannels]; |
452 | for(Int_t i=0; i<kNChannels; i++){ | |
3d155339 | 453 | if(hPedOutOfTimehg[i]->GetEntries() == 0){ |
a5747d8f | 454 | printf("\n WARNING! Empty histos for out-of-time channels!!!\n\n"); |
3d155339 | 455 | return -1; |
456 | } | |
b3b2b9dd | 457 | hPedOutOfTimehg[i]->Fit("gaus","Q"); |
458 | ADCootfunchg[i] = hPedOutOfTimehg[i]->GetFunction("gaus"); | |
442e1b18 | 459 | MeanPedOOT[i] = (Double_t) ADCootfunchg[i]->GetParameter(1); |
460 | MeanPedWidthOOT[i] = (Double_t) ADCootfunchg[i]->GetParameter(2); | |
218f916a | 461 | fprintf(fileShuttle,"\t%f\t%f\n",MeanPedOOT[i],MeanPedWidthOOT[i]); |
057ed567 | 462 | //printf("\t MeanPedOOThg[%d] = %f\n",i, MeanPedOOT[i]); |
b3b2b9dd | 463 | } |
464 | TF1 *ADCootfunclg[kNChannels]; | |
465 | for(Int_t i=0; i<kNChannels; i++){ | |
466 | hPedOutOfTimelg[i]->Fit("gaus","Q"); | |
467 | ADCootfunclg[i] = hPedOutOfTimelg[i]->GetFunction("gaus"); | |
442e1b18 | 468 | MeanPedOOT[i+kNChannels] = (Double_t) ADCootfunclg[i]->GetParameter(1); |
469 | MeanPedWidthOOT[i+kNChannels] = (Double_t) ADCootfunclg[i]->GetParameter(2); | |
218f916a | 470 | fprintf(fileShuttle,"\t%f\t%f\n",MeanPedOOT[i+kNChannels],MeanPedWidthOOT[i+kNChannels]); |
057ed567 | 471 | //printf("\t MeanPedOOTlg[%d] = %f\n",i+kNChannels, MeanPedOOT[i+kNChannels]); |
f3ca8792 | 472 | } |
63604324 | 473 | |
b3b2b9dd | 474 | // --- Correlations |
65bae83c | 475 | // NB -> The correlations are NOT fitted since at the moment |
476 | // (Sptember 2009) they are NOT correlations and the DA would fail!!! | |
477 | /* Float_t CorrCoeff0[2*kNChannels], CorrCoeff1[2*kNChannels]; | |
b3b2b9dd | 478 | TProfile *hPedCorrProfhg[kNChannels], *hPedCorrProflg[kNChannels]; |
479 | TF1 *ffunchg[kNChannels], *ffunclg[kNChannels]; | |
f3ca8792 | 480 | char namhist4[50]; |
b3b2b9dd | 481 | for(int i=0;i<kNChannels;i++) { |
482 | sprintf(namhist4,"ADCHRvsOOT%d_Prof",i); | |
483 | hPedCorrProfhg[i] = hPedCorrhg[i]->ProfileX(namhist4,-1,-1,"S"); | |
484 | hPedCorrProfhg[i]->SetName(namhist4); | |
485 | hPedCorrProfhg[i]->Fit("pol1","Q"); | |
486 | ffunchg[i] = hPedCorrProfhg[i]->GetFunction("pol1"); | |
442e1b18 | 487 | CorrCoeff0[i] = (Double_t) ffunchg[i]->GetParameter(0); |
488 | CorrCoeff1[i] = (Double_t) ffunchg[i]->GetParameter(1); | |
9d0b658a | 489 | fprintf(fileShuttle,"\t%f\t%f\n",CorrCoeff0[i],CorrCoeff1[i]); |
ea628de7 | 490 | //printf("\t CorrCoeff0[%d] = %f, CorrCoeff1[%d] = %f\n",i, CorrCoeff0[i], i, CorrCoeff1[i]); |
b3b2b9dd | 491 | } |
492 | for(int i=0;i<kNChannels;i++) { | |
493 | sprintf(namhist4,"ADCLRvsOOT%d_Prof",i); | |
494 | hPedCorrProflg[i] = hPedCorrlg[i]->ProfileX(namhist4,-1,-1,"S"); | |
495 | hPedCorrProflg[i]->SetName(namhist4); | |
496 | hPedCorrProflg[i]->Fit("pol1","Q"); | |
497 | ffunclg[i] = hPedCorrProflg[i]->GetFunction("pol1"); | |
442e1b18 | 498 | CorrCoeff0[i+kNChannels] = (Double_t) ffunclg[i]->GetParameter(0); |
499 | CorrCoeff1[i+kNChannels] = (Double_t) ffunclg[i]->GetParameter(1); | |
9d0b658a | 500 | fprintf(fileShuttle,"\t%f\t%f\n",CorrCoeff0[i+kNChannels],CorrCoeff1[i+kNChannels]); |
ea628de7 | 501 | //printf("\t CorrCoeff0[%d] = %f, CorrCoeff1[%d] = %f\n", |
502 | // i+kNChannels, CorrCoeff0[i+kNChannels], i+kNChannels, CorrCoeff1[i+kNChannels]); | |
f3ca8792 | 503 | } |
65bae83c | 504 | */ |
f3ca8792 | 505 | // |
506 | fclose(fileShuttle); | |
b3b2b9dd | 507 | // |
65448375 | 508 | /* report progress */ |
509 | daqDA_progressReport(80); | |
510 | // | |
511 | TFile *histofile = new TFile(PEDHISTO_FILE,"RECREATE"); | |
512 | histofile->cd(); | |
513 | for(int k=0; k<kNChannels; k++){ | |
514 | hPedhg[k]->Write(); | |
515 | hPedOutOfTimehg[k]->Write(); | |
516 | hPedCorrhg[k]->Write(); | |
517 | hPedlg[k]->Write(); | |
518 | hPedOutOfTimelg[k]->Write(); | |
519 | hPedCorrlg[k]->Write(); | |
520 | } | |
521 | // | |
522 | histofile->Close(); | |
523 | // | |
039a08c3 | 524 | for(Int_t j=0; j<kNChannels; j++){ |
525 | delete hPedhg[j]; | |
526 | delete hPedOutOfTimehg[j]; | |
527 | delete hPedCorrhg[j]; | |
528 | delete hPedlg[j]; | |
529 | delete hPedOutOfTimelg[j]; | |
530 | delete hPedCorrlg[j]; | |
531 | } | |
f3ca8792 | 532 | |
533 | /* write report */ | |
534 | fprintf(fp,"Run #%s, received %d physics events out of %d\n",getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total); | |
535 | ||
536 | /* close result file */ | |
537 | fclose(fp); | |
039a08c3 | 538 | |
539 | /* report progress */ | |
540 | daqDA_progressReport(90); | |
541 | ||
442e1b18 | 542 | /* store the result files on FES */ |
63604324 | 543 | // [1] File with mapping |
a8eaff64 | 544 | status = daqDA_FES_storeFile(MAPDATA_FILE, "MAPPING"); |
442e1b18 | 545 | if(status){ |
65448375 | 546 | printf("Failed to export mapping data file to DAQ FES\n"); |
442e1b18 | 547 | return -1; |
548 | } | |
63604324 | 549 | // [2] File with pedestal data |
a8eaff64 | 550 | status = daqDA_FES_storeFile(PEDDATA_FILE, "PEDESTALDATA"); |
039a08c3 | 551 | if(status){ |
65448375 | 552 | printf("Failed to export pedestal data file to DAQ FES\n"); |
553 | return -1; | |
554 | } | |
555 | // [3] File with pedestal histos | |
a8eaff64 | 556 | status = daqDA_FES_storeFile(PEDHISTO_FILE, "PEDESTALHISTOS"); |
65448375 | 557 | if(status){ |
558 | printf("Failed to export pedestal histos file to DAQ FES\n"); | |
63604324 | 559 | return -1; |
560 | } | |
561 | ||
562 | /* store the result files on DB */ | |
a8eaff64 | 563 | status = daqDA_DB_storeFile(PEDDATA_FILE, PEDDATA_FILE); |
63604324 | 564 | if(status){ |
65448375 | 565 | printf("Failed to store pedestal data file to DAQ DB\n"); |
039a08c3 | 566 | return -1; |
567 | } | |
568 | ||
569 | /* report progress */ | |
570 | daqDA_progressReport(100); | |
f3ca8792 | 571 | |
f3ca8792 | 572 | return status; |
573 | } |