#include <PomwigHadronizer.h>

Inheritance diagram for gen::PomwigHadronizer:

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 ()

	PomwigHadronizer (const edm::ParameterSet &params)

bool	readSettings (int)

bool	residualDecay ()

void	statistics ()

	~PomwigHadronizer () override

Public Member Functions inherited from gen::BaseHadronizer
	BaseHadronizer (edm::ParameterSet const &ps)

void	cleanLHE ()

void	generateLHE (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine, unsigned int ncpu)

edm::Event &	getEDMEvent () const

std::unique_ptr< HepMC::GenEvent >	getGenEvent ()

std::unique_ptr< GenEventInfoProduct >	getGenEventInfo ()

virtual std::unique_ptr< GenLumiInfoHeader >	getGenLumiInfoHeader () const

GenRunInfoProduct &	getGenRunInfo ()

std::unique_ptr< lhef::LHEEvent >	getLHEEvent ()

const std::shared_ptr< lhef::LHERunInfo > &	getLHERunInfo () const

const std::string &	gridpackPath () const

int	randomIndex () const

const std::string &	randomInitConfigDescription () const

void	randomizeIndex (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine)

void	resetEvent (std::unique_ptr< HepMC::GenEvent > event)

void	resetEventInfo (std::unique_ptr< GenEventInfoProduct > eventInfo)

virtual bool	select (HepMC::GenEvent *) const

void	setEDMEvent (edm::Event &event)

void	setLHEEvent (std::unique_ptr< lhef::LHEEvent > event)

void	setLHERunInfo (std::unique_ptr< lhef::LHERunInfo > runInfo)

void	setRandomEngine (CLHEP::HepRandomEngine *v)

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

virtual	~BaseHadronizer () noexcept(false)

Public Member Functions inherited from gen::Herwig6Instance
bool	callWithTimeout (unsigned int secs, void(*fn)())

bool	give (const std::string &line)

	Herwig6Instance ()

void	openParticleSpecFile (const std::string fileName)

void	setHerwigRandomEngine (CLHEP::HepRandomEngine *v)

	~Herwig6Instance () override

Public Member Functions inherited from gen::FortranInstance
template<typename T >
T	call (T(&fn)())

template<typename T , typename A >
T	call (T(&fn)(A), A a)

template<typename T , typename A1 , typename A2 >
T	call (T(&fn)(A1, A2), A1 a1, A2 a2)

void	call (void(&fn)())

template<typename A >
void	call (void(&fn)(A), A a)

template<typename A1 , typename A2 >
void	call (void(&fn)(A1, A2), A1 a1, A2 a2)

virtual void	enter ()

	FortranInstance ()

virtual void	leave ()

virtual void	upEvnt ()

virtual void	upInit ()

virtual bool	upVeto ()

virtual	~FortranInstance () noexcept(false)

Private Member Functions
void	clear ()

void	doSetRandomEngine (CLHEP::HepRandomEngine *v) override

std::vector< std::string > const &	doSharedResources () const override

bool	initializeDPDF ()

Private Attributes
double	comEnergy

HepMC::IO_HERWIG	conv

int	diffTopology

bool	doMPInteraction

bool	doPDGConvert

int	h1fit

int	hepmcVerbosity

int	herwigVerbosity

int	maxEventsToPrint

bool	needClear

int	numTrials

gen::ParameterCollector	parameters

bool	printCards

double	survivalProbability

bool	useJimmy

Static Private Attributes
static const std::vector< std::string >	theSharedResources

Additional Inherited Members
Static Public Member Functions inherited from gen::FortranInstance
template<typename T >
static T *	getInstance ()

Static Public Attributes inherited from gen::FortranInstance
static const std::string	kFortranInstance = "FortranInstance"

Protected Member Functions inherited from gen::BaseHadronizer
std::unique_ptr< HepMC::GenEvent > &	event ()

std::unique_ptr< GenEventInfoProduct > &	eventInfo ()

lhef::LHEEvent *	lheEvent ()

lhef::LHERunInfo *	lheRunInfo ()

GenRunInfoProduct &	runInfo ()

Protected Member Functions inherited from gen::Herwig6Instance
virtual bool	hwwarn (const std::string &fn, int code)

Protected Attributes inherited from gen::BaseHadronizer
std::string	lheFile_

int	randomIndex_

Detailed Description

Definition at line 20 of file PomwigHadronizer.h.

Constructor & Destructor Documentation

◆ PomwigHadronizer()

