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->SetRequireTPCRefit(kTRUE);
73 esdTrackCutsL->SetMaxDCAToVertexXY(2.4);
74 esdTrackCutsL->SetMaxDCAToVertexZ(3.2);
75 esdTrackCutsL->SetDCAToVertex2D(kTRUE);
76 esdTrackCutsL->SetRequireSigmaToVertex(kFALSE);
77 esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);
80 AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
81 trackFilter->AddCuts(esdTrackCutsL);
83 AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
84 esdfilter->SetTrackFilter(trackFilter);
85 esdfilter->SetDebugLevel(10);
86 mgr->AddTask(esdfilter);
92 AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
93 jetana->SetDebugLevel(10);
94 jetana->SetConfigFile("ConfigJetAnalysisAOD.C");
97 // Create containers for input/output
98 AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
100 AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
101 coutput1->SetSpecialOutput();
103 AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
104 AliAnalysisManager::kOutputContainer, "histos.root");
106 mgr->ConnectInput (esdfilter, 0, cinput1 );
107 mgr->ConnectOutput (esdfilter, 0, coutput1);
109 mgr->ConnectInput (jetana, 0, cinput1);
110 mgr->ConnectOutput(jetana, 0, coutput1);
111 mgr->ConnectOutput(jetana, 1, coutput2);
119 mgr->StartAnalysis("proof","/PWG4/arian/jetjetAbove_50_real");