[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / macros / runEMCalJetAnalysis.C

// runEMCalJetAnalysis.C
// =====================
// This macro can be used to run a jet analysis within the EMCal Jet Framework.
//
// Examples:
// -> Analyze ESDs from the pA pilot run on the AliEn grid with your task in AnaClass.cxx/.h
//     dataType = "ESD", useGrid = kTRUE, pattern = "*ESDs/pass2/*ESDs.root", addCXXs = "AnaClass.cxx", 
//     addHs = "AnaClass.h", gridDir = "/alice/data/2012/LHC12g", gridMode = "full", runNumbers = "188359 188362"
//     
// -> Analyze AODs (up to 96 files) locally given in files_aod.txt
//     dataType = "AOD", useGrid = kFALSE, numLocalFiles = 96
//
// MERGING ON ALIEN
// ++++++++++++++++
// If you run on the grid, you can monitor the jobs with alimonitor.cern.ch. When enough of them are in DONE state,
// you have to merge the output. This can be done automatically, if you just change the gridMode to "terminate" and
// give the EXACT name of the task whose output should be merged in uniqueName.
// 
//
// Authors: R. Haake, S. Aiola

#include <ctime>
#include "TGrid.h"

AliAnalysisGrid* CreateAlienHandler(const char* uniqueName, const char* gridDir, const char* gridMode, const char* runNumbers, 
                                     const char* pattern, TString additionalCode, TString additionalHeaders, Int_t maxFilesPerWorker, 
                                     Int_t workerTTL, Bool_t isMC);
                                    
//______________________________________________________________________________
void runEMCalJetAnalysis(
         const char*    dataType            = "ESD",                       // set the analysis type, AOD, ESD or sESD
         Bool_t         useGrid             = kFALSE,                      // local or grid
         const char*    gridMode            = "test",                      // set the grid run mode (can be "full", "test", "offline", "submit" or "terminate")
         const char*    pattern             = "*ESDs/pass2/*ESDs.root",    // file pattern (here one can specify subdirs like passX etc.) (used on grid)
         const char*    gridDir             = "/alice/data/2012/LHC12g",   // dir on alien, where the files live (used on grid)
         const char*    runNumbers          = "188359 188362",             // considered run numbers (used on grid)
         UInt_t         numLocalFiles       = 50,                          // number of files analyzed locally  
         const char*    runPeriod           = "LHC12g",                    // set the run period (used on grid)
         const char*    uniqueName          = "EMCalJF_LEGOTrainTest",     // sets base string for the name of the task on the grid
         UInt_t         pSel                = AliVEvent::kAny,             // used event selection for every task except for the analysis tasks
         Bool_t         useTender           = kFALSE,                      // trigger, if tender task should be used
         Bool_t         isMC                = kFALSE,                      // trigger, if MC handler should be used

         // Here you have to specify additional code files you want to use but that are not in aliroot
         const char*    addCXXs             = "",
         const char*    addHs               = "",

         // These two settings depend on the dataset and your quotas on the AliEN services
         Int_t          maxFilesPerWorker   = 4,
         Int_t          workerTTL           = 7200

         )
{

  // Some pre-settings and constants

  enum AlgoType {kKT, kANTIKT};
  enum JetType  {kFULLJETS, kCHARGEDJETS, kNEUTRALJETS};
  gSystem->SetFPEMask();
  gSystem->Setenv("ETRAIN_ROOT", ".");
  gSystem->Setenv("ETRAIN_PERIOD", runPeriod);
  // change this objects to strings
  TString usedData(dataType);
  TString additionalCXXs(addCXXs);
  TString additionalHs(addHs);
  cout << dataType << " analysis chosen" << endl;
  if (useGrid)  
  {
    cout << "-- using AliEn grid.\n";
    if (usedData == "sESD") 
    {
      cout << "Skimmed ESD analysis not available on the grid!" << endl;
      return;
    }
  }
  else
    cout << "-- using local analysis.\n";
  

  // Load necessary libraries
  LoadLibs();

  // Create analysis manager
  AliAnalysisManager* mgr = new AliAnalysisManager(uniqueName);

  // Check type of input and create handler for it
  TString localFiles("-1");
  if(usedData == "AOD")
  {
    localFiles = "files_aod.txt";
    gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/train/AddAODHandler.C");
    AliAODInputHandler* aodH = AddAODHandler();
  }
  else if((usedData == "ESD") || (usedData == "sESD"))
  {
    if (usedData == "ESD")
      localFiles = "files_esd.txt";
    else
      localFiles = "files_sesd.txt";
    
    gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/train/AddESDHandler.C");
    AliESDInputHandler* esdH = AddESDHandler();
  }
  else
  {
    cout << "Data type not recognized! You have to specify ESD, AOD, or sESD!\n";
  }
  
  if(!useGrid)
    cout << "Using " << localFiles.Data() << " as input file list.\n";

  // Create MC handler, if MC is demanded
  if (isMC && (usedData != "AOD"))
  {
    AliMCEventHandler* mcH = new AliMCEventHandler();
    mcH->SetPreReadMode(AliMCEventHandler::kLmPreRead);
    mcH->SetReadTR(kTRUE);
    mgr->SetMCtruthEventHandler(mcH); 
  }
  
  // ################# Now: Add some basic tasks

  // Physics selection task
  gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEmcalPhysicsSelection.C");
  AliPhysicsSelectionTask *physSelTask = AddTaskEmcalPhysicsSelection(kTRUE, kTRUE, pSel, 5, 5, 10, kTRUE, -1, -1, -1, -1);

  if (!physSelTask) 
  {
    cout << "no physSelTask"; 
    return; 
  }

  // Centrality task
  if (usedData == "ESD") {
    gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");
    AliCentralitySelectionTask *centralityTask = AddTaskCentrality(kTRUE);
  }

  // Compatibility task, only needed for skimmed ESD
  if (usedData == "sESD") {
    gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEmcalCompat.C");
    AliEmcalCompatTask *comptask = AddTaskEmcalCompat();
  }

  // Setup task
  gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEmcalSetup.C");
  AliEmcalSetupTask *setupTask = AddTaskEmcalSetup();
  setupTask->SetGeoPath("$ALICE_ROOT/OADB/EMCAL");
  
  UInt_t tpc = AliAnalysisTaskEmcal::kTPC;

  // Tender Supplies
  if (useTender)
  {
    gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEMCALTender.C");
    AliAnalysisTaskSE *tender = AddTaskEMCALTender(runPeriod, kFALSE, kTRUE, kTRUE, kTRUE, kTRUE, kTRUE, kFALSE, kTRUE, kTRUE, kTRUE,
                                                   kFALSE, AliEMCALRecParam::kClusterizerNxN);
    if (usedData != "AOD" && !useGrid) {
      AliTender *alitender = dynamic_cast<AliTender*>(tender);
      alitender->SetDefaultCDBStorage("local://$ALICE_ROOT/OCDB"); 
    }
  }

  // ################# Now: Call jet preparation macro (picotracks, hadronic corrected caloclusters, ...) 
  
  // Jet preparation
  gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskJetPreparation.C");
  AddTaskJetPreparation(dataType);

  // ################# Now: Add jet finders+analyzers

  gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C");
  AliEmcalJetTask* jetFinderTask = AddTaskEmcalJet("PicoTracks", "CaloClustersCorr", kANTIKT, 0.2, kCHARGEDJETS, 0.150, 0.300);

  // Here you can put in your AddTaskMacro for your task
  gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetSample.C");
  AliAnalysisTaskEmcalJetSample* anaTask = AddTaskEmcalJetSample("PicoTracks", "CaloClustersCorr", jetFinderTask->GetName(), "");


  // Set the physics selection for all given tasks
  TObjArray *toptasks = mgr->GetTasks();
  for (Int_t i=0; i<toptasks->GetEntries(); ++i) 
  {
    AliAnalysisTaskSE *task = dynamic_cast<AliAnalysisTaskSE*>(toptasks->At(i));
    if (!task)
      continue;
    if (task->InheritsFrom("AliPhysicsSelectionTask"))
      continue;
    ::Info("setPSel", "Set physics selection for %s (%s)", task->GetName(), task->ClassName());
    task->SelectCollisionCandidates(pSel);
  }

  mgr->SetUseProgressBar(1, 25);
        
  if (!mgr->InitAnalysis()) 
    return;
  mgr->PrintStatus();

  if (useGrid) 
  {  // GRID CALCULATION

    AliAnalysisGrid *plugin = CreateAlienHandler(uniqueName, gridDir, gridMode, runNumbers, pattern, additionalCXXs, additionalHs, maxFilesPerWorker, workerTTL, isMC);
    mgr->SetGridHandler(plugin);

    // start analysis
    cout << "Starting GRID Analysis...";
    mgr->SetDebugLevel(0);
    mgr->StartAnalysis("grid");
  }
  else
  {  // LOCAL CALCULATION

    TChain* chain = 0;
    if (usedData == "AOD") 
    {
      gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/CreateAODChain.C");
      chain = CreateAODChain(localFiles.Data(), numLocalFiles);
    }
    else
    {  // ESD or skimmed ESD
      gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/CreateESDChain.C");
      chain = CreateESDChain(localFiles.Data(), numLocalFiles);
    }
    
    // start analysis
    cout << "Starting LOCAL Analysis...";
    mgr->SetDebugLevel(2);
    mgr->StartAnalysis("local", chain);
  }
}

//______________________________________________________________________________
void LoadLibs()
{
  // Load common libraries (better too many than too few)
  gSystem->Load("libTree");
  gSystem->Load("libVMC");
  gSystem->Load("libGeom");
  gSystem->Load("libGui");
  gSystem->Load("libXMLParser");
  gSystem->Load("libMinuit");
  gSystem->Load("libMinuit2");
  gSystem->Load("libProof");
  gSystem->Load("libPhysics");
  gSystem->Load("libSTEERBase");
  gSystem->Load("libESD");
  gSystem->Load("libAOD");
  gSystem->Load("libOADB");
  gSystem->Load("libANALYSIS");
  gSystem->Load("libCDB");
  gSystem->Load("libRAWDatabase");
  gSystem->Load("libSTEER");
  gSystem->Load("libEVGEN");
  gSystem->Load("libANALYSISalice");
  gSystem->Load("libCORRFW");
  gSystem->Load("libTOFbase");
  //gSystem->Load("libTOFrec");
  gSystem->Load("libRAWDatabase.so");
  gSystem->Load("libRAWDatarec.so");
  gSystem->Load("libTPCbase.so");
  gSystem->Load("libTPCrec.so");
  gSystem->Load("libITSbase.so");
  gSystem->Load("libITSrec.so");
  gSystem->Load("libTRDbase.so");
  gSystem->Load("libTENDER.so");
  gSystem->Load("libSTAT.so");
  gSystem->Load("libTRDrec.so");
  gSystem->Load("libHMPIDbase.so");
  gSystem->Load("libPWGPP.so");
  gSystem->Load("libPWGHFbase");
  gSystem->Load("libPWGDQdielectron");
  gSystem->Load("libPWGHFhfe");
  gSystem->Load("libEMCALUtils");
  gSystem->Load("libPHOSUtils");
  gSystem->Load("libPWGCaloTrackCorrBase");
  gSystem->Load("libEMCALraw");
  gSystem->Load("libEMCALbase");
  gSystem->Load("libEMCALrec");
  gSystem->Load("libTRDbase");
  gSystem->Load("libVZERObase");
  gSystem->Load("libVZEROrec");
  gSystem->Load("libTENDER");
  gSystem->Load("libTENDERSupplies");
  gSystem->Load("libPWGEMCAL");
  gSystem->Load("libPWGGAEMCALTasks");
  gSystem->Load("libPWGTools");
  gSystem->Load("libPWGCFCorrelationsBase");
  gSystem->Load("libPWGCFCorrelationsDPhi");

  //load CGAL, Fastjet and SISCone
  gSystem->Load("libCGAL");
  gSystem->Load("libfastjet");
  gSystem->Load("libSISConePlugin");

  gSystem->Load("libJETAN");
  gSystem->Load("libFASTJETAN");
  gSystem->Load("libPWGJEEMCALJetTasks");


  // include paths
  gSystem->AddIncludePath("-Wno-deprecated");
  gSystem->AddIncludePath("-I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/EMCAL");
  gSystem->AddIncludePath("-I$ALICE_ROOT/PWGDQ/dielectron -I$ALICE_ROOT/PWGHF/hfe");
  gSystem->AddIncludePath("-I$ALICE_ROOT/JETAN -I$ALICE_ROOT/JETAN/fastjet");
}

AliAnalysisGrid* CreateAlienHandler(const char* uniqueName, const char* gridDir, const char* gridMode, const char* runNumbers, 
                                     const char* pattern, TString additionalCode, TString additionalHeaders, Int_t maxFilesPerWorker, 
                                     Int_t workerTTL, Bool_t isMC)
{
  TDatime currentTime;
  TString tmpName(uniqueName);

  // Only add current date and time when not in terminate mode! In this case the exact name has to be supplied by the user
  if(strcmp(gridMode, "terminate"))
  {
    tmpName += "_";
    tmpName += currentTime.GetDate();
    tmpName += "_";
    tmpName += currentTime.GetTime();
  }

  TString tmpAdditionalLibs("");
  tmpAdditionalLibs = Form("libTree.so libVMC.so libGeom.so libGui.so libXMLParser.so libMinuit.so libMinuit2.so libProof.so libPhysics.so libSTEERBase.so libESD.so libAOD.so libOADB.so libANALYSIS.so libCDB.so libRAWDatabase.so libSTEER.so libANALYSISalice.so libCORRFW.so libTOFbase.so libRAWDatabase.so libRAWDatarec.so libTPCbase.so libTPCrec.so libITSbase.so libITSrec.so libTRDbase.so libTENDER.so libSTAT.so libTRDrec.so libHMPIDbase.so libPWGPP.so libPWGHFbase.so libPWGDQdielectron.so libPWGHFhfe.so libEMCALUtils.so libPHOSUtils.so libPWGCaloTrackCorrBase.so libEMCALraw.so libEMCALbase.so libEMCALrec.so libTRDbase.so libVZERObase.so libVZEROrec.so libTENDER.so libTENDERSupplies.so libPWGEMCAL.so libPWGGAEMCALTasks.so libPWGTools.so libPWGCFCorrelationsBase.so libPWGCFCorrelationsDPhi.so  libCGAL.so libJETAN.so libfastjet.so libSISConePlugin.so libFASTJETAN.so libPWGJE.so libPWGmuon.so libPWGJEEMCALJetTasks.so %s %s",additionalCode.Data(),additionalHeaders.Data());


  TString macroName("");
  TString execName("");
  TString jdlName("");
  macroName = Form("%s.C", tmpName.Data());
  execName = Form("%s.sh", tmpName.Data());
  jdlName = Form("%s.jdl", tmpName.Data());

  AliAnalysisAlien *plugin = new AliAnalysisAlien();
  plugin->SetOverwriteMode();
  plugin->SetRunMode(gridMode);
     
  // Here you can set the (Ali)ROOT version you want to use
  plugin->SetAPIVersion("V1.1x");
  plugin->SetROOTVersion("v5-34-02-1");
  plugin->SetAliROOTVersion("v5-04-21-AN-1");

  plugin->SetGridDataDir(gridDir); // e.g. "/alice/sim/LHC10a6"
  plugin->SetDataPattern(pattern); //dir structure in run directory
  if (!isMC)
   plugin->SetRunPrefix("000");

  plugin->AddRunList(runNumbers);

  plugin->SetGridWorkingDir(Form("work/%s",tmpName.Data()));
  plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output

  plugin->SetAnalysisSource(additionalCode.Data());
  plugin->SetAdditionalLibs(tmpAdditionalLibs.Data());
  plugin->AddExternalPackage("boost::v1_43_0");
  plugin->AddExternalPackage("cgal::v3.6");
  plugin->AddExternalPackage("fastjet::v2.4.2");

  plugin->SetDefaultOutputs(kTRUE);
  //plugin->SetMergeExcludes("");
  plugin->SetAnalysisMacro(macroName.Data());
  plugin->SetSplitMaxInputFileNumber(maxFilesPerWorker);
  plugin->SetExecutable(execName.Data());
  plugin->SetTTL(workerTTL);
  plugin->SetInputFormat("xml-single");
  plugin->SetJDLName(jdlName.Data());
  plugin->SetPrice(1);      
  plugin->SetSplitMode("se");

  return plugin;
}
Commit	Line	Data
8628b70c	1	// runEMCalJetAnalysis.C
	2	// =====================
	3	// This macro can be used to run a jet analysis within the EMCal Jet Framework.
	4	//
	5	// Examples:
	6	// -> Analyze ESDs from the pA pilot run on the AliEn grid with your task in AnaClass.cxx/.h
	7	// dataType = "ESD", useGrid = kTRUE, pattern = "ESDs/pass2/ESDs.root", addCXXs = "AnaClass.cxx",
	8	// addHs = "AnaClass.h", gridDir = "/alice/data/2012/LHC12g", gridMode = "full", runNumbers = "188359 188362"
	9	//
	10	// -> Analyze AODs (up to 96 files) locally given in files_aod.txt
	11	// dataType = "AOD", useGrid = kFALSE, numLocalFiles = 96
	12	//
	13	// MERGING ON ALIEN
	14	// ++++++++++++++++
	15	// If you run on the grid, you can monitor the jobs with alimonitor.cern.ch. When enough of them are in DONE state,
	16	// you have to merge the output. This can be done automatically, if you just change the gridMode to "terminate" and
	17	// give the EXACT name of the task whose output should be merged in uniqueName.
	18	//
	19	//
	20	// Authors: R. Haake, S. Aiola
	21
	22	#include <ctime>
	23	#include "TGrid.h"
	24
	25	AliAnalysisGrid* CreateAlienHandler(const char* uniqueName, const char* gridDir, const char* gridMode, const char* runNumbers,
	26	const char* pattern, TString additionalCode, TString additionalHeaders, Int_t maxFilesPerWorker,
	27	Int_t workerTTL, Bool_t isMC);
	28
	29	//______________________________________________________________________________
	30	void runEMCalJetAnalysis(
	31	const char* dataType = "ESD", // set the analysis type, AOD, ESD or sESD
	32	Bool_t useGrid = kFALSE, // local or grid
	33	const char* gridMode = "test", // set the grid run mode (can be "full", "test", "offline", "submit" or "terminate")
	34	const char* pattern = "ESDs/pass2/ESDs.root", // file pattern (here one can specify subdirs like passX etc.) (used on grid)
	35	const char* gridDir = "/alice/data/2012/LHC12g", // dir on alien, where the files live (used on grid)
	36	const char* runNumbers = "188359 188362", // considered run numbers (used on grid)
	37	UInt_t numLocalFiles = 50, // number of files analyzed locally
	38	const char* runPeriod = "LHC12g", // set the run period (used on grid)
	39	const char* uniqueName = "EMCalJF_LEGOTrainTest", // sets base string for the name of the task on the grid
	40	UInt_t pSel = AliVEvent::kAny, // used event selection for every task except for the analysis tasks
	41	Bool_t useTender = kFALSE, // trigger, if tender task should be used
95449deb	42	Bool_t isMC = kFALSE, // trigger, if MC handler should be used
8628b70c	43
	44	// Here you have to specify additional code files you want to use but that are not in aliroot
	45	const char* addCXXs = "",
	46	const char* addHs = "",
	47
	48	// These two settings depend on the dataset and your quotas on the AliEN services
	49	Int_t maxFilesPerWorker = 4,
	50	Int_t workerTTL = 7200
	51
	52	)
	53	{
	54
	55	// Some pre-settings and constants
	56
	57	enum AlgoType {kKT, kANTIKT};
	58	enum JetType {kFULLJETS, kCHARGEDJETS, kNEUTRALJETS};
	59	gSystem->SetFPEMask();
	60	gSystem->Setenv("ETRAIN_ROOT", ".");
	61	gSystem->Setenv("ETRAIN_PERIOD", runPeriod);
	62	// change this objects to strings
	63	TString usedData(dataType);
	64	TString additionalCXXs(addCXXs);
	65	TString additionalHs(addHs);
	66	cout << dataType << " analysis chosen" << endl;
	67	if (useGrid)
	68	{
	69	cout << "-- using AliEn grid.\n";
	70	if (usedData == "sESD")
	71	{
	72	cout << "Skimmed ESD analysis not available on the grid!" << endl;
	73	return;
	74	}
	75	}
	76	else
	77	cout << "-- using local analysis.\n";
	78
	79
	80	// Load necessary libraries
	81	LoadLibs();
	82
	83	// Create analysis manager
	84	AliAnalysisManager* mgr = new AliAnalysisManager(uniqueName);
	85
	86	// Check type of input and create handler for it
	87	TString localFiles("-1");
	88	if(usedData == "AOD")
	89	{
	90	localFiles = "files_aod.txt";
	91	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/train/AddAODHandler.C");
	92	AliAODInputHandler* aodH = AddAODHandler();
	93	}
	94	else if((usedData == "ESD") \|\| (usedData == "sESD"))
	95	{
	96	if (usedData == "ESD")
	97	localFiles = "files_esd.txt";
	98	else
	99	localFiles = "files_sesd.txt";
	100
	101	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/train/AddESDHandler.C");
	102	AliESDInputHandler* esdH = AddESDHandler();
	103	}
	104	else
	105	{
	106	cout << "Data type not recognized! You have to specify ESD, AOD, or sESD!\n";
107	}
108
109	if(!useGrid)
110	cout << "Using " << localFiles.Data() << " as input file list.\n";
111
112	// Create MC handler, if MC is demanded
12eb4ac1	113	if (isMC && (usedData != "AOD"))
8628b70c	114	{
	115	AliMCEventHandler* mcH = new AliMCEventHandler();
	116	mcH->SetPreReadMode(AliMCEventHandler::kLmPreRead);
	117	mcH->SetReadTR(kTRUE);
	118	mgr->SetMCtruthEventHandler(mcH);
	119	}
	120
	121	// ################# Now: Add some basic tasks
	122
	123	// Physics selection task
	124	gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEmcalPhysicsSelection.C");
	125	AliPhysicsSelectionTask *physSelTask = AddTaskEmcalPhysicsSelection(kTRUE, kTRUE, pSel, 5, 5, 10, kTRUE, -1, -1, -1, -1);
	126
	127	if (!physSelTask)
	128	{
	129	cout << "no physSelTask";
	130	return;
	131	}
	132
	133	// Centrality task
	134	if (usedData == "ESD") {
	135	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");
	136	AliCentralitySelectionTask *centralityTask = AddTaskCentrality(kTRUE);
	137	}
	138
	139	// Compatibility task, only needed for skimmed ESD
	140	if (usedData == "sESD") {
	141	gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEmcalCompat.C");
	142	AliEmcalCompatTask *comptask = AddTaskEmcalCompat();
	143	}
	144
	145	// Setup task
	146	gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEmcalSetup.C");
	147	AliEmcalSetupTask *setupTask = AddTaskEmcalSetup();
7030f36f	148	setupTask->SetGeoPath("$ALICE_ROOT/OADB/EMCAL");
8628b70c	149
	150	UInt_t tpc = AliAnalysisTaskEmcal::kTPC;
	151
	152	// Tender Supplies
	153	if (useTender)
	154	{
	155	gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEMCALTender.C");
	156	AliAnalysisTaskSE *tender = AddTaskEMCALTender(runPeriod, kFALSE, kTRUE, kTRUE, kTRUE, kTRUE, kTRUE, kFALSE, kTRUE, kTRUE, kTRUE,
7030f36f	157	kFALSE, AliEMCALRecParam::kClusterizerNxN);
8628b70c	158	if (usedData != "AOD" && !useGrid) {
8628b70c	159	AliTender alitender = dynamic_cast<AliTender>(tender);
7030f36f	160	alitender->SetDefaultCDBStorage("local://$ALICE_ROOT/OCDB");
8628b70c	161	}
	162	}
	163
	164	// ################# Now: Call jet preparation macro (picotracks, hadronic corrected caloclusters, ...)
	165
	166	// Jet preparation
	167	gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskJetPreparation.C");
	168	AddTaskJetPreparation(dataType);
	169
	170	// ################# Now: Add jet finders+analyzers
	171
	172	gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C");
	173	AliEmcalJetTask* jetFinderTask = AddTaskEmcalJet("PicoTracks", "CaloClustersCorr", kANTIKT, 0.2, kCHARGEDJETS, 0.150, 0.300);
	174
	175	// Here you can put in your AddTaskMacro for your task
	176	gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJetSample.C");
95449deb	177	AliAnalysisTaskEmcalJetSample* anaTask = AddTaskEmcalJetSample("PicoTracks", "CaloClustersCorr", jetFinderTask->GetName(), "");
8628b70c	178
	179
	180	// Set the physics selection for all given tasks
	181	TObjArray *toptasks = mgr->GetTasks();
	182	for (Int_t i=0; i<toptasks->GetEntries(); ++i)
	183	{
	184	AliAnalysisTaskSE task = dynamic_cast<AliAnalysisTaskSE>(toptasks->At(i));
	185	if (!task)
	186	continue;
	187	if (task->InheritsFrom("AliPhysicsSelectionTask"))
	188	continue;
	189	::Info("setPSel", "Set physics selection for %s (%s)", task->GetName(), task->ClassName());
	190	task->SelectCollisionCandidates(pSel);
	191	}
	192
	193	mgr->SetUseProgressBar(1, 25);
	194
	195	if (!mgr->InitAnalysis())
	196	return;
	197	mgr->PrintStatus();
	198
	199	if (useGrid)
	200	{ // GRID CALCULATION
	201
	202	AliAnalysisGrid *plugin = CreateAlienHandler(uniqueName, gridDir, gridMode, runNumbers, pattern, additionalCXXs, additionalHs, maxFilesPerWorker, workerTTL, isMC);
	203	mgr->SetGridHandler(plugin);
	204
	205	// start analysis
	206	cout << "Starting GRID Analysis...";
	207	mgr->SetDebugLevel(0);
	208	mgr->StartAnalysis("grid");
	209	}
	210	else
	211	{ // LOCAL CALCULATION
	212
	213	TChain* chain = 0;
	214	if (usedData == "AOD")
	215	{
	216	gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/CreateAODChain.C");
	217	chain = CreateAODChain(localFiles.Data(), numLocalFiles);
	218	}
	219	else
	220	{ // ESD or skimmed ESD
	221	gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/CreateESDChain.C");
	222	chain = CreateESDChain(localFiles.Data(), numLocalFiles);
	223	}
	224
	225	// start analysis
	226	cout << "Starting LOCAL Analysis...";
	227	mgr->SetDebugLevel(2);
	228	mgr->StartAnalysis("local", chain);
	229	}
	230	}
	231
	232	//______________________________________________________________________________
	233	void LoadLibs()
	234	{
	235	// Load common libraries (better too many than too few)
	236	gSystem->Load("libTree");
	237	gSystem->Load("libVMC");
	238	gSystem->Load("libGeom");
	239	gSystem->Load("libGui");
	240	gSystem->Load("libXMLParser");
	241	gSystem->Load("libMinuit");
242	gSystem->Load("libMinuit2");
243	gSystem->Load("libProof");
244	gSystem->Load("libPhysics");
245	gSystem->Load("libSTEERBase");
246	gSystem->Load("libESD");
247	gSystem->Load("libAOD");
248	gSystem->Load("libOADB");
249	gSystem->Load("libANALYSIS");
250	gSystem->Load("libCDB");
251	gSystem->Load("libRAWDatabase");
252	gSystem->Load("libSTEER");
253	gSystem->Load("libEVGEN");
254	gSystem->Load("libANALYSISalice");
255	gSystem->Load("libCORRFW");
256	gSystem->Load("libTOFbase");
257	//gSystem->Load("libTOFrec");
258	gSystem->Load("libRAWDatabase.so");
259	gSystem->Load("libRAWDatarec.so");
260	gSystem->Load("libTPCbase.so");
261	gSystem->Load("libTPCrec.so");
262	gSystem->Load("libITSbase.so");
263	gSystem->Load("libITSrec.so");
264	gSystem->Load("libTRDbase.so");
265	gSystem->Load("libTENDER.so");
266	gSystem->Load("libSTAT.so");
267	gSystem->Load("libTRDrec.so");
268	gSystem->Load("libHMPIDbase.so");
269	gSystem->Load("libPWGPP.so");
270	gSystem->Load("libPWGHFbase");
271	gSystem->Load("libPWGDQdielectron");
272	gSystem->Load("libPWGHFhfe");
273	gSystem->Load("libEMCALUtils");
274	gSystem->Load("libPHOSUtils");
275	gSystem->Load("libPWGCaloTrackCorrBase");
276	gSystem->Load("libEMCALraw");
277	gSystem->Load("libEMCALbase");
278	gSystem->Load("libEMCALrec");
279	gSystem->Load("libTRDbase");
280	gSystem->Load("libVZERObase");
281	gSystem->Load("libVZEROrec");
282	gSystem->Load("libTENDER");
283	gSystem->Load("libTENDERSupplies");
284	gSystem->Load("libPWGEMCAL");
285	gSystem->Load("libPWGGAEMCALTasks");
286	gSystem->Load("libPWGTools");
287	gSystem->Load("libPWGCFCorrelationsBase");
288	gSystem->Load("libPWGCFCorrelationsDPhi");
289
290	//load CGAL, Fastjet and SISCone
291	gSystem->Load("libCGAL");
292	gSystem->Load("libfastjet");
293	gSystem->Load("libSISConePlugin");
294
295	gSystem->Load("libJETAN");
296	gSystem->Load("libFASTJETAN");
297	gSystem->Load("libPWGJEEMCALJetTasks");
298
299
300	// include paths
301	gSystem->AddIncludePath("-Wno-deprecated");
302	gSystem->AddIncludePath("-I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/EMCAL");
303	gSystem->AddIncludePath("-I$ALICE_ROOT/PWGDQ/dielectron -I$ALICE_ROOT/PWGHF/hfe");
304	gSystem->AddIncludePath("-I$ALICE_ROOT/JETAN -I$ALICE_ROOT/JETAN/fastjet");
305	}
306
307	AliAnalysisGrid* CreateAlienHandler(const char* uniqueName, const char* gridDir, const char* gridMode, const char* runNumbers,
308	const char* pattern, TString additionalCode, TString additionalHeaders, Int_t maxFilesPerWorker,
309	Int_t workerTTL, Bool_t isMC)
310	{
311	TDatime currentTime;
312	TString tmpName(uniqueName);
313
314	// Only add current date and time when not in terminate mode! In this case the exact name has to be supplied by the user
315	if(strcmp(gridMode, "terminate"))
316	{
317	tmpName += "_";
318	tmpName += currentTime.GetDate();
319	tmpName += "_";
320	tmpName += currentTime.GetTime();
321	}
322
323	TString tmpAdditionalLibs("");
324	tmpAdditionalLibs = Form("libTree.so libVMC.so libGeom.so libGui.so libXMLParser.so libMinuit.so libMinuit2.so libProof.so libPhysics.so libSTEERBase.so libESD.so libAOD.so libOADB.so libANALYSIS.so libCDB.so libRAWDatabase.so libSTEER.so libANALYSISalice.so libCORRFW.so libTOFbase.so libRAWDatabase.so libRAWDatarec.so libTPCbase.so libTPCrec.so libITSbase.so libITSrec.so libTRDbase.so libTENDER.so libSTAT.so libTRDrec.so libHMPIDbase.so libPWGPP.so libPWGHFbase.so libPWGDQdielectron.so libPWGHFhfe.so libEMCALUtils.so libPHOSUtils.so libPWGCaloTrackCorrBase.so libEMCALraw.so libEMCALbase.so libEMCALrec.so libTRDbase.so libVZERObase.so libVZEROrec.so libTENDER.so libTENDERSupplies.so libPWGEMCAL.so libPWGGAEMCALTasks.so libPWGTools.so libPWGCFCorrelationsBase.so libPWGCFCorrelationsDPhi.so libCGAL.so libJETAN.so libfastjet.so libSISConePlugin.so libFASTJETAN.so libPWGJE.so libPWGmuon.so libPWGJEEMCALJetTasks.so %s %s",additionalCode.Data(),additionalHeaders.Data());
325
326
327	TString macroName("");
328	TString execName("");
329	TString jdlName("");
330	macroName = Form("%s.C", tmpName.Data());
331	execName = Form("%s.sh", tmpName.Data());
332	jdlName = Form("%s.jdl", tmpName.Data());
333
334	AliAnalysisAlien *plugin = new AliAnalysisAlien();
335	plugin->SetOverwriteMode();
336	plugin->SetRunMode(gridMode);
337
338	// Here you can set the (Ali)ROOT version you want to use
339	plugin->SetAPIVersion("V1.1x");
340	plugin->SetROOTVersion("v5-34-02-1");
341	plugin->SetAliROOTVersion("v5-04-21-AN-1");
342
343	plugin->SetGridDataDir(gridDir); // e.g. "/alice/sim/LHC10a6"
344	plugin->SetDataPattern(pattern); //dir structure in run directory
345	if (!isMC)
346	plugin->SetRunPrefix("000");
347
348	plugin->AddRunList(runNumbers);
349
350	plugin->SetGridWorkingDir(Form("work/%s",tmpName.Data()));
351	plugin->SetGridOutputDir("output"); // In this case will be $HOME/work/output
352
353	plugin->SetAnalysisSource(additionalCode.Data());
354	plugin->SetAdditionalLibs(tmpAdditionalLibs.Data());
355	plugin->AddExternalPackage("boost::v1_43_0");
356	plugin->AddExternalPackage("cgal::v3.6");
357	plugin->AddExternalPackage("fastjet::v2.4.2");
358
359	plugin->SetDefaultOutputs(kTRUE);
360	//plugin->SetMergeExcludes("");
361	plugin->SetAnalysisMacro(macroName.Data());
362	plugin->SetSplitMaxInputFileNumber(maxFilesPerWorker);
363	plugin->SetExecutable(execName.Data());
364	plugin->SetTTL(workerTTL);
365	plugin->SetInputFormat("xml-single");
366	plugin->SetJDLName(jdlName.Data());
367	plugin->SetPrice(1);
368	plugin->SetSplitMode("se");
369
370	return plugin;
371	}