ThePEGHadronizer Class Reference

Inheritance diagram for ThePEGHadronizer:

Public Member Functions
const char *	classname () const
bool	decay ()
bool	declareSpecialSettings (const std::vector< std::string >)
bool	declareStableParticles (const std::vector< int > &pdgIds)
void	finalizeEvent ()
bool	generatePartonsAndHadronize ()
bool	hadronize ()
bool	initializeForExternalPartons ()
bool	initializeForInternalPartons ()
bool	readSettings (int)
bool	residualDecay ()
void	statistics ()
	ThePEGHadronizer (const edm::ParameterSet &params)
virtual	~ThePEGHadronizer ()
Public Member Functions inherited from ThePEGInterface
void	setPEGRandomEngine (CLHEP::HepRandomEngine *)
	ThePEGInterface (const edm::ParameterSet &params)
virtual	~ThePEGInterface ()
Public Member Functions inherited from gen::BaseHadronizer
	BaseHadronizer (edm::ParameterSet const &ps)
edm::Event &	getEDMEvent () const
HepMC::GenEvent *	getGenEvent ()
GenEventInfoProduct *	getGenEventInfo ()
GenRunInfoProduct &	getGenRunInfo ()
const boost::shared_ptr < lhef::LHERunInfo > &	getLHERunInfo () const
void	resetEvent (HepMC::GenEvent *event)
void	resetEventInfo (GenEventInfoProduct *eventInfo)
virtual bool	select (HepMC::GenEvent *) const
void	setEDMEvent (edm::Event &event)
void	setLHEEvent (lhef::LHEEvent *event)
void	setLHERunInfo (lhef::LHERunInfo *runInfo)
void	setRandomEngine (CLHEP::HepRandomEngine *v)
std::vector< std::string > const &	sharedResources () const
	~BaseHadronizer ()

Private Member Functions
virtual void	doSetRandomEngine (CLHEP::HepRandomEngine *v) override

Private Attributes
unsigned int	eventsToPrint
const std::string	handlerDirectory_
boost::shared_ptr< lhef::LHEProxy >	proxy_
ThePEG::EventPtr	thepegEvent

Additional Inherited Members
Protected Member Functions inherited from ThePEGInterface
std::string	dataFile (const std::string &fileName) const
std::string	dataFile (const edm::ParameterSet &pset, const std::string &paramName) const
void	flushRandomNumberGenerator ()
void	initGenerator ()
void	initRepository (const edm::ParameterSet &params) const
Protected Member Functions inherited from gen::BaseHadronizer
std::auto_ptr< HepMC::GenEvent > &	event ()
std::auto_ptr < GenEventInfoProduct > &	eventInfo ()
lhef::LHEEvent *	lheEvent ()
lhef::LHERunInfo *	lheRunInfo ()
GenRunInfoProduct &	runInfo ()
Static Protected Member Functions inherited from ThePEGInterface
static void	clearAuxiliary (HepMC::GenEvent *hepmc, HepMC::PdfInfo *pdf)
static std::auto_ptr < HepMC::GenEvent >	convert (const ThePEG::EventPtr &event)
static void	fillAuxiliary (HepMC::GenEvent *hepmc, HepMC::PdfInfo *pdf, const ThePEG::EventPtr &event)
static double	pthat (const ThePEG::EventPtr &event)
Protected Attributes inherited from ThePEGInterface
ThePEG::EGPtr	eg_
std::auto_ptr < HepMC::IO_BaseClass >	iobc_

Detailed Description

Definition at line 30 of file ThePEGHadronizer.cc.

Constructor & Destructor Documentation

ThePEGHadronizer::ThePEGHadronizer

(

const edm::ParameterSet &

params

)

Definition at line 63 of file ThePEGHadronizer.cc.

References ThePEGInterface::initRepository().

                                                              :
    ThePEGInterface(pset),
    BaseHadronizer(pset),
    eventsToPrint(pset.getUntrackedParameter<unsigned int>("eventsToPrint", 0)),
    handlerDirectory_(pset.getParameter<std::string>("eventHandlers"))
{  
    initRepository(pset);

}

ThePEGHadronizer::~ThePEGHadronizer ( )

virtual

Definition at line 73 of file ThePEGHadronizer.cc.

{
}

Member Function Documentation

const char* ThePEGHadronizer::classname ( ) const

inline

