#include <BeamHaloSource.h>

Inheritance diagram for edm::BeamHaloSource:

Public Member Functions
	BeamHaloSource (const ParameterSet &, const InputSourceDescription &)
	Constructor. More...

virtual	~BeamHaloSource ()
	Destructor. More...

Public Member Functions inherited from edm::GeneratedInputSource
	GeneratedInputSource (ParameterSet const &pset, InputSourceDescription const &desc)

virtual	~GeneratedInputSource ()

Public Member Functions inherited from edm::ConfigurableInputSource
	ConfigurableInputSource (ParameterSet const &pset, InputSourceDescription const &desc, bool realData=true)

EventNumber_t	event () const

unsigned int	eventCreationDelay () const

LuminosityBlockNumber_t	luminosityBlock () const

unsigned int	numberEventsInLumi () const

unsigned int	numberEventsInRun () const

unsigned int	numberEventsInThisLumi () const

unsigned int	numberEventsInThisRun () const

TimeValue_t	presentTime () const

RunNumber_t	run () const

unsigned int	timeBetweenEvents () const

virtual	~ConfigurableInputSource ()

Public Member Functions inherited from edm::InputSource
boost::shared_ptr < ActivityRegistry >	actReg () const
	Accessor for Activity Registry. More...

void	closeFile (boost::shared_ptr< FileBlock >)
	close current file More...

void	doBeginJob ()
	Called by framework at beginning of job. More...

void	doBeginLumi (LuminosityBlockPrincipal &lbp)
	Called by framework at beginning of lumi block. More...

void	doBeginRun (RunPrincipal &rp)
	Called by framework at beginning of run. More...

void	doEndJob ()
	Called by framework at end of job. More...

void	doEndLumi (LuminosityBlockPrincipal &lbp)
	Called by framework at end of lumi block. More...

void	doEndRun (RunPrincipal &rp)
	Called by framework at end of run. More...

void	doPostForkReacquireResources (boost::shared_ptr< multicore::MessageReceiverForSource >)

void	doPreForkReleaseResources ()
	Called by the framework before forking the process. More...

ProcessingController::ForwardState	forwardState () const

bool	goToEvent (EventID const &eventID)

	InputSource (ParameterSet const &, InputSourceDescription const &)
	Constructor. More...

void	issueReports (EventID const &eventID)
	issue an event report More...

LuminosityBlockNumber_t	luminosityBlock () const
	Accessor for current luminosity block number. More...

boost::shared_ptr < LuminosityBlockAuxiliary >	luminosityBlockAuxiliary () const
	Called by the framework to merge or insert lumi in principal cache. More...

int	markLumi ()
	Mark lumi as read. More...

int	markRun ()
	Mark run as read. More...

int	maxEvents () const

int	maxLuminosityBlocks () const

ModuleDescription const &	moduleDescription () const
	Accessor for 'module' description. More...

ItemType	nextItemType ()

bool const	primary () const
	Accessor for primary input source flag. More...

ProcessConfiguration const &	processConfiguration () const
	Accessor for Process Configuration. More...

std::string const &	processGUID () const
	Accessor for global process identifier. More...

ProcessHistoryID const &	processHistoryID () const
	Accessor for the input process history ID of the current run. More...

ProcessingMode	processingMode () const
	RunsLumisAndEvents (default), RunsAndLumis, or Runs. More...

boost::shared_ptr < ProductRegistry const >	productRegistry () const
	Accessor for product registry. More...

bool	randomAccess () const

void	readAndCacheLumi ()
	Read next luminosity block. More...

void	readAndCacheRun ()
	Read next run. More...

EventPrincipal *	readEvent (boost::shared_ptr< LuminosityBlockPrincipal > lbCache)

EventPrincipal *	readEvent (EventID const &)
	Read a specific event. More...

boost::shared_ptr< FileBlock >	readFile ()
	Read next file. More...

