from marta

[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / macros / AddTaskJetResponseMaker.C

diff --git a/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C b/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C
index 3458981125a2d7461a669092d20e7708f3d6930b..edb74a84ed715f8bd347e14284211c50bf986541 100644 (file)
--- a/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C
+++ b/PWGJE/EMCALJetTasks/macros/AddTaskJetResponseMaker.C
@@ -5,20 +5,23 @@ AliJetResponseMaker* AddTaskJetResponseMaker(
   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    jetradius          = 0.2,
+  Double_t    jetradius2         = 0.2,
   Double_t    jetptcut           = 1,
   Double_t    jetareacut         = 0.557,
   Double_t    jetBiasTrack       = 5,
-  Double_t    jetBiasClus        = 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
@@ -46,16 +49,19 @@ AliJetResponseMaker* AddTaskJetResponseMaker(
   //-------------------------------------------------------
 
   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::kUser) 
     name += "USER";
-  if (ptHardBin != -999) {
-    name += "_PtHard";
-    name += ptHardBin;
-  }
+
+  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)
@@ -68,11 +74,12 @@ AliJetResponseMaker* AddTaskJetResponseMaker(
   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->SetJetRadius(jetradius);
+  jetTask->SetJet2Radius(jetradius2);
   jetTask->SetJetPtCut(jetptcut);
   jetTask->SetPercAreaCut(jetareacut);
   jetTask->SetPtBiasJetTrack(jetBiasTrack);
@@ -80,6 +87,7 @@ AliJetResponseMaker* AddTaskJetResponseMaker(
   jetTask->SetMatching(matching, maxDistance1, maxDistance2);
   jetTask->SetVzRange(-10,10);
   jetTask->SetPtHardBin(ptHardBin);
+  jetTask->SetCentRange(minCent,maxCent);
 
   if (biggerMatrix) 
     jetTask->SetHistoBins(1000,0,500);