626 |
|
_cosdphiMetgg = TMath::Cos(_dphiMetgg); |
627 |
|
_dphiPhPh = MathUtils::DeltaPhi((phHard->Mom()).Phi(),(phSoft->Mom()).Phi()); |
628 |
|
_mass = (phHard->Mom()+phSoft->Mom()).M(); |
629 |
< |
_masserr = 0.5*_mass*TMath::Sqrt(phHard->EnergyErr()*phHard->EnergyErr()/phHard->E()/phHard->E() + phSoft->EnergyErr()*phSoft->EnergyErr()/phSoft->E()/phSoft->E()); |
630 |
< |
_masserrsmeared = 0.5*_mass*TMath::Sqrt(phHard->EnergyErrSmeared()*phHard->EnergyErrSmeared()/phHard->E()/phHard->E() + phSoft->EnergyErrSmeared()*phSoft->EnergyErrSmeared()/phSoft->E()/phSoft->E()); |
629 |
> |
_masserr = 0.5*TMath::Sqrt(phHard->EnergyErr()*phHard->EnergyErr() + phSoft->EnergyErr()*phSoft->EnergyErr()); |
630 |
> |
_masserrsmeared = 0.5*TMath::Sqrt(phHard->EnergyErrSmeared()*phHard->EnergyErrSmeared() + phSoft->EnergyErrSmeared()*phSoft->EnergyErrSmeared()); |
631 |
|
_ptgg = (phHard->Mom()+phSoft->Mom()).Pt(); |
632 |
|
_etagg = (phHard->Mom()+phSoft->Mom()).Eta(); |
633 |
|
_phigg = (phHard->Mom()+phSoft->Mom()).Phi(); |
635 |
|
_evtcat = PhotonTools::DiphotonR9EtaPtCat(phHard,phSoft); |
636 |
|
|
637 |
|
const Double_t dz = sqrt(2.0)*fBeamspotWidth; |
638 |
< |
Double_t deltamvtx = _mass*VertexTools::DeltaMassVtx(phHard->CaloPos().X(), |
638 |
> |
Double_t deltamvtx = VertexTools::DeltaMassVtx(phHard->CaloPos().X(), |
639 |
|
phHard->CaloPos().Y(), |
640 |
|
phHard->CaloPos().Z(), |
641 |
|
phSoft->CaloPos().X(), |