#include <MCatNLOSource.h>

Inheritance diagram for MCatNLOSource:

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

	~MCatNLOSource () override

Public Member Functions inherited from edm::ProducerSourceFromFiles
InputFileCatalog &	catalog ()

std::vector< std::string >	fileNames (unsigned iCatalog) const

std::vector< std::string > const &	logicalFileNames () const

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

	~ProducerSourceFromFiles () override

Public Member Functions inherited from edm::ProducerSourceBase
	ProducerSourceBase (ParameterSet const &pset, InputSourceDescription const &desc, bool realData)

	~ProducerSourceBase () noexcept(false) override

Public Member Functions inherited from edm::IDGeneratorSourceBase< PuttableSourceBase >
EventNumber_t	event () const

unsigned int	eventCreationDelay () const

EventID const &	eventID () const

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

LuminosityBlockNumber_t	luminosityBlock () const

unsigned int	numberEventsInLumi () const

unsigned int	numberEventsInRun () const

unsigned int	numberEventsInThisLumi () const

unsigned int	numberEventsInThisRun () const

TimeValue_t	presentTime () const

RunNumber_t	run () const

unsigned int	timeBetweenEvents () const

	~IDGeneratorSourceBase () noexcept(false) override

Public Member Functions inherited from edm::PuttableSourceBase
bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

	PuttableSourceBase (ParameterSet const &, InputSourceDescription const &)

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

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

std::shared_ptr< BranchIDListHelper > &	branchIDListHelper ()

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

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 (InputSource const &)=delete

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

void	issueReports (EventID const &eventID, StreamID streamID)
	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 &	processHistoryRegistry ()

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

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

std::shared_ptr< ProductRegistry > &	productRegistry ()

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

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

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

void	readEvent (EventPrincipal &ep, StreamContext &)
	Read next 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...

bool	readProcessBlock ()
	Read next process block. 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 > &	thinnedAssociationsHelper ()

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

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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

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
template<typename T >
T	call (T(&fn)())

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

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

void	call (void(&fn)())

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

template<typename A1 , typename A2 >
void	call (void(&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	beginRun (edm::Run &run) override

void	endJob () override

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

void	nextEvent ()

void	produce (edm::Event &event) override

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

Private Attributes
std::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

std::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::IDGeneratorSourceBase< PuttableSourceBase >
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::ProducerSourceFromFiles
void	incrementFileIndex ()

Protected Member Functions inherited from edm::IDGeneratorSourceBase< PuttableSourceBase >
void	doReadEvent (EventPrincipal &eventPrincipal, F &&f)

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

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Detailed Description

Definition at line 40 of file MCatNLOSource.h.

Constructor & Destructor Documentation

◆ MCatNLOSource()

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

explicit

Definition at line 40 of file MCatNLOSource.cc.

     : 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(0);
  
   // 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 = std::make_unique<std::ifstream>(fileName.c_str());
  
   produces<LHEEventProduct>();
   produces<LHERunInfoProduct, edm::Transition::BeginRun>();
 }

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

◆ ~MCatNLOSource()

MCatNLOSource::~MCatNLOSource ( )

override

Definition at line 68 of file MCatNLOSource.cc.

68 {}

Member Function Documentation

◆ beginRun()

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

overrideprivatevirtual

Reimplemented from edm::PuttableSourceBase.

Definition at line 88 of file MCatNLOSource.cc.

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

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, lhe2HepMCConverter_cff::LHERunInfoProduct, makeConfigLine(), mcatnloupinit_(), mcpars_, eostools::move(), processCode, lhef::CommonBlocks::readHEPRUP(), MCPARS_::rmass, MCPARS_::rmasss, edm::IDGeneratorSourceBase< PuttableSourceBase >::run(), runInfo, and wrapper.

◆ endJob()

void MCatNLOSource::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::InputSource.

Definition at line 70 of file MCatNLOSource.cc.

70 { reader.reset(); }

References reader.

◆ hwwarn()

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

overrideprivatevirtual

Reimplemented from gen::Herwig6Instance.

Definition at line 180 of file MCatNLOSource.cc.

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

◆ nextEvent()

void MCatNLOSource::nextEvent ( )

private

Definition at line 72 of file MCatNLOSource.cc.

72 { return; }

◆ produce()

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

overrideprivatevirtual

Implements edm::ProducerSourceBase.

Definition at line 166 of file MCatNLOSource.cc.

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

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

◆ setRunAndEventInfo()

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

overrideprivatevirtual

Implements edm::IDGeneratorSourceBase< PuttableSourceBase >.

Definition at line 144 of file MCatNLOSource.cc.

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

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

Member Data Documentation

◆ event

std::shared_ptr<lhef::LHEEvent> MCatNLOSource::event

private

Definition at line 74 of file MCatNLOSource.h.

Referenced by Types.EventID::cppID(), looper.Looper::process(), and core.AutoHandle.AutoHandle::product().

◆ fileName

std::string MCatNLOSource::fileName

private

Name of the input file.

Definition at line 56 of file MCatNLOSource.h.

Referenced by beginRun(), MCatNLOSource(), and readConfig.fileINI::read().

◆ ihpro

int MCatNLOSource::ihpro

private

Definition at line 67 of file MCatNLOSource.h.

Referenced by produce(), and setRunAndEventInfo().

◆ inputFile

std::unique_ptr<std::ifstream> MCatNLOSource::inputFile

private

Pointer to the input file.

Definition at line 59 of file MCatNLOSource.h.

◆ nEvents

unsigned int MCatNLOSource::nEvents

private

Number of events.

Definition at line 65 of file MCatNLOSource.h.

Referenced by looper.Looper::loop().

◆ processCode

int MCatNLOSource::processCode

private

Definition at line 69 of file MCatNLOSource.h.

Referenced by beginRun(), and setRunAndEventInfo().

◆ reader

std::unique_ptr<std::ifstream> MCatNLOSource::reader

private

Definition at line 71 of file MCatNLOSource.h.

Referenced by endJob(), and MCatNLOSource().

◆ runInfo

std::shared_ptr<lhef::LHERunInfo> MCatNLOSource::runInfo

private

Definition at line 73 of file MCatNLOSource.h.

Referenced by beginRun().

◆ skipEvents

unsigned int MCatNLOSource::skipEvents

private

Number of events to skip.

Definition at line 62 of file MCatNLOSource.h.

Referenced by setRunAndEventInfo().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MCatNLOSource()

◆ ~MCatNLOSource()

Member Function Documentation

◆ beginRun()

◆ endJob()

◆ hwwarn()

◆ nextEvent()

◆ produce()

◆ setRunAndEventInfo()

Member Data Documentation

◆ event

◆ fileName

◆ ihpro

◆ inputFile

◆ nEvents

◆ processCode

◆ reader

◆ runInfo

◆ skipEvents