X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ANALYSIS%2FrunAODFilterMC.C;h=e8d5d8e9ee6c8e833102a29b8bbb4c8317e6f4e9;hb=4010aa177a28d6461a17d169277d49cafc625b5c;hp=091470ef86cf322348a90a30a6aa9f9ceda1f76f;hpb=d2740fba2f4b5099b4b42fb8dc17a6f4e78aa501;p=u%2Fmrichter%2FAliRoot.git diff --git a/ANALYSIS/runAODFilterMC.C b/ANALYSIS/runAODFilterMC.C index 091470ef86c..e8d5d8e9ee6 100644 --- a/ANALYSIS/runAODFilterMC.C +++ b/ANALYSIS/runAODFilterMC.C @@ -8,8 +8,10 @@ void runAODFilterMC() gSystem->Load("libVMC.so"); bool bKineFilter = true; - + // TProof::Mgr("alicecaf")->SetROOTVersion("v5-21-01-alice_dbg"); TProof::Open("alicecaf"); + // gProof->SetParallel(1); + char *dataset = "/COMMON/COMMON/LHC08c11_10TeV_0.5T"; // char *dataset = "/PWG4/kleinb/LHC08q_jetjet100"; // gProof->ClearPackages(); @@ -67,11 +69,16 @@ void runAODFilterMC() esdTrackCutsL->SetRequireTPCRefit(kTRUE); esdTrackCutsL->SetDCAToVertexZ(3.0); esdTrackCutsL->SetDCAToVertexXY(3.0); + esdTrackCutsL->SetDCAToVertex2D(kTRUE); esdTrackCutsL->SetRequireSigmaToVertex(kFALSE); - esdTrackCutsL->SetAcceptKingDaughters(kFALSE); + esdTrackCutsL->SetAcceptKinkDaughters(kFALSE); + + AliESDtrackCuts* esdTrackCutsITSsa = new AliESDtrackCuts("AliESDtrackCuts", "ITS stand-alone"); + esdTrackCutsITSsa->SetRequireITSStandAlone(kTRUE); AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter"); trackFilter->AddCuts(esdTrackCutsL); + trackFilter->AddCuts(esdTrackCutsITSsa); AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter"); esdfilter->SetTrackFilter(trackFilter); @@ -81,17 +88,17 @@ void runAODFilterMC() // // Create containers for input/output - AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(), - AliAnalysisManager::kInputContainer); - - AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(), - AliAnalysisManager::kOutputContainer, "default"); + AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer(); + AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer(); coutput1->SetSpecialOutput(); if(bKineFilter){ mgr->ConnectInput (kinefilter, 0, cinput1 ); mgr->ConnectOutput (kinefilter, 0, coutput1 ); + AliAnalysisDataContainer *coutputEx = mgr->CreateContainer("cFilterList", TList::Class(), + AliAnalysisManager::kOutputContainer,"pyxsec_hists.root"); + mgr->ConnectOutput (kinefilter, 1, coutputEx ); } mgr->ConnectInput (esdfilter, 0, cinput1 ); @@ -105,3 +112,4 @@ void runAODFilterMC() mgr->StartAnalysis("proof",dataset,10000); } +