boost::shared_ptr < LuminosityBlockAuxiliary >	readLuminosityBlockAuxiliary ()
	Read next luminosity block Auxilary. More...

boost::shared_ptr< RunAuxiliary >	readRunAuxiliary ()
	Read next run Auxiliary. More...

void	registerProducts ()
	Register any produced products. More...

int	remainingEvents () const

int	remainingLuminosityBlocks () const

void	repeat ()
	Reset the remaining number of events/lumis to the maximum number. More...

ProcessingController::ReverseState	reverseState () const

void	rewind ()
	Begin again at the first event. More...

RunNumber_t	run () const
	Accessor for current run number. More...

boost::shared_ptr< RunAuxiliary >	runAuxiliary () const
	Called by the framework to merge or insert run in principal cache. More...

void	setLuminosityBlockNumber_t (LuminosityBlockNumber_t lb)
	Set the luminosity block ID. More...

void	setRunNumber (RunNumber_t r)
	Set the run number. More...

void	skipEvents (int offset)

Timestamp const &	timestamp () const
	Accessor for the current time, as seen by the input source. More...

void	wakeUp ()
	Wake up the input source. More...

virtual	~InputSource ()
	Destructor. More...

Private Member Functions
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 ()

virtual bool	produce (Event &e)

Private Attributes
float	EG_MAX_

float	EG_MIN_

HepMC::GenEvent *	evt

std::string	G3FNAME_

int	GENMOD_

int	IW_HAD_

int	IW_MUO_

int	LHC_B1_

int	LHC_B2_

Additional Inherited Members
Public Types inherited from edm::InputSource
enum	ItemType { IsInvalid, IsStop, IsFile, IsRun, IsLumi, IsEvent, IsRepeat }

enum	ProcessingMode { Runs, RunsAndLumis, RunsLumisAndEvents }

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::ConfigurableInputSource
static void	fillDescription (ParameterSetDescription &desc)

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

static void	fillDescription (ParameterSetDescription &desc)

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ConfigurableInputSource
void	reallyReadEvent ()

void	setEventNumber (EventNumber_t e)

void	setTime (TimeValue_t t)

Protected Member Functions inherited from edm::InputSource
void	decreaseRemainingEventsBy (int iSkipped)

EventPrincipal *const	eventPrincipalCache ()

boost::shared_ptr < LuminosityBlockPrincipal > const	luminosityBlockPrincipal () const

PrincipalCache const &	principalCache () const

PrincipalCache &	principalCache ()

ProductRegistry &	productRegistryUpdate () const

void	reset () const

void	resetLuminosityBlockAuxiliary () const

void	resetRunAuxiliary () const

boost::shared_ptr < RunPrincipal > const	runPrincipal () const

void	setLuminosityBlockAuxiliary (LuminosityBlockAuxiliary *lbp)

void	setLumiPrematurelyRead ()

void	setRunAuxiliary (RunAuxiliary *rp)

void	setRunPrematurelyRead ()

void	setTimestamp (Timestamp const &theTime)
	To set the current time, as seen by the input source. More...

ItemType	state () const

Detailed Description

Definition at line 22 of file BeamHaloSource.h.

Constructor & Destructor Documentation

BeamHaloSource::BeamHaloSource	(	const ParameterSet &	pset,
		const InputSourceDescription &	desc
	)

Constructor.

Definition at line 56 of file BeamHaloSource.cc.

References _BeamHalo_randomEngine, call_bh_set_parameters(), call_ki_bhg_init(), gather_cfg::cout, edm::RandomNumberGenerator::getEngine(), edm::ParameterSet::getUntrackedParameter(), edm::RandomNumberGenerator::mySeed(), and edm::ConfigurableInputSource::numberEventsInRun().

                                                                  :
     GeneratedInputSource(pset, desc), evt(0)
 {
 
    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]  = numberEventsInRun();
    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);
 
 
 // -- Seed for randomnumbers
     Service<RandomNumberGenerator> rng;
     _BeamHalo_randomEngine = &(rng->getEngine());
     long seed = (long)(rng->mySeed());
 
 
 // -- initialisation
    call_ki_bhg_init(seed);
 
 
   produces<HepMCProduct>();
   cout << "BeamHaloSource: starting event generation ... " << endl;
 
 }

