70 std::vector<TauolaParticle*> tau_list =
findParticles(pdg_id);
71 std::vector<TauolaParticle*> stable_tau_list;
73 for(
int i=0; i<(int) tau_list.size(); i++){
75 if(!tau_list.at(i)->hasDaughters())
76 stable_tau_list.push_back(tau_list.at(i));
79 std::vector<TauolaParticle*> t = tau_list.at(i)->getDaughters();
81 if(t.size()==1)
continue;
82 if(t.size()==2 && (abs(t[0]->getPdgID())==15 || abs(t[1]->getPdgID())==15) )
continue;
83 Log::Warning()<<
"Particle with pdg code "<<tau_list.at(i)->getPdgID()
84 <<
" already has daughters" <<endl;
88 return stable_tau_list;