Skip to content
Merged
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
31 changes: 17 additions & 14 deletions PWGHF/TableProducer/candidateCreatorSigmac0plusplus.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -85,6 +85,7 @@ struct HfCandidateCreatorSigmac0plusplus {
Configurable<bool> applyGlobalTrkWoDcaCutsSoftPi{"applyGlobalTrkWoDcaCutsSoftPi", false, "Switch on the application of the global-track w/o dca cuts for soft pion BEFORE ALL OTHER CUSTOM CUTS"};
Configurable<float> softPiEtaMax{"softPiEtaMax", 0.9f, "Soft pion max value for pseudorapidity (abs vale)"};
Configurable<float> softPiChi2Max{"softPiChi2Max", 36.f, "Soft pion max value for chi2 ITS"};
Configurable<int> softPiApplyCustomITSHitMap{"softPiApplyCustomITSHitMap", true, "Flag to enable/disable the application of the custom ITS hitmap requirement for the candidate soft pion"};
Configurable<int> softPiItsHitMap{"softPiItsHitMap", 127, "Soft pion ITS hitmap"};
Configurable<int> softPiItsHitsMin{"softPiItsHitsMin", 1, "Minimum number of ITS layers crossed by the soft pion among those in \"softPiItsHitMap\""};
Configurable<float> softPiDcaXYMax{"softPiDcaXYMax", 0.065, "Soft pion max dcaXY (cm)"};
Expand Down Expand Up @@ -164,22 +165,24 @@ struct HfCandidateCreatorSigmac0plusplus {
// ITS chi2
softPiCuts.SetMaxChi2PerClusterITS(softPiChi2Max);
// ITS hitmap
std::set<uint8_t> setSoftPiItsHitMap; // = {};
constexpr std::size_t NLayersIts = 7;
for (std::size_t idItsLayer = 0u; idItsLayer < NLayersIts; idItsLayer++) {
if (TESTBIT(softPiItsHitMap, idItsLayer)) {
setSoftPiItsHitMap.insert(static_cast<uint8_t>(idItsLayer));
if (softPiApplyCustomITSHitMap) {
std::set<uint8_t> setSoftPiItsHitMap; // = {};
constexpr std::size_t NLayersIts = 7;
for (std::size_t idItsLayer = 0u; idItsLayer < NLayersIts; idItsLayer++) {
if (TESTBIT(softPiItsHitMap, idItsLayer)) {
setSoftPiItsHitMap.insert(static_cast<uint8_t>(idItsLayer));
}
}
LOG(info) << "### ITS hitmap for soft pion";
LOG(info) << " >>> setSoftPiItsHitMap.size(): " << setSoftPiItsHitMap.size();
LOG(info) << " >>> Custom ITS hitmap dfchecked: ";
for (const auto& it : setSoftPiItsHitMap) {
LOG(info) << " Layer " << static_cast<int>(it) << " ";
}
LOG(info) << "############";
softPiCuts.SetRequireITSRefit();
softPiCuts.SetRequireHitsInITSLayers(softPiItsHitsMin, setSoftPiItsHitMap);
}
LOG(info) << "### ITS hitmap for soft pion";
LOG(info) << " >>> setSoftPiItsHitMap.size(): " << setSoftPiItsHitMap.size();
LOG(info) << " >>> Custom ITS hitmap dfchecked: ";
for (const auto it : setSoftPiItsHitMap) {
LOG(info) << " Layer " << static_cast<int>(it) << " ";
}
LOG(info) << "############";
softPiCuts.SetRequireITSRefit();
softPiCuts.SetRequireHitsInITSLayers(softPiItsHitsMin, setSoftPiItsHitMap);

/// CCDB for dcaXY, dcaZ recalculation of soft pions reassigned to another collision
ccdb->setURL(ccdbUrl);
Expand Down
Loading