/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/IOMC/EventVertexGenerators/interface/GaussEvtVtxGenerator.h

Go to the documentation of this file.

#ifndef IOMC_GaussEvtVtxGenerator_H
#define IOMC_GaussEvtVtxGenerator_H

#include "IOMC/EventVertexGenerators/interface/BaseEvtVtxGenerator.h"

namespace CLHEP {
   class RandGaussQ;
}

class GaussEvtVtxGenerator : public BaseEvtVtxGenerator 
{
public:
  GaussEvtVtxGenerator(const edm::ParameterSet & p);
  virtual ~GaussEvtVtxGenerator();

  //virtual CLHEP::Hep3Vector* newVertex();
  virtual HepMC::FourVector* newVertex() ;

  virtual TMatrixD* GetInvLorentzBoost() {
          return 0;
  }

   
  void sigmaX(double s=1.0);
  void sigmaY(double s=1.0);
  void sigmaZ(double s=1.0);

  void meanX(double m=0) { fMeanX=m; }
  void meanY(double m=0) { fMeanY=m; }
  void meanZ(double m=0) { fMeanZ=m; }
  
private:
  GaussEvtVtxGenerator(const GaussEvtVtxGenerator &p);
  GaussEvtVtxGenerator&  operator = (const GaussEvtVtxGenerator & rhs );
private:
  double fSigmaX, fSigmaY, fSigmaZ;
  double fMeanX,  fMeanY,  fMeanZ;
  CLHEP::RandGaussQ*  fRandom ;
  double fTimeOffset;
};

#endif