adding some cuts for systematic study
[u/mrichter/AliRoot.git] / T0 / T0Cosmicda.cxx
1 /*
2 T0 DA for online calibration
3  
4 Contact: Michal.Oledzki@cern.ch
5 Link: http://users.jyu.fi/~mioledzk/
6 Run Type: PHYSICS
7 DA Type: MON
8 Number of events needed: 400000 
9 Input Files: inCosmicLaser.dat, external parameters
10 Output Files: daCosmicLaser.root, to be exported to the DAQ FXS
11 Trigger types used: PHYSICS_EVENT
12
13 */
14
15 #define FILE_OUT "daCosmic.root"
16 //
17 //author Michal Oledzki
18 //DA for old cosmic runs
19 //
20 //
21 #define FILE_IN "inCosmic.dat"
22 #include <daqDA.h>
23 #include <event.h>
24 #include <monitor.h>
25  
26 #include <Riostream.h>
27 #include <stdio.h>
28 #include <stdlib.h>
29
30 //AliRoot
31 #include <AliRawReaderDate.h>
32 #include <AliRawReader.h>
33 #include <AliT0RawReader.h>
34
35 //ROOT
36 #include "TROOT.h"
37 #include "TPluginManager.h"
38 #include "TFile.h"
39 #include "TKey.h"
40 #include "TH2S.h"
41 #include "TObject.h"
42 #include "TBenchmark.h"
43 #include "TRandom.h"
44 #include "TMath.h"
45 #include "TCanvas.h"
46 #include "TString.h"
47 #include "TH1.h"
48 #include "TF1.h"
49 #include "TSpectrum.h"
50 #include "TVirtualFitter.h"
51 //#include "TProfile.h"
52 int kCcbx;
53 float kCclx,kCcmx;
54 /* Main routine
55       Arguments: 
56       1- monitoring data source
57 */
58 int main(int argc, char **argv) {
59   int status;
60   
61   /* magic line */
62   gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
63                                         "*",
64                                         "TStreamerInfo",
65                                         "RIO",
66                                         "TStreamerInfo()"); 
67   
68   if(daqDA_DB_getFile(FILE_IN, FILE_IN)){
69      printf("Couldn't get input file >>inCosmic.dat<< from DAQ_DB !!!\n");
70      return -1;
71   }
72   
73  
74   FILE *inp;
75   char c;
76   inp = fopen(FILE_IN, "r");  
77   if(!inp){
78         printf("Input file >>inCosmic.dat<< not found !!!\n");
79         return -1;
80   }
81
82   while((c=getc(inp))!=EOF) {
83     switch(c) {
84       case 'a': {fscanf(inp, "%d", &kCcbx ); break;} //N of X bins hCFD1_CFD
85       case 'b': {fscanf(inp, "%f", &kCclx ); break;} //Low x hCFD1_CFD
86       case 'c': {fscanf(inp, "%f", &kCcmx ); break;} //High x hCFD1_CFD
87     }
88   }
89   fclose(inp);
90
91   if (argc!=2) {
92     printf("Wrong number of arguments\n");
93     return -1;
94   }
95
96
97   /* define data source : this is argument 1 */  
98   status=monitorSetDataSource( argv[1] );
99   if (status!=0) {
100     printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
101     return -1;
102   }
103
104
105   /* declare monitoring program */
106   status=monitorDeclareMp( __FILE__ );
107   if (status!=0) {
108     printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
109     return -1;
110   }
111
112
113   /* define wait event timeout - 1s max */
114   monitorSetNowait();
115   monitorSetNoWaitNetworkTimeout(1000);
116   
117
118   /* log start of process */
119   printf("T0 monitoring program started\n");  
120
121   // Allocation of histograms - start
122   TH1F *hCFD1minCFD[24]; 
123
124    for(Int_t ic=0; ic<24; ic++) {
125       hCFD1minCFD[ic] = new TH1F(Form("CFD1-CFD%d",ic+1),"CFD-CFD",kCcbx,kCclx,kCcmx);
126     }
127
128   // Allocation of histograms - end
129
130   Int_t iev=0;
131   /* main loop (infinite) */
132   for(;;) {
133     struct eventHeaderStruct *event;
134     eventTypeType eventT;
135   
136     /* check shutdown condition */
137     if (daqDA_checkShutdown()) {break;}
138     
139     /* get next event (blocking call until timeout) */
140     status=monitorGetEventDynamic((void **)&event);
141     if (status==(int)MON_ERR_EOF) {
142       printf ("End of File detected\n");
143       break; /* end of monitoring file has been reached */
144     }
145     
146     if (status!=0) {
147       printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
148       break;
149     }
150
151     /* retry if got no event */
152     if (event==NULL) {
153       continue;
154     }
155  
156 //    iev++; 
157
158     /* use event - here, just write event id to result file */
159     eventT=event->eventType;
160    
161     switch (event->eventType){
162
163       case START_OF_RUN:
164         break;
165
166       case END_OF_RUN:
167         break;
168
169       case CALIBRATION_EVENT:
170 //      case PHYSICS_EVENT:
171       iev++; 
172
173       if(iev==1){
174            printf("First event - %i\n",iev);
175       }
176
177       // Initalize raw-data reading and decoding
178       AliRawReader *reader = new AliRawReaderDate((void*)event);
179           
180       // Enable the following two lines in case of real-data
181           reader->RequireHeader(kTRUE);
182       AliT0RawReader *start = new AliT0RawReader(reader, kTRUE);
183
184       // Read raw data
185       Int_t allData[105][5];
186       for(Int_t i0=0;i0<105;i0++)
187         for(Int_t j0=0;j0<5;j0++)
188                 allData[i0][j0] = 0;
189     
190        if(start->Next()){
191        for (Int_t i=0; i<105; i++) {
192         for(Int_t iHit=0;iHit<5;iHit++){
193           allData[i][iHit]= start->GetData(i,iHit);
194         }
195        }
196       }
197         else
198         printf("No T0 data found!!!\n");
199
200       // Fill the histograms
201         
202       for (Int_t ik = 0; ik<24; ik++)
203          for (Int_t iHt=0; iHt<5; iHt++){
204                 if(allData[ik+1][iHt]!=0 ){
205                   if(ik<12){
206                          hCFD1minCFD[ik]->Fill(allData[ik+1][iHt]-allData[1][iHt]);
207                   }
208                   if(ik>11){
209                          hCFD1minCFD[ik]->Fill(allData[ik+45][iHt]-allData[57][iHt]);
210                   }
211                 }
212         }
213         
214      delete start;
215         start = 0x0;
216      reader->Reset();
217       // End of fill histograms
218
219     }
220
221     /* free resources */
222     free(event);
223     
224     /* exit when last event received, no need to wait for TERM signal */
225     if (eventT==END_OF_RUN) {
226       printf("EOR event detected\n");
227       printf("Number of events processed - %i\n ",iev);
228       break;
229     }
230   }
231   printf("After loop, before writing histos\n");
232   // write a file with the histograms
233   TFile *hist = new TFile(FILE_OUT,"RECREATE");
234
235   for(Int_t j=0;j<24;j++){
236      hCFD1minCFD[j]->Write();
237     }
238   hist->Close();
239   delete hist;
240
241   status=0;
242
243   /* export file to FXS */
244   if (daqDA_FES_storeFile(FILE_OUT, "COSMIC")) {
245     status=-2;
246   }
247
248   return status;
249 }
250