gen::BaseHadronizer Class Reference

#include <BaseHadronizer.h>

Inheritance diagram for gen::BaseHadronizer:

Public Member Functions
	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)

Protected Member Functions
std::unique_ptr< HepMC::GenEvent > &	event ()
std::unique_ptr< GenEventInfoProduct > &	eventInfo ()
lhef::LHEEvent *	lheEvent ()
lhef::LHERunInfo *	lheRunInfo ()
GenRunInfoProduct &	runInfo ()

Protected Attributes
std::string	lheFile_
int	randomIndex_

Private Member Functions
virtual void	doSetRandomEngine (CLHEP::HepRandomEngine *v)
virtual std::vector< std::string > const &	doSharedResources () const

Private Attributes
edm::Event *	edmEvent_
std::unique_ptr< HepMC::GenEvent >	genEvent_
std::unique_ptr< GenEventInfoProduct >	genEventInfo_
GenRunInfoProduct	genRunInfo_
std::vector< std::string >	gridpackPaths_
std::unique_ptr< lhef::LHEEvent >	lheEvent_
std::shared_ptr< lhef::LHERunInfo >	lheRunInfo_
std::vector< std::string >	randomInitConfigDescriptions_
std::vector< double >	randomInitWeights_

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

Detailed Description

Definition at line 47 of file BaseHadronizer.h.

Constructor & Destructor Documentation

gen::BaseHadronizer::BaseHadronizer

(

edm::ParameterSet const &

ps

)

