TOF/TOFpulserda.cxx

   1 /*
   2
   3 TOF DA for online calibration from pulser data
   4
   5 Contact: Chiara.Zampolli@bo.infn.it
   6 Link: www.bo.infn.it/~zampolli
   7 Run Type: CALIBRATION, PULSER - to be properly defined
   8 DA Type: LDC
   9 Number of events needed: 10000
  10 Input Files: TOF<nrun>.raw, where <nrun> is the run number
  11 Output Files: TOFoutPulserLDC_<LDCid>.root, where <LDCid> is the id of the LDC which is read (2 characters field, e.g. TOFoutPulserLDC_03.root),to be exported to the DAQ FXS
  12 Trigger types used: PHYSICS_EVENT, for the time being
  13
  14 */
  15
  16 // DATE
  17 #include "event.h"
  18 #include "monitor.h"
  19 #include "daqDA.h"
  20
  21 #include <stdio.h>
  22 #include <stdlib.h>
  23
  24 //AliRoot
  25 #include <AliTOFRawStream.h>
  26 #include <AliRawReaderDate.h>
  27 #include <AliRawReader.h>
  28 #include <AliTOFGeometry.h>
  29 #include <AliDAQ.h>
  30 #include <AliTOFHitData.h>
  31 #include <AliTOFHitDataBuffer.h>
  32 #include <AliTOFDecoder.h>
  33
  34 //ROOT
  35 #include <TFile.h>
  36 #include <TKey.h>
  37 #include <TH1I.h>
  38 #include <TObject.h>
  39 #include <TMath.h>
  40 #include <TSystem.h>
  41
  42 /* Main routine
  43       Arguments: list of DATE raw data files
  44 */
  45 int main(int argc, char **argv) {
  46
  47   AliTOFGeometry * geom = new AliTOFGeometry();
  48
  49   static const Int_t size = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
  50   TH1F::AddDirectory(0);
  51   TH1I * htofPulser = new TH1I("hTOFpulser","histo with signals on TOF during pulser", size,-0.5,size-0.5);
  52   for (Int_t ibin =1;ibin<=size;ibin++){
  53     htofPulser->SetBinContent(ibin,-1);
  54   }
  55
  56   UInt_t ldcId=0;
  57
  58   int status;
  59
  60   /* log start of process */
  61   printf("TOF DA started\n");
  62
  63   /* check that we got some arguments = list of files */
  64   if (argc<2) {
  65     printf("Wrong number of arguments\n");
  66     return -1;
  67   }
  68
  69   /* open result file */
  70   FILE *fp=NULL;
  71   fp=fopen("./result.txt","a");
  72   if (fp==NULL) {
  73     printf("Failed to open file\n");
  74     return -1;
  75   }
  76
  77   /* init some counters */
  78   int nevents_physics=0;
  79   int nevents_total=0;
  80
  81   /* read the data files */
  82   int n;
  83   for (n=1;n<argc;n++) {
  84
  85     status=monitorSetDataSource( argv[n] );
  86     if (status!=0) {
  87       printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
  88       return -1;
  89     }
  90
  91     /* read the file */
  92     for(;;) {
  93       struct eventHeaderStruct *event;
  94       eventTypeType eventT;
  95
  96       /* get next event */
  97       status=monitorGetEventDynamic((void **)&event);
  98       if (status==MON_ERR_EOF) break; /* end of monitoring file has been reached */
  99       if (status!=0) {
 100         printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
 101         return -1;
 102       }
 103
 104       /* retry if got no event */
 105       if (event==NULL) {
 106         break;
 107       }
 108
 109       /* use event - here, just write event id to result file */
 110       eventT=event->eventType;
 111
 112       if (eventT==PHYSICS_EVENT) {
 113         //printf ("event %i \n", nevents_physics);
 114         Int_t nPDBEntriesToT = 0;
 115         Int_t nDBEntriesToT = 0;
 116         AliTOFHitData *HitData;
 117         Int_t dummy = -1;
 118         Int_t Volume[5];
 119         AliRawReader *rawReader = new AliRawReaderDate((void*)event);
 120         AliTOFRawStream *rawStreamTOF = new AliTOFRawStream(rawReader);
 121         AliTOFDecoder * decoderTOF = new AliTOFDecoder();
 122         AliTOFHitDataBuffer *DataBuffer[72];
 123         AliTOFHitDataBuffer *PackedDataBuffer[72];
 124         for (Int_t i=0;i<AliDAQ::NumberOfDdls("TOF");i++){
 125           DataBuffer[i]=new AliTOFHitDataBuffer();
 126           PackedDataBuffer[i]=new AliTOFHitDataBuffer();
 127         }
 128         Int_t currentEquipment;
 129         Int_t currentDDL;
 130         UChar_t *data = 0x0;
 131         while(rawReader->ReadHeader()){
 132           ldcId = rawReader->GetLDCId();
 133           //printf ("ldcId = %i \n",ldcId);
 134
 135           //memory leak prevention (actually data should be always 0x0 here)
 136           if (data != 0x0)delete [] data;
 137
 138           //get equipment infos
 139           currentEquipment = rawReader->GetEquipmentId();
 140           currentDDL = rawReader->GetDDLID();
 141           Int_t nchDDL = 0;
 142           if (currentDDL%2==0) {
 143             nchDDL = 2160;
 144           }
 145           else {
 146             nchDDL = 2208;
 147           }
 148           Int_t * array = new Int_t[nchDDL];
 149           decoderTOF->GetArrayDDL(array, currentDDL);
 150
 151           for (Int_t i=0;i<nchDDL;i++){
 152             if (htofPulser->GetBinContent(array[i]+1)<0) htofPulser->SetBinContent(array[i]+1,0);
 153           }
 154
 155           //printf(" Equipment = %i, and DDL = %i \n", currentEquipment,currentDDL);
 156           const Int_t kDataSize = rawReader->GetDataSize();
 157           const Int_t kDataWords = kDataSize / 4;
 158           data = new UChar_t[kDataSize];
 159           decoderTOF->SetDataBuffer(DataBuffer[currentDDL]);
 160           decoderTOF->SetPackedDataBuffer(PackedDataBuffer[currentDDL]);
 161           //start decoding
 162           if (!rawReader->ReadNext(data, kDataSize)){
 163             rawReader->AddMajorErrorLog(AliTOFRawStream::kDDLdataReading);
 164             printf("Error while reading DDL data. Go to next equipment \n");
 165             delete [] data;
 166             data = 0x0;
 167             continue;
 168           }
 169           if (decoderTOF->Decode((UInt_t *)data, kDataWords) == kTRUE) {
 170             rawReader->AddMajorErrorLog(AliTOFRawStream::kDDLDecoder,Form("DDL # = %d",currentDDL));
 171             printf("Error while decoding DDL # %d: decoder returned with errors \n", currentDDL);
 172           }
 173
 174           Int_t nDBEntries = DataBuffer[currentDDL]->GetEntries();
 175           Int_t nPDBEntries = PackedDataBuffer[currentDDL]->GetEntries();
 176           nPDBEntriesToT+=nPDBEntries;
 177           nDBEntriesToT+=nDBEntries;
 178           /* reset buffer */
 179           DataBuffer[currentDDL]->Reset();
 180
 181           /* read data buffer hits */
 182           for (Int_t iHit = 0; iHit < nPDBEntries; iHit++){
 183             HitData = PackedDataBuffer[currentDDL]->GetHit(iHit);
 184             /* add volume information */
 185             HitData->SetDDLID(currentDDL);
 186             rawStreamTOF->EquipmentId2VolumeId(HitData, Volume);
 187             if (Volume[0]==-1 ||
 188                 Volume[1]==-1 ||
 189                 Volume[2]==-1 ||
 190                 Volume[3]==-1 ||
 191                 Volume[4]==-1) continue;
 192             else {
 193               dummy = Volume[3];
 194               Volume[3] = Volume[4];
 195               Volume[4] = dummy;
 196               Int_t index = geom->GetIndex(Volume);
 197               //printf ("index = %i \n",index);
 198               htofPulser->Fill(index); //channel index start from 0, bin index from 1
 199               //debugging printings
 200               //printf("sector %i, plate %i, strip %i, padz %i, padx %i \n",Volume[0],Volume[1],Volume[2],Volume[3],Volume[4]);
 201             }
 202           }
 203           /* reset buffer */
 204           PackedDataBuffer[currentDDL]->Reset();
 205           delete [] data;
 206           data = 0x0;
 207           delete [] array;
 208         }
 209
 210         //printf(" Packed Hit Buffer Entries = %i \n",nPDBEntriesToT);
 211         //printf(" Hit Buffer Entries = %i \n",nDBEntriesToT);
 212
 213         delete decoderTOF;
 214         decoderTOF=0x0;
 215
 216         for (Int_t i=0;i<72;i++){
 217           delete DataBuffer[i];
 218           delete PackedDataBuffer[i];
 219         }
 220
 221         delete rawStreamTOF;
 222         rawStreamTOF = 0x0;
 223
 224         delete rawReader;
 225         rawReader = 0x0;
 226         nevents_physics++;
 227       }
 228       nevents_total++;
 229
 230       /* free resources */
 231       free(event);
 232     }
 233   }
 234
 235   delete geom;
 236   geom = 0x0;
 237
 238   //write the Run level file
 239   char filename[100];
 240   sprintf(filename,"TOFoutPulserLDC_%02i.root",ldcId);
 241   TFile * fileRun = new TFile (filename,"RECREATE");
 242   htofPulser->Write();
 243   fileRun->Close();
 244
 245   /* write report */
 246   fprintf(fp,"Run #%s, received %d physics events out of %d\n",getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total);
 247
 248   /* close result file */
 249   fclose(fp);
 250
 251   /* store the result file on FES */
 252   status=daqDA_FES_storeFile(filename,filename);
 253   if (status) {
 254     status = -2;
 255   }
 256
 257   return status;
 258 }