Correction.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 22 Oct 2007 17:30:47 +0000 (17:30 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 22 Oct 2007 17:30:47 +0000 (17:30 +0000)
JETAN/JetAnalysisManagerLoc.C

index bf6fb3a..fc966a9 100644 (file)
@@ -1,37 +1,74 @@
 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("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("libPWG0base.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");
-    //
+
+    /////////////////////////////////////////////////////////////////////////////////// 
     // 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();
-    mgr->SetMCtruthEventHandler(mcHandler);
+    
 
+    AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
+    mgr->SetInputEventHandler  (inpHandler);
+    mgr->SetOutputEventHandler (aodHandler);
+    mgr->SetMCtruthEventHandler(mcHandler);
+    mgr-> SetDebugLevel(10);
+    /////////////////////////////////////////////////////////////////////////////////// 
+    
+    
     //
-    // Jet analysis
+    // Set of cuts
+    // 
+    // standard
+    AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Loose");
+    esdTrackCutsL->SetMinNClustersTPC(50);
+    esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
+    esdTrackCutsL->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
+    esdTrackCutsL->SetRequireTPCRefit(kTRUE);
+    esdTrackCutsL->SetMinNsigmaToVertex(3);
+    esdTrackCutsL->SetRequireSigmaToVertex(kTRUE);
+    esdTrackCutsL->SetAcceptKingDaughters(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->SetMinNsigmaToVertex(2);
+    esdTrackCutsH->SetRequireSigmaToVertex(kTRUE);
+    esdTrackCutsH->SetAcceptKingDaughters(kFALSE);
+    //
+    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);
     mgr->AddTask(jetana);
@@ -42,14 +79,18 @@ void JetAnalysisManagerLoc()
                                                             AliAnalysisManager::kInputContainer);
 
     AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
-                                                             AliAnalysisManager::kOutputContainer, "aod.root");
-
-    AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TH1F::Class(),
+                                                             AliAnalysisManager::kOutputContainer, "default");
+    AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
                                                              AliAnalysisManager::kOutputContainer, "histos.root");
 
-    mgr->ConnectInput  (jetana,  0, cinput1  );
-    mgr->ConnectOutput (jetana,  0, coutput1 );
-    mgr->ConnectOutput (jetana,  1, coutput2 );
+    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 );
+
+
     //
     // Run the analysis
     //