PWG1/CalibMacros/Pass0/AddTaskTRDCalib.C

   1 //=============================================================================
   2 //
   3 // *** AddTaskTRDCalib
   4 //
   5 // This macros setup the TRD calibration task
   6 //
   7 //=============================================================================
   8
   9 AliAnalysisTask  *AddTaskTRDCalib(Int_t runNumber)
  10 {
  11   gSystem->Load("libTRDcalib");
  12   // pointer to the analysis manager
  13   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  14   if (!mgr) {
  15     Error("AddTaskTRDCalib", "No analysis manager to connect to.");
  16     return NULL;
  17   }
  18
  19   // check the input handler
  20   if (!mgr->GetInputEventHandler()) {
  21     ::Error("AddTask", "This task requires an input event handler");
  22     return NULL;
  23   }
  24
  25   ////////////////////////////////////////////
  26   // Number of timebins
  27   ///////////////////////////////////////////
  28   AliTRDcalibDB *calib = AliTRDcalibDB::Instance();
  29   Int_t nbOfTimeBins = calib->GetNumberOfTimeBinsDCS();
  30   if(nbOfTimeBins < 0) nbOfTimeBins = 27;
  31   ////////////////////////////////////////////
  32   //
  33   /////////////////////////////////////////////
  34   Int_t versiongain, subversiongain, versionvdrift, subversionvdrift;
  35
  36   /////////////////////////
  37   // The TRD calib Task
  38   /////////////////////////
  39   AliTRDCalibTask *calibTask = new AliTRDCalibTask();
  40   calibTask->SetHisto2d(kTRUE);
  41   calibTask->SetVector2d(kFALSE);
  42   calibTask->SetVdriftLinear(kTRUE);
  43   calibTask->SetExbAlt(kTRUE);
  44   calibTask->SetNz(0,0);
  45   calibTask->SetNrphi(0,0);
  46   calibTask->SetNz(0,1);
  47   calibTask->SetNrphi(0,1);
  48   calibTask->SetNz(0,2);
  49   calibTask->SetNrphi(0,2);
  50   calibTask->SetLow(0);
  51   calibTask->SetHigh(30);
  52   calibTask->SetFillZero(kFALSE);
  53   // now
  54   calibTask->AddSelectedTriggerClass("C0OB0-ABCE-NOPF-ALL");
  55   calibTask->AddSelectedTriggerClass("CTRDCO2-ABCE-NOPF-CENT");
  56   calibTask->AddSelectedTriggerClass("CTRDCO2-ABCE-NOPF-TRD");
  57   calibTask->AddSelectedTriggerClass("CTRDCO2-ABCE-NOPF-ALL");
  58   calibTask->SetReject(kTRUE);
  59   // before
  60   //calibTask->AddSelectedTriggerClass("CINT1B-ABCE-NOPF-ALL");
  61   //calibTask->AddSelectedTriggerClass("CINT1WU-B-NOPF-ALL");
  62   //calibTask->AddSelectedTriggerClass("CINT7WU-B-NOPF-ALL");
  63   //calibTask->AddSelectedTriggerClass("CINT7WU-I-NOPF-ALL");
  64   //calibTask->SetReject(kFALSE);
  65   //calibTask->SetDebug(2);
  66   calibTask->SetNbTimeBins(nbOfTimeBins);
  67   //calibTask->SetMaxEvent(10);
  68   //calibTask->SetThresholdP(1.0);
  69   calibTask->SetRequirePrimaryVertex(kTRUE);
  70   calibTask->SetMinNbOfContributors(1);
  71   calibTask->SetMaxCluster(100.0);
  72   calibTask->SetNbMaxCluster(2);
  73   //calibTask->SetLimitChargeIntegration(kTRUE);
  74
  75
  76   /////////////////////////////
  77   // Track cuts
  78   /////////////////////////////
  79   AliESDtrackCuts *trackCuts = new AliESDtrackCuts("trackcuts","trackcuts");
  80   trackCuts->SetMinNClustersTPC(50);
  81   trackCuts->SetMaxChi2PerClusterTPC(3.5);
  82   //trackCuts->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
  83   trackCuts->SetRequireTPCRefit(kTRUE);
  84   //trackCuts->SetRequireITSRefit(kTRUE);
  85   //trackCuts->SetMinNsigmaToVertex(10);
  86   trackCuts->SetRequireSigmaToVertex(kFALSE);
  87   trackCuts->SetAcceptKinkDaughters(kFALSE);
  88   trackCuts->SetMaxDCAToVertexZ(30.0);
  89   trackCuts->SetMaxDCAToVertexXY(3.0);
  90   trackCuts->SetDCAToVertex2D(kFALSE);
  91
  92   calibTask->SetESDtrackCuts(trackCuts);
  93
  94   mgr->AddTask(calibTask);
  95
  96   AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  97   //AliAnalysisDataContainer *cinput = mgr->GetCommonOutputContainer();
  98
  99   if (!cinput) cinput = mgr->CreateContainer("cchain",TChain::Class(),
 100                                       AliAnalysisManager::kInputContainer);
 101
 102   AliAnalysisDataContainer *coutput =mgr->CreateContainer("TRDCalib",TList::Class(), AliAnalysisManager::kOutputContainer, "AliESDfriends_v1.root");
 103
 104
 105   mgr->ConnectInput(calibTask,0,cinput);
 106   mgr->ConnectOutput(calibTask,1,coutput);
 107   return calibTask;
 108
 109 }