#include <RandomXiThetaGunProducer.h>

Inheritance diagram for edm::RandomXiThetaGunProducer:

Public Member Functions
	RandomXiThetaGunProducer (const ParameterSet &)

	~RandomXiThetaGunProducer () override=default

Public Member Functions inherited from edm::one::EDProducer<>
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	generateParticle (double z_sign, double mass, unsigned int barcode, HepMC::GenVertex *vtx) const

void	produce (Event &, const EventSetup &) override

Private Attributes
double	energy_

CLHEP::HepRandomEngine *	engine_

unsigned int	nParticlesSector45_

unsigned int	nParticlesSector56_

unsigned int	particleId_

double	theta_x_mean_

double	theta_x_sigma_

double	theta_y_mean_

double	theta_y_sigma_

unsigned int	verbosity_

double	xi_max_

double	xi_min_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 19 of file RandomXiThetaGunProducer.h.

Constructor & Destructor Documentation

RandomXiThetaGunProducer::RandomXiThetaGunProducer ( const ParameterSet & iConfig )

Definition at line 24 of file RandomXiThetaGunProducer.cc.

                                                                                  :
   verbosity_         (iConfig.getUntrackedParameter<unsigned int>("verbosity", 0)),
   particleId_        (iConfig.getParameter<unsigned int>("particleId")),
   energy_            (iConfig.getParameter<double>("energy")),
   xi_min_            (iConfig.getParameter<double>("xi_min")),
   xi_max_            (iConfig.getParameter<double>("xi_max")),
   theta_x_mean_      (iConfig.getParameter<double>("theta_x_mean")),
   theta_x_sigma_     (iConfig.getParameter<double>("theta_x_sigma")),
   theta_y_mean_      (iConfig.getParameter<double>("theta_y_mean")),
   theta_y_sigma_     (iConfig.getParameter<double>("theta_y_sigma")),
   nParticlesSector45_(iConfig.getParameter<unsigned int>("nParticlesSector45")),
   nParticlesSector56_(iConfig.getParameter<unsigned int>("nParticlesSector56")),
   engine_(nullptr)
 {
   produces<HepMCProduct>("unsmeared");
 }

edm::RandomXiThetaGunProducer::~RandomXiThetaGunProducer ( )

overridedefault

Member Function Documentation

void RandomXiThetaGunProducer::generateParticle	(	double	z_sign,
		double	mass,
		unsigned int	barcode,
		HepMC::GenVertex *	vtx
	)		const

private

Definition at line 79 of file RandomXiThetaGunProducer.cc.

References MillePedeFileConverter_cfg::e, energy_, engine_, GenParticle::GenParticle, AlCaHLTBitMon_ParallelJobs::p, particleId_, mathSSE::sqrt(), theta_x_mean_, theta_x_sigma_, theta_y_mean_, theta_y_sigma_, verbosity_, protons_cff::xi, xi_max_, and xi_min_.

Referenced by produce().

 {
   const double xi = CLHEP::RandFlat::shoot(engine_, xi_min_, xi_max_);
   const double theta_x = CLHEP::RandGauss::shoot(engine_, theta_x_mean_, theta_x_sigma_);
   const double theta_y = CLHEP::RandGauss::shoot(engine_, theta_y_mean_, theta_y_sigma_);
 
   if (verbosity_)
     LogInfo("RandomXiThetaGunProducer") << "xi = " << xi << ", theta_x = " << theta_x
       << ", theta_y" << theta_y << ", z_sign = " << z_sign;
 
   const double cos_theta = sqrt(1. - theta_x*theta_x - theta_y*theta_y);
 
   const double p_nom = sqrt(energy_*energy_ - mass*mass);
   const double p = p_nom * (1. - xi);
   const double e = sqrt(p*p + mass*mass);
 
   HepMC::FourVector momentum(
     p * theta_x,
     p * theta_y,
     z_sign * p * cos_theta,
     e
   );
 
   HepMC::GenParticle* particle = new HepMC::GenParticle(momentum, particleId_, 1);
   particle->suggest_barcode(barcode);
   vtx->add_particle_out(particle);
 }

void RandomXiThetaGunProducer::produce	(	edm::Event &	e,
		const EventSetup &	es
	)

overrideprivate

Definition at line 43 of file RandomXiThetaGunProducer.cc.

References engine_, edm::EventID::event(), generateParticle(), edm::EventSetup::getData(), edm::RandomNumberGenerator::getEngine(), mps_fire::i, edm::EventBase::id(), ResonanceBuilder::mass, eostools::move(), nParticlesSector45_, nParticlesSector56_, convertSQLitetoXML_cfg::output, source_particleGun_cfi::ParticleID, particleId_, edm::Event::put(), edm::Event::streamID(), and extraflags_cff::vtx.

 {
   // get conditions
   edm::Service<edm::RandomNumberGenerator> rng;
   engine_ = &rng->getEngine(e.streamID());
 
   ESHandle<HepPDT::ParticleDataTable> pdgTable;
   es.getData(pdgTable);
 
   // prepare HepMC event
   HepMC::GenEvent *fEvt = new HepMC::GenEvent();
   fEvt->set_event_number(e.id().event());
 
   HepMC::GenVertex *vtx = new HepMC::GenVertex(HepMC::FourVector(0., 0., 0., 0.));
   fEvt->add_vertex(vtx);
 
   const HepPDT::ParticleData *pData = pdgTable->particle(HepPDT::ParticleID(particleId_));
   double mass = pData->mass().value();
 
   // generate particles
   unsigned int barcode = 0;
 
   for (unsigned int i = 0; i < nParticlesSector45_; ++i)
     generateParticle(+1., mass, ++barcode, vtx);
 
   for (unsigned int i = 0; i < nParticlesSector56_; ++i)
     generateParticle(-1., mass, ++barcode, vtx);
 
   // save output
   std::unique_ptr<HepMCProduct> output(new HepMCProduct()) ;
   output->addHepMCData(fEvt);
   e.put(std::move(output), "unsmeared");
 }