|
35 | 35 |
|
36 | 36 | #include "CCDB/BasicCCDBManager.h" |
37 | 37 | #include "CCDB/CcdbApi.h" |
38 | | -#include "MathUtils/BetheBlochAleph.h" |
39 | 38 | #include "Framework/ASoA.h" |
40 | 39 | #include "Framework/ASoAHelpers.h" |
41 | 40 | #include "Framework/AnalysisDataModel.h" |
|
45 | 44 | #include "Framework/Logger.h" |
46 | 45 | #include "Framework/RunningWorkflowInfo.h" |
47 | 46 | #include "Framework/runDataProcessing.h" |
| 47 | +#include "MathUtils/BetheBlochAleph.h" |
48 | 48 | #include "ReconstructionDataFormats/DCA.h" |
49 | 49 | #include "ReconstructionDataFormats/PID.h" |
50 | 50 | #include "ReconstructionDataFormats/Track.h" |
@@ -3704,42 +3704,42 @@ struct AntinucleiInJets { |
3704 | 3704 | registryMC.fill(HIST("genEventsCoalescenceCorr"), 2.5); |
3705 | 3705 |
|
3706 | 3706 | // Build deuterons |
3707 | | - for (size_t iP = 0; iP < protonCandidates.size(); ++iP) { |
3708 | | - if (protonCandidates[iP].used) |
3709 | | - continue; |
3710 | | - |
3711 | | - for (size_t iN = 0; iN < neutronCandidates.size(); ++iN) { |
3712 | | - if (neutronCandidates[iN].used) |
3713 | | - continue; |
3714 | | - |
3715 | | - // Physics consistency check |
3716 | | - if (protonCandidates[iP].pdgCode * neutronCandidates[iN].pdgCode < 0) |
3717 | | - continue; |
3718 | | - |
3719 | | - if (passDeuteronCoalescence(protonCandidates[iP], neutronCandidates[iN], coalescenceMomentum, mRand)) { |
3720 | | - neutronCandidates[iN].used = true; |
3721 | | - protonCandidates[iP].used = true; |
3722 | | - |
3723 | | - int sign = (protonCandidates[iP].pdgCode > 0) ? +1 : -1; |
3724 | | - int deuteronPdg = sign * o2::constants::physics::Pdg::kDeuteron; |
3725 | | - |
3726 | | - double pxDeut = protonCandidates[iP].px + neutronCandidates[iN].px; |
3727 | | - double pyDeut = protonCandidates[iP].py + neutronCandidates[iN].py; |
3728 | | - double pzDeut = protonCandidates[iP].pz + neutronCandidates[iN].pz; |
3729 | | - double energyDeut = std::sqrt(pxDeut * pxDeut + pyDeut * pyDeut + pzDeut * pzDeut + massDeut * massDeut); |
3730 | | - LorentzVector pd(pxDeut, pyDeut, pzDeut, energyDeut); |
3731 | | - if (pd.Eta() >= minEta && pd.Eta() <= maxEta && (0.5 * pd.Pt()) >= MinPtPerNucleon) { |
3732 | | - // Store Deuteron |
3733 | | - finalDeuterons.push_back({pxDeut, pyDeut, pzDeut, deuteronPdg, protonCandidates[iP].mcIndex, false}); |
3734 | | - } |
3735 | | - |
3736 | | - break; |
3737 | | - } |
| 3707 | + for (size_t iP = 0; iP < protonCandidates.size(); ++iP) { |
| 3708 | + if (protonCandidates[iP].used) |
| 3709 | + continue; |
| 3710 | + |
| 3711 | + for (size_t iN = 0; iN < neutronCandidates.size(); ++iN) { |
| 3712 | + if (neutronCandidates[iN].used) |
| 3713 | + continue; |
| 3714 | + |
| 3715 | + // Physics consistency check |
| 3716 | + if (protonCandidates[iP].pdgCode * neutronCandidates[iN].pdgCode < 0) |
| 3717 | + continue; |
| 3718 | + |
| 3719 | + if (passDeuteronCoalescence(protonCandidates[iP], neutronCandidates[iN], coalescenceMomentum, mRand)) { |
| 3720 | + neutronCandidates[iN].used = true; |
| 3721 | + protonCandidates[iP].used = true; |
| 3722 | + |
| 3723 | + int sign = (protonCandidates[iP].pdgCode > 0) ? +1 : -1; |
| 3724 | + int deuteronPdg = sign * o2::constants::physics::Pdg::kDeuteron; |
| 3725 | + |
| 3726 | + double pxDeut = protonCandidates[iP].px + neutronCandidates[iN].px; |
| 3727 | + double pyDeut = protonCandidates[iP].py + neutronCandidates[iN].py; |
| 3728 | + double pzDeut = protonCandidates[iP].pz + neutronCandidates[iN].pz; |
| 3729 | + double energyDeut = std::sqrt(pxDeut * pxDeut + pyDeut * pyDeut + pzDeut * pzDeut + massDeut * massDeut); |
| 3730 | + LorentzVector pd(pxDeut, pyDeut, pzDeut, energyDeut); |
| 3731 | + if (pd.Eta() >= minEta && pd.Eta() <= maxEta && (0.5 * pd.Pt()) >= MinPtPerNucleon) { |
| 3732 | + // Store Deuteron |
| 3733 | + finalDeuterons.push_back({pxDeut, pyDeut, pzDeut, deuteronPdg, protonCandidates[iP].mcIndex, false}); |
3738 | 3734 | } |
| 3735 | + |
| 3736 | + break; |
| 3737 | + } |
3739 | 3738 | } |
3740 | | - |
3741 | | - // Add unused protons to final vectors |
3742 | | - for (const auto& proton : protonCandidates) { |
| 3739 | + } |
| 3740 | + |
| 3741 | + // Add unused protons to final vectors |
| 3742 | + for (const auto& proton : protonCandidates) { |
3743 | 3743 | if (!proton.used) { |
3744 | 3744 | finalProtons.push_back(proton); |
3745 | 3745 | } |
|
0 commit comments