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)
edm::Event &	getEDMEvent () const
HepMC::GenEvent *	getGenEvent ()
GenEventInfoProduct *	getGenEventInfo ()
virtual GenLumiInfoHeader *	getGenLumiInfoHeader () const
GenRunInfoProduct &	getGenRunInfo ()
const boost::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 (HepMC::GenEvent *event)
void	resetEventInfo (GenEventInfoProduct *eventInfo)
virtual bool	select (HepMC::GenEvent *) const
void	setEDMEvent (edm::Event &event)
void	setLHEEvent (lhef::LHEEvent *event)
void	setLHERunInfo (lhef::LHERunInfo *runInfo)
void	setRandomEngine (CLHEP::HepRandomEngine *v)
std::vector< std::string > const &	sharedResources () const
	~BaseHadronizer ()

Protected Member Functions
std::auto_ptr< HepMC::GenEvent > &	event ()
std::auto_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::auto_ptr< HepMC::GenEvent >	genEvent_
std::auto_ptr < GenEventInfoProduct >	genEventInfo_
GenRunInfoProduct	genRunInfo_
std::vector< std::string >	gridpackPaths_
std::auto_ptr< lhef::LHEEvent >	lheEvent_
boost::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.)) );

  }

gen::BaseHadronizer::~BaseHadronizer ( )

inline

Definition at line 50 of file BaseHadronizer.h.

{}

Member Function Documentation

void gen::BaseHadronizer::cleanLHE

(

)

Definition at line 147 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, gen::HydjetHadronizer, SherpaHadronizer, gen::AMPTHadronizer, gen::Cascade2Hadronizer, gen::PyquenHadronizer, ThePEGHadronizer, gen::HijingHadronizer, and gen::PomwigHadronizer.

Definition at line 96 of file BaseHadronizer.h.

Referenced by setRandomEngine().

{ }

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 98 of file BaseHadronizer.h.

References theSharedResources.

Referenced by sharedResources().

{ return theSharedResources; }

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

inlineprotected

Definition at line 86 of file BaseHadronizer.h.

References genEvent_.

Referenced by gen::Py8GunBase::finalizeEvent(), gen::Py8PtGun::generatePartonsAndHadronize(), gen::Py8JetGun::generatePartonsAndHadronize(), gen::Py8EGun::generatePartonsAndHadronize(), gen::Py8GunBase::residualDecay(), and setEDMEvent().

{ return genEvent_; }

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

inlineprotected

Definition at line 87 of file BaseHadronizer.h.

References genEventInfo_.

{ return genEventInfo_; }

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

Definition at line 78 of file BaseHadronizer.cc.

References createfilelist::args, constexpr, relativeConstraints::empty, Exception, edm::FileInPath::fullPath(), edm::pset::Registry::getMapped(), gridpackPath(), edm::LuminosityBlockBase::id(), edm::pset::Registry::instance(), lheFile_, edm::LuminosityBlockID::luminosityBlock(), cmsHarvester::nevents, NULL, EcalCondDB::outfilename, 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";
    
    char *args[5];
    args[0] = strdup(script.fullPath().c_str());
    args[1] = strdup(gridpackPath().c_str());
    args[2] = strdup(nevStream.str().c_str());
    args[3] = strdup(randomStream.str().c_str());
    args[4] = NULL;
    
    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],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;
    
    for (int iarg=0; iarg<4; ++iarg) {
      delete[] args[iarg];
    }
    
  }

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

inline

Definition at line 69 of file BaseHadronizer.h.

References edmEvent_.

{ return *edmEvent_; }

HepMC::GenEvent* gen::BaseHadronizer::getGenEvent ( )

inline

Definition at line 54 of file BaseHadronizer.h.

References genEvent_.

{ return genEvent_.release(); }

GenEventInfoProduct* gen::BaseHadronizer::getGenEventInfo ( )

inline

Definition at line 55 of file BaseHadronizer.h.

References genEventInfo_.

{ return genEventInfo_.release(); }

GenLumiInfoHeader * gen::BaseHadronizer::getGenLumiInfoHeader ( ) const

virtual

Reimplemented in Pythia8Hadronizer.

Definition at line 47 of file BaseHadronizer.cc.

References randomIndex_, randomInitConfigDescription(), GenLumiInfoHeader::setConfigDescription(), and GenLumiInfoHeader::setRandomConfigIndex().

                                                                {
    
    GenLumiInfoHeader *genLumiInfoHeader = new 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.

References genRunInfo_.

{ return genRunInfo_; }

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

inline

Definition at line 62 of file BaseHadronizer.h.

References lheRunInfo_.

{ return lheRunInfo_; }

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

inline

Definition at line 78 of file BaseHadronizer.h.

References gridpackPaths_, bookConverter::max, and randomIndex_.

Referenced by generateLHE().

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

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

inlineprotected

Definition at line 89 of file BaseHadronizer.h.

References lheEvent_.

{ return lheEvent_.get(); }

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

inlineprotected

Definition at line 90 of file BaseHadronizer.h.

References lheRunInfo_.

{ return lheRunInfo_.get(); }

int gen::BaseHadronizer::randomIndex ( ) const

inline

Definition at line 76 of file BaseHadronizer.h.

References randomIndex_.

{ return randomIndex_; }

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

inline

Definition at line 77 of file BaseHadronizer.h.

References randomIndex_, and randomInitConfigDescriptions_.

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_.size()>0) {
      //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 ( HepMC::GenEvent *  event )

inline

Definition at line 58 of file BaseHadronizer.h.

References genEvent_.

{ genEvent_.reset(event); }

void gen::BaseHadronizer::resetEventInfo ( GenEventInfoProduct *  eventInfo )

inline

Definition at line 59 of file BaseHadronizer.h.

References genEventInfo_.

{ genEventInfo_.reset(eventInfo); }

GenRunInfoProduct& gen::BaseHadronizer::runInfo ( )

inlineprotected

Definition at line 85 of file BaseHadronizer.h.

References genRunInfo_.

Referenced by BaseHadronizer(), and gen::Py8GunBase::statistics().

{ return genRunInfo_; }

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

inlinevirtual

Reimplemented in gen::PyquenHadronizer.

Definition at line 70 of file BaseHadronizer.h.

{ return true;}

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

inline

Definition at line 68 of file BaseHadronizer.h.

References edmEvent_, and event().

{ edmEvent_ = &event; }

void gen::BaseHadronizer::setLHEEvent ( lhef::LHEEvent *  event )

inline

Definition at line 65 of file BaseHadronizer.h.

References lheEvent_.

{ lheEvent_.reset(event); }

void gen::BaseHadronizer::setLHERunInfo ( lhef::LHERunInfo *  runInfo )

inline

Definition at line 64 of file BaseHadronizer.h.

References lheRunInfo_.

{ lheRunInfo_.reset(runInfo); }

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

inline

Definition at line 72 of file BaseHadronizer.h.

References doSetRandomEngine().

{ doSetRandomEngine(v); }

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

inline

Definition at line 74 of file BaseHadronizer.h.

References doSharedResources().

{ return doSharedResources(); }

Member Data Documentation

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

private

Definition at line 107 of file BaseHadronizer.h.

Referenced by getEDMEvent(), and setEDMEvent().

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

private

Definition at line 101 of file BaseHadronizer.h.

Referenced by event(), getGenEvent(), and resetEvent().

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

private

Definition at line 102 of file BaseHadronizer.h.

Referenced by eventInfo(), getGenEventInfo(), and resetEventInfo().

GenRunInfoProduct gen::BaseHadronizer::genRunInfo_

private

Definition at line 100 of file BaseHadronizer.h.

Referenced by getGenRunInfo(), and runInfo().

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

private

Definition at line 113 of file BaseHadronizer.h.

Referenced by BaseHadronizer(), and gridpackPath().

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

private

Definition at line 105 of file BaseHadronizer.h.

Referenced by lheEvent(), and setLHEEvent().

std::string gen::BaseHadronizer::lheFile_

protected

Definition at line 92 of file BaseHadronizer.h.

Referenced by cleanLHE(), and generateLHE().

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

private

Definition at line 104 of file BaseHadronizer.h.

Referenced by getLHERunInfo(), lheRunInfo(), and setLHERunInfo().

int gen::BaseHadronizer::randomIndex_

protected

Definition at line 91 of file BaseHadronizer.h.

Referenced by getGenLumiInfoHeader(), gridpackPath(), randomIndex(), randomInitConfigDescription(), and randomizeIndex().

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

private

Definition at line 112 of file BaseHadronizer.h.

Referenced by BaseHadronizer(), and randomInitConfigDescription().

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

private

Definition at line 111 of file BaseHadronizer.h.

Referenced by BaseHadronizer(), and randomizeIndex().

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

staticprivate

Definition at line 109 of file BaseHadronizer.h.

Referenced by doSharedResources().