#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 ()(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 21 of file PomwigHadronizer.h.

Constructor & Destructor Documentation

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

Definition at line 94 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 115 of file PomwigHadronizer.cc.

References clear().

 {
     clear();
 }

Member Function Documentation

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

inline

Definition at line 41 of file PomwigHadronizer.h.

References findQualityFiles::v.

41 { return "PomwigHadronizer"; }

void gen::PomwigHadronizer::clear ( void )

private

Definition at line 125 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 358 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().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;
 }

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

inline

Definition at line 32 of file PomwigHadronizer.h.

References SelectingProcedure_cff::decay, and myMessageLogger_cff::statistics.

32 { return true; }

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

Definition at line 289 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 120 of file PomwigHadronizer.cc.

References gen::Herwig6Instance::setHerwigRandomEngine().

 {
         setHerwigRandomEngine(v);
 }

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

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 46 of file PomwigHadronizer.h.

References hitfit::clear().

46 { return theSharedResources; }

gen::PomwigHadronizer::theSharedResources

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

Definition: PomwigHadronizer.h:51

void gen::PomwigHadronizer::finalizeEvent ( )

Definition at line 349 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 314 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 309 of file PomwigHadronizer.cc.

 {
    return false;
 }

bool gen::PomwigHadronizer::initializeDPDF ( )

private

Definition at line 243 of file PomwigHadronizer.cc.

References printConversionInfo::aux, GetRecoTauVFromDQM_MC_cff::cl, edm::errors::Configuration, diffTopology, Exception, connectstrParser::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 138 of file PomwigHadronizer.cc.

 {
    return false;
 }

bool gen::PomwigHadronizer::initializeForInternalPartons ( )

Definition at line 229 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 143 of file PomwigHadronizer.cc.

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.

 {
 
    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 389 of file PomwigHadronizer.cc.

 {
     return true;
 }

void gen::PomwigHadronizer::statistics ( )

Definition at line 298 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 61 of file PomwigHadronizer.h.

Referenced by readSettings().

HepMC::IO_HERWIG gen::PomwigHadronizer::conv

private

Definition at line 72 of file PomwigHadronizer.h.

Referenced by decay().

int gen::PomwigHadronizer::diffTopology

private

Definition at line 63 of file PomwigHadronizer.h.

Referenced by initializeDPDF(), and readSettings().

bool gen::PomwigHadronizer::doMPInteraction

private

Definition at line 67 of file PomwigHadronizer.h.

Referenced by generatePartonsAndHadronize().

bool gen::PomwigHadronizer::doPDGConvert

private

Definition at line 70 of file PomwigHadronizer.h.

Referenced by decay(), and PomwigHadronizer().

int gen::PomwigHadronizer::h1fit

private

Definition at line 64 of file PomwigHadronizer.h.

Referenced by initializeDPDF().

int gen::PomwigHadronizer::hepmcVerbosity

private

Definition at line 57 of file PomwigHadronizer.h.

int gen::PomwigHadronizer::herwigVerbosity

private

Definition at line 56 of file PomwigHadronizer.h.

Referenced by readSettings().

int gen::PomwigHadronizer::maxEventsToPrint

private

Definition at line 58 of file PomwigHadronizer.h.

Referenced by readSettings().

bool gen::PomwigHadronizer::needClear

private

Definition at line 53 of file PomwigHadronizer.h.

Referenced by clear(), and readSettings().

int gen::PomwigHadronizer::numTrials

private

Definition at line 68 of file PomwigHadronizer.h.

gen::ParameterCollector gen::PomwigHadronizer::parameters

private

Definition at line 55 of file PomwigHadronizer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

bool gen::PomwigHadronizer::printCards

private

Definition at line 59 of file PomwigHadronizer.h.

double gen::PomwigHadronizer::survivalProbability

private

Definition at line 62 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 51 of file PomwigHadronizer.h.

bool gen::PomwigHadronizer::useJimmy

private

Definition at line 66 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