gen::PomwigHadronizer::PomwigHadronizer ( const edm::ParameterSet & params )

Definition at line 89 of file PomwigHadronizer.cc.

       : BaseHadronizer(params),
         needClear(false),
         parameters(params.getParameter<edm::ParameterSet>("HerwigParameters")),
         herwigVerbosity(params.getUntrackedParameter<int>("herwigVerbosity", 0)),
         hepmcVerbosity(params.getUntrackedParameter<int>("hepmcVerbosity", 0)),
         maxEventsToPrint(params.getUntrackedParameter<int>("maxEventsToPrint", 0)),
         printCards(params.getUntrackedParameter<bool>("printCards", false)),
         comEnergy(params.getParameter<double>("comEnergy")),
         survivalProbability(params.getParameter<double>("survivalProbability")),
         diffTopology(params.getParameter<int>("diffTopology")),
         h1fit(params.getParameter<int>("h1fit")),
         useJimmy(params.getParameter<bool>("useJimmy")),
         doMPInteraction(params.getParameter<bool>("doMPInteraction")),
         numTrials(params.getUntrackedParameter<int>("numTrialsMPI", 100)),
         doPDGConvert(false) {
     if (params.exists("doPDGConvert"))
       doPDGConvert = params.getParameter<bool>("doPDGConvert");
   }

References doPDGConvert, and CalibrationSummaryClient_cfi::params.

◆ ~PomwigHadronizer()

gen::PomwigHadronizer::~PomwigHadronizer ( )

override

Definition at line 109 of file PomwigHadronizer.cc.

109 { clear(); }

References clear().

Member Function Documentation

◆ classname()

const char* gen::PomwigHadronizer::classname ( ) const

inline

Definition at line 39 of file PomwigHadronizer.h.

39 { return "PomwigHadronizer"; }

◆ clear()

void gen::PomwigHadronizer::clear ( void )

private

Definition at line 113 of file PomwigHadronizer.cc.

                                {
     if (!needClear)
       return;
  
     // teminate elementary process
     call(hwefin);
     if (useJimmy)
       call(jmefin);
  
     needClear = false;
   }

References gen::FortranInstance::call(), jmefin, needClear, and useJimmy.

Referenced by readSettings(), and ~PomwigHadronizer().

◆ decay()

bool gen::PomwigHadronizer::decay ( )

Definition at line 333 of file PomwigHadronizer.cc.

                                {
     // hadron decays
  
     InstanceWrapper wrapper(this);  // safe guard
  
     hwdhad();  // unstable particle decays
     hwdhvy();  // heavy flavour decays
     hwmevt();  // soft underlying event
  
     hwufne();  // finalize event
  
     if (hwevnt.IERROR)
       return false;
  
     event().reset(new HepMC::GenEvent);
     if (!conv.fill_next_event(event().get()))
       throw cms::Exception("PomwigError") << "HepMC Conversion problems in event." << std::endl;
  
     // do particle ID conversion Herwig->PDG, if requested
     if (doPDGConvert) {
       for (HepMC::GenEvent::particle_iterator part = event()->particles_begin(); part != event()->particles_end();
            ++part) {
         if ((*part)->pdg_id() != HepPID::translateHerwigtoPDT((*part)->pdg_id()))
           (*part)->set_pdg_id(HepPID::translateHerwigtoPDT((*part)->pdg_id()));
       }
     }
  
     return true;
   }

References conv, doPDGConvert, gen::BaseHadronizer::event(), get, and wrapper.

◆ declareSpecialSettings()

bool gen::PomwigHadronizer::declareSpecialSettings ( const std::vector< std::string > & )

inline

Definition at line 30 of file PomwigHadronizer.h.

30 { return true; }

◆ declareStableParticles()

bool gen::PomwigHadronizer::declareStableParticles ( const std::vector< int > & pdgIds )

Definition at line 271 of file PomwigHadronizer.cc.

                                                                             {
     for (std::vector<int>::const_iterator iter = pdgIds.begin(); iter != pdgIds.end(); ++iter)
       if (!markStable(*iter))
         return false;
     return true;
   }

References CosmicGenFilterHelix_cfi::pdgIds.

◆ doSetRandomEngine()

void gen::PomwigHadronizer::doSetRandomEngine ( CLHEP::HepRandomEngine * v )

overrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 111 of file PomwigHadronizer.cc.

111 { setHerwigRandomEngine(v); }

References gen::Herwig6Instance::setHerwigRandomEngine(), and gen::v.

◆ doSharedResources()

