#include <LHEEvent.h>

Public Types
typedef LHEEventProduct::PDF	PDF

typedef LHEEventProduct::WGT	WGT

Public Member Functions
void	addComment (const std::string &line)

void	addWeight (const WGT &wgt)

std::unique_ptr< HepMC::GenEvent >	asHepMCEvent () const

void	attempted ()

void	count (LHERunInfo::CountMode count, double weight=1.0, double matchWeight=1.0)

void	fillEventInfo (HepMC::GenEvent *hepmc) const

void	fillPdfInfo (HepMC::PdfInfo *info) const

const std::vector< std::string > &	getComments () const

const HEPEUP *	getHEPEUP () const

const HEPRUP *	getHEPRUP () const

const PDF *	getPDF () const

const int	getReadAttempts ()

const std::shared_ptr< LHERunInfo > &	getRunInfo () const

	LHEEvent (const std::shared_ptr< LHERunInfo > &runInfo, const HEPEUP &hepeup)

	LHEEvent (const std::shared_ptr< LHERunInfo > &runInfo, const HEPEUP &hepeup, const LHEEventProduct::PDF *pdf, const std::vector< std::string > &comments)

	LHEEvent (const std::shared_ptr< LHERunInfo > &runInfo, const LHEEventProduct &product)

	LHEEvent (const std::shared_ptr< LHERunInfo > &runInfo, std::istream &in)

int	npLO () const

int	npNLO () const

double	originalXWGTUP () const

void	removeResonances (const std::vector< int > &ids)

const std::vector< float > &	scales () const

void	setNpLO (int n)

void	setNpNLO (int n)

void	setPDF (std::unique_ptr< PDF > pdf)

void	setScales (const std::vector< float > &scales)

const std::vector< WGT > &	weights () const

	~LHEEvent ()

Static Public Member Functions
static const HepMC::GenVertex *	findSignalVertex (const HepMC::GenEvent *event, bool status3=true)

static void	fixHepMCEventTimeOrdering (HepMC::GenEvent *event)

static void	removeParticle (lhef::HEPEUP &hepeup, int index)

Private Member Functions
HepMC::GenParticle *	makeHepMCParticle (unsigned int i) const

Static Private Member Functions
static bool	checkHepMCTree (const HepMC::GenEvent *event)

Private Attributes
std::vector< std::string >	comments

bool	counted

HEPEUP	hepeup

int	npLO_

int	npNLO_

double	originalXWGTUP_

std::unique_ptr< PDF >	pdf

int	readAttemptCounter

const std::shared_ptr< LHERunInfo >	runInfo

std::vector< float >	scales_

std::vector< WGT >	weights_

Detailed Description

Definition at line 23 of file LHEEvent.h.

Member Typedef Documentation

◆ PDF

typedef LHEEventProduct::PDF lhef::LHEEvent::PDF

Definition at line 34 of file LHEEvent.h.

◆ WGT

typedef LHEEventProduct::WGT lhef::LHEEvent::WGT

Definition at line 35 of file LHEEvent.h.

Constructor & Destructor Documentation

◆ LHEEvent() [1/4]

lhef::LHEEvent::LHEEvent	(	const std::shared_ptr< LHERunInfo > &	runInfo,
		std::istream &	in
	)

