Mapping file put in the CDB
[u/mrichter/AliRoot.git] / PMD / PMDPEDESTALda.cxx
1 /*
2 PMD DA for online calibration
3
4 contact: basanta@phy.iitb.ac.in
5 Link:
6 Reference Run:/afs/cern.ch/user/b/bnandi/public/pedestaldata/run37820.raw
7 Run Type: PEDESTAL
8 DA Type: LDC
9 Number of events needed: 1000
10 Input Files: 
11 Output Files: pmd_ped.root, to be exported to the DAQ FXS, pedestal230*.ped
12 Trigger types used: PHYSICS_EVENT
13
14 */
15 extern "C" {
16 #include <daqDA.h>
17 }
18
19 #include "event.h"
20 #include "monitor.h"
21 //#include "daqDA.h"
22
23 #include <Riostream.h>
24 #include <stdio.h>
25 #include <stdlib.h>
26
27 //AliRoot
28 #include "AliRawReaderDate.h"
29 #include "AliPMDCalibPedestal.h"
30
31 //ROOT
32 #include "TFile.h"
33 #include "TH1F.h"
34 #include "TBenchmark.h"
35 #include "TTree.h"
36 #include "TROOT.h"
37 #include "TPluginManager.h"
38
39 /* Main routine
40       Arguments: 
41       1- monitoring data source
42 */
43 int main(int argc, char **argv) {
44   
45     /* magic line from Rene */
46     gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
47                                           "*",
48                                           "TStreamerInfo",
49                                           "RIO",
50                                           "TStreamerInfo()");
51
52     
53     AliPMDCalibPedestal calibped;
54
55     TTree *ped = NULL;
56
57     //    TH1F::AddDirectory(0);
58   
59       
60     // decoding the events
61   
62     int status;
63     
64     if (argc!=2) {
65         printf("Wrong number of arguments\n");
66         return -1;
67     }
68
69     
70     /* define data source : this is argument 1 */  
71     status=monitorSetDataSource( argv[1] );
72     if (status!=0) {
73         printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
74         return -1;
75     }
76     
77     /* declare monitoring program */
78     status=monitorDeclareMp( __FILE__ );
79     if (status!=0) {
80         printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
81         return -1;
82     }
83     
84     /* define wait event timeout - 1s max */
85     monitorSetNowait();
86     monitorSetNoWaitNetworkTimeout(1000);
87     
88     /* log start of process */
89     printf("PMD PED DA - started generating mean and rms of each channel\n");  
90     
91     /* init some counters */
92
93
94     struct eventHeaderStruct *event;
95     eventTypeType eventT = 0;
96     Int_t iev=0;
97     
98     /* main loop (infinite) */
99     for(;;) {
100         
101         /* check shutdown condition */
102         if (daqDA_checkShutdown()) {break;}
103         
104         /* get next event (blocking call until timeout) */
105         status=monitorGetEventDynamic((void **)&event);
106         if (status==MON_ERR_EOF) {
107             printf ("End of File detected\n");
108             break; /* end of monitoring file has been reached */
109         }
110         
111         if (status!=0) {
112             printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
113             break;
114         }
115         
116         /* retry if got no event */
117         if (event==NULL) {
118             continue;
119         }
120         
121         iev++; 
122         
123         /* use event - here, just write event id to result file */
124         eventT=event->eventType;
125         switch (event->eventType){
126             
127             /* START OF RUN */
128             case START_OF_RUN:
129                 break;
130                 /* END START OF RUN */
131                 
132                 /* END OF RUN */
133             case END_OF_RUN:
134                 break;
135                 
136             case PHYSICS_EVENT:
137
138                 AliRawReader *rawReader = new AliRawReaderDate((void*)event);
139                 TObjArray *pmdddlcont = new TObjArray();
140                 calibped.ProcessEvent(rawReader,pmdddlcont);
141                 
142                 delete pmdddlcont;
143                 pmdddlcont = 0;
144                 delete rawReader;
145                 rawReader = 0x0;
146         }
147         
148         /* free resources */
149         free(event);
150         
151         /* exit when last event received, no need to wait for TERM signal */
152
153         // either END_OF_RUN or 1000 events
154
155         if (iev == 1000)
156           {
157             printf("-- 1000 pedestal events processed : terminating --\n");
158             printf("-- eventT is set to END_OF_RUN --\n");
159             eventT = END_OF_RUN;
160             break;
161           }
162
163     }
164
165     printf(" Total number of events processed = %i \n",iev);
166
167     ped = new TTree("ped","PMD Pedestal tree");
168     
169     if (eventT==END_OF_RUN) {
170       printf("EOR event detected\n");
171       calibped.Analyse(ped);
172     }
173     
174     TFile * pedRun = new TFile ("PMD_PED.root","RECREATE"); 
175     ped->Write();
176     pedRun->Close();
177
178     delete ped;
179     ped = 0;
180
181 /* store the pedestal file in database */
182
183     status = daqDA_DB_storeFile("PMD_PED.root","PMD_PED.root");
184
185     if (!status)
186       {
187         printf("--- PEDESTAL root FILE IS IN THE DAQ DATABASE\n");
188       }
189     else
190       {
191         printf("--- Storing the root file into the database failed\n");
192       }
193
194 /* store the pedestal file on FES */
195  
196     status = daqDA_FES_storeFile("PMD_PED.root","PMD_PED.root");
197     if (status) {
198         status = -2;
199     }
200
201
202
203
204     return status;
205 }