TRD/TRDVDRIFTda.cxx

   1 /*
   2
   3 Contact: r.bailhache@gsi.de
   4 Link:
   5 Reference run: 25909
   6 Run Type: PHYSICS STANDALONE
   7 DA Type: MON
   8 Number of events needed: many
   9 Input Files:  TRD raw files
  10 Output Files: trdCalibrationv.root
  11 Trigger types used: PHYSICS_EVENT
  12
  13 */
  14
  15 #define RESULT_FILE  "trdVdrift.root"
  16 #define FILE_ID "VDRIFT"
  17
  18
  19 extern "C" {
  20 #include <daqDA.h>
  21 }
  22
  23 #include "event.h"
  24 #include "monitor.h"
  25 #include <stdio.h>
  26 #include <stdlib.h>
  27
  28
  29 //
  30 //Root includes
  31 //
  32 #include <TObjArray.h>
  33 #include <TString.h>
  34 #include <TVectorF.h>
  35 #include <TDirectory.h>
  36 #include <TSystem.h>
  37 #include <TFile.h>
  38 #include "TROOT.h"
  39 #include "TPluginManager.h"
  40 //
  41 //AliRoot includes
  42 //
  43 #include "AliRawReader.h"
  44 #include "AliRawReaderDate.h"
  45 //#include "AliTRDrawFastStream.h"
  46 //#include "AliTRDrawStreamBase.h"
  47 #include "AliLog.h"
  48
  49 //
  50 //AMORE
  51 //
  52 #include <AmoreDA.h>
  53
  54 //
  55 // TRD calibration algorithm includes
  56 //
  57 #include "AliTRDCalibraFillHisto.h"
  58
  59 //functions, implementation below
  60 //void SendToAmoreDB(TObject *o, unsigned long32 runNb);
  61
  62
  63 /* Main routine
  64       Arguments:
  65       1- monitoring data source
  66 */
  67 int main(int argc, char **argv) {
  68
  69   /* magic line from Rene */
  70   gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
  71                                         "*",
  72                                         "TStreamerInfo",
  73                                         "RIO",
  74                                         "TStreamerInfo()");
  75   int status;
  76
  77   if (argc!=2) {
  78     printf("Wrong number of arguments\n");
  79     return -1;
  80   }
  81
  82
  83   /* define data source : this is argument 1 */
  84   status=monitorSetDataSource( argv[1] );
  85   if (status!=0) {
  86     printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
  87     return -1;
  88   }
  89
  90
  91   /* declare monitoring program */
  92   status=monitorDeclareMp( __FILE__ );
  93   if (status!=0) {
  94     printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
  95     return -1;
  96   }
  97
  98
  99   /* define wait event timeout - 1s max */
 100   monitorSetNowait();
 101   monitorSetNoWaitNetworkTimeout(1000);
 102
 103
 104   /* log start of process */
 105   printf("TRD DA VDRIFT started\n");
 106
 107
 108   /* init some counters */
 109   int nevents_physics=0;
 110   int nevents_total=0;
 111   //unsigned long32 runNb=0;      //run number
 112
 113   //Instance of AliCDBManager: needed by AliTRDrawFastStream
 114   //AliCDBManager *man = AliCDBManager::Instance();
 115   //man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
 116   //man->SetRun(0);
 117   //Instance of AliTRDCalibraFillHisto
 118   AliTRDCalibraFillHisto *calibra      = AliTRDCalibraFillHisto::Instance();
 119   calibra->SetNumberClustersProcent(0.9);
 120
 121   // everythings are okey for vdrift
 122   Bool_t passvdrift  = kTRUE;    // if timebin okey
 123   Int_t  nbvdrift    = 0;     // number of events with entries for vdrift
 124
 125    // setting
 126   //AliTRDrawFastStream::DisableSkipData();
 127   AliLog::SetGlobalLogLevel(AliLog::kFatal);
 128
 129   /* main loop (infinite) */
 130   for(;;) {
 131   //for(Int_t k = 0; k < 20; k++) {
 132     struct eventHeaderStruct *event;
 133     eventTypeType eventT;
 134
 135     /* check shutdown condition */
 136     if (daqDA_checkShutdown()) {break;}
 137
 138     /* get next event (blocking call until timeout) */
 139     status=monitorGetEventDynamic((void **)&event);
 140     if (status==MON_ERR_EOF) {
 141       printf ("End of File detected\n");
 142       break; /* end of monitoring file has been reached */
 143     }
 144
 145     if (status!=0) {
 146       printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
 147       break;
 148     }
 149
 150     /* retry if got no event */
 151     if (event==NULL) {
 152       continue;
 153     }
 154
 155     if( ((Int_t)nevents_total)%1000 == 0 ) printf(" event number %d (physic event number %d) will be processed\n",(Int_t) nevents_total,(Int_t) nevents_physics);
 156
 157
 158
 159     /* use event - here, just write event id to result file */
 160     eventT=event->eventType;
 161
 162     //
 163     // vdrift calibration: run only for physics events
 164     //
 165     if ((eventT==PHYSICS_EVENT) && (passvdrift)) {
 166       //if (eventT==PHYSICS_EVENT) {
 167       AliRawReader *rawReader = new AliRawReaderDate((void*)event);
 168       rawReader->Select("TRD");
 169
 170       Int_t result = calibra->ProcessEventDAQ3((AliRawReader *)rawReader);
 171       if(!result) passvdrift = kFALSE;
 172       else nbvdrift += (Int_t) result/2;
 173
 174
 175       delete rawReader;
 176
 177     }
 178
 179     // get the run number
 180     //runNb = event->eventRunNb;
 181
 182     if(eventT==PHYSICS_EVENT){
 183
 184       nevents_physics++;
 185
 186     }
 187
 188     /*
 189       if( ((Int_t)nevents_physics)%100 == 0 ) {
 190       SendToAmoreDB(&calibra);
 191       }
 192     */
 193
 194
 195     nevents_total++;
 196
 197
 198     /* free resources */
 199     free(event);
 200
 201     /* exit when last event received, no need to wait for TERM signal */
 202     if (eventT==END_OF_RUN) {
 203       printf("EOR event detected\n");
 204       break;
 205     }
 206   }
 207
 208
 209   /* report progress */
 210   printf("%d events processed and %d used\n",nevents_total,nbvdrift);
 211
 212   //
 213   // Write a local file and put it on the FX
 214   //
 215   TFile  *fileTRD     = new TFile(RESULT_FILE,"recreate");
 216
 217   if((nbvdrift > 0) && passvdrift){
 218     //Double_t *stat = calibra->StatH((TH2 *)(calibra->GetPH2d()),1);
 219     //if((stat[6] < 0.20) && (stat[5] > 3000.0)) {
 220     // write the histogram in any case to see if problems
 221     calibra->GetPH2d()->Write();
 222     //}
 223   }
 224   fileTRD->Close();
 225   status=0;
 226   // Export the file in any case to see if problems
 227   if(daqDA_FES_storeFile(RESULT_FILE,FILE_ID)) status = -2;
 228
 229   delete fileTRD;
 230
 231   return status;
 232 }
 233
 234 /*
 235   void SendToAmoreDB(TObject *o)
 236   {
 237   //AMORE
 238   const char *amoreDANameorig=gSystem->Getenv("AMORE_DA_NAME");
 239   gSystem->Setenv("AMORE_DA_NAME","TRD-dataQA");
 240
 241   amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender);
 242   Int_t statusDA=0;
 243   statusDA+=amoreDA.Send("DataQA-VDRIFT",o);
 244   if ( statusDA!=0 )
 245   printf("Waring: Failed to write one of the calib objects to the AMORE database\n");
 246   // reset env var
 247   if (amoreDANameorig) gSystem->Setenv("AMORE_DA_NAME",amoreDANameorig);
 248   }
 249 */