1 void JetAnalysisManagerCAF()
4 if (gApplication) gApplication->InitializeGraphics();
8 // TProof::Reset("proof://morsch@lxb6046.cern.ch");
9 TProof::Open("proof://morsch@lxb6046.cern.ch");
10 // gProof->SetParallel(1);
11 // gProof->ClearPackage("STEERBase");
12 // gProof->ClearPackage("ESD");
13 // gProof->ClearPackage("AOD");
14 // gProof->ClearPackage("JETAN");
15 // gProof->ClearPackage("ANALYSIS");
16 // gProof->ClearPackage("ANALYSISalice");
18 gProof->ShowEnabledPackages();
20 // Enable the STEERBase Package
21 gProof->UploadPackage("STEERBase.par");
22 gProof->EnablePackage("STEERBase");
23 // Enable the ESD Package
24 gProof->UploadPackage("ESD.par");
25 gProof->EnablePackage("ESD");
26 // Enable the AOD Package
27 gProof->UploadPackage("AOD.par");
28 gProof->EnablePackage("AOD");
29 // Enable the Analysis Package
30 gProof->UploadPackage("ANALYSIS.par");
31 gProof->EnablePackage("ANALYSIS");
33 gProof->UploadPackage("ANALYSISalice.par");
34 gProof->EnablePackage("ANALYSISalice");
36 // Enable the JETAN Package
37 gProof->UploadPackage("JETAN.par");
38 gProof->EnablePackage("JETAN");
41 gProof->ShowEnabledPackages();
47 // TChain* chain = CreateESDChain("test.txt", 200);
49 AliESDInputHandler* inpHandler = new AliESDInputHandler();
51 // Create the analysis manager
53 AliAODHandler* aodHandler = new AliAODHandler();
54 aodHandler->SetOutputFileName("jets.root");
56 AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager");
57 mgr->SetOutputEventHandler(aodHandler);
58 mgr->SetInputEventHandler(inpHandler);
59 mgr-> SetDebugLevel(10);
69 AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Loose");
70 esdTrackCutsL->SetMinNClustersTPC(50);
71 esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
72 esdTrackCutsL->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
73 esdTrackCutsL->SetRequireTPCRefit(kTRUE);
74 esdTrackCutsL->SetMinNsigmaToVertex(3);
75 esdTrackCutsL->SetRequireSigmaToVertex(kTRUE);
76 esdTrackCutsL->SetAcceptKingDaughters(kFALSE);
79 AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
80 trackFilter->AddCuts(esdTrackCutsL);
82 AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
83 esdfilter->SetTrackFilter(trackFilter);
84 esdfilter->SetDebugLevel(10);
85 mgr->AddTask(esdfilter);
91 AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
92 jetana->SetDebugLevel(10);
93 jetana->SetConfigFile("ConfigJetAnalysisAOD.C");
96 // Create containers for input/output
97 AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
99 AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
100 coutput1->SetSpecialOutput();
102 AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
103 AliAnalysisManager::kOutputContainer, "histos.root");
105 mgr->ConnectInput (esdfilter, 0, cinput1 );
106 mgr->ConnectOutput (esdfilter, 0, coutput1);
108 mgr->ConnectInput (jetana, 0, cinput1);
109 mgr->ConnectOutput(jetana, 0, coutput1);
110 mgr->ConnectOutput(jetana, 1, coutput2);
118 mgr->StartAnalysis("proof","/PWG4/arian/jetjetAbove_50_real");