#include <HepMCFileReader.h>

Public Member Functions
HepMC::GenEvent *	fillCurrentEventData ()

virtual void	getStatsFromTuple (int &mo1, int &mo2, int &da1, int &da2, int &status, int &pid, int j) const

virtual void	initialize (const std::string &filename)

bool	isInitialized () const

void	printEvent () const

virtual bool	printHepMcEvent () const

virtual bool	readCurrentEvent ()

virtual void	ReadStats ()

virtual bool	setEvent (int event)

virtual	~HepMCFileReader ()

Static Public Member Functions
static HepMCFileReader *	instance ()

Protected Member Functions
	HepMCFileReader ()

Private Member Functions
int	find_in_map (const std::map< HepMC::GenParticle , int > &m, HepMC::GenParticle p) const

HepMC::IO_BaseClass const *	input () const

HepMC::IO_BaseClass *&	input ()

int	rdstate () const

Private Attributes
edm::propagate_const< HepMC::GenEvent * >	evt_

std::vector< HepMC::GenParticle * >	index_to_particle

edm::propagate_const< HepMC::IO_BaseClass * >	input_

std::map< HepMC::GenParticle *, int >	particle_to_index

Static Private Attributes
static HepMCFileReader *	instance_ = nullptr

Detailed Description

This class is used by the implementation of DaqEventFactory present in this package to read in the full event raw data from a flat binary file. WARNING: If you want to use this class for other purposes you must always invoke the method initialize before starting using the interface it exposes.

Author: G. Bruno - CERN, EP Division

Definition at line 27 of file HepMCFileReader.h.

Constructor & Destructor Documentation

◆ HepMCFileReader()

HepMCFileReader::HepMCFileReader ( )

protected

Definition at line 38 of file HepMCFileReader.cc.

References instance_.

                                  : evt_(nullptr), input_(nullptr) {
   // Default constructor.
   if (instance_ == nullptr) {
     instance_ = this;
   } else {
     edm::LogError("HepMCFileReader") << "Constructing a new instance";
   }
 }

◆ ~HepMCFileReader()

HepMCFileReader::~HepMCFileReader ( )

virtual

Definition at line 48 of file HepMCFileReader.cc.

References edm::propagate_const< T >::get(), input_, and instance_.

                                   {
   edm::LogInfo("HepMCFileReader") << "Destructing HepMCFileReader";
 
   instance_ = nullptr;
   delete input_.get();
 }

Member Function Documentation

◆ fillCurrentEventData()

HepMC::GenEvent * HepMCFileReader::fillCurrentEventData ( )

Definition at line 108 of file HepMCFileReader.cc.

References evt_, and readCurrentEvent().

                                                      {
   readCurrentEvent();
   return evt_;
 }

◆ find_in_map()

int HepMCFileReader::find_in_map	(	const std::map< HepMC::GenParticle *, int > &	m,
		HepMC::GenParticle *	p
	)		const

private

Definition at line 252 of file HepMCFileReader.cc.

References visualization-live-secondInstance_cfg::m, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by getStatsFromTuple().

                                                                                                     {
   std::map<HepMC::GenParticle *, int>::const_iterator iter = m.find(p);
   return (iter == m.end()) ? 0 : iter->second;
 }

◆ getStatsFromTuple()

void HepMCFileReader::getStatsFromTuple	(	int &	mo1,
		int &	mo2,
		int &	da1,
		int &	da2,
		int &	status,
		int &	pid,
		int	j
	)		const

virtual

Definition at line 205 of file HepMCFileReader.cc.

References class-composition::children, gather_cfg::cout, evt_, find_in_map(), index_to_particle, dqmiolumiharvest::j, particle_to_index, and mps_update::status.

Referenced by printEvent().

                                                                                                                   {
   if (!evt_) {
     cout << "HepMCFileReader: Got no event :-(  Game over already  ?" << endl;
   } else {
     status = index_to_particle[j]->status();
     pid = index_to_particle[j]->pdg_id();
     if (index_to_particle[j]->production_vertex()) {
       //HepLorentzVector p = index_to_particle[j]->
       //production_vertex()->position();
 
       int num_mothers = index_to_particle[j]->production_vertex()->particles_in_size();
       int first_mother =
           find_in_map(particle_to_index, *(index_to_particle[j]->production_vertex()->particles_in_const_begin()));
       int last_mother = first_mother + num_mothers - 1;
       if (first_mother == 0)
         last_mother = 0;
       mo1 = first_mother;
       mo2 = last_mother;
     } else {
       mo1 = 0;
       mo2 = 0;
     }
     if (!index_to_particle[j]->end_vertex()) {
       //find # of 1. daughter
       int first_daughter =
           find_in_map(particle_to_index, *(index_to_particle[j]->end_vertex()->particles_begin(HepMC::children)));
       //cout <<"first_daughter "<< first_daughter <<  "num_daughters " << num_daughters << endl;
       HepMC::GenVertex::particle_iterator ic;
       int last_daughter = 0;
       //find # of last daughter
       for (ic = index_to_particle[j]->end_vertex()->particles_begin(HepMC::children);
            ic != index_to_particle[j]->end_vertex()->particles_end(HepMC::children);
            ++ic)
         last_daughter = find_in_map(particle_to_index, *ic);
 
       if (first_daughter == 0)
         last_daughter = 0;
       da1 = first_daughter;
       da2 = last_daughter;
     } else {
       da1 = 0;
       da2 = 0;
     }
   }
 }

◆ initialize()

void HepMCFileReader::initialize ( const std::string & filename )

virtual

Definition at line 56 of file HepMCFileReader.cc.

References Exception, corrVsCorr::filename, edm::propagate_const< T >::get(), recoMuon::in, input_, isInitialized(), and rdstate().

                                                        {
   if (isInitialized()) {
     edm::LogError("HepMCFileReader") << "Was already initialized... reinitializing";
     delete input_.get();
   }
 
   edm::LogInfo("HepMCFileReader") << "Opening file" << filename << "using HepMC::IO_GenEvent";
   input_ = new HepMC::IO_GenEvent(filename, std::ios::in);
 
   if (rdstate() == std::ios::failbit) {
     throw cms::Exception("FileNotFound", "HepMCFileReader::initialize()") << "File " << filename << " was not found.\n";
   }
 }

◆ input() [1/2]

HepMC::IO_BaseClass const* HepMCFileReader::input ( ) const

inlineprivate

Definition at line 53 of file HepMCFileReader.h.

References edm::get_underlying_safe(), and input_.

Referenced by rdstate().

53 { return get_underlying_safe(input_); }

HepMCFileReader::input_

edm::propagate_const< HepMC::IO_BaseClass * > input_

Definition: HepMCFileReader.h:58

edm::get_underlying_safe

constexpr std::shared_ptr< T > & get_underlying_safe(propagate_const< std::shared_ptr< T >> &iP)

Definition: get_underlying_safe.h:41

◆ input() [2/2]

HepMC::IO_BaseClass*& HepMCFileReader::input ( )

inlineprivate

Definition at line 54 of file HepMCFileReader.h.

References edm::get_underlying_safe(), and input_.

54 { return get_underlying_safe(input_); }

HepMCFileReader::input_

edm::propagate_const< HepMC::IO_BaseClass * > input_

Definition: HepMCFileReader.h:58

edm::get_underlying_safe

constexpr std::shared_ptr< T > & get_underlying_safe(propagate_const< std::shared_ptr< T >> &iP)

Definition: get_underlying_safe.h:41

◆ instance()

HepMCFileReader * HepMCFileReader::instance ( )

static

Definition at line 28 of file HepMCFileReader.cc.

                                            {
   // Implement a HepMCFileReader singleton.
 
   if (instance_ == nullptr) {
     instance_ = new HepMCFileReader();
   }
   return instance_;
 }

◆ isInitialized()

bool HepMCFileReader::isInitialized ( ) const

inline

Definition at line 71 of file HepMCFileReader.h.

References input_.

Referenced by initialize().

71 { return input_ != nullptr; }

HepMCFileReader::input_

edm::propagate_const< HepMC::IO_BaseClass * > input_

Definition: HepMCFileReader.h:58

◆ printEvent()

void HepMCFileReader::printEvent ( ) const

