#include <MCatNLOSource.h>

Inheritance diagram for MCatNLOSource:

Public Member Functions
	MCatNLOSource (const edm::ParameterSet &params, const edm::InputSourceDescription &desc)

virtual	~MCatNLOSource ()

Public Member Functions inherited from edm::ProducerSourceFromFiles
bool	noFiles () const override

	ProducerSourceFromFiles (ParameterSet const &pset, InputSourceDescription const &desc, bool realData)

	~ProducerSourceFromFiles () override

Public Member Functions inherited from edm::ProducerSourceBase
EventNumber_t	event () const

unsigned int	eventCreationDelay () const

EventID const &	eventID () const

LuminosityBlockNumber_t	luminosityBlock () const

unsigned int	numberEventsInLumi () const

unsigned int	numberEventsInRun () const

unsigned int	numberEventsInThisLumi () const

unsigned int	numberEventsInThisRun () const

TimeValue_t	presentTime () const

	ProducerSourceBase (ParameterSet const &pset, InputSourceDescription const &desc, bool realData)

RunNumber_t	run () const

unsigned int	timeBetweenEvents () const

	~ProducerSourceBase () noexcept(false) override

Public Member Functions inherited from edm::PuttableSourceBase
	PuttableSourceBase (ParameterSet const &, InputSourceDescription const &)

void	registerProducts () final
	Register any produced products. More...

Public Member Functions inherited from edm::InputSource
std::shared_ptr< ActivityRegistry >	actReg () const
	Accessor for Activity Registry. More...

std::shared_ptr< BranchIDListHelper const >	branchIDListHelper () const
	Accessors for branchIDListHelper. More...

std::shared_ptr< BranchIDListHelper > &	branchIDListHelper ()

void	closeFile (FileBlock *, bool cleaningUpAfterException)
	close current file More...

void	doBeginJob ()
	Called by framework at beginning of job. More...

void	doEndJob ()
	Called by framework at end of job. More...

ProcessingController::ForwardState	forwardState () const

bool	goToEvent (EventID const &eventID)

	InputSource (ParameterSet const &, InputSourceDescription const &)
	Constructor. More...

	InputSource (InputSource const &)=delete

void	issueReports (EventID const &eventID)
	issue an event report More...

LuminosityBlockNumber_t	luminosityBlock () const
	Accessor for current luminosity block number. More...

std::shared_ptr< LuminosityBlockAuxiliary >	luminosityBlockAuxiliary () const
	Called by the framework to merge or insert lumi in principal cache. More...

int	maxEvents () const

int	maxLuminosityBlocks () const

ModuleDescription const &	moduleDescription () const
	Accessor for 'module' description. More...

ItemType	nextItemType ()
	Advances the source to the next item. More...

InputSource &	operator= (InputSource const &)=delete

ProcessConfiguration const &	processConfiguration () const
	Accessor for Process Configuration. More...

std::string const &	processGUID () const
	Accessor for global process identifier. More...

ProcessHistoryRegistry const &	processHistoryRegistry () const
	Accessors for process history registry. More...

ProcessHistoryRegistry &	processHistoryRegistry ()

ProcessingMode	processingMode () const
	RunsLumisAndEvents (default), RunsAndLumis, or Runs. More...

std::shared_ptr< ProductRegistry const >	productRegistry () const
	Accessors for product registry. More...

std::shared_ptr< ProductRegistry > &	productRegistry ()

bool	randomAccess () const

void	readAndMergeLumi (LuminosityBlockPrincipal &lbp)
	Read next luminosity block (same as a prior lumi) More...

void	readAndMergeRun (RunPrincipal &rp)
	Read next run (same as a prior run) More...

void	readEvent (EventPrincipal &ep, StreamContext &)
	Read next event. More...

bool	readEvent (EventPrincipal &ep, EventID const &, StreamContext &)
	Read a specific event. More...

std::unique_ptr< FileBlock >	readFile ()
	Read next file. More...

void	readLuminosityBlock (LuminosityBlockPrincipal &lumiPrincipal, HistoryAppender &historyAppender)
	Read next luminosity block (new lumi) More...

std::shared_ptr< LuminosityBlockAuxiliary >	readLuminosityBlockAuxiliary ()
	Read next luminosity block Auxilary. More...

