#include <JetMatchingHook.h>

Inheritance diagram for JetMatchingHook:

Public Member Functions
virtual void	beforeHadronization (lhef::LHEEvent *lhee)

bool	canVetoPartonLevel () override

bool	doVetoPartonLevel (const Pythia8::Event &event) override

virtual void	init (lhef::LHERunInfo *runInfo)

bool	initAfterBeams () override

	JetMatchingHook (const edm::ParameterSet &, Pythia8::Info *)

void	resetMatchingStatus ()

void	setEventNumber (int ievt)

	~JetMatchingHook () override

Protected Member Functions
	JetMatchingHook ()

void	setLHEEvent (lhef::LHEEvent *lhee)

void	setLHERunInfo (lhef::LHERunInfo *lheri)

Protected Attributes
lhef::LHEEvent *	fEventBlock

int	fEventNumber

Pythia8::Info *	fInfoPtr

bool	fIsInitialized

Py8toJetInput *	fJetInputFill

gen::JetMatching *	fJetMatching

lhef::LHERunInfo *	fRunBlock

Detailed Description

Definition at line 18 of file JetMatchingHook.h.

Constructor & Destructor Documentation

JetMatchingHook::JetMatchingHook	(	const edm::ParameterSet &	,
		Pythia8::Info *
	)

JetMatchingHook::~JetMatchingHook ( )

override

Definition at line 56 of file JetMatchingHook.cc.

References fJetMatching.

                                   {
   if (fJetMatching)
     delete fJetMatching;
 }

JetMatchingHook::JetMatchingHook ( )

inlineprotected

Definition at line 53 of file JetMatchingHook.h.

53 : UserHooks() {}

Member Function Documentation

void JetMatchingHook::beforeHadronization ( lhef::LHEEvent * lhee )

virtual

Definition at line 72 of file JetMatchingHook.cc.

References gen::JetMatching::beforeHadronisation(), gen::JetMatching::beforeHadronisationExec(), fJetMatching, and setLHEEvent().

Referenced by resetMatchingStatus().

                                                             {
   setLHEEvent(lhee);
   fJetMatching->beforeHadronisation(lhee);
 
   // here we'll have to adjust, if needed, for "massless" particles
   // from earlier Madgraph version(s)
   // also, we'll have to setup elements of the Py6 fortran array
   // VINT(357), VINT(358), VINT(360) and VINT(390)
   // if ( fJetMatching->getMatchingScheme() == "Madgraph" )
   // {
   //
   // }
 
   fJetMatching->beforeHadronisationExec();
 
   return;
 }

bool JetMatchingHook::canVetoPartonLevel ( )

inlineoverride

Definition at line 30 of file JetMatchingHook.h.

References doVetoPartonLevel().

30 { return true; }

bool JetMatchingHook::doVetoPartonLevel ( const Pythia8::Event & event )

override

Definition at line 90 of file JetMatchingHook.cc.

References fEventBlock, Py8toJetInput::fillJetAlgoInput(), fJetInputFill, fJetMatching, gen::JetMatching::getPartonList(), hepeup_, gen::JetMatching::isMatchingDone(), gen::JetMatching::match(), and HEPEUP_::nup.

Referenced by canVetoPartonLevel().

 {
   // event.list();
 
   // extract "hardest" event - the output will go into workEvent,
   // which is a data mamber of base class UserHooks
   //
   subEvent(event, true);
 
   if (!hepeup_.nup || fJetMatching->isMatchingDone()) {
     return true;
   }
 
   //
   // bool jmtch = fJetMatching->match( 0, 0, true ); // true if veto-ed, false if accepted (not veto-ed)
   std::vector<fastjet::PseudoJet> jetInput =
       fJetInputFill->fillJetAlgoInput(event, workEvent, fEventBlock, fJetMatching->getPartonList());
   bool jmtch = fJetMatching->match(fEventBlock, &jetInput);
   if (jmtch) {
     return true;
   }
 
   // Do not veto events that got this far
   //
   return false;
 }

void JetMatchingHook::init ( lhef::LHERunInfo * runInfo )

virtual

Definition at line 61 of file JetMatchingHook.cc.

References ALCARECOTkAlBeamHalo_cff::etaMax, Exception, fJetInputFill, fJetMatching, fRunBlock, gen::JetMatching::getJetEtaMax(), gen::JetMatching::init(), Py8toJetInput::setJetEtaMax(), and setLHERunInfo().

Referenced by setEventNumber().

                                                   {
   setLHERunInfo(runInfo);
   if (!fRunBlock) {
     throw cms::Exception("JetMatching") << "Invalid RunInfo" << std::endl;
   }
   fJetMatching->init(runInfo);
   double etaMax = fJetMatching->getJetEtaMax();
   fJetInputFill->setJetEtaMax(etaMax);
   return;
 }

bool JetMatchingHook::initAfterBeams ( )

inlineoverride

Definition at line 39 of file JetMatchingHook.h.

References fIsInitialized, fJetMatching, and gen::JetMatching::initAfterBeams().

                                  {
     if (fIsInitialized)
       return true;
     fJetMatching->initAfterBeams();
     fIsInitialized = true;
     return true;
   }

void JetMatchingHook::resetMatchingStatus ( )

inline

Definition at line 46 of file JetMatchingHook.h.

References beforeHadronization(), fJetMatching, and gen::JetMatching::resetMatchingStatus().

                              {
     fJetMatching->resetMatchingStatus();
     return;
   }

void JetMatchingHook::setEventNumber ( int ievt )

inline

Definition at line 33 of file JetMatchingHook.h.

References fEventNumber, and init().

                                 {
     fEventNumber = ievt;
     return;
   }

void JetMatchingHook::setLHEEvent ( lhef::LHEEvent * lhee )

inlineprotected

Definition at line 63 of file JetMatchingHook.h.

References fEventBlock, lhef::CommonBlocks::fillHEPEUP(), and lhef::LHEEvent::getHEPEUP().

Referenced by beforeHadronization().

                                        {
     fEventBlock = lhee;
     if (fEventBlock == nullptr)
       return;
     const lhef::HEPEUP* hepeup = fEventBlock->getHEPEUP();
     lhef::CommonBlocks::fillHEPEUP(hepeup);
     return;
   }

void JetMatchingHook::setLHERunInfo ( lhef::LHERunInfo * lheri )

inlineprotected

Definition at line 55 of file JetMatchingHook.h.

References lhef::CommonBlocks::fillHEPRUP(), fRunBlock, and lhef::LHERunInfo::getHEPRUP().

Referenced by init().

                                             {
     fRunBlock = lheri;
     if (fRunBlock == nullptr)
       return;
     const lhef::HEPRUP* heprup = fRunBlock->getHEPRUP();
     lhef::CommonBlocks::fillHEPRUP(heprup);
     return;
   }