Inheritance diagram for AlpgenSource:

Public Member Functions
	AlpgenSource (const edm::ParameterSet &params, const edm::InputSourceDescription &desc)
	Constructor. More...

	~AlpgenSource () override
	Destructor. More...

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

	PuttableSourceBase (const PuttableSourceBase &)=delete

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

void	fillProcessBlockHelper ()
	Fill the ProcessBlockHelper with info for the current file. 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...

bool	nextProcessBlock (ProcessBlockPrincipal &)
	Next process block, return false if there is none, sets the processName in the principal. More...

InputSource &	operator= (InputSource const &)=delete

std::shared_ptr< ProcessBlockHelper > &	processBlockHelper ()

std::shared_ptr< ProcessBlockHelper const >	processBlockHelper () const
	Accessors for processBlockHelper. More...

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::shared_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	readProcessBlock (ProcessBlockPrincipal &)
	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

Private Member Functions
void	beginRun (edm::Run &run) override

template<typename T >
T	getParameter (AlpgenHeader::Parameter index) const
	Function to get parameter by name from AlpgenHeader. More...

template<typename T >
T	getParameter (AlpgenHeader::Parameter index, const T &defValue) const
	Function to get parameter by name from AlpgenHeader, w/ default. More...

unsigned int	processID () const

void	produce (edm::Event &event) override

bool	readAlpgenEvent (lhef::HEPEUP &hepeup)
	Read an event and put it into the HEPEUP. More...

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

std::string	slhaMassLine (int pdgId, AlpgenHeader::Masses mass, const std::string &comment) const

Private Attributes
std::string	extraHeaderFileName_
	Name of the extra header file. More...

std::string	extraHeaderName_
	Name given to the extra header. More...

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

AlpgenHeader	header
	Alpgen _unw.par file as an AlpgenHeader. More...

std::unique_ptr< lhef::HEPEUP >	hepeup_

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

LHERunInfoProduct::Header	lheAlpgenUnwParHeader
	Alpgen _unw.par file as a LHE header. More...

unsigned long	nEvent_
	Number of events. More...

unsigned long	skipEvents_
	Number of events to skip. More...

bool	writeAlpgenParFile

bool	writeAlpgenWgtFile
	configuration flags More...

bool	writeExtraHeader

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

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

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

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 31 of file AlpgenSource.cc.

Constructor & Destructor Documentation

◆ AlpgenSource()

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

Constructor.

Definition at line 94 of file AlpgenSource.cc.

     : edm::ProducerSourceFromFiles(params, desc, false),
       skipEvents_(params.getUntrackedParameter<unsigned int>("skipEvents", 0)),
       nEvent_(0),
       lheAlpgenUnwParHeader("AlpgenUnwParFile"),
       extraHeaderFileName_(params.getUntrackedParameter<std::string>("extraHeaderFileName", "")),
       extraHeaderName_(params.getUntrackedParameter<std::string>("extraHeaderName", "")),
       writeAlpgenWgtFile(params.getUntrackedParameter<bool>("writeAlpgenWgtFile", true)),
       writeAlpgenParFile(params.getUntrackedParameter<bool>("writeAlpgenParFile", true)),
       writeExtraHeader(params.getUntrackedParameter<bool>("writeExtraHeader", false)) {
   std::vector<std::string> allFileNames = fileNames(0);
  
   // Only one filename
   if (allFileNames.size() != 1)
     throw cms::Exception("Generator|AlpgenInterface") << "AlpgenSource needs exactly one file specified "
                                                          "for now."
                                                       << std::endl;
  
   fileName_ = allFileNames[0];
  
   // Strip the "file:" prefix
   if (fileName_.find("file:") != 0)
     throw cms::Exception("Generator|AlpgenInterface") << "AlpgenSource only supports the file: scheme "
                                                          "for now."
                                                       << std::endl;
   fileName_.erase(0, 5);
  
   // Open the _unw.par file to store additional
   // informations in the LHERunInfoProduct
   std::ifstream reader((fileName_ + "_unw.par").c_str());
   if (!reader.good())
     throw cms::Exception("Generator|AlpgenInterface")
         << "AlpgenSource was unable to open the file \"" << fileName_ << "_unw.par\"." << std::endl;
  
   // A full copy of the _unw.par file in the LHE header.
   char buffer[256];
   while (reader.getline(buffer, sizeof buffer))
     lheAlpgenUnwParHeader.addLine(std::string(buffer) + "\n");
  
   // Parse that header to setup an Alpgen header,
   // which will be used in the production itself.
   if (!header.parse(lheAlpgenUnwParHeader.begin(), lheAlpgenUnwParHeader.end()))
     throw cms::Exception("Generator|AlpgenInterface") << "AlpgenSource was unable to parse the Alpgen "
                                                       << "unweighted parameter file." << std::endl;
  
   // Declare the products.
   produces<LHERunInfoProduct, edm::Transition::BeginRun>();
   produces<LHEEventProduct>();
 }

