Algorithm:

• Find two highest-pt B-quarks (B-jets)
• Among the rest of the particles (jets), find two highest-pt quarks (quark jets), i.e. KFJET=0..6. Call them "tagging quarks (jets)"
• In the Eta of both tagging quarks (jets) histo, save Eta even if only one quark (jet) was found. Thus, we have slighly less than 200k entries.
• In the Minimum |Eta| histo, save min |Eta| out of the two tagging quarks (jets) - i.e., we have 1 entry per good event. Bad events are saved with |Eta| = -1 or -2.

PYTHIA	ATLFAST
Pythia with GenStat=3 (hard scattered particles): -1 means that the event failed because we couldn't find two b-quarks. -->I DON'T UNDERSTAND WHY PYTHIA STORES EVENTS LIKE THIS! -2 means that the event failed because we couldn't find two tagging-quarks	Atlfast: -1 means that the event failed because we couldn't find two b-jets. -->Because b-tagging extends to Eta<2.5, ~40% of all events fail this cut. -2 means that the event failed because we couldn't find at least two tagging-jets
Pythia with GenStat=3 (hard scattered particles): #(tagging quarks) = #(total quarks) - 2 (because we already identified the two b's from H->bbar)	Atlfast: #(tagging jets) = #(total quark-jets) - 2 (because we already identified the two b-jets from H->bbar)

---

Here is the # of b-quarks per event on Pythia level. I have no idea why Pythia saved events with less than 2 b-quarks!!!