// AddTaskFullpAJets.C 2013-02-07 cyaldo

AliAnalysisTaskFullpAJets *AddTaskFullpAJets()
{
    const char *usedTracks="PicoTracks";
    const char *usedClusters="CaloClusters";
    const char *outClusName="CaloClustersCorr";
    const Double_t minTrackPt=0.15;
    const Double_t minClusterPt=0.30;
    const Double_t jetRadius=0.4;
    
    // Some constants for the jet finders
    const Int_t cKT                 = 0;
    const Int_t cANTIKT             = 1;
    const Int_t cFULLJETS           = 0;
    const Int_t cCHARGEDJETS        = 1;
    const Int_t cNEUTRALJETS        = 2;
    
    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
    if (!mgr)
    {
        Error("AddTaskJetCommon","No analysis manager found.");
        return 0;
    }

    // Determine the int of the jet radius for naming purposes
    Int_t drjet=Int_t(100*jetRadius);
    if (drjet%10 == 0)
    {
        drjet/=10;
    }
    
    TString taskName = Form("AnalysisFullpAJetsR%d",drjet);
    TString listName = Form("ListR%d",drjet);
    TString fileName = taskName + ".root";
    
    // Jet finders (RECONSTRUCTED DATA)
    TString tmpTaskName("");
    AliEmcalJetTask* jetFinderTask = NULL;

    gROOT->LoadMacro("$ALICE_ROOT/PWGJE/EMCALJetTasks/macros/AddTaskEmcalJet.C");
    
    // ########## CHARGED JETS ##########
    jetFinderTask = AddTaskEmcalJet(usedTracks,"",cANTIKT,jetRadius,cCHARGEDJETS,minTrackPt,minClusterPt);
    RequestMemory(jetFinderTask,250*1024);//more memory

    // ########## FULL JETS ##########
    // last two settings are for min pt tracks/clusters
    // anti-kT
    jetFinderTask = AddTaskEmcalJet(usedTracks,outClusName,cANTIKT,jetRadius,cFULLJETS,minTrackPt,minClusterPt);
    RequestMemory(jetFinderTask,250*1024);//more memory

    // kT
    jetFinderTask = AddTaskEmcalJet(usedTracks,outClusName,cKT,jetRadius,cFULLJETS,minTrackPt,minClusterPt);
    RequestMemory(jetFinderTask,250*1024);//more memory

    // Add User Task
    AliAnalysisTaskFullpAJets *task = new AliAnalysisTaskFullpAJets(taskName);
    mgr->AddTask(task);
    //task->SetR_JET(drjet);
    task->SetRjet(drjet);
    AliAnalysisDataContainer *coutput = mgr->CreateContainer(listName,TList::Class(),AliAnalysisManager::kOutputContainer,fileName);
    mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
    mgr->ConnectOutput(task,1,coutput);
    RequestMemory(task,250*1024);//more memory

    return task;

}