[PWGJE] New code for hadronjet correlations of Jan Janis#16387
Conversation
|
O2 linter results: ❌ 1 errors, |
|
Error while checking build/O2Physics/staging for e4c372b at 2026-05-26 10:49: Full log here. |
[PWGJE] Please consider the following formatting changes to AliceO2Group#16387
|
Error while checking build/O2Physics/o2 for e4c372b at 2026-05-26 10:52: Full log here. |
|
Please test the code properly before opening the PR for review. |
|
Please restore the commit history of this PR and reopen it. |
[PWGJE] Please consider the following formatting changes to AliceO2Group#16387
[PWGJE] Please consider the following formatting changes to AliceO2Group#16387
| #include "Common/DataModel/Multiplicity.h" | ||
|
|
||
| #include <CommonConstants/MathConstants.h> | ||
| #include <Framework/ASoA.h> |
There was a problem hiding this comment.
Could you please check if all these includes are needed? They were recently cleaned up so you could discuss with vit if these are an updated list
|
|
||
| using FilteredCollPartLevelMB = soa::Filtered<soa::Join<aod::JetMcCollisions, aod::BkgChargedMcRhos>>::iterator; | ||
|
|
||
| using CollPartLevelMBEfficiency = aod::JetMcCollisions::iterator; |
There was a problem hiding this comment.
for places where these are only used once you can directly use the table instead of this alias
| using namespace o2::framework::expressions; | ||
|
|
||
| using FilteredColl = soa::Filtered<soa::Join<aod::JetCollisions, aod::BkgChargedRhos>>::iterator; | ||
| using FilteredEventMultiplicity = soa::Filtered<soa::Join<aod::JetCollisions, aod::ZDCMults>>::iterator; |
There was a problem hiding this comment.
the ZDCMult table would not allow you to use the JE derived data. is the information needed?
|
|
||
| using FilteredTracksMatch = soa::Filtered<soa::Join<aod::JetTracks, aod::McTrackLabels>>; | ||
|
|
||
| using ParticlesFullTable = soa::Join<aod::JetParticles, aod::JMcParticlePIs>; |
There was a problem hiding this comment.
you should remove the PIs tables evrywhere since you dont actually access the original tables
|
|
||
| auto mcpjetsPerMCCollision = mcpjets.sliceBy(partJetsPerCollision, collision.mcCollisionId()); | ||
|
|
||
| fillMatchedHistograms(tracks, mcpjetsPerMCCollision, mcdjets); |
There was a problem hiding this comment.
does this not need a weight?
| float detJetPt = jetMatched.pt(); | ||
|
|
||
| spectra.fill(HIST("hJetPt_resolution"), (partJetPt - detJetPt) / partJetPt, partJetPt, weight); | ||
| spectra.fill(HIST("hJetPhi_resolution"), partJet.phi() - jetMatched.phi(), partJetPt, weight); |
There was a problem hiding this comment.
you are not wrapping phi. You should use the constrainAngle function in recodecay
| continue; | ||
|
|
||
| auto collTracks = tracks.sliceBy(tracksPerJCollision, collision.globalIndex()); | ||
| spectra.fill(HIST("hNTracksPerCollision"), collTracks.size()); |
There was a problem hiding this comment.
does this need a weight?
| while (bCloseJet) { | ||
| iii++; | ||
| if (iii > maxRandomConeIterations) { | ||
| break; |
There was a problem hiding this comment.
if this loop breaks previous values of randJetPhi are still used in line 495
| } | ||
|
|
||
| randJetPhi = rand.Uniform(-constants::math::PI, constants::math::PI); | ||
| randJetEta = rand.Uniform(-0.5, 0.5); |
There was a problem hiding this comment.
should it be a configurable?
|
|
||
| // Z coordinate of collision vertex | ||
| spectra.add("hVertexZ_NoCut", "z vertex of collisions w/o cut", kTH1F, {{100, -20., 20., "#it{z}_{vertex}"}}); | ||
| spectra.add("hVertexZ_Cut", "z vertex of collisions w. cut", kTH1F, {{100, -20., 20., "#it{z}_{vertex}"}}); |
There was a problem hiding this comment.
some histograms like this are never filled
This is a code of Jan Janis who is undergrad student of Filip Krizek. Jan will use this code to study hadron jet correlations.