Definition at line 12 of file BaseHadronizer.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gridpackPaths_, randomInitConfigDescriptions_, randomInitWeights_, runInfo(), GenRunInfoProduct::setExternalXSecLO(), GenRunInfoProduct::setExternalXSecNLO(), GenRunInfoProduct::setFilterEfficiency(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                            :
  randomIndex_(-1),
  gridpackPaths_(1)
  {
    
    if (ps.exists("RandomizedParameters")) {
      std::vector<edm::ParameterSet> randomizedParameters = ps.getParameter<std::vector<edm::ParameterSet> >("RandomizedParameters");
      randomInitWeights_.resize(randomizedParameters.size());
      randomInitConfigDescriptions_.resize(randomizedParameters.size());
      gridpackPaths_.resize(randomizedParameters.size());
      for (unsigned int irand = 0; irand<randomizedParameters.size(); ++irand) {
        randomInitWeights_[irand] = randomizedParameters[irand].getParameter<double>("ConfigWeight");
        if (randomizedParameters[irand].exists("ConfigDescription")) {
          randomInitConfigDescriptions_[irand] = randomizedParameters[irand].getParameter<std::string>("ConfigDescription");
        }
        if (randomizedParameters[irand].exists("GridpackPath")) {
          gridpackPaths_[irand] = randomizedParameters[irand].getParameter<std::string>("GridpackPath");
        }
      }
    }
    else {
      if (ps.exists("GridpackPath")) {
        gridpackPaths_[0] = ps.getParameter<std::string>("GridpackPath");
      }
    }
    
    runInfo().setFilterEfficiency(
        ps.getUntrackedParameter<double>("filterEfficiency", -1.) );
    runInfo().setExternalXSecLO(
        GenRunInfoProduct::XSec(ps.getUntrackedParameter<double>("crossSection", -1.)) );
    runInfo().setExternalXSecNLO(
        GenRunInfoProduct::XSec(ps.getUntrackedParameter<double>("crossSectionNLO", -1.)) );

  }

virtual gen::BaseHadronizer::~BaseHadronizer ( )

inlinevirtual

Definition at line 50 of file BaseHadronizer.h.

{}

Member Function Documentation

void gen::BaseHadronizer::cleanLHE

(

)

Definition at line 146 of file BaseHadronizer.cc.

References lheFile_, and MatrixUtil::remove().

                                {
    if (lheFile_.empty()) {
      return;
    }
    
    std::remove(lheFile_.c_str());
  }

virtual void gen::BaseHadronizer::doSetRandomEngine ( CLHEP::HepRandomEngine *  v )

inlineprivatevirtual

Reimplemented in gen::ReggeGribovPartonMCHadronizer, Herwig6Hadronizer, gen::Hydjet2Hadronizer, Pythia8Hadronizer, gen::Pythia6Hadronizer, gen::ExhumeHadronizer, Herwig7Hadronizer, SherpaHadronizer, gen::HydjetHadronizer, gen::AMPTHadronizer, gen::Cascade2Hadronizer, gen::PyquenHadronizer, gen::HijingHadronizer, and gen::PomwigHadronizer.

Definition at line 97 of file BaseHadronizer.h.

{ }

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

inlineprivatevirtual

Reimplemented in Herwig6Hadronizer, gen::Hydjet2Hadronizer, Pythia8Hadronizer, gen::Pythia6Hadronizer, gen::ExhumeHadronizer, gen::HydjetHadronizer, gen::Cascade2Hadronizer, gen::PyquenHadronizer, gen::HijingHadronizer, and gen::PomwigHadronizer.

Definition at line 99 of file BaseHadronizer.h.

{ return theSharedResources; }

std::unique_ptr<HepMC::GenEvent>& gen::BaseHadronizer::event ( )

inlineprotected

Definition at line 87 of file BaseHadronizer.h.

Referenced by gen::PomwigHadronizer::decay(), gen::PomwigHadronizer::finalizeEvent(), gen::Cascade2Hadronizer::finalizeEvent(), gen::Py8GunBase::finalizeEvent(), Herwig7Hadronizer::finalizeEvent(), SherpaHadronizer::finalizeEvent(), gen::ExhumeHadronizer::finalizeEvent(), gen::Pythia6Hadronizer::finalizeEvent(), Pythia8Hadronizer::finalizeEvent(), Herwig6Hadronizer::finalizeEvent(), gen::Py8PtAndDxyGun::generatePartonsAndHadronize(), gen::Py8PtGun::generatePartonsAndHadronize(), gen::Py8EGun::generatePartonsAndHadronize(), gen::Py8JetGun::generatePartonsAndHadronize(), gen::HijingHadronizer::generatePartonsAndHadronize(), gen::PomwigHadronizer::generatePartonsAndHadronize(), gen::PyquenHadronizer::generatePartonsAndHadronize(), gen::AMPTHadronizer::generatePartonsAndHadronize(), gen::HydjetHadronizer::generatePartonsAndHadronize(), gen::Cascade2Hadronizer::generatePartonsAndHadronize(), Herwig7Hadronizer::generatePartonsAndHadronize(), gen::ExhumeHadronizer::generatePartonsAndHadronize(), gen::Pythia6Hadronizer::generatePartonsAndHadronize(), gen::Hydjet2Hadronizer::generatePartonsAndHadronize(), Pythia8Hadronizer::generatePartonsAndHadronize(), gen::ReggeGribovPartonMCHadronizer::generatePartonsAndHadronize(), Herwig7Hadronizer::hadronize(), gen::Pythia6Hadronizer::hadronize(), Pythia8Hadronizer::hadronize(), Herwig6Hadronizer::hadronize(), gen::Cascade2Hadronizer::imposeProperTime(), gen::Pythia6Hadronizer::imposeProperTime(), Herwig6Hadronizer::pythiaStatusCodes(), gen::Cascade2Hadronizer::residualDecay(), gen::Py8GunBase::residualDecay(), gen::Pythia6Hadronizer::residualDecay(), and Pythia8Hadronizer::residualDecay().

{ return genEvent_; }

std::unique_ptr<GenEventInfoProduct>& gen::BaseHadronizer::eventInfo ( )

inlineprotected

Definition at line 88 of file BaseHadronizer.h.

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

{ return genEventInfo_; }

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

Definition at line 78 of file BaseHadronizer.cc.

References writedatasetfile::args, begin, constexpr, relativeConstraints::empty, Exception, edm::FileInPath::fullPath(), edm::pset::Registry::getMapped(), gridpackPath(), edm::LuminosityBlockBase::id(), edm::pset::Registry::instance(), lheFile_, edm::LuminosityBlockID::luminosityBlock(), sysUtil::pid, edm::LuminosityBlock::processHistory(), edm::ProcessHistory::rbegin(), edm::LuminosityBlockID::run(), and mps_update::status.

                                                                                                                   {
        
    if (gridpackPath().empty()) {
      return;
    }
    
    //get random seed from HepRandomEngine state plus run and lumi section numbers
    //(random engine state will be the same for every lumi section in a job)
    std::vector<long unsigned int> seeds = rengine->put();
    seeds.push_back(lumi.id().run());
    seeds.push_back(lumi.id().luminosityBlock());
    std::seed_seq seedseq(seeds.begin(),seeds.end());
    std::array<unsigned int,1> lheseed;
    seedseq.generate(lheseed.begin(),lheseed.end());
    
    constexpr unsigned int maxseed = 30081*30081; //madgraph cannot handle seeds larger than this
    unsigned int seedval = lheseed[0]%(maxseed+1);

    unsigned int nevents = edm::pset::Registry::instance()->getMapped(lumi.processHistory().rbegin()->parameterSetID())->getParameter<edm::ParameterSet>("@main_input").getUntrackedParameter<unsigned int>("numberEventsInLuminosityBlock");
    
    std::ostringstream nevStream;
    nevStream << nevents;
    
    std::ostringstream randomStream;
    randomStream << seedval;
   
    edm::FileInPath script("GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh");
    const char *outfilename = "cmsgrid_final.lhe";


    
    std::array<std::string,5> argStrs;
    argStrs[0]=script.fullPath();
    argStrs[1]=gridpackPath();
    argStrs[2]=nevStream.str();
    argStrs[3]=randomStream.str();
    argStrs[4]=std::to_string(ncpu);
    std::array<char*,6>args{ { &argStrs[0][0], &argStrs[1][0],&argStrs[2][0],&argStrs[3][0], &argStrs[4][0], nullptr } };
    
    pid_t pid = fork();

    if (pid == -1) {
        // error, failed to fork()
      throw cms::Exception("BaseHadronizer::generateLHE") << "Unable to fork a child";
    } 
    else if (pid==0) {
      //child
      execvp(args[0],std::begin(args));
      _exit(1);   // exec never returns
    }
    else {
      //parent
      int status;
      waitpid(pid, &status, 0);
      if (status) {
        throw cms::Exception("BaseHadronizer::generateLHE") << "Failed to execute script";
      }
    }
    FILE* lhef = std::fopen(outfilename, "r");
    if (!lhef) {
      throw cms::Exception("BaseHadronizer::generateLHE") << "Output file " << outfilename << " not found.";
    }
    std::fclose(lhef);

    lheFile_ = outfilename;
    
  }

edm::Event& gen::BaseHadronizer::getEDMEvent ( ) const

inline

Definition at line 70 of file BaseHadronizer.h.

Referenced by gen::PyquenHadronizer::generatePartonsAndHadronize(), gen::HydjetHadronizer::generatePartonsAndHadronize(), and gen::Hydjet2Hadronizer::generatePartonsAndHadronize().

{ return *edmEvent_; }

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

inline

Definition at line 54 of file BaseHadronizer.h.

References eostools::move().

{ return std::move(genEvent_); }

std::unique_ptr<GenEventInfoProduct> gen::BaseHadronizer::getGenEventInfo ( )

inline

Definition at line 55 of file BaseHadronizer.h.

References eostools::move().

{ return std::move(genEventInfo_); }

std::unique_ptr< GenLumiInfoHeader > gen::BaseHadronizer::getGenLumiInfoHeader ( ) const

virtual

Reimplemented in Pythia8Hadronizer, Herwig7Hadronizer, and SherpaHadronizer.

Definition at line 47 of file BaseHadronizer.cc.

References genWeightsTable_cfi::genLumiInfoHeader, randomIndex_, and randomInitConfigDescription().

                                                                              {
    
    auto genLumiInfoHeader = std::make_unique<GenLumiInfoHeader>();
    
    //fill information on randomized configs for parameter scans
    genLumiInfoHeader->setRandomConfigIndex(randomIndex_);
    if (randomIndex_>=0) {
      genLumiInfoHeader->setConfigDescription(randomInitConfigDescription());      
    }
    
    return genLumiInfoHeader;
    
  }

GenRunInfoProduct& gen::BaseHadronizer::getGenRunInfo ( )

inline

Definition at line 53 of file BaseHadronizer.h.

{ return genRunInfo_; }

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

inline

Definition at line 57 of file BaseHadronizer.h.

References eostools::move().

{ return std::move(lheEvent_);}

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

inline

Definition at line 63 of file BaseHadronizer.h.

{ return lheRunInfo_; }

const std::string& gen::BaseHadronizer::gridpackPath ( ) const

inline

Definition at line 79 of file BaseHadronizer.h.

References SiStripPI::max.

Referenced by generateLHE().

{ return gridpackPaths_[std::max(randomIndex_,0)]; }

lhef::LHEEvent* gen::BaseHadronizer::lheEvent ( )

inlineprotected

Definition at line 90 of file BaseHadronizer.h.

Referenced by gen::Pythia6Hadronizer::finalizeEvent(), Pythia8Hadronizer::finalizeEvent(), Herwig6Hadronizer::finalizeEvent(), Pythia8Hadronizer::generatePartonsAndHadronize(), Herwig7Hadronizer::hadronize(), gen::Pythia6Hadronizer::hadronize(), Pythia8Hadronizer::hadronize(), Herwig6Hadronizer::hadronize(), and Herwig6Hadronizer::upEvnt().

{ return lheEvent_.get(); }

lhef::LHERunInfo* gen::BaseHadronizer::lheRunInfo ( )

inlineprotected

Definition at line 91 of file BaseHadronizer.h.

Referenced by gen::Pythia6Hadronizer::finalizeEvent(), Herwig6Hadronizer::finalizeEvent(), Herwig6Hadronizer::initialize(), gen::Pythia6Hadronizer::initializeForExternalPartons(), Pythia8Hadronizer::initializeForExternalPartons(), Herwig6Hadronizer::readSettings(), and Herwig6Hadronizer::upInit().

{ return lheRunInfo_.get(); }

int gen::BaseHadronizer::randomIndex ( ) const

inline

Definition at line 77 of file BaseHadronizer.h.

Referenced by gen::Py8InterfaceBase::readSettings().

{ return randomIndex_; }

const std::string& gen::BaseHadronizer::randomInitConfigDescription ( ) const

inline

Definition at line 78 of file BaseHadronizer.h.

Referenced by getGenLumiInfoHeader().

{ return randomInitConfigDescriptions_[randomIndex_]; }

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

Definition at line 61 of file BaseHadronizer.cc.

References edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), randomIndex_, randomInitWeights_, and edm::LuminosityBlockID::run().

                                                                                                   {
    if (!randomInitWeights_.empty()) {
      //randomly select from a list of provided configuration sets (for parameter scans)
      
      //seeds std 32-bit mersene twister with HepRandomEngine state plus run and lumi section numbers
      //(random engine state will be the same for every lumi section in a job)
      std::vector<long unsigned int> seeds = rengine->put();
      seeds.push_back(lumi.id().run());
      seeds.push_back(lumi.id().luminosityBlock());
      std::seed_seq seedseq(seeds.begin(),seeds.end());
      std::mt19937 randgen(seedseq);
      std::discrete_distribution<int> randdist(randomInitWeights_.begin(),randomInitWeights_.end());
      
      randomIndex_ = randdist(randgen);
    }
  }

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

inline

Definition at line 59 of file BaseHadronizer.h.

References eostools::move().

Referenced by SherpaHadronizer::generatePartonsAndHadronize().

{ genEvent_ = std::move(event); }

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

inline

Definition at line 60 of file BaseHadronizer.h.

References eostools::move().

{ genEventInfo_ = std::move(eventInfo); }

GenRunInfoProduct& gen::BaseHadronizer::runInfo ( )

inlineprotected

Definition at line 86 of file BaseHadronizer.h.

Referenced by BaseHadronizer(), gen::PomwigHadronizer::statistics(), gen::Cascade2Hadronizer::statistics(), Herwig7Hadronizer::statistics(), gen::Py8GunBase::statistics(), SherpaHadronizer::statistics(), gen::ExhumeHadronizer::statistics(), gen::Pythia6Hadronizer::statistics(), Pythia8Hadronizer::statistics(), and Herwig6Hadronizer::statistics().

{ return genRunInfo_; }

virtual bool gen::BaseHadronizer::select ( HepMC::GenEvent *  ) const

inlinevirtual

Reimplemented in gen::PyquenHadronizer.

Definition at line 71 of file BaseHadronizer.h.

{ return true;}

void gen::BaseHadronizer::setEDMEvent ( edm::Event &  event )

inline

Definition at line 69 of file BaseHadronizer.h.

References event().

{ edmEvent_ = &event; }

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

inline

Definition at line 66 of file BaseHadronizer.h.

References eostools::move().

{ lheEvent_ = std::move(event); }

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

inline

Definition at line 65 of file BaseHadronizer.h.

References eostools::move().

{ lheRunInfo_ = std::move(runInfo); }

void gen::BaseHadronizer::setRandomEngine ( CLHEP::HepRandomEngine *  v )

inline

Definition at line 73 of file BaseHadronizer.h.

Referenced by CMS_SHERPA_RNG::CMS_SHERPA_RNG().

{ doSetRandomEngine(v); }

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

inline

Definition at line 75 of file BaseHadronizer.h.

{ return doSharedResources(); }

Member Data Documentation

edm::Event* gen::BaseHadronizer::edmEvent_

private

Definition at line 108 of file BaseHadronizer.h.

std::unique_ptr<HepMC::GenEvent> gen::BaseHadronizer::genEvent_

private

Definition at line 102 of file BaseHadronizer.h.

std::unique_ptr<GenEventInfoProduct> gen::BaseHadronizer::genEventInfo_

private

Definition at line 103 of file BaseHadronizer.h.

GenRunInfoProduct gen::BaseHadronizer::genRunInfo_

private

Definition at line 101 of file BaseHadronizer.h.

std::vector<std::string> gen::BaseHadronizer::gridpackPaths_

private

Definition at line 114 of file BaseHadronizer.h.

Referenced by BaseHadronizer().

std::unique_ptr<lhef::LHEEvent> gen::BaseHadronizer::lheEvent_

private

Definition at line 106 of file BaseHadronizer.h.

std::string gen::BaseHadronizer::lheFile_

protected

Definition at line 93 of file BaseHadronizer.h.

Referenced by cleanLHE(), generateLHE(), and Pythia8Hadronizer::initializeForInternalPartons().

std::shared_ptr<lhef::LHERunInfo> gen::BaseHadronizer::lheRunInfo_

private

Definition at line 105 of file BaseHadronizer.h.

int gen::BaseHadronizer::randomIndex_

protected

Definition at line 92 of file BaseHadronizer.h.

Referenced by getGenLumiInfoHeader(), and randomizeIndex().

std::vector<std::string> gen::BaseHadronizer::randomInitConfigDescriptions_

private

Definition at line 113 of file BaseHadronizer.h.

Referenced by BaseHadronizer().

std::vector<double> gen::BaseHadronizer::randomInitWeights_

private

Definition at line 112 of file BaseHadronizer.h.

Referenced by BaseHadronizer(), and randomizeIndex().

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

staticprivate

Definition at line 110 of file BaseHadronizer.h.