TOF + macros to each detector folder
[u/mrichter/AliRoot.git] / ITS / ITSSPDVertexDiamondda.cxx
CommitLineData
eb35e591 1/*
2e0441d1 2Contact: cvetan.cheshkov@cern.ch
3Link: http://alisoft.cern.ch/viewvc/trunk/ITS/ITSSPDVertexDiamondda.cxx?root=AliRoot&view=log , /afs/cern.ch/user/c/cheshkov/public/ITS/VD_da_test.date , /afs/cern.ch/user/c/cheshkov/public/08000058338016.30.root.date.gz
4Reference Run: 58338
5Run Type: PHYSICS
6DA Type: MON
7Number of events needed: 100
8Input Files: GRP/Geometry/Data , ITS/Align/Data , spd_noisy_ocdb , spd_dead_ocdb , TRIGGER/SPD/PITConditions (all the files are taken from SPD daqDetDB)
9Output Files: SPDVertexDiamondDA.root
10Trigger types used: PHYSICS, SPD-F0
11 */
eb35e591 12
97a75e92 13#define OUTPUT_FILE "SPDVertexDiamondDA.root"
eb35e591 14
15extern "C" {
16#include "daqDA.h"
17}
18
19#include "event.h"
20#include "monitor.h"
21
2e0441d1 22#define ALI_AMORE
23
55c5e86d 24#ifdef ALI_AMORE
25#include <AmoreDA.h>
308c2f7c 26//int amore::da::Updated(char const*) {}
55c5e86d 27#endif
28
eb35e591 29#include <TPluginManager.h>
30#include <TROOT.h>
cf0b1cea 31#include <TH1.h>
32#include <TH2.h>
694ea942 33#include <TPaveText.h>
6bf5b296 34#include <TSystem.h>
ee0a8262 35#include <TGeoGlobalMagField.h>
13cf0903 36#include <TTimeStamp.h>
eb35e591 37
e1f6be44 38#include "AliLog.h"
ee0a8262 39#include "AliMagF.h"
eb35e591 40#include "AliRawReaderDate.h"
eb35e591 41#include "AliCDBManager.h"
cf0b1cea 42#include "AliITSMeanVertexer.h"
cc7ca4a6 43#include <cstdlib>
eb35e591 44
45int main(int argc, char **argv) {
46
2e0441d1 47 gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
48 "*",
49 "TStreamerInfo",
50 "RIO",
51 "TStreamerInfo()");
52
53 int status;
54 if (argc<2) {
55 printf("Wrong number of arguments\n");
56 return -1;
57 }
58 // get mean vertex quality cuts
59 status = daqDA_DB_getFile("ITSSPD_VertexQualityTuning_DA.config","./ITSSPD_VertexQualityTuning_DA.config");
60 if (status) {
61 printf("Failed to get config file (ITSSPD_VertexQualityTuning_DA.config) from DAQ DB, status=%d\n", status);
62 return -1;
63 }
64 /* open the config file and retrieve running parameters */
65 Float_t errX = -1;
66 Float_t r = -1;
67 UInt_t minClInner = 999 ;
68 UInt_t maxClInner = 999 ;
69 Int_t nEvFirstLoop = 0;
70 Int_t nEvAUTOSAVE = 0;
d3d3859c 71 Int_t zFiducialRegion=0;
13cf0903 72 Int_t timeWindowExport=0;
73 Int_t timeErrWindowExport=0;
74
d3d3859c 75
2e0441d1 76 char name[10][10];
77
78 FILE *fpConfig = fopen("ITSSPD_VertexQualityTuning_DA.config","r");
13cf0903 79 fscanf(fpConfig,"%s %f\n %s %f\n %s %d\n %s %d \n %s %d \n %s %d \n %s %d \n %s %d \n %s %d",&name[0], &errX, &name[1], &r, &name[2], &minClInner,&name[3], &maxClInner, &name[4],&nEvFirstLoop,&name[5],&nEvAUTOSAVE,&name[6],&zFiducialRegion, &name[7], &timeWindowExport, &name[8], &timeErrWindowExport);
80
2e0441d1 81 fclose(fpConfig);
82
13cf0903 83 printf("\n\n Mean Vertex quality cuts : \n- errX = %f\n- r = %f\n- minSPD0 = %d maxSPD0 = %d\n- nEventsFirstLoop = %d nEventsAUTOSAVE = %d \n- zFiducialRegion = %d\n- timeWindowExport = %d \n- timeErrWindowExport = %d \n\n\n",errX,r,minClInner,maxClInner,nEvFirstLoop,nEvAUTOSAVE,zFiducialRegion, timeWindowExport, timeErrWindowExport);
2e0441d1 84
85 /* define data source : this is argument 1 */
86 status=monitorSetDataSource( argv[1] );
87 if (status!=0) {
88 printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
89 return -1;
90 }
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 /* define wait event timeout - 1s max */
100 monitorSetNowait();
101 monitorSetNoWaitNetworkTimeout(1000);
102
103 /* log start of process */
104 printf("Vertex-Diamond SPD DA started\n");
105
106 /* init some counters */
107 int nevents_with_vertex = 0;
108 int nevents_physics=0;
109 int nevents_total=0;
110
111 struct eventHeaderStruct *event;
112 eventTypeType eventT;
113
114 // Get run number
115 if (getenv("DATE_RUN_NUMBER")==0) {
116 printf("DATE_RUN_NUMBER not properly set.\n");
117 return -1;
118 }
119 int runNr = atoi(getenv("DATE_RUN_NUMBER"));
120
121 // Get the necessary OCDB files from the DAQ detector DB
122 if (gSystem->AccessPathName("localOCDB/GRP/Geometry/Data",kFileExists)) {
123 if (gSystem->mkdir("localOCDB/GRP/Geometry/Data",kTRUE) != 0) {
124 printf("Failed to create directory: localOCDB/GRP/Geometry/Data");
125 return -1;
694ea942 126 }
2e0441d1 127 }
128 status = daqDA_DB_getFile("GRP/Geometry/Data","localOCDB/GRP/Geometry/Data/Run0_999999999_v0_s0.root");
129 if (status) {
130 printf("Failed to get geometry file (GRP/Geometry/Data) from DAQdetDB, status=%d\n", status);
131 return -1;
132 }
133
134 if (gSystem->AccessPathName("localOCDB/ITS/Align/Data",kFileExists)) {
135 if (gSystem->mkdir("localOCDB/ITS/Align/Data",kTRUE) != 0) {
136 printf("Failed to create directory: localOCDB/ITS/Align/Data");
137 return -1;
eb35e591 138 }
2e0441d1 139 }
140 status = daqDA_DB_getFile("ITS/Align/Data","localOCDB/ITS/Align/Data/Run0_999999999_v0_s0.root");
141 if (status) {
142 printf("Failed to get its-alignment file (ITS/Align/Data) from DAQdetDB, status=%d\n", status);
143 return -1;
144 }
145
146 if (gSystem->AccessPathName("localOCDB/ITS/Calib/SPDNoisy",kFileExists)) {
147 if (gSystem->mkdir("localOCDB/ITS/Calib/SPDNoisy",kTRUE) != 0) {
148 printf("Failed to create directory: localOCDB/ITS/Calib/SPDNoisy");
149 return -1;
6bf5b296 150 }
2e0441d1 151 }
152 status = daqDA_DB_getFile("spd_noisy_ocdb","localOCDB/ITS/Calib/SPDNoisy/Run0_999999999_v0_s0.root");
153 if (status) {
154 printf("Failed to get spd file (spd_noisy_ocdb) from DAQdetDB, status=%d\n", status);
155 return -1;
156 }
157
158 if (gSystem->AccessPathName("localOCDB/ITS/Calib/SPDDead",kFileExists)) {
159 if (gSystem->mkdir("localOCDB/ITS/Calib/SPDDead",kTRUE) != 0) {
160 printf("Failed to create directory: localOCDB/ITS/Calib/SPDDead");
161 return -1;
6bf5b296 162 }
2e0441d1 163 }
164 status = daqDA_DB_getFile("spd_dead_ocdb","localOCDB/ITS/Calib/SPDDead/Run0_999999999_v0_s0.root");
165 if (status) {
166 printf("Failed to get spd file (spd_dead_ocdb) from DAQdetDB, status=%d\n", status);
167 return -1;
168 }
169
170 if (gSystem->AccessPathName("localOCDB/ITS/Calib/SPDSparseDead",kFileExists)) {
171 if (gSystem->mkdir("localOCDB/ITS/Calib/SPDSparseDead",kTRUE) != 0) {
172 printf("Failed to create directory: localOCDB/ITS/Calib/SPDSparseDead");
173 return -1;
6bf5b296 174 }
2e0441d1 175 }
176
177 status = daqDA_DB_getFile("spd_sparsedead_ocdb","localOCDB/ITS/Calib/SPDSparseDead/Run0_999999999_v0_s0.root");
178 if (status) {
179 printf("Failed to get spd file (spd_sparsedead_ocdb) from DAQdetDB, status=%d\n", status);
180 return -1;
181 }
182
183 if (gSystem->AccessPathName("localOCDB/TRIGGER/SPD/PITConditions",kFileExists)) {
184 if (gSystem->mkdir("localOCDB/TRIGGER/SPD/PITConditions",kTRUE) != 0) {
185 printf("Failed to create directory: localOCDB/TRIGGER/SPD/PITConditions");
186 return -1;
6bf5b296 187 }
2e0441d1 188 }
189 status = daqDA_DB_getFile("TRIGGER/SPD/PITConditions","localOCDB/TRIGGER/SPD/PITConditions/Run0_999999999_v0_s0.root");
190 if (status) {
191 printf("Failed to get spd trigger file (TRIGGER/SPD/PITConditions) from DAQdetDB, status=%d\n", status);
192 return -1;
193 }
194
195 status = daqDA_DB_getFile("mfchebKGI_sym.root","localOCDB/mfchebKGI_sym.root");
196 if (status) {
197 printf("Failed to get spd file (mfchebKGI_sym.root) from DAQdetDB, status=%d\n", status);
198 return -1;
199 }
200
201 // Global initializations
202
203 // The B filed is required in AliITSClusterFinderV2SPD
204 // for the Lorentz angle correction. B set to 0.
205 TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform,AliMagF::kBeamTypeAA,-1,2,15,"localOCDB/mfchebKGI_sym.root"));
206 AliLog::SetGlobalLogLevel(AliLog::kError);
207 AliCDBManager *man = AliCDBManager::Instance();
208 man->SetDefaultStorage("local://localOCDB");
209 man->SetRun(runNr);
210
211 // Init mean vertexer
212 AliITSMeanVertexer *mv = new AliITSMeanVertexer();
4f40f207 213
214 // Optionally change the Z fiducial region here (default +/- 40 cm
d3d3859c 215 mv->SetZFiducialRegion(zFiducialRegion);
4f40f207 216
2e0441d1 217 if (!mv->Init()) {
218 printf("Initialization of mean vertexer object failed ! Check the log for details");
219 return -1;
220 }
221
222 mv->SetCutOnErrX(errX);
223 mv->SetCutOnR(r);
224 mv->SetCutOnCls(minClInner,maxClInner);
225
e0672fdb 226 gSystem->Exec(Form("rm %s",OUTPUT_FILE));
227
13cf0903 228 TTimeStamp *timeStamp;
229 timeStamp->Set();
230 Int_t t1 = timeStamp->GetSec();
231 Int_t t2=0;
232
233// Initialization of AMORE sender
2e0441d1 234#ifdef ALI_AMORE
235 amore::da::AmoreDA vtxAmore(amore::da::AmoreDA::kSender);
236#endif
a33265ea 237 Int_t amore_status = 0;
2e0441d1 238 /* main loop (infinite) */
239 for(;;) {
6bf5b296 240
2e0441d1 241 /* check shutdown condition */
242 if (daqDA_checkShutdown()) {break;}
69ebd58d 243
2e0441d1 244 /* get next event (blocking call until timeout) */
245 status=monitorGetEventDynamic((void **)&event);
246 if (status==MON_ERR_EOF) {
247 printf ("End of File detected\n");
248 break; /* end of monitoring file has been reached */
69ebd58d 249 }
250
2e0441d1 251 if (status!=0) {
252 printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
253 break;
2a16d0d6 254 }
bdffed96 255
2e0441d1 256 /* retry if got no event */
257 if (event==NULL) {
258 continue;
37add1d1 259 }
eb35e591 260
2e0441d1 261 nevents_total++;
262 eventT=event->eventType;
263 switch (event->eventType){
264
265 /* START OF RUN */
266 case START_OF_RUN:
267 break;
268 /* END START OF RUN */
269
270 /* END OF RUN */
271 case END_OF_RUN:
272 break;
273
274 case PHYSICS_EVENT:
275 nevents_physics++;
276 AliRawReader *rawReader = new AliRawReaderDate((void*)event);
277
278 // Run mean-vertexer reco
279 if (mv->Reconstruct(rawReader)) nevents_with_vertex++;
280 if(nevents_physics < nEvFirstLoop) continue;
281 // Auto save
282 if ((nevents_physics%nEvAUTOSAVE) == 0) {
283 TH2F *histo = ((TH2F*)mv->GetVertexXY());
284 histo->SetStats(kFALSE);
285 histo->SetTitle("");
286 histo->GetListOfFunctions()->SetOwner(kTRUE);
287 if(histo->GetListOfFunctions()->GetEntries()<1) histo->GetListOfFunctions()->Add(new TPaveText(-5,4.5,5.,6.2,"br"));
288 for(Int_t i=0; i<histo->GetListOfFunctions()->GetEntries(); i++){
289 TString funcName = histo->GetListOfFunctions()->At(i)->ClassName();
290 if(funcName.Contains("TPaveText")){
291 TPaveText *p = (TPaveText*)histo->GetListOfFunctions()->At(i);
292 p->Clear();
293 p->AddText(Form("%f events with vertex (%i out of %i processed events)",((Double_t)nevents_with_vertex)/((Double_t)nevents_physics),nevents_with_vertex,nevents_physics));
294 p->AddText(Form("%f events with good vertex (%i out of %i events with vertex)",histo->GetEntries()/((Double_t)nevents_with_vertex),(Int_t)histo->GetEntries(),nevents_with_vertex));
295 }
296 }
97a75e92 297 mv->WriteVertices(OUTPUT_FILE);
13cf0903 298
299 timeStamp->Set();
300 t2 = timeStamp->GetSec();
301
302
303
55c5e86d 304#ifdef ALI_AMORE
13cf0903 305
306 if (TMath::Abs((t2-t1)-timeWindowExport) < timeErrWindowExport){
307 t1=t2;
308 // send the histos to AMORE pool
309 amore_status=vtxAmore.Send(mv->GetVertexXY()->GetName(),mv->GetVertexXY());
310 if(amore_status) printf("AMORE XY send status: %d\n",amore_status);
311 TH1D *hProj = mv->GetVertexXY()->ProjectionX();
312 amore_status=vtxAmore.Send(Form("%s_ProjX",mv->GetVertexXY()->GetName()),hProj);
313 if(amore_status) printf("AMORE X send status: %d\n",amore_status);
314 if(hProj) delete hProj;
315 hProj = mv->GetVertexXY()->ProjectionY();
316 amore_status=vtxAmore.Send(Form("%s_ProjY",mv->GetVertexXY()->GetName()),hProj);
317 if(amore_status) printf("AMORE Y send status: %d\n",amore_status);
318 if(hProj) hProj->Delete();
319 amore_status=vtxAmore.Send(mv->GetVertexZ()->GetName(),mv->GetVertexZ());
320 if(amore_status) printf("AMORE Z send status: %d\n",amore_status);
321 }
322 if (TMath::Abs(t2-t1)> timeWindowExport) t1=t2;
323
55c5e86d 324#endif
eb35e591 325 }
326
2e0441d1 327 delete rawReader;
328 }
eb35e591 329
2e0441d1 330 /* free resources */
331 free(event);
6bf5b296 332
2e0441d1 333 /* exit when last event received, no need to wait for TERM signal */
334 if (eventT==END_OF_RUN) {
335 printf("EOR event detected\n");
336 break;
eb35e591 337 }
2e0441d1 338 }
339
340 TH2F *histo = ((TH2F*)mv->GetVertexXY());
341 histo->SetStats(kFALSE);
342 histo->SetTitle("");
343 histo->GetListOfFunctions()->SetOwner(kTRUE);
344 if(histo->GetListOfFunctions()->GetEntries()<1) histo->GetListOfFunctions()->Add(new TPaveText(-5,4.5,5.,6.2,"br"));
345 for(Int_t i=0; i<histo->GetListOfFunctions()->GetEntries(); i++){
346 TString funcName = histo->GetListOfFunctions()->At(i)->ClassName();
347 if(funcName.Contains("TPaveText")){
348 TPaveText *p = (TPaveText*)histo->GetListOfFunctions()->At(i);
349 p->Clear();
350 p->AddText(Form("%f events with vertex (%i out of %i processed events)",((Double_t)nevents_with_vertex)/((Double_t)nevents_physics),nevents_with_vertex,nevents_physics));
351 p->AddText(Form("%f events with good vertex (%i out of %i events with vertex)",histo->GetEntries()/((Double_t)nevents_with_vertex),(Int_t)histo->GetEntries(),nevents_with_vertex));
694ea942 352 }
2e0441d1 353 }
354 mv->WriteVertices(OUTPUT_FILE);
eb35e591 355
55c5e86d 356#ifdef ALI_AMORE
13cf0903 357 if (TMath::Abs((t2-t1)-timeWindowExport) < timeErrWindowExport){
358 t1=t2;
359 // send the histos to AMORE pool
360 amore_status=vtxAmore.Send(mv->GetVertexXY()->GetName(),mv->GetVertexXY());
361 if(amore_status) printf("AMORE XY send status: %d\n",amore_status);
362 TH1D *hProj = mv->GetVertexXY()->ProjectionX();
363 amore_status=vtxAmore.Send(Form("%s_ProjX",mv->GetVertexXY()->GetName()),hProj);
364 if(amore_status) printf("AMORE X send status: %d\n",amore_status);
365 if(hProj) delete hProj;
366 hProj = mv->GetVertexXY()->ProjectionY();
367 amore_status=vtxAmore.Send(Form("%s_ProjY",mv->GetVertexXY()->GetName()),hProj);
368 if(amore_status) printf("AMORE Y send status: %d\n",amore_status);
369 if(hProj) hProj->Delete();
370 amore_status=vtxAmore.Send(mv->GetVertexZ()->GetName(),mv->GetVertexZ());
371 if(amore_status) printf("AMORE Z send status: %d\n",amore_status);
372 }
373 if (TMath::Abs(t2-t1)> timeWindowExport) t1=t2;
55c5e86d 374
13cf0903 375#endif
376
2e0441d1 377 delete mv;
eb35e591 378
2e0441d1 379 /* write report */
380 printf("Run #%s, received %d events with vertex, out of %d physics and out of %d total events\n",getenv("DATE_RUN_NUMBER"),nevents_with_vertex,nevents_physics,nevents_total);
eb35e591 381
2e0441d1 382 status=0;
eb35e591 383
2e0441d1 384 /* export file to FXS */
385 if (daqDA_FES_storeFile(OUTPUT_FILE, "VertexDiamond")) {
386 status=-2;
387 }
388
389 return status;
eb35e591 390}