void	readRun (RunPrincipal &runPrincipal, HistoryAppender &historyAppender)
	Read next run (new run) More...

std::shared_ptr< RunAuxiliary >	readRunAuxiliary ()
	Read next run Auxiliary. More...

ProcessHistoryID const &	reducedProcessHistoryID () const

int	remainingEvents () const

int	remainingLuminosityBlocks () const

void	repeat ()
	Reset the remaining number of events/lumis to the maximum number. More...

std::pair< SharedResourcesAcquirer , std::recursive_mutex >	resourceSharedWithDelayedReader ()
	Returns nullptr if no resource shared between the Source and a DelayedReader. More...

ProcessingController::ReverseState	reverseState () const

void	rewind ()
	Begin again at the first event. More...

RunNumber_t	run () const
	Accessor for current run number. More...

std::shared_ptr< RunAuxiliary >	runAuxiliary () const
	Called by the framework to merge or insert run in principal cache. More...

void	setLuminosityBlockNumber_t (LuminosityBlockNumber_t lb)
	Set the luminosity block ID. More...

void	setRunNumber (RunNumber_t r)
	Set the run number. More...

void	skipEvents (int offset)

std::shared_ptr< ThinnedAssociationsHelper const >	thinnedAssociationsHelper () const
	Accessors for thinnedAssociationsHelper. More...

std::shared_ptr< ThinnedAssociationsHelper > &	thinnedAssociationsHelper ()

Timestamp const &	timestamp () const
	Accessor for the current time, as seen by the input source. More...

virtual	~InputSource () noexcept(false)
	Destructor. More...

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () 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)

virtual	~Herwig6Instance ()

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
virtual void	beginRun (edm::Run &run)

virtual void	endJob ()

virtual bool	hwwarn (const std::string &fn, int code)

void	nextEvent ()

virtual void	produce (edm::Event &event)

virtual bool	setRunAndEventInfo (edm::EventID &, edm::TimeValue_t &, edm::EventAuxiliary::ExperimentType &)

Private Attributes
boost::shared_ptr< lhef::LHEEvent >	event

std::string	fileName
	Name of the input file. More...

int	ihpro

std::unique_ptr< std::ifstream >	inputFile
	Pointer to the input file. More...

unsigned int	nEvents
	Number of events. More...

int	processCode

std::unique_ptr< std::ifstream >	reader

boost::shared_ptr< lhef::LHERunInfo >	runInfo

unsigned int	skipEvents
	Number of events to skip. More...

Additional Inherited Members
Public Types inherited from edm::InputSource
enum	ItemType { IsInvalid, IsStop, IsFile, IsRun, IsLumi, IsEvent, IsRepeat, IsSynchronize }

enum	ProcessingMode { Runs, RunsAndLumis, RunsLumisAndEvents }

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::ProducerSourceFromFiles
static void	fillDescription (ParameterSetDescription &desc)

Static Public Member Functions inherited from edm::ProducerSourceBase
static void	fillDescription (ParameterSetDescription &desc)

Static Public Member Functions inherited from edm::InputSource
static const std::string &	baseType ()

static void	fillDescription (ParameterSetDescription &desc)

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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

Public Attributes inherited from edm::InputSource
signalslot::Signal< void(StreamContext const &, ModuleCallingContext const &)>	postEventReadFromSourceSignal_

signalslot::Signal< void(StreamContext const &, ModuleCallingContext const &)>	preEventReadFromSourceSignal_

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

Protected Member Functions inherited from edm::PuttableSourceBase
void	beginJob () override

Protected Member Functions inherited from edm::InputSource
void	decreaseRemainingEventsBy (int iSkipped)

bool	eventCached () const

bool	newLumi () const

bool	newRun () const

ProcessHistoryRegistry &	processHistoryRegistryForUpdate ()

ProductRegistry &	productRegistryUpdate ()

void	reset () const

void	resetEventCached ()

void	resetLuminosityBlockAuxiliary (bool isNewLumi=true) const

void	resetNewLumi ()

void	resetNewRun ()

void	resetRunAuxiliary (bool isNewRun=true) const

