diff --git a/PWGHF/TableProducer/treeCreatorOmegacSt.cxx b/PWGHF/TableProducer/treeCreatorOmegacSt.cxx index e5727d54bcb..92947c752d8 100644 --- a/PWGHF/TableProducer/treeCreatorOmegacSt.cxx +++ b/PWGHF/TableProducer/treeCreatorOmegacSt.cxx @@ -17,16 +17,16 @@ /// \author Ruiqi Yin #include "PWGHF/DataModel/TrackIndexSkimmingTables.h" -#include "PWGHF/Utils/utilsTrkCandHf.h" #include "PWGHF/Utils/utilsEvSelHf.h" +#include "PWGHF/Utils/utilsTrkCandHf.h" #include "PWGLF/DataModel/LFStrangenessTables.h" #include "Common/Core/RecoDecay.h" #include "Common/Core/Zorro.h" #include "Common/Core/ZorroSummary.h" #include "Common/Core/trackUtilities.h" -#include "Common/DataModel/CollisionAssociationTables.h" #include "Common/DataModel/Centrality.h" +#include "Common/DataModel/CollisionAssociationTables.h" #include "Common/DataModel/EventSelection.h" #include "Common/DataModel/PIDResponseTOF.h" #include "Common/DataModel/PIDResponseTPC.h" @@ -278,13 +278,13 @@ struct HfTreeCreatorOmegacSt { // whether to require TOF in addition to TPC Configurable useTofPid{"useTofPid", false, "Require TOF PID together with TPC PID when true; otherwise use TPC only"}; - //pT cut for CharmBaryon + // pT cut for CharmBaryon Configurable minPtCharmedBaryon{"minPtCharmedBaryon", 0.f, "Minimum pT for the charmed baryon"}; Configurable maxPtCharmedBaryon{"maxPtCharmedBaryon", 50.f, "Maximum pT for the charmed baryon"}; // CPA cuts (defaults: >= 0.9) Configurable minCpaCharmedBaryon{"minCpaCharmedBaryon", 0.9f, "Minimum CPA for charmed baryon"}; - Configurable minCpaXYCharmedBaryon{"minCpaXYCharmedBaryon", 0.9f, "Minimum CPA XY for charmed baryon"};//remove maybe??? + Configurable minCpaXYCharmedBaryon{"minCpaXYCharmedBaryon", 0.9f, "Minimum CPA XY for charmed baryon"}; // remove maybe??? Configurable minCpaCasc{"minCpaCasc", 0.9f, "Minimum CPA for cascade"}; Configurable minCpaXYCasc{"minCpaXYCasc", 0.9f, "Minimum CPA XY for cascade"}; @@ -494,15 +494,15 @@ struct HfTreeCreatorOmegacSt { const auto matCorr = static_cast(materialCorrectionType.value); for (const auto& collision : collisions) { - // Event selection & centrality using HF helper - float centrality{-1.f}; - const auto rejectionMask = hfEvSel.getHfCollisionRejectionMask(collision, centrality, ccdb, registry); - if (rejectionMask != 0) { - continue; // reject collisions failing HF event selection - } - // record event-level monitoring (centrality, vertex, etc.) - hfEvSel.fillHistograms(collision, rejectionMask, centrality); - const auto bc = collision.template bc_as(); + // Event selection & centrality using HF helper + float centrality{-1.f}; + const auto rejectionMask = hfEvSel.getHfCollisionRejectionMask(collision, centrality, ccdb, registry); + if (rejectionMask != 0) { + continue; // reject collisions failing HF event selection + } + // record event-level monitoring (centrality, vertex, etc.) + hfEvSel.fillHistograms(collision, rejectionMask, centrality); + const auto bc = collision.template bc_as(); if (runNumber != bc.runNumber()) { if (skimmedProcessing) { if (runNumber == 0) { @@ -696,7 +696,7 @@ struct HfTreeCreatorOmegacSt { (track.tpcNClsCrossedRows() >= TpcNclsFindableFraction * track.tpcNClsFindable()) && (track.tpcChi2NCl() <= TpcChi2NclMax) && (track.itsChi2NCl() <= ItsChi2NclMax)) { - //TPC and TOF PID for pion/kaon from chambaryon + // TPC and TOF PID for pion/kaon from chambaryon const bool passTPCpid = (std::abs(track.tpcNSigmaPi()) < maxNSigmaPion) || (std::abs(track.tpcNSigmaKa()) < maxNSigmaKaon); const float tofPiAbs = std::abs(track.tofNSigmaPi()); const float tofKaAbs = std::abs(track.tofNSigmaKa()); @@ -799,10 +799,10 @@ struct HfTreeCreatorOmegacSt { //--- do the MC Rec match if (mcParticles) { auto arrayDaughters = std::array{ - trackId.template track_as(), // bachelor <- charm baryon - casc.template bachelor_as(), // bachelor <- cascade - v0.template posTrack_as(), // p <- lambda - v0.template negTrack_as()}; // pi <- lambda + trackId.template track_as(), // bachelor <- charm baryon + casc.template bachelor_as(), // bachelor <- cascade + v0.template posTrack_as(), // p <- lambda + v0.template negTrack_as()}; // pi <- lambda auto arrayDaughtersCasc = std::array{ casc.template bachelor_as(), // bachelor <- cascade @@ -896,7 +896,7 @@ struct HfTreeCreatorOmegacSt { bachelor.tpcNSigmaKa(), bachelor.tofNSigmaKa(), momenta[0][0], // cascade momentum - momenta[0][1], + momenta[0][1], momenta[0][2], static_cast(trackCasc.sign() > 0), momenta[1][0], // pion/kaon momentum @@ -950,46 +950,45 @@ struct HfTreeCreatorOmegacSt { } } } - - // No centrality selection - void processDataNocent(Collisions const& collisions, - soa::SmallGroups const& trackedCascades, - aod::TrackAssoc const& trackIndices, - aod::Cascades const&, - aod::V0s const&, - TracksExt const&, - aod::BCsWithTimestamps const&) - { - fillTable(collisions, trackedCascades, trackIndices); - } - PROCESS_SWITCH(HfTreeCreatorOmegacSt, processDataNocent, "Process data (No centrality)", true); - - // FT0C centrality selection - void processDataFT0C(CollisionsFT0C const& collisions, - soa::SmallGroups const& trackedCascades, - aod::TrackAssoc const& trackIndices, - aod::Cascades const&, - aod::V0s const&, - TracksExt const&, - aod::BCsWithTimestamps const&) - { - fillTable(collisions, trackedCascades, trackIndices); - } - PROCESS_SWITCH(HfTreeCreatorOmegacSt, processDataFT0C, "Process data (FT0C centrality)", false); - - // FT0M centrality selection - void processDataFT0M(CollisionsFT0M const& collisions, - soa::SmallGroups const& trackedCascades, - aod::TrackAssoc const& trackIndices, - aod::Cascades const&, - aod::V0s const&, - TracksExt const&, - aod::BCsWithTimestamps const&) - { - fillTable(collisions, trackedCascades, trackIndices); - } - PROCESS_SWITCH(HfTreeCreatorOmegacSt, processDataFT0M, "Process data (FT0M centrality)", false); + // No centrality selection + void processDataNocent(Collisions const& collisions, + soa::SmallGroups const& trackedCascades, + aod::TrackAssoc const& trackIndices, + aod::Cascades const&, + aod::V0s const&, + TracksExt const&, + aod::BCsWithTimestamps const&) + { + fillTable(collisions, trackedCascades, trackIndices); + } + PROCESS_SWITCH(HfTreeCreatorOmegacSt, processDataNocent, "Process data (No centrality)", true); + + // FT0C centrality selection + void processDataFT0C(CollisionsFT0C const& collisions, + soa::SmallGroups const& trackedCascades, + aod::TrackAssoc const& trackIndices, + aod::Cascades const&, + aod::V0s const&, + TracksExt const&, + aod::BCsWithTimestamps const&) + { + fillTable(collisions, trackedCascades, trackIndices); + } + PROCESS_SWITCH(HfTreeCreatorOmegacSt, processDataFT0C, "Process data (FT0C centrality)", false); + + // FT0M centrality selection + void processDataFT0M(CollisionsFT0M const& collisions, + soa::SmallGroups const& trackedCascades, + aod::TrackAssoc const& trackIndices, + aod::Cascades const&, + aod::V0s const&, + TracksExt const&, + aod::BCsWithTimestamps const&) + { + fillTable(collisions, trackedCascades, trackIndices); + } + PROCESS_SWITCH(HfTreeCreatorOmegacSt, processDataFT0M, "Process data (FT0M centrality)", false); // MC reco: No centrality selection void processMcRecNocent(Collisions const& collisions,