Merge branch 'master' of https://git.cern.ch/reps/AliRoot

[u/mrichter/AliRoot.git] / ANALYSIS / runAODFilterMC.C

diff --git a/ANALYSIS/runAODFilterMC.C b/ANALYSIS/runAODFilterMC.C
index 091470ef86cf322348a90a30a6aa9f9ceda1f76f..73ed697e416541291699208a076735a12fdfca95 100644 (file)
--- a/ANALYSIS/runAODFilterMC.C
+++ b/ANALYSIS/runAODFilterMC.C
@@ -2,14 +2,16 @@ void runAODFilterMC()
 {
   // PROOF example
 
-      gSystem->Load("libTree.so");
-      gSystem->Load("libPhysics.so");
-      gSystem->Load("libGeom.so");
-      gSystem->Load("libVMC.so");
+      gSystem->Load("libTree");
+      gSystem->Load("libPhysics");
+      gSystem->Load("libGeom");
+      gSystem->Load("libVMC");
 
       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);
 
 }
+