Definition at line 49 of file ThePEGHadronizer.cc.

{ return "ThePEGHadronizer"; }

bool ThePEGHadronizer::decay

(

)

Definition at line 219 of file ThePEGHadronizer.cc.

{
    return true;
}

bool ThePEGHadronizer::declareSpecialSettings ( const std::vector< std::string >  )

inline

Definition at line 39 of file ThePEGHadronizer.cc.

{ return true; }

bool ThePEGHadronizer::declareStableParticles

(

const std::vector< int > &

pdgIds

)

Definition at line 101 of file ThePEGHadronizer.cc.

{
    return false;
}

virtual void ThePEGHadronizer::doSetRandomEngine ( CLHEP::HepRandomEngine *  v )

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 53 of file ThePEGHadronizer.cc.

References ThePEGInterface::setPEGRandomEngine().

{ setPEGRandomEngine(v); }

void ThePEGHadronizer::finalizeEvent

(

)

Definition at line 197 of file ThePEGHadronizer.cc.

References ThePEGInterface::clearAuxiliary(), gen::BaseHadronizer::event(), gen::BaseHadronizer::eventInfo(), eventsToPrint, ThePEGInterface::fillAuxiliary(), ThePEGInterface::iobc_, ThePEGInterface::pthat(), and thepegEvent.

{
    HepMC::PdfInfo pdf;
    clearAuxiliary(event().get(), &pdf);
    fillAuxiliary(event().get(), &pdf, thepegEvent);
    event()->set_pdf_info(pdf);

    eventInfo().reset(new GenEventInfoProduct(event().get()));
    eventInfo()->setBinningValues(
            std::vector<double>(1, pthat(thepegEvent)));

    if (eventsToPrint) {
        eventsToPrint--;
        event()->print();
    }

    if (iobc_.get())
        iobc_->write_event(event().get());

    edm::LogInfo("Generator|ThePEGHadronizer") << "Event produced";
}

bool ThePEGHadronizer::generatePartonsAndHadronize

(

)

Definition at line 113 of file ThePEGHadronizer.cc.

References ThePEGInterface::convert(), ThePEGInterface::eg_, gen::BaseHadronizer::event(), cppFunctionSkipper::exception, ThePEGInterface::flushRandomNumberGenerator(), and thepegEvent.

{
    edm::LogInfo("Generator|ThePEGHadronizer") << "Start production";

    flushRandomNumberGenerator();

        try {
                thepegEvent = eg_->shoot();
        } catch (std::exception& exc) {
                edm::LogWarning("Generator|ThePEGHadronizer") << "EGPtr::shoot() thrown an exception, event skipped: " << exc.what();
                return false;
        } catch (...) {
                edm::LogWarning("Generator|ThePEGHadronizer") << "EGPtr::shoot() thrown an unknown exception, event skipped";
                return false;
        }        
        
    if (!thepegEvent) {
        edm::LogWarning("Generator|ThePEGHadronizer") << "thepegEvent not initialized";
        return false;
    }

    event() = convert(thepegEvent);
    if (!event().get()) {
        edm::LogWarning("Generator|ThePEGHadronizer") << "genEvent not initialized";
        return false;
    }

    return true;
}

bool ThePEGHadronizer::hadronize

(

)

Definition at line 143 of file ThePEGHadronizer.cc.

References ThePEGInterface::clearAuxiliary(), ThePEGInterface::convert(), lhef::LHEEvent::count(), ThePEGInterface::eg_, gen::BaseHadronizer::event(), cppFunctionSkipper::exception, ThePEGInterface::fillAuxiliary(), lhef::LHEEvent::fillPdfInfo(), ThePEGInterface::flushRandomNumberGenerator(), lhef::LHEEvent::getHEPEUP(), lhef::LHERunInfo::kAccepted, lhef::LHERunInfo::kSelected, gen::BaseHadronizer::lheEvent(), proxy_, thepegEvent, and lhef::HEPEUP::XWGTUP.

