Skip to content

[PWGEM] Adding histograms for accepted pi0 and eta into gamma gamma decay cha…#15465

Draft
amaringarcia wants to merge 1 commit intoAliceO2Group:masterfrom
amaringarcia:master
Draft

[PWGEM] Adding histograms for accepted pi0 and eta into gamma gamma decay cha…#15465
amaringarcia wants to merge 1 commit intoAliceO2Group:masterfrom
amaringarcia:master

Conversation

@amaringarcia
Copy link
Contributor

…nnel

BinnedAcceptedPi0/Eta are not stored

@github-actions github-actions bot added pwgem photon Photon and neutral meson PAG label labels Mar 20, 2026
@github-actions github-actions bot changed the title Adding histograms for accepted pi0 and eta into gamma gamma decay cha… [PWGEM] Adding histograms for accepted pi0 and eta into gamma gamma decay cha… Mar 20, 2026
auto lDaughter0 = lDaughters.begin();
if (lNDaughters == 2) {
auto lDaughter1 = lDaughters.iteratorAt(1);
if ( std::fabs(lDaughter0.eta())< cfg_max_eta_photon && std::fabs(lDaughter1.eta()) < cfg_max_eta_photon){
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Maybe we can add a proper function to make the acceptance checks to allow the check for the different photon detection methods? So we can have comparable acceptance between Run2 and Run3

Comment on lines +180 to +182
auto lDaughter0 = lDaughters.begin();
if (lNDaughters == 2) {
auto lDaughter1 = lDaughters.iteratorAt(1);
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Spawning iterators all the time in the loop is very expensive. One should create two iterators once before the loop and use setCursor(index) instead to move the iterator to the correct position.

if ((mcParticle.isPhysicalPrimary() || mcParticle.producedByGenerator()) && std::fabs(mcParticle.y()) < 0.9f && mcParticle.pt() < 20.f) {
auto binNumber = hBinFinder->FindBin(mcParticle.pt(), std::fabs(mcParticle.y())); // caution: pack
size_t lNDaughters = 0;
size_t mesonAccepted = 0;
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is currently used as a bool, please check if this is the intended use. If so, might be best to rename it to isMesonAccepted and make it a bool for clarity.

@vkucera vkucera marked this pull request as draft March 20, 2026 12:20
@vkucera
Copy link
Collaborator

vkucera commented Mar 20, 2026

@amaringarcia Please check your changes for formatting before you open a PR for review.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

photon Photon and neutral meson PAG label pwgem

Development

Successfully merging this pull request may close these issues.

3 participants