References LHERunInfoProduct::Header::addLine(), LHERunInfoProduct::Header::begin(), edmScanValgrind::buffer, LHERunInfoProduct::Header::end(), Exception, fileName_, edm::FromFiles::fileNames(), header, lheAlpgenUnwParHeader, AlpgenHeader::parse(), DQM::reader, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~AlpgenSource()

AlpgenSource::~AlpgenSource ( )

override

Destructor.

Definition at line 144 of file AlpgenSource.cc.

144 {}

Member Function Documentation

◆ beginRun()

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

overrideprivatevirtual

Reimplemented from edm::PuttableSourceBase.

Definition at line 153 of file AlpgenSource.cc.

                                        {
   // At this point, the lheUnwParHeader has the full contents of the _unw.par
   // file. So we can get the HEPRUP information from the LHE header itself.
   lhef::HEPRUP heprup;
  
   // Get basic run information.
   // Beam identity.
   heprup.IDBMUP.first = 2212;
   switch (getParameter<int>(AlpgenHeader::ih2)) {
     case 1:
       heprup.IDBMUP.second = 2212;
       break;
     case -1:
       heprup.IDBMUP.second = -2212;
       break;
     default:
       throw cms::Exception("Generator|AlpgenInterface") << "AlpgenSource was unable to understand the ih2 "
                                                         << "parameter." << std::endl;
   }
  
   // Beam energy.
   heprup.EBMUP.second = heprup.EBMUP.first = getParameter<double>(AlpgenHeader::ebeam);
  
   // PDF info. Initially, Alpgen doesn't fill it.
   heprup.PDFGUP.first = -1;
   heprup.PDFGUP.second = -1;
   heprup.PDFSUP.first = -1;
   heprup.PDFSUP.second = -1;
  
   // Unweighted events.
   heprup.IDWTUP = 3;
  
   // Only one process.
   heprup.resize(1);
  
   // Cross section and error.
   heprup.XSECUP[0] = header.xsec;
   heprup.XERRUP[0] = header.xsecErr;
  
   // Maximum weight.
   heprup.XMAXUP[0] = header.xsec;
  
   // Process code for Pythia.
   heprup.LPRUP[0] = processID();
  
   // Comments on top.
   LHERunInfoProduct::Header comments;
   comments.addLine("\n");
   comments.addLine("\tExtracted by AlpgenInterface\n");
  
   // Add SLHA header containing particle masses from Alpgen.
   // Pythia6Hadronisation will feed the masses to Pythia automatically.
   LHERunInfoProduct::Header slha("slha");
   slha.addLine("\n# SLHA header containing masses from Alpgen\n");
   slha.addLine("Block MASS      #  Mass spectrum (kinematic masses)\n");
   slha.addLine("#       PDG   Mass\n");
   slha.addLine(slhaMassLine(4, AlpgenHeader::mc, "charm    pole mass"));
   slha.addLine(slhaMassLine(5, AlpgenHeader::mb, "bottom   pole mass"));
   slha.addLine(slhaMassLine(6, AlpgenHeader::mt, "top      pole mass"));
   slha.addLine(slhaMassLine(23, AlpgenHeader::mz, "Z        mass"));
   slha.addLine(slhaMassLine(24, AlpgenHeader::mw, "W        mass"));
   slha.addLine(slhaMassLine(25, AlpgenHeader::mh, "H        mass"));
  
   char buffer[512];
  
   // We also add the information on weighted events.
   LHERunInfoProduct::Header lheAlpgenWgtHeader("AlpgenWgtFile");
   if (writeAlpgenWgtFile) {
     std::ifstream wgtascii((fileName_ + ".wgt").c_str());
     while (wgtascii.getline(buffer, 512)) {
       lheAlpgenWgtHeader.addLine(std::string(buffer) + "\n");
     }
   }
  
   LHERunInfoProduct::Header lheAlpgenParHeader("AlpgenParFile");
   if (writeAlpgenParFile) {
     std::ifstream parascii((fileName_ + ".par").c_str());
     while (parascii.getline(buffer, 512)) {
       lheAlpgenParHeader.addLine(std::string(buffer) + "\n");
     }
   }
  
   // If requested by the user, we also add any specific header provided.
   // Nota bene: the header is put in the LHERunInfoProduct AS IT WAS GIVEN.
   // That means NO CROSS-CHECKS WHATSOEVER. Use with care.
   LHERunInfoProduct::Header extraHeader(extraHeaderName_);
   if (writeExtraHeader) {
     std::ifstream extraascii(extraHeaderFileName_.c_str());
     while (extraascii.getline(buffer, 512)) {
       extraHeader.addLine(std::string(buffer) + "\n");
     }
   }
  
   // Build the final Run info object. Backwards-compatible order.
   std::unique_ptr<LHERunInfoProduct> runInfo(new LHERunInfoProduct(heprup));
   runInfo->addHeader(comments);
   runInfo->addHeader(lheAlpgenUnwParHeader);
   if (writeAlpgenWgtFile)
     runInfo->addHeader(lheAlpgenWgtHeader);
   if (writeAlpgenParFile)
     runInfo->addHeader(lheAlpgenParHeader);
   runInfo->addHeader(slha);
   if (writeExtraHeader)
     runInfo->addHeader(extraHeader);
   run.put(std::move(runInfo));
  
   // Open the .unw file in the heap, and set the global pointer to it.
   inputFile_ = std::make_unique<std::ifstream>((fileName_ + ".unw").c_str());
   if (!inputFile_->good())
     throw cms::Exception("Generator|AlpgenInterface")
         << "AlpgenSource was unable to open the file \"" << fileName_ << ".unw\"." << std::endl;
 }