Definition at line 115 of file HepMCFileReader.cc.

References gather_cfg::cout, evt_, alignBH_cfg::fixed, g, GenParticle::GenParticle, getStatsFromTuple(), index_to_particle, dqmiodumpmetadata::n, mathSSE::sqrt(), and mps_update::status.

                                        {
   int mo1 = 0, mo2 = 0, da1 = 0, da2 = 0, status = 0, pid = 0;
   if (evt_ != nullptr) {
     cout << "---#-------pid--st---Mo1---Mo2---Da1---Da2------px------py------pz-------E-";
     cout << "------m---------x---------y---------z---------t-";
     cout << endl;
     cout.setf(ios::right, ios::adjustfield);
     for (int n = 1; n <= evt_->particles_size(); n++) {
       HepMC::GenParticle *g = index_to_particle[n];
       getStatsFromTuple(mo1, mo2, da1, da2, status, pid, n);
       cout << setw(4) << n << setw(8) << pid << setw(5) << status << setw(6) << mo1 << setw(6) << mo2 << setw(6) << da1
            << setw(6) << da2;
       cout.setf(ios::fixed, ios::floatfield);
       cout.setf(ios::right, ios::adjustfield);
       cout << setw(10) << setprecision(2) << g->momentum().x();
       cout << setw(8) << setprecision(2) << g->momentum().y();
       cout << setw(10) << setprecision(2) << g->momentum().z();
       cout << setw(8) << setprecision(2) << g->momentum().t();
       cout << setw(8) << setprecision(2) << g->generatedMass();
       // tau=L/(gamma*beta*c)
       if (g->production_vertex() != nullptr && g->end_vertex() != nullptr && status == 2) {
         cout << setw(10) << setprecision(2) << g->production_vertex()->position().x();
         cout << setw(10) << setprecision(2) << g->production_vertex()->position().y();
         cout << setw(10) << setprecision(2) << g->production_vertex()->position().z();
 
         double xm = g->production_vertex()->position().x();
         double ym = g->production_vertex()->position().y();
         double zm = g->production_vertex()->position().z();
         double xd = g->end_vertex()->position().x();
         double yd = g->end_vertex()->position().y();
         double zd = g->end_vertex()->position().z();
         double decl = sqrt((xd - xm) * (xd - xm) + (yd - ym) * (yd - ym) + (zd - zm) * (zd - zm));
         double labTime = decl / c_light;
         // convert lab time to proper time
         double properTime = labTime / g->momentum().rho() * (g->generatedMass());
         // set the proper time in nanoseconds
         cout << setw(8) << setprecision(2) << properTime;
       } else {
         cout << setw(10) << setprecision(2) << 0.0;
         cout << setw(10) << setprecision(2) << 0.0;
         cout << setw(10) << setprecision(2) << 0.0;
         cout << setw(8) << setprecision(2) << 0.0;
       }
       cout << endl;
     }
   } else {
     cout << " HepMCFileReader: No event available !" << endl;
   }
 }

◆ printHepMcEvent()

bool HepMCFileReader::printHepMcEvent ( ) const

virtual

Definition at line 101 of file HepMCFileReader.cc.

References evt_.

                                             {
   if (evt_ != nullptr)
     evt_->print();
   return true;
 }

◆ rdstate()

int HepMCFileReader::rdstate ( ) const

private

Definition at line 71 of file HepMCFileReader.cc.