Definition at line 35 of file LHEEvent.cc.

       : runInfo(runInfo),
         weights_(0),
         counted(false),
         readAttemptCounter(0),
         npLO_(-99),
         npNLO_(-99)
  
   {
     hepeup.NUP = 0;
     hepeup.XPDWUP.first = hepeup.XPDWUP.second = 0.0;
  
     in >> hepeup.NUP >> hepeup.IDPRUP >> hepeup.XWGTUP >> hepeup.SCALUP >> hepeup.AQEDUP >> hepeup.AQCDUP;
     if (!in.good())
       throw cms::Exception("InvalidFormat") << "Les Houches file contained invalid"
                                                " event header."
                                             << std::endl;
  
     // store the original value of XWGTUP for the user
     originalXWGTUP_ = hepeup.XWGTUP;
  
     int idwtup = runInfo->getHEPRUP()->IDWTUP;
     if (idwtup >= 0 && hepeup.XWGTUP < 0) {
       edm::LogWarning("Generator|LHEInterface") << "Non-allowed negative event weight encountered." << std::endl;
       hepeup.XWGTUP = std::abs(hepeup.XWGTUP);
     }
  
     if (std::abs(idwtup) == 3 && std::abs(hepeup.XWGTUP) != 1.) {
       edm::LogInfo("Generator|LHEInterface") << "Event weight not set to one for abs(IDWTUP) == 3" << std::endl;
       hepeup.XWGTUP = hepeup.XWGTUP > 0. ? 1.0 : -1.0;
     }
  
     hepeup.resize();
  
     for (int i = 0; i < hepeup.NUP; i++) {
       in >> hepeup.IDUP[i] >> hepeup.ISTUP[i] >> hepeup.MOTHUP[i].first >> hepeup.MOTHUP[i].second >>
           hepeup.ICOLUP[i].first >> hepeup.ICOLUP[i].second >> hepeup.PUP[i][0] >> hepeup.PUP[i][1] >>
           hepeup.PUP[i][2] >> hepeup.PUP[i][3] >> hepeup.PUP[i][4] >> hepeup.VTIMUP[i] >> hepeup.SPINUP[i];
       if (!in.good())
         throw cms::Exception("InvalidFormat") << "Les Houches file contained invalid event"
                                                  " in particle line "
                                               << (i + 1) << "." << std::endl;
     }
  
     while (skipWhitespace(in) == '#') {
       std::string line;
       std::getline(in, line);
       std::istringstream ss(line);
       std::string tag;
       ss >> tag;
       if (tag == "#pdf") {
         pdf.reset(new PDF);
         ss >> pdf->id.first >> pdf->id.second >> pdf->x.first >> pdf->x.second >> pdf->scalePDF >> pdf->xPDF.first >>
             pdf->xPDF.second;
         if (ss.bad()) {
           edm::LogWarning("Generator|LHEInterface") << "Les Houches event contained"
                                                        " unparseable PDF information."
                                                     << std::endl;
           pdf.reset();
         } else
           continue;
       }
       comments.push_back(line + "\n");
     }
  
     if (!in.eof())
       edm::LogWarning("Generator|LHEInterface") << "Les Houches file contained spurious"
                                                    " content after event data."
                                                 << std::endl;
   }

References funct::abs(), lhef::HEPEUP::AQCDUP, lhef::HEPEUP::AQEDUP, comments, hepeup, mps_fire::i, lhef::HEPEUP::ICOLUP, lhef::HEPEUP::IDPRUP, lhef::HEPEUP::IDUP, recoMuon::in, lhef::HEPEUP::ISTUP, mps_splice::line, lhef::HEPEUP::MOTHUP, lhef::HEPEUP::NUP, originalXWGTUP_, pdf, lhef::HEPEUP::PUP, lhef::HEPEUP::resize(), runInfo, lhef::HEPEUP::SCALUP, skipWhitespace(), lhef::HEPEUP::SPINUP, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, lhef::HEPEUP::VTIMUP, lhef::HEPEUP::XPDWUP, and lhef::HEPEUP::XWGTUP.

◆ LHEEvent() [2/4]

lhef::LHEEvent::LHEEvent	(	const std::shared_ptr< LHERunInfo > &	runInfo,
		const HEPEUP &	hepeup
	)

Definition at line 106 of file LHEEvent.cc.

107 : runInfo(runInfo), hepeup(hepeup), counted(false), readAttemptCounter(0), npLO_(-99), npNLO_(-99) {}

◆ LHEEvent() [3/4]

lhef::LHEEvent::LHEEvent	(	const std::shared_ptr< LHERunInfo > &	runInfo,
		const HEPEUP &	hepeup,
		const LHEEventProduct::PDF *	pdf,
		const std::vector< std::string > &	comments
	)

Definition at line 109 of file LHEEvent.cc.

       : runInfo(runInfo),
         hepeup(hepeup),
         pdf(pdf ? new PDF(*pdf) : nullptr),
         comments(comments),
         counted(false),
         readAttemptCounter(0),
         npLO_(-99),
         npNLO_(-99) {}

◆ LHEEvent() [4/4]

lhef::LHEEvent::LHEEvent	(	const std::shared_ptr< LHERunInfo > &	runInfo,
		const LHEEventProduct &	product
	)

Definition at line 122 of file LHEEvent.cc.

       : runInfo(runInfo),
         hepeup(product.hepeup()),
         pdf(product.pdf() ? new PDF(*product.pdf()) : nullptr),
         weights_(product.weights()),
         comments(product.comments_begin(), product.comments_end()),
         counted(false),
         readAttemptCounter(0),
         originalXWGTUP_(product.originalXWGTUP()),
         scales_(product.scales()),
         npLO_(product.npLO()),
         npNLO_(product.npNLO()) {}

◆ ~LHEEvent()

lhef::LHEEvent::~LHEEvent ( )

Definition at line 135 of file LHEEvent.cc.

135 {}

Member Function Documentation

◆ addComment()

void lhef::LHEEvent::addComment ( const std::string & line )

inline

Definition at line 59 of file LHEEvent.h.

59 { comments.push_back(line); }

References comments, and mps_splice::line.

◆ addWeight()

void lhef::LHEEvent::addWeight ( const WGT & wgt )

inline

Definition at line 44 of file LHEEvent.h.

44 { weights_.push_back(wgt); }

References weights_.

◆ asHepMCEvent()

std::unique_ptr< HepMC::GenEvent > lhef::LHEEvent::asHepMCEvent ( ) const

Definition at line 225 of file LHEEvent.cc.

                                                             {
     std::unique_ptr<HepMC::GenEvent> hepmc(new HepMC::GenEvent);
  
     hepmc->set_signal_process_id(hepeup.IDPRUP);
     hepmc->set_event_scale(hepeup.SCALUP);
     hepmc->set_alphaQED(hepeup.AQEDUP);
     hepmc->set_alphaQCD(hepeup.AQCDUP);
  
     unsigned int nup = hepeup.NUP;  // particles in event
  
     // any particles in HEPEUP block?
     if (!nup) {
       edm::LogWarning("Generator|LHEInterface") << "Les Houches Event does not contain any partons. "
                                                 << "Not much to convert.";
       return hepmc;
     }
  
     // stores (pointers to) converted particles
     std::vector<HepMC::GenParticle *> genParticles;
     std::vector<HepMC::GenVertex *> genVertices;
  
     // I. convert particles
     for (unsigned int i = 0; i < nup; i++)
       genParticles.push_back(makeHepMCParticle(i));
  
     // II. loop again to build vertices
     for (unsigned int i = 0; i < nup; i++) {
       unsigned int mother1 = hepeup.MOTHUP.at(i).first;
       unsigned int mother2 = hepeup.MOTHUP.at(i).second;
       double cTau = hepeup.VTIMUP.at(i);  // decay time
  
       // current particle has a mother? --- Sorry, parent! We're PC.
       if (mother1) {
         mother1--;  // FORTRAN notation!
         if (mother2)
           mother2--;
         else
           mother2 = mother1;
  
         HepMC::GenParticle *in_par = genParticles.at(mother1);
         HepMC::GenVertex *current_vtx = in_par->end_vertex();  // vertex of first mother
  
         if (!current_vtx) {
           current_vtx = new HepMC::GenVertex(HepMC::FourVector(0, 0, 0, cTau));
  
           // add vertex to event
           genVertices.push_back(current_vtx);
         }
  
         for (unsigned int j = mother1; j <= mother2; j++)  // set mother-daughter relations
           if (!genParticles.at(j)->end_vertex())
             current_vtx->add_particle_in(genParticles.at(j));
  
         // connect THIS outgoing particle to current vertex
         current_vtx->add_particle_out(genParticles.at(i));
       }
     }
  
     checkHepMCTree(hepmc.get());
  
     // III. restore color flow
     // ok, nobody knows how to do it so far...
  
     // IV. fill run information
     const HEPRUP *heprup = getHEPRUP();
  
     // set beam particles
     HepMC::GenParticle *b1 = new HepMC::GenParticle(
         HepMC::FourVector(0.0, 0.0, +heprup->EBMUP.first, heprup->EBMUP.first), heprup->IDBMUP.first);
     HepMC::GenParticle *b2 = new HepMC::GenParticle(
         HepMC::FourVector(0.0, 0.0, -heprup->EBMUP.second, heprup->EBMUP.second), heprup->IDBMUP.second);
     b1->set_status(3);
     b2->set_status(3);
  
     HepMC::GenVertex *v1 = new HepMC::GenVertex();
     HepMC::GenVertex *v2 = new HepMC::GenVertex();
     v1->add_particle_in(b1);
     v2->add_particle_in(b2);
  
     hepmc->add_vertex(v1);
     hepmc->add_vertex(v2);
     hepmc->set_beam_particles(b1, b2);
  
     // first two particles have to be the hard partons going into the interaction
     if (genParticles.size() >= 2) {
       if (!genParticles.at(0)->production_vertex() && !genParticles.at(1)->production_vertex()) {
         v1->add_particle_out(genParticles.at(0));
         v2->add_particle_out(genParticles.at(1));
       } else
         edm::LogWarning("Generator|LHEInterface") << "Initial partons do already have a"
                                                      " production vertex. "
                                                   << std::endl
                                                   << "Beam particles not connected.";
     } else
       edm::LogWarning("Generator|LHEInterface") << "Can't find any initial partons to be"
                                                    " connected to the beam particles.";
  
     for (std::vector<HepMC::GenVertex *>::const_iterator iter = genVertices.begin(); iter != genVertices.end(); ++iter)
       hepmc->add_vertex(*iter);
  
     // do some more consistency checks
     for (unsigned int i = 0; i < nup; i++) {
       if (!genParticles.at(i)->parent_event()) {
         edm::LogWarning("Generator|LHEInterface") << "Not all LHE particles could be stored"
                                                      "  stored in the HepMC event. "
                                                   << std::endl
                                                   << "Check the mother-daughter relations"
                                                      " in the given LHE input file.";
         break;
       }
     }
  
     hepmc->set_signal_process_vertex(const_cast<HepMC::GenVertex *>(findSignalVertex(hepmc.get(), false)));
  
     return hepmc;
   }

