Inheritance diagram for LHEProducer:

Public Member Functions
	LHEProducer (const edm::ParameterSet &params)

virtual	~LHEProducer ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Protected Member Functions
virtual void	beginJob ()

virtual bool	beginRun (edm::Run &run, const edm::EventSetup &es)

virtual void	endJob ()

virtual bool	endRun (edm::Run &run, const edm::EventSetup &es)

virtual bool	filter (edm::Event &event, const edm::EventSetup &es)

Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Private Member Functions
double	matching (const HepMC::GenEvent *event, bool shower) const

void	onBeforeHadronisation ()

void	onInit ()

bool	showeredEvent (const boost::shared_ptr< HepMC::GenEvent > &event)

Private Attributes
BranchingRatios	branchingRatios

unsigned int	eventsToPrint

double	extCrossSect

double	extFilterEff

std::auto_ptr< Hadronisation >	hadronisation

unsigned int	index

std::auto_ptr< JetMatching >	jetMatching

bool	matchingDone

bool	matchSummary

boost::shared_ptr< LHEEvent >	partonLevel

std::vector< int >	removeResonances

boost::shared_ptr< LHERunInfo >	runInfo

double	weight

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

typedef WorkerT< EDFilter >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 33 of file LHEProducer.cc.

Constructor & Destructor Documentation

LHEProducer::LHEProducer ( const edm::ParameterSet & params )

explicit

Definition at line 69 of file LHEProducer.cc.