References input(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by initialize().

                                    {
   // work around a HepMC IO_ inheritence shortfall
 
   HepMC::IO_GenEvent const *p = dynamic_cast<HepMC::IO_GenEvent const *>(input());
   if (p)
     return p->rdstate();
 
   return std::ios::failbit;
 }

◆ readCurrentEvent()

bool HepMCFileReader::readCurrentEvent ( )

virtual

Definition at line 82 of file HepMCFileReader.cc.

References evt_, input_, and ReadStats().

Referenced by fillCurrentEventData().

                                        {
   evt_ = input_->read_next_event();
   if (evt_) {
     edm::LogInfo("HepMCFileReader") << "| --- Event Nr. " << evt_->event_number() << " with " << evt_->particles_size()
                                     << " particles --- !";
     ReadStats();
     //  printHepMcEvent();
     //          printEvent();
   } else {
     edm::LogInfo("HepMCFileReader") << "Got no event" << endl;
   }
 
   return evt_ != nullptr;
 }

◆ ReadStats()

void HepMCFileReader::ReadStats ( )

virtual

Definition at line 166 of file HepMCFileReader.cc.

References evt_, index_to_particle, LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, particle_to_index, and findQualityFiles::v.

Referenced by readCurrentEvent().

                                 {
   unsigned int particle_counter = 0;
   index_to_particle.reserve(evt_->particles_size() + 1);
   index_to_particle[0] = nullptr;
   for (HepMC::GenEvent::vertex_const_iterator v = evt_->vertices_begin(); v != evt_->vertices_end(); ++v) {
     // making a list of incoming particles of the vertices
     // so that the mother indices in HEPEVT can be filled properly
     for (HepMC::GenVertex::particles_in_const_iterator p1 = (*v)->particles_in_const_begin();
          p1 != (*v)->particles_in_const_end();
          ++p1) {
       ++particle_counter;
       //particle_counter can be very large for heavy ions
       if (particle_counter >= index_to_particle.size()) {
         //make it large enough to hold up to this index
         index_to_particle.resize(particle_counter + 1);
       }
       index_to_particle[particle_counter] = *p1;
       particle_to_index[*p1] = particle_counter;
     }
     // daughters are entered only if they aren't a mother of
     // another vertex
     for (HepMC::GenVertex::particles_out_const_iterator p2 = (*v)->particles_out_const_begin();
          p2 != (*v)->particles_out_const_end();
          ++p2) {
       if (!(*p2)->end_vertex()) {
         ++particle_counter;
         //particle_counter can be very large for heavy ions
         if (particle_counter >= index_to_particle.size()) {
           //make it large enough to hold up to this index
           index_to_particle.resize(particle_counter + 1);
         }
         index_to_particle[particle_counter] = *p2;
         particle_to_index[*p2] = particle_counter;
       }
     }
   }
 }

◆ setEvent()

bool HepMCFileReader::setEvent ( int event )

virtual

Definition at line 98 of file HepMCFileReader.cc.

98 { return true; }

Member Data Documentation

◆ evt_

edm::propagate_const<HepMC::GenEvent *> HepMCFileReader::evt_

private

Definition at line 57 of file HepMCFileReader.h.

Referenced by fillCurrentEventData(), getStatsFromTuple(), printEvent(), printHepMcEvent(), readCurrentEvent(), and ReadStats().

◆ index_to_particle

std::vector<HepMC::GenParticle *> HepMCFileReader::index_to_particle

private

Definition at line 65 of file HepMCFileReader.h.

Referenced by getStatsFromTuple(), printEvent(), and ReadStats().

◆ input_

edm::propagate_const<HepMC::IO_BaseClass *> HepMCFileReader::input_

private

Definition at line 58 of file HepMCFileReader.h.

Referenced by initialize(), input(), isInitialized(), readCurrentEvent(), and ~HepMCFileReader().

◆ instance_

HepMCFileReader * HepMCFileReader::instance_ = nullptr

staticprivate

Definition at line 60 of file HepMCFileReader.h.

Referenced by HepMCFileReader(), and ~HepMCFileReader().

◆ particle_to_index

std::map<HepMC::GenParticle *, int> HepMCFileReader::particle_to_index

private

Definition at line 66 of file HepMCFileReader.h.

Referenced by getStatsFromTuple(), and ReadStats().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HepMCFileReader()

◆ ~HepMCFileReader()

Member Function Documentation

◆ fillCurrentEventData()

◆ find_in_map()

◆ getStatsFromTuple()

◆ initialize()

◆ input() [1/2]

◆ input() [2/2]

◆ instance()

◆ isInitialized()

◆ printEvent()

◆ printHepMcEvent()

◆ rdstate()

◆ readCurrentEvent()

◆ ReadStats()

◆ setEvent()

Member Data Documentation

◆ evt_

◆ index_to_particle

◆ input_

◆ instance_

◆ particle_to_index