#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
void	call (void(&fn)())

template<typename T >
T	call (T(&fn)())

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

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

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

template<typename T , typename A1 , typename A2 >
T	call (T(&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

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

Definition at line 91 of file PomwigHadronizer.cc.

References doPDGConvert, edm::ParameterSet::exists(), and edm::ParameterSet::getParameter().

       : 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");
   }

gen::PomwigHadronizer::~PomwigHadronizer ( )

override

Definition at line 111 of file PomwigHadronizer.cc.

References clear().

111 { clear(); }

gen::PomwigHadronizer::clear

void clear()

Definition: PomwigHadronizer.cc:115

Member Function Documentation

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

inline

Definition at line 39 of file PomwigHadronizer.h.

39 { return "PomwigHadronizer"; }

void gen::PomwigHadronizer::clear ( void )

private

Definition at line 115 of file PomwigHadronizer.cc.

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

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

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

bool gen::PomwigHadronizer::decay ( )

Definition at line 335 of file PomwigHadronizer.cc.

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

                                {
     // 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() = std::make_unique<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;
   }

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

inline

Definition at line 30 of file PomwigHadronizer.h.

30 { return true; }

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

Definition at line 273 of file PomwigHadronizer.cc.

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

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

overrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 113 of file PomwigHadronizer.cc.

References gen::Herwig6Instance::setHerwigRandomEngine().

113 { setHerwigRandomEngine(v); }

gen::v

double v[5][pyjets_maxn]

Definition: Cascade2Hadronizer.cc:76

gen::Herwig6Instance::setHerwigRandomEngine

void setHerwigRandomEngine(CLHEP::HepRandomEngine *v)

Definition: Herwig6Instance.h:38

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

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 43 of file PomwigHadronizer.h.

References theSharedResources.

43 { return theSharedResources; }

gen::PomwigHadronizer::theSharedResources

static const std::vector< std::string > theSharedResources

Definition: PomwigHadronizer.h:48

void gen::PomwigHadronizer::finalizeEvent ( )

Definition at line 327 of file PomwigHadronizer.cc.

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

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

bool gen::PomwigHadronizer::generatePartonsAndHadronize ( )

Definition at line 292 of file PomwigHadronizer.cc.

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

                                                      {
     // 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;
   }

bool gen::PomwigHadronizer::hadronize ( )

Definition at line 290 of file PomwigHadronizer.cc.

290 { return false; }

bool gen::PomwigHadronizer::initializeDPDF ( )

private

Definition at line 225 of file PomwigHadronizer.cc.

References printConversionInfo::aux, haddnano::cl, edm::errors::Configuration, diffTopology, Exception, validate-o2o-wbm::f2, h1fit, qcd_1994, qcd_2006, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by initializeForInternalPartons().

                                         {
     // 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;
   }

bool gen::PomwigHadronizer::initializeForExternalPartons ( )

Definition at line 127 of file PomwigHadronizer.cc.

127 { return false; }

bool gen::PomwigHadronizer::initializeForInternalPartons ( )

Definition at line 212 of file PomwigHadronizer.cc.

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

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

bool gen::PomwigHadronizer::readSettings ( int )

Definition at line 129 of file PomwigHadronizer.cc.

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

                                          {
     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;
   }

bool gen::PomwigHadronizer::residualDecay ( )

Definition at line 365 of file PomwigHadronizer.cc.

365 { return true; }

void gen::PomwigHadronizer::statistics ( )

Definition at line 280 of file PomwigHadronizer.cc.

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

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

Member Data Documentation

double gen::PomwigHadronizer::comEnergy

private

Definition at line 58 of file PomwigHadronizer.h.

Referenced by readSettings().

HepMC::IO_HERWIG gen::PomwigHadronizer::conv

private

Definition at line 69 of file PomwigHadronizer.h.

Referenced by decay().

int gen::PomwigHadronizer::diffTopology

private

Definition at line 60 of file PomwigHadronizer.h.

Referenced by initializeDPDF(), and readSettings().

bool gen::PomwigHadronizer::doMPInteraction

private

Definition at line 64 of file PomwigHadronizer.h.

Referenced by generatePartonsAndHadronize().

bool gen::PomwigHadronizer::doPDGConvert

private

Definition at line 67 of file PomwigHadronizer.h.

Referenced by decay(), and PomwigHadronizer().

int gen::PomwigHadronizer::h1fit

private

Definition at line 61 of file PomwigHadronizer.h.

Referenced by initializeDPDF().

int gen::PomwigHadronizer::hepmcVerbosity

private

Definition at line 54 of file PomwigHadronizer.h.

int gen::PomwigHadronizer::herwigVerbosity

private

Definition at line 53 of file PomwigHadronizer.h.

Referenced by readSettings().

int gen::PomwigHadronizer::maxEventsToPrint

private

Definition at line 55 of file PomwigHadronizer.h.

Referenced by readSettings().

bool gen::PomwigHadronizer::needClear

private

Definition at line 50 of file PomwigHadronizer.h.

Referenced by clear(), and readSettings().

int gen::PomwigHadronizer::numTrials

private

Definition at line 65 of file PomwigHadronizer.h.

gen::ParameterCollector gen::PomwigHadronizer::parameters

private

Definition at line 52 of file PomwigHadronizer.h.

Referenced by readSettings().

bool gen::PomwigHadronizer::printCards

private

Definition at line 56 of file PomwigHadronizer.h.

double gen::PomwigHadronizer::survivalProbability

private

Definition at line 59 of file PomwigHadronizer.h.

Referenced by statistics().

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

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

Member Function Documentation

Member Data Documentation