9 #include "CLHEP/Random/RandGaussZiggurat.h"
12 #include "CLHEP/Vector/LorentzVector.h"
13 #include "CLHEP/Units/PhysicalConstants.h"
33 return StatusCode::SUCCESS;
40 return StatusCode::SUCCESS;
45 return StatusCode::SUCCESS;
56 const CLHEP::HepRotation rotx(CLHEP::HepRep3x3(1.0,0.0,0.0, 0.0,cosX,-sinX, 0.0,sinX,cosX));
59 const CLHEP::HepRotation roty(CLHEP::HepRep3x3(cosY,0.0,sinY, 0.0,1.0,0.0, -sinY,0.0,cosY));
62 const CLHEP::HepRotation rotz(CLHEP::HepRep3x3(cosZ,-sinZ,0.0, sinZ,cosZ,0.0, 0.0,0.0,1.0));
64 const CLHEP::HepRotation rot = rotx*roty*rotz;
66 return StatusCode::SUCCESS;
73 CLHEP::HepLorentzRotation
transform = CLHEP::HepLorentzRotation();
76 for(
auto particleIter: ge) {
79 return StatusCode::SUCCESS;
83 const CLHEP::HepLorentzRotation&
transform)
const
86 const HepMC::FourVector
mom =
p->momentum();
87 CLHEP::HepLorentzVector hv(
mom.px(),
mom.py(),
mom.pz(),
mom.e());
91 p->set_momentum(HepMC::FourVector(hv.px(),hv.py(),hv.pz(),hv.e()));