Simulation::GenEventRotator Class Reference

This tool takes a HepMC::GenEvent and applies rotations See https://twiki.cern.ch/twiki/bin/viewauth/Atlas/AtlasG4EventFilters for more details. More...

#include <GenEventRotator.h>

Inheritance diagram for Simulation::GenEventRotator:

Collaboration diagram for Simulation::GenEventRotator:

Public Member Functions
	GenEventRotator (const std::string &t, const std::string &n, const IInterface *p)
	Constructor with parameters.
StatusCode	initialize () override final
	Athena algtool's Hooks.
StatusCode	finalize () override final
	Athena algtool's Hooks.
StatusCode	initializeAthenaEvent ()
StatusCode	manipulate (HepMC::GenEvent &ge, const EventContext &ctx) const override final
	modifies the given GenEvent

Private Member Functions
StatusCode	initializeGenEvent (CLHEP::HepLorentzRotation &transform, const EventContext &ctx) const
	calculate the transformations that we want to apply to the particles in the current GenEvent
void	rotateParticle (HepMC::GenParticlePtr &p, const CLHEP::HepLorentzRotation &transform) const
	apply rotations to individual GenParticles

Private Attributes
Gaudi::Property< double >	m_xangle {this, "xAngle", 0.0, "rotation about x-axis"}
Gaudi::Property< double >	m_yangle {this, "yAngle", 0.0, "rotation about y-axis"}
Gaudi::Property< double >	m_zangle {this, "zAngle", 0.0, "rotation about z-axis"}

Detailed Description

This tool takes a HepMC::GenEvent and applies rotations See https://twiki.cern.ch/twiki/bin/viewauth/Atlas/AtlasG4EventFilters for more details.

Definition at line 28 of file GenEventRotator.h.

Constructor & Destructor Documentation

GenEventRotator()

Simulation::GenEventRotator::GenEventRotator	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p )

Constructor with parameters.

Constructor.

Definition at line 22 of file GenEventRotator.cxx.

    : base_class(t,n,p)
  {
  }

Member Function Documentation

finalize()

StatusCode Simulation::GenEventRotator::finalize ( )

finaloverride

Athena algtool's Hooks.

Definition at line 37 of file GenEventRotator.cxx.

  {
    ATH_MSG_VERBOSE("Finalizing ...");
    return StatusCode::SUCCESS;
  }

initialize()

StatusCode Simulation::GenEventRotator::initialize ( )

finaloverride

Athena algtool's Hooks.

Definition at line 30 of file GenEventRotator.cxx.

  {
    ATH_MSG_VERBOSE("Initializing ...");
    return StatusCode::SUCCESS;
  }

initializeAthenaEvent()

StatusCode Simulation::GenEventRotator::initializeAthenaEvent

(

)

Definition at line 43 of file GenEventRotator.cxx.

  {
    return StatusCode::SUCCESS;
  }

initializeGenEvent()

StatusCode Simulation::GenEventRotator::initializeGenEvent ( CLHEP::HepLorentzRotation & transform, const EventContext & ctx ) const

private

calculate the transformations that we want to apply to the particles in the current GenEvent

Definition at line 48 of file GenEventRotator.cxx.

  {
    // Reset the transformation
    transform = CLHEP::HepLorentzRotation(); //TODO drop this
 
    // Setting up three rotation matrices, via 3x3 rep. rather than even more messy EulerAngles
    const double sinX = sin(m_xangle*CLHEP::deg);
    const double cosX = cos(m_xangle*CLHEP::deg);
    const CLHEP::HepRotation rotx(CLHEP::HepRep3x3(1.0,0.0,0.0, 0.0,cosX,-sinX, 0.0,sinX,cosX));
    const double sinY = sin(m_yangle*CLHEP::deg);
    const double cosY = cos(m_yangle*CLHEP::deg);
    const CLHEP::HepRotation roty(CLHEP::HepRep3x3(cosY,0.0,sinY, 0.0,1.0,0.0, -sinY,0.0,cosY));
    const double sinZ = sin(m_zangle*CLHEP::deg);
    const double cosZ = cos(m_zangle*CLHEP::deg);
    const CLHEP::HepRotation rotz(CLHEP::HepRep3x3(cosZ,-sinZ,0.0, sinZ,cosZ,0.0, 0.0,0.0,1.0));
    // Combine the three rotations - order of trf execution is unimportant as they are orthogonal.
    const CLHEP::HepRotation rot = rotx*roty*rotz;
    transform.set(rot);
    return StatusCode::SUCCESS;
  }

manipulate()

StatusCode Simulation::GenEventRotator::manipulate ( HepMC::GenEvent & ge, const EventContext & ctx ) const

finaloverride

modifies the given GenEvent

modifies (displaces) the given GenEvent

Definition at line 70 of file GenEventRotator.cxx.

  {
    // Obtain the transformation
    CLHEP::HepLorentzRotation transform = CLHEP::HepLorentzRotation();
    ATH_CHECK( initializeGenEvent(transform, ctx) );
 
    for(auto particleIter:  ge) {
      rotateParticle(particleIter, transform);
    }
    return StatusCode::SUCCESS;
  }

rotateParticle()

void Simulation::GenEventRotator::rotateParticle ( HepMC::GenParticlePtr & p, const CLHEP::HepLorentzRotation & transform ) const

private

apply rotations to individual GenParticles

Definition at line 82 of file GenEventRotator.cxx.

  {
    // Apply the same transformation for EVERY HepMC::GenParticle
    const HepMC::FourVector mom = p->momentum();
    CLHEP::HepLorentzVector hv(mom.px(), mom.py(), mom.pz(), mom.e()); //TODO check units
    ATH_MSG_VERBOSE("initial momentum " << hv );
    hv.transform(transform);
    ATH_MSG_VERBOSE("transformed momentum " << hv);
    p->set_momentum(HepMC::FourVector(hv.px(),hv.py(),hv.pz(),hv.e())); //TODO check units
  }

Member Data Documentation

m_xangle

Gaudi::Property<double> Simulation::GenEventRotator::m_xangle {this, "xAngle", 0.0, "rotation about x-axis"}

private

Definition at line 49 of file GenEventRotator.h.

{this, "xAngle", 0.0, "rotation about x-axis"};

m_yangle

Gaudi::Property<double> Simulation::GenEventRotator::m_yangle {this, "yAngle", 0.0, "rotation about y-axis"}

private

Definition at line 50 of file GenEventRotator.h.

{this, "yAngle", 0.0, "rotation about y-axis"};

m_zangle

Gaudi::Property<double> Simulation::GenEventRotator::m_zangle {this, "zAngle", 0.0, "rotation about z-axis"}

private

Definition at line 51 of file GenEventRotator.h.

{this, "zAngle", 0.0, "rotation about z-axis"};

---

The documentation for this class was generated from the following files:

• GenEventRotator.h
• GenEventRotator.cxx