]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGDQ/dielectron/macrosLMEE/AddTask_Remi_LMEEPbPb2011AOD_Peri.C
Split: removed dirs now in AliPhysics
[u/mrichter/AliRoot.git] / PWGDQ / dielectron / macrosLMEE / AddTask_Remi_LMEEPbPb2011AOD_Peri.C
diff --git a/PWGDQ/dielectron/macrosLMEE/AddTask_Remi_LMEEPbPb2011AOD_Peri.C b/PWGDQ/dielectron/macrosLMEE/AddTask_Remi_LMEEPbPb2011AOD_Peri.C
deleted file mode 100644 (file)
index 1c22082..0000000
+++ /dev/null
@@ -1,133 +0,0 @@
-AliAnalysisTask *AddTask_Remi_LMEEPbPb2011AOD_Peri(Char_t* outputFileName="LMEEoutput.root", Bool_t runAll=kFALSE,Bool_t setMC=kFALSE,Bool_t getFromAlien=kFALSE, Bool_t PIDbaseline=kFALSE, Bool_t rejOnly=kTRUE) {
-  
-  Bool_t bESDANA=kFALSE; //Autodetect via InputHandler
-  //get the current analysis manager
-  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
-  if (!mgr) {
-    Error("AddTask_Remi_LMEEPbPb2011", "No analysis manager found.");
-    return 0;
-  }
-  
-  
-  //  create task and add it to the manager
-  //   gSystem->AddIncludePath("$ALICE_ROOT/PWGDQ/dielectron/macrosLMEE");
-  
-  
-  TString configBasePath("$TRAIN_ROOT/cbaumann_dielectron/");
-  TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
-  //  if (trainRoot.IsNull()) configBasePath= "/home/tanizaki/nfs/LMee_Deflection/ver1/";
-  if (trainRoot.IsNull()) configBasePath= "$ALICE_ROOT/PWGDQ/dielectron/macrosLMEE/";
-  
-  
-  if (getFromAlien &&
-      (!gSystem->Exec("alien_cp alien:///alice/cern.ch/user/c/cbaumann/PWGDQ/dielectron/macrosLMEE/ConfigRemiLMEEPbPb2011AOD.C .")) &&
-      (!gSystem->Exec("alien_cp alien:///alice/cern.ch/user/c/cbaumann/PWGDQ/dielectron/macrosLMEE/LMEECutLibRemi.C ."))
-      ) {
-    configBasePath=Form("%s/",gSystem->pwd());
-  }
-  
-  
-  TString configFile("ConfigRemiLMEEPbPb2011AOD.C");
-  TString configLMEECutLib("LMEECutLibRemi.C");
-  
-  TString configFilePath(configBasePath+configFile);
-  TString configLMEECutLibPath(configBasePath+configLMEECutLib);
-  
-  //AOD Usage currently tested with separate task, to be merged
-  if (mgr->GetInputEventHandler()->IsA()==AliAODInputHandler::Class()){
-    ::Info("AddTaskLMEEPbPb2011", "no dedicated AOD configuration");
-  }
-  else if (mgr->GetInputEventHandler()->IsA()==AliESDInputHandler::Class()){
-    ::Info("AddTaskLMEEPbPb2011AOD","switching on ESD specific code");
-    bESDANA=kTRUE;
-  }
-  
-  
-  //Do we have an MC handler?
-  Bool_t hasMC=setMC;
-  if (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0)
-    hasMC=kTRUE;
-  
-  
-  
-  if (!gROOT->GetListOfGlobalFunctions()->FindObject(configLMEECutLib.Data()))
-    gROOT->LoadMacro(configLMEECutLibPath.Data());
-  if (!gROOT->GetListOfGlobalFunctions()->FindObject(configFile.Data()))
-    gROOT->LoadMacro(configFilePath.Data());
-  
-  
-  LMEECutLibRemi* cutlib = new LMEECutLibRemi();
-  AliAnalysisTaskMultiDielectron *task=new AliAnalysisTaskMultiDielectron("MultiDiEData");
-  if (!hasMC) task->UsePhysicsSelection();
-  task->SetTriggerMask(AliVEvent::kMB+AliVEvent::kCentral+AliVEvent::kSemiCentral);
-  //  task->SelectCollisionCandidates(AliVEvent::kMB+AliVEvent::kCentral+AliVEvent::kSemiCentral);
-  //  task->SetRejectPileup();
-  task->SelectCollisionCandidates(AliVEvent::kAny);  
-  task->SetEventFilter(cutlib->GetEventCuts(LMEECutLibRemi::kPbPb2011pidITSTPCTOF)); //
-  
-  //load dielectron configuration file
-  
-  //add dielectron analysis with different cuts to the task
-  
-  
-  
-    AliDielectron *lowmass4=ConfigRemiLMEEPbPb2011AOD(4,hasMC,bESDANA);
-    task->AddDielectron(lowmass4);
-    printf("add: %s\n",lowmass4->GetName());
-    AliDielectron *lowmass8=ConfigRemiLMEEPbPb2011AOD(8,hasMC,bESDANA);
-    task->AddDielectron(lowmass8);
-    printf("add: %s\n",lowmass8->GetName());
-    AliDielectron *lowmass12=ConfigRemiLMEEPbPb2011AOD(12,hasMC,bESDANA);
-    task->AddDielectron(lowmass12);
-    printf("add: %s\n",lowmass12->GetName());
-    AliDielectron *lowmass16=ConfigRemiLMEEPbPb2011AOD(16,hasMC,bESDANA);
-    task->AddDielectron(lowmass16);
-    printf("add: %s\n",lowmass16->GetName());
-
-    AliDielectron *lowmass20=ConfigRemiLMEEPbPb2011AOD(20,hasMC,bESDANA);
-    task->AddDielectron(lowmass20);
-    printf("add: %s\n",lowmass20->GetName());
-    AliDielectron *lowmass24=ConfigRemiLMEEPbPb2011AOD(24,hasMC,bESDANA);
-    task->AddDielectron(lowmass24);
-    printf("add: %s\n",lowmass24->GetName());
-    AliDielectron *lowmass28=ConfigRemiLMEEPbPb2011AOD(28,hasMC,bESDANA);
-    task->AddDielectron(lowmass28);
-    printf("add: %s\n",lowmass28->GetName());
-
-  
-  mgr->AddTask(task);
-
-  //create output container
-  AliAnalysisDataContainer *coutput1 =
-       mgr->CreateContainer("Remi_LMEEPbPb2011Peri_tree",
-               TTree::Class(),
-               AliAnalysisManager::kExchangeContainer,
-               outputFileName);
-
-  AliAnalysisDataContainer *cOutputHist1 =
-       mgr->CreateContainer("Remi_LMEEPbPb2011Peri_out",
-               TList::Class(),
-               AliAnalysisManager::kOutputContainer,
-               outputFileName);
-
-  AliAnalysisDataContainer *cOutputHist2 =
-       mgr->CreateContainer("Remi_LMEEPbPb2011Peri_CF",
-               TList::Class(),
-               AliAnalysisManager::kOutputContainer,
-               outputFileName);
-
-  AliAnalysisDataContainer *cOutputHist3 =
-       mgr->CreateContainer("Remi_EventStatPbPb2011Peri",
-               TH1D::Class(),
-               AliAnalysisManager::kOutputContainer,
-               outputFileName);
-
-
-  mgr->ConnectInput(task,  0, mgr->GetCommonInputContainer());
-  mgr->ConnectOutput(task, 0, coutput1 );
-  mgr->ConnectOutput(task, 1, cOutputHist1);
-  mgr->ConnectOutput(task, 2, cOutputHist2);
-  mgr->ConnectOutput(task, 3, cOutputHist3);
-
-  return task;
-}