References LHERunInfoProduct::Header::addLine(), edmScanValgrind::buffer, AlpgenHeader::ebeam, lhef::HEPRUP::EBMUP, Exception, extraHeaderFileName_, extraHeaderName_, fileName_, header, lhef::HEPRUP::IDBMUP, lhef::HEPRUP::IDWTUP, AlpgenHeader::ih2, inputFile_, lheAlpgenUnwParHeader, lhe2HepMCConverter_cff::LHERunInfoProduct, lhef::HEPRUP::LPRUP, AlpgenHeader::mb, AlpgenHeader::mc, AlpgenHeader::mh, eostools::move(), AlpgenHeader::mt, AlpgenHeader::mw, AlpgenHeader::mz, lhef::HEPRUP::PDFGUP, lhef::HEPRUP::PDFSUP, processID(), lhef::HEPRUP::resize(), edm::IDGeneratorSourceBase< PuttableSourceBase >::run(), submitPVValidationJobs::runInfo, slhaMassLine(), AlCaHLTBitMon_QueryRunRegistry::string, writeAlpgenParFile, writeAlpgenWgtFile, writeExtraHeader, lhef::HEPRUP::XERRUP, lhef::HEPRUP::XMAXUP, AlpgenHeader::xsec, AlpgenHeader::xsecErr, and lhef::HEPRUP::XSECUP.

◆ getParameter() [1/2]

template<typename T >

T AlpgenSource::getParameter ( AlpgenHeader::Parameter index ) const

private

Function to get parameter by name from AlpgenHeader.

Definition at line 267 of file AlpgenSource.cc.

                                                               {
   std::map<AlpgenHeader::Parameter, double>::const_iterator pos = header.params.find(index);
   if (pos == header.params.end())
     throw cms::Exception("Generator|AlpgenInterface")
         << "Requested Alpgen parameter \"" << AlpgenHeader::parameterName(index)
         << "\" "
            "not found in Alpgen parameter file."
         << std::endl;
  
   return T(pos->second);
 }

References header, AlpgenHeader::parameterName(), and AlpgenHeader::params.

◆ getParameter() [2/2]

template<typename T >

T AlpgenSource::getParameter	(	AlpgenHeader::Parameter	index,
		const T &	defValue
	)		const

private

Function to get parameter by name from AlpgenHeader, w/ default.

Definition at line 280 of file AlpgenSource.cc.

                                                                                  {
   std::map<AlpgenHeader::Parameter, double>::const_iterator pos = header.params.find(index);
   if (pos == header.params.end())
     return defValue;
   else
     return T(pos->second);
 }

References header, and AlpgenHeader::params.

◆ processID()

unsigned int AlpgenSource::processID ( ) const

private

The Alpgen process ID. This is defined as processID() = 100*X + 10*Y + Z, where = ihrd, Y = ihvy, Z = njets

Definition at line 288 of file AlpgenSource.cc.

                                            {
   // return 661; // The original, old thing.
   // digits #XYZ: X = ihrd, Y = ihvy, Z = njets
   return header.ihrd * 100 + getParameter<unsigned int>(AlpgenHeader::ihvy, 0) * 10 +
          getParameter<unsigned int>(AlpgenHeader::njets, 0);
 }

References header, AlpgenHeader::ihrd, AlpgenHeader::ihvy, and AlpgenHeader::njets.

Referenced by beginRun(), and readAlpgenEvent().

◆ produce()

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

overrideprivatevirtual

Implements edm::ProducerSourceBase.

Definition at line 422 of file AlpgenSource.cc.

                                           {
   // Here are the Alpgen routines for filling up the rest
   // of the LHE Event Record. The .unw file has the information
   // in a compressed way, e.g. it doesn't list the W boson -
   // one has to reconstruct it from the e nu pair.
   lhef::HEPEUP &hepeup = *hepeup_;
  
   switch (header.ihrd) {
     case 1:
     case 2:
     case 3:
     case 4:
     case 10:
     case 14:
     case 15:
       alpgen::fixEventWZ(hepeup);
       break;
     case 5:
       alpgen::fixEventMultiBoson(hepeup);
       break;
     case 6:
       alpgen::fixEventTTbar(hepeup);
       break;
     case 8:
       alpgen::fixEventHiggsTTbar(hepeup);
       break;
     case 13:
       alpgen::fixEventSingleTop(hepeup, header.masses[AlpgenHeader::mb], int(header.params[AlpgenHeader::itopprc]));
       break;
     case 7:
     case 9:
     case 11:
     case 12:
     case 16:
       // No fixes needed.
       break;
  
     default:
       throw cms::Exception("Generator|AlpgenInterface") << "Unrecognized IHRD process code" << std::endl;
   }
  
   // Create the LHEEventProduct and put it into the Event.
   std::unique_ptr<LHEEventProduct> lheEvent(new LHEEventProduct(hepeup));
   event.put(std::move(lheEvent));
  
   hepeup_.reset();
 }

References Exception, alpgen::fixEventHiggsTTbar(), alpgen::fixEventMultiBoson(), alpgen::fixEventSingleTop(), alpgen::fixEventTTbar(), alpgen::fixEventWZ(), header, hepeup_, AlpgenHeader::ihrd, AlpgenHeader::itopprc, lhe2HepMCConverter_cff::LHEEventProduct, AlpgenHeader::masses, AlpgenHeader::mb, eostools::move(), and AlpgenHeader::params.

◆ readAlpgenEvent()

bool AlpgenSource::readAlpgenEvent ( lhef::HEPEUP & hepeup )

private

Read an event and put it into the HEPEUP.

