X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=JETAN%2FJetAnalysisManager.C;h=16bba38d64922d6bf116bcab2a5883775b10e042;hb=f11779ec42ceb686acee93673b486b28acc28f98;hp=bac6a188de7595a8ce4b550f327ee2d1f8186406;hpb=f8355f0993fee45acb48548b7d6bddf06ef65f7c;p=u%2Fmrichter%2FAliRoot.git diff --git a/JETAN/JetAnalysisManager.C b/JETAN/JetAnalysisManager.C index bac6a188de7..16bba38d649 100644 --- a/JETAN/JetAnalysisManager.C +++ b/JETAN/JetAnalysisManager.C @@ -3,41 +3,58 @@ void JetAnalysisManager() // // Load relevant libraries // - gSystem->Load("libEG.so"); - gSystem->Load("libGeom.so"); - gSystem->Load("libANALYSIS_NEW.so"); - gSystem->Load("libESD.so"); - gSystem->Load("libJETAN.so"); + gSystem->Load("libTree"); + gSystem->Load("libNetx"); + gSystem->Load("libProof"); + gSystem->Load("libProofPlayer"); + gSystem->Load("libGeom"); + gSystem->Load("libEG"); + + gSystem->Load("libANALYSIS"); + gSystem->Load("libESD"); + gSystem->Load("libJETAN"); // // Connect to alien // TGrid::Connect("alien://"); - // - // Prepare the input event chain - // + AliTagAnalysis *TagAna = new AliTagAnalysis(); - // create an EventTagCut object AliEventTagCuts *EvCuts = new AliEventTagCuts(); AliRunTagCuts *RuCuts = new AliRunTagCuts(); - TAlienCollection* coll = TAlienCollection::Open("tags/tag100.xml"); - TGridResult* TagResult = coll->GetGridResult(""); - TagAna->ChainGridTags(TagResult); - TChain* chain1 = 0x0; - chain1 = TagAna->QueryTags(RuCuts, EvCuts); + //EvCuts->SetNChargedAbove1GeVRange(1, 1000); + //EvCuts->SetMultiplicityRange(11,120); + //EvCuts->SetNPionRange(2,10000); + TAlienCollection* coll = TAlienCollection::Open("tag100.xml"); + TGridResult* TagResult = coll->GetGridResult("", 0, 0); + TagResult->Print(); + TagAna->ChainGridTags(TagResult); + + ////////////////////////////////////////////////////////////////// + //Get the chain + printf("*******************************\n"); + printf("*** Getting the Chain ***\n"); + printf("*******************************\n"); + TChain* chain1 = 0x0; + chain1 = TagAna->QueryTags(RuCuts, EvCuts); + chain1->ls(); + // // Make the analysis manager // - AliAnalysisManager *mgr = new AliAnalysisManager(); - AliAnalysisTask *jetana = new AliAnalysisTaskJets("JetAnalysis"); + AliAnalysisManager *mgr = new AliAnalysisManager("Manager", "Manager"); + mgr-> SetDebugLevel(10); + + AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis"); + jetana->SetDebugLevel(10); mgr->AddTask(jetana); // Create containers for input/output AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain1",TChain::Class(), AliAnalysisManager::kInputContainer); + AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist1", TH1::Class(), AliAnalysisManager::kOutputContainer); - printf("Connectiong I/O \n"); - + mgr->ConnectInput (jetana,0,cinput1); mgr->ConnectOutput(jetana,0,coutput1); cinput1->SetData(chain1); @@ -48,6 +65,6 @@ void JetAnalysisManager() if (mgr->InitAnalysis()) { mgr->PrintStatus(); - chain1->Process(mgr); + mgr->StartAnalysis("local", chain1); } }