@@ -331,25 +331,26 @@ bool MCSignal::CheckProng(int i, bool checkSources, const T& track)
331331 break ;
332332 }
333333 }
334- } /* else { // check generation in time
335- if (!currentMCParticle.has_daughters()) {
336- return false;
337- }
338- const auto& daughtersSlice = currentMCParticle.template daughters_as<P>();
339- for (auto& d : daughtersSlice) {
340- if (!fProngs[i].fExcludePDGInHistory[k] && fProngs[i].ComparePDG(d.pdgCode(), fProngs[i].fPDGInHistory[k], true, fProngs[i].fExcludePDGInHistory[k])) {
341- pdgInHistory.emplace_back(d.pdgCode());
342- break;
343- }
344- if (fProngs[i].fExcludePDGInHistory[k] && !fProngs[i].ComparePDG(d.pdgCode(), fProngs[i].fPDGInHistory[k], true, fProngs[i].fExcludePDGInHistory[k])) {
345- return false;
346- }
347- ith++;
348- if (ith > 10) { // need error message. Given pdg code was not found within 10 generations of the particles decay chain.
349- break;
350- }
351- }
352- }*/
334+ }
335+ // else { // check generation in time
336+ // if (!currentMCParticle.has_daughters()) {
337+ // return false;
338+ // }
339+ // const auto& daughtersSlice = currentMCParticle.template daughters_as<P>();
340+ // for (auto& d : daughtersSlice) {
341+ // if (!fProngs[i].fExcludePDGInHistory[k] && fProngs[i].ComparePDG(d.pdgCode(), fProngs[i].fPDGInHistory[k], true, fProngs[i].fExcludePDGInHistory[k])) {
342+ // pdgInHistory.emplace_back(d.pdgCode());
343+ // break;
344+ // }
345+ // if (fProngs[i].fExcludePDGInHistory[k] && !fProngs[i].ComparePDG(d.pdgCode(), fProngs[i].fPDGInHistory[k], true, fProngs[i].fExcludePDGInHistory[k])) {
346+ // return false;
347+ // }
348+ // ith++;
349+ // if (ith > 10) { // need error message. Given pdg code was not found within 10 generations of the particles decay chain.
350+ // break;
351+ // }
352+ // }
353+ // }
353354 }
354355 if (pdgInHistory.size () != nIncludedPDG) { // vector has as many entries as mothers (daughters) defined for prong
355356 return false ;
0 commit comments