BeamHaloSource::~BeamHaloSource ( )

virtual

Destructor.

Definition at line 50 of file BeamHaloSource.cc.

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

Member Function Documentation

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

private

Definition at line 147 of file BeamHaloSource.cc.

References BHSETPARAM.

Referenced by BeamHaloSource().

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

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

private

Definition at line 157 of file BeamHaloSource.cc.

References KI_BHG_FILL.

Referenced by produce().

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

bool BeamHaloSource::call_ki_bhg_init ( long & seed )

private

Definition at line 152 of file BeamHaloSource.cc.

References KI_BHG_INIT.

Referenced by BeamHaloSource().

                                                 {
     KI_BHG_INIT(seed);
     return true;
 }

bool BeamHaloSource::call_ki_bhg_stat ( int & iret )

private

Definition at line 162 of file BeamHaloSource.cc.

References KI_BHG_STAT.

                                                {
     KI_BHG_STAT(iret);
     return true;
 }

void BeamHaloSource::clear ( void )

private

Definition at line 100 of file BeamHaloSource.cc.

Referenced by python.Vispa.Views.BoxDecayView.BoxDecayView::closeEvent(), python.Vispa.Views.LineDecayView.LineDecayView::setDataObjects(), and python.Vispa.Views.BoxDecayView.BoxDecayView::updateContent().

101 {

102 }

bool BeamHaloSource::produce ( Event & e )

privatevirtual

Implements edm::ConfigurableInputSource.

Definition at line 104 of file BeamHaloSource.cc.

References call_ki_bhg_fill(), edm::ConfigurableInputSource::event(), evt, configurableAnalysis::GenParticle, L1TEmulatorMonitor_cff::p, edm::Event::put(), CommonMethods::weight(), ExpressReco_HICollisions_FallBack::weights, and wrapper.

                                       {
     // cout << "in produce " << endl;
 
   //        auto_ptr<HepMCProduct> bare_product(new HepMCProduct());
 
     // cout << "apres autoptr " << endl;
 
     int iret=0;
         float weight = 0;
         call_ki_bhg_fill(iret, weight);
 
     if( iret < 0 ) return false;
 
     // 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(event());
 
     HepMC::WeightContainer& weights = evt -> weights();
     weights.push_back(weight);
     //  evt->print();
   std::auto_ptr<HepMCProduct> CMProduct(new HepMCProduct());
   if (evt) CMProduct->addHepMCData(evt );
   e.put(CMProduct);
     
     return true;
 }

Member Data Documentation

float edm::BeamHaloSource::EG_MAX_

private

Definition at line 49 of file BeamHaloSource.h.

float edm::BeamHaloSource::EG_MIN_

private

Definition at line 48 of file BeamHaloSource.h.

HepMC::GenEvent* edm::BeamHaloSource::evt

private

Definition at line 41 of file BeamHaloSource.h.

Referenced by produce().

std::string edm::BeamHaloSource::G3FNAME_

private

Definition at line 50 of file BeamHaloSource.h.

int edm::BeamHaloSource::GENMOD_

private

Definition at line 43 of file BeamHaloSource.h.

int edm::BeamHaloSource::IW_HAD_

private

Definition at line 47 of file BeamHaloSource.h.

int edm::BeamHaloSource::IW_MUO_

private

Definition at line 46 of file BeamHaloSource.h.

int edm::BeamHaloSource::LHC_B1_

private

Definition at line 44 of file BeamHaloSource.h.

int edm::BeamHaloSource::LHC_B2_

private

Definition at line 45 of file BeamHaloSource.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation