Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
10 changes: 5 additions & 5 deletions EventFiltering/PWGHF/HFFilter.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -171,7 +171,7 @@ struct HfFilter { // Main struct for HF triggers
std::array<std::shared_ptr<TH1>, kNCharmParticles> hCharmDeuteronKstarDistr{};
std::array<std::shared_ptr<TH2>, nTotBeautyParts> hMassVsPtB{};
std::array<std::shared_ptr<TH2>, kNCharmParticles + 19> hMassVsPtC{}; // +9 for resonances (D*+, D*0, Ds*+, Ds1+, Ds2*+, Xic+* right sign, Xic+* wrong sign, Xic0* right sign, Xic0* wrong sign) +2 for SigmaC (SigmaC++, SigmaC0) +2 for SigmaCK pairs (SigmaC++K-, SigmaC0K0s) +3 for charm baryons (Xi+Pi, Xi+Ka, Xi+Pi+Pi) + JPsi
std::array<std::shared_ptr<TH2>, 4> hPrDePID; // proton TPC, proton TOF, deuteron TPC, deuteron TOF
std::array<std::shared_ptr<TH2>, 4> hPrDePID; // proton TPC, proton TOF, deuteron TPC, deuteron TOF
std::array<std::shared_ptr<TH1>, kNCharmParticles> hBDTScoreBkg{};
std::array<std::shared_ptr<TH1>, kNCharmParticles> hBDTScorePrompt{};
std::array<std::shared_ptr<TH1>, kNCharmParticles> hBDTScoreNonPrompt{};
Expand Down Expand Up @@ -421,7 +421,7 @@ struct HfFilter { // Main struct for HF triggers
auto cand2ProngsThisColl = cand2Prongs.sliceBy(hf2ProngPerCollision, thisCollId);
for (const auto& cand2Prong : cand2ProngsThisColl) { // start loop over 2 prongs

int8_t preselD0 = TESTBIT(cand2Prong.hfflag(), o2::aod::hf_cand_2prong::DecayType::D0ToPiK); // check if it's a D0
int8_t preselD0 = TESTBIT(cand2Prong.hfflag(), o2::aod::hf_cand_2prong::DecayType::D0ToPiK); // check if it's a D0
int8_t preselJPsiToMuMu = TESTBIT(cand2Prong.hfflag(), o2::aod::hf_cand_2prong::DecayType::JpsiToMuMu); // check if it's a JPsi
if (preselD0 == 0 && preselJPsiToMuMu == 0) {
continue;
Expand Down Expand Up @@ -494,7 +494,7 @@ struct HfFilter { // Main struct for HF triggers
int8_t selD0InMass{0};
double massD0Cand{-1.}, massD0BarCand{-1.};
if (isD0SignalTagged) {
// single D0
// single D0
keepEvent[kSingleCharm2P] = true;
if (isD0BeautyTagged) {
keepEvent[kSingleNonPromptCharm2P] = true;
Expand All @@ -509,7 +509,7 @@ struct HfFilter { // Main struct for HF triggers
// multi-charm selection
if (isD0CharmTagged) {
indicesDau2Prong.push_back(std::vector<int64_t>{trackPos.globalIndex(), trackNeg.globalIndex()});
}
}

if (applyOptimisation) {
optimisationTreeCharm(thisCollId, o2::constants::physics::Pdg::kD0, pt2Prong, scores[0], scores[1], scores[2]);
Expand Down Expand Up @@ -710,7 +710,7 @@ struct HfFilter { // Main struct for HF triggers
std::array<float, 3> pVecPosVtx{}, pVecNegVtx{}, pVecThirdVtx{}, pVecFourthVtx{};
// 3-prong vertices
if (!keepEvent[kBtoJPsiKa] || !keepEvent[kBtoJPsiPi]) {
if(df3.process(trackParPos, trackParNeg, trackParThird) != 0) {
if (df3.process(trackParPos, trackParNeg, trackParThird) != 0) {
df3.propagateTracksToVertex();
const auto& secondaryVertexBto3tracks = df3.getPCACandidate();
df3.getTrack(0).getPxPyPzGlo(pVecPosVtx);
Expand Down
6 changes: 3 additions & 3 deletions EventFiltering/PWGHF/HFFilterHelpers.h
Original file line number Diff line number Diff line change
Expand Up @@ -301,8 +301,8 @@ static const std::vector<std::string> labelsRowsCutsPt = {"Minimum", "Maximum"};

// PID cuts
constexpr float cutsNsigma[4][8] = {
{3., 3., 3., 5., 3., 3., 5., 3.}, // TPC proton from Lc, pi/K from D0, K from 3-prong, femto selected proton, pi/K from Xic/Omegac, K from Xic*->SigmaC-Kaon, femto selected deuteron, K/p from beauty->JPsiX
{3., 3., 3., 2.5, 3., 3., 5., 3.}, // TOF proton from Lc, pi/K from D0, K from 3-prong, femto selected proton, pi/K from Xic/Omegac, K from Xic*->SigmaC-Kaon, femto selected deuteron, K/p from beauty->JPsiX
{3., 3., 3., 5., 3., 3., 5., 3.}, // TPC proton from Lc, pi/K from D0, K from 3-prong, femto selected proton, pi/K from Xic/Omegac, K from Xic*->SigmaC-Kaon, femto selected deuteron, K/p from beauty->JPsiX
{3., 3., 3., 2.5, 3., 3., 5., 3.}, // TOF proton from Lc, pi/K from D0, K from 3-prong, femto selected proton, pi/K from Xic/Omegac, K from Xic*->SigmaC-Kaon, femto selected deuteron, K/p from beauty->JPsiX
{999., 999., 999., 2.5, 999., 999., 5., 999.}, // Sum in quadrature of TPC and TOF (used only for femto selected proton and deuteron for pT < 4 GeV/c)
{999., 999., 999., 999., 999., 999., -4., 999.} // ITS used only for femto selected deuteron for less than pt threshold
};
Expand Down Expand Up @@ -2043,7 +2043,7 @@ inline bool HfFilterHelper::isSelectedBhadronInMassRange(T1 const& ptCand, T2 co
/// \param tracksDauNoMu is the array of tracks for the daughters that are no muons
/// \param primVtx is the primary vertex
/// \param secVtx is the secondary vertex
/// \param activateQA is the flag to enable the
/// \param activateQA is the flag to enable the
/// \param hMassVsPt is the array of histograms for QA
/// \return true if the beauty candidate passes all cuts
template <int Nprongs, typename T1, typename T2, typename T3, typename T4, typename H2>
Expand Down
Loading