#include <BeamHaloProducer.h>

Inheritance diagram for edm::BeamHaloProducer:

Public Member Functions
	BeamHaloProducer (const ParameterSet &)
	Constructor. More...

void	setRandomEngine (CLHEP::HepRandomEngine *v)

	~BeamHaloProducer () override
	Destructor. More...

Public Member Functions inherited from edm::one::EDProducer< EndRunProducer, one::WatchLuminosityBlocks, one::SharedResources >
	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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginLuminosityBlock (LuminosityBlock const &, EventSetup const &) override

bool	call_bh_set_parameters (int ival, float fval, const std::string cval_string)

bool	call_ki_bhg_fill (int &iret, float &weight)

bool	call_ki_bhg_init (long &seed)

bool	call_ki_bhg_stat (int &iret)

void	clear ()

void	endLuminosityBlock (LuminosityBlock const &, EventSetup const &) override

void	endRunProduce (Run &r, const EventSetup &es) override

void	produce (Event &e, const EventSetup &es) override

Private Attributes
float	EG_MAX_

float	EG_MIN_

HepMC::GenEvent *	evt

std::string	G3FNAME_

int	GENMOD_

bool	isInitialized_

int	IW_HAD_

int	IW_MUO_

int	LHC_B1_

int	LHC_B2_

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 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 23 of file BeamHaloProducer.h.

Constructor & Destructor Documentation

BeamHaloProducer::BeamHaloProducer ( const ParameterSet & pset )

Constructor.

Definition at line 59 of file BeamHaloProducer.cc.

References call_bh_set_parameters(), gather_cfg::cout, objects.autophobj::float, edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                              :
   evt(nullptr),
   isInitialized_(false)
 {
 
    int iparam[8];
    float fparam[4];
    std::string cparam;
  // -- from bhgctrl.inc
    iparam[0]  = pset.getUntrackedParameter<int>("GENMOD");
    iparam[1]  = pset.getUntrackedParameter<int>("LHC_B1");
    iparam[2]  = pset.getUntrackedParameter<int>("LHC_B2");
    iparam[3]  = pset.getUntrackedParameter<int>("IW_MUO");
    iparam[4]  = pset.getUntrackedParameter<int>("IW_HAD");
    iparam[5]  = 9999999;
    iparam[6]  = pset.getUntrackedParameter<int>("OFFSET",0);
    iparam[7]  = pset.getUntrackedParameter<int>("shift_bx");
    
    fparam[0]  = (float)pset.getUntrackedParameter<double>("EG_MIN");
    fparam[1]  = (float)pset.getUntrackedParameter<double>("EG_MAX");
 
    fparam[2] = (float)pset.getUntrackedParameter<double>("BXNS");
    fparam[3]  = (float)pset.getUntrackedParameter<double>("W0",1.0);
 
    cparam     = pset.getUntrackedParameter<std::string>("G3FNAME","input.txt");
    call_bh_set_parameters(iparam,fparam,cparam);
 
   produces<HepMCProduct>("unsmeared");
   produces<GenEventInfoProduct>();
   produces<GenRunInfoProduct, Transition::EndRun>();
 
   usesResource("BeamHaloProducer");
 
   cout << "BeamHaloProducer: starting event generation ... " << endl;
 }

BeamHaloProducer::~BeamHaloProducer ( )

override

Destructor.

Definition at line 53 of file BeamHaloProducer.cc.

                                     {
     int iret=0;
     call_ki_bhg_stat(iret);
 }

Member Function Documentation

void BeamHaloProducer::beginLuminosityBlock	(	LuminosityBlock const &	lumi,
		EventSetup const &
	)

overrideprivate

Definition at line 104 of file BeamHaloProducer.cc.

References call_ki_bhg_init(), edm::LuminosityBlock::index(), isInitialized_, and SurveyInfoScenario_cff::seed.

 {
   if(!isInitialized_) {
     isInitialized_ = true;
     RandomEngineSentry<BeamHaloProducer> randomEngineSentry(this, lumi.index());
 
     // -- initialisation
     long seed = 1; // This seed is not actually used
     call_ki_bhg_init(seed);
   }
 }

bool BeamHaloProducer::call_bh_set_parameters	(	int *	ival,
		float *	fval,
		const std::string	cval_string
	)

private

Definition at line 175 of file BeamHaloProducer.cc.

References BHSETPARAM.

Referenced by BeamHaloProducer().

                                                                                                  {
   BHSETPARAM(ival,fval,cval_string.c_str(),cval_string.length());
     return true;
 }

bool BeamHaloProducer::call_ki_bhg_fill	(	int &	iret,
		float &	weight
	)

private

Definition at line 185 of file BeamHaloProducer.cc.

References KI_BHG_FILL.

Referenced by produce().

                                                                 {
     KI_BHG_FILL(iret,weight);
     return true;
 }

bool BeamHaloProducer::call_ki_bhg_init ( long & seed )

private

Definition at line 180 of file BeamHaloProducer.cc.

References KI_BHG_INIT.

Referenced by beginLuminosityBlock().

                                                   {
     KI_BHG_INIT(seed);
     return true;
 }

bool BeamHaloProducer::call_ki_bhg_stat ( int & iret )

private

