@@ -182,6 +182,7 @@ struct CorrelationTask {
182182 if (doprocessSame2Prong2Prong || doprocessSame2Prong2ProngML) {
183183 registry.add (" invMassTwoPart" , " 2D 2-prong invariant mass (GeV/c^2)" , {HistType::kTHnSparseF , {axisSpecMass, axisSpecMass, axisPtTrigger, axisPtAssoc, axisMultiplicity}});
184184 registry.add (" invMassTwoPartDPhi" , " 2D 2-prong invariant mass (GeV/c^2)" , {HistType::kTHnSparseF , {axisSpecMass, axisSpecMass, axisPtTrigger, axisPtAssoc, axisDeltaPhi}});
185+ registry.add (" invMassTwoPartDEta" , " 2D 2-prong invariant mass (GeV/c^2)" , {HistType::kTHnSparseF , {axisSpecMass, axisSpecMass, axisPtTrigger, axisPtAssoc, axisDeltaEta}});
185186 }
186187 }
187188 if (doprocessSameDerivedMultSet) {
@@ -394,6 +395,9 @@ struct CorrelationTask {
394395 continue ;
395396 registry.fill (HIST (" invMassTwoPart" ), track1.invMass (), track2.invMass (), track1.pt (), track2.pt (), multiplicity);
396397 registry.fill (HIST (" invMassTwoPartDPhi" ), track1.invMass (), track2.invMass (), track1.pt (), track2.pt (), TVector2::Phi_0_2pi (track1.phi () - track2.phi () + TMath::Pi () / 2.0 ) - TMath::Pi () / 2.0 );
398+ if (std::abs (track1.phi () - track2.phi ()) < constants::math::PI * 0.5 ) {
399+ registry.fill (HIST (" invMassTwoPartDEta" ), track1.invMass (), track2.invMass (), track1.pt (), track2.pt (), track1.eta () - track2.eta ());
400+ }
397401 }
398402 }
399403 }
@@ -637,7 +641,7 @@ struct CorrelationTask {
637641 continue ; // skip particles that do not match the decay mask
638642 }
639643 if (cfgV0RapidityMax > 0 ) {
640- auto [t, y] = getV0Rapidity (track1 );
644+ auto [t, y] = getV0Rapidity (track2 );
641645 if (t && std::abs (y) > cfgV0RapidityMax)
642646 continue ; // V0s are not allowed to be outside the rapidity range
643647 }
0 commit comments