#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
virtual bool	noFiles () const override

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

virtual	~ProducerSourceFromFiles ()

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

virtual	~ProducerSourceBase ()

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

boost::shared_ptr < BranchIDListHelper >	branchIDListHelper () const
	Accessor for branchIDListHelper. More...

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

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

void	doBeginLumi (LuminosityBlockPrincipal &lbp, ProcessContext const *)
	Called by framework at beginning of lumi block. More...

void	doBeginRun (RunPrincipal &rp, ProcessContext const *)
	Called by framework at beginning of run. More...

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

void	doEndLumi (LuminosityBlockPrincipal &lbp, bool cleaningUpAfterException, ProcessContext const *)
	Called by framework at end of lumi block. More...

void	doEndRun (RunPrincipal &rp, bool cleaningUpAfterException, ProcessContext const *)
	Called by framework at end of run. More...

void	doPostForkReacquireResources (boost::shared_ptr< multicore::MessageReceiverForSource >)

void	doPreForkReleaseResources ()
	Called by the framework before forking the process. 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...

boost::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

bool	primary () const
	Accessor for primary input source flag. More...

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

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

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

ProcessHistoryRegistry &	processHistoryRegistryForUpdate ()
	Non-const accessor for process history registry. More...

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

boost::shared_ptr < ProductRegistry const >	productRegistry () const
	Accessor 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...

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

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

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

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

ProcessHistoryID const &	reducedProcessHistoryID () const

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

int	remainingEvents () const

int	remainingLuminosityBlocks () const

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

SharedResourcesAcquirer *	resourceSharedWithDelayedReader () const
	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...

boost::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)

bool	skipForForking ()

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

virtual	~InputSource ()
	Destructor. More...

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 ()

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 &)

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 }

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 ()

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

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

bool	eventCached () const

boost::shared_ptr < LuminosityBlockPrincipal > const	luminosityBlockPrincipal () const

bool	newLumi () const

bool	newRun () const

ProcessHistoryRegistry &	processHistoryRegistryUpdate () const

ProductRegistry &	productRegistryUpdate () const

void	reset () const

void	resetEventCached ()

void	resetLuminosityBlockAuxiliary (bool isNewLumi=true) const

void	resetNewLumi ()

void	resetNewRun ()

void	resetRunAuxiliary (bool isNewRun=true) const

boost::shared_ptr < RunPrincipal > const	runPrincipal () 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::InRun>();
 }

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

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, i, makeConfigLine(), mcatnloupinit_(), mcpars_, 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::auto_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(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.

 {
   // 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(), 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::auto_ptr<LHEEventProduct> lhEvent(new LHEEventProduct(hepeup,hepeup.XWGTUP));
   lhEvent->addComment(makeConfigLine("#IHPRO", ihpro));
   event.put(lhEvent);
 }

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

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

Member Data Documentation

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

private

Definition at line 79 of file MCatNLOSource.h.

Referenced by Types.EventID::cppID().

std::string MCatNLOSource::fileName

private

Name of the input file.

Definition at line 61 of file MCatNLOSource.h.

Referenced by beginRun(), and MCatNLOSource().

int MCatNLOSource::ihpro

private

Definition at line 72 of file MCatNLOSource.h.

Referenced by produce(), and setRunAndEventInfo().

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

private

Pointer to the input file.

Definition at line 64 of file MCatNLOSource.h.

unsigned int MCatNLOSource::nEvents

private

Number of events.

Definition at line 70 of file MCatNLOSource.h.

int MCatNLOSource::processCode

private

Definition at line 74 of file MCatNLOSource.h.

Referenced by beginRun(), and setRunAndEventInfo().

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

private

Definition at line 76 of file MCatNLOSource.h.

Referenced by endJob(), and MCatNLOSource().

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

private

Definition at line 78 of file MCatNLOSource.h.

Referenced by beginRun().

unsigned int MCatNLOSource::skipEvents

private

Number of events to skip.

Definition at line 67 of file MCatNLOSource.h.

Referenced by setRunAndEventInfo().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation