@@ -62,10 +62,10 @@ struct HfCandidateSelectorToOmegaKa {
6262 // Produces<aod::HfMlSelOmegaKaToOmegaKa> hfMlSelToOmegaKa;
6363
6464 // LF analysis selections
65- Configurable<double > radiusCascMin{" radiusCascMin" , 0.5 , " Min cascade radius" };
65+ Configurable<double > radiusCascMin{" radiusCascMin" , 0.1 , " Min cascade radius" };
6666 Configurable<double > radiusV0Min{" radiusV0Min" , 1.1 , " Min V0 radius" };
67- Configurable<double > cosPAV0Min{" cosPAV0Min" , 0.97 , " Min valueCosPA V0" };
68- Configurable<double > cosPACascMin{" cosPACascMin" , 0.97 , " Min value CosPA cascade" };
67+ Configurable<double > cosPAV0Min{" cosPAV0Min" , 0.95 , " Min valueCosPA V0" };
68+ Configurable<double > cosPACascMin{" cosPACascMin" , 0.95 , " Min value CosPA cascade" };
6969 Configurable<double > dcaCascDauMax{" dcaCascDauMax" , 1.0 , " Max DCA cascade daughters" };
7070 Configurable<double > dcaV0DauMax{" dcaV0DauMax" , 1.0 , " Max DCA V0 daughters" };
7171 Configurable<float > dcaBachToPvMin{" dcaBachToPvMin" , 0.04 , " DCA Bach To PV" };
@@ -92,8 +92,8 @@ struct HfCandidateSelectorToOmegaKa {
9292 Configurable<double > impactParameterXYCascMin{" impactParameterXYCascMin" , 0 ., " Min dcaxy cascade track to PV" };
9393 Configurable<double > impactParameterXYCascMax{" impactParameterXYCascMax" , 10 ., " Max dcaxy cascade track to PV" };
9494
95- Configurable<double > ptCandMin{" ptCandMin" , 0 ., " Lower bound of candidate pT" };
96- Configurable<double > ptCandMax{" ptCandMax" , 50 ., " Upper bound of candidate pT" };
95+ Configurable<double > ptCandMin{" ptCandMin" , 1 ., " Lower bound of candidate pT" };
96+ Configurable<double > ptCandMax{" ptCandMax" , 12 ., " Upper bound of candidate pT" };
9797
9898 Configurable<double > dcaCharmBaryonDauMax{" dcaCharmBaryonDauMax" , 2.0 , " Max DCA charm baryon daughters" };
9999
@@ -248,7 +248,6 @@ struct HfCandidateSelectorToOmegaKa {
248248 registry.add (" hSelMassCasc" , " hSelMassCasc;status;entries" , {HistType::kTH1D , {axisSel}});
249249 registry.add (" hSelMassCharmBaryon" , " hSelMassCharmBaryon;status;entries" , {HistType::kTH1D , {axisSel}});
250250 registry.add (" hSelDcaXYToPvKaFromCasc" , " hSelDcaXYToPvKaFromCasc;status;entries" , {HistType::kTH1D , {axisSel}});
251- registry.add (" hSelPtOmegaKa" , " hSelPtOmegaKa;status;entries" , {HistType::kTH1D , {axisSel}});
252251 registry.add (" hSelCompetingCasc" , " hSelCompetingCasc;status;entries" , {HistType::kTH1D , {axisSel}});
253252 registry.add (" hSelV0_Casc_OmegaKaldl" , " hSelV0_Casc_OmegaKaldl;status;entries" , {HistType::kTH1D , {axisSel}});
254253 registry.add (" hSelctauOmegaKa" , " hSelctauOmegaKa;status;entries" , {HistType::kTH1D , {axisSel}});
@@ -265,21 +264,18 @@ struct HfCandidateSelectorToOmegaKa {
265264 bool selectionTopol (const T1& hfCandOmegaKa)
266265 {
267266 auto candpT = hfCandOmegaKa.kfPtOmegaKa ();
267+ // check that the candidate pT is within the analysis range
268+ if (candpT < ptCandMin || candpT > ptCandMax) {
269+ return false ;
270+ }
268271 auto kaPtFromOmegaKa = hfCandOmegaKa.kfPtKaFromOmegaKa ();
269272 int const pTBin = findBin (binsPt, candpT);
270273 if (pTBin == -1 ) {
271274 return false ;
272275 }
273276
274277 // check that the candidate pT is within the analysis range
275- if (candpT <= ptCandMin || candpT >= ptCandMax) {
276- registry.fill (HIST (" hSelPtOmegaKa" ), 0 );
277- return false ;
278- }
279- registry.fill (HIST (" hSelPtOmegaKa" ), 1 );
280-
281- // check that the candidate pT is within the analysis range
282- if (kaPtFromOmegaKa < cuts->get (pTBin, " pT ka from OmegaKa" )) {
278+ if (kaPtFromOmegaKa < cuts->get (pTBin, " pT Ka from Omegac" )) {
283279 registry.fill (HIST (" hSelPtKaFromCharm" ), 0 );
284280 return false ;
285281 }
@@ -329,10 +325,13 @@ struct HfCandidateSelectorToOmegaKa {
329325 } else if (signDecay < 0 ) {
330326 registry.fill (HIST (" hSelSignDec" ), 0 ); // particle decay
331327 }
332-
333328 // pt-dependent selection
334329 if (!selectionTopol (candidate)) {
335330 resultSelections = false ;
331+ hfSelToOmegaKaKf (statusPidLambda, statusPidCascade, statusPidCharmBaryon, statusInvMassLambda, statusInvMassCascade, statusInvMassCharmBaryon, resultSelections, infoTpcStored, infoTofStored,
332+ trackKaFromCharm.tpcNSigmaKa (), trackKaFromCasc.tpcNSigmaKa (), trackPiFromLam.tpcNSigmaPi (), trackPrFromLam.tpcNSigmaPr (),
333+ trackKaFromCharm.tofNSigmaKa (), trackKaFromCasc.tofNSigmaKa (), trackPiFromLam.tofNSigmaPi (), trackPrFromLam.tofNSigmaPr ());
334+ continue ;
336335 }
337336
338337 // eta selection
0 commit comments