PHOS/PHOSLEDda.cxx

   1 /*
   2 contact: Boris.Polishchuk@cern.ch
   3 link: see comments in the $ALICE_ROOT/PHOS/AliPHOSRcuDA1.cxx
   4 reference run: /alice/data/2009/LHC09b_PHOS/000075883/raw/09000075883017.20.root
   5 run type: LED
   6 DA type: MON
   7 number of events needed: 1000
   8 input files: Mod0RCU0.data Mod0RCU1.data Mod0RCU2.data Mod0RCU3.data Mod1RCU0.data Mod1RCU1.data Mod1RCU2.data Mod1RCU3.data Mod2RCU0.data Mod2RCU1.data Mod2RCU2.data Mod2RCU3.data Mod3RCU0.data Mod3RCU1.data Mod3RCU2.data Mod3RCU3.data Mod4RCU0.data Mod4RCU1.data Mod4RCU2.data Mod4RCU3.data zs.txt
   9 Output files: PHOS_LED.root.
  10 Trigger types used: CALIBRATION_EVENT
  11 */
  12
  13
  14 #include "event.h"
  15 #include "monitor.h"
  16 extern "C" {
  17 #include "daqDA.h"
  18 }
  19
  20 #include <stdio.h>
  21 #include <stdlib.h>
  22
  23 #include <TSystem.h>
  24 #include <TROOT.h>
  25 #include <TPluginManager.h>
  26
  27 #include "AliRawReader.h"
  28 #include "AliRawReaderDate.h"
  29 #include "AliPHOSRcuDA1.h"
  30 #include "AliPHOSRawFitterv0.h"
  31 #include "AliCaloAltroMapping.h"
  32 #include "AliCaloRawStreamV3.h"
  33 #include "AliLog.h"
  34
  35
  36 /* Main routine
  37       Arguments:
  38       1- monitoring data source
  39 */
  40 int main(int argc, char **argv) {
  41
  42   gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
  43                                         "*",
  44                                         "TStreamerInfo",
  45                                         "RIO",
  46                                         "TStreamerInfo()");
  47
  48   AliLog::SetGlobalDebugLevel(0) ;
  49   AliLog::SetGlobalLogLevel(AliLog::kFatal);
  50
  51   int status;
  52
  53   if (argc!=2) {
  54     printf("Wrong number of arguments\n");
  55     return -1;
  56   }
  57
  58   /* Retrieve ZS parameters from DAQ DB */
  59   const char* zsfile = "zs.txt";
  60   int failZS = daqDA_DB_getFile(zsfile, zsfile);
  61
  62   short offset,threshold;
  63
  64   if(!failZS) {
  65     FILE *f = fopen(zsfile,"r");
  66     int scan = fscanf(f,"%d %d",&offset,&threshold);
  67   }
  68
  69   /* Retrieve mapping files from DAQ DB */
  70   const char* mapFiles[20] = {
  71     "Mod0RCU0.data",
  72     "Mod0RCU1.data",
  73     "Mod0RCU2.data",
  74     "Mod0RCU3.data",
  75     "Mod1RCU0.data",
  76     "Mod1RCU1.data",
  77     "Mod1RCU2.data",
  78     "Mod1RCU3.data",
  79     "Mod2RCU0.data",
  80     "Mod2RCU1.data",
  81     "Mod2RCU2.data",
  82     "Mod2RCU3.data",
  83     "Mod3RCU0.data",
  84     "Mod3RCU1.data",
  85     "Mod3RCU2.data",
  86     "Mod3RCU3.data",
  87     "Mod4RCU0.data",
  88     "Mod4RCU1.data",
  89     "Mod4RCU2.data",
  90     "Mod4RCU3.data"
  91   };
  92
  93   for(Int_t iFile=0; iFile<20; iFile++) {
  94     int failed = daqDA_DB_getFile(mapFiles[iFile], mapFiles[iFile]);
  95     if(failed) {
  96       printf("Cannot retrieve file %s from DAQ DB. Exit.\n",mapFiles[iFile]);
  97       return -1;
  98     }
  99   }
 100
 101   /* Open mapping files */
 102   AliAltroMapping *mapping[20];
 103   TString path = "./";
 104
 105   path += "Mod";
 106   TString path2;
 107   TString path3;
 108   Int_t iMap = 0;
 109
 110   for(Int_t iMod = 0; iMod < 5; iMod++) {
 111     path2 = path;
 112     path2 += iMod;
 113     path2 += "RCU";
 114
 115     for(Int_t iRCU=0; iRCU<4; iRCU++) {
 116       path3 = path2;
 117       path3 += iRCU;
 118       path3 += ".data";
 119       mapping[iMap] = new AliCaloAltroMapping(path3.Data());
 120       iMap++;
 121     }
 122   }
 123
 124   /* define data source : this is argument 1 */
 125   status=monitorSetDataSource( argv[1] );
 126   if (status!=0) {
 127     printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
 128     return -1;
 129   }
 130
 131
 132   /* declare monitoring program */
 133   status=monitorDeclareMp( __FILE__ );
 134   if (status!=0) {
 135     printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
 136     return -1;
 137   }
 138
 139
 140   /* define wait event timeout - 1s max */
 141   monitorSetNowait();
 142   monitorSetNoWaitNetworkTimeout(1000);
 143
 144    /* init some counters */
 145   int nevents_physics=0;
 146   int nevents_total=0;
 147
 148   AliRawReader *rawReader = NULL;
 149   AliPHOSRcuDA1* dAs[5];
 150
 151   for(Int_t iMod=0; iMod<5; iMod++) {
 152     dAs[iMod] = 0;
 153   }
 154
 155   Float_t e[64][56][2];
 156   Float_t t[64][56][2];
 157
 158   for(Int_t iX=0; iX<64; iX++) {
 159     for(Int_t iZ=0; iZ<56; iZ++) {
 160       for(Int_t iGain=0; iGain<2; iGain++) {
 161         e[iX][iZ][iGain] = 0.;
 162         t[iX][iZ][iGain] = 0.;
 163       }
 164     }
 165   }
 166
 167   Int_t cellX    = -1;
 168   Int_t cellZ    = -1;
 169   Int_t nBunches =  0;
 170   Int_t nFired   = -1;
 171   Int_t sigStart, sigLength;
 172   Int_t caloFlag;
 173
 174
 175   TH1I fFiredCells("fFiredCells","Number of fired cells per event",100,0,1000);
 176
 177   /* main loop (infinite) */
 178   for(;;) {
 179     struct eventHeaderStruct *event;
 180     eventTypeType eventT;
 181
 182     /* check shutdown condition */
 183     if (daqDA_checkShutdown()) {break;}
 184
 185     /* get next event (blocking call until timeout) */
 186     status=monitorGetEventDynamic((void **)&event);
 187     if (status==MON_ERR_EOF) {
 188       printf ("End of File detected\n");
 189       break; /* end of monitoring file has been reached */
 190     }
 191
 192     if (status!=0) {
 193       printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
 194       break;
 195     }
 196
 197     /* retry if got no event */
 198     if (event==NULL) {
 199       continue;
 200     }
 201
 202
 203     /* use event - here, just write event id to result file */
 204     eventT=event->eventType;
 205
 206     if (eventT==PHYSICS_EVENT) {
 207
 208       nFired = 0;
 209
 210       rawReader = new AliRawReaderDate((void*)event);
 211       AliCaloRawStreamV3 stream(rawReader,"PHOS",mapping);
 212       AliPHOSRawFitterv0 fitter;
 213       fitter.SubtractPedestals(kTRUE); // assume that data is non-ZS
 214
 215       if(!failZS) {
 216         fitter.SubtractPedestals(kFALSE);
 217         fitter.SetAmpOffset(offset);
 218         fitter.SetAmpThreshold(threshold);
 219       }
 220
 221       while (stream.NextDDL()) {
 222         while (stream.NextChannel()) {
 223
 224           /* Retrieve ZS parameters from data*/
 225           if(failZS) {
 226             short value = stream.GetAltroCFG1();
 227             bool ZeroSuppressionEnabled = (value >> 15) & 0x1;
 228             bool AutomaticBaselineSubtraction = (value >> 14) & 0x1;
 229             if(ZeroSuppressionEnabled) {
 230               offset = (value >> 10) & 0xf;
 231               threshold = value & 0x3ff;
 232               fitter.SubtractPedestals(kFALSE);
 233               fitter.SetAmpOffset(offset);
 234               fitter.SetAmpThreshold(threshold);
 235             }
 236           }
 237
 238           cellX    = stream.GetCellX();
 239           cellZ    = stream.GetCellZ();
 240           caloFlag = stream.GetCaloFlag();  // 0=LG, 1=HG, 2=TRU
 241
 242           if(caloFlag!=0 && caloFlag!=1) continue; //TRU data!
 243
 244           // In case of oscillating signals with ZS, a channel can have several bunches
 245           nBunches = 0;
 246           while (stream.NextBunch()) {
 247             nBunches++;
 248             if (nBunches > 1) continue;
 249             sigStart  = stream.GetStartTimeBin();
 250             sigLength = stream.GetBunchLength();
 251             fitter.SetChannelGeo(stream.GetModule(),cellX,cellZ,caloFlag);
 252             fitter.Eval(stream.GetSignals(),sigStart,sigLength);
 253           } // End of NextBunch()
 254
 255           if(nBunches != 1) continue;
 256
 257           e[cellX][cellZ][caloFlag] = fitter.GetEnergy();
 258           t[cellX][cellZ][caloFlag] = fitter.GetTime();
 259
 260           if(caloFlag==1 && fitter.GetEnergy()>40)
 261             nFired++;
 262         }
 263
 264         if(stream.GetModule()<0 || stream.GetModule()>4) continue;
 265
 266         if(dAs[stream.GetModule()])
 267           dAs[stream.GetModule()]->FillHistograms(e,t);
 268         else {
 269           dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1,0);
 270           dAs[stream.GetModule()]->FillHistograms(e,t);
 271         }
 272
 273         for(Int_t iX=0; iX<64; iX++) {
 274           for(Int_t iZ=0; iZ<56; iZ++) {
 275             for(Int_t iGain=0; iGain<2; iGain++) {
 276               e[iX][iZ][iGain] = 0.;
 277               t[iX][iZ][iGain] = 0.;
 278             }
 279           }
 280         }
 281
 282       }
 283
 284       fFiredCells.Fill(nFired);
 285
 286       delete rawReader;
 287       nevents_physics++;
 288     }
 289
 290     nevents_total++;
 291
 292     /* free resources */
 293     free(event);
 294
 295     /* exit when last event received, no need to wait for TERM signal */
 296     if (eventT==END_OF_RUN) {
 297       printf("EOR event detected\n");
 298       break;
 299     }
 300   }
 301
 302   for(Int_t i = 0; i < 20; i++) delete mapping[i];
 303
 304   /* Be sure that all histograms are saved */
 305
 306   const TH2F* h2=0;
 307   const TH1F* h1=0;
 308   char localfile[128];
 309
 310   Int_t nGood=0;    // >10 entries in peak
 311   Int_t nMax=-111;  // max. number of entries in peak
 312   Int_t iXmax=-1;
 313   Int_t iZmax=-1;
 314   Int_t iModMax=-1;
 315
 316   sprintf(localfile,"PHOS_LED.root");
 317   TFile* f = new TFile(localfile,"recreate");
 318
 319   for(Int_t iMod=0; iMod<5; iMod++) {
 320     if(!dAs[iMod]) continue;
 321
 322     printf("DA1 for module %d detected.\n",iMod);
 323
 324     for(Int_t iX=0; iX<64; iX++) {
 325       for(Int_t iZ=0; iZ<56; iZ++) {
 326
 327         h1 = dAs[iMod]->GetHgLgRatioHistogram(iX,iZ); // High Gain/Low Gain ratio
 328         if(h1) {
 329           if(h1->GetMaximum()>10.) nGood++;
 330           if(h1->GetMaximum()>nMax) {
 331             nMax = (Int_t)h1->GetMaximum(); iXmax=iX; iZmax=iZ; iModMax=iMod;
 332           }
 333           h1->Write();
 334         }
 335
 336         for(Int_t iGain=0; iGain<2; iGain++) {
 337           h2 = dAs[iMod]->GetTimeEnergyHistogram(iX,iZ,iGain); // Time vs Energy
 338           if(h2) h2->Write();
 339         }
 340
 341       }
 342     }
 343
 344     fFiredCells.Write();
 345   }
 346
 347   f->Close();
 348
 349   /* Store output files to the File Exchange Server */
 350   daqDA_FES_storeFile(localfile,"LED");
 351
 352   printf("%d physics events of %d total processed.\n",nevents_physics,nevents_total);
 353   printf("%d histograms has >10 entries in maximum, max. is %d entries ",nGood,nMax);
 354   printf("(module,iX,iZ)=(%d,%d,%d)",iModMax,iXmax,iZmax);
 355   printf("\n");
 356
 357   return status;
 358 }