References lhef::HEPEUP::AQCDUP, lhef::HEPEUP::AQEDUP, b1, b2, checkHepMCTree(), pdg::cTau(), findSignalVertex(), GenParticle::GenParticle, genParticles2HepMC_cfi::genParticles, getHEPRUP(), hepeup, mps_fire::i, lhef::HEPEUP::IDPRUP, dqmiolumiharvest::j, makeHepMCParticle(), lhef::HEPEUP::MOTHUP, lhef::HEPEUP::NUP, lhef::HEPEUP::SCALUP, and lhef::HEPEUP::VTIMUP.

◆ attempted()

void lhef::LHEEvent::attempted ( )

inline

Definition at line 66 of file LHEEvent.h.

                      {
       readAttemptCounter++;
       return;
     }

References readAttemptCounter.

◆ checkHepMCTree()

bool lhef::LHEEvent::checkHepMCTree ( const HepMC::GenEvent * event )

staticprivate

Definition at line 355 of file LHEEvent.cc.

                                                           {
     double px = 0, py = 0, pz = 0, E = 0;
  
     for (HepMC::GenEvent::particle_const_iterator iter = event->particles_begin(); iter != event->particles_end();
          iter++) {
       int status = (*iter)->status();
       HepMC::FourVector fv = (*iter)->momentum();
  
       // incoming particles
       if (status == 3 && *iter != event->beam_particles().first && *iter != event->beam_particles().second) {
         px -= fv.px();
         py -= fv.py();
         pz -= fv.pz();
         E -= fv.e();
       }
  
       // outgoing particles
       if (status == 1) {
         px += fv.px();
         py += fv.py();
         pz += fv.pz();
         E += fv.e();
       }
     }
  
     if (px * px + py * py + pz * pz + E * E > 0.1) {
       edm::LogWarning("Generator|LHEInterface")
           << "Energy-momentum badly conserved. " << std::setprecision(3) << "sum p_i  = [" << std::setw(7) << E << ", "
           << std::setw(7) << px << ", " << std::setw(7) << py << ", " << std::setw(7) << pz << "]";
  
       return false;
     }
  
     return true;
   }

References multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, and mps_update::status.

Referenced by asHepMCEvent().

◆ count()

void lhef::LHEEvent::count	(	LHERunInfo::CountMode	count,
		double	weight = `1.0`,
		double	matchWeight = `1.0`
	)

Definition at line 182 of file LHEEvent.cc.

                                                                                   {
     if (counted)
       edm::LogWarning("Generator|LHEInterface") << "LHEEvent::count called twice on same event!" << std::endl;
  
     runInfo->count(hepeup.IDPRUP, mode, hepeup.XWGTUP, weight, matchWeight);
  
     counted = true;
   }

References counted, hepeup, lhef::HEPEUP::IDPRUP, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, runInfo, and lhef::HEPEUP::XWGTUP.

Referenced by gen::Pythia6Hadronizer::hadronize(), Pythia8Hadronizer::hadronize(), and Herwig6Hadronizer::hadronize().

◆ fillEventInfo()

void lhef::LHEEvent::fillEventInfo ( HepMC::GenEvent * hepmc ) const

Definition at line 218 of file LHEEvent.cc.

                                                          {
     event->set_signal_process_id(hepeup.IDPRUP);
     event->set_event_scale(hepeup.SCALUP);
     event->set_alphaQED(hepeup.AQEDUP);
     event->set_alphaQCD(hepeup.AQCDUP);
   }

References lhef::HEPEUP::AQCDUP, lhef::HEPEUP::AQEDUP, hepeup, lhef::HEPEUP::IDPRUP, and lhef::HEPEUP::SCALUP.

Referenced by gen::Pythia6Hadronizer::finalizeEvent(), and Herwig6Hadronizer::finalizeEvent().

◆ fillPdfInfo()

void lhef::LHEEvent::fillPdfInfo ( HepMC::PdfInfo * info ) const

Definition at line 191 of file LHEEvent.cc.

                                                      {
     if (pdf.get()) {
       info->set_id1(pdf->id.first);
       info->set_id2(pdf->id.second);
       info->set_x1(pdf->x.first);
       info->set_x2(pdf->x.second);
       info->set_pdf1(pdf->xPDF.first);
       info->set_pdf2(pdf->xPDF.second);
       info->set_scalePDF(pdf->scalePDF);
     } else if (hepeup.NUP >= 2) {
       const HEPRUP *heprup = getHEPRUP();
       info->set_id1(hepeup.IDUP[0] == 21 ? 0 : hepeup.IDUP[0]);
       info->set_id2(hepeup.IDUP[1] == 21 ? 0 : hepeup.IDUP[1]);
       info->set_x1(std::abs(hepeup.PUP[0][2] / heprup->EBMUP.first));
       info->set_x2(std::abs(hepeup.PUP[1][2] / heprup->EBMUP.second));
       info->set_pdf1(-1.0);
       info->set_pdf2(-1.0);
       info->set_scalePDF(hepeup.SCALUP);
     } else {
       info->set_x1(-1.0);
       info->set_x2(-1.0);
       info->set_pdf1(-1.0);
       info->set_pdf2(-1.0);
       info->set_scalePDF(hepeup.SCALUP);
     }
   }

References funct::abs(), lhef::HEPRUP::EBMUP, getHEPRUP(), hepeup, lhef::HEPEUP::IDUP, info(), lhef::HEPEUP::NUP, pdf, lhef::HEPEUP::PUP, and lhef::HEPEUP::SCALUP.

Referenced by gen::Pythia6Hadronizer::finalizeEvent(), and Herwig6Hadronizer::finalizeEvent().

◆ findSignalVertex()

const HepMC::GenVertex * lhef::LHEEvent::findSignalVertex	(	const HepMC::GenEvent *	event,
		bool	status3 = `true`
	)

static

Definition at line 391 of file LHEEvent.cc.

                                                                                            {
     double largestMass2 = -9.0e-30;
     const HepMC::GenVertex *vertex = nullptr;
     for (HepMC::GenEvent::vertex_const_iterator iter = event->vertices_begin(); iter != event->vertices_end(); ++iter) {
       if ((*iter)->particles_in_size() < 2)
         continue;
       if ((*iter)->particles_out_size() < 1 ||
           ((*iter)->particles_out_size() == 1 &&
            (!(*(*iter)->particles_out_const_begin())->end_vertex() ||
             !(*(*iter)->particles_out_const_begin())->end_vertex()->particles_out_size())))
         continue;
  
       double px = 0.0, py = 0.0, pz = 0.0, E = 0.0;
       bool hadStatus3 = false;
       for (HepMC::GenVertex::particles_out_const_iterator iter2 = (*iter)->particles_out_const_begin();
            iter2 != (*iter)->particles_out_const_end();
            ++iter2) {
         hadStatus3 = hadStatus3 || (*iter2)->status() == 3;
         px += (*iter2)->momentum().px();
         py += (*iter2)->momentum().py();
         pz += (*iter2)->momentum().pz();
         E += (*iter2)->momentum().e();
       }
       if (status3 && !hadStatus3)
         continue;
  
       double mass2 = E * E - (px * px + py * py + pz * pz);
       if (mass2 > largestMass2) {
         vertex = *iter;
         largestMass2 = mass2;
       }
     }
  
     return vertex;
   }

References multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, and bphysicsOniaDQM_cfi::vertex.

Referenced by asHepMCEvent().

◆ fixHepMCEventTimeOrdering()

void lhef::LHEEvent::fixHepMCEventTimeOrdering ( HepMC::GenEvent * event )

static

Definition at line 451 of file LHEEvent.cc.

                                                                {
     std::set<const HepMC::GenVertex *> visited;
     HepMC::FourVector zeroTime;
     for (HepMC::GenEvent::vertex_iterator iter = event->vertices_begin(); iter != event->vertices_end(); ++iter)
       fixSubTree(*iter, zeroTime, visited);
   }

References lhef::fixSubTree(), and class-composition::visited.

Referenced by gen::PomwigHadronizer::finalizeEvent(), and Herwig6Hadronizer::finalizeEvent().

◆ getComments()

const std::vector<std::string>& lhef::LHEEvent::getComments ( ) const

inline

Definition at line 41 of file LHEEvent.h.

41 { return comments; }

References comments.

Referenced by Herwig6Hadronizer::upEvnt().

◆ getHEPEUP()

const HEPEUP* lhef::LHEEvent::getHEPEUP ( ) const

inline

Definition at line 38 of file LHEEvent.h.

38 { return &hepeup; }

References hepeup.

Referenced by gen::JetMatchingMGFastJet::beforeHadronisation(), gen::TauolappInterface::decay(), gen::FortranCallback::fillEvent(), Py8toJetInput::fillJetAlgoInput(), Py8toJetInputHEPEVT::fillJetAlgoInput(), gen::JetMatchingMGFastJet::match(), and JetMatchingHook::setLHEEvent().

◆ getHEPRUP()

const HEPRUP* lhef::LHEEvent::getHEPRUP ( ) const

inline

Definition at line 39 of file LHEEvent.h.

39 { return runInfo->getHEPRUP(); }

References runInfo.

Referenced by asHepMCEvent(), and fillPdfInfo().

◆ getPDF()

const PDF* lhef::LHEEvent::getPDF ( ) const

inline

Definition at line 40 of file LHEEvent.h.

40 { return pdf.get(); }

References pdf.

◆ getReadAttempts()

const int lhef::LHEEvent::getReadAttempts ( )

inline

Definition at line 42 of file LHEEvent.h.

42 { return readAttemptCounter; }

References readAttemptCounter.

◆ getRunInfo()

const std::shared_ptr<LHERunInfo>& lhef::LHEEvent::getRunInfo ( ) const

inline

Definition at line 37 of file LHEEvent.h.

37 { return runInfo; }

References runInfo.

◆ makeHepMCParticle()

HepMC::GenParticle * lhef::LHEEvent::makeHepMCParticle ( unsigned int i ) const

private

Definition at line 342 of file LHEEvent.cc.

                                                                     {
     HepMC::GenParticle *particle = new HepMC::GenParticle(
         HepMC::FourVector(hepeup.PUP.at(i)[0], hepeup.PUP.at(i)[1], hepeup.PUP.at(i)[2], hepeup.PUP.at(i)[3]),
         hepeup.IDUP.at(i));
  
     int status = hepeup.ISTUP.at(i);
  
     particle->set_generated_mass(hepeup.PUP.at(i)[4]);
     particle->set_status(status > 0 ? (status == 2 ? 3 : status) : 3);
  
     return particle;
   }

References GenParticle::GenParticle, hepeup, mps_fire::i, lhef::HEPEUP::IDUP, lhef::HEPEUP::ISTUP, lhef::HEPEUP::PUP, and mps_update::status.

Referenced by asHepMCEvent().

◆ npLO()

int lhef::LHEEvent::npLO ( ) const

inline

Definition at line 53 of file LHEEvent.h.

53 { return npLO_; }

References npLO_.

◆ npNLO()

int lhef::LHEEvent::npNLO ( ) const

inline

Definition at line 54 of file LHEEvent.h.

54 { return npNLO_; }

References npNLO_.

◆ originalXWGTUP()

double lhef::LHEEvent::originalXWGTUP ( ) const

inline

Definition at line 47 of file LHEEvent.h.

47 { return originalXWGTUP_; }

References originalXWGTUP_.

◆ removeParticle()

void lhef::LHEEvent::removeParticle	(	lhef::HEPEUP &	hepeup,
		int	index
	)

static

Definition at line 137 of file LHEEvent.cc.

                                                          {
     index--;
     if (index < 0 || index >= hepeup.NUP) {
       edm::LogError("Generator|LHEInterface")
           << "removeParticle: Index " << (index + 1) << " out of bounds." << std::endl;
       return;
     }
  
     std::pair<int, int> mo = hepeup.MOTHUP[index];
  
     hepeup.IDUP.erase(hepeup.IDUP.begin() + index);
     hepeup.ISTUP.erase(hepeup.ISTUP.begin() + index);
     hepeup.MOTHUP.erase(hepeup.MOTHUP.begin() + index);
     hepeup.ICOLUP.erase(hepeup.ICOLUP.begin() + index);
     hepeup.PUP.erase(hepeup.PUP.begin() + index);
     hepeup.VTIMUP.erase(hepeup.VTIMUP.begin() + index);
     hepeup.SPINUP.erase(hepeup.SPINUP.begin() + index);
  
     hepeup.NUP--;
  
     index++;
     for (int i = 0; i < hepeup.NUP; i++) {
       if (hepeup.MOTHUP[i].first == index) {
         if (hepeup.MOTHUP[i].second > 0)
           edm::LogError("Generator|LHEInterface")
               << "removeParticle: Particle " << (i + 2) << " has two mothers." << std::endl;
         hepeup.MOTHUP[i] = mo;
       }
  
       if (hepeup.MOTHUP[i].first > index)
         hepeup.MOTHUP[i].first--;
       if (hepeup.MOTHUP[i].second > index)
         hepeup.MOTHUP[i].second--;
     }
   }

References hepeup, mps_fire::i, lhef::HEPEUP::ICOLUP, lhef::HEPEUP::IDUP, lhef::HEPEUP::ISTUP, lhef::HEPEUP::MOTHUP, lhef::HEPEUP::NUP, lhef::HEPEUP::PUP, lhef::HEPEUP::SPINUP, and lhef::HEPEUP::VTIMUP.

Referenced by removeResonances().

◆ removeResonances()

void lhef::LHEEvent::removeResonances ( const std::vector< int > & ids )

Definition at line 173 of file LHEEvent.cc.

                                                            {
     for (int i = 1; i <= hepeup.NUP; i++) {
       int id = std::abs(hepeup.IDUP[i - 1]);
       if (hepeup.MOTHUP[i - 1].first > 0 && hepeup.MOTHUP[i - 1].second > 0 &&
           std::find(ids.begin(), ids.end(), id) != ids.end())
         removeParticle(hepeup, i--);
     }
   }

References funct::abs(), spr::find(), hepeup, mps_fire::i, triggerObjects_cff::id, lhef::HEPEUP::IDUP, lhef::HEPEUP::MOTHUP, lhef::HEPEUP::NUP, and removeParticle().

◆ scales()

const std::vector<float>& lhef::LHEEvent::scales ( ) const

inline

Definition at line 50 of file LHEEvent.h.

50 { return scales_; }

References scales_.

Referenced by setScales().

◆ setNpLO()

void lhef::LHEEvent::setNpLO ( int n )

inline

Definition at line 56 of file LHEEvent.h.

56 { npLO_ = n; }

References dqmiodumpmetadata::n, and npLO_.

◆ setNpNLO()

void lhef::LHEEvent::setNpNLO ( int n )

inline

Definition at line 57 of file LHEEvent.h.

