@@ -260,11 +260,11 @@ struct HfTreeCreatorOmegacSt {
260260 Configurable<int > selectedChannel{" selectedChannel" , 0 , " Selected decay channel for Xic or Omegac: 0=all, 1=Xic->Xi+Pi, 2=Xic->Omega+K, 3=Omegac->Omega+Pi, 4=Omegac->Omega+K, 5=Omegac->Xi+Pi" };
261261 Configurable<float > maxMatchingChi2TrackedCascade{" maxMatchingChi2TrackedCascade" , 2000 ., " Max matching chi2 for tracked cascades" };
262262 Configurable<bool > recalculateMasses{" recalculateMasses" , true , " Recalculate Xi/Omega masses" };
263- Configurable<float > maxNSigmaBachelor{ " maxNSigmaBachelor " , 5 ., " Max Nsigma for bachelor of tracked cascade" };
264- Configurable<float > maxNSigmaV0Pr{ " maxNSigmaV0Pr " , 5 ., " Max Nsigma for proton from V0 from tracked cascade" };
265- Configurable<float > maxNSigmaV0Pi{ " maxNSigmaV0Pi " , 5 ., " Max Nsigma for pion from V0 from tracked cascade" };
266- Configurable<float > maxNSigmaPion{ " maxNSigmaPion " , 5 ., " Max Nsigma for pion to be paired with Omega" };
267- Configurable<float > maxNSigmaKaon{ " maxNSigmaKaon " , 5 ., " Max Nsigma for kaon to be paired with Omega" };
263+ Configurable<float > nSigmaBachelorMax{ " nSigmaBachelorMax " , 5 ., " Max Nsigma for bachelor of tracked cascade" };
264+ Configurable<float > nSigmaV0PrMax{ " nSigmaV0PrMax " , 5 ., " Max Nsigma for proton from V0 from tracked cascade" };
265+ Configurable<float > nSigmaV0PiMax{ " nSigmaV0PiMax " , 5 ., " Max Nsigma for pion from V0 from tracked cascade" };
266+ Configurable<float > nSigmaPionMax{ " nSigmaPionMax " , 5 ., " Max Nsigma for pion to be paired with Omega" };
267+ Configurable<float > nSigmaKaonMax{ " nSigmaKaonMax " , 5 ., " Max Nsigma for kaon to be paired with Omega" };
268268 Configurable<bool > bzOnly{" bzOnly" , true , " Use B_z instead of full field map" };
269269 Configurable<std::string> cfgTriggersOfInterest{" cfgTriggersOfInterest" , " fTrackedOmega,fHfCharmBarToXiBach" , " Triggers of interest, comma separated for Zorro" };
270270
@@ -279,28 +279,28 @@ struct HfTreeCreatorOmegacSt {
279279 Configurable<bool > useTofPid{" useTofPid" , false , " Require TOF PID together with TPC PID when true; otherwise use TPC only" };
280280
281281 // pT cut for CharmBaryon
282- Configurable<float > minPtCharmedBaryon{ " minPtCharmedBaryon " , 0 .f , " Minimum pT for the charmed baryon" };
283- Configurable<float > maxPtCharmedBaryon{ " maxPtCharmedBaryon " , 50 .f , " Maximum pT for the charmed baryon" };
282+ Configurable<float > ptCharmedBaryonMin{ " ptCharmedBaryonMin " , 0 .f , " Minimum pT for the charmed baryon" };
283+ Configurable<float > ptCharmedBaryonMax{ " ptCharmedBaryonMax " , 50 .f , " Maximum pT for the charmed baryon" };
284284
285285 // CPA cuts (defaults: >= 0.9)
286- Configurable<float > minCpaCharmedBaryon{ " minCpaCharmedBaryon " , 0 .9f , " Minimum CPA for charmed baryon" };
287- Configurable<float > minCpaXYCharmedBaryon{ " minCpaXYCharmedBaryon " , 0 .9f , " Minimum CPA XY for charmed baryon" }; // remove maybe???
288- Configurable<float > minCpaCasc{ " minCpaCasc " , 0 .9f , " Minimum CPA for cascade" };
289- Configurable<float > minCpaXYCasc{ " minCpaXYCasc " , 0 .9f , " Minimum CPA XY for cascade" };
286+ Configurable<float > cpaCharmedBaryonMin{ " cpaCharmedBaryonMin " , 0 .9f , " Minimum CPA for charmed baryon" };
287+ Configurable<float > cpaXYCharmedBaryonMin{ " cpaXYCharmedBaryonMin " , 0 .9f , " Minimum CPA XY for charmed baryon" }; // remove maybe???
288+ Configurable<float > cpaCascMin{ " cpaCascMin " , 0 .9f , " Minimum CPA for cascade" };
289+ Configurable<float > cpaXYCascMin{ " cpaXYCascMin " , 0 .9f , " Minimum CPA XY for cascade" };
290290
291291 // Chi2 cuts (defaults: <= 10)
292- Configurable<float > maxChi2TopCharmedBaryon{ " maxChi2TopCharmedBaryon " , 10 .f , " Maximum chi2 for top charmed baryon fit" };
293- Configurable<float > maxChi2TopCasc{ " maxChi2TopCasc " , 10 .f , " Maximum topologyChi2 of the tracked cascade" };
292+ Configurable<float > chi2TopCharmedBaryonMax{ " chi2TopCharmedBaryonMax " , 10 .f , " Maximum chi2 for top charmed baryon fit" };
293+ Configurable<float > chi2TopCascMax{ " chi2TopCascMax " , 10 .f , " Maximum topologyChi2 of the tracked cascade" };
294294
295295 // Impact parameter cuts (defaults: absolute value < 10)
296- Configurable<float > maxImpactCascY{ " maxImpactCascY " , 10 .f , " Max abs Max abs impactParameter of cascade(Y)" };
297- Configurable<float > maxImpactCascZ{ " maxImpactCascZ " , 10 .f , " Max abs DCA Z of cascade" };
298- Configurable<float > maxImpactCharmBachelorY{ " maxImpactCharmBachelorY " , 10 .f , " Max abs impactParameter of charm bachelor(Y)" };
299- Configurable<float > maxImpactCharmBachelorZ{ " maxImpactCharmBachelorZ " , 10 .f , " Max abs DCA Z of charm bachelor" };
300- Configurable<float > minPtCascade{ " minPtCascade " , 0 .f , " Minimum pT requirement for the cascade track" };
301- Configurable<float > minPtPionOrKaon{ " minPtPionOrKaon " , 0 .f , " Minimum pT requirement for the charm baryon daughter pion/kaon track" };
302- Configurable<float > minImpProduct{ " minImpProduct " , -1 .f , " Lower bound for DCA impact product" };
303- Configurable<float > maxImpProduct{ " maxImpProduct " , 0 .f , " Upper bound for DCA impact product" };
296+ Configurable<float > impactParameterCascYMax{ " impactParameterCascYMax " , 10 .f , " Max abs Max abs impactParameter of cascade(Y)" };
297+ Configurable<float > absDcaCascZMax{ " absDcaCascZMax " , 10 .f , " Max abs DCA Z of cascade" };
298+ Configurable<float > impactParameterCharmBachelorYMax{ " impactParameterCharmBachelorYMax " , 10 .f , " Max abs impactParameter of charm bachelor(Y)" };
299+ Configurable<float > absDcaCharmBachelorZMax{ " absDcaCharmBachelorZMax " , 10 .f , " Max abs DCA Z of charm bachelor" };
300+ Configurable<float > ptCascadeMin{ " ptCascadeMin " , 0 .f , " Minimum pT requirement for the cascade track" };
301+ Configurable<float > ptPionOrKaonMin{ " ptPionOrKaonMin " , 0 .f , " Minimum pT requirement for the charm baryon daughter pion/kaon track" };
302+ Configurable<float > impProductMin{ " impProductMin " , -1 .f , " Lower bound for DCA impact product" };
303+ Configurable<float > impProductMax{ " impProductMax " , 0 .f , " Upper bound for DCA impact product" };
304304
305305 SliceCache cache;
306306 Service<o2::ccdb::BasicCCDBManager> ccdb;
@@ -410,12 +410,12 @@ struct HfTreeCreatorOmegacSt {
410410 bool isMatched = false ;
411411 static constexpr std::size_t NDaughters{2u };
412412 enum ChannelType {
413- anyChannel = 0 ,
414- xicToXiPi = 1 ,
415- xicToOmegaKa = 2 ,
416- omegacToOmegaPi = 3 ,
417- omegacToOmegaKa = 4 ,
418- omegacToXiPi = 5
413+ AnyChannel = 0 ,
414+ XicToXiPi = 1 ,
415+ XicToOmegaKa = 2 ,
416+ OmegacToOmegaPi = 3 ,
417+ OmegacToOmegaKa = 4 ,
418+ OmegacToXiPi = 5
419419 };
420420
421421 void processMc (aod::McCollisions const &,
@@ -659,16 +659,16 @@ struct HfTreeCreatorOmegacSt {
659659 (std::abs (massXi - o2::constants::physics::MassXiMinus) < massWindowTrackedXi)) {
660660 // combine configurable TPC and TOF nsigma requirements for bachelor and V0 daughters
661661 // PID: decide based on `useTofPid` switch
662- const bool tpcBachelor = (std::abs (bachelor.tpcNSigmaKa ()) < maxNSigmaBachelor ) || (std::abs (bachelor.tpcNSigmaPi ()) < maxNSigmaBachelor );
662+ const bool tpcBachelor = (std::abs (bachelor.tpcNSigmaKa ()) < nSigmaBachelorMax ) || (std::abs (bachelor.tpcNSigmaPi ()) < nSigmaBachelorMax );
663663 const float tofBachelorPiAbs = std::abs (bachelor.tofNSigmaPi ());
664664 const float tofBachelorKaAbs = std::abs (bachelor.tofNSigmaKa ());
665665 const bool tofBachelorPiPass = (tofBachelorPiAbs > TofWoSignalRange[0 ] && tofBachelorPiAbs < TofWoSignalRange[1 ]) || (tofBachelorPiAbs < nSigmaAbsTofBachelorPiMax);
666666 const bool tofBachelorKaPass = (tofBachelorKaAbs > TofWoSignalRange[0 ] && tofBachelorKaAbs < TofWoSignalRange[1 ]) || (tofBachelorKaAbs < nSigmaAbsTofBachelorKaMax);
667667 const bool tofBachelorPass = tofBachelorPiPass || tofBachelorKaPass;
668668 const bool bachelorPass = useTofPid.value ? (tpcBachelor && tofBachelorPass) : tpcBachelor;
669669
670- const bool tpcV0Pr = (std::abs (v0TrackPr.tpcNSigmaPr ()) < maxNSigmaV0Pr );
671- const bool tpcV0Pi = (std::abs (v0TrackPi.tpcNSigmaPi ()) < maxNSigmaV0Pi );
670+ const bool tpcV0Pr = (std::abs (v0TrackPr.tpcNSigmaPr ()) < nSigmaV0PrMax );
671+ const bool tpcV0Pi = (std::abs (v0TrackPi.tpcNSigmaPi ()) < nSigmaV0PiMax );
672672 const float tofV0PrAbs = std::abs (v0TrackPr.tofNSigmaPr ());
673673 const float tofV0PiAbs = std::abs (v0TrackPi.tofNSigmaPi ());
674674 const bool tofV0PrPass = (tofV0PrAbs > TofWoSignalRange[0 ] && tofV0PrAbs < TofWoSignalRange[1 ]) || (tofV0PrAbs < nSigmaAbsTofV0PrMax);
@@ -713,7 +713,7 @@ struct HfTreeCreatorOmegacSt {
713713 (track.tpcChi2NCl () <= TpcChi2NclMax) &&
714714 (track.itsChi2NCl () <= ItsChi2NclMax)) {
715715 // TPC and TOF PID for pion/kaon from chambaryon
716- const bool passTPCpid = (std::abs (track.tpcNSigmaPi ()) < maxNSigmaPion ) || (std::abs (track.tpcNSigmaKa ()) < maxNSigmaKaon );
716+ const bool passTPCpid = (std::abs (track.tpcNSigmaPi ()) < nSigmaPionMax ) || (std::abs (track.tpcNSigmaKa ()) < nSigmaKaonMax );
717717 const float tofPiAbs = std::abs (track.tofNSigmaPi ());
718718 const float tofKaAbs = std::abs (track.tofNSigmaKa ());
719719 const bool tofPiPass = (tofPiAbs > TofWoSignalRange[0 ] && tofPiAbs < TofWoSignalRange[1 ]) || (tofPiAbs < nSigmaAbsTofTrackPiMax);
@@ -787,22 +787,22 @@ struct HfTreeCreatorOmegacSt {
787787 const bool anyChannelPass = massOmegacToOmegaKPass || massOmegacToOmegaPiPass || massXicToXiPiPass || massXicToOmegaKaPass || massOmegacToXiPiPass;
788788 bool passSelectedChannel = true ;
789789 switch (selectedChannel.value ) {
790- case anyChannel :
790+ case AnyChannel :
791791 passSelectedChannel = anyChannelPass;
792792 break ;
793- case xicToXiPi :
793+ case XicToXiPi :
794794 passSelectedChannel = massXicToXiPiPass;
795795 break ;
796- case xicToOmegaKa :
796+ case XicToOmegaKa :
797797 passSelectedChannel = massXicToOmegaKaPass;
798798 break ;
799- case omegacToOmegaPi :
799+ case OmegacToOmegaPi :
800800 passSelectedChannel = massOmegacToOmegaPiPass;
801801 break ;
802- case omegacToOmegaKa :
802+ case OmegacToOmegaKa :
803803 passSelectedChannel = massOmegacToOmegaKPass;
804804 break ;
805- case omegacToXiPi :
805+ case OmegacToXiPi :
806806 passSelectedChannel = massOmegacToXiPiPass;
807807 break ;
808808 default :
@@ -867,31 +867,31 @@ struct HfTreeCreatorOmegacSt {
867867 }
868868 }
869869 // apply configurable quality cuts (CPA, chi2, impact parameters)
870- if (cpaCharmedBaryon < minCpaCharmedBaryon )
870+ if (cpaCharmedBaryon < cpaCharmedBaryonMin )
871871 continue ;
872- if (cpaXYCharmedBaryon < minCpaXYCharmedBaryon )
872+ if (cpaXYCharmedBaryon < cpaXYCharmedBaryonMin )
873873 continue ;
874- if (cpaCasc < minCpaCasc )
874+ if (cpaCasc < cpaCascMin )
875875 continue ;
876- if (cpaXYCasc < minCpaXYCasc )
876+ if (cpaXYCasc < cpaXYCascMin )
877877 continue ;
878- if (chi2TopCharmedBaryon > maxChi2TopCharmedBaryon )
878+ if (chi2TopCharmedBaryon > chi2TopCharmedBaryonMax )
879879 continue ;
880- if (trackedCascade.topologyChi2 () > maxChi2TopCasc )
880+ if (trackedCascade.topologyChi2 () > chi2TopCascMax )
881881 continue ;
882- if (std::abs (impactParameterCasc.getY ()) > maxImpactCascY || std::abs (impactParameterCasc.getZ ()) > maxImpactCascZ )
882+ if (std::abs (impactParameterCasc.getY ()) > impactParameterCascYMax || std::abs (impactParameterCasc.getZ ()) > absDcaCascZMax )
883883 continue ;
884- if (std::abs (impactParameterPion.getY ()) > maxImpactCharmBachelorY || std::abs (impactParameterPion.getZ ()) > maxImpactCharmBachelorZ )
884+ if (std::abs (impactParameterPion.getY ()) > impactParameterCharmBachelorYMax || std::abs (impactParameterPion.getZ ()) > absDcaCharmBachelorZMax )
885885 continue ;
886- if (ptCascade < minPtCascade )
886+ if (ptCascade < ptCascadeMin )
887887 continue ;
888- if (ptPionOrKaon < minPtPionOrKaon )
888+ if (ptPionOrKaon < ptPionOrKaonMin )
889889 continue ;
890- if (ptCharmedBaryon < minPtCharmedBaryon )
890+ if (ptCharmedBaryon < ptCharmedBaryonMin )
891891 continue ;
892- if (ptCharmedBaryon > maxPtCharmedBaryon )
892+ if (ptCharmedBaryon > ptCharmedBaryonMax )
893893 continue ;
894- if (impactProduct < minImpProduct || impactProduct > maxImpProduct )
894+ if (impactProduct < impProductMin || impactProduct > impProductMax )
895895 continue ;
896896 registry.fill (HIST (" hDecayLength" ), decayLength * 1e4 );
897897 registry.fill (HIST (" hDecayLengthScaled" ), decayLength * o2::constants::physics::MassOmegaC0 / RecoDecay::p (momenta[0 ], momenta[1 ]) * 1e4 );
0 commit comments