void	setEventCached ()
	Called by the framework to merge or ached() const {return eventCached_;}. More...

void	setLuminosityBlockAuxiliary (LuminosityBlockAuxiliary *lbp)

void	setNewLumi ()

void	setNewRun ()

void	setRunAuxiliary (RunAuxiliary *rp)

void	setTimestamp (Timestamp const &theTime)
	To set the current time, as seen by the input source. More...

ItemType	state () const

Detailed Description

Definition at line 43 of file MCatNLOSource.h.

Constructor & Destructor Documentation

MCatNLOSource::MCatNLOSource	(	const edm::ParameterSet &	params,
		const edm::InputSourceDescription &	desc
	)

explicit

Definition at line 40 of file MCatNLOSource.cc.

References fileName, edm::FromFiles::fileNames(), and reader.

                                                             :
     ProducerSourceFromFiles(params, desc, false),
         gen::Herwig6Instance(),
     skipEvents(params.getUntrackedParameter<unsigned int>("skipEvents", 0)),
         nEvents(0),
         ihpro(0),
     processCode(params.getParameter<int>("processCode"))
 {
         std::vector<std::string> allFileNames = fileNames();
 
         // Only one filename
     if (allFileNames.size() != 1)
       throw cms::Exception("Generator|MCatNLOInterface")
         << "MCatNLOSource needs exactly one file specified. " <<std::endl;
 
     fileName = allFileNames[0];
 
     // Strip the "file:" prefix 
     if (fileName.find("file:") != 0)
       throw cms::Exception("Generator|MCatNLOInterface") << "MCatNLOSource only supports the file: scheme. "<<std::endl;
     fileName.erase(0, 5);
 
     // open input file
     reader.reset(new std::ifstream(fileName.c_str()));
 
     produces<LHEEventProduct>();
     produces<LHERunInfoProduct, edm::Transition::BeginRun>();
 }

MCatNLOSource::~MCatNLOSource ( )

virtual

Definition at line 70 of file MCatNLOSource.cc.

71 {

72 }

Member Function Documentation

void MCatNLOSource::beginRun ( edm::Run & run )

privatevirtual

Reimplemented from edm::PuttableSourceBase.

Definition at line 100 of file MCatNLOSource.cc.

References funct::abs(), LHERunInfoProduct::Header::addLine(), MCPARS_::emmax, MCPARS_::emmaxs, MCPARS_::emmin, MCPARS_::emmins, fileName, MCPARS_::gammax, MCPARS_::gammaxs, MCPARS_::gamw, MCPARS_::gamws, MCPARS_::gamz, MCPARS_::gamzs, mps_fire::i, makeConfigLine(), mcatnloupinit_(), mcpars_, eostools::move(), processCode, edm::Run::put(), lhef::CommonBlocks::readHEPRUP(), MCPARS_::rmass, MCPARS_::rmasss, runInfo, and wrapper.

 {
   InstanceWrapper wrapper(this);
 
   // call UPINIT privided by MC@NLO (v3.4)
   mcatnloupinit_(&processCode, fileName.c_str(), fileName.length());
 
   // fill HEPRUP common block and store in edm::Run
   lhef::HEPRUP heprup;
   lhef::CommonBlocks::readHEPRUP(&heprup);
 
   // make sure we write a valid LHE header, Herwig6Hadronizer
   // will interpret it correctly and set up LHAPDF
   heprup.PDFGUP.first = 0;
   heprup.PDFGUP.second = 0;
 
   std::unique_ptr<LHERunInfoProduct> runInfo(new LHERunInfoProduct(heprup));
 
   LHERunInfoProduct::Header hw6header("herwig6header");
   hw6header.addLine("\n");
   hw6header.addLine("# Herwig6 parameters\n");
   hw6header.addLine(makeConfigLine("IPROC", processCode));
   // add lines for parameter that have been touched by UPINIT
   if(mcpars_.emmins) 
     hw6header.addLine(makeConfigLine("EMMIN", mcpars_.emmin));
   if(mcpars_.emmaxs) 
     hw6header.addLine(makeConfigLine("EMMAX", mcpars_.emmax));
   if(mcpars_.gammaxs) 
     hw6header.addLine(makeConfigLine("GAMMAX",mcpars_.gammax));
   if(mcpars_.gamzs) 
     hw6header.addLine(makeConfigLine("GAMZ",mcpars_.gamz));
   if(mcpars_.gamws) 
     hw6header.addLine(makeConfigLine("GAMW",mcpars_.gamw));
   for(unsigned int i=0; i<1000; ++i) {
     if(mcpars_.rmasss[i])
       hw6header.addLine(makeConfigLine("RMASS",i+1,mcpars_.rmass[i]));
   }
 
   // other needed MC@NLO defaults (from mcatnlo_hwdriver.f v3.4)
   hw6header.addLine(makeConfigLine("SOFTME", false));
   hw6header.addLine(makeConfigLine("NOWGT", false));
   hw6header.addLine(makeConfigLine("NEGWTS", true));
   if(abs(processCode)==1705 || abs(processCode)==11705)
     hw6header.addLine(makeConfigLine("PSPLT",2,0.5));
   double wgtmax_=1.000001;
   hw6header.addLine(makeConfigLine("WGTMAX", wgtmax_));
   hw6header.addLine(makeConfigLine("AVABW", wgtmax_));
   hw6header.addLine(makeConfigLine("RLTIM",6, 1.0E-23));
   hw6header.addLine(makeConfigLine("RLTIM",12, 1.0E-23));
   
 
   runInfo->addHeader(hw6header);
 
   run.put(std::move(runInfo));
 
   return;
 }

void MCatNLOSource::endJob ( void )

privatevirtual

Reimplemented from edm::InputSource.

Definition at line 74 of file MCatNLOSource.cc.

References reader.

 {
   reader.reset();
 }

bool MCatNLOSource::hwwarn	(	const std::string &	fn,
		int	code
	)

privatevirtual

Reimplemented from gen::Herwig6Instance.

Definition at line 194 of file MCatNLOSource.cc.

References DEFINE_FWK_INPUT_SOURCE.

 {
   // dummy ignoring useless HERWIG warnings
   return true;
 }

void MCatNLOSource::nextEvent ( )

private

Definition at line 79 of file MCatNLOSource.cc.

 {
   return;
 }

void MCatNLOSource::produce ( edm::Event & event )

privatevirtual

Implements edm::ProducerSourceBase.

Definition at line 179 of file MCatNLOSource.cc.

References lhef::HEPEUP::IDPRUP, ihpro, lhef::HEPRUP::LPRUP, makeConfigLine(), eostools::move(), lhef::CommonBlocks::readHEPEUP(), lhef::CommonBlocks::readHEPRUP(), wrapper, and lhef::HEPEUP::XWGTUP.

 {
   InstanceWrapper wrapper(this);
 
   // fill HEPRUP common block and store in edm::Run
   lhef::HEPRUP heprup;
   lhef::HEPEUP hepeup;
   lhef::CommonBlocks::readHEPRUP(&heprup);
   lhef::CommonBlocks::readHEPEUP(&hepeup);
   hepeup.IDPRUP = heprup.LPRUP[0];
   std::unique_ptr<LHEEventProduct> lhEvent(new LHEEventProduct(hepeup,hepeup.XWGTUP));
   lhEvent->addComment(makeConfigLine("#IHPRO", ihpro));
   event.put(std::move(lhEvent));
 }

bool MCatNLOSource::setRunAndEventInfo	(	edm::EventID &	,
		edm::TimeValue_t &	,
		edm::EventAuxiliary::ExperimentType &
	)

privatevirtual

Implements edm::ProducerSourceBase.

Definition at line 158 of file MCatNLOSource.cc.

References ihpro, mcatnloupevnt_(), processCode, skipEvents, and wrapper.

 {
   InstanceWrapper wrapper(this);
 
   int lastEventDone=0;
   ihpro=0;
   // skip events if asked to...
 
   while(skipEvents>0) {
     skipEvents--;
     mcatnloupevnt_(&processCode,&lastEventDone,&ihpro);
     if(lastEventDone) return false;
   }
 
   // call UPINIT privided by MC@NLO (v3.4)
   mcatnloupevnt_(&processCode,&lastEventDone,&ihpro);
 
   if(lastEventDone) return false;
   return true;
 }