Skip to content

Commit 5bf3e44

Browse files
YazhenLinYazhenLin
authored
Merge branch 'AliceO2Group:master' into master
2 parents 58eff73 + 2e678d5 commit 5bf3e44

File tree

11 files changed

+877
-329
lines changed

11 files changed

+877
-329
lines changed

‎ALICE3/TableProducer/alice3-decayfinder.cxx‎

Lines changed: 9 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -66,7 +66,8 @@ using std::array;
6666
// #define bitcheck(var, nbit) ((var) & (static_cast<uint32_t>(1) << (nbit)))
6767

6868
// For MC association in pre-selection
69-
using Alice3TracksWPid = soa::Join<aod::Tracks, aod::TracksCov, aod::Alice3DecayMaps, aod::McTrackLabels, aod::TracksDCA, aod::UpgradeTrkPids, aod::UpgradeTofs, aod::UpgradeRichs>;
69+
using Alice3TracksWPid = soa::Join<aod::Tracks, aod::TracksCov, aod::Alice3DecayMaps, aod::McTrackLabels, aod::TracksDCA, aod::UpgradeTofs, aod::UpgradeRichs>;
70+
using Alice3TracksWTrkPid = soa::Join<Alice3TracksWPid, aod::UpgradeTrkPids>;
7071

7172
struct alice3decayFinder {
7273
SliceCache cache;
@@ -185,18 +186,17 @@ struct alice3decayFinder {
185186
((aod::a3DecayMap::decayMap & trackSelectionKaMinusFromD) == trackSelectionKaMinusFromD) && aod::track::signed1Pt < 0.0f && nabs(aod::track::dcaXY) > kaFromD_dcaXYconstant + kaFromD_dcaXYpTdep* nabs(aod::track::signed1Pt);
186187

187188
// partitions for Lc baryons
188-
Partition<Alice3TracksWPid> tracksPiPlusFromLc =
189+
Partition<Alice3TracksWTrkPid> tracksPiPlusFromLc =
189190
((aod::a3DecayMap::decayMap & trackSelectionPiPlusFromLc) == trackSelectionPiPlusFromLc) && aod::track::signed1Pt > 0.0f && nabs(aod::track::dcaXY) > piFromLc_dcaXYconstant + piFromLc_dcaXYpTdep* nabs(aod::track::signed1Pt);
190-
Partition<Alice3TracksWPid> tracksKaPlusFromLc =
191+
Partition<Alice3TracksWTrkPid> tracksKaPlusFromLc =
191192
((aod::a3DecayMap::decayMap & trackSelectionKaPlusFromLc) == trackSelectionKaPlusFromLc) && aod::track::signed1Pt > 0.0f && nabs(aod::track::dcaXY) > kaFromLc_dcaXYconstant + kaFromLc_dcaXYpTdep* nabs(aod::track::signed1Pt);
192-
Partition<Alice3TracksWPid> tracksPrPlusFromLc =
193+
Partition<Alice3TracksWTrkPid> tracksPrPlusFromLc =
193194
((aod::a3DecayMap::decayMap & trackSelectionPrPlusFromLc) == trackSelectionPrPlusFromLc) && aod::track::signed1Pt > 0.0f && nabs(aod::track::dcaXY) > prFromLc_dcaXYconstant + prFromLc_dcaXYpTdep* nabs(aod::track::signed1Pt);
194-
// partitions for Lc baryons
195-
Partition<Alice3TracksWPid> tracksPiMinusFromLc =
195+
Partition<Alice3TracksWTrkPid> tracksPiMinusFromLc =
196196
((aod::a3DecayMap::decayMap & trackSelectionPiMinusFromLc) == trackSelectionPiMinusFromLc) && aod::track::signed1Pt < 0.0f && nabs(aod::track::dcaXY) > piFromLc_dcaXYconstant + piFromLc_dcaXYpTdep* nabs(aod::track::signed1Pt);
197-
Partition<Alice3TracksWPid> tracksKaMinusFromLc =
197+
Partition<Alice3TracksWTrkPid> tracksKaMinusFromLc =
198198
((aod::a3DecayMap::decayMap & trackSelectionKaMinusFromLc) == trackSelectionKaMinusFromLc) && aod::track::signed1Pt < 0.0f && nabs(aod::track::dcaXY) > kaFromLc_dcaXYconstant + kaFromLc_dcaXYpTdep* nabs(aod::track::signed1Pt);
199-
Partition<Alice3TracksWPid> tracksPrMinusFromLc =
199+
Partition<Alice3TracksWTrkPid> tracksPrMinusFromLc =
200200
((aod::a3DecayMap::decayMap & trackSelectionPrMinusFromLc) == trackSelectionPrMinusFromLc) && aod::track::signed1Pt < 0.0f && nabs(aod::track::dcaXY) > prFromLc_dcaXYconstant + prFromLc_dcaXYpTdep* nabs(aod::track::signed1Pt);
201201

202202
// Helper struct to pass candidate information
@@ -1068,7 +1068,7 @@ struct alice3decayFinder {
10681068

10691069
void processFindLc(aod::Collision const& collision,
10701070
aod::McParticles const& mcParticles,
1071-
Alice3TracksWPid const& tracks)
1071+
Alice3TracksWTrkPid const& tracks)
10721072
{
10731073
LOG(debug) << "Processing Lc candidates for collision " << collision.globalIndex() << " with " << tracks.size() << " tracks";
10741074
for (auto const& track : tracks) {

‎PWGCF/FemtoUniverse/Core/FemtoUniverseDetaDphiStar.h‎

Lines changed: 70 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -187,8 +187,9 @@ class FemtoUniverseDetaDphiStar
187187
}
188188
}
189189

190-
template <typename t1>
191-
void init_kT(HistogramRegistry* registry, t1& ktbins, std::vector<float> ldeltaphistarcutmin, std::vector<float> ldeltaphistarcutmax, std::vector<float> ldeltaetacutmin, std::vector<float> ldeltaetacutmax)
190+
template <typename t1, typename t2>
191+
// void init_kT(HistogramRegistry* registry, t1& ktbins, std::vector<float> ldeltaphistarcutmin, std::vector<float> ldeltaphistarcutmax, std::vector<float> ldeltaetacutmin, std::vector<float> ldeltaetacutmax, std::vector<float> ldeltaphistarcutFractionmin, std::vector<float> ldeltaphistarcutFractionmax, std::vector<float> ldeltaetacutFractionmin, std::vector<float> ldeltaetacutFractionmax)
192+
void init_kT(HistogramRegistry* registry, t1& ktbins, t2& ldeltaphistarcutmin, t2& ldeltaphistarcutmax, t2& ldeltaetacutmin, t2& ldeltaetacutmax, t2& ldeltaphistarcutFractionmin, t2& ldeltaphistarcutFractionmax, t2& ldeltaetacutFractionmin, t2& ldeltaetacutFractionmax)
192193
{
193194
mHistogramRegistry = registry;
194195
ktBins = ktbins;
@@ -198,6 +199,11 @@ class FemtoUniverseDetaDphiStar
198199
cutDeltaEtaMaxVector = ldeltaetacutmax;
199200
cutDeltaEtaMinVector = ldeltaetacutmin;
200201

202+
cutDeltaPhiStarFractionMaxVector = ldeltaphistarcutFractionmax;
203+
cutDeltaPhiStarFractionMinVector = ldeltaphistarcutFractionmin;
204+
cutDeltaEtaFractionMaxVector = ldeltaetacutFractionmax;
205+
cutDeltaEtaFractionMinVector = ldeltaetacutFractionmin;
206+
201207
if constexpr (kPartOneType == o2::aod::femtouniverseparticle::ParticleType::kTrack && kPartTwoType == o2::aod::femtouniverseparticle::ParticleType::kTrack) {
202208
std::string dirName = static_cast<std::string>(DirNames[0]);
203209
for (int j = 1; j < static_cast<int>(ktBins.size() - 1); j++) {
@@ -208,6 +214,8 @@ class FemtoUniverseDetaDphiStar
208214
histdetadphimixedbeforekT[j] = mHistogramRegistry->add<TH2>((dirName + histFolderkT + "detadphidetadphiBeforeMixed").c_str(), "; #Delta #eta; #Delta #phi", kTH2F, {{100, -0.15, 0.15}, {100, -0.15, 0.15}});
209215
histdetadphisameafterkT[j] = mHistogramRegistry->add<TH2>((dirName + histFolderkT + "detadphidetadphiAfterSame").c_str(), "; #Delta #eta; #Delta #phi", kTH2F, {{100, -0.15, 0.15}, {100, -0.15, 0.15}});
210216
histdetadphimixedafterkT[j] = mHistogramRegistry->add<TH2>((dirName + histFolderkT + "detadphidetadphiAfterMixed").c_str(), "; #Delta #eta; #Delta #phi", kTH2F, {{100, -0.15, 0.15}, {100, -0.15, 0.15}});
217+
histdetadphisameafterFractionkT[j] = mHistogramRegistry->add<TH2>((dirName + histFolderkT + "detadphidetadphiAfterFractionSame").c_str(), "; #Delta #eta; #Delta #phi", kTH2F, {{100, -0.15, 0.15}, {100, -0.15, 0.15}});
218+
histdetadphimixedafterFractionkT[j] = mHistogramRegistry->add<TH2>((dirName + histFolderkT + "detadphidetadphiAfterFractionMIxed").c_str(), "; #Delta #eta; #Delta #phi", kTH2F, {{100, -0.15, 0.15}, {100, -0.15, 0.15}});
211219
}
212220
}
213221
}
@@ -713,6 +721,55 @@ class FemtoUniverseDetaDphiStar
713721
}
714722
}
715723

724+
/// Check if pair is close or not
725+
template <typename Part>
726+
bool isClosePairFractionkT(Part const& part1, Part const& part2, uint8_t ChosenEventType, float ktval, float lmagfield, float DistMax, float FracMax)
727+
{
728+
/// Track-Track combination
729+
// check if provided particles are in agreement with the class instantiation
730+
if (part1.partType() != o2::aod::femtouniverseparticle::ParticleType::kTrack || part2.partType() != o2::aod::femtouniverseparticle::ParticleType::kTrack) {
731+
LOG(fatal) << "FemtoUniverseDetaDphiStar: passed arguments don't agree with FemtoUniverseDetaDphiStar instantiation! Please provide kTrack,kTrack candidates.";
732+
return false;
733+
}
734+
735+
int ktbinval = 1;
736+
if (ktval >= ktBins[1] && ktval < ktBins[2]) {
737+
ktbinval = 1;
738+
} else if (ktval >= ktBins[2] && ktval < ktBins[3]) {
739+
ktbinval = 2;
740+
} else if (ktval >= ktBins[3] && ktval < ktBins[4]) {
741+
ktbinval = 3;
742+
} else if (ktval >= ktBins[4] && ktval < ktBins[5]) {
743+
ktbinval = 4;
744+
}
745+
746+
magfield = lmagfield;
747+
748+
auto deta = part1.eta() - part2.eta();
749+
auto dphiAvg = averagePhiStar(part1, part2, 0);
750+
auto distfrac = averagePhiStarFrac(part1, part2, DistMax);
751+
auto DeltaPhiStarMax = static_cast<float>(cutDeltaPhiStarFractionMaxVector[ktbinval]);
752+
auto DeltaPhiStarMin = static_cast<float>(cutDeltaPhiStarFractionMinVector[ktbinval]);
753+
auto DeltaEtaMax = static_cast<float>(cutDeltaEtaFractionMaxVector[ktbinval]);
754+
auto DeltaEtaMin = static_cast<float>(cutDeltaEtaFractionMinVector[ktbinval]);
755+
756+
double outerVal = std::pow(dphiAvg, 2) / std::pow(DeltaPhiStarMax, 2) + std::pow(deta, 2) / std::pow(DeltaEtaMax, 2);
757+
double innerVal = std::pow(dphiAvg, 2) / std::pow(DeltaPhiStarMin, 2) + std::pow(deta, 2) / std::pow(DeltaEtaMin, 2);
758+
759+
if ((innerVal >= 1.0) && (outerVal < 1.0) && (distfrac > FracMax)) {
760+
return true;
761+
} else {
762+
if (ChosenEventType == femto_universe_container::EventType::same) {
763+
histdetadphisameafterFractionkT[ktbinval]->Fill(deta, dphiAvg);
764+
} else if (ChosenEventType == femto_universe_container::EventType::mixed) {
765+
histdetadphimixedafterFractionkT[ktbinval]->Fill(deta, dphiAvg);
766+
} else {
767+
LOG(fatal) << "FemtoUniverseDetaDphiStar: passed arguments don't agree with FemtoUniverseDetaDphiStar's type of events! Please provide same or mixed.";
768+
}
769+
return false;
770+
}
771+
}
772+
716773
/// Check if pair is close or not
717774
template <typename Part>
718775
void ClosePairqLCMS(Part const& part1, Part const& part2, float lmagfield, uint8_t ChosenEventType, double qlcms) // add typename Parts and variable parts for adding MClabels
@@ -787,10 +844,15 @@ class FemtoUniverseDetaDphiStar
787844
float cutDeltaEtaMax;
788845
float cutDeltaEtaMin;
789846

790-
std::vector<float> cutDeltaPhiStarMaxVector;
791-
std::vector<float> cutDeltaPhiStarMinVector;
792-
std::vector<float> cutDeltaEtaMaxVector;
793-
std::vector<float> cutDeltaEtaMinVector;
847+
std::vector<double> cutDeltaPhiStarMaxVector;
848+
std::vector<double> cutDeltaPhiStarMinVector;
849+
std::vector<double> cutDeltaEtaMaxVector;
850+
std::vector<double> cutDeltaEtaMinVector;
851+
852+
std::vector<double> cutDeltaPhiStarFractionMaxVector;
853+
std::vector<double> cutDeltaPhiStarFractionMinVector;
854+
std::vector<double> cutDeltaEtaFractionMaxVector;
855+
std::vector<double> cutDeltaEtaFractionMinVector;
794856

795857
float magfield;
796858
bool plotForEveryRadii = false;
@@ -805,6 +867,8 @@ class FemtoUniverseDetaDphiStar
805867
std::array<std::shared_ptr<TH2>, 4> histdetadphimixedbeforekT{};
806868
std::array<std::shared_ptr<TH2>, 4> histdetadphisameafterkT{};
807869
std::array<std::shared_ptr<TH2>, 4> histdetadphimixedafterkT{};
870+
std::array<std::shared_ptr<TH2>, 4> histdetadphisameafterFractionkT{};
871+
std::array<std::shared_ptr<TH2>, 4> histdetadphimixedafterFractionkT{};
808872

809873
std::array<std::array<std::shared_ptr<TH2>, 9>, 7> histdetadpiRadii{};
810874

‎PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrackSpherHarMultKtExtended.cxx‎

Lines changed: 26 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -89,14 +89,11 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
8989
Configurable<bool> confIsCPRatITS{"confIsCPRatITS", false, "Close Pair check at ITS"};
9090
Configurable<bool> confCPRPlotPerRadii{"confCPRPlotPerRadii", false, "Plot CPR per radii"};
9191
Configurable<float> confCPRdeltaPhiCutMax{"confCPRdeltaPhiCutMax", 0.0, "Delta Phi max cut for Close Pair Rejection"};
92-
Configurable<float> confCPRdeltaPhiCutMin{"confCPRdeltaPhiCutMin", 0., "Delta Phi min cut for Close Pair Rejection"};
93-
Configurable<float> confCPRdeltaEtaCutMax{"confCPRdeltaEtaCutMax", 0., "Delta Eta max cut for Close Pair Rejection"};
94-
Configurable<float> confCPRdeltaEtaCutMin{"confCPRdeltaEtaCutMin", 0., "Delta Eta min cut for Close Pair Rejection"};
95-
Configurable<std::vector<float>> confCPRdeltaPhiCutMaxVector{"confCPRdeltaPhiCutMaxVector", std::vector<float>{0.0, 0.0, 0.0, 0.0}, "Delta Phi max cut for Close Pair Rejection"};
96-
Configurable<std::vector<float>> confCPRdeltaPhiCutMinVector{"confCPRdeltaPhiCutMinVector", std::vector<float>{0.0, 0.0, 0.0, 0.0}, "Delta Phi min cut for Close Pair Rejection"};
97-
Configurable<std::vector<float>> confCPRdeltaEtaCutMaxVector{"confCPRdeltaEtaCutMaxVector", std::vector<float>{0.0, 0.0, 0.0, 0.0}, "Delta Eta max cut for Close Pair Rejection"};
98-
Configurable<std::vector<float>> confCPRdeltaEtaCutMinVector{"confCPRdeltaEtaCutMinVector", std::vector<float>{0.0, 0.0, 0.0, 0.0}, "Delta Eta min cut for Close Pair Rejection"};
92+
Configurable<float> confCPRdeltaPhiCutMin{"confCPRdeltaPhiCutMin", 0.0, "Delta Phi min cut for Close Pair Rejection"};
93+
Configurable<float> confCPRdeltaEtaCutMax{"confCPRdeltaEtaCutMax", 0.0, "Delta Eta max cut for Close Pair Rejection"};
94+
Configurable<float> confCPRdeltaEtaCutMin{"confCPRdeltaEtaCutMin", 0.0, "Delta Eta min cut for Close Pair Rejection"};
9995
Configurable<bool> confIsCPRkT{"confIsCPRkT", true, "kT dependent deltaEta-deltaPhi cut for Close Pair Rejection"};
96+
Configurable<bool> confIsCPRFraction{"confIsCPRFraction", true, "Close Pair Rejection with fraction of merged points"};
10097
Configurable<float> confTrkDCAxyMax{"confTrkDCAxyMax", 0.2, "Max DCA in xy "};
10198
Configurable<float> confTrkDCAzMax{"confTrkDCAzMax", 0.2, "Max DCA in z "};
10299
Configurable<int> confTrkTPCcRowsMin{"confTrkTPCcRowsMin", 80, "Min of TPC crossed rows"};
@@ -120,6 +117,17 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
120117
Configurable<bool> confisIdenLCMS{"confisIdenLCMS", true, "Choosing identical or non-identical pairs in LCMS"};
121118
Configurable<bool> confIsWeight{"confIsWeight", true, "Fill quantum weight"};
122119
Configurable<bool> confisIdenPRF{"confisIdenPRF", false, "Choosing identical or non-identical pairs in PRF"};
120+
121+
ConfigurableAxis confCPRdeltaEtaCutMaxVec{"confCPRdeltaEtaCutMaxVec", {VARIABLE_WIDTH, 0.0f, 0.0f, 0.0f, 0.0f}, "Bins for max delta eta cut in CPR"};
122+
ConfigurableAxis confCPRdeltaEtaCutMinVec{"confCPRdeltaEtaCutMinVec", {VARIABLE_WIDTH, 0.0f, 0.0f, 0.0f, 0.0f}, "Bins for min delta eta cut in CPR"};
123+
ConfigurableAxis confCPRdeltaPhiCutMaxVec{"confCPRdeltaPhiCutMaxVec", {VARIABLE_WIDTH, 0.0f, 0.0f, 0.0f, 0.0f}, "Bins for max delta phi cut in CPR"};
124+
ConfigurableAxis confCPRdeltaPhiCutMinVec{"confCPRdeltaPhiCutMinVec", {VARIABLE_WIDTH, 0.0f, 0.0f, 0.0f, 0.0f}, "Bins for min delta phi cut in CPR"};
125+
126+
ConfigurableAxis confCPRdeltaEtaCutFractionMaxVec{"confCPRdeltaEtaCutFractionMaxVec", {VARIABLE_WIDTH, 0.0f, 0.0f, 0.0f, 0.0f}, "Bins for max delta eta cut in CPR Fraction"};
127+
ConfigurableAxis confCPRdeltaEtaCutFractionMinVec{"confCPRdeltaEtaCutFractionMinVec", {VARIABLE_WIDTH, 0.0f, 0.0f, 0.0f, 0.0f}, "Bins for min delta eta cut in CPR Fraction"};
128+
ConfigurableAxis confCPRdeltaPhiCutFractionMaxVec{"confCPRdeltaPhiCutFractionMaxVec", {VARIABLE_WIDTH, 0.0f, 0.0f, 0.0f, 0.0f}, "Bins for max delta phi cut in CPR Fraction"};
129+
ConfigurableAxis confCPRdeltaPhiCutFractionMinVec{"confCPRdeltaPhiCutFractionMinVec", {VARIABLE_WIDTH, 0.0f, 0.0f, 0.0f, 0.0f}, "Bins for min delta phi cut in CPR Fraction"};
130+
123131
} twotracksconfigs;
124132

