HLTANALYSIS/JET/macros/tasks/JetAnalysisManagerHLTMC.C

   1 void JetAnalysisManagerHLTMC() {
   2
   3   if ( getenv("FASTJET") ) {
   4     gSystem->Load("libCGAL");
   5     gSystem->Load("libfastjet");
   6     gSystem->Load("libSISConePlugin");
   7   }
   8
   9   gSystem->Load("libTree");
  10   gSystem->Load("libPhysics");
  11   gSystem->Load("libGeom");
  12   gSystem->Load("libVMC");
  13   gSystem->Load("libANALYSIS");
  14   gSystem->Load("libSTEERBase");
  15   gSystem->Load("libAOD");
  16   gSystem->Load("libESD");
  17   gSystem->Load("libANALYSISalice");
  18   gSystem->Load("libJETAN");
  19
  20   gSystem->Load("libHLTbase");
  21   gSystem->Load("libAliHLTUtil");
  22   gSystem->Load("libAliHLTJET");
  23
  24   Int_t debugLevel = 10;
  25
  26   // --------------------------------------------------------------------------------
  27
  28   if (gApplication) gApplication->InitializeGraphics();
  29
  30   // --------------------------------------------------------------------------------
  31   //
  32   // Create the chain
  33   //
  34   // --------------------------------------------------------------------------------
  35
  36   gROOT->LoadMacro("${ALICE_ROOT}/PWG0/CreateESDChain.C");
  37   TChain* chain = new TChain("TE");
  38
  39   chain->Add("/home/jthaeder/jet/data/HEAD_2009-10-26/FastGen/kPythia6Jets125_150_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_200ev/galice.root");
  40
  41   // --------------------------------------------------------------------------------
  42   //
  43   // Create the analysis manager
  44   //
  45   // --------------------------------------------------------------------------------
  46
  47   // Output
  48   AliAODHandler* aodHandler = new AliAODHandler();
  49   aodHandler->SetOutputFileName("aod.root");
  50
  51   // MC Truth
  52   AliMCEventHandler* mcHandler = new AliMCEventHandler();
  53   mcHandler->SetReadTR(kFALSE); // for fastgen
  54
  55   AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
  56   mgr->SetOutputEventHandler (aodHandler);
  57   mgr->SetMCtruthEventHandler(mcHandler);
  58   mgr->SetDebugLevel(debugLevel);
  59
  60   // --------------------------------------------------------------------------------
  61   //
  62   // Set Configfiles
  63   //
  64   // --------------------------------------------------------------------------------
  65
  66   // -- HLT FFSC - MC
  67 #if FFSC
  68   AliAnalysisTaskJets *taskFFSC = new AliAnalysisTaskJets("JetAnalysisHLTMC_FFSC");
  69   taskFFSC->SetConfigFile("./tasks/ConfigJetAnalysisHLTMC.C");
  70   taskFFSC->SetNonStdBranch("jetsHLTMC_FFSC");
  71   taskFFSC->SetDebugLevel(debugLevel);
  72   mgr->AddTask(taskFFSC);
  73 #endif
  74
  75   // -- HLT FastJet - Kt
  76 #if FASTJET_KT
  77   if ( getenv("FASTJET") ) {
  78     AliAnalysisTaskJets *taskKt = new AliAnalysisTaskJets("JetAnalysisHLTMC_Kt");
  79     taskKt->SetConfigFile("./tasks/ConfigJetAnalysisHLTMCKt.C");
  80     taskKt->SetNonStdBranch("jetsHLTMC_Kt");
  81     taskKt->SetDebugLevel(debugLevel);
  82     mgr->AddTask(taskKt);
  83   }
  84 #endif
  85
  86 #if FASTJET_ANTIKT
  87   // -- HLT FastJet - AntiKt
  88   if ( getenv("FASTJET") ) {
  89     AliAnalysisTaskJets *taskAntiKt = new AliAnalysisTaskJets("JetAnalysisHLTMC_AntiKt");
  90     taskAntiKt->SetConfigFile("./tasks/ConfigJetAnalysisHLTMCAntiKt.C");
  91     taskAntiKt->SetNonStdBranch("jetsHLTMC_AntiKt");
  92     taskAntiKt->SetDebugLevel(debugLevel);
  93     mgr->AddTask(taskAntiKt);
  94   }
  95 #endif
  96
  97   // --------------------------------------------------------------------------------
  98   //
  99   // Create containers for input/output
 100   //
 101   // --------------------------------------------------------------------------------
 102
 103   AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
 104   AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
 105
 106 #if FFSC
 107   AliAnalysisDataContainer *coutputFFSC = mgr->CreateContainer("histos", TList::Class(),
 108                                                                AliAnalysisManager::kOutputContainer, "histos_MC_FFSC.root");
 109   mgr->ConnectInput  (taskFFSC,     0, cinput1  );
 110   mgr->ConnectOutput (taskFFSC,     0, coutput1 );
 111   mgr->ConnectOutput (taskFFSC,     1, coutputFFSC );
 112 #endif
 113
 114 #if FASTJET_KT
 115   if ( getenv("FASTJET") ) {
 116     AliAnalysisDataContainer *coutputKt = mgr->CreateContainer("histos", TList::Class(),
 117                                                                AliAnalysisManager::kOutputContainer, "histos_MC_KT.root");
 118     mgr->ConnectInput  (taskKt,     0, cinput1  );
 119     mgr->ConnectOutput (taskKt,     0, coutput1 );
 120     mgr->ConnectOutput (taskKt,     1, coutputKt );
 121   }
 122 #endif
 123
 124 #if FASTJET_ANTIKT
 125   if ( getenv("FASTJET") ) {
 126     AliAnalysisDataContainer *coutputAntiKt = mgr->CreateContainer("histos", TList::Class(),
 127                                                                    AliAnalysisManager::kOutputContainer, "histos_MC_ANTIKT.root");
 128     mgr->ConnectInput  (taskAntiKt,     0, cinput1  );
 129     mgr->ConnectOutput (taskAntiKt,     0, coutput1 );
 130     mgr->ConnectOutput (taskAntiKt,     1, coutputAntiKt );
 131   }
 132 #endif
 133
 134   //  AliAnalysisTaskKineFilter *kinefilter = new AliAnalysisTaskKineFilter("Kine Filter");
 135   //  mgr->AddTask(kinefilter);
 136   //  mgr->ConnectInput  (kinefilter,     0, cinput1  );
 137   //  mgr->ConnectOutput (kinefilter,     0, coutput1 );
 138
 139   // --------------------------------------------------------------------------------
 140   //
 141   // Run the analysis
 142   //
 143   // --------------------------------------------------------------------------------
 144
 145   mgr->InitAnalysis();
 146   mgr->PrintStatus();
 147   mgr->StartAnalysis("local",chain);
 148 }