void
PrintAllParticles ()
{
int i;
fflush(stdout);
printf(" PARTICLE POSITIONS\n\n");
for (i = 0; i < Total_Particles; i++) {
PrintParticle(Particle_List[i]);
printf("\n");
}
}
void GenTopEvent::Print() const{
std::cout << "top" << std::endl;
PrintParticle(GetTop());
std::cout << "top decay quark" << std::endl;
PrintParticle(GetTopDecayQuark());
std::cout << "higgs" << std::endl;
PrintParticle(GetHiggs());
std::cout << "higgs decay products" << std::endl;
PrintParticles(GetHiggsDecayProducts());
std::cout << "topbar" << std::endl;
PrintParticle(GetTopBar());
std::cout << "topbar decay quark" << std::endl;
PrintParticle(GetTopBarDecayQuark());
std::cout << "wplus" << std::endl;
PrintParticle(GetWplus());
std::cout << "wplus decay products" << std::endl;
PrintParticles(GetWplusDecayProducts());
std::cout << "wminus" << std::endl;
PrintParticle(GetWminus());
std::cout << "wminus decay products" << std::endl;
PrintParticles(GetWminusDecayProducts());
std::cout << "all top hads" << std::endl;
PrintParticles(GetAllTopHads());
std::cout << "all W hads" << std::endl;
PrintParticles(GetAllWhads());
std::cout << "all tophad bs" << std::endl;
PrintParticles(GetAllTopHadDecayQuarks());
std::cout << "all tophad W quarks" << std::endl;
PrintParticles(GetAllWQuarks());
std::cout << "all tophad W anti quarks" << std::endl;
PrintParticles(GetAllWAntiQuarks());
std::cout << "all top leps" << std::endl;
PrintParticles(GetAllTopLeps());
std::cout << "all w leps" << std::endl;
PrintParticles(GetAllWleps());
std::cout << "all top lep bs" << std::endl;
PrintParticles(GetAllTopLepDecayQuarks());
std::cout << "all leptons" << std::endl;
PrintParticles(GetAllLeptons());
std::cout << "all neutrinos" << std::endl;
PrintParticles(GetAllNeutrinos());
std::cout << "============================" << std::endl;
}
void GenTopEvent::PrintParticles(std::vector<reco::GenParticle> ps) const{
for(auto p=ps.begin();p!=ps.end();p++){
PrintParticle(*p);
}
}