PWGPP/CalibMacros/CPass1/makeOCDB.C

   1 /*
   2   macro to extract the OCDB entries
   3
   4   input: CalibObjects.root
   5   ouput: TimeGain and TimeVdrift calibration objects for TPC and TRD
   6
   7   Example:
   8   .L $ALICE_ROOT/PWGPP/CalibMacros/CPass1/makeOCDB.C
   9   makeOCDB("105160");
  10
  11 */
  12
  13 void printCalibStat(Int_t run, const char * fname,  TTreeSRedirector * pcstream);
  14
  15 void makeOCDB(Int_t runNumber, TString  targetOCDBstorage="", TString sourceOCDBstorage="raw://")
  16 {
  17   //
  18   // extract OCDB entries for detectors participating in the calibration for the current run
  19   //
  20
  21   gROOT->Macro("$ALICE_ROOT/PWGPP/CalibMacros/CPass1/LoadLibraries.C");
  22   gROOT->LoadMacro("$ALICE_ROOT/PWGPP/CalibMacros/CPass1/ConfigCalibTrain.C");
  23
  24   // switch off log info
  25   AliLog::SetClassDebugLevel("AliESDEvent",0);
  26
  27   // config GRP
  28   printf("runNumber from runCalibTrain = %d\n",runNumber);
  29   ConfigCalibTrain(runNumber, sourceOCDBstorage.Data());
  30
  31   // check the presence of the detectors
  32   AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/GRP/Data");
  33   AliGRPObject* grpData = dynamic_cast<AliGRPObject*>(entry->GetObject());
  34   if (!grpData) {printf("Failed to get GRP data for run",runNumber); return;}
  35   Int_t activeDetectors = grpData->GetDetectorMask();
  36   TString detStr = AliDAQ::ListOfTriggeredDetectors(activeDetectors);
  37   printf("Detectors in the data:\n%s\n",detStr.Data());
  38
  39   // Steering Tasks - set output storage
  40   // DefaultStorage set already before - in ConfigCalibTrain.C
  41   //targetOCDBstorage+="?se=ALICE::CERN::SE";
  42
  43   // for tests, the target OCDB has to be different from raw://, but since the calibration procedure may need to read the OCDB,
  44   // one cannot set the target OCDB as new default storage, since it could it that it does not contain all the necessary entries.
  45   // This is true only if the target OCDB storage is different from the one used as default source and to configure the train
  46   if (targetOCDBstorage != sourceOCDBstorage) AliCDBManager::Instance()->SetSpecificStorage("*/*/*", targetOCDBstorage.Data());
  47
  48   if (gSystem->AccessPathName("TPC", kFileExists)==0) {
  49     AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/Correction","local://");
  50   }
  51   // set OCDB storage
  52   if (targetOCDBstorage.Length()==0) targetOCDBstorage+="local://"+gSystem->GetFromPipe("pwd")+"/OCDB";
  53   TString ocdbLocal = "local://"+gSystem->GetFromPipe("pwd")+"/OCDB";
  54
  55   // TPC part
  56   AliTPCPreprocessorOffline* procesTPC=0;
  57   if (detStr.Contains("TPC")){
  58     Printf("\n******* Calibrating TPC *******");
  59     procesTPC = new AliTPCPreprocessorOffline;
  60     // switch on parameter validation
  61     procesTPC->SetTimeGainRange(0.5,4.0);
  62     procesTPC->SwitchOnValidation();
  63     // Make timegain calibration
  64     //proces.CalibTimeGain("CalibObjects.root", runNumber,AliCDBRunRange::Infinity(),targetOCDBstorage);
  65     //procesTPC->CalibTimeGain("CalibObjects.root", runNumber,runNumber,ocdbLocal);
  66     // Make vdrift calibration
  67     //proces.CalibTimeVdrift("CalibObjects.root",runNumber,AliCDBRunRange::Infinity(),targetOCDBstorage);
  68     //procesTPC->CalibTimeVdrift("CalibObjects.root",runNumber,runNumber,ocdbLocal);
  69   }
  70
  71   // TOF part
  72   AliTOFAnalysisTaskCalibPass0 *procesTOF=0;
  73   if ( detStr.Contains("TOF") && detStr.Contains("TPC")){
  74     procesTOF = new AliTOFAnalysisTaskCalibPass0;
  75     Printf("\n******* Calibrating TOF *******");
  76     procesTOF->ProcessOutput("CalibObjects.root", targetOCDBstorage);
  77   }
  78
  79   // T0 part
  80   AliT0PreprocessorOffline *procesT0= 0;
  81   if ( detStr.Contains("T0")) {
  82     Printf("\n******* Calibrating T0 *******");
  83     procesT0 = new AliT0PreprocessorOffline;
  84     // Make  calibration of channels offset
  85     procesT0->setDArun(177000);
  86     procesT0->Process("CalibObjects.root",runNumber, runNumber, targetOCDBstorage);
  87   }
  88
  89   //TRD part
  90   AliTRDPreprocessorOffline *procesTRD = 0;
  91   if ( detStr.Contains("TRD")){
  92     Printf("\n******* Calibrating TRD *******");
  93     procesTRD = new  AliTRDPreprocessorOffline;
  94     procesTRD->SetLinearFitForVdrift(kTRUE);
  95     procesTRD->SetMinStatsVdriftT0PH(600*10);
  96     procesTRD->SetMinStatsVdriftLinear(50);
  97     procesTRD->SetMinStatsGain(600);
  98     procesTRD->SetLimitValidateNoData(60);
  99     procesTRD->SetLimitValidateBadCalib(60);
 100     procesTRD->SetAlternativeDriftVelocityFit(kTRUE);
 101     procesTRD->Init("CalibObjects.root");
 102     Int_t versionVdriftUsed = procesTRD->GetVersionVdriftUsed();
 103     Int_t subversionVdriftUsed = procesTRD->GetSubVersionVdriftUsed();
 104     Int_t versionGainUsed = procesTRD->GetVersionGainUsed();
 105     Int_t subversionGainUsed = procesTRD->GetSubVersionGainUsed();
 106     Int_t versionExBUsed = procesTRD->GetVersionExBUsed();
 107     Int_t subversionExBUsed = procesTRD->GetSubVersionExBUsed();
 108     printf("version and subversion vdrift %d and %d\n",versionVdriftUsed,subversionVdriftUsed);
 109     printf("version and subversion gain %d and %d\n",versionGainUsed,subversionGainUsed);
 110     printf("version and subversion exb %d and %d\n",versionExBUsed,subversionExBUsed);
 111     procesTRD->Process("CalibObjects.root",runNumber,runNumber,targetOCDBstorage);
 112   }
 113
 114   // switched OFF at CPass1 in any case
 115   /*
 116   //Mean Vertex
 117   AliMeanVertexPreprocessorOffline * procesMeanVtx=0;
 118   if ( detStr.Contains("ITSSPD")) {
 119     Printf("\n******* Calibrating MeanVertex *******");
 120     procesMeanVtx = new AliMeanVertexPreprocessorOffline;
 121     procesMeanVtx->ProcessOutput("CalibObjects.root", targetOCDBstorage, runNumber);
 122   }
 123   */
 124
 125   //
 126   // Print calibration status into the stdout
 127   //
 128   Int_t trdStatus = (procesTRD) ?  procesTRD->GetStatus():0;
 129   Int_t tofStatus = (procesTOF) ?  procesTOF->GetStatus():0;
 130   Int_t t0Status  = (procesT0)  ?  procesT0->GetStatus():0;
 131   Int_t tpcStatus = (procesTPC) ? ((procesTPC->ValidateTimeDrift() || procesTPC->ValidateTimeGain())==kFALSE):0;
 132   //
 133   printf("\n\n\n\n");
 134   printf("CPass1 calibration status\n");
 135   printf("TRD calibration status=%d\n",trdStatus);
 136   printf("TOF calibration status=%d\n",tofStatus);
 137   printf("TPC calibration status=%d\n",tpcStatus);
 138   printf("T0  calibration status=%d\n",t0Status);
 139
 140   TTreeSRedirector *pcstream = new TTreeSRedirector("cpassStat.root","recreate");
 141   printCalibStat(runNumber, "CalibObjects.root",pcstream);
 142   delete pcstream;
 143   return;
 144 }
 145
 146
 147 // function to print statistics used to calibrate the various detectors
 148
 149 void printCalibStat(Int_t run, const char * fname,  TTreeSRedirector * pcstream){
 150
 151   //
 152   // Dump the statistical information about all histograms in the calibration files
 153   //    into the statistical tree, print on the screen (log files) as well
 154   //
 155   //
 156   // 1. Default dump for all histograms
 157   //    Information to dump:
 158   //    stat =Entries, Mean, MeanError,  RMS, MaxBin
 159   //    Branch naming convention:
 160   //    <detName>_<hisName><statName>
 161   //
 162   // 2. Detector statistical information  - to be implemented by expert
 163   //                                      - First version implemented by MI
 164   //
 165   //
 166
 167   TFile *fin = TFile::Open(fname);
 168   if (!fin) return;
 169   const Double_t kMaxHis=10000;
 170
 171   TList * keyList = fin->GetListOfKeys();
 172   Int_t nkeys=keyList->GetEntries();
 173   Double_t *hisEntries = new Double_t[kMaxHis];
 174   Double_t *hisMean = new Double_t[kMaxHis];
 175   Double_t *hisMeanError = new Double_t[kMaxHis];
 176   Double_t *hisRMS = new Double_t[kMaxHis];
 177   Double_t *hisMaxBin = new Double_t[kMaxHis];
 178   Int_t counter=0;
 179
 180   if (pcstream) (*pcstream)<<"calibStatAll"<<"run="<<run;
 181   for (Int_t ikey=0; ikey<nkeys; ikey++){
 182     TObject * object = fin->Get(keyList->At(ikey)->GetName());
 183     if (!object) continue;
 184     if (object->InheritsFrom("TCollection")==0) continue;
 185     TSeqCollection *collection  = (TSeqCollection*)object;
 186     Int_t nentries= collection->GetEntries();
 187     for (Int_t ihis=0; ihis<nentries; ihis++){
 188       TObject * ohis = collection->At(ihis);
 189       if (!ohis) continue;
 190       if (ohis->InheritsFrom("TH1")==0) continue;
 191       TH1* phis = (TH1*)ohis;
 192       hisEntries[counter]=phis->GetEntries();
 193       Int_t idim=1;
 194       if (ohis->InheritsFrom("TH2")) idim=2;
 195       if (ohis->InheritsFrom("TH3")) idim=3;
 196       hisMean[counter]=phis->GetMean(idim);
 197       hisMeanError[counter]=phis->GetMeanError(idim);
 198       hisRMS[counter]=phis->GetRMS(idim);
 199       hisMaxBin[counter]=phis->GetBinCenter(phis->GetMaximumBin());
 200       if (pcstream) (*pcstream)<<"calibStatAll"<<
 201                       Form("%s_%sEntries=",keyList->At(ikey)->GetName(), phis->GetName())<<hisEntries[counter]<<
 202                       Form("%s_%sMean=",keyList->At(ikey)->GetName(), phis->GetName())<<hisMean[counter]<<
 203                       Form("%s_%sMeanError=",keyList->At(ikey)->GetName(), phis->GetName())<<hisMeanError[counter]<<
 204                       Form("%s_%sRMS=",keyList->At(ikey)->GetName(), phis->GetName())<<hisRMS[counter]<<
 205                       Form("%s_%sMaxBin=",keyList->At(ikey)->GetName(), phis->GetName())<<hisMaxBin[counter];
 206       //printf("Histo:\t%s_%s\t%f\t%d\n",keyList->At(ikey)->GetName(), phis->GetName(), hisEntries[counter],idim);
 207       counter++;
 208     }
 209     delete object;
 210   }
 211
 212   //
 213   // Expert dump example (MI first iteration):
 214   //
 215   // 0.)  TOF dump
 216   //
 217
 218   Int_t tofEvents=0;
 219   Int_t tofTracks=0;
 220   TList * TOFCalib = (TList*)fin->Get("TOFHistos");
 221   if (TOFCalib) {
 222     TH1 *histoEvents = (TH1*)TOFCalib->FindObject("hHistoVertexTimestamp");
 223     TH1 *histoTracks = (TH1*)TOFCalib->FindObject("hHistoDeltatTimestamp");
 224     if (histoEvents && histoTracks){
 225       tofEvents = TMath::Nint(histoEvents->GetEntries());
 226       tofTracks = TMath::Nint(histoTracks->GetEntries());
 227     }
 228     delete TOFCalib;
 229   }
 230   printf("Monalisa TOFevents\t%d\n",tofEvents);
 231   if (pcstream) (*pcstream)<<"calibStatAll"<<"TOFevents="<<tofEvents;
 232   printf("Monalisa TOFtracks\t%d\n",tofTracks);
 233   if (pcstream) (*pcstream)<<"calibStatAll"<<"TOFtracks="<<tofTracks;
 234
 235   //
 236   // 1.)  TPC  dump - usefull events/tracks  for the calibration
 237   //
 238   Int_t tpcEvents=0;
 239   Int_t tpcTracks=0;
 240   TList * TPCCalib = (TList*)fin->Get("TPCCalib");
 241   if (TPCCalib) {
 242     AliTPCcalibTime  * calibTime = (AliTPCcalibTime  *)TPCCalib->FindObject("calibTime");
 243     if (calibTime){
 244       tpcEvents = TMath::Nint(calibTime->GetTPCVertexHisto(0)->GetEntries());
 245       tpcTracks = TMath::Nint(calibTime->GetResHistoTPCITS(0)->GetEntries());
 246     }
 247   }
 248   printf("Monalisa TPCevents\t%d\n",tpcEvents);
 249   if (pcstream) (*pcstream)<<"calibStatAll"<<"TPCevents="<<tpcEvents;
 250   printf("Monalisa TPCtracks\t%d\n",tpcTracks);
 251   if (pcstream) (*pcstream)<<"calibStatAll"<<"TPCtracks="<<tpcTracks;
 252
 253   //
 254   // 2. TRD dump
 255   //
 256   Int_t trdEvents=0;
 257   Int_t trdTracks=0;
 258   TList * TRDCalib = (TList*)fin->Get("TRDCalib");
 259   if (TRDCalib) {
 260     TH1  *histoEvents = (TH1*)TRDCalib->FindObject("NEventsInput_AliTRDCalibTask");
 261     TH1  *histoTracks = (TH1*)TRDCalib->FindObject("AbsoluteGain_AliTRDCalibTask");
 262     if (histoEvents && histoTracks){
 263       trdEvents= TMath::Nint(histoEvents->GetEntries());
 264       trdTracks= TMath::Nint(histoTracks->GetEntries());
 265     }
 266     delete TRDCalib;
 267   }
 268   printf("Monalisa TRDevents\t%d\n",trdEvents);
 269   if (pcstream) (*pcstream)<<"calibStatAll"<<"TRDevents="<<trdEvents;
 270   printf("Monalisa TRDtracks\t%d\n",trdTracks);
 271   if (pcstream) (*pcstream)<<"calibStatAll"<<"TRDtracks="<<trdTracks;
 272
 273   //
 274   // 3. T0 dump
 275   //
 276   Int_t T0Events=0;
 277   TList * T0Calib = (TList*)fin->Get("T0Calib");
 278   if (T0Calib) {
 279     TH1  *histoEvents = (TH1*) T0Calib->FindObject("fTzeroORAplusORC");
 280     if (histoEvents){
 281       T0Events= TMath::Nint(histoEvents->GetEntries());
 282     }
 283     delete T0Calib;
 284   }
 285   printf("Monalisa T0events\t%d\n",T0Events);
 286   if (pcstream) (*pcstream)<<"calibStatAll"<<"T0events="<<T0Events;
 287
 288   //
 289   // 4. Mean vertex -   dump
 290   // Not present in CPass1
 291   /*
 292     Int_t meanVertexEvents=0;
 293   TList * meanVertexCalib = (TList*)fin->Get("MeanVertex");
 294   if (meanVertexCalib) {
 295     TH1  *histoEvents = (TH1*) meanVertexCalib->FindObject("hTRKVertexX");
 296     if (histoEvents){
 297       meanVertexEvents = TMath::Nint(histoEvents->GetEntries());
 298     }
 299     delete meanVertexCalib;
 300   }
 301   printf("Monalisa MeanVertexevents\t%d\n",meanVertexEvents);
 302   if (pcstream) (*pcstream)<<"calibStatAll"<<"MeanVertexevents="<<meanVertexEvents;
 303   */
 304
 305   //
 306   // 5. SDD dump
 307   //
 308   Int_t sddEvents=0;
 309   Int_t sddTracks=0;
 310   TList * SDDCalib = (TList*)fin->Get("clistSDDCalib");
 311   if (SDDCalib) {
 312     TH1  *histoEvents = (TH1*) SDDCalib->FindObject("hNEvents");
 313     if (histoEvents ){
 314       sddEvents = TMath::Nint(histoEvents->GetBinContent(4));
 315       sddTracks = TMath::Nint(histoEvents->GetBinContent(5));
 316     }
 317     delete SDDCalib;
 318   }
 319   printf("Monalisa SDDevents\t%d\n",sddEvents);
 320   if (pcstream) (*pcstream)<<"calibStatAll"<<"SDDevents="<<sddEvents;
 321   printf("Monalisa SDDtracks\t%d\n",sddTracks);
 322   if (pcstream) (*pcstream)<<"calibStatAll"<<"SDDtracks="<<sddTracks;
 323
 324   //
 325   if (pcstream) (*pcstream)<<"calibStatAll"<<"\n";
 326   delete fin;
 327
 328 }
 329
 330