From 435283af94ed6621a60d817062fb20a7fdf8ae6b Mon Sep 17 00:00:00 2001 From: mvl Date: Mon, 1 Dec 2014 12:28:32 +0100 Subject: [PATCH] Update from Tomas --- PWGJE/EMCALJetTasks/macros/AddTaskJetJTJT.C | 16 +++++++++++++++- 1 file changed, 15 insertions(+), 1 deletion(-) diff --git a/PWGJE/EMCALJetTasks/macros/AddTaskJetJTJT.C b/PWGJE/EMCALJetTasks/macros/AddTaskJetJTJT.C index 440738b0427..ba91c655ed0 100644 --- a/PWGJE/EMCALJetTasks/macros/AddTaskJetJTJT.C +++ b/PWGJE/EMCALJetTasks/macros/AddTaskJetJTJT.C @@ -1,17 +1,20 @@ // $Id$ AliAnalysisTaskJetJTJT* AddTaskJetJTJT( + const char *runPeriod = "LHC13b", const char *ntracks = "Tracks", const char *nclusters = "CaloClusters", const char *njets = "Jets", const char *nrho = "Rho", + Int_t trigger = AliVEvent::kMB, Int_t nCentBins = 1, Double_t jetradius = 0.2, Double_t jetptcut = 1, Double_t jetareacut = 0.6, const char *type = "EMCAL", Int_t leadhadtype = 0, - const char *taskname = "AliAnalysisTaskJetJTJT" + const char *taskname = "AliAnalysisTaskJetJTJT", + Int_t debug = 0 ) { // Get the pointer to the existing analysis manager via the static access method. @@ -48,6 +51,10 @@ AliAnalysisTaskJetJTJT* AddTaskJetJTJT( name += "_"; name += type; } + name += "_R"; + name += jetradius*10; + name += "_T"; + name += trigger; TString tracksName = "PicoTracks"; TString clustersName = "EmcCaloClusters"; @@ -66,7 +73,12 @@ AliAnalysisTaskJetJTJT* AddTaskJetJTJT( AliAnalysisTaskJetJTJT* jtTask = new AliAnalysisTaskJetJTJT(name); jtTask->SetCentRange(0.,100.); jtTask->SetNCentBins(nCentBins); + if(debug > 1){ + cout << "SetTrackArrayName: " << ntracks << endl; + } jtTask->SetTrackArrayName(ntracks); + jtTask->setDebug(debug); + jtTask->setRunPeriod(runPeriod); Double_t borders[5] = {0,10,20,40,100}; Double_t triggpt[8] = {0,5,10,20,40,80,100,150}; Double_t triggpta[2] = {0,100}; @@ -74,6 +86,7 @@ AliAnalysisTaskJetJTJT* AddTaskJetJTJT( jtTask->setCentBinBorders(5,borders); jtTask->setTriggPtBorders(8,triggpt); jtTask->setAssocPtBorders(2,triggpta); + jtTask->SelectCollisionCandidates(trigger); AliParticleContainer *trackCont = jtTask->AddParticleContainer(ntracks); trackCont->SetClassName("AliVTrack"); @@ -101,6 +114,7 @@ AliAnalysisTaskJetJTJT* AddTaskJetJTJT( AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer() ; TString contname(name); contname += "_histos"; + cout << "Create container " << contname << endl; AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(contname.Data(), TList::Class(),AliAnalysisManager::kOutputContainer, Form("%s", AliAnalysisManager::GetCommonFileName())); -- 2.39.3