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
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
bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

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

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

Constructor & Destructor Documentation

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

Constructor.

Definition at line 94 of file AlpgenSource.cc.

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

                                                       :
   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();
 
   // 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>();
 }

AlpgenSource::~AlpgenSource ( )

override

Destructor.

Definition at line 147 of file AlpgenSource.cc.

148 {

149 }

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from edm::PuttableSourceBase.

Definition at line 161 of file AlpgenSource.cc.

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

 {
   // 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_.reset(new 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;
 
 }

template<typename T >

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

private

Function to get parameter by name from AlpgenHeader.

Definition at line 280 of file AlpgenSource.cc.

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

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

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

References header, and AlpgenHeader::params.

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

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

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

Referenced by beginRun(), and readAlpgenEvent().

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

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

overrideprivatevirtual

Implements edm::ProducerSourceBase.

Definition at line 433 of file AlpgenSource.cc.

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

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

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

private

Read an event and put it into the HEPEUP.

Definition at line 314 of file AlpgenSource.cc.

References edmScanValgrind::buffer, GenHFHadronMatcher_cff::flavour, mps_fire::i, lhef::HEPEUP::ICOLUP, lhef::HEPEUP::IDPRUP, lhef::HEPEUP::IDUP, inputFile_, lhef::HEPEUP::ISTUP, eostools::ls(), ResonanceBuilder::mass, lhef::HEPEUP::MOTHUP, nPart(), processID(), lhef::HEPEUP::PUP, lhef::HEPEUP::resize(), lhef::HEPEUP::SCALUP, mathSSE::sqrt(), and lhef::HEPEUP::XWGTUP.

Referenced by setRunAndEventInfo().

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

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

overrideprivatevirtual

Implements edm::ProducerSourceBase.

Definition at line 406 of file AlpgenSource.cc.

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

 {
   // The LHE Event Record
   hepeup_.reset(new 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;
 }

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

References header, ResonanceBuilder::mass, and AlpgenHeader::masses.

Referenced by beginRun().

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

Member Data Documentation

std::string AlpgenSource::extraHeaderFileName_

private

Name of the extra header file.

Definition at line 83 of file AlpgenSource.cc.

Referenced by beginRun().

std::string AlpgenSource::extraHeaderName_

private

Name given to the extra header.

Definition at line 86 of file AlpgenSource.cc.

Referenced by beginRun().

std::string AlpgenSource::fileName_

private

Name of the input file.

Definition at line 64 of file AlpgenSource.cc.

Referenced by AlpgenSource(), and beginRun().

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

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

private

Definition at line 80 of file AlpgenSource.cc.

Referenced by produce(), and setRunAndEventInfo().

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

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

unsigned long AlpgenSource::nEvent_

private

Number of events.

Definition at line 73 of file AlpgenSource.cc.

Referenced by setRunAndEventInfo().

unsigned long AlpgenSource::skipEvents_

private

Number of events to skip.

Definition at line 70 of file AlpgenSource.cc.

Referenced by setRunAndEventInfo().

bool AlpgenSource::writeAlpgenParFile

private

Definition at line 90 of file AlpgenSource.cc.

Referenced by beginRun().

bool AlpgenSource::writeAlpgenWgtFile

private

configuration flags

Definition at line 89 of file AlpgenSource.cc.

Referenced by beginRun().

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

Member Function Documentation

Member Data Documentation