void JetAnalysisManagerLoc()
{
- gSystem->Load("libTree.so");
- gSystem->Load("libGeom.so");
- gSystem->Load("libVMC.so");
- gSystem->Load("libANALYSIS.so");
- gSystem->Load("libSTEERBase.so");
- gSystem->Load("libAOD.so");
- gSystem->Load("libESD.so");
- gSystem->Load("libJETAN.so");
+ gSystem->Load("libTree.so");
+ gSystem->Load("libPhysics.so");
+ gSystem->Load("libGeom.so");
+ gSystem->Load("libVMC.so");
+ gSystem->Load("libANALYSIS.so");
+ gSystem->Load("libSTEERBase.so");
+ gSystem->Load("libAOD.so");
+ gSystem->Load("libESD.so");
+ gSystem->Load("libANALYSISalice.so");
+ gSystem->Load("libJETAN.so");
+
//
if (gApplication) gApplication->InitializeGraphics();
// Create the chain
//
+ gROOT->LoadMacro("CreateESDChain.C");
TChain* chain = new TChain("esdTree");
- chain->Add("/home/morsch/AliRoot/data/data_jets102/AliESDs.root");
- chain->Add("/home/morsch/AliRoot/data/data_jets103/AliESDs.root");
- //
+ chain->Add("~/alice/data/highpt/kPythia6Jets125_150/030/AliESDs.root");
+
+ ///////////////////////////////////////////////////////////////////////////////////
// Create the analysis manager
//
- AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager");
- mgr->SetDebugLevel(10);
- //
- // Common output service
- AliAODHandler* aodHandler = new AliAODHandler();
+ // Input
+ AliESDInputHandler* inpHandler = new AliESDInputHandler();
+ // Output
+ AliAODHandler* aodHandler = new AliAODHandler();
aodHandler->SetOutputFileName("aod.root");
- mgr->SetOutputEventHandler(aodHandler);
- //
- // Common MC truth services
+ // MC Truth
AliMCEventHandler* mcHandler = new AliMCEventHandler();
+
+
+ AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager");
+ mgr->SetInputEventHandler (inpHandler);
+ mgr->SetOutputEventHandler (aodHandler);
mgr->SetMCtruthEventHandler(mcHandler);
+ mgr->SetDebugLevel(10);
+ ///////////////////////////////////////////////////////////////////////////////////
+
+
+ //
+ // Set of cuts
+ //
+ // standard
+ AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("Standard Track Cuts", "ESD Track Cuts");
+ esdTrackCutsL->SetMinNClustersTPC(50);
+ esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
+ esdTrackCutsL->SetRequireTPCRefit(kTRUE);
+ esdTrackCutsL->SetMaxDCAToVertexXY(2.4);
+ esdTrackCutsL->SetMaxDCAToVertexZ(3.2);
+ esdTrackCutsL->SetDCAToVertex2D(kTRUE);
+ esdTrackCutsL->SetRequireSigmaToVertex(kFALSE);
+ esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);
+
+ //
+ // hard
+ AliESDtrackCuts* esdTrackCutsH = new AliESDtrackCuts("AliESDtrackCuts", "Hard");
+ esdTrackCutsH->SetMinNClustersTPC(100);
+ esdTrackCutsH->SetMaxChi2PerClusterTPC(2.0);
+ esdTrackCutsH->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
+ esdTrackCutsH->SetRequireTPCRefit(kTRUE);
+ esdTrackCutsH->SetMaxNsigmaToVertex(2);
+ esdTrackCutsH->SetRequireSigmaToVertex(kTRUE);
+ esdTrackCutsH->SetAcceptKinkDaughters(kFALSE);
//
- // Jet analysis
+ AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
+ trackFilter->AddCuts(esdTrackCutsL);
+ trackFilter->AddCuts(esdTrackCutsH);
+ //
+ AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
+ esdfilter->SetTrackFilter(trackFilter);
+ esdfilter->SetDebugLevel(10);
+ mgr->AddTask(esdfilter);
+
AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
jetana->SetDebugLevel(10);
+
+
+
+ AliAnalysisTaskJets *jetanaMC = new AliAnalysisTaskJets("JetAnalysisMC");
+ jetanaMC->SetDebugLevel(10);
+ jetanaMC->SetConfigFile("ConfigJetAnalysisMC.C");
+ jetanaMC->SetNonStdBranch("jetsMC");
+ mgr->AddTask(jetanaMC);
mgr->AddTask(jetana);
//
// Create containers for input/output
- AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
- AliAnalysisManager::kInputContainer);
+ AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
+ AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
+ AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
+ AliAnalysisManager::kOutputContainer, "histos.root");
- AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
- AliAnalysisManager::kOutputContainer, "aod.root");
+ AliAnalysisDataContainer *coutputMC2 = mgr->CreateContainer("histosMC", TList::Class(),
+ AliAnalysisManager::kOutputContainer, "histosMC.root");
+
+ mgr->ConnectInput (esdfilter, 0, cinput1 );
+ mgr->ConnectOutput (esdfilter, 0, coutput1 );
+
+ mgr->ConnectInput (jetana, 0, cinput1 );
+ mgr->ConnectOutput (jetana, 0, coutput1 );
+ mgr->ConnectOutput (jetana, 1, coutput2 );
+
+ mgr->ConnectInput (jetanaMC, 0, cinput1 );
+ mgr->ConnectOutput (jetanaMC, 0, coutput1 );
+ mgr->ConnectOutput (jetanaMC, 1, coutputMC2 );
- AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TH1F::Class(),
- AliAnalysisManager::kOutputContainer, "histos.root");
- mgr->ConnectInput (jetana, 0, cinput1 );
- mgr->ConnectOutput (jetana, 0, coutput1 );
- mgr->ConnectOutput (jetana, 1, coutput2 );
//
// Run the analysis
//