DefaultSmearer.h

Go to the documentation of this file.

// DefaultSmearer class
//
// Namespace Atlfast::
//
//-----------------------------------------------------------------
//
// This class provides a default implementation of the ISmearer
// interface. It manages the generation of random numbers and
// implements the smear() function required by ISmear
//
// Prototype- needs tidying up!
//
// Authors: H.T.Phillips, P. Clarke, E. Richter-Was, P. Sherwood, R. Steward
//
//
//

#ifndef ATLFAST_DEFAULTSMEARER_H
#define ATLFAST_DEFAULTSMEARER_H

//Atlfast includes
#include "AtlfastAlgs/ISmearer.h"

#include "CLHEP/Random/RandomEngine.h"
#include "CLHEP/Random/RandGauss.h"
#include "CLHEP/Random/RandFlat.h"

namespace Atlfast {

  class DefaultSmearer : virtual public ISmearer {

  public:
    //----------------------------------------
    // Constructors and Destructors
    //----------------------------------------
    DefaultSmearer (const int aseed );
    
    virtual ~DefaultSmearer();
    
    //----------------------------------------
    // The Smear Method
    //----------------------------------------
    virtual HepLorentzVector smear(const HepMC::GenParticle& );
    virtual HepLorentzVector smear( const HepLorentzVector& );
    
    virtual int setSmearParameters( const std::vector<double>& smearValues );
    virtual int setSmearParamSchema( const int smearSchema);
    
    //------------------------------------------------
    inline RandGauss* randGauss( ) { return m_randGauss;}
    //------------------------------------------------
    inline RandFlat* randFlat( ) { return m_randFlat;}

  private:
    //----------------------------------------
    // Private Methods
    //----------------------------------------

    //----------------------------------------
    // Data Members
    //---------------------------------------

    // For now, we manage the random number engine
    // ourselves. In the future we can use the framework
    // random number service (Gaudi 4 has one for example)
    // WE SHOULD USE THIS SERVICE ASAP!!
 
    //HepRandomEngine* m_randomEngine;
    RandGauss* m_randGauss;

    RandFlat* m_randFlat;

    // And the random number seed which I put here as a member
    // variable because I hope to allow it to be written with
    // jobOptions eventually...
    // Although right now it is not used anywhere other than the
    // constructor.
    int m_randSeed;

  };


} // end of namespace bracket

#endif

---