For Pythia with tune don't switch off MI in ConfigHeavyFlavor
[u/mrichter/AliRoot.git] / ITS / ITSSPDVertexDiamondda.cxx
1 /*
2 Contact: cvetan.cheshkov@cern.ch
3 Link: 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
4 Reference Run: 58338
5 Run Type: PHYSICS
6 DA Type: MON
7 Number of events needed: 100
8 Input Files: GRP/Geometry/Data , ITS/Align/Data , spd_noisy_ocdb , spd_dead_ocdb , TRIGGER/SPD/PITConditions (all the files are taken from SPD daqDetDB)
9 Output Files: SPDVertexDiamondDA.root
10 Trigger types used: PHYSICS, SPD-F0 
11 */
12
13 #define OUTPUT_FILE "SPDVertexDiamondDA.root"
14 #define N_EVENTS_AUTOSAVE 50
15
16 extern "C" {
17 #include "daqDA.h"
18 }
19
20 #include "event.h"
21 #include "monitor.h"
22
23 #ifdef ALI_AMORE
24 #include <AmoreDA.h>
25 //int amore::da::Updated(char const*) {}
26 #endif
27
28 #include <TPluginManager.h>
29 #include <TROOT.h>
30 #include <TH1.h>
31 #include <TH2.h>
32 #include <TSystem.h>
33
34 #include "AliLog.h"
35 #include "AliRawReaderDate.h"
36 #include "AliCDBManager.h"
37 #include "AliITSMeanVertexer.h"
38
39 int main(int argc, char **argv) {
40
41   gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
42                                         "*",
43                                         "TStreamerInfo",
44                                         "RIO",
45                                         "TStreamerInfo()"); 
46
47   int status;
48   if (argc<2) {
49     printf("Wrong number of arguments\n");
50     return -1;
51   }
52
53   /* define data source : this is argument 1 */  
54   status=monitorSetDataSource( argv[1] );
55   if (status!=0) {
56     printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
57     return -1;
58   }
59
60   /* declare monitoring program */
61   status=monitorDeclareMp( __FILE__ );
62   if (status!=0) {
63     printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
64     return -1;
65   }
66
67   /* define wait event timeout - 1s max */
68   monitorSetNowait();
69   monitorSetNoWaitNetworkTimeout(1000);
70   
71   /* log start of process */
72   printf("Vertex-Diamond SPD DA started\n");  
73
74   /* init some counters */
75   int nevents_with_vertex = 0;
76   int nevents_physics=0;
77   int nevents_total=0;
78
79   struct eventHeaderStruct *event;
80   eventTypeType eventT;
81
82   // Get run number
83   if (getenv("DATE_RUN_NUMBER")==0) {
84     printf("DATE_RUN_NUMBER not properly set.\n");
85     return -1;
86   }
87   int runNr = atoi(getenv("DATE_RUN_NUMBER"));
88
89   // Get the necessary OCDB files from the DAQ detector DB
90   if (gSystem->AccessPathName("localOCDB/GRP/Geometry/Data",kFileExists)) {
91     if (gSystem->mkdir("localOCDB/GRP/Geometry/Data",kTRUE) != 0) {
92       printf("Failed to create directory: localOCDB/GRP/Geometry/Data");
93       return -1;
94     }
95   }
96   status = daqDA_DB_getFile("GRP/Geometry/Data","localOCDB/GRP/Geometry/Data/Run0_999999999_v0_s0.root");
97   if (status) {
98     printf("Failed to get geometry file (GRP/Geometry/Data) from DAQdetDB, status=%d\n", status);
99     return -1;
100   }
101
102   if (gSystem->AccessPathName("localOCDB/ITS/Align/Data",kFileExists)) {
103     if (gSystem->mkdir("localOCDB/ITS/Align/Data",kTRUE) != 0) {
104       printf("Failed to create directory: localOCDB/ITS/Align/Data");
105       return -1;
106     }
107   }
108   status = daqDA_DB_getFile("ITS/Align/Data","localOCDB/ITS/Align/Data/Run0_999999999_v0_s0.root");
109   if (status) {
110     printf("Failed to get its-alignment file (ITS/Align/Data) from DAQdetDB, status=%d\n", status);
111     return -1;
112   }
113
114   if (gSystem->AccessPathName("localOCDB/ITS/Calib/SPDNoisy",kFileExists)) {
115     if (gSystem->mkdir("localOCDB/ITS/Calib/SPDNoisy",kTRUE) != 0) {
116       printf("Failed to create directory: localOCDB/ITS/Calib/SPDNoisy");
117       return -1;
118     }
119   }
120   status = daqDA_DB_getFile("spd_noisy_ocdb","localOCDB/ITS/Calib/SPDNoisy/Run0_999999999_v0_s0.root");
121   if (status) {
122     printf("Failed to get spd file (spd_noisy_ocdb) from DAQdetDB, status=%d\n", status);
123     return -1;
124   }
125
126   if (gSystem->AccessPathName("localOCDB/ITS/Calib/SPDDead",kFileExists)) {
127     if (gSystem->mkdir("localOCDB/ITS/Calib/SPDDead",kTRUE) != 0) {
128       printf("Failed to create directory: localOCDB/ITS/Calib/SPDDead");
129       return -1;
130     }
131   }
132   status = daqDA_DB_getFile("spd_dead_ocdb","localOCDB/ITS/Calib/SPDDead/Run0_999999999_v0_s0.root");
133   if (status) {
134     printf("Failed to get spd file (spd_dead_ocdb) from DAQdetDB, status=%d\n", status);
135     return -1;
136   }
137
138   if (gSystem->AccessPathName("localOCDB/TRIGGER/SPD/PITConditions",kFileExists)) {
139     if (gSystem->mkdir("localOCDB/TRIGGER/SPD/PITConditions",kTRUE) != 0) {
140       printf("Failed to create directory: localOCDB/TRIGGER/SPD/PITConditions");
141       return -1;
142     }
143   }
144   status = daqDA_DB_getFile("TRIGGER/SPD/PITConditions","localOCDB/TRIGGER/SPD/PITConditions/Run0_999999999_v0_s0.root");
145   if (status) {
146     printf("Failed to get spd trigger file (TRIGGER/SPD/PITConditions) from DAQdetDB, status=%d\n", status);
147     return -1;
148   }
149
150   // Global initializations
151   AliLog::SetGlobalLogLevel(AliLog::kError);
152   AliCDBManager *man = AliCDBManager::Instance();
153   man->SetDefaultStorage("local://localOCDB");
154   man->SetRun(runNr);
155
156   // Init mean vertexer
157   AliITSMeanVertexer *mv = new AliITSMeanVertexer();
158   if (!mv->Init()) {
159     printf("Initialization of mean vertexer object failed ! Check the log for details");
160     return -1;
161   }
162
163   // Initialization of AMORE sender
164 #ifdef ALI_AMORE
165   amore::da::AmoreDA vtxAmore(amore::da::AmoreDA::kSender);
166 #endif
167   /* main loop (infinite) */
168   for(;;) {
169     
170     /* check shutdown condition */
171     if (daqDA_checkShutdown()) {break;}
172     
173     /* get next event (blocking call until timeout) */
174     status=monitorGetEventDynamic((void **)&event);
175     if (status==MON_ERR_EOF) {
176       printf ("End of File detected\n");
177       break; /* end of monitoring file has been reached */
178     }
179     
180     if (status!=0) {
181       printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
182       break;
183     }
184     
185     /* retry if got no event */
186     if (event==NULL) {
187       continue;
188     }
189
190     nevents_total++;
191     eventT=event->eventType;
192     switch (event->eventType){
193       
194       /* START OF RUN */
195     case START_OF_RUN:
196       break;
197       /* END START OF RUN */
198       
199     /* END OF RUN */
200     case END_OF_RUN:
201       break;
202       
203     case PHYSICS_EVENT:
204       nevents_physics++;
205       AliRawReader *rawReader = new AliRawReaderDate((void*)event);
206
207       // Run mean-vertexer reco
208       if (mv->Reconstruct(rawReader)) nevents_with_vertex++;
209
210       // Auto save
211       if ((nevents_physics%N_EVENTS_AUTOSAVE) == 0) {
212         mv->WriteVertices(OUTPUT_FILE);
213
214 #ifdef ALI_AMORE
215       // send the histos to AMORE pool
216         printf("AMORE send status: %d\n",vtxAmore.Send(mv->GetVertexXY()->GetName(),mv->GetVertexXY()));
217         printf("AMORE send status: %d\n",vtxAmore.Send(mv->GetVertexZ()->GetName(),mv->GetVertexZ()));
218 #endif
219       }
220
221       delete rawReader;
222     }
223     
224     /* free resources */
225     free(event);
226     
227     /* exit when last event received, no need to wait for TERM signal */
228     if (eventT==END_OF_RUN) {
229       printf("EOR event detected\n");
230       break;
231     }
232   }
233
234   mv->WriteVertices(OUTPUT_FILE);
235
236 #ifdef ALI_AMORE
237   // send the histos to AMORE pool
238   printf("AMORE send status: %d\n",vtxAmore.Send(mv->GetVertexXY()->GetName(),mv->GetVertexXY()));
239   printf("AMORE send status: %d\n",vtxAmore.Send(mv->GetVertexZ()->GetName(),mv->GetVertexZ()));
240 #endif
241
242   delete mv;
243
244   /* write report */
245   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);
246
247   status=0;
248
249   /* export file to FXS */
250   if (daqDA_FES_storeFile(OUTPUT_FILE, "VertexDiamond")) {
251     status=-2;
252   }
253   
254   return status;
255 }