Definition at line 295 of file AlpgenSource.cc.

                                                      {
   char buffer[512];
   double dummy;
   int nPart;
   double sWgtRes;
   double sQ;
  
   inputFile_->getline(buffer, sizeof buffer);
   if (!inputFile_->good())
     return false;
  
   std::istringstream ls(buffer);
  
   // Event number and process don't matter (or do they?)
   ls >> dummy >> dummy;
  
   // Number of particles in the record, sample's average weight and Q scale
   ls >> nPart >> sWgtRes >> sQ;
  
   if (ls.bad() || nPart < 1 || nPart > 1000)
     return false;
  
   // Make room for the particles listed in the Alpgen file
   hepeup.resize(nPart);
  
   // Scales, weights and process ID.
   hepeup.SCALUP = sQ;
   hepeup.XWGTUP = sWgtRes;
   hepeup.IDPRUP = processID();
  
   // Incoming lines
   for (int i = 0; i != 2; ++i) {
     inputFile_->getline(buffer, sizeof buffer);
     std::istringstream ls(buffer);
     int flavour;
     ls >> flavour;
     int colour1;
     ls >> colour1;
     int colour2;
     ls >> colour2;
     double zMomentum;
     ls >> zMomentum;
  
     if (inputFile_->bad())
       return false;
  
     // Setting the HEPEUP of the incoming lines.
     hepeup.IDUP[i] = flavour;
     hepeup.ISTUP[i] = -1;
     hepeup.MOTHUP[i].first = 0;
     hepeup.MOTHUP[i].second = 0;
     hepeup.PUP[i][0] = 0.;
     hepeup.PUP[i][1] = 0.;
     hepeup.PUP[i][2] = zMomentum;
     hepeup.PUP[i][3] = std::fabs(zMomentum);
     hepeup.PUP[i][4] = 0.;
     if (colour1)
       colour1 += 500;
     if (colour2)
       colour2 += 500;
     hepeup.ICOLUP[i].first = colour1;
     hepeup.ICOLUP[i].second = colour2;
   }
  
   // Outgoing lines
   for (int i = 2; i != nPart; ++i) {
     inputFile_->getline(buffer, sizeof buffer);
     std::istringstream ls(buffer);
     int flavour;
     ls >> flavour;
     int colour1;
     ls >> colour1;
     int colour2;
     ls >> colour2;
     double px, py, pz, mass;
     ls >> px >> py >> pz >> mass;
     double energy = std::sqrt(px * px + py * py + pz * pz + mass * mass);
  
     if (inputFile_->bad())
       return false;
  
     // Setting the HEPEUP of the outgoing lines.
     hepeup.IDUP[i] = flavour;
     hepeup.ISTUP[i] = 1;
     hepeup.MOTHUP[i].first = 1;
     hepeup.MOTHUP[i].second = 2;
     hepeup.PUP[i][0] = px;
     hepeup.PUP[i][1] = py;
     hepeup.PUP[i][2] = pz;
     hepeup.PUP[i][3] = energy;
     hepeup.PUP[i][4] = mass;
     if (colour1)
       colour1 += 500;
     if (colour2)
       colour2 += 500;
     hepeup.ICOLUP[i].first = colour1;
     hepeup.ICOLUP[i].second = colour2;
   }
  
   return true;
 }

References edmScanValgrind::buffer, HCALHighEnergyHPDFilter_cfi::energy, mps_fire::i, lhef::HEPEUP::ICOLUP, lhef::HEPEUP::IDPRUP, lhef::HEPEUP::IDUP, inputFile_, lhef::HEPEUP::ISTUP, eostools::ls(), EgHLTOffHistBins_cfi::mass, lhef::HEPEUP::MOTHUP, nPart(), processID(), lhef::HEPEUP::PUP, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, lhef::HEPEUP::resize(), lhef::HEPEUP::SCALUP, mathSSE::sqrt(), and lhef::HEPEUP::XWGTUP.

Referenced by setRunAndEventInfo().

◆ setRunAndEventInfo()

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

overrideprivatevirtual

Implements edm::IDGeneratorSourceBase< PuttableSourceBase >.

Definition at line 397 of file AlpgenSource.cc.

                                                                                                        {
   // The LHE Event Record
   hepeup_ = std::make_unique<lhef::HEPEUP>();
  
   lhef::HEPEUP &hepeup = *hepeup_;
  
   // Read the .unw file until it is over.
   for (;;) {
     if (!readAlpgenEvent(hepeup)) {
       if (inputFile_->eof())
         return false;
  
       throw cms::Exception("Generator|AlpgenInterface")
           << "AlpgenSource is not able read event no. " << nEvent_ << std::endl;
     }
  
     nEvent_++;
     if (skipEvents_ > 0)
       skipEvents_--;
     else
       break;
   }
   return true;
 }

