| 86 |
|
_np=0; |
| 87 |
|
_npc=0; |
| 88 |
|
|
| 89 |
< |
vector<fastjet::PseudoJet> sortedbJets; //create a pseudojet |
| 90 |
< |
sortedbJets = sorted_by_pt(ufs->bJetList()); //psuedojet is filled with bjets that are sorted by pt |
| 91 |
< |
|
| 89 |
> |
vector<TLorentzVector> sortedbJets = ufs->bJetList(); |
| 90 |
|
// below is to work out what taulist and gentaulist is |
| 91 |
|
int countertau_fin(0), countertau_nf(0), countermu_fin(0), countermu_nf(0); |
| 92 |
|
for (int i = 0; i < pythia->event.size(); ++i){ |
| 191 |
|
|
| 192 |
|
//checks the number of bjets with pt over cut |
| 193 |
|
for(unsigned int bi =0; bi<ufs->bJetList().size(); bi++){ |
| 194 |
< |
if(abs(sortedbJets[bi].perp())>bjet_ptcut) bjet_counter25++; //this counts the bjets that pass cutt |
| 195 |
< |
bjet_pt->Fill(abs(sortedbJets[bi].perp())); |
| 194 |
> |
if(abs(sortedbJets[bi].Et())>bjet_ptcut) bjet_counter25++; //this counts the bjets that pass cutt |
| 195 |
> |
bjet_pt->Fill(abs(sortedbJets[bi].Et())); |
| 196 |
|
} |
| 197 |
|
|
| 198 |
|
//implement the cut on event of btagging |
