Phys DA , works on DAQ tst machine

[u/mrichter/AliRoot.git] / JETAN / JetAnalysisManagerCAF.C

diff --git a/JETAN/JetAnalysisManagerCAF.C b/JETAN/JetAnalysisManagerCAF.C
index 3916ea09f0be119d495164b08b5a07f5439fab57..d387af51ee74d09043af3f3bc0832df0e9d3caf4 100644 (file)
--- a/JETAN/JetAnalysisManagerCAF.C
+++ b/JETAN/JetAnalysisManagerCAF.C
@@ -2,30 +2,41 @@ void JetAnalysisManagerCAF()
 {
     //
     if (gApplication) gApplication->InitializeGraphics();
-    gROOT->LoadMacro("CreateESDChain.C");
     //
     // Connect to proof
-
-    TProof::Reset("proof://morsch@lxb6046.cern.ch"); 
-    TProof::Open("proof://morsch@lxb6046.cern.ch"); 
-
-    //gProof->SetParallel(3);
-    gProof->ClearPackage("ESD");
-    gProof->ClearPackage("JETAN");
-    gProof->ClearPackage("ANALYSIS");
+    
+//    TProof::Reset("proof://morsch@lxb6046.cern.ch"); 
+    TProof::Open("proof://morsch@lxb6046.cern.ch");
+//    gProof->SetParallel(1);
+//    gProof->ClearPackage("STEERBase");
+//    gProof->ClearPackage("ESD");
+//    gProof->ClearPackage("AOD");
+//    gProof->ClearPackage("JETAN");
+//    gProof->ClearPackage("ANALYSIS");
+//    gProof->ClearPackage("ANALYSISalice");
     
     gProof->ShowEnabledPackages();
+
+    // Enable the STEERBase Package
+    gProof->UploadPackage("STEERBase.par");
+    gProof->EnablePackage("STEERBase");
     // Enable the ESD Package
     gProof->UploadPackage("ESD.par");
     gProof->EnablePackage("ESD");
+     // Enable the AOD Package
+    gProof->UploadPackage("AOD.par");
+    gProof->EnablePackage("AOD");
      // Enable the Analysis Package
     gProof->UploadPackage("ANALYSIS.par");
     gProof->EnablePackage("ANALYSIS");
+
+    gProof->UploadPackage("ANALYSISalice.par");
+    gProof->EnablePackage("ANALYSISalice");
+
     // Enable the JETAN Package
     gProof->UploadPackage("JETAN.par");
     gProof->EnablePackage("JETAN");
-    // Load Configuration macro
-    //gProof->Load("ConfigJetAnalysis.C");
+
     //
     gProof->ShowEnabledPackages();
     
@@ -33,32 +44,77 @@ void JetAnalysisManagerCAF()
     //
     // Create the chain
     //
-    TChain* chain = CreateESDChain("test.txt", 200);
+    // TChain* chain = CreateESDChain("test.txt", 200);
+ // Input 
+    AliESDInputHandler* inpHandler = new AliESDInputHandler();
     //
-    // Make the analysis manager
+    // Create the analysis manager
     //
-    AliAnalysisManager *mgr     = new AliAnalysisManager("Jet Manager", "Jet Manager");
     AliAODHandler* aodHandler   = new AliAODHandler();
-    mgr->SetEventHandler(aodHandler);
-    mgr-> SetDebugLevel(10);
+    aodHandler->SetOutputFileName("jets.root");
     
+    AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
+    mgr->SetOutputEventHandler(aodHandler);
+    mgr->SetInputEventHandler(inpHandler);
+    mgr-> SetDebugLevel(10);
+
+
+    //
+    //  ESD Filter Task
+    //
+    //
+    // Set of cuts
+    // 
+    // standard
+    AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Loose");
+    esdTrackCutsL->SetMinNClustersTPC(50);
+    esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
+    esdTrackCutsL->SetRequireTPCRefit(kTRUE);
+    esdTrackCutsL->SetMaxDCAToVertexXY(2.4);
+    esdTrackCutsL->SetMaxDCAToVertexZ(3.2);
+    esdTrackCutsL->SetDCAToVertex2D(kTRUE);
+    esdTrackCutsL->SetRequireSigmaToVertex(kFALSE);
+    esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);
+    //
+    //
+    AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
+    trackFilter->AddCuts(esdTrackCutsL);
+    //
+    AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
+    esdfilter->SetTrackFilter(trackFilter);
+    esdfilter->SetDebugLevel(10);
+    mgr->AddTask(esdfilter);
+
+//
+//   Jet Finder Task
+//
+
     AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
     jetana->SetDebugLevel(10);
+    jetana->SetConfigFile("ConfigJetAnalysisAOD.C");
     mgr->AddTask(jetana);
     //
     // Create containers for input/output
-    AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(), 
-                                                            AliAnalysisManager::kInputContainer);
+    AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
+
+    AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
+    coutput1->SetSpecialOutput();
+    
+    AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
+                                                             AliAnalysisManager::kOutputContainer, "histos.root");
 
-    AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
-                                                             AliAnalysisManager::kOutputContainer, "aod.root");
+    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
     //    
     mgr->InitAnalysis();
     mgr->PrintStatus();
-    mgr->StartAnalysis("proof",chain);
+    mgr->StartAnalysis("proof","/PWG4/arian/jetjetAbove_50_real");
 }