| 291 |
|
- met->Px()*dilepton->Px() - met->Py()*dilepton->Py() |
| 292 |
|
- leptons->At(0)->Px()*leptons->At(1)->Px() - leptons->At(0)->Py()*leptons->At(1)->Py()); |
| 293 |
|
} |
| 294 |
< |
else if(nsel == 7){ // Use of M mass and mnu == 0, no fancy stuff |
| 295 |
< |
enell = TMath::Sqrt(dilepton->Pt()*dilepton->Pt() + dilepton->Mass()*dilepton->Mass()); |
| 296 |
< |
enenn = TMath::Sqrt(met->Pt()*met->Pt() + 0.0*0.0); |
| 297 |
< |
enex = dilepton->Px() + met->Px(); |
| 298 |
< |
eney = dilepton->Py() + met->Py(); |
| 299 |
< |
mtHiggs = (enell+enenn)*(enell+enenn) - enex*enex - eney*eney; |
| 294 |
> |
else if(nsel == 7){ // Use of M mass and mnu == 0 |
| 295 |
> |
double deltaPhiDileptonMet = MathUtils::DeltaPhi(dilepton->Phi(), |
| 296 |
> |
met->Phi()); |
| 297 |
> |
mtHiggs = 2.0*dilepton->Pt()*met->Pt()*(1.0 - cos(deltaPhiDileptonMet)); |
| 298 |
|
} |
| 299 |
|
|
| 300 |
|
if(nsel >= 0 && nsel <= 4){ |
