Add jet embedding MC (M. Verweij)

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

void JetAnalysisManager()
{
    //
    // Load relevant libraries
    //
    gSystem->Load("libTree.so");
    gSystem->Load("libNetx.so");
    gSystem->Load("libProof.so");
    gSystem->Load("libProofPlayer.so");
    gSystem->Load("libGeom.so");
    gSystem->Load("libEG");

    gSystem->Load("libANALYSIS.so");
    gSystem->Load("libESD.so"); 
    gSystem->Load("libJETAN.so");
    //
    // Connect to alien
    //
    TGrid::Connect("alien://"); 
    
    AliTagAnalysis *TagAna = new AliTagAnalysis(); 
    AliEventTagCuts *EvCuts = new AliEventTagCuts();
    AliRunTagCuts   *RuCuts = new AliRunTagCuts();
    //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("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);

    mgr->ConnectInput (jetana,0,cinput1);
    mgr->ConnectOutput(jetana,0,coutput1);
    cinput1->SetData(chain1);

//
// Run the analysis
//    

    if (mgr->InitAnalysis()) {
	mgr->PrintStatus();
	mgr->StartAnalysis("local", chain1);
    }
}