References branchingRatios, SurfaceDeformationFactory::create(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hadronisation, i, jetMatching, matchSummary, onBeforeHadronisation(), onInit(), removeResonances, lhef::BranchingRatios::set(), and showeredEvent().

                                                       :
     eventsToPrint(params.getUntrackedParameter<unsigned int>("eventsToPrint", 0)),
     extCrossSect(params.getUntrackedParameter<double>("crossSection", -1.0)),
     extFilterEff(params.getUntrackedParameter<double>("filterEfficiency", -1.0)),
     matchSummary(false)
 {
     hadronisation = Hadronisation::create(
         params.getParameter<edm::ParameterSet>("hadronisation"));
 
     if (params.exists("removeResonances"))
         removeResonances =
             params.getParameter<std::vector<int> >(
                             "removeResonances");
 
     std::set<std::string> matchingCapabilities;
     if (params.exists("jetMatching")) {
         edm::ParameterSet jetParams =
             params.getUntrackedParameter<edm::ParameterSet>(
                                 "jetMatching");
         jetMatching = JetMatching::create(jetParams);
 
         matchingCapabilities = jetMatching->capabilities();
         hadronisation->matchingCapabilities(matchingCapabilities);
     }
 
     produces<edm::HepMCProduct>();
     produces<GenEventInfoProduct>();
     produces<GenRunInfoProduct, edm::InRun>();
 
     if (jetMatching.get()) {
         if (params.getUntrackedParameter<bool>(
                     "preferShowerVetoCallback", true))
             hadronisation->onShoweredEvent().connect(
                 sigc::mem_fun(*this,
                               &LHEProducer::showeredEvent));
         hadronisation->onInit().connect(
                 sigc::mem_fun(*this, &LHEProducer::onInit));
         hadronisation->onBeforeHadronisation().connect(
             sigc::mem_fun(*this,
                           &LHEProducer::onBeforeHadronisation));
 
         matchSummary = matchingCapabilities.count("matchSummary");
         if (matchSummary) {
             produces< std::vector<double> >("matchDeltaR");
             produces< std::vector<double> >("matchDeltaPRel");
         }
     }
 
     // force total branching ratio for QCD/QED to 1
     for(int i = 0; i < 6; i++)
         branchingRatios.set(i);
     for(int i = 9; i < 23; i++)
         branchingRatios.set(i);
 }

LHEProducer::~LHEProducer ( )

virtual

Definition at line 124 of file LHEProducer.cc.

125 {

126 }

Member Function Documentation

void LHEProducer::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDFilter.

Definition at line 128 of file LHEProducer.cc.

References hadronisation.

 {
     hadronisation->init();
 }

bool LHEProducer::beginRun	(	edm::Run &	run,
		const edm::EventSetup &	es
	)

protectedvirtual

Reimplemented from edm::EDFilter.

Definition at line 139 of file LHEProducer.cc.

References edm::Run::getByLabel(), index, and runInfo.

 {
     edm::Handle<LHERunInfoProduct> product;
     run.getByLabel("source", product);
 
     runInfo.reset(new LHERunInfo(*product));
     index = 0;
 
     return true;
 }

void LHEProducer::endJob ( void )

protectedvirtual

Reimplemented from edm::EDFilter.

Definition at line 133 of file LHEProducer.cc.

References hadronisation, and jetMatching.

 {
     hadronisation.reset();
     jetMatching.reset();
 }

bool LHEProducer::endRun	(	edm::Run &	run,
		const edm::EventSetup &	es
	)

protectedvirtual

Reimplemented from edm::EDFilter.

Definition at line 150 of file LHEProducer.cc.

References PYTHIA6_MinBias_2360GeV_cff_py_GEN_FASTSIM::crossSection, lhef::LHERunInfo::XSec::error, extCrossSect, extFilterEff, hadronisation, edm::Run::put(), runInfo, and lhef::LHERunInfo::XSec::value.

 {
     hadronisation->statistics();
 
     LHERunInfo::XSec crossSection;
     if (runInfo) {
         crossSection = runInfo->xsec();
         runInfo->statistics();
     }
 
     std::auto_ptr<GenRunInfoProduct> runInfo(new GenRunInfoProduct);
 
     runInfo->setInternalXSec(
             GenRunInfoProduct::XSec(crossSection.value,
                                     crossSection.error));
     runInfo->setExternalXSecLO(extCrossSect);
     runInfo->setFilterEfficiency(extFilterEff);
 
     run.put(runInfo);
 
     runInfo.reset();
 
     return true;
 }

bool LHEProducer::filter	(	edm::Event &	event,
		const edm::EventSetup &	es
	)

protectedvirtual

Implements edm::EDFilter.

Definition at line 175 of file LHEProducer.cc.

References beamvalidation::br, branchingRatios, eventsToPrint, lhef::BranchingRatios::getFactor(), hadronisation, index, info, jetMatching, matching(), matchingDone, matchSummary, partonLevel, removeResonances, query::result, runInfo, and weight.

 {
     std::auto_ptr<edm::HepMCProduct> result(new edm::HepMCProduct);
 
     edm::Handle<LHEEventProduct> product;
     event.getByLabel("source", product);
 
     partonLevel.reset(new LHEEvent(runInfo, *product));
     if (!removeResonances.empty())
         partonLevel->removeResonances(removeResonances);
 
     if (product->pdf())
         partonLevel->setPDF(
             std::auto_ptr<LHEEvent::PDF>(
                 new LHEEvent::PDF(*product->pdf())));
 
     hadronisation->setEvent(partonLevel);
 
     double br = branchingRatios.getFactor(hadronisation.get(),
                                           partonLevel);
 
     matchingDone = false;
     weight = 1.0;
     std::auto_ptr<HepMC::GenEvent> hadronLevel(hadronisation->hadronize());
 
     if (!hadronLevel.get()) {
         if (matchingDone) {
             if (weight == 0.0)
                 partonLevel->count(LHERunInfo::kSelected, br);
             else
                 partonLevel->count(LHERunInfo::kKilled,
                                    br, weight);
         } else
             partonLevel->count(LHERunInfo::kTried, br);
     }
 
     if (!matchingDone && jetMatching.get() && hadronLevel.get())
         weight = matching(hadronLevel.get(), false);
 
     if (weight == 0.0) {
         edm::LogInfo("Generator|LHEInterface")
             << "Event got rejected by the "
                "jet matching." << std::endl;
 
         if (hadronLevel.get()) {
             partonLevel->count(LHERunInfo::kSelected);
             hadronLevel.reset();
         }
     }
 
     if (!hadronLevel.get()) {
         event.put(result);
         std::auto_ptr<GenEventInfoProduct> info(
                         new GenEventInfoProduct);
         event.put(info);
         return false;
     }
 
     partonLevel->count(LHERunInfo::kAccepted, br, weight);
 
     hadronLevel->set_event_number(++index);
 
     if (eventsToPrint) {
         eventsToPrint--;
         hadronLevel->print();
     }
 
     std::auto_ptr<GenEventInfoProduct> info(
                 new GenEventInfoProduct(hadronLevel.get()));
     result->addHepMCData(hadronLevel.release());
     event.put(result);
     event.put(info);
 
     if (jetMatching.get() && matchSummary) {
         std::auto_ptr< std::vector<double> > matchDeltaR(
                         new std::vector<double>);
         std::auto_ptr< std::vector<double> > matchDeltaPRel(
                         new std::vector<double>);
 
         typedef std::vector<JetMatching::JetPartonMatch> Matches;
         Matches matches = jetMatching->getMatchSummary();
 
         for(Matches::const_iterator iter = matches.begin();
             iter != matches.end(); ++iter) {
             if (!iter->isMatch())
                 continue;
 
             matchDeltaR->push_back(iter->delta);
             matchDeltaPRel->push_back(iter->jet.rho() /
                                       iter->parton.rho() - 1.0);
         }
 
         event.put(matchDeltaR, "matchDeltaR");
         event.put(matchDeltaPRel, "matchDeltaPRel");
     }
 
     return true;
 }

double LHEProducer::matching	(	const HepMC::GenEvent *	event,
		bool	shower
	)		const

private

Definition at line 274 of file LHEProducer.cc.

References event(), jetMatching, and partonLevel.

Referenced by filter(), and showeredEvent().

 {
     if (!jetMatching.get())
         return 1.0;
 
     return jetMatching->match(partonLevel->asHepMCEvent().get(),
                               event, shower);
 }

void LHEProducer::onBeforeHadronisation ( )

private

Definition at line 295 of file LHEProducer.cc.

References jetMatching, and partonLevel.

Referenced by LHEProducer().

 {
     jetMatching->beforeHadronisation(partonLevel);
 }

void LHEProducer::onInit ( )

private

Definition at line 290 of file LHEProducer.cc.

References jetMatching, and runInfo.

Referenced by LHEProducer().

 {
     jetMatching->init(runInfo);
 }

bool LHEProducer::showeredEvent ( const boost::shared_ptr< HepMC::GenEvent > & event )

private

Definition at line 283 of file LHEProducer.cc.

References matching(), matchingDone, and weight.

Referenced by LHEProducer().

 {
     weight = matching(event.get(), true);
     matchingDone = true;
     return weight == 0.0;
 }

Member Data Documentation

BranchingRatios LHEProducer::branchingRatios

private

Definition at line 66 of file LHEProducer.cc.

Referenced by filter(), and LHEProducer().

unsigned int LHEProducer::eventsToPrint

private

Definition at line 52 of file LHEProducer.cc.

Referenced by filter().

double LHEProducer::extCrossSect

private

Definition at line 57 of file LHEProducer.cc.

Referenced by endRun().

double LHEProducer::extFilterEff

private

Definition at line 58 of file LHEProducer.cc.

Referenced by endRun().

std::auto_ptr<Hadronisation> LHEProducer::hadronisation

private

Definition at line 54 of file LHEProducer.cc.

Referenced by beginJob(), endJob(), endRun(), filter(), and LHEProducer().

unsigned int LHEProducer::index

private

Definition at line 63 of file LHEProducer.cc.

Referenced by beginRun(), and filter().

std::auto_ptr<JetMatching> LHEProducer::jetMatching

private

Definition at line 55 of file LHEProducer.cc.

Referenced by endJob(), filter(), LHEProducer(), matching(), onBeforeHadronisation(), and onInit().

bool LHEProducer::matchingDone

private

Definition at line 64 of file LHEProducer.cc.

Referenced by filter(), and showeredEvent().

bool LHEProducer::matchSummary

private

Definition at line 59 of file LHEProducer.cc.

Referenced by filter(), and LHEProducer().

boost::shared_ptr<LHEEvent> LHEProducer::partonLevel

private

Definition at line 61 of file LHEProducer.cc.

Referenced by filter(), matching(), and onBeforeHadronisation().

std::vector<int> LHEProducer::removeResonances

private

Definition at line 53 of file LHEProducer.cc.

Referenced by filter(), and LHEProducer().

boost::shared_ptr<LHERunInfo> LHEProducer::runInfo

private

Definition at line 62 of file LHEProducer.cc.

Referenced by beginRun(), endRun(), filter(), and onInit().

double LHEProducer::weight

private

Definition at line 65 of file LHEProducer.cc.

Referenced by filter(), and showeredEvent().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation