new macro for EMCal electron
authorssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 19 May 2012 09:02:52 +0000 (09:02 +0000)
committerssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 19 May 2012 09:02:52 +0000 (09:02 +0000)
PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C [new file with mode: 0644]

diff --git a/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C b/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C
new file mode 100644 (file)
index 0000000..de876cb
--- /dev/null
@@ -0,0 +1,60 @@
+AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC){\r
+  //\r
+  // HFE standard task configuration\r
+  //\r
+\r
+  Bool_t kAnalyseTaggedTracks = kTRUE;\r
+  \r
+  AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");\r
+  hfecuts->CreateStandardCuts();\r
+  hfecuts->SetMinNClustersTPC(100);\r
+  hfecuts->SetMinRatioTPCclusters(0.6);\r
+  hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);\r
+  hfecuts->SetMinNClustersITS(3);\r
+  hfecuts->SetCutITSpixel(AliHFEextraCuts::kFirst);\r
+  hfecuts->SetCheckITSLayerStatus(kFALSE);\r
+  hfecuts->SetVertexRange(10.);\r
+  hfecuts->SetTOFPIDStep(kFALSE);\r
+  hfecuts->SetPtRange(2, 50);\r
+  hfecuts->SetMaxImpactParam(1,2);\r
+  \r
+  AliAnalysisTaskHFECal *task = new AliAnalysisTaskHFECal("HFE v2");\r
+  printf("task ------------------------ %p\n ", task);\r
+  task->SetHFECuts(hfecuts);\r
+  task->SetInvariantMassCut(0.05);\r
+\r
+  // Define PID\r
+  AliHFEpid *pid = task->GetPID();\r
+  if(useMC) pid->SetHasMCData(kTRUE);\r
+  pid->AddDetector("TPC", 0);\r
+  pid->AddDetector("EMCAL", 1);\r
+  // change E/p cuts\r
+  AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);\r
+  emcpid->SetEoPMax(1.2);\r
+  emcpid->SetEoPMim(0.9);\r
+\r
+  Double_t params[4];\r
+  char *cutmodel;\r
+  if(useMC){\r
+         // Monte-Carlo needs modelling of the falling mean with momentum at low momentum\r
+         // for high momentum it is consistent with a flat -0.94\r
+         cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";\r
+         Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};\r
+         for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];\r
+  } else {\r
+         // Data is consistent with a flat 0.12\r
+         cutmodel = "pol0";\r
+         //params[0] = -0.0015;\r
+         //params[0] = -3.0;\r
+         //params[0] = -0.05; //sigma min\r
+         params[0] = -1.0; //sigma min\r
+  }\r
+  pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); \r
+\r
+  printf("*************************************\n");\r
+  printf("Configuring standard Task:\n");\r
+//  task->PrintStatus();\r
+  pid->PrintStatus();\r
+  printf("*************************************\n"); \r
+  return task;\r
+}\r