Definition at line 190 of file BeamHaloProducer.cc.

References KI_BHG_STAT.

                                                  {
     KI_BHG_STAT(iret);
     return true;
 }

void BeamHaloProducer::clear ( void )

private

Definition at line 96 of file BeamHaloProducer.cc.

97 {

98 }

void edm::BeamHaloProducer::endLuminosityBlock	(	LuminosityBlock const &	,
		EventSetup const &
	)

inlineoverrideprivate

Definition at line 44 of file BeamHaloProducer.h.

References hitfit::clear().

44 { }

void BeamHaloProducer::endRunProduce	(	Run &	r,
		const EventSetup &	es
	)

overrideprivate

Definition at line 166 of file BeamHaloProducer.cc.

References eostools::move(), and edm::Run::put().

 {
    // 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
    unique_ptr<GenRunInfoProduct> genRunInfo( new GenRunInfoProduct() );
    run.put(std::move(genRunInfo));
 }

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

overrideprivate

Definition at line 116 of file BeamHaloProducer.cc.

References call_ki_bhg_fill(), edm::EventID::event(), edm::errors::EventCorruption, evt, Exception, particlelevel_cff::genEventInfo, GenParticle::GenParticle, edm::EventBase::id(), eostools::move(), AlCaHLTBitMon_ParallelJobs::p, edm::Event::put(), edm::Event::streamID(), create_public_lumi_plots::weights, and wrapper.

                                                                {
 
   RandomEngineSentry<BeamHaloProducer> randomEngineSentry(this, e.streamID());
 
     // cout << "in produce " << endl;
 
   //        unique_ptr<HepMCProduct> bare_product(new HepMCProduct());
 
     // cout << "apres autoptr " << endl;
 
     int iret=0;
         float weight = 0;
         call_ki_bhg_fill(iret, weight);
 
 // Throw an exception if call_ki_bhg_fill(...) fails.  Use the EventCorruption
 // exception since it maps onto SkipEvent which is what we want to do here.
 
         if( iret < 0 )
           throw edm::Exception(edm::errors::EventCorruption)
             << "BeamHaloProducer: function call_ki_bhg_fill returned " << iret << endl;
 
     // cout << "apres fortran " << endl;
 
 
         // HepMC::GenEvent* evt = conv.getGenEventfromHEPEVT();
     //  HepMC::GenEvent* evt = conv.read_next_event();  seems to be broken (?)
   evt = new HepMC::GenEvent();
 
   for (int theindex = 1; theindex<=wrapper.number_entries(); theindex++) {
   HepMC::GenVertex* Vtx = new  HepMC::GenVertex(HepMC::FourVector(wrapper.x(theindex),wrapper.y(theindex),wrapper.z(theindex),wrapper.t(theindex)));
   HepMC::FourVector p(wrapper.px(theindex),wrapper.py(theindex),wrapper.pz(theindex),wrapper.e(theindex));
   HepMC::GenParticle* Part = 
     new HepMC::GenParticle(p,wrapper.id(theindex),wrapper.status(theindex));
   Vtx->add_particle_out(Part); 
   evt->add_vertex(Vtx);
   }
 
   evt->set_event_number(e.id().event());
 
     HepMC::WeightContainer& weights = evt -> weights();
     weights.push_back(weight);
     //  evt->print();
   std::unique_ptr<HepMCProduct> CMProduct(new HepMCProduct());
   if (evt) CMProduct->addHepMCData(evt);
   e.put(std::move(CMProduct), "unsmeared");
 
   unique_ptr<GenEventInfoProduct> genEventInfo(new GenEventInfoProduct(evt));
   e.put(std::move(genEventInfo));
 }

void BeamHaloProducer::setRandomEngine ( CLHEP::HepRandomEngine * v )

Definition at line 100 of file BeamHaloProducer.cc.

References _BeamHalo_randomEngine, and findQualityFiles::v.

                                                               {
   _BeamHalo_randomEngine = v;
 }

Member Data Documentation

float edm::BeamHaloProducer::EG_MAX_

private

Definition at line 56 of file BeamHaloProducer.h.

float edm::BeamHaloProducer::EG_MIN_

private

Definition at line 55 of file BeamHaloProducer.h.

HepMC::GenEvent* edm::BeamHaloProducer::evt

private

Definition at line 48 of file BeamHaloProducer.h.

Referenced by produce().

std::string edm::BeamHaloProducer::G3FNAME_

private

Definition at line 57 of file BeamHaloProducer.h.

int edm::BeamHaloProducer::GENMOD_

private

Definition at line 50 of file BeamHaloProducer.h.

bool edm::BeamHaloProducer::isInitialized_

private

Definition at line 59 of file BeamHaloProducer.h.

Referenced by beginLuminosityBlock().

int edm::BeamHaloProducer::IW_HAD_

private

Definition at line 54 of file BeamHaloProducer.h.

int edm::BeamHaloProducer::IW_MUO_

private

Definition at line 53 of file BeamHaloProducer.h.

int edm::BeamHaloProducer::LHC_B1_

private

Definition at line 51 of file BeamHaloProducer.h.

int edm::BeamHaloProducer::LHC_B2_

private

Definition at line 52 of file BeamHaloProducer.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation