#include <FlatEGunASCIIWriter.h>

Inheritance diagram for edm::FlatEGunASCIIWriter:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

void	beginRun (const edm::Run &, const EventSetup &) override

	FlatEGunASCIIWriter (const edm::ParameterSet &)

	~FlatEGunASCIIWriter () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
int	fCurrentEvent

HepMC::GenEvent *	fEvt

double	fMaxE

double	fMaxEta

double	fMaxPhi

double	fMinE

double	fMinEta

double	fMinPhi

std::string	fOutFileName

HepMC::IO_GenEvent *	fOutStream

std::vector< int >	fPartIDs

ESHandle< HepPDT::ParticleDataTable >	fPDGTable

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 33 of file FlatEGunASCIIWriter.h.

Constructor & Destructor Documentation

◆ FlatEGunASCIIWriter()

FlatEGunASCIIWriter::FlatEGunASCIIWriter ( const edm::ParameterSet & pset )

explicit

Definition at line 37 of file FlatEGunASCIIWriter.cc.

     : fEvt(nullptr),
       fOutFileName(pset.getUntrackedParameter<string>("OutFileName", "FlatEGunHepMC.dat")),
       fCurrentEvent(0) {
   ParameterSet pgun_params = pset.getParameter<ParameterSet>("PGunParameters");
   fPartIDs = pgun_params.getParameter<vector<int> >("PartID");
   fMinEta = pgun_params.getParameter<double>("MinEta");
   fMaxEta = pgun_params.getParameter<double>("MaxEta");
   fMinPhi = pgun_params.getParameter<double>("MinPhi");
   fMaxPhi = pgun_params.getParameter<double>("MaxPhi");
   fMinE = pgun_params.getParameter<double>("MinE");
   fMaxE = pgun_params.getParameter<double>("MaxE");
 }

References fMaxE, fMaxEta, fMaxPhi, fMinE, fMinEta, fMinPhi, fPartIDs, edm::ParameterSet::getParameter(), and muonDTDigis_cfi::pset.

◆ ~FlatEGunASCIIWriter()

FlatEGunASCIIWriter::~FlatEGunASCIIWriter ( )

override

Definition at line 51 of file FlatEGunASCIIWriter.cc.

                                           {
   if (fEvt != nullptr)
     delete fEvt;
   if (fOutStream)
     delete fOutStream;
 }

References fEvt, and fOutStream.

Member Function Documentation

◆ analyze()

void FlatEGunASCIIWriter::analyze	(	const edm::Event &	,
		const edm::EventSetup &
	)

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 70 of file FlatEGunASCIIWriter.cc.

                                                                    {
   // clean up GenEvent memory : also deletes all vtx/part in it
   //
   if (fEvt != nullptr)
     delete fEvt;
  
   // here re-create fEvt (memory)
   fEvt = new HepMC::GenEvent();
  
   // now actualy, cook up the event from PDGTable and gun parameters
   //
   //HepMC::GenVertex* Vtx = new HepMC::GenVertex( CLHEP::HepLorentzVector(0.,0.,0.) );
   HepMC::GenVertex* Vtx = new HepMC::GenVertex(HepMC::FourVector(0., 0., 0.));
  
   // loop over particles
   //
   for (unsigned int ip = 0; ip < fPartIDs.size(); ip++) {
     double energy = CLHEP::RandFlat::shoot(fMinE, fMaxE);
     double eta = CLHEP::RandFlat::shoot(fMinEta, fMaxEta);
     double phi = CLHEP::RandFlat::shoot(fMinPhi, fMaxPhi);
     const HepPDT::ParticleData* PData = fPDGTable->particle(HepPDT::ParticleID(abs(fPartIDs[ip])));
     double mass = PData->mass().value();
     double mom2 = energy * energy - mass * mass;
     double mom = 0.;
     if (mom2 > 0.)
       mom = sqrt(mom2);
     double theta = 2. * atan(exp(-eta));
     double px = mom * sin(theta) * cos(phi);
     double py = mom * sin(theta) * sin(phi);
     double pz = mom * cos(theta);
     //CLHEP::Hep3Vector p(px,py,pz) ;
     //HepMC::GenParticle* Part =
     //    new HepMC::GenParticle(CLHEP::HepLorentzVector(p,energy),fPartIDs[ip],1);
     HepMC::FourVector p(px, py, pz, energy);
     HepMC::GenParticle* Part = new HepMC::GenParticle(p, fPartIDs[ip], 1);
     Vtx->add_particle_out(Part);
   }
   fEvt->add_vertex(Vtx);
   fEvt->set_event_number(fCurrentEvent + 1);
   fEvt->set_signal_process_id(20);
  
   // for testing purpose only
   // fEvt->print() ;
  
   fOutStream->write_event(fEvt);
  
   fCurrentEvent++;
  
   return;
 }

