]>
Commit | Line | Data |
---|---|---|
0ec8015e | 1 | /********************************************************************************* |
2 | - Contact: Brigitte Cheynis b.cheynis@ipnl.in2p3.fr | |
3 | - Link: http | |
4 | - Raw data test file : | |
5 | - Reference run number : | |
6 | - Run Type: CHANNEL_DELAY_TUNING | |
7 | - DA Type: LDC | |
8 | - Number of events needed: 500 | |
9 | - Input Files: argument list | |
10 | - Output Files: local file V0_ChannelDelayTuning.dat | |
11 | FXS file V0_ChannelDelayTuning.dat | |
12 | - Trigger types used: PHYSICS_EVENT | |
13 | **********************************************************************************/ | |
14 | ||
15 | ||
16 | /********************************************************************************** | |
17 | * * | |
18 | * VZERO Detector Algorithm used for tuning FEE parameters * | |
19 | * * | |
20 | * This program reads data on the LDC * | |
21 | * It cumulates fBB and fBG flags, populates local "./V0_ChannelDelayTuning.dat" * | |
22 | * file, exports it to the FES, and stores it to DAQ DB * | |
23 | * We have 128 channels instead of 64 as expected for V0 due to the two sets of * | |
24 | * charge integrators which are used by the FEE ... * | |
25 | * The program reports about its processing progress. * | |
26 | * * | |
27 | ***********************************************************************************/ | |
28 | ||
29 | // DATE | |
30 | #include "event.h" | |
31 | #include "monitor.h" | |
32 | #include "daqDA.h" | |
33 | ||
34 | //AliRoot | |
35 | #include <AliVZERORawStream.h> | |
36 | #include <AliRawReaderDate.h> | |
37 | #include <AliRawReader.h> | |
38 | #include <AliDAQ.h> | |
39 | ||
40 | // standard | |
41 | #include <stdio.h> | |
42 | #include <stdlib.h> | |
43 | ||
44 | //ROOT | |
45 | #include "TROOT.h" | |
46 | #include "TPluginManager.h" | |
47 | #include <TFile.h> | |
48 | #include <TH1F.h> | |
49 | #include <TMath.h> | |
50 | ||
51 | ||
52 | /* Main routine --- Arguments: list of DATE raw data files */ | |
53 | ||
54 | int main(int argc, char **argv) { | |
55 | ||
56 | /* magic line from Cvetan */ | |
57 | gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", | |
58 | "*", | |
59 | "TStreamerInfo", | |
60 | "RIO", | |
61 | "TStreamerInfo()"); | |
62 | int status; | |
073bb639 | 63 | status = daqDA_DB_getFile("V00da_results","./V0_Pedestals.dat"); |
0ec8015e | 64 | if (status) { |
073bb639 | 65 | printf("Failed to get Pedestal file (V0_Pedestals.dat) from DAQ DB, status=%d\n", status); |
0ec8015e | 66 | return -1; |
67 | } | |
68 | ||
69 | Float_t MeanPed[128], SigPed[128], fdump; | |
70 | ||
71 | /* open the pedestal file and retrieve pedestal mean and sigma */ | |
073bb639 | 72 | FILE *fpPed = fopen("V0_Pedestals.dat","r"); |
0ec8015e | 73 | for(int i=0;i<128;i++){ |
74 | fscanf(fpPed,"%f %f %f %f \n",&MeanPed[i],&SigPed[i],&fdump,&fdump); | |
75 | // printf("%.3f %.3f \n",MeanPed[i],SigPed[i]); | |
76 | } | |
77 | fclose(fpPed); | |
78 | ||
79 | //______________________________________________________________________________ | |
80 | // Get running parameters from V00_ChannelDelayTuning_DA.config file | |
81 | ||
82 | Int_t kNbEventSkipped; // = 100; number of events skipped - to be tuned | |
83 | Float_t kSigmaCut; // = 3.0; number of sigmas for threshold cut - to be tuned | |
84 | ||
85 | status = daqDA_DB_getFile("V00_ChannelDelayTuning_DA.config","./V00_ChannelDelayTuning_DA.config"); | |
86 | if (status) { | |
87 | printf("Failed to get config file (V00_ChannelDelayTuning_DA.config) from DAQ DB, status=%d\n", status); | |
88 | return -1; | |
89 | } | |
90 | /* open the config file and retrieve running parameters */ | |
91 | FILE *fpConfig = fopen("V00_ChannelDelayTuning_DA.config","r"); | |
92 | fscanf(fpConfig,"%f %d",&kSigmaCut,&kNbEventSkipped); | |
93 | fclose(fpConfig); | |
94 | ||
95 | printf("Number of events skipped = %d ; Number of sigmas for threshold cut = %f\n",kNbEventSkipped,kSigmaCut); | |
96 | //______________________________________________________________________________ | |
97 | ||
073bb639 | 98 | Int_t BBFlag[64]; |
99 | Int_t BGFlag[64]; | |
100 | Int_t ChargeEoI = 0; | |
0ec8015e | 101 | Bool_t Integrator = 0; |
102 | Int_t NHit[64]; | |
103 | for(Int_t i=0; i<64; i++) { | |
104 | BBFlag[i] = 0; | |
105 | BGFlag[i] = 0; | |
073bb639 | 106 | NHit[i] = 0; |
0ec8015e | 107 | } |
108 | ||
109 | /* log start of process */ | |
110 | printf("VZERO DA program started - Channel Delay Tuning \n"); | |
111 | ||
112 | /* check that we got some arguments */ | |
113 | if (argc<2) { | |
114 | printf("Wrong number of arguments\n"); | |
115 | return -1;} | |
116 | ||
117 | /* open result file to be exported to FES */ | |
118 | FILE *fp=NULL; | |
119 | fp=fopen("./V0_ChannelDelayTuning.dat","a"); | |
120 | if (fp==NULL) { | |
121 | printf("Failed to open result file\n"); | |
122 | return -1;} | |
123 | ||
124 | /* open log file to inform user */ | |
125 | FILE *flog=NULL; | |
126 | flog=fopen("./V00log.txt","w"); | |
127 | if (flog==NULL) { | |
128 | printf("Failed to open log file\n"); | |
129 | return -1; } | |
130 | ||
131 | /* report progress */ | |
132 | daqDA_progressReport(10); | |
133 | ||
134 | /* init counters on events */ | |
135 | int nevents_physics=0; | |
136 | int nevents_total=0; | |
137 | int iteration; | |
138 | sscanf(argv[1],"%d",&iteration); | |
139 | ||
140 | /* read the n data files */ | |
141 | for (int n=2; n<argc; n++) { | |
142 | ||
143 | /* read the data */ | |
144 | status=monitorSetDataSource( argv[n] ); | |
145 | if (status!=0) { | |
146 | printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status)); | |
147 | return -1; } | |
148 | ||
149 | /* report progress */ | |
150 | daqDA_progressReport(10+50*n/argc); | |
151 | ||
152 | /* read the data file */ | |
153 | for(;;) { | |
154 | struct eventHeaderStruct *event; | |
155 | eventTypeType eventT; | |
156 | ||
157 | /* get next event */ | |
158 | status=monitorGetEventDynamic((void **)&event); | |
159 | if (status==MON_ERR_EOF) break; /* end of monitoring file has been reached */ | |
160 | if (status!=0) { | |
161 | printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status)); | |
162 | return -1; } | |
163 | ||
164 | /* retry if got no event */ | |
165 | if (event==NULL) break; | |
166 | ||
167 | /* decode event */ | |
168 | eventT=event->eventType; | |
169 | ||
170 | switch (event->eventType){ | |
171 | ||
172 | case START_OF_RUN: | |
173 | printf("START Of Run detected\n"); | |
174 | break; | |
175 | ||
176 | case END_OF_RUN: | |
177 | printf("END Of Run detected\n"); | |
178 | break; | |
179 | ||
180 | case PHYSICS_EVENT: | |
181 | nevents_physics++; | |
182 | ||
183 | AliRawReader *rawReader = new AliRawReaderDate((void*)event); | |
184 | ||
185 | AliVZERORawStream* rawStream = new AliVZERORawStream(rawReader); | |
186 | rawStream->Next(); | |
187 | if(nevents_physics > kNbEventSkipped){ | |
188 | for(Int_t i=0; i<64; i++) { | |
073bb639 | 189 | //ChargeEoI= rawStream->GetADC(i); takes EoI as central clock |
190 | // Look for the maximum in the LHC clock train instead of central clock 10 | |
191 | ChargeEoI = 0; | |
192 | Int_t iClock = 0; | |
193 | for(size_t iEvent=0; iEvent<21; iEvent++){ | |
194 | if(rawStream->GetPedestal(i,iEvent)>ChargeEoI) | |
195 | {ChargeEoI= rawStream->GetPedestal(i,iEvent); | |
196 | iClock = iEvent;} | |
197 | } | |
198 | Integrator = rawStream->GetIntegratorFlag(i,iClock); | |
0ec8015e | 199 | Float_t Threshold = MeanPed[i + 64 * Integrator] + kSigmaCut * SigPed[i + 64 * Integrator]; |
200 | ||
201 | if((float)ChargeEoI>Threshold) { | |
202 | NHit[i]+=1; | |
073bb639 | 203 | if(rawStream->GetBBFlag(i,iClock)) BBFlag[i]+=1; |
204 | if(rawStream->GetBGFlag(i,iClock)) BGFlag[i]+=1; | |
0ec8015e | 205 | } |
206 | } | |
207 | } | |
208 | delete rawStream; | |
209 | rawStream = 0x0; | |
210 | delete rawReader; | |
211 | rawReader = 0x0; | |
212 | } // end of switch on event type | |
213 | ||
214 | nevents_total++; | |
215 | /* free resources */ | |
216 | free(event); | |
217 | ||
218 | } // end of loop over events | |
219 | } // end of loop over data files | |
220 | ||
221 | //________________________________________________________________________ | |
222 | // Computes mean values, dumps them into the output text file | |
223 | Float_t fBB, fBG; | |
224 | for(Int_t i=0; i<64; i++) { | |
225 | if(NHit[i] > 0) { | |
226 | fBB = (float)BBFlag[i]/(float)NHit[i]; | |
227 | fBG = (float)BGFlag[i]/(float)NHit[i]; | |
228 | }else{ | |
229 | fBB = fBG = 0.; | |
230 | } | |
231 | fprintf(fp," %d %d %f %f\n",iteration,i,fBB,fBG); | |
232 | fprintf(flog," it %d ch %d BB %f BG %f BBFlag %d BGFlag %d Hit %d\n",iteration,i,fBB,fBG,BBFlag[i],BGFlag[i],NHit[i]); | |
233 | } | |
234 | ||
235 | //________________________________________________________________________ | |
236 | ||
237 | /* write report */ | |
238 | fprintf(flog,"Run #%s, received %d physics events out of %d\n",getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total); | |
239 | printf("Run #%s, received %d physics events out of %d\n",getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total); | |
240 | ||
241 | /* close result and log files */ | |
242 | fclose(fp); | |
243 | fclose(flog); | |
244 | ||
245 | /* report progress */ | |
246 | daqDA_progressReport(90); | |
247 | ||
248 | /* export result file to FES */ | |
249 | status=daqDA_FES_storeFile("./V0_ChannelDelayTuning.dat","V00da_ChannelDelayTuning"); | |
250 | if (status) { | |
251 | printf("Failed to export file : %d\n",status); | |
252 | return -1; } | |
253 | ||
254 | /* store result file on Online DB */ | |
255 | // status=daqDA_DB_storeFile("./V0_ChannelDelayTuning.dat","V00da_ChannelDelayTuning"); | |
256 | // if (status) { | |
257 | // printf("Failed to store file into Online DB: %d\n",status); | |
258 | // return -1; } | |
259 | ||
260 | /* report progress */ | |
261 | daqDA_progressReport(100); | |
262 | ||
263 | return status; | |
264 | } |