std::vector<std::string> const& gen::PomwigHadronizer::doSharedResources ( ) const

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 43 of file PomwigHadronizer.h.

43 { return theSharedResources; }

References theSharedResources.

◆ finalizeEvent()

void gen::PomwigHadronizer::finalizeEvent ( )

Definition at line 325 of file PomwigHadronizer.cc.

                                        {
     lhef::LHEEvent::fixHepMCEventTimeOrdering(event().get());
  
     event()->set_signal_process_id(hwproc.IPROC);
  
     event()->weights().push_back(hwevnt.EVWGT);
   }

References gen::BaseHadronizer::event(), lhef::LHEEvent::fixHepMCEventTimeOrdering(), and get.

◆ generatePartonsAndHadronize()

bool gen::PomwigHadronizer::generatePartonsAndHadronize ( )

Definition at line 290 of file PomwigHadronizer.cc.

                                                      {
     // hard process generation, parton shower, hadron formation
  
     InstanceWrapper wrapper(this);  // safe guard
  
     event().reset();
  
     // call herwig routines to create HEPEVT
  
     hwuine();  // initialize event
  
     if (callWithTimeout(10, hwepro)) {  // process event and PS
       // We hung for more than 10 seconds
       int error = 199;
       hwwarn_("HWHGUP", &error);
     }
  
     hwbgen();  // parton cascades
  
     // call jimmy ... only if event is not killed yet by HERWIG
     if (useJimmy && doMPInteraction && !hwevnt.IERROR && call_hwmsct())
       return false;
  
     hwdhob();  // heavy quark decays
     hwcfor();  // cluster formation
     hwcdec();  // cluster decays
  
     // if event *not* killed by HERWIG, return true
     if (!hwevnt.IERROR)
       return true;
  
     hwufne();  // finalize event
     return false;
   }

References gen::Herwig6Instance::callWithTimeout(), doMPInteraction, relativeConstraints::error, gen::BaseHadronizer::event(), gen::hwwarn_(), useJimmy, and wrapper.

◆ hadronize()

bool gen::PomwigHadronizer::hadronize ( )

Definition at line 288 of file PomwigHadronizer.cc.

288 { return false; }

◆ initializeDPDF()

bool gen::PomwigHadronizer::initializeDPDF ( )

private

Definition at line 223 of file PomwigHadronizer.cc.

                                         {
     // Initialize H1 pomeron/reggeon
  
     if (diffTopology == 3)
       return true;
  
     if ((diffTopology != 0) && (diffTopology != 1) && (diffTopology != 2))
       return false;
  
     int nstru = hwpram.NSTRU;
     int ifit = h1fit;
     if ((nstru == 9) || (nstru == 10)) {
       if ((ifit <= 0) || (ifit >= 7)) {
         throw edm::Exception(edm::errors::Configuration, "PomwigError")
             << " Attempted to set non existant H1 1997 fit index. Has to be 1...6";
       }
       std::string aux((nstru == 9) ? "Pomeron" : "Reggeon");
       edm::LogVerbatim("") << "   H1 1997 pdf's: " << aux << "\n"
                            << "   IFIT = " << ifit;
       double xp = 0.1;
       double Q2 = 75.0;
       double xpq[13];
       qcd_1994(xp, Q2, xpq, ifit);
     } else if ((nstru >= 12) && (nstru <= 15)) {
       bool isPom = (nstru == 12) || (nstru == 14);
       bool isFitA = (nstru == 12) || (nstru == 13);
       ifit = (isFitA) ? 1 : 2;
       std::string aux_0((isFitA) ? "A" : "B");
       std::string aux_1((isPom) ? "Pomeron" : "Reggeon");
       edm::LogVerbatim("") << "   H1 2006 Fit " << aux_0 << " " << aux_1 << "\n"
                            << "   IFIT = " << ifit;
       double xp = 0.1;
       double Q2 = 75.0;
       double xpq[13];
       double f2[2];
       double fl[2];
       double c2[2];
       double cl[2];
       qcd_2006(xp, Q2, ifit, xpq, f2, fl, c2, cl);
     } else {
       throw edm::Exception(edm::errors::Configuration, "PomwigError")
           << " Only running Pomeron H1 1997 (NSTRU=9), H1 2006 fit A (NSTRU=12) and H1 2006 fit B (NSTRU=14) or "
              "Reggeon H1 1997 (NSTRU=10), H1 2006 fit A (NSTRU=13) and H1 2006 fit B (NSTRU=15)";
     }
  
     return true;
   }

References printConversionInfo::aux, GetRecoTauVFromDQM_MC_cff::cl, edm::errors::Configuration, diffTopology, Exception, DeadROC_duringRun::f2, h1fit, qcd_1994, qcd_2006, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by initializeForInternalPartons().

◆ initializeForExternalPartons()

bool gen::PomwigHadronizer::initializeForExternalPartons ( )

Definition at line 125 of file PomwigHadronizer.cc.

125 { return false; }

◆ initializeForInternalPartons()

bool gen::PomwigHadronizer::initializeForInternalPartons ( )

Definition at line 210 of file PomwigHadronizer.cc.

                                                       {
     call(hwuinc);
  
     hwusta("PI0     ", 1);
  
     if (!initializeDPDF())
       return false;
  
     call(hweini);
  
     return true;
   }

References gen::FortranInstance::call(), and initializeDPDF().

◆ readSettings()

bool gen::PomwigHadronizer::readSettings ( int )

Definition at line 127 of file PomwigHadronizer.cc.

                                          {
     clear();
  
     edm::LogVerbatim("") << "----------------------------------------------\n"
                          << "Initializing PomwigHadronizer\n"
                          << "----------------------------------------------\n";
  
     // Call hwudat to set up HERWIG block data
     hwudat();
  
     // Setting basic parameters ...
     hwproc.PBEAM1 = comEnergy / 2.;
     hwproc.PBEAM2 = comEnergy / 2.;
     // Choose beam particles for POMWIG depending on topology
     switch (diffTopology) {
       case 0:  //DPE
         hwbmch.PART1[0] = 'E';
         hwbmch.PART1[1] = '-';
         hwbmch.PART2[0] = 'E';
         hwbmch.PART2[1] = '-';
         break;
       case 1:  //SD survive PART1
         hwbmch.PART1[0] = 'E';
         hwbmch.PART1[1] = '-';
         hwbmch.PART2[0] = 'P';
         hwbmch.PART2[1] = ' ';
         break;
       case 2:  //SD survive PART2
         hwbmch.PART1[0] = 'P';
         hwbmch.PART1[1] = ' ';
         hwbmch.PART2[0] = 'E';
         hwbmch.PART2[1] = '-';
         break;
       case 3:  //Non diffractive
         hwbmch.PART1[0] = 'P';
         hwbmch.PART1[1] = ' ';
         hwbmch.PART2[0] = 'P';
         hwbmch.PART2[1] = ' ';
         break;
       default:
         throw edm::Exception(edm::errors::Configuration, "PomwigError")
             << " Invalid Diff. Topology. Must be DPE(diffTopology = 0), SD particle 1 (diffTopology = 1), SD particle "
                "2 (diffTopology = 2) and Non diffractive (diffTopology = 3)";
         break;
     }
     for (int i = 2; i < 8; ++i) {
       hwbmch.PART1[i] = ' ';
       hwbmch.PART2[i] = ' ';
     }
  
     // initialize other common blocks ...
     call(hwigin);
  
     hwevnt.MAXER = 100000000;  // O(inf)
     hwpram.LWSUD = 0;          // don't write Sudakov form factors
     hwdspn.LWDEC = 0;          // don't write three/four body decays
                                // (no fort.77 and fort.88 ...)a
  
     std::memset(hwprch.AUTPDF, ' ', sizeof hwprch.AUTPDF);
     for (unsigned int i = 0; i < 2; i++) {
       hwpram.MODPDF[i] = -111;
       std::memcpy(hwprch.AUTPDF[i], "HWLHAPDF", 8);
     }
  
     hwevnt.MAXPR = maxEventsToPrint;
     hwpram.IPRINT = herwigVerbosity;
  
     edm::LogVerbatim("") << "------------------------------------\n"
                          << "Reading HERWIG parameters\n"
                          << "------------------------------------\n";
  
     for (gen::ParameterCollector::const_iterator line = parameters.begin(); line != parameters.end(); ++line) {
       edm::LogVerbatim("") << "   " << *line;
       if (!give(*line))
         throw edm::Exception(edm::errors::Configuration)
             << "Herwig 6 did not accept the following: \"" << *line << "\"." << std::endl;
     }
  
     needClear = true;
  
     return true;
   }

References gen::FortranInstance::call(), clear(), comEnergy, edm::errors::Configuration, diffTopology, Exception, gen::Herwig6Instance::give(), herwigVerbosity, hwdspn, hwprch, mps_fire::i, mps_splice::line, maxEventsToPrint, and needClear.

◆ residualDecay()

bool gen::PomwigHadronizer::residualDecay ( )

Definition at line 363 of file PomwigHadronizer.cc.

363 { return true; }

◆ statistics()

void gen::PomwigHadronizer::statistics ( )

Definition at line 278 of file PomwigHadronizer.cc.

                                     {
     double RNWGT = 1. / hwevnt.NWGTS;
     double AVWGT = hwevnt.WGTSUM * RNWGT;
  
     double xsec = 1.0e3 * AVWGT;
     xsec = survivalProbability * xsec;
  
     runInfo().setInternalXSec(xsec);
   }

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

Member Data Documentation

◆ comEnergy

double gen::PomwigHadronizer::comEnergy

private

Definition at line 58 of file PomwigHadronizer.h.

Referenced by readSettings().

◆ conv

HepMC::IO_HERWIG gen::PomwigHadronizer::conv

private

Definition at line 69 of file PomwigHadronizer.h.

Referenced by decay().

◆ diffTopology

int gen::PomwigHadronizer::diffTopology

private

Definition at line 60 of file PomwigHadronizer.h.

Referenced by initializeDPDF(), and readSettings().

◆ doMPInteraction

bool gen::PomwigHadronizer::doMPInteraction

private

Definition at line 64 of file PomwigHadronizer.h.

Referenced by generatePartonsAndHadronize().

◆ doPDGConvert

bool gen::PomwigHadronizer::doPDGConvert

private

Definition at line 67 of file PomwigHadronizer.h.

Referenced by decay(), and PomwigHadronizer().

◆ h1fit

int gen::PomwigHadronizer::h1fit

private

Definition at line 61 of file PomwigHadronizer.h.

Referenced by initializeDPDF().

◆ hepmcVerbosity

int gen::PomwigHadronizer::hepmcVerbosity

private

Definition at line 54 of file PomwigHadronizer.h.

◆ herwigVerbosity

int gen::PomwigHadronizer::herwigVerbosity

private

Definition at line 53 of file PomwigHadronizer.h.

Referenced by readSettings().

◆ maxEventsToPrint

int gen::PomwigHadronizer::maxEventsToPrint

private

Definition at line 55 of file PomwigHadronizer.h.

Referenced by readSettings().

◆ needClear

bool gen::PomwigHadronizer::needClear

private

Definition at line 50 of file PomwigHadronizer.h.

Referenced by clear(), and readSettings().

◆ numTrials

int gen::PomwigHadronizer::numTrials

private

Definition at line 65 of file PomwigHadronizer.h.

◆ parameters

gen::ParameterCollector gen::PomwigHadronizer::parameters

private

Definition at line 52 of file PomwigHadronizer.h.

◆ printCards

bool gen::PomwigHadronizer::printCards

private

Definition at line 56 of file PomwigHadronizer.h.

◆ survivalProbability

double gen::PomwigHadronizer::survivalProbability

private

Definition at line 59 of file PomwigHadronizer.h.

Referenced by statistics().

◆ theSharedResources

const std::vector< std::string > gen::PomwigHadronizer::theSharedResources

staticprivate

Initial value:

= {edm::SharedResourceNames::kHerwig6,

gen::FortranInstance::kFortranInstance}

Definition at line 48 of file PomwigHadronizer.h.

Referenced by doSharedResources().

◆ useJimmy

bool gen::PomwigHadronizer::useJimmy

private

Definition at line 63 of file PomwigHadronizer.h.

Referenced by clear(), and generatePartonsAndHadronize().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PomwigHadronizer()

◆ ~PomwigHadronizer()

Member Function Documentation

◆ classname()

◆ clear()

◆ decay()

◆ declareSpecialSettings()

◆ declareStableParticles()

◆ doSetRandomEngine()

◆ doSharedResources()

◆ finalizeEvent()

◆ generatePartonsAndHadronize()

◆ hadronize()

◆ initializeDPDF()

◆ initializeForExternalPartons()

◆ initializeForInternalPartons()

◆ readSettings()

◆ residualDecay()

◆ statistics()

Member Data Documentation

◆ comEnergy

◆ conv

◆ diffTopology

◆ doMPInteraction

◆ doPDGConvert

◆ h1fit

◆ hepmcVerbosity

◆ herwigVerbosity

◆ maxEventsToPrint

◆ needClear

◆ numTrials

◆ parameters

◆ printCards

◆ survivalProbability

◆ theSharedResources

◆ useJimmy