Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 14:21:07

0001 
0002 
0003 #include "IOMC/EventVertexGenerators/interface/GaussEvtVtxGenerator.h"
0004 #include "FWCore/Utilities/interface/Exception.h"
0005 
0006 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0007 
0008 #include "CLHEP/Random/RandGaussQ.h"
0009 #include "CLHEP/Units/GlobalSystemOfUnits.h"
0010 #include "CLHEP/Units/GlobalPhysicalConstants.h"
0011 //#include "CLHEP/Vector/ThreeVector.h"
0012 #include "HepMC/SimpleVector.h"
0013 
0014 GaussEvtVtxGenerator::GaussEvtVtxGenerator(const edm::ParameterSet& p) : BaseEvtVtxGenerator(p) {
0015   fMeanX = p.getParameter<double>("MeanX") * cm;
0016   fMeanY = p.getParameter<double>("MeanY") * cm;
0017   fMeanZ = p.getParameter<double>("MeanZ") * cm;
0018   fSigmaX = p.getParameter<double>("SigmaX") * cm;
0019   fSigmaY = p.getParameter<double>("SigmaY") * cm;
0020   fSigmaZ = p.getParameter<double>("SigmaZ") * cm;
0021   fTimeOffset = p.getParameter<double>("TimeOffset") * ns * c_light;
0022 
0023   if (fSigmaX < 0) {
0024     throw cms::Exception("Configuration") << "Error in GaussEvtVtxGenerator: "
0025                                           << "Illegal resolution in X (SigmaX is negative)";
0026   }
0027   if (fSigmaY < 0) {
0028     throw cms::Exception("Configuration") << "Error in GaussEvtVtxGenerator: "
0029                                           << "Illegal resolution in Y (SigmaY is negative)";
0030   }
0031   if (fSigmaZ < 0) {
0032     throw cms::Exception("Configuration") << "Error in GaussEvtVtxGenerator: "
0033                                           << "Illegal resolution in Z (SigmaZ is negative)";
0034   }
0035 }
0036 
0037 GaussEvtVtxGenerator::~GaussEvtVtxGenerator() {}
0038 
0039 HepMC::FourVector GaussEvtVtxGenerator::newVertex(CLHEP::HepRandomEngine* engine) const {
0040   double X, Y, Z, T;
0041   X = CLHEP::RandGaussQ::shoot(engine, fMeanX, fSigmaX);
0042   Y = CLHEP::RandGaussQ::shoot(engine, fMeanY, fSigmaY);
0043   Z = CLHEP::RandGaussQ::shoot(engine, fMeanZ, fSigmaZ);
0044   T = CLHEP::RandGaussQ::shoot(engine, fTimeOffset, fSigmaZ);
0045 
0046   return HepMC::FourVector(X, Y, Z, T);
0047 }
0048 
0049 void GaussEvtVtxGenerator::sigmaX(double s) {
0050   if (s >= 0) {
0051     fSigmaX = s;
0052   } else {
0053     throw cms::Exception("LogicError") << "Error in GaussEvtVtxGenerator::sigmaX: "
0054                                        << "Illegal resolution in X (negative)";
0055   }
0056 }
0057 
0058 void GaussEvtVtxGenerator::sigmaY(double s) {
0059   if (s >= 0) {
0060     fSigmaY = s;
0061   } else {
0062     throw cms::Exception("LogicError") << "Error in GaussEvtVtxGenerator::sigmaY: "
0063                                        << "Illegal resolution in Y (negative)";
0064   }
0065 }
0066 
0067 void GaussEvtVtxGenerator::sigmaZ(double s) {
0068   if (s >= 0) {
0069     fSigmaZ = s;
0070   } else {
0071     throw cms::Exception("LogicError") << "Error in GaussEvtVtxGenerator::sigmaZ: "
0072                                        << "Illegal resolution in Z (negative)";
0073   }
0074 }