|
|
|
#include <JetMatchingHook.h>
Public Member Functions | |
| void | beforeHadronization (lhef::LHEEvent *lhee) |
| virtual bool | canVetoPartonLevel () |
| virtual bool | doVetoPartonLevel (const Pythia8::Event &event) |
| void | init (lhef::LHERunInfo *runInfo) |
| JetMatchingHook (const edm::ParameterSet &, Pythia8::Info *) | |
| void | resetMatchingStatus () |
| void | setEventNumber (int ievt) |
| virtual | ~JetMatchingHook () |
Protected Member Functions | |
| void | setLHEEvent (lhef::LHEEvent *lhee) |
| void | setLHERunInfo (lhef::LHERunInfo *lheri) |
Private Member Functions | |
| int | getAncestor (int, const Pythia8::Event &) |
| void | setHEPEVT (const Pythia8::Event &) |
Private Attributes | |
| lhef::LHEEvent * | fEventBlock |
| int | fEventNumber |
| Pythia8::Info * | fInfoPtr |
| gen::JetMatching * | fJetMatching |
| lhef::LHERunInfo * | fRunBlock |
Definition at line 15 of file JetMatchingHook.h.
| JetMatchingHook::JetMatchingHook | ( | const edm::ParameterSet & | , |
| Pythia8::Info * | |||
| ) |
Definition at line 23 of file JetMatchingHook.cc.
References Exception, fInfoPtr, fJetMatching, and edm::ParameterSet::getParameter().
: UserHooks(), fRunBlock(0), fEventBlock(0), fEventNumber(0), fInfoPtr(info), fJetMatching(0) { assert( fInfoPtr ); std::string scheme = ps.getParameter<std::string>("scheme"); if ( scheme == "Madgraph" ) { fJetMatching = new JetMatchingMadgraph(ps); } else if ( scheme == "MLM" || scheme == "Alpgen" ) { throw cms::Exception("JetMatching") << "Port of " << scheme << "scheme \"" << "\"" " for parton-shower matching is still in progress." << std::endl; } else throw cms::Exception("InvalidJetMatching") << "Unknown scheme \"" << scheme << "\"" " specified for parton-shower matching." << std::endl; }
| JetMatchingHook::~JetMatchingHook | ( | ) | [virtual] |
Definition at line 50 of file JetMatchingHook.cc.
References fJetMatching.
{
if ( fJetMatching ) delete fJetMatching;
}
| void JetMatchingHook::beforeHadronization | ( | lhef::LHEEvent * | lhee | ) |
Definition at line 70 of file JetMatchingHook.cc.
References gen::JetMatching::beforeHadronisation(), gen::JetMatching::beforeHadronisationExec(), fJetMatching, and setLHEEvent().
Referenced by Pythia8Hadronizer::hadronize().
{
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;
}
| virtual bool JetMatchingHook::canVetoPartonLevel | ( | ) | [inline, virtual] |
Definition at line 23 of file JetMatchingHook.h.
{ return true; }
| bool JetMatchingHook::doVetoPartonLevel | ( | const Pythia8::Event & | event | ) | [virtual] |
Definition at line 91 of file JetMatchingHook.cc.
References fJetMatching, hepeup_, gen::JetMatching::isMatchingDone(), gen::JetMatching::match(), HEPEUP_::nup, and setHEPEVT().
{
// extract "hardest" event - the output will go into workEvent,
// which is a data mamber of base class UserHooks
//
subEvent(event,true);
setHEPEVT( event ); // here we pass in "full" event, not workEvent !
// std::cout << " NPartons= " << hepeup_.nup << std::endl;
if ( !hepeup_.nup || fJetMatching->isMatchingDone() )
{
return true;
}
// Note from JY:
// match(...)input agrs here are reserved for future development and are irrelevat at this point
// just for future references, they're:
// const HepMC::GenEvent* partonLevel,
// const HepMC::GenEvent *finalState,
// bool showeredFinalState
//
bool jmtch = fJetMatching->match( 0, 0, true ); // true if veto-ed, false if accepted (not veto-ed)
if ( jmtch )
{
return true;
}
// Do not veto events that got this far
//
return false;
}
| int JetMatchingHook::getAncestor | ( | int | , |
| const Pythia8::Event & | |||
| ) | [private] |
Definition at line 234 of file JetMatchingHook.cc.
References abs, pos, and ntuplemaker::status.
Referenced by setHEPEVT().
{
int parentId = fullEvent[pos].mother1();
int parentPrevId = 0;
int counter = pos;
while ( parentId > 0 )
{
if ( parentId == fullEvent[counter].mother2() ) // carbon copy, keep walking up
{
parentPrevId = parentId;
counter = parentId;
parentId = fullEvent[parentPrevId].mother1();
continue;
}
// we get here if not a carbon copy
// let's check if it's a normal process, etc.
//
if ( (parentId < parentPrevId) || parentId < fullEvent[counter].mother2() ) // normal process
{
// first of all, check if hard block
if ( abs(fullEvent[counter].status()) == 22 || abs(fullEvent[counter].status()) == 23 )
{
// yes, it's the hard block
// we got what we want, and can exit now !
parentId = counter;
break;
}
else
{
parentPrevId = parentId;
parentId = fullEvent[parentPrevId].mother1();
}
}
else if ( parentId > parentPrevId || parentId > pos ) // "circular"/"forward-pointing parent" - intermediate process
{
parentId = -1;
break;
}
// additional checks... although we shouldn't be geeting here all that much...
//
if ( abs(fullEvent[parentId].status()) == 22 || abs(fullEvent[parentId].status())== 23 ) // hard block
{
break;
}
if ( abs(fullEvent[parentId].status()) < 22 ) // incoming
{
parentId = -1;
break;
}
}
return parentId;
}
| void JetMatchingHook::init | ( | lhef::LHERunInfo * | runInfo | ) |
Definition at line 55 of file JetMatchingHook.cc.
References Exception, fJetMatching, fRunBlock, gen::JetMatching::init(), and setLHERunInfo().
Referenced by Pythia8Hadronizer::initializeForExternalPartons().
{
setLHERunInfo( runInfo );
if ( !fRunBlock )
{
throw cms::Exception("JetMatching")
<< "Invalid RunInfo" << std::endl;
}
fJetMatching->init( runInfo );
return;
}
| void JetMatchingHook::resetMatchingStatus | ( | ) | [inline] |
Definition at line 29 of file JetMatchingHook.h.
References fJetMatching, and gen::JetMatching::resetMatchingStatus().
Referenced by Pythia8Hadronizer::hadronize().
{ fJetMatching->resetMatchingStatus(); return; }
| void JetMatchingHook::setEventNumber | ( | int | ievt | ) | [inline] |
Definition at line 26 of file JetMatchingHook.h.
References fEventNumber.
{ fEventNumber = ievt; return ; }
| void JetMatchingHook::setHEPEVT | ( | const Pythia8::Event & | ) | [private] |
Definition at line 127 of file JetMatchingHook.cc.
References abs, event(), fEventBlock, fEventNumber, getAncestor(), lhef::LHEEvent::getHEPEUP(), lhef::HEPEUP::IDUP, getHLTprescales::index, lhef::HEPEUP::MOTHUP, lhef::HEPEUP::NUP, and lhef::HEPEUP::PUP.
Referenced by doVetoPartonLevel().
{
HepMC::HEPEVT_Wrapper::zero_everything();
// service container for further mother-daughters links
//
std::vector<int> Py8PartonIdx; // position of original (LHE) partons in Py8::Event
Py8PartonIdx.clear();
std::vector<int> HEPEVTPartonIdx; // position of LHE partons in HEPEVT (incl. ME-generated decays)
HEPEVTPartonIdx.clear();
// general counter
//
int index = 0;
int Py8PartonCounter = 0;
int HEPEVTPartonCounter = 0;
// find the fisrt parton that comes from LHE (ME-generated)
// skip the incoming particles/partons
for ( int iprt=1; iprt<event.size(); iprt++ )
{
const Particle& part = event[iprt];
if ( abs(part.status()) < 22 ) continue; // below 10 is "service"
// 11-19 are beam particles; below 10 is "service"
// 21 is incoming partons
Py8PartonCounter = iprt;
break;
}
const lhef::HEPEUP& hepeup = *fEventBlock->getHEPEUP();
// start the counter from 2, because we don't want the incoming particles/oartons !
for ( int iprt=2; iprt<hepeup.NUP; iprt++ )
{
index++;
HepMC::HEPEVT_Wrapper::set_id( index, hepeup.IDUP[iprt] );
HepMC::HEPEVT_Wrapper::set_status( index, 2 );
HepMC::HEPEVT_Wrapper::set_momentum( index, hepeup.PUP[iprt][0], hepeup.PUP[iprt][1], hepeup.PUP[iprt][2], hepeup.PUP[iprt][4] );
HepMC::HEPEVT_Wrapper::set_mass( index, hepeup.PUP[iprt][4] );
// --> FIXME HepMC::HEPEVT_Wrapper::set_position( index, part.xProd(), part.yProd(), part.zProd(), part.tProd() );
HepMC::HEPEVT_Wrapper::set_parents( index, 0, 0 ); // NO, not anymore to the "system particle"
HepMC::HEPEVT_Wrapper::set_children( index, 0, 0 );
if ( hepeup.MOTHUP[iprt].first > 2 && hepeup.MOTHUP[iprt].second > 2 ) // decay from LHE, will NOT show at the start of Py8 event !!!
{
HEPEVTPartonCounter++;
continue;
}
Py8PartonIdx.push_back( Py8PartonCounter );
Py8PartonCounter++;
HEPEVTPartonIdx.push_back( HEPEVTPartonCounter);
HEPEVTPartonCounter++;
}
HepMC::HEPEVT_Wrapper::set_number_entries( index );
// now that the initial partons are in, attach parton-level from Pythia8
// do NOT reset index as we need to *add* more particles sequentially
//
for ( int iprt=1; iprt<workEvent.size(); iprt++ ) // from 0-entry (system) or from 1st entry ???
{
const Particle& part = workEvent[iprt];
if ( part.status() != 62 ) continue;
index++;
HepMC::HEPEVT_Wrapper::set_id( index, part.id() );
// HepMC::HEPEVT_Wrapper::set_status( index, event.statusHepMC(iprt) );
HepMC::HEPEVT_Wrapper::set_status( index, 1 );
HepMC::HEPEVT_Wrapper::set_momentum( index, part.px(), part.py(), part.pz(), part.e() );
HepMC::HEPEVT_Wrapper::set_mass( index, part.m() );
HepMC::HEPEVT_Wrapper::set_position( index, part.xProd(), part.yProd(), part.zProd(), part.tProd() );
HepMC::HEPEVT_Wrapper::set_parents( index, 0, 0 ); // just set to 0 like in Py6...
// although for some, mother will need to be re-set properly !
HepMC::HEPEVT_Wrapper::set_children( index, 0, 0 );
// now refine mother-daughters links, where applicable
int parentId = getAncestor( part.daughter1(), event );
if ( parentId <= 0 ) continue;
for ( int idx=0; idx<(int)Py8PartonIdx.size(); idx++ )
{
if ( parentId == Py8PartonIdx[idx] )
{
int idx1 = HEPEVTPartonIdx[idx];
HepMC::HEPEVT_Wrapper::set_parents( index, idx1+1, idx1+1 );
break;
}
}
}
HepMC::HEPEVT_Wrapper::set_number_entries( index );
HepMC::HEPEVT_Wrapper::set_event_number( fEventNumber ); // well, if you know it... well, it's one of the counters...
// HepMC::HEPEVT_Wrapper::print_hepevt();
return;
}
| void JetMatchingHook::setLHEEvent | ( | lhef::LHEEvent * | lhee | ) | [inline, protected] |
Definition at line 41 of file JetMatchingHook.h.
References fEventBlock, lhef::CommonBlocks::fillHEPEUP(), and lhef::LHEEvent::getHEPEUP().
Referenced by beforeHadronization().
{
fEventBlock=lhee;
if ( fEventBlock == 0 ) return;
const lhef::HEPEUP* hepeup = fEventBlock->getHEPEUP();
lhef::CommonBlocks::fillHEPEUP(hepeup);
return;
}
| void JetMatchingHook::setLHERunInfo | ( | lhef::LHERunInfo * | lheri | ) | [inline, protected] |
Definition at line 34 of file JetMatchingHook.h.
References lhef::CommonBlocks::fillHEPRUP(), fRunBlock, and lhef::LHERunInfo::getHEPRUP().
Referenced by init().
{
fRunBlock=lheri;
if ( fRunBlock == 0 ) return;
const lhef::HEPRUP* heprup = fRunBlock->getHEPRUP();
lhef::CommonBlocks::fillHEPRUP(heprup);
return;
}
lhef::LHEEvent* JetMatchingHook::fEventBlock [private] |
Definition at line 52 of file JetMatchingHook.h.
Referenced by setHEPEVT(), and setLHEEvent().
int JetMatchingHook::fEventNumber [private] |
Definition at line 53 of file JetMatchingHook.h.
Referenced by setEventNumber(), and setHEPEVT().
Pythia8::Info* JetMatchingHook::fInfoPtr [private] |
Definition at line 54 of file JetMatchingHook.h.
Referenced by JetMatchingHook().
gen::JetMatching* JetMatchingHook::fJetMatching [private] |
Definition at line 55 of file JetMatchingHook.h.
Referenced by beforeHadronization(), doVetoPartonLevel(), init(), JetMatchingHook(), resetMatchingStatus(), and ~JetMatchingHook().
lhef::LHERunInfo* JetMatchingHook::fRunBlock [private] |
Definition at line 51 of file JetMatchingHook.h.
Referenced by init(), and setLHERunInfo().
1.7.3