Example for analysis on kine tree.

diff --git a/JETAN/JetAnalysisManagerKine.C b/JETAN/JetAnalysisManagerKine.C
new file mode 100644 (file)
index 0000000..610e210
--- /dev/null
+++ b/JETAN/JetAnalysisManagerKine.C
@@ -0,0 +1,66 @@
+void JetAnalysisManagerKine()
+{
+      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("TE");
+    chain->Add("/home/morsch/AliRoot/newtrunk/FASTSIM/jets/galice.root");
+    /////////////////////////////////////////////////////////////////////////////////// 
+    // Create the analysis manager
+    //
+    // Output
+    AliAODHandler* aodHandler = new AliAODHandler();
+    aodHandler->SetOutputFileName("aod.root");
+    // MC Truth
+    AliMCEventHandler* mcHandler = new AliMCEventHandler();
+    AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
+    mgr->SetOutputEventHandler (aodHandler);
+    mgr->SetMCtruthEventHandler(mcHandler);
+    mgr->SetDebugLevel(10);
+
+    AliAnalysisTaskKineFilter *kinefilter = new AliAnalysisTaskKineFilter("Kine Filter");
+    mgr->AddTask(kinefilter);
+    
+    
+    AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
+    jetana->SetDebugLevel(10);
+    mgr->AddTask(jetana);
+
+    //
+    // 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 *coutput2 = mgr->CreateContainer("histos", TList::Class(),
+                                                             AliAnalysisManager::kOutputContainer, "histos.root");
+
+    mgr->ConnectInput  (kinefilter,     0, cinput1  );
+    mgr->ConnectOutput (kinefilter,     0, coutput1 );
+
+    mgr->ConnectInput  (jetana,     0, cinput1  );
+    mgr->ConnectOutput (jetana,     0, coutput1 );
+    mgr->ConnectOutput (jetana,     1, coutput2 );
+
+
+    //
+    // Run the analysis
+    //    
+    mgr->InitAnalysis();
+    mgr->PrintStatus();
+    mgr->StartAnalysis("local",chain);
+}