References funct::abs(), funct::cos(), HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, JetChargeProducer_cfi::exp, fCurrentEvent, fEvt, fMaxE, fMaxEta, fMaxPhi, fMinE, fMinEta, fMinPhi, fOutStream, fPartIDs, fPDGTable, GenParticle::GenParticle, EgHLTOffHistBins_cfi::mass, AlCaHLTBitMon_ParallelJobs::p, LHEGenericFilter_cfi::ParticleID, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, funct::sin(), mathSSE::sqrt(), and theta().

◆ beginJob()

void FlatEGunASCIIWriter::beginJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 60 of file FlatEGunASCIIWriter.cc.

                                    {
   fOutStream = new HepMC::IO_GenEvent(fOutFileName);
   if (fOutStream->rdstate() == std::ios::failbit) {
     throw cms::Exception("FileNotOpen", "FlatEGunASCIIWriter::beginJob()")
         << "File " << fOutFileName << " was not open.\n";
   }
  
   return;
 }

References Exception, fOutFileName, and fOutStream.

◆ beginRun()

void FlatEGunASCIIWriter::beginRun	(	const edm::Run &	r,
		const EventSetup &	es
	)

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 58 of file FlatEGunASCIIWriter.cc.

58 { es.getData(fPDGTable); }

References fPDGTable, and edm::EventSetup::getData().

Member Data Documentation

◆ fCurrentEvent

int edm::FlatEGunASCIIWriter::fCurrentEvent

private

Definition at line 65 of file FlatEGunASCIIWriter.h.

Referenced by analyze().

◆ fEvt

HepMC::GenEvent* edm::FlatEGunASCIIWriter::fEvt

private

Definition at line 57 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and ~FlatEGunASCIIWriter().

◆ fMaxE

double edm::FlatEGunASCIIWriter::fMaxE

private

Definition at line 54 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and FlatEGunASCIIWriter().

◆ fMaxEta

double edm::FlatEGunASCIIWriter::fMaxEta

private

Definition at line 50 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and FlatEGunASCIIWriter().

◆ fMaxPhi

double edm::FlatEGunASCIIWriter::fMaxPhi

private

Definition at line 52 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and FlatEGunASCIIWriter().

◆ fMinE

double edm::FlatEGunASCIIWriter::fMinE

private

Definition at line 53 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and FlatEGunASCIIWriter().

◆ fMinEta

double edm::FlatEGunASCIIWriter::fMinEta

private

Definition at line 49 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and FlatEGunASCIIWriter().

◆ fMinPhi

double edm::FlatEGunASCIIWriter::fMinPhi

private

Definition at line 51 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and FlatEGunASCIIWriter().

◆ fOutFileName

std::string edm::FlatEGunASCIIWriter::fOutFileName

private

Definition at line 62 of file FlatEGunASCIIWriter.h.

Referenced by beginJob().

◆ fOutStream

HepMC::IO_GenEvent* edm::FlatEGunASCIIWriter::fOutStream

private

Definition at line 63 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), beginJob(), and ~FlatEGunASCIIWriter().

◆ fPartIDs

std::vector<int> edm::FlatEGunASCIIWriter::fPartIDs

private

Definition at line 48 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and FlatEGunASCIIWriter().

◆ fPDGTable

ESHandle<HepPDT::ParticleDataTable> edm::FlatEGunASCIIWriter::fPDGTable

private

Definition at line 60 of file FlatEGunASCIIWriter.h.

Referenced by analyze(), and beginRun().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FlatEGunASCIIWriter()

◆ ~FlatEGunASCIIWriter()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

Member Data Documentation

◆ fCurrentEvent

◆ fEvt

◆ fMaxE

◆ fMaxEta

◆ fMaxPhi

◆ fMinE

◆ fMinEta

◆ fMinPhi

◆ fOutFileName

◆ fOutStream

◆ fPartIDs

◆ fPDGTable