{
    edm::LogInfo("Generator|ThePEGHadronizer") << "Start production";

    flushRandomNumberGenerator();
    
    //need to copy lhe event here unfortunately because of interface mismatch
    proxy_->loadEvent(boost::shared_ptr<lhef::LHEEvent>(new lhef::LHEEvent(*lheEvent())));

        //dummy for now
        double mergeweight = 1.0;
        
    try {
        thepegEvent = eg_->shoot();
    } catch (std::exception& exc) {
        edm::LogWarning("Generator|ThePEGHadronizer") << "EGPtr::shoot() thrown an exception, event skipped: " << exc.what();
                lheEvent()->count( lhef::LHERunInfo::kSelected, 1.0, mergeweight );
        return false;
    } catch (...) {
        edm::LogWarning("Generator|ThePEGHadronizer") << "EGPtr::shoot() thrown an unknown exception, event skipped";
                lheEvent()->count( lhef::LHERunInfo::kSelected, 1.0, mergeweight );
        return false;
    }

    if (!thepegEvent) {
        edm::LogWarning("Generator|ThePEGHadronizer") << "thepegEvent not initialized";
                lheEvent()->count( lhef::LHERunInfo::kSelected, 1.0, mergeweight );
        return false;
    }

    event() = convert(thepegEvent);
    if (!event().get()) {
        edm::LogWarning("Generator|ThePEGHadronizer") << "genEvent not initialized";
                lheEvent()->count( lhef::LHERunInfo::kSelected, 1.0, mergeweight );
        return false;
    }
    
    //Fill LHE weight (since it's not otherwise propagated)
    event()->weights()[0] *= lheEvent()->getHEPEUP()->XWGTUP;

    HepMC::PdfInfo pdf;
    clearAuxiliary(event().get(), &pdf);
    lheEvent()->fillPdfInfo(&pdf);
    fillAuxiliary(event().get(), &pdf, thepegEvent);
    event()->set_pdf_info(pdf);

        // update LHE matching statistics
        //
        lheEvent()->count( lhef::LHERunInfo::kAccepted, 1.0, mergeweight );        
        
    return true;

}

bool ThePEGHadronizer::initializeForExternalPartons

(

)

Definition at line 83 of file ThePEGHadronizer.cc.

References lhef::LHEProxy::create(), gen::BaseHadronizer::getLHERunInfo(), handlerDirectory_, ThePEGInterface::initGenerator(), proxy_, and contentValuesCheck::ss.

{
    proxy_ = lhef::LHEProxy::create();  
  
    std::ostringstream ss;
    ss << proxy_->getID();

    std::ostringstream logstream;
    ThePEG::Repository::exec("set " + handlerDirectory_ +
                             "/LHEReader:ProxyID " + ss.str(), logstream);
    edm::LogInfo("Generator|LHEInterface") << logstream.str();  
    
    proxy_->loadRunInfo(getLHERunInfo());   
    
    initGenerator();
    return true;
}

bool ThePEGHadronizer::initializeForInternalPartons

(

)

Definition at line 77 of file ThePEGHadronizer.cc.

References ThePEGInterface::initGenerator().

{
    initGenerator();
    return true;
}

bool ThePEGHadronizer::readSettings ( int  )

inline

Definition at line 35 of file ThePEGHadronizer.cc.

{ return true; }

bool ThePEGHadronizer::residualDecay

(

)

Definition at line 224 of file ThePEGHadronizer.cc.

{
    return true;
}

void ThePEGHadronizer::statistics

(

)

Definition at line 106 of file ThePEGHadronizer.cc.

References ThePEGInterface::eg_, gen::BaseHadronizer::runInfo(), and GenRunInfoProduct::setInternalXSec().

{
    runInfo().setInternalXSec(GenRunInfoProduct::XSec(
        eg_->integratedXSec() / ThePEG::picobarn,
        eg_->integratedXSecErr() / ThePEG::picobarn));
}

Member Data Documentation

unsigned int ThePEGHadronizer::eventsToPrint

private

Definition at line 55 of file ThePEGHadronizer.cc.

Referenced by finalizeEvent().

const std::string ThePEGHadronizer::handlerDirectory_

private

Definition at line 60 of file ThePEGHadronizer.cc.

Referenced by initializeForExternalPartons().

boost::shared_ptr<lhef::LHEProxy> ThePEGHadronizer::proxy_

private

Definition at line 59 of file ThePEGHadronizer.cc.

Referenced by hadronize(), and initializeForExternalPartons().

ThePEG::EventPtr ThePEGHadronizer::thepegEvent

private

Definition at line 57 of file ThePEGHadronizer.cc.

Referenced by finalizeEvent(), generatePartonsAndHadronize(), and hadronize().