JETAN/JetAnalysisManagerNeutralsLoc.C

   1 void JetAnalysisManagerLoc()
   2 {
   3       gSystem->Load("libTree.so");
   4       gSystem->Load("libPhysics.so");
   5       gSystem->Load("libGeom.so");
   6       gSystem->Load("libVMC.so");
   7       gSystem->Load("libANALYSIS.so");
   8       gSystem->Load("libSTEERBase.so");
   9       gSystem->Load("libAOD.so");
  10       gSystem->Load("libESD.so");
  11       gSystem->Load("libANALYSISalice.so");
  12       gSystem->Load("libfastjet.so");
  13       gSystem->Load("libsiscone.so");
  14       gSystem->Load("libSISConePlugin.so");
  15       gSystem->Load("libJETAN.so");
  16
  17
  18     if (gApplication) gApplication->InitializeGraphics();
  19     //
  20     // Create the chain
  21     //
  22     gROOT->LoadMacro("CreateESDChain.C");
  23     TChain* chain = new TChain("esdTree");
  24     chain->Add("~/alice/data/highpt/kPythia6Jets125_150/030/AliESDs.root");
  25
  26     ///////////////////////////////////////////////////////////////////////////////////
  27     // Create the analysis manager
  28     //
  29     // Input
  30     AliESDInputHandler* inpHandler = new AliESDInputHandler();
  31     // Output
  32     AliAODHandler* aodHandler = new AliAODHandler();
  33     aodHandler->SetOutputFileName("aod.root");
  34     // MC Truth
  35     AliMCEventHandler* mcHandler = new AliMCEventHandler();
  36
  37
  38     AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
  39     mgr->SetInputEventHandler  (inpHandler);
  40     mgr->SetOutputEventHandler (aodHandler);
  41     mgr->SetMCtruthEventHandler(mcHandler);
  42     mgr->SetDebugLevel(10);
  43     ///////////////////////////////////////////////////////////////////////////////////
  44
  45
  46     //
  47     // Set of cuts
  48     //
  49     // standard
  50     AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Loose");
  51     esdTrackCutsL->SetMinNClustersTPC(50);
  52     esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
  53     esdTrackCutsL->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
  54     esdTrackCutsL->SetRequireTPCRefit(kTRUE);
  55     esdTrackCutsL->SetMaxNsigmaToVertex(3);
  56     esdTrackCutsL->SetRequireSigmaToVertex(kTRUE);
  57     esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);
  58     //
  59     // hard
  60     AliESDtrackCuts* esdTrackCutsH = new AliESDtrackCuts("AliESDtrackCuts", "Hard");
  61     esdTrackCutsH->SetMinNClustersTPC(100);
  62     esdTrackCutsH->SetMaxChi2PerClusterTPC(2.0);
  63     esdTrackCutsH->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
  64     esdTrackCutsH->SetRequireTPCRefit(kTRUE);
  65     esdTrackCutsH->SetMaxNsigmaToVertex(2);
  66     esdTrackCutsH->SetRequireSigmaToVertex(kTRUE);
  67     esdTrackCutsH->SetAcceptKinkDaughters(kFALSE);
  68     //
  69     AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
  70     trackFilter->AddCuts(esdTrackCutsL);
  71     trackFilter->AddCuts(esdTrackCutsH);
  72     //
  73     AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
  74     esdfilter->SetTrackFilter(trackFilter);
  75     esdfilter->SetDebugLevel(10);
  76     mgr->AddTask(esdfilter);
  77
  78     //
  79     // Jet analysis
  80     //
  81     // AliAnalysisTaskJets *jetana = new AliAnalysisTaskJetsV2("JetAnalysis");         // Charged only
  82     AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis",chain); // Charged + neutral
  83     jetana->SetDebugLevel(10);
  84
  85     AliAnalysisTaskJets *jetanaMC = new AliAnalysisTaskJets("JetAnalysisMC");
  86     jetanaMC->SetDebugLevel(10);
  87     jetanaMC->SetConfigFile("ConfigJetAnalysisMC.C");
  88     jetanaMC->SetNonStdBranch("jetsMC");
  89     mgr->AddTask(jetanaMC);
  90     mgr->AddTask(jetana);
  91
  92     //
  93     // Create containers for input/output
  94     AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
  95                                                              AliAnalysisManager::kInputContainer);
  96
  97     AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
  98                                                               AliAnalysisManager::kOutputContainer, "default");
  99     AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
 100                                                               AliAnalysisManager::kOutputContainer, "histos.root");
 101
 102     AliAnalysisDataContainer *coutputMC2 = mgr->CreateContainer("histosMC", TList::Class(),
 103                                                               AliAnalysisManager::kOutputContainer, "histosMC.root");
 104
 105     mgr->ConnectInput  (esdfilter,  0, cinput1  );
 106     mgr->ConnectOutput (esdfilter,  0, coutput1 );
 107
 108     mgr->ConnectInput  (jetana,     0, cinput1  );
 109     mgr->ConnectOutput (jetana,     0, coutput1 );
 110     mgr->ConnectOutput (jetana,     1, coutput2 );
 111
 112     mgr->ConnectInput  (jetanaMC,     0, cinput1  );
 113     mgr->ConnectOutput (jetanaMC,     0, coutput1 );
 114     mgr->ConnectOutput (jetanaMC,     1, coutputMC2 );
 115
 116
 117     //
 118     // Run the analysis
 119     //
 120     mgr->InitAnalysis();
 121     mgr->PrintStatus();
 122     mgr->StartAnalysis("local",chain);
 123 }