---

DefaultSmearer.cxx

Go to the documentation of this file.

// DefaultSmearer.cxx
//
// Implementation of Default Smearer class
//
// Namespace Atlfast::
//
//--------------------------------------------------------------------
//
// Authors: H.T. Phillips, P. Clarke, E. Richter-Was, P. Sherwood, R. Steward
//
//
#include "AtlfastAlgs/DefaultSmearer.h"
// CLHEP includes
#include "CLHEP/Vector/LorentzVector.h"
#include "CLHEP/Random/JamesRandom.h"
#include "CLHEP/Random/RandGauss.h"





namespace Atlfast {

  DefaultSmearer::DefaultSmearer(const int aseed) {
    // instantiate the random number engine based on the seed and store it
    HepRandomEngine* randomEngine = new HepJamesRandom( aseed );
    m_randGauss    = new RandGauss( randomEngine );
    m_randSeed     = aseed;
  }


  DefaultSmearer::~DefaultSmearer() {
    // we check to see if our random engine and gauss are still set and if so kill them
    //    if (m_randomEngine) {
    //      delete m_randomEngine;
    //    }
    if (m_randGauss) {
      delete m_randGauss;
    }
  }

  HepLorentzVector DefaultSmearer::smear( const HepLorentzVector& avec) {
    // do the actual business... to fake it up, just smear in each variable ofthe four vector with width 1.0
    // for testing purposes only

    HepLorentzVector bvec;
    // set each component by smearing around the central value with width 1.0
    bvec.setE( m_randGauss->fire(avec.e(), 1.0) );
    bvec.setPx(m_randGauss->fire(avec.px(),1.0) );
    bvec.setPy(m_randGauss->fire(avec.py(),1.0) );
    bvec.setPz(m_randGauss->fire(avec.pz(),1.0) );

    return bvec;
  }



} // end of namespace bracket

---