18 double& wideJetBroadening,
double& totalJetBroadening,
19 CLHEP::Hep3Vector
thrust,
bool useThreeD )
21 if(abs(
thrust.mag()-1)>0.01)
43 CLHEP::Hep3Vector
c( (*itr)->px(), (*itr)->py(),
z );
47 etau += (*itr)->eta() *
c.mag();
48 phiu += (*itr)->phi() *
c.mag();
51 etad += (*itr)->eta() *
c.mag();
52 phid += (*itr)->phi() *
c.mag();
56 const double inv_Qu = 1. / Qu;
57 const double inv_Qd = 1. / Qd;
76 CLHEP::Hep3Vector
c( (*itr)->px(), (*itr)->py(),
z );
81 sqrt((etau-(*itr)->eta())*(etau-(*itr)->eta())+
82 (phiu-(*itr)->phi())*(phiu-(*itr)->phi()));
85 sqrt((etad-(*itr)->eta())*(etad-(*itr)->eta())+
86 (phid-(*itr)->phi())*(phid-(*itr)->phi()));
93 totalJetBroadening=Bu+Bd;