Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

[PWGJE] Adding extra options for rho estimation #9664

Merged
merged 1 commit into from
Feb 3, 2025
Merged
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 8 additions & 5 deletions PWGJE/TableProducer/rhoEstimator.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -47,17 +47,19 @@
Configurable<float> trackPtMax{"trackPtMax", 1000.0, "maximum track pT"};
Configurable<float> trackEtaMin{"trackEtaMin", -0.9, "minimum track eta"};
Configurable<float> trackEtaMax{"trackEtaMax", 0.9, "maximum track eta"};
Configurable<float> trackPhiMin{"trackPhiMin", -999, "minimum track phi"};
Configurable<float> trackPhiMax{"trackPhiMax", 999, "maximum track phi"};
Configurable<float> trackPhiMin{"trackPhiMin", -99.0, "minimum track phi"};
Configurable<float> trackPhiMax{"trackPhiMax", 99.0, "maximum track phi"};
Configurable<double> trackingEfficiency{"trackingEfficiency", 1.0, "tracking efficiency applied to jet finding"};
Configurable<std::string> trackSelections{"trackSelections", "globalTracks", "set track selections"};

Configurable<std::string> particleSelections{"particleSelections", "PhysicalPrimary", "set particle selections"};

Configurable<int> jetAlgorithm{"jetAlgorithm", 0, "jet clustering algorithm. 0 = kT, 1 = C/A, 2 = Anti-kT"};
Configurable<int> jetRecombScheme{"jetRecombScheme", 0, "jet recombination scheme. 0 = E-scheme, 1 = pT-scheme, 2 = pT2-scheme"};
Configurable<float> bkgjetR{"bkgjetR", 0.2, "jet resolution parameter for determining background density"};
Configurable<float> bkgEtaMin{"bkgEtaMin", -0.9, "minimim pseudorapidity for determining background density"};
Configurable<float> bkgEtaMax{"bkgEtaMax", 0.9, "maximum pseudorapidity for determining background density"};
Configurable<float> bkgPhiMin{"bkgPhiMin", 0., "minimim phi 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.0, "minimim phi for determining background density"};
Configurable<float> bkgPhiMax{"bkgPhiMax", 99.0, "maximum phi for determining background density"};
Configurable<bool> doSparse{"doSparse", false, "perfom sparse estimation"};

Expand All @@ -74,10 +76,11 @@
trackSelection = jetderiveddatautilities::initialiseTrackSelection(static_cast<std::string>(trackSelections));
particleSelection = static_cast<std::string>(particleSelections);

bkgSub.setJetAlgorithmAndScheme(static_cast<fastjet::JetAlgorithm>(static_cast<int>(jetAlgorithm)), static_cast<fastjet::RecombinationScheme>(static_cast<int>(jetRecombScheme)));
bkgSub.setJetBkgR(bkgjetR);
bkgSub.setEtaMinMax(bkgEtaMin, bkgEtaMax);
if (bkgPhiMax > 98.0) {
bkgPhiMax_ = 2.0 * M_PI;

Check warning on line 83 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[external-pi]

Consider using the PI constant (and its multiples and fractions) defined in o2::constants::math.

Check warning on line 83 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pi-multiple-fraction]

Consider using multiples/fractions of PI defined in o2::constants::math.
}
bkgSub.setPhiMinMax(bkgPhiMin, bkgPhiMax_);
}
Expand Down Expand Up @@ -114,7 +117,7 @@
void processD0Collisions(aod::JetCollision const&, soa::Filtered<aod::JetTracks> const& tracks, aod::CandidatesD0Data const& candidates)
{
inputParticles.clear();
for (auto& candidate : candidates) {

Check warning on line 120 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
inputParticles.clear();
jetfindingutilities::analyseTracks(inputParticles, tracks, trackSelection, trackingEfficiency, std::optional{candidate});

Expand All @@ -127,7 +130,7 @@
void processD0McCollisions(aod::JetMcCollision const&, soa::Filtered<aod::JetParticles> const& particles, aod::CandidatesD0MCP const& candidates)
{
inputParticles.clear();
for (auto& candidate : candidates) {

Check warning on line 133 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
inputParticles.clear();
jetfindingutilities::analyseParticles<true>(inputParticles, particleSelection, 1, particles, pdgDatabase, std::optional{candidate});

Expand All @@ -140,7 +143,7 @@
void processLcCollisions(aod::JetCollision const&, soa::Filtered<aod::JetTracks> const& tracks, aod::CandidatesLcData const& candidates)
{
inputParticles.clear();
for (auto& candidate : candidates) {

Check warning on line 146 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
inputParticles.clear();
jetfindingutilities::analyseTracks(inputParticles, tracks, trackSelection, trackingEfficiency, std::optional{candidate});

Expand All @@ -153,7 +156,7 @@
void processLcMcCollisions(aod::JetMcCollision const&, soa::Filtered<aod::JetParticles> const& particles, aod::CandidatesLcMCP const& candidates)
{
inputParticles.clear();
for (auto& candidate : candidates) {

Check warning on line 159 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
inputParticles.clear();
jetfindingutilities::analyseParticles<true>(inputParticles, particleSelection, 1, particles, pdgDatabase, std::optional{candidate});

Expand All @@ -166,7 +169,7 @@
void processBplusCollisions(aod::JetCollision const&, soa::Filtered<aod::JetTracks> const& tracks, aod::CandidatesBplusData const& candidates)
{
inputParticles.clear();
for (auto& candidate : candidates) {

Check warning on line 172 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
inputParticles.clear();
jetfindingutilities::analyseTracks(inputParticles, tracks, trackSelection, trackingEfficiency, std::optional{candidate});

Expand All @@ -179,7 +182,7 @@
void processBplusMcCollisions(aod::JetMcCollision const&, soa::Filtered<aod::JetParticles> const& particles, aod::CandidatesBplusMCP const& candidates)
{
inputParticles.clear();
for (auto& candidate : candidates) {

Check warning on line 185 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
inputParticles.clear();
jetfindingutilities::analyseParticles<true>(inputParticles, particleSelection, 1, particles, pdgDatabase, std::optional{candidate});

Expand All @@ -192,7 +195,7 @@
void processDielectronCollisions(aod::JetCollision const&, soa::Filtered<aod::JetTracks> const& tracks, aod::CandidatesDielectronData const& candidates)
{
inputParticles.clear();
for (auto& candidate : candidates) {

Check warning on line 198 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
inputParticles.clear();
jetfindingutilities::analyseTracks(inputParticles, tracks, trackSelection, trackingEfficiency, std::optional{candidate});

Expand All @@ -205,7 +208,7 @@
void processDielectronMcCollisions(aod::JetMcCollision const&, soa::Filtered<aod::JetParticles> const& particles, aod::CandidatesDielectronMCP const& candidates)
{
inputParticles.clear();
for (auto& candidate : candidates) {

Check warning on line 211 in PWGJE/TableProducer/rhoEstimator.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
inputParticles.clear();
jetfindingutilities::analyseParticles<true>(inputParticles, particleSelection, 1, particles, pdgDatabase, std::optional{candidate});

Expand Down
Loading