#include <FlatRandomEGunSource.h>

Inheritance diagram for edm::FlatRandomEGunSource:

Public Member Functions
	FlatRandomEGunSource (const ParameterSet &, const InputSourceDescription &)

virtual	~FlatRandomEGunSource ()

Public Member Functions inherited from edm::BaseFlatGunSource
	BaseFlatGunSource (const ParameterSet &, const InputSourceDescription &)

void	beginRun (edm::Run &r, const edm::EventSetup &)

virtual	~BaseFlatGunSource ()

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 >, bool cleaningUpAfterException)
	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, bool cleaningUpAfterException)
	Called by framework at end of lumi block. More...

void	doEndRun (RunPrincipal &rp, bool cleaningUpAfterException)
	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	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...

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 (bool merge, HistoryAppender &historyAppender)
	Read next luminosity block. More...

void	readAndCacheRun (bool merge, HistoryAppender &historyAppender)
	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...

ProcessHistoryID const &	reducedProcessHistoryID () const

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...

virtual	~InputSource ()
	Destructor. More...

Protected Attributes
double	fMaxE

double	fMinE

Protected Attributes inherited from edm::BaseFlatGunSource
bool	fAddAntiParticle

HepMC::GenEvent *	fEvt

double	fMaxEta

double	fMaxPhi

double	fMinEta

double	fMinPhi

std::vector< int >	fPartIDs

ESHandle < HepPDT::ParticleDataTable >	fPDGTable

CLHEP::HepRandomEngine &	fRandomEngine

CLHEP::RandFlat *	fRandomGenerator

int	fVerbosity

Private Member Functions
virtual bool	produce (Event &e)

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)

static void	prevalidate (ConfigurationDescriptions &)

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 *	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	setRunAuxiliary (RunAuxiliary *rp)

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 15 of file FlatRandomEGunSource.h.

Constructor & Destructor Documentation

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

Definition at line 19 of file FlatRandomEGunSource.cc.

References gather_cfg::cout, fMaxE, fMinE, edm::ParameterSet::getUntrackedParameter(), and edm::InputSource::remainingEvents().

                                                                                :
    BaseFlatGunSource(pset, desc)
 {
 
    ParameterSet defpset ;
    // ParameterSet pgun_params = pset.getParameter<ParameterSet>("PGunParameters") ;
    ParameterSet pgun_params = 
       pset.getUntrackedParameter<ParameterSet>("PGunParameters",defpset) ;
   
    // doesn't seem necessary to check if pset is empty - if this
    // is the case, default values will be taken for params
    fMinE = pgun_params.getUntrackedParameter<double>("MinE",0.99);
    fMaxE = pgun_params.getUntrackedParameter<double>("MaxE",1.01); 
   
    produces<HepMCProduct>();
 
    cout << "Internal FlatRandomEGun is initialzed" << endl ;
    cout << "It is going to generate " << remainingEvents() << "events" << endl ;
    
 }

FlatRandomEGunSource::~FlatRandomEGunSource ( )

virtual

Definition at line 41 of file FlatRandomEGunSource.cc.

 {
    // no need to cleanup fEvt since it's done in HepMCProduct
 }

Member Function Documentation

bool FlatRandomEGunSource::produce ( Event & e )

privatevirtual

Implements edm::ConfigurableInputSource.

Definition at line 46 of file FlatRandomEGunSource.cc.

 {
 
    if ( fVerbosity > 0 )
    {
       cout << " FlatRandomEGunSource : Begin New Event Generation" << endl ; 
    }
    
    // event loop (well, another step in it...)
           
    // no need to clean up GenEvent memory - done in HepMCProduct
 
    // here re-create fEvt (memory)
    //
    fEvt = new HepMC::GenEvent() ;
    
    // now actualy, cook up the event from PDGTable and gun parameters
    //
 
    // 1st, primary vertex
    //
    HepMC::GenVertex* Vtx = new HepMC::GenVertex( HepMC::FourVector(0.,0.,0.));
    
    // loop over particles
    //
    int barcode = 1;
    for (unsigned int ip=0; ip<fPartIDs.size(); ip++)
    {
        double energy = fRandomGenerator->fire(fMinE, fMaxE) ;
        double eta    = fRandomGenerator->fire(fMinEta, fMaxEta) ;
        double phi    = fRandomGenerator->fire(fMinPhi, fMaxPhi) ;
        int PartID = fPartIDs[ip] ;
        const HepPDT::ParticleData* 
           PData = fPDGTable->particle(HepPDT::ParticleID(abs(PartID))) ;
        double mass   = PData->mass().value() ;
        double mom2   = energy*energy - mass*mass ;
        double mom    = 0. ;
        if (mom2 > 0.) 
        {
           mom = sqrt(mom2) ;
        }
        else
        {
           mom = 0. ;
        }
        double theta  = 2.*atan(exp(-eta)) ;
        double px     = mom*sin(theta)*cos(phi) ;
        double py     = mom*sin(theta)*sin(phi) ;
        double pz     = mom*cos(theta) ;
 
        HepMC::FourVector p(px,py,pz,energy) ;
        HepMC::GenParticle* Part = 
            new HepMC::GenParticle(p,PartID,1);
        Part->suggest_barcode( barcode ) ;
        barcode++ ;
        Vtx->add_particle_out(Part);
        
        if ( fAddAntiParticle )
        {
           HepMC::FourVector ap(-px,-py,-pz,energy) ;
       int APartID = -PartID ;
       if ( PartID == 22 || PartID == 23 )
       {
          APartID = PartID ;
       }
       HepMC::GenParticle* APart =
          new HepMC::GenParticle(ap,APartID,1);
       APart->suggest_barcode( barcode ) ;
       barcode++ ;
       Vtx->add_particle_out(APart) ;
        }
        
    }
    fEvt->add_vertex(Vtx) ;
    fEvt->set_event_number(event()) ;
    fEvt->set_signal_process_id(20) ;  
    
    
    if ( fVerbosity > 0 )
    {
       fEvt->print() ;  
    }  
 
    auto_ptr<HepMCProduct> BProduct(new HepMCProduct()) ;
    BProduct->addHepMCData( fEvt );
    e.put(BProduct);
     
    if ( fVerbosity > 0 )
    {
       // for testing purpose only
       //fEvt->print() ;  // for some strange reason, it prints NO event info
       // after it's been put into edm::Event...
       cout << " FlatRandomEGunSource : Event Generation Done " << endl;
    }
       
    return true;
 }

Member Data Documentation

double edm::FlatRandomEGunSource::fMaxE

protected

Definition at line 31 of file FlatRandomEGunSource.h.

Referenced by FlatRandomEGunSource(), and produce().

double edm::FlatRandomEGunSource::fMinE

protected

Definition at line 30 of file FlatRandomEGunSource.h.

Referenced by FlatRandomEGunSource(), and produce().

Public Member Functions

Protected Attributes

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation