]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGHF/hfe/macros/configs/PbPb/ConfigHFEElecHadronCorl.C
end-of-line normalization
[u/mrichter/AliRoot.git] / PWGHF / hfe / macros / configs / PbPb / ConfigHFEElecHadronCorl.C
index d459c79b4c17df1215365c29fad0cef83b8312f8..dc38d3a9cf9c7eac33bd5605eccca4f5cb2e2c6e 100644 (file)
@@ -1,68 +1,68 @@
-AliAnalysisTaskElecHadronCorrel* ConfigHFEElecHadronCorl(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->SetTPCiter1(kTRUE);\r
-  hfecuts->SetMinNClustersTPC(120);\r
-  hfecuts->SetMinRatioTPCclusters(0.6);\r
-  //hfecuts->SetTPCmodes(AliHFEextraCuts::kCrossedRows, AliHFEextraCuts::kFoundOverCR);\r
-  hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);\r
-  //hfecuts->SetUseMixedVertex(kTRUE);\r
-\r
-  hfecuts->SetMinNClustersITS(3);\r
-  hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny); //shud be put back\r
-  hfecuts->SetCheckITSLayerStatus(kFALSE); // shud be put back\r
-//  hfecuts->UnsetVertexRequirement();\r
-  hfecuts->SetVertexRange(10.);\r
-  //hfecuts->SetSigmaToVertex(10);\r
-  hfecuts->SetTOFPIDStep(kFALSE);\r
-//  hfecuts->SetQAOn();\r
- hfecuts->SetPtRange(0, 30);\r
- hfecuts->SetMaxImpactParam(1,2);\r
-\r
-  AliAnalysisTaskElecHadronCorrel *task = new AliAnalysisTaskElecHadronCorrel("HFE-hadron correlations");\r
-  printf("task ------------------------ %p\n ", task);\r
-  task->SetHFECuts(hfecuts);\r
-  task->SetInvariantMassCut(0.05);\r
-//  task->SetRemovePileUp(kTRUE);\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.8);\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
+AliAnalysisTaskElecHadronCorrel* ConfigHFEElecHadronCorl(Bool_t useMC){
+  //
+  // HFE standard task configuration
+  //
+
+  Bool_t kAnalyseTaggedTracks = kTRUE;
+  
+  AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");
+  hfecuts->CreateStandardCuts();
+  //hfecuts->SetTPCiter1(kTRUE);
+  hfecuts->SetMinNClustersTPC(120);
+  hfecuts->SetMinRatioTPCclusters(0.6);
+  //hfecuts->SetTPCmodes(AliHFEextraCuts::kCrossedRows, AliHFEextraCuts::kFoundOverCR);
+  hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);
+  //hfecuts->SetUseMixedVertex(kTRUE);
+
+  hfecuts->SetMinNClustersITS(3);
+  hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny); //shud be put back
+  hfecuts->SetCheckITSLayerStatus(kFALSE); // shud be put back
+//  hfecuts->UnsetVertexRequirement();
+  hfecuts->SetVertexRange(10.);
+  //hfecuts->SetSigmaToVertex(10);
+  hfecuts->SetTOFPIDStep(kFALSE);
+//  hfecuts->SetQAOn();
+ hfecuts->SetPtRange(0, 30);
+ hfecuts->SetMaxImpactParam(1,2);
+
+  AliAnalysisTaskElecHadronCorrel *task = new AliAnalysisTaskElecHadronCorrel("HFE-hadron correlations");
+  printf("task ------------------------ %p\n ", task);
+  task->SetHFECuts(hfecuts);
+  task->SetInvariantMassCut(0.05);
+//  task->SetRemovePileUp(kTRUE);
+
+  // Define PID
+  AliHFEpid *pid = task->GetPID();
+  if(useMC) pid->SetHasMCData(kTRUE);
+  pid->AddDetector("TPC", 0);
+  pid->AddDetector("EMCAL", 1);
+  // change E/p cuts
+  AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);
+  emcpid->SetEoPMax(1.2);
+  emcpid->SetEoPMim(0.8);
+
+  Double_t params[4];
+  char *cutmodel;
+  if(useMC){
+         // Monte-Carlo needs modelling of the falling mean with momentum at low momentum
+         // for high momentum it is consistent with a flat -0.94
+         cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";
+         Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};
+         for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];
+  } else {
+         // Data is consistent with a flat 0.12
+         cutmodel = "pol0";
+         //params[0] = -0.0015;
+         //params[0] = -3.0;
+         //params[0] = -0.05; //sigma min
+         params[0] = -1.0; //sigma min
+  }
+  pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); 
+
+  printf("*************************************\n");
+  printf("Configuring standard Task:\n");
+//  task->PrintStatus();
+  pid->PrintStatus();
+  printf("*************************************\n"); 
+  return task;
+}