Added some more scripts

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

diff --git a/ANALYSIS/runAODFilterMC.C b/ANALYSIS/runAODFilterMC.C
index b7e656d56ee75255b205de63c09cd8a77bd0edd8..e8d5d8e9ee6c8e833102a29b8bbb4c8317e6f4e9 100644 (file)
--- a/ANALYSIS/runAODFilterMC.C
+++ b/ANALYSIS/runAODFilterMC.C
@@ -69,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);
@@ -83,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  );
@@ -107,3 +112,4 @@ void runAODFilterMC()
       mgr->StartAnalysis("proof",dataset,10000);
 
 }
+