]>
Commit | Line | Data |
---|---|---|
9c754ce7 | 1 | /* |
2 | TPC DA for online calibration | |
3 | ||
4 | Contact: Haavard.Helstrup@cern.ch | |
5 | Link: | |
6 | Run Type: PULSER_RUN | |
7 | DA Type: LDC | |
8 | Number of events needed: 100 | |
9 | Input Files: | |
10 | Output Files: tpcPulser.root, to be exported to the DAQ FXS | |
11 | Trigger types used: CALIBRATION_EVENT | |
12 | ||
13 | */ | |
14 | ||
ec2624ea | 15 | /* |
16 | ||
17 | TPCda_pulser.cxx - calibration algorithm for TPC pulser events | |
18 | ||
19 | 10/06/2007 sylvain.chapeland@cern.ch : first version - clean skeleton based on DAQ DA case1 | |
20 | 30/09/2007 haavard.helstrup@cern.ch : created pulser DA based on pedestal code | |
21 | ||
22 | contact: marian.ivanov@cern.ch | |
23 | ||
24 | ||
25 | This process reads RAW data from the files provided as command line arguments | |
26 | and save results in a file (named from RESULT_FILE define - see below). | |
27 | ||
28 | */ | |
29 | ||
30 | #define RESULT_FILE "tpcPulser.root" | |
31 | ||
32 | ||
33 | #include <daqDA.h> | |
34 | #include "event.h" | |
35 | #include "monitor.h" | |
36 | #include <stdio.h> | |
37 | #include <stdlib.h> | |
38 | ||
39 | // | |
40 | //Root includes | |
41 | // | |
42 | #include <TFile.h> | |
43 | #include "TROOT.h" | |
44 | #include "TPluginManager.h" | |
45 | // | |
46 | //AliRoot includes | |
47 | // | |
48 | #include "AliRawReader.h" | |
49 | #include "AliRawReaderDate.h" | |
50 | #include "AliTPCRawStream.h" | |
51 | #include "AliTPCROC.h" | |
52 | #include "AliTPCCalROC.h" | |
53 | #include "AliTPCCalPad.h" | |
54 | #include "AliMathBase.h" | |
55 | #include "TTreeStream.h" | |
56 | ||
57 | // | |
58 | // TPC calibration algorithm includes | |
59 | // | |
60 | #include "AliTPCCalibPulser.h" | |
61 | ||
62 | ||
63 | ||
64 | ||
65 | /* Main routine | |
66 | Arguments: list of DATE raw data files | |
67 | */ | |
68 | int main(int argc, char **argv) { | |
69 | ||
70 | gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", | |
71 | "*", | |
72 | "TStreamerInfo", | |
73 | "RIO", | |
74 | "TStreamerInfo()"); | |
75 | ||
76 | ||
77 | int i,status; | |
78 | AliTPCCalibPulser calibPulser; // pedestal and noise calibration | |
79 | ||
80 | if (argc<2) { | |
81 | printf("Wrong number of arguments\n"); | |
82 | return -1; | |
83 | } | |
84 | ||
85 | ||
86 | /* log start of process */ | |
87 | printf("TPC Pulser DA started - %s\n",__FILE__); | |
88 | ||
ca273e23 | 89 | /* set time bin range */ |
86dd6fa4 | 90 | calibPulser.SetRangeTime(400,500); |
ec2624ea | 91 | |
92 | /* declare monitoring program */ | |
93 | status=monitorDeclareMp( __FILE__ ); | |
94 | if (status!=0) { | |
95 | printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status)); | |
96 | return -1; | |
97 | } | |
98 | ||
99 | ||
100 | /* loop over RAW data files */ | |
101 | int nevents=0; | |
102 | for(i=1;i<argc;i++) { | |
103 | ||
104 | /* define data source : this is argument i */ | |
105 | printf("Processing file %s\n", argv[i]); | |
106 | status=monitorSetDataSource( argv[i] ); | |
107 | if (status!=0) { | |
108 | printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status)); | |
109 | return -1; | |
110 | } | |
111 | ||
112 | /* read until EOF */ | |
113 | while (true) { | |
114 | struct eventHeaderStruct *event; | |
115 | ||
116 | /* check shutdown condition */ | |
117 | if (daqDA_checkShutdown()) {break;} | |
118 | ||
119 | /* get next event (blocking call until timeout) */ | |
120 | status=monitorGetEventDynamic((void **)&event); | |
121 | if (status==MON_ERR_EOF) { | |
122 | printf ("End of File %d detected\n",i); | |
123 | break; /* end of monitoring file has been reached */ | |
124 | } | |
125 | ||
126 | if (status!=0) { | |
127 | printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status)); | |
128 | break; | |
129 | } | |
130 | ||
131 | /* retry if got no event */ | |
132 | if (event==NULL) { | |
133 | continue; | |
134 | } | |
135 | nevents++; | |
136 | ||
137 | // Pulser calibration | |
138 | ||
139 | AliRawReader *rawReader = new AliRawReaderDate((void*)event); | |
86dd6fa4 | 140 | // calibPulser.ProcessEvent(rawReader); |
141 | calibPulser.ProcessEventFast(rawReader); | |
ec2624ea | 142 | delete rawReader; |
143 | ||
144 | /* free resources */ | |
145 | free(event); | |
146 | } | |
147 | } | |
148 | ||
149 | calibPulser.Analyse(); | |
150 | printf ("%d events processed\n",nevents); | |
151 | ||
152 | TFile * fileTPC = new TFile (RESULT_FILE,"recreate"); | |
153 | calibPulser.Write("calibPulser"); | |
154 | delete fileTPC; | |
155 | printf("Wrote %s\n",RESULT_FILE); | |
156 | ||
f2c72763 | 157 | /* store the result file on FES */ |
158 | ||
159 | status=daqDA_FES_storeFile(RESULT_FILE,RESULT_FILE); | |
160 | if (status) { | |
161 | status = -2; | |
162 | } | |
163 | ||
ec2624ea | 164 | return status; |
165 | } |