Merge branch 'feature-movesplit'

[u/mrichter/AliRoot.git] / PWGHF / hfe / macros / AddTaskHFECal.C

diff --git a/PWGHF/hfe/macros/AddTaskHFECal.C b/PWGHF/hfe/macros/AddTaskHFECal.C
index ce542c5d092ad374282a825e713ff9b133c13c37..31ddac5b989fdb64957ce04b50594553861d87d2 100644 (file)
--- a/PWGHF/hfe/macros/AddTaskHFECal.C
+++ b/PWGHF/hfe/macros/AddTaskHFECal.C
@@ -1,4 +1,4 @@
-AliAnalysisTask *AddTaskHFECal(Bool_t MassConst)
+AliAnalysisTask *AddTaskHFECal(Bool_t MassConst, Bool_t MassWidthCut, Bool_t MassNonlinear)
 {
   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   if (!mgr) {
@@ -23,18 +23,21 @@ AliAnalysisTask *AddTaskHFECal(Bool_t MassConst)
     MCthere=kTRUE;
   }
   
-  
+  Bool_t MassCal = kTRUE;
+  Double_t masscut = 0.05;
+  if(!MassConst)masscut = 0.1;
   //analysis task 
   gROOT->LoadMacro("$ALICE_ROOT/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C");
 
-  AliAnalysisTaskHFECal *hfetaskCent = ConfigHFECal(MCthere,MassConst);
-  AliAnalysisTaskHFECal *hfetaskTrig= ConfigHFECal(MCthere,MassConst);
-  AliAnalysisTaskHFECal *hfetaskTrig2= ConfigHFECal(MCthere,MassConst);
+  AliAnalysisTaskHFECal *hfetaskCent = ConfigHFECal(MCthere,MassConst,MassWidthCut,MassCal,MassNonlinear,0.5,0.1,masscut,-5,1);
+  AliAnalysisTaskHFECal *hfetaskTrig= ConfigHFECal(MCthere,MassConst,MassWidthCut,MassCal,MassNonlinear,0.5,0.1,masscut,-5,1);
+  AliAnalysisTaskHFECal *hfetaskTrig2= ConfigHFECal(MCthere,MassConst,MassWidthCut,MassCal,MassNonlinear,0.5,0.1,masscut,-5,1);
  
   mgr->AddTask(hfetaskCent);
   mgr->AddTask(hfetaskTrig);
   mgr->AddTask(hfetaskTrig2);
   
+  
   // central trigger
   hfetaskCent->SelectCollisionCandidates(AliVEvent::kCentral);
   
@@ -45,7 +48,7 @@ AliAnalysisTask *AddTaskHFECal(Bool_t MassConst)
   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("HFE_Results_EMCalCentral", TList::Class(),AliAnalysisManager::kOutputContainer, containerName.Data());
   mgr->ConnectInput(hfetaskCent, 0, cinput);
   mgr->ConnectOutput(hfetaskCent, 1, coutput1);
-  
+
   //L1 gamma trigger
   hfetaskTrig->SelectCollisionCandidates(AliVEvent::kEMCEGA);
   
@@ -57,10 +60,10 @@ AliAnalysisTask *AddTaskHFECal(Bool_t MassConst)
   mgr->ConnectInput(hfetaskTrig, 0, cinput);
   mgr->ConnectOutput(hfetaskTrig, 1, coutput1);
   
-
-  //Jet trigger
+  
   // semi-central
-  hfetaskTrig2->SelectCollisionCandidates(AliVEvent::kSemiCentral);
+  //hfetaskTrig2->SelectCollisionCandidates(AliVEvent::kSemiCentral | AliVEvent::kCentral);
+  hfetaskTrig2->SelectCollisionCandidates(AliVEvent::kAny);
   
   TString containerName3 = mgr->GetCommonFileName();
   containerName3 += ":PWGHF_hfeCalSemiCentral";
@@ -69,24 +72,23 @@ AliAnalysisTask *AddTaskHFECal(Bool_t MassConst)
   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("HFE_Results_SemiCentral", TList::Class(),AliAnalysisManager::kOutputContainer, containerName3.Data());
   mgr->ConnectInput(hfetaskTrig2, 0, cinput);
   mgr->ConnectOutput(hfetaskTrig2, 1, coutput1);
-  
 
+ 
   //if(MCthere)
     //{
     //MB trigger
-    AliAnalysisTaskHFECal *hfetaskMB = ConfigHFECal(MCthere,MassConst);
+    AliAnalysisTaskHFECal *hfetaskMB = ConfigHFECal(MCthere,MassConst,MassWidthCut,MassCal,MassNonlinear,0.5,0.1,masscut,-5,1);
     mgr->AddTask(hfetaskMB);
     hfetaskMB->SelectCollisionCandidates(AliVEvent::kMB);
 
-    TString containerName3 = mgr->GetCommonFileName();
-    containerName3 += ":PWGHF_hfeCalkMB";
+    TString containerName4 = mgr->GetCommonFileName();
+    containerName4 += ":PWGHF_hfeCalkMB";
 
      AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
-     AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("HFE_Results_EMCalMB", TList::Class(),AliAnalysisManager::kOutputContainer, containerName3.Data());
+     AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("HFE_Results_EMCalMB", TList::Class(),AliAnalysisManager::kOutputContainer, containerName4.Data());
      mgr->ConnectInput(hfetaskMB, 0, cinput);
      mgr->ConnectOutput(hfetaskMB, 1, coutput1);
      //}
-  
 
   return NULL;
 }