PWGHF/hfe/macros/configs/PbPb/ConfigHFE_FLOW_TPCEMCal_EP.C

   1 AliAnalysisTaskFlowTPCEMCalEP* ConfigHFE_FLOW_TPCEMCal_EP(Bool_t useMC){
   2   //
   3   // HFE standard task configuration
   4   //
   5
   6   Bool_t kAnalyseTaggedTracks = kTRUE;
   7
   8   AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");
   9   hfecuts->CreateStandardCuts();
  10   hfecuts->SetMinNClustersTPC(100);
  11   hfecuts->SetMinRatioTPCclusters(0.6);
  12   hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);
  13   hfecuts->SetMinNClustersITS(3);
  14   hfecuts->SetCutITSpixel(AliHFEextraCuts::kFirst);
  15   hfecuts->SetCheckITSLayerStatus(kFALSE);
  16   hfecuts->SetVertexRange(10.);
  17   hfecuts->SetTOFPIDStep(kFALSE);
  18   hfecuts->SetPtRange(2, 50);
  19   hfecuts->SetMaxImpactParam(1,2);
  20
  21   AliAnalysisTaskFlowTPCEMCalEP *task = new AliAnalysisTaskFlowTPCEMCalEP("HFE v2");
  22   printf("task ------------------------ %p\n ", task);
  23   task->SetHFECuts(hfecuts);
  24   task->SetInvariantMassCut(0.05);
  25
  26   // Define PID
  27   AliHFEpid *pid = task->GetPID();
  28   if(useMC) pid->SetHasMCData(kTRUE);
  29   pid->AddDetector("TPC", 0);
  30   pid->AddDetector("EMCAL", 1);
  31   // change E/p cuts
  32   AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);
  33   emcpid->SetEoPMax(1.2);
  34   emcpid->SetEoPMim(0.9);
  35
  36   Double_t params[4];
  37   char *cutmodel;
  38   if(useMC){
  39           // Monte-Carlo needs modelling of the falling mean with momentum at low momentum
  40           // for high momentum it is consistent with a flat -0.94
  41           cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";
  42           Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};
  43           for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];
  44   } else {
  45           // Data is consistent with a flat 0.12
  46           cutmodel = "pol0";
  47           //params[0] = -0.0015;
  48           //params[0] = -3.0;
  49           //params[0] = -0.05; //sigma min
  50           params[0] = -1.0; //sigma min
  51   }
  52   pid->ConfigureTPCdefaultCut(cutmodel, params,3.0);
  53
  54   printf("*************************************\n");
  55   printf("Configuring standard Task:\n");
  56 //  task->PrintStatus();
  57   pid->PrintStatus();
  58   printf("*************************************\n");
  59   return task;
  60 }