/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/IOMC/ParticleGuns/src/FlatBaseThetaGunProducer.cc

Go to the documentation of this file.

#include <iostream>
#include <memory>

#include "FWCore/Framework/interface/Run.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/RandomNumberGenerator.h"
#include "FWCore/Utilities/interface/Exception.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/GenRunInfoProduct.h"

#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"

#include "IOMC/ParticleGuns/interface/FlatBaseThetaGunProducer.h"


using namespace edm;

namespace {

  CLHEP::HepRandomEngine& getEngineReference() {

    edm::Service<edm::RandomNumberGenerator> rng;
    if(!rng.isAvailable()) {
      throw cms::Exception("Configuration")
        << "The RandomNumberProducer module requires the RandomNumberGeneratorService\n"
        "which appears to be absent.  Please add that service to your configuration\n"
        "or remove the modules that require it.";
    }

    // The Service has already instantiated an engine.  Make contact with it.
    return (rng->getEngine());
  }
}

FlatBaseThetaGunProducer::FlatBaseThetaGunProducer(const edm::ParameterSet& pset) :
   fEvt(0),
   fRandomEngine(getEngineReference()),
   fRandomGenerator(0) {

  edm::ParameterSet pgun_params = pset.getParameter<edm::ParameterSet>("PGunParameters") ;
  
  fPartIDs  = pgun_params.getParameter<std::vector<int> >("PartID");  
  fMinTheta = pgun_params.getParameter<double>("MinTheta");
  fMaxTheta = pgun_params.getParameter<double>("MaxTheta");
  fMinPhi   = pgun_params.getParameter<double>("MinPhi");
  fMaxPhi   = pgun_params.getParameter<double>("MaxPhi");
  fVerbosity = pset.getUntrackedParameter<int>( "Verbosity",0 ) ;

  // The Service has already instantiated an engine.  Use it.
   fRandomGenerator = new CLHEP::RandFlat(fRandomEngine) ;
   fAddAntiParticle = pset.getParameter<bool>("AddAntiParticle") ;

   produces<GenRunInfoProduct, InRun>();
}

FlatBaseThetaGunProducer::~FlatBaseThetaGunProducer() {
//  if ( fRandomGenerator != NULL ) delete fRandomGenerator;
}


void FlatBaseThetaGunProducer::beginRun( edm::Run &r, const edm::EventSetup& es ) {
   es.getData( fPDGTable ) ;
}

void FlatBaseThetaGunProducer::endRun( Run &run, const EventSetup& es )
{
   // just create an empty product
   // to keep the EventContent definitions happy
   // later on we might put the info into the run info that this is a PGun
   std::auto_ptr<GenRunInfoProduct> genRunInfo( new GenRunInfoProduct() );
   run.put( genRunInfo );
}