Use standard cuts.

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

void JetAnalysisManagerCAF()
{
    //
    if (gApplication) gApplication->InitializeGraphics();
    //
    // Connect to proof
    
//    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");

    //
    gProof->ShowEnabledPackages();
    
    //
    //
    // Create the chain
    //
    // TChain* chain = CreateESDChain("test.txt", 200);
 // Input 
    AliESDInputHandler* inpHandler = new AliESDInputHandler();
    //
    // Create the analysis manager
    //
    AliAODHandler* aodHandler   = new AliAODHandler();
    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->GetCommonInputContainer();

    AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
    coutput1->SetSpecialOutput();
    
    AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
							      AliAnalysisManager::kOutputContainer, "histos.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","/PWG4/arian/jetjetAbove_50_real");
}