Improving the performance of the bkg estimators (#13796)

Author: nzardosh

PWGJE/TableProducer/eventwiseConstituentSubtractor.cxx

@@ -80,6 +80,12 @@ struct eventWiseConstituentSubtractorTask {
   Configurable<float> rMax{"rMax", 0.24, "maximum distance of subtraction"};
   Configurable<float> eventEtaMax{"eventEtaMax", 0.9, "maximum pseudorapidity of event"};
   Configurable<bool> doRhoMassSub{"doRhoMassSub", true, "perfom mass subtraction as well"};
+  Configurable<double> ghostRapMax{"ghostRapMax", 0.9, "Ghost rapidity max"};
+  Configurable<int> ghostRepeat{"ghostRepeat", 1, "Ghost tiling repeats"};
+  Configurable<double> ghostArea{"ghostArea", 0.005, "Area per ghost"};
+  Configurable<double> ghostGridScatter{"ghostGridScatter", 1.0, "Grid scatter"};
+  Configurable<double> ghostKtScatter{"ghostKtScatter", 0.1, "kT scatter"};
+  Configurable<double> ghostMeanPt{"ghostMeanPt", 1e-100, "Mean ghost pT"};
 
   JetBkgSubUtils eventWiseConstituentSubtractor;
   std::vector<fastjet::PseudoJet> inputParticles;
@@ -98,6 +104,9 @@ struct eventWiseConstituentSubtractorTask {
     eventWiseConstituentSubtractor.setDoRhoMassSub(doRhoMassSub);
     eventWiseConstituentSubtractor.setConstSubAlphaRMax(alpha, rMax);
     eventWiseConstituentSubtractor.setMaxEtaEvent(eventEtaMax);
+    fastjet::GhostedAreaSpec ghostAreaSpec(ghostRapMax, ghostRepeat, ghostArea,
+                                           ghostGridScatter, ghostKtScatter, ghostMeanPt);
+    eventWiseConstituentSubtractor.setGhostAreaSpec(ghostAreaSpec);
 
     if (applyTrackingEfficiency) {
       if (trackingEfficiencyPtBinning->size() < 2) {

PWGJE/TableProducer/rhoEstimator.cxx

@@ -86,9 +86,15 @@ struct RhoEstimatorTask {
     Configurable<float> bkgjetR{"bkgjetR", 0.2, "jet resolution parameter for determining background density"};
     Configurable<float> bkgEtaMin{"bkgEtaMin", -0.7, "minimim pseudorapidity for determining background density"};
     Configurable<float> bkgEtaMax{"bkgEtaMax", 0.7, "maximum pseudorapidity for determining background density"};
-    Configurable<float> bkgPhiMin{"bkgPhiMin", -99., "minimim phi for determining background density"};
-    Configurable<float> bkgPhiMax{"bkgPhiMax", 99., "maximum phi for determining background density"};
+    Configurable<float> bkgPhiMin{"bkgPhiMin", -6.283, "minimim phi for determining background density"};
+    Configurable<float> bkgPhiMax{"bkgPhiMax", 6.283, "maximum phi for determining background density"};
     Configurable<bool> doSparse{"doSparse", false, "perfom sparse estimation"};
+    Configurable<double> ghostRapMax{"ghostRapMax", 0.9, "Ghost rapidity max"};
+    Configurable<int> ghostRepeat{"ghostRepeat", 1, "Ghost tiling repeats"};
+    Configurable<double> ghostArea{"ghostArea", 0.005, "Area per ghost"};
+    Configurable<double> ghostGridScatter{"ghostGridScatter", 1.0, "Grid scatter"};
+    Configurable<double> ghostKtScatter{"ghostKtScatter", 0.1, "kT scatter"};
+    Configurable<double> ghostMeanPt{"ghostMeanPt", 1e-100, "Mean ghost pT"};
 
     Configurable<float> thresholdTriggerTrackPtMin{"thresholdTriggerTrackPtMin", 0.0, "Minimum trigger track pt to accept event"};
     Configurable<float> thresholdClusterEnergyMin{"thresholdClusterEnergyMin", 0.0, "Minimum cluster energy to accept event"};
@@ -134,6 +140,11 @@ struct RhoEstimatorTask {
       bkgPhiMin_ = -2.0 * M_PI;
     }
     bkgSub.setPhiMinMax(bkgPhiMin_, bkgPhiMax_);
+
+    fastjet::GhostedAreaSpec ghostAreaSpec(config.ghostRapMax, config.ghostRepeat, config.ghostArea,
+                                           config.ghostGridScatter, config.ghostKtScatter, config.ghostMeanPt);
+    bkgSub.setGhostAreaSpec(ghostAreaSpec);
+
     eventSelectionBits = jetderiveddatautilities::initialiseEventSelectionBits(static_cast<std::string>(config.eventSelections));
     triggerMaskBits = jetderiveddatautilities::initialiseTriggerMaskBits(config.triggerMasks);
 

PWGJE/Tasks/jetTutorial.cxx

@@ -215,6 +215,7 @@ struct JetTutorialTask {
 
   using JetMCPTable = soa::Filtered<soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents, aod::ChargedMCParticleLevelJetsMatchedToChargedMCDetectorLevelJets>>;
   void processMCMatchedCharged(soa::Filtered<aod::JetCollisionsMCD>::iterator const& collision,
+                               aod::JetMcCollisions const&,
                                soa::Filtered<soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents, aod::ChargedMCDetectorLevelJetsMatchedToChargedMCParticleLevelJets>> const& mcdjets,
                                JetMCPTable const&,
                                aod::JetTracks const&,
@@ -226,8 +227,8 @@ struct JetTutorialTask {
     for (const auto& mcdjet : mcdjets) {
       for (auto& mcpjet : mcdjet.template matchedJetGeo_as<JetMCPTable>()) {
         registry.fill(HIST("h_matched_jets_pt"), mcpjet.pt(), mcdjet.pt(), collision.mcCollision().weight());
-        registry.fill(HIST("h_matched_jets_pt"), mcpjet.phi(), mcdjet.phi(), collision.mcCollision().weight());
-        registry.fill(HIST("h_matched_jets_pt"), mcpjet.eta(), mcdjet.eta(), collision.mcCollision().weight());
+        registry.fill(HIST("h_matched_jets_phi"), mcpjet.phi(), mcdjet.phi(), collision.mcCollision().weight());
+        registry.fill(HIST("h_matched_jets_eta"), mcpjet.eta(), mcdjet.eta(), collision.mcCollision().weight());
       }
     }
   }