125133
using FemtoFullParticles = soa::Join<aod::FDParticles, aod::FDExtParticles>;
@@ -489,7 +497,7 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
489497
pairCleaner.init(&qaRegistry);
490498
if (twotracksconfigs.confIsCPR.value) {
491499
pairCloseRejection.init(&resultRegistry, &qaRegistry, twotracksconfigs.confCPRdeltaPhiCutMin.value, twotracksconfigs.confCPRdeltaPhiCutMax.value, twotracksconfigs.confCPRdeltaEtaCutMin.value, twotracksconfigs.confCPRdeltaEtaCutMax.value, twotracksconfigs.confCPRChosenRadii.value, twotracksconfigs.confCPRPlotPerRadii.value);
492-
pairCloseRejection.init_kT(&resultRegistry, confKtKstarBins, twotracksconfigs.confCPRdeltaPhiCutMinVector, twotracksconfigs.confCPRdeltaPhiCutMaxVector, twotracksconfigs.confCPRdeltaEtaCutMinVector, twotracksconfigs.confCPRdeltaEtaCutMaxVector);
500+
pairCloseRejection.init_kT(&resultRegistry, confKtKstarBins, twotracksconfigs.confCPRdeltaPhiCutMinVec.value, twotracksconfigs.confCPRdeltaPhiCutMaxVec.value, twotracksconfigs.confCPRdeltaEtaCutMinVec.value, twotracksconfigs.confCPRdeltaEtaCutMaxVec.value, twotracksconfigs.confCPRdeltaPhiCutFractionMinVec.value, twotracksconfigs.confCPRdeltaPhiCutFractionMaxVec.value, twotracksconfigs.confCPRdeltaEtaCutFractionMinVec.value, twotracksconfigs.confCPRdeltaEtaCutFractionMaxVec.value);
493501
}
494502

495503
vPIDPartOne = trackonefilter.confPIDPartOne.value;
@@ -633,6 +641,11 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
633641
if (pairCloseRejection.isClosePairkT(p1, p2, femto_universe_container::EventType::same, kT, twotracksconfigs.confIsCircularCut, twotracksconfigs.confCPRDphiAvgOrDist, magFieldTesla, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax)) {
634642
continue;
635643
}
644+
if (twotracksconfigs.confIsCPRFraction) {
645+
if (pairCloseRejection.isClosePairFractionkT(p1, p2, femto_universe_container::EventType::same, kT, magFieldTesla, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax)) {
646+
continue;
647+
}
648+
}
636649
} else {
637650
if (pairCloseRejection.isClosePairFrac(part1, part2, magFieldTesla, femto_universe_container::EventType::same, twotracksconfigs.confCPRDphiAvgOrDist, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax, twotracksconfigs.confIsCircularCut)) {
638651
continue;
@@ -996,6 +1009,11 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
9961009
if (pairCloseRejection.isClosePairkT(p1, p2, femto_universe_container::EventType::mixed, kT, twotracksconfigs.confIsCircularCut, twotracksconfigs.confCPRDphiAvgOrDist, magFieldTesla, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax)) {
9971010
continue;
9981011
}
1012+
if (twotracksconfigs.confIsCPRFraction) {
1013+
if (pairCloseRejection.isClosePairFractionkT(p1, p2, femto_universe_container::EventType::mixed, kT, magFieldTesla, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax)) {
1014+
continue;
1015+
}
1016+
}
9991017
} else {
10001018
if (pairCloseRejection.isClosePairFrac(part1, part2, magFieldTesla, femto_universe_container::EventType::mixed, twotracksconfigs.confCPRDphiAvgOrDist, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax, twotracksconfigs.confIsCircularCut)) {
10011019
continue;

0 commit comments

Comments
 (0)