status code for b-tagging

[antara001]

Hi Delphes Team,

I am trying to study events where a Higgs decays to b-quarks. So all my events have two b-quarks at the parton level. I used MG5 for generating the MC events and showered them with Pythia8. Next I have performed jet clustering and detector effects using Delphes (with delphes_card_CMS.tcl)

I have two queries in it:

1. I am doing b-tagging by comparing b-quark that is the closest in DR to the reconstructed jet.

So, for this I used particles in GenParticle class: if(particle->Status==3 && abs(particle->PID)==5) nb++;
and cout nb, but it is all 0. Instead, printing their status, I found the b-quarks have statuses : 23, 51, 52, 71. Could you tell why this is so and what these statuses mean?

2. I was also trying to do the same by looping over the jet constituents, if I could find a match for the b-jet by accessing the GenParticle in each reconstructed jet. However, here I could not find any b-quark within the jet constituents. I hope it is a correct approach.

Am I missing something here? Is this expected or are there some other factors playing a role here?

It will be really helpful if you can help me with understanding this.

Thanks and regards,
Antara