// $Id$
AliJetResponseMaker* AddTaskJetResponseMaker(
- const char *ntracks = "Tracks",
- const char *nclusters = "CaloClusters",
- const char *njets = "Jets",
- const char *nmcjets = "MCJets",
- const char *nmctracks = "MCParticles",
- Double_t jetradius = 0.4,
+ const char *ntracks1 = "Tracks",
+ const char *nclusters1 = "CaloClusters",
+ const char *njets1 = "Jets",
+ const char *nrho1 = "Rho",
+ Double_t jetradius1 = 0.2,
+ const char *ntracks2 = "MCParticles",
+ const char *nclusters2 = "",
+ const char *njets2 = "MCJets",
+ const char *nrho2 = "",
+ Double_t jetradius2 = 0.2,
Double_t jetptcut = 1,
- Double_t jetareacut = 0.8,
- Double_t ptcut = 0.15,
- Double_t jetBiasTrack = 10,
- Double_t jetBiasClus = 10,
- Double_t maxDistance = 0.25,
+ Double_t jetareacut = 0.557,
+ Double_t jetBiasTrack = 5,
+ Double_t jetBiasClus = 1000,
+ UInt_t matching = AliJetResponseMaker::kGeometrical,
+ Double_t maxDistance1 = 0.25,
+ Double_t maxDistance2 = 0.25,
UInt_t type = AliAnalysisTaskEmcal::kTPC,
Int_t ptHardBin = -999,
+ Double_t minCent = -999,
+ Double_t maxCent = -999,
const char *taskname = "AliJetResponseMaker",
+ Bool_t biggerMatrix = kFALSE,
AliJetResponseMaker* address = 0
)
{
// Init the task and do settings
//-------------------------------------------------------
- TString name(taskname);
- name += "_";
- name += njets;
- name += "_Track";
- name += jetBiasTrack;
- name += "_Clus";
- name += jetBiasClus;
- name += "_R0";
- name += floor(jetradius*10+0.5);
- name += "_PtCut";
- name += floor(ptcut*1000+0.5);
- name += "_";
- if (type == AliAnalysisTaskEmcal::kTPC)
+ TString name(Form("%s_%s_%s_Track%d_Clus%d_",taskname,njets1,njets2,(Int_t)floor(jetBiasTrack),(Int_t)floor(jetBiasClus)));
+
+ if (type == AliAnalysisTaskEmcal::kTPC)
name += "TPC";
else if (type == AliAnalysisTaskEmcal::kEMCAL)
name += "EMCAL";
- else if (type == AliAnalysisTaskEmcal::kTPCSmall)
- name += "TPCSmall";
- else if (type == AliAnalysisTaskEmcal::kEMCALOnly)
- name += "EMCALOnly";
- if (ptHardBin != -999) {
- name += "_PtHard";
- name += ptHardBin;
- }
+ else if (type == AliAnalysisTaskEmcal::kUser)
+ name += "USER";
+
+ if (minCent != -999 && maxCent != -999)
+ name += Form("_Cent%d_%d", (Int_t)floor(minCent), (Int_t)floor(maxCent));
+
+ if (ptHardBin != -999)
+ name += Form("_PtHard%d", ptHardBin);
AliJetResponseMaker* jetTask = address;
if (jetTask)
jetTask = new AliJetResponseMaker(name);
jetTask->SetAnaType(type);
- jetTask->SetTracksName(ntracks);
- jetTask->SetClusName(nclusters);
- jetTask->SetJetsName(njets);
- jetTask->SetMCJetsName(nmcjets);
- jetTask->SetMCTracksName(nmctracks);
- jetTask->SetPtCut(ptcut);
- jetTask->SetJetRadius(jetradius);
+ jetTask->SetTracksName(ntracks1);
+ jetTask->SetClusName(nclusters1);
+ jetTask->SetJetsName(njets1);
+ jetTask->SetRhoName(nrho1);
+ jetTask->SetJetRadius(jetradius1);
+ jetTask->SetTracks2Name(ntracks2);
+ jetTask->SetClus2Name(nclusters2);
+ jetTask->SetJets2Name(njets2);
+ jetTask->SetRho2Name(nrho2);
+ jetTask->SetJet2Radius(jetradius2);
jetTask->SetJetPtCut(jetptcut);
jetTask->SetPercAreaCut(jetareacut);
jetTask->SetPtBiasJetTrack(jetBiasTrack);
jetTask->SetPtBiasJetClus(jetBiasClus);
- jetTask->SetMaxDistance(maxDistance);
+ jetTask->SetMatching(matching, maxDistance1, maxDistance2);
jetTask->SetVzRange(-10,10);
jetTask->SetPtHardBin(ptHardBin);
+ jetTask->SetCentRange(minCent,maxCent);
+
+ if (biggerMatrix)
+ jetTask->SetHistoBins(1000,0,500);
//-------------------------------------------------------
// Final settings, pass to manager and set the containers