@@ -44,13 +44,9 @@ struct Filter2Prong {
4444 O2_DEFINE_CONFIGURABLE (cfgYMax, float , -1 .0f , " Maximum candidate rapidity" )
4545 O2_DEFINE_CONFIGURABLE (cfgImPart1Mass, float , o2::constants::physics::MassKPlus, " Daughter particle 1 mass in GeV" )
4646 O2_DEFINE_CONFIGURABLE (cfgImPart2Mass, float , o2::constants::physics::MassKMinus, " Daughter particle 2 mass in GeV" )
47- // O2_DEFINE_CONFIGURABLE(cfgImPart1PID, int, o2::track::PID::Kaon, "PID of daughter particle 1 (O2 PID ID)")
48- // O2_DEFINE_CONFIGURABLE(cfgImPart2PID, int, o2::track::PID::Kaon, "PID of daughter particle 2 (O2 PID ID)")
49- O2_DEFINE_CONFIGURABLE (cfgMomDepPID, bool , 1 , " Use mommentum dependent PID for Phi meson" )
5047 O2_DEFINE_CONFIGURABLE (cfgImCutPt, float , 0 .2f , " Minimal pT for candidates" )
5148 O2_DEFINE_CONFIGURABLE (cfgImMinInvMass, float , 0 .95f , " Minimum invariant mass (GeV)" )
5249 O2_DEFINE_CONFIGURABLE (cfgImMaxInvMass, float , 1 .07f , " Maximum invariant mass (GeV)" )
53- // O2_DEFINE_CONFIGURABLE(cfgImSigmaFormula, std::string, "(z < 0.5 && x < 3.0) || (z >= 0.5 && x < 2.5 && y < 3.0)", "pT dependent daughter track sigma pass condition (x = TPC sigma, y = TOF sigma, z = pT)")
5450 O2_DEFINE_CONFIGURABLE (cfgDoPhi, bool , false , " Store phi information" )
5551 O2_DEFINE_CONFIGURABLE (cfgDoV0, bool , true , " Store V0s candidates" )
5652 O2_DEFINE_CONFIGURABLE (tpcNClsCrossedRowsTrackMin, float , 70 , " Minimum number of crossed rows in TPC" )
@@ -62,6 +58,7 @@ struct Filter2Prong {
6258 O2_DEFINE_CONFIGURABLE (ITSPIDPthreshold, float , 1.0 , " Momentum threshold for ITS PID (GeV/c) (only used if ITSPIDSelection is true)" )
6359 O2_DEFINE_CONFIGURABLE (ITSPIDNsigma, float , 3.0 , " PID nsigma for ITS" )
6460 O2_DEFINE_CONFIGURABLE (ConfDaughPIDCuts, float , 4.0 , " PID nsigma for V0s" )
61+ O2_DEFINE_CONFIGURABLE (cfgMomDepPID, bool , true , " Use momentum dependent PID cuts for V0s" )
6562 O2_DEFINE_CONFIGURABLE (massK0Min, float , 0.4 , " Minimum mass for K0" )
6663 O2_DEFINE_CONFIGURABLE (massK0Max, float , 0.6 , " Maximum mass for K0" )
6764 O2_DEFINE_CONFIGURABLE (massLambdaMin, float , 1.0 , " Minimum mass for lambda" )
@@ -108,12 +105,8 @@ struct Filter2Prong {
108105 using PIDTrack = soa::Join<aod::Tracks, aod::TracksExtra, aod::TrackSelection, aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr, aod::pidTOFPi, aod::pidTOFKa, aod::pidTOFPr, aod::pidTOFbeta>;
109106 using ResoV0s = aod::V0Datas;
110107
111- // std::unique_ptr<TFormula> sigmaFormula;
112-
113108 void init (InitContext&)
114109 {
115- // if (doprocessDataInvMass)
116- // sigmaFormula = std::make_unique<TFormula>("sigmaFormula", cfgImSigmaFormula.value.c_str());
117110 }
118111
119112 template <class HFCandidatesType >
@@ -394,6 +387,7 @@ struct Filter2Prong {
394387 if (removefaketrack && isFakeKaon (p1)) { // Check if the track is a fake kaon
395388 continue ;
396389 }
390+
397391 for (const auto & cftrack2 : cftracks) { // Loop over second track
398392 if (cftrack2.globalIndex () == cftrack1.globalIndex ()) // Skip if it's the same track as the first one
399393 continue ;
@@ -414,6 +408,7 @@ struct Filter2Prong {
414408 if (!selectionPair (p1, p2)) {
415409 continue ;
416410 }
411+
417412 ROOT::Math::PtEtaPhiMVector vec1 (p1.pt (), p1.eta (), p1.phi (), cfgImPart1Mass);
418413 ROOT::Math::PtEtaPhiMVector vec2 (p2.pt (), p2.eta (), p2.phi (), cfgImPart2Mass);
419414 ROOT::Math::PtEtaPhiMVector s = vec1 + vec2;
0 commit comments