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"
#include "CLHEP/Random/RandFlat.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_randFlat     = new RandFlat( *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_randFlat) {
      delete m_randFlat;
    }
    
    if (m_randGauss) {
      delete m_randGauss;
    }
    
    //    if (m_randomEngine) {
    //      delete m_randomEngine;
    //    }
    
  }
  
  HepLorentzVector DefaultSmearer::smear(const HepMC::GenParticle& particle){
    return smear(particle.momentum());
  }
  
  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;
  }

  //cct: implement the setSmearParameters method
  int DefaultSmearer::setSmearParameters (const std::vector<double>& /*smearValues*/){
    return 0;
  }
  //cct: implement the setSmearParamSchema method
  int DefaultSmearer::setSmearParamSchema ( const int /*smearSchema*/){
    return 0;
  }
  
} // end of namespace bracket

---