57 { npNLO_ = n; }

References dqmiodumpmetadata::n, and npNLO_.

◆ setPDF()

void lhef::LHEEvent::setPDF ( std::unique_ptr< PDF > pdf )

inline

Definition at line 45 of file LHEEvent.h.

45 { this->pdf = std::move(pdf); }

References eostools::move(), and pdf.

◆ setScales()

void lhef::LHEEvent::setScales ( const std::vector< float > & scales )

inline

Definition at line 51 of file LHEEvent.h.

51 { scales_ = scales; }

References scales(), and scales_.

◆ weights()

const std::vector<WGT>& lhef::LHEEvent::weights ( ) const

inline

Definition at line 48 of file LHEEvent.h.

48 { return weights_; }

References weights_.

Member Data Documentation

◆ comments

std::vector<std::string> lhef::LHEEvent::comments

private

Definition at line 89 of file LHEEvent.h.

Referenced by addComment(), getComments(), and LHEEvent().

◆ counted

bool lhef::LHEEvent::counted

private

Definition at line 90 of file LHEEvent.h.

Referenced by count().

◆ hepeup

HEPEUP lhef::LHEEvent::hepeup

private

Definition at line 86 of file LHEEvent.h.

Referenced by asHepMCEvent(), count(), fillEventInfo(), fillPdfInfo(), getHEPEUP(), LHEEvent(), makeHepMCParticle(), removeParticle(), and removeResonances().

◆ npLO_

int lhef::LHEEvent::npLO_

private

Definition at line 94 of file LHEEvent.h.

Referenced by npLO(), and setNpLO().

◆ npNLO_

int lhef::LHEEvent::npNLO_

private

Definition at line 95 of file LHEEvent.h.

Referenced by npNLO(), and setNpNLO().

◆ originalXWGTUP_

double lhef::LHEEvent::originalXWGTUP_

private

Definition at line 92 of file LHEEvent.h.

Referenced by LHEEvent(), and originalXWGTUP().

◆ pdf

std::unique_ptr<PDF> lhef::LHEEvent::pdf

private

Definition at line 87 of file LHEEvent.h.

Referenced by fillPdfInfo(), getPDF(), LHEEvent(), and setPDF().

◆ readAttemptCounter

int lhef::LHEEvent::readAttemptCounter

private

Definition at line 91 of file LHEEvent.h.

Referenced by attempted(), and getReadAttempts().

◆ runInfo

const std::shared_ptr<LHERunInfo> lhef::LHEEvent::runInfo

private

Definition at line 84 of file LHEEvent.h.

Referenced by count(), getHEPRUP(), getRunInfo(), and LHEEvent().

◆ scales_

std::vector<float> lhef::LHEEvent::scales_

private

Definition at line 93 of file LHEEvent.h.

Referenced by scales(), and setScales().

◆ weights_

std::vector<WGT> lhef::LHEEvent::weights_

private

Definition at line 88 of file LHEEvent.h.

Referenced by addWeight(), and weights().

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ PDF

◆ WGT

Constructor & Destructor Documentation

◆ LHEEvent() [1/4]

◆ LHEEvent() [2/4]

◆ LHEEvent() [3/4]

◆ LHEEvent() [4/4]

◆ ~LHEEvent()

Member Function Documentation

◆ addComment()

◆ addWeight()

◆ asHepMCEvent()

◆ attempted()

◆ checkHepMCTree()

◆ count()

◆ fillEventInfo()

◆ fillPdfInfo()

◆ findSignalVertex()

◆ fixHepMCEventTimeOrdering()

◆ getComments()

◆ getHEPEUP()

◆ getHEPRUP()

◆ getPDF()

◆ getReadAttempts()

◆ getRunInfo()

◆ makeHepMCParticle()

◆ npLO()

◆ npNLO()

◆ originalXWGTUP()

◆ removeParticle()

◆ removeResonances()

◆ scales()

◆ setNpLO()

◆ setNpNLO()

◆ setPDF()

◆ setScales()

◆ weights()

Member Data Documentation

◆ comments

◆ counted

◆ hepeup

◆ npLO_

◆ npNLO_

◆ originalXWGTUP_

◆ pdf

◆ readAttemptCounter

◆ runInfo

◆ scales_

◆ weights_