PWG1/CalibMacros/Pass0/AddTaskT0Calib.C

   1
   2
   3 void    readCDB (TObject *task1);
   4 //_____________________________________________________________________________
   5 AliAnalysisTask  *AddTaskT0Calib(Int_t runNumber)
   6 {
   7   //
   8   // add calibration task
   9   //
  10   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  11   if (!mgr) {
  12     ::Error("AddTaskT0Calib", "No analysis manager to connect to.");
  13     return NULL;
  14   }
  15
  16   // check the input handler
  17   if (!mgr->GetInputEventHandler()) {
  18     ::Error("AddTaskT0Calib", "This task requires an input event handler");
  19     return NULL;
  20   }
  21
  22   // set TPC OCDB parameters
  23   //ConfigOCDB(runNumber);
  24
  25   // setup task
  26   AliT0CalibOffsetChannelsTask  *task1 = new AliT0CalibOffsetChannelsTask("CalibObjectsTrain1");
  27   readCDB(task1);
  28   mgr->AddTask(task1);
  29
  30   //  AliT0AnalysisTaskQA * task2 = new AliT0AnalysisTaskQA("QA task");
  31   //    mgr->AddTask(task2);
  32
  33   AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
  34   if (!cinput1) cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
  35                                       AliAnalysisManager::kInputContainer);
  36   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("T0Calib",TObjArray::Class(), AliAnalysisManager::kOutputContainer, "AliESDfriends_v1.root");
  37
  38   mgr->ConnectInput(task1,0,cinput1);
  39   mgr->ConnectOutput(task1,1,coutput1);
  40   return task1;
  41 }
  42 //_____________________________________________________________________________
  43 void    readCDB (TObject *task1) {
  44
  45   Float_t zero_timecdb[24]={0};
  46   Float_t *timecdb = zero_timecdb;
  47   Float_t cfdvalue[24][5];
  48   Float_t zero_shiftcdb[4]={0};
  49   Float_t *shiftcdb = zero_shiftcdb;
  50   AliT0CalibOffsetChannelsTask *mytask = (AliT0CalibOffsetChannelsTask*)task1;
  51
  52   //  AliCDBManager* man = AliCDBManager::Instance();
  53   //  man->SetDefaultStorage("raw://");
  54   // man->SetRun(158124);
  55   AliCDBEntry *entryCalib1 = man->Get("T0/Calib/TimeDelay");
  56   if(!entryCalib1) {
  57     AliError::(Form("Cannot find any AliCDBEntry for [Calib, TimeDelay]!"));
  58     return;
  59   }
  60   else
  61     {
  62       AliT0CalibTimeEq *clb = (AliT0CalibTimeEq*)entryCalib1->GetObject();
  63       timecdb = clb->GetTimeEq();
  64       for(Int_t i=0; i<24; i++)
  65         for (Int_t i0=0; i0<5; i0++){
  66           cfdvalue[i][i0] = clb->GetCFDvalue(i, i0);
  67         }
  68     }
  69   for (Int_t i=0; i<24; i++) {
  70     mytask->SetCFDvalue(i,cfdvalue[i][0]);
  71     mytask->SetTimeEq(i,timecdb[i]);
  72   }
  73
  74   AliCDBEntry *entryCalib2 = man->Get("T0/Calib/TimeAdjust");
  75   if(!entryCalib2) {
  76      AliError(Form("Cannot find any AliCDBEntry for [Calib, TimeAdjust]!"));
  77   }
  78  else
  79     {
  80       AliT0CalibSeasonTimeShift *clb1 = (AliT0CalibSeasonTimeShift*)entryCalib2->GetObject();
  81       shiftcdb = clb1->GetT0Means();
  82     }
  83
  84   for (Int_t i=0; i<4; i++)  mytask->SetT0Means(i,shiftcdb[i]);
  85 }