References Exception, hepeup_, inputFile_, nEvent_, readAlpgenEvent(), and skipEvents_.

◆ slhaMassLine()

std::string AlpgenSource::slhaMassLine	(	int	pdgId,
		AlpgenHeader::Masses	mass,
		const std::string &	comment
	)		const

private

Converts the AlpgenHeader::Masses to a std::string formatted as a slhaMassLine to facilitate passing them to Alpgen.

Definition at line 146 of file AlpgenSource.cc.

                                                                                                        {
   std::ostringstream ss;
   ss << std::setw(9) << pdgId << "     " << std::scientific << std::setprecision(9) << header.masses[mass] << "   # "
      << comment << std::endl;
   return ss.str();
 }

References header, EgHLTOffHistBins_cfi::mass, AlpgenHeader::masses, EgammaValidation_cff::pdgId, and contentValuesCheck::ss.

Referenced by beginRun().

Member Data Documentation

◆ extraHeaderFileName_

std::string AlpgenSource::extraHeaderFileName_

private

Name of the extra header file.

Definition at line 83 of file AlpgenSource.cc.

Referenced by beginRun().

◆ extraHeaderName_

std::string AlpgenSource::extraHeaderName_

private

Name given to the extra header.

Definition at line 86 of file AlpgenSource.cc.

Referenced by beginRun().

◆ fileName_

std::string AlpgenSource::fileName_

private

Name of the input file.

Definition at line 64 of file AlpgenSource.cc.

Referenced by AlpgenSource(), and beginRun().

◆ header

AlpgenHeader AlpgenSource::header

private

Alpgen _unw.par file as an AlpgenHeader.

Definition at line 78 of file AlpgenSource.cc.

Referenced by AlpgenSource(), beginRun(), getParameter(), processID(), produce(), and slhaMassLine().

◆ hepeup_

std::unique_ptr<lhef::HEPEUP> AlpgenSource::hepeup_

private

Definition at line 80 of file AlpgenSource.cc.

Referenced by produce(), and setRunAndEventInfo().

◆ inputFile_

std::unique_ptr<std::ifstream> AlpgenSource::inputFile_

private

Pointer to the input file.

Definition at line 67 of file AlpgenSource.cc.

Referenced by beginRun(), readAlpgenEvent(), and setRunAndEventInfo().

◆ lheAlpgenUnwParHeader

LHERunInfoProduct::Header AlpgenSource::lheAlpgenUnwParHeader

private

Alpgen _unw.par file as a LHE header.

Definition at line 76 of file AlpgenSource.cc.

Referenced by AlpgenSource(), and beginRun().

◆ nEvent_

unsigned long AlpgenSource::nEvent_

private

Number of events.

Definition at line 73 of file AlpgenSource.cc.

Referenced by setRunAndEventInfo().

◆ skipEvents_

unsigned long AlpgenSource::skipEvents_

private

Number of events to skip.

Definition at line 70 of file AlpgenSource.cc.

Referenced by setRunAndEventInfo().

◆ writeAlpgenParFile

bool AlpgenSource::writeAlpgenParFile

private

Definition at line 90 of file AlpgenSource.cc.

Referenced by beginRun().

◆ writeAlpgenWgtFile

bool AlpgenSource::writeAlpgenWgtFile

private

configuration flags

Definition at line 89 of file AlpgenSource.cc.

Referenced by beginRun().

◆ writeExtraHeader

bool AlpgenSource::writeExtraHeader

private

Definition at line 91 of file AlpgenSource.cc.

Referenced by beginRun().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AlpgenSource()

◆ ~AlpgenSource()

Member Function Documentation

◆ beginRun()

◆ getParameter() [1/2]

◆ getParameter() [2/2]

◆ processID()

◆ produce()

◆ readAlpgenEvent()

◆ setRunAndEventInfo()

◆ slhaMassLine()

Member Data Documentation

◆ extraHeaderFileName_

◆ extraHeaderName_

◆ fileName_

◆ header

◆ hepeup_

◆ inputFile_

◆ lheAlpgenUnwParHeader

◆ nEvent_

◆ skipEvents_

◆ writeAlpgenParFile

◆ writeAlpgenWgtFile

◆ writeExtraHeader