[PWGHF] Add hadron and electron PID cuts in HFE–hadron correlation#14767
[PWGHF] Add hadron and electron PID cuts in HFE–hadron correlation#14767rashigupt wants to merge 18 commits intoAliceO2Group:masterfrom
Conversation
Two separate columns are included in the analysis table: one corresponding to the LSmass selection and the other to the ULmass selection, and the common mass column has been removed.
We define the mass column in the analysis by splitting it into two separate columns, LSmass and ULmass.
…ass. We use the mass column correctly by separating it into LSmass and ULmass.
Add PID cut for find the efficiency
Add PID cut for find Hadron efficiency
|
Please do not write a full story of your motivations in the PR title. You can give all the details in the PR description. The PR title is supposed to be a concise summary of what changed and where. |
|
Thank you for the suggestion. I have updated the PR title to make it concise and focused on the actual change. |
The title has not changed. |
| // Electron-hadron Correlation | ||
| template <typename TracksType, typename ElectronType, typename CollisionType, typename BcType> | ||
| void fillCorrelation(CollisionType const& collision, ElectronType const& electrons, TracksType const& tracks, BcType const&) | ||
| template <bool isMc, typename TracksType, typename ElectronType, typename CollisionType, typename BcType, typename McParticles> |
There was a problem hiding this comment.
| template <bool isMc, typename TracksType, typename ElectronType, typename CollisionType, typename BcType, typename McParticles> | |
| template <bool IsMc, typename TracksType, typename ElectronType, typename CollisionType, typename BcType, typename McParticles> |
McParticles is an existing table type. Do not use it as a template parameter.
There was a problem hiding this comment.
Thankyou for the suggestion
| if (std::abs(mcParticle.pdgCode()) != kElectron) { | ||
| registry.fill(HIST("hMCRecptEleHadron"), hTrack.pt()); | ||
| } | ||
| if (std::abs(mcParticle.pdgCode()) != kMuonMinus) { | ||
| registry.fill(HIST("hMCRecptMuonHadron"), hTrack.pt()); | ||
| } | ||
| if (std::abs(mcParticle.pdgCode()) != kPiPlus) { | ||
| registry.fill(HIST("hMCRecptPionHadron"), hTrack.pt()); | ||
| } | ||
| if (std::abs(mcParticle.pdgCode()) != kKPlus) { | ||
| registry.fill(HIST("hMCRecptKaonHadron"), hTrack.pt()); | ||
| } | ||
| if (std::abs(mcParticle.pdgCode()) != kProton) { | ||
| registry.fill(HIST("hMCRecptProtonHadron"), hTrack.pt()); | ||
| } |
There was a problem hiding this comment.
These conditions are mutually exclusive. Use else if or a switch statement.
There was a problem hiding this comment.
Thankyou for the guidance
| if (std::abs(particleMc.pdgCode()) != kElectron) { | ||
| registry.fill(HIST("hMCgenptEleHadron"), particleMc.pt()); | ||
| } | ||
| if (std::abs(particleMc.pdgCode()) != kMuonMinus) { | ||
| registry.fill(HIST("hMCgenptMuonHadron"), particleMc.pt()); | ||
| } | ||
| if (std::abs(particleMc.pdgCode()) != kPiPlus) { | ||
| registry.fill(HIST("hMCgenptPionHadron"), particleMc.pt()); | ||
| } | ||
| if (std::abs(particleMc.pdgCode()) != kKPlus) { |
There was a problem hiding this comment.
Thank you for the guidance.
| if (std::abs(mcParticle.pdgCode()) != kElectron) { | ||
| registry.fill(HIST("hMCRecptEleHadron"), hTrack.pt()); | ||
| } |
There was a problem hiding this comment.
Why do you fill the histogram with particles which are not electrons?
There was a problem hiding this comment.
Apologies for the mistake. Thank you for pointing it out!
…ter PID for correlation
Electron and hadron efficiencies are evaluated by applying PID cuts at both the MC generation and reconstruction levels: