#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::ExternalInputSource
InputFileCatalog &	catalog ()

	ExternalInputSource (ParameterSet const &pset, InputSourceDescription const &desc, bool realData=true)

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

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

virtual	~ExternalInputSource ()

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

EventNumber_t	event () const

unsigned int	eventCreationDelay () 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

RunNumber_t	run () const

unsigned int	timeBetweenEvents () const

virtual	~ConfigurableInputSource ()

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

void	closeFile (boost::shared_ptr< FileBlock >, bool cleaningUpAfterException)
	close current file More...

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

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

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

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

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

void	doEndRun (RunPrincipal &rp, bool cleaningUpAfterException)
	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...

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	markLumi ()
	Mark lumi as read. More...

int	markRun ()
	Mark run as read. More...

int	maxEvents () const

int	maxLuminosityBlocks () const

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

ItemType	nextItemType ()

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

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	readAndCacheLumi (bool merge, HistoryAppender &historyAppender)
	Read next luminosity block. More...

void	readAndCacheRun (bool merge, HistoryAppender &historyAppender)
	Read next run. More...

EventPrincipal *	readEvent (boost::shared_ptr< LuminosityBlockPrincipal > lbCache)

EventPrincipal *	readEvent (EventID const &)
	Read a specific event. More...

boost::shared_ptr< FileBlock >	readFile ()
	Read next file. More...

boost::shared_ptr < LuminosityBlockAuxiliary >	readLuminosityBlockAuxiliary ()
	Read next luminosity block Auxilary. 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...

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)

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 (CLHEP::HepRandomEngine *randomEngine=0)

	Herwig6Instance (int dummy)

void	openParticleSpecFile (const std::string fileName)

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

Protected Member Functions
virtual void	beginRun (edm::Run &run)

virtual void	endJob ()

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

void	nextEvent ()

virtual bool	produce (edm::Event &event)

Protected Member Functions inherited from edm::ExternalInputSource
void	incrementFileIndex ()

Protected Member Functions inherited from edm::ConfigurableInputSource
void	reallyReadEvent ()

void	setEventNumber (EventNumber_t e)

void	setTime (TimeValue_t t)

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

EventPrincipal *	eventPrincipalCache ()

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

PrincipalCache const &	principalCache () const

PrincipalCache &	principalCache ()

ProductRegistry &	productRegistryUpdate () const

void	reset () const

void	resetLuminosityBlockAuxiliary () const

void	resetRunAuxiliary () const

boost::shared_ptr < RunPrincipal > const	runPrincipal () const

void	setLuminosityBlockAuxiliary (LuminosityBlockAuxiliary *lbp)

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 Attributes
boost::shared_ptr< lhef::LHEEvent >	event

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

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

unsigned int	nEvents
	Number of events. More...

int	processCode

std::auto_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 }

enum	ProcessingMode { Runs, RunsAndLumis, RunsLumisAndEvents }

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

Static Public Member Functions inherited from edm::ConfigurableInputSource
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 ()

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 41 of file MCatNLOSource.cc.

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

                                                             :
     ExternalInputSource(params, desc, false),
         gen::Herwig6Instance(0),
     skipEvents(params.getUntrackedParameter<unsigned int>("skipEvents", 0)),
         nEvents(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 )

protectedvirtual

Reimplemented from edm::ConfigurableInputSource.

Definition at line 100 of file MCatNLOSource.cc.

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

protectedvirtual

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
	)

protectedvirtual

Reimplemented from gen::Herwig6Instance.

Definition at line 190 of file MCatNLOSource.cc.

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

void MCatNLOSource::nextEvent ( )

protectedvirtual

Reimplemented from edm::ConfigurableInputSource.

Definition at line 79 of file MCatNLOSource.cc.

 {
   return;
 }

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

protectedvirtual

Implements edm::ConfigurableInputSource.

Definition at line 158 of file MCatNLOSource.cc.

References lhef::HEPEUP::IDPRUP, lhef::HEPRUP::LPRUP, makeConfigLine(), mcatnloupevnt_(), processCode, lhef::CommonBlocks::readHEPEUP(), lhef::CommonBlocks::readHEPRUP(), skipEvents, wrapper, and lhef::HEPEUP::XWGTUP.

 {
   InstanceWrapper wrapper(this);
 
   int lastEventDone=0;
   int 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;
 
   // 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);
 
   return true;
 }

Member Data Documentation

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

protected

Definition at line 76 of file MCatNLOSource.h.

std::string MCatNLOSource::fileName

protected

Name of the input file.

Definition at line 60 of file MCatNLOSource.h.

Referenced by beginRun(), and MCatNLOSource().

std::auto_ptr<std::ifstream> MCatNLOSource::inputFile

protected

Pointer to the input file.

Definition at line 63 of file MCatNLOSource.h.

unsigned int MCatNLOSource::nEvents

protected

Number of events.

Definition at line 69 of file MCatNLOSource.h.

int MCatNLOSource::processCode

protected

Definition at line 71 of file MCatNLOSource.h.

Referenced by beginRun(), and produce().

std::auto_ptr<std::ifstream> MCatNLOSource::reader

protected

Definition at line 73 of file MCatNLOSource.h.

Referenced by endJob(), and MCatNLOSource().

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

protected

Definition at line 75 of file MCatNLOSource.h.

Referenced by beginRun().

unsigned int MCatNLOSource::skipEvents

protected

Number of events to skip.

Definition at line 66 of file MCatNLOSource.h.

Referenced by produce().

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation