RHStopTracer Class Reference

#include <RHStopTracer.h>

Inheritance diagram for RHStopTracer:

Classes
struct	StopPoint

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override
	RHStopTracer (edm::ParameterSet const &p)
void	update (const BeginOfRun *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfEvent *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...
	~RHStopTracer () override
Public Member Functions inherited from SimProducer
const SimProducer &	operator= (const SimProducer &)=delete
void	registerProducts (edm::ProducesCollector producesCollector)
	SimProducer ()
	SimProducer (const SimProducer &)=delete
Public Member Functions inherited from SimWatcher
virtual void	beginRun (edm::EventSetup const &)
bool	isMT () const
const SimWatcher &	operator= (const SimWatcher &)=delete
virtual void	registerConsumes (edm::ConsumesCollector)
	SimWatcher ()
	SimWatcher (const SimWatcher &)=delete
virtual	~SimWatcher ()
Public Member Functions inherited from Observer< const BeginOfRun *>
	Observer ()
void	slotForUpdate (const BeginOfRun * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfTrack *>
	Observer ()
void	slotForUpdate (const BeginOfTrack * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const EndOfTrack *>
	Observer ()
void	slotForUpdate (const EndOfTrack * iT)
virtual	~Observer ()

Private Attributes
int	maxPdgId
int	minPdgId
std::vector< StopPoint >	mStopPoints
bool	mStopRegular
double	mTraceEnergy
std::string	mTraceParticleName
int	otherPdgId
std::regex	rePartName

Additional Inherited Members
Protected Member Functions inherited from SimProducer
template<class T >
void	produces ()
template<class T >
void	produces (const std::string &instanceName)
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 18 of file RHStopTracer.h.

Constructor & Destructor Documentation

RHStopTracer()

RHStopTracer::RHStopTracer

(

edm::ParameterSet const &

p

)

Definition at line 20 of file RHStopTracer.cc.

References maxPdgId, minPdgId, mStopRegular, mTraceEnergy, mTraceParticleName, otherPdgId, AlCaHLTBitMon_ParallelJobs::p, rePartName, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                   {
  edm::ParameterSet parameters = p.getParameter<edm::ParameterSet>("RHStopTracer");
  mStopRegular = parameters.getUntrackedParameter<bool>("stopRegularParticles", false);
  mTraceEnergy = parameters.getUntrackedParameter<double>("traceEnergy", 1.e20);
  mTraceParticleName = parameters.getParameter<std::string>("traceParticle");
  minPdgId = parameters.getUntrackedParameter<int>("minPdgId", 1000000);
  maxPdgId = parameters.getUntrackedParameter<int>("maxPdgId", 2000000);
  otherPdgId = parameters.getUntrackedParameter<int>("otherPdgId", 17);
  produces<std::vector<std::string> >("StoppedParticlesName");
  produces<std::vector<float> >("StoppedParticlesX");
  produces<std::vector<float> >("StoppedParticlesY");
  produces<std::vector<float> >("StoppedParticlesZ");
  produces<std::vector<float> >("StoppedParticlesTime");
  produces<std::vector<int> >("StoppedParticlesPdgId");
  produces<std::vector<float> >("StoppedParticlesMass");
  produces<std::vector<float> >("StoppedParticlesCharge");

  mTraceEnergy *= CLHEP::GeV;
  rePartName = mTraceParticleName;

  edm::LogInfo("SimG4CoreCustomPhysics") << "RHStopTracer::RHStopTracer " << mTraceParticleName
                                         << " Eth(GeV)= " << mTraceEnergy;
}

~RHStopTracer()

RHStopTracer::~RHStopTracer ( )

override

Definition at line 44 of file RHStopTracer.cc.

{}

Member Function Documentation

produce()

void RHStopTracer::produce ( edm::Event &  fEvent, const edm::EventSetup &    )

overridevirtual

Implements SimProducer.

Definition at line 105 of file RHStopTracer.cc.

References CosmicGenFilterHelix_cfi::charges, hcaldqm::fEvent, LogDebug, particleFlowDisplacedVertex_cfi::masses, eostools::move(), mStopPoints, and names.

                                                                 {
  LogDebug("SimG4CoreCustomPhysics") << "RHStopTracer::produce->";

  std::unique_ptr<std::vector<std::string> > names(new std::vector<std::string>);
  std::unique_ptr<std::vector<float> > xs(new std::vector<float>);
  std::unique_ptr<std::vector<float> > ys(new std::vector<float>);
  std::unique_ptr<std::vector<float> > zs(new std::vector<float>);
  std::unique_ptr<std::vector<float> > ts(new std::vector<float>);
  std::unique_ptr<std::vector<int> > ids(new std::vector<int>);
  std::unique_ptr<std::vector<float> > masses(new std::vector<float>);
  std::unique_ptr<std::vector<float> > charges(new std::vector<float>);

  std::vector<StopPoint>::const_iterator stopPoint = mStopPoints.begin();
  for (; stopPoint != mStopPoints.end(); ++stopPoint) {
    names->push_back(stopPoint->name);
    xs->push_back(stopPoint->x);
    ys->push_back(stopPoint->y);
    zs->push_back(stopPoint->z);
    ts->push_back(stopPoint->t);
    ids->push_back(stopPoint->id);
    masses->push_back(stopPoint->mass);
    charges->push_back(stopPoint->charge);
  }
  fEvent.put(std::move(names), "StoppedParticlesName");
  fEvent.put(std::move(xs), "StoppedParticlesX");
  fEvent.put(std::move(ys), "StoppedParticlesY");
  fEvent.put(std::move(zs), "StoppedParticlesZ");
  fEvent.put(std::move(ts), "StoppedParticlesTime");
  fEvent.put(std::move(ids), "StoppedParticlesPdgId");
  fEvent.put(std::move(masses), "StoppedParticlesMass");
  fEvent.put(std::move(charges), "StoppedParticlesCharge");
  mStopPoints.clear();
}

update() [1/4]

void RHStopTracer::update ( const BeginOfRun *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 46 of file RHStopTracer.cc.

References LogDebug.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                {
  LogDebug("SimG4CoreCustomPhysics") << "RHStopTracer::update-> begin of the run " << (*fRun)()->GetRunID();
}

update() [2/4]

void RHStopTracer::update ( const BeginOfEvent *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 50 of file RHStopTracer.cc.

References LogDebug.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                    {
  LogDebug("SimG4CoreCustomPhysics") << "RHStopTracer::update-> begin of the event " << (*fEvent)()->GetEventID();
}

update() [3/4]

void RHStopTracer::update ( const BeginOfTrack *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 54 of file RHStopTracer.cc.

References funct::abs(), LogDebug, muonTagProbeFilters_cff::matched, maxPdgId, minPdgId, mStopRegular, mTraceEnergy, otherPdgId, EgammaValidation_cff::pdgid, rePartName, AlCaHLTBitMon_QueryRunRegistry::string, and HLT_2022v15_cff::track.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                    {
  const G4Track* track = (*fTrack)();
  const G4ParticleDefinition* part = track->GetDefinition();
  const std::string& stringPartName = part->GetParticleName();
  bool matched = false;
  int pdgid = std::abs(part->GetPDGEncoding());
  if ((pdgid > minPdgId && pdgid < maxPdgId) || pdgid == otherPdgId)
    matched = std::regex_match(stringPartName, rePartName);
  if (matched || track->GetKineticEnergy() > mTraceEnergy) {
    LogDebug("SimG4CoreCustomPhysics") << "RHStopTracer::update-> new track: ID/Name/pdgId/mass/charge/Parent: "
                                       << track->GetTrackID() << '/' << part->GetParticleName() << '/'
                                       << part->GetPDGEncoding() << '/' << part->GetPDGMass() / GeV << " GeV/"
                                       << part->GetPDGCharge() << '/' << track->GetParentID()
                                       << " Position: " << track->GetPosition() << ' '
                                       << " R/phi: " << track->GetPosition().perp() << '/' << track->GetPosition().phi()
                                       << "   4vec " << track->GetMomentum();
  } else if (mStopRegular) {  // kill regular particles
    const_cast<G4Track*>(track)->SetTrackStatus(fStopAndKill);
  }
}

update() [4/4]

void RHStopTracer::update ( const EndOfTrack *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack *>.

Definition at line 75 of file RHStopTracer.cc.

References funct::abs(), LogDebug, muonTagProbeFilters_cff::matched, maxPdgId, minPdgId, mStopPoints, mTraceEnergy, otherPdgId, EgammaValidation_cff::pdgid, rePartName, AlCaHLTBitMon_QueryRunRegistry::string, and HLT_2022v15_cff::track.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                  {
  const G4Track* track = (*fTrack)();
  const G4ParticleDefinition* part = track->GetDefinition();
  const std::string& stringPartName = part->GetParticleName();
  bool matched = false;
  int pdgid = std::abs(part->GetPDGEncoding());
  if ((pdgid > minPdgId && pdgid < maxPdgId) || pdgid == otherPdgId)
    matched = std::regex_match(stringPartName, rePartName);
  if (matched || track->GetKineticEnergy() > mTraceEnergy) {
    LogDebug("SimG4CoreCustomPhysics") << "RHStopTracer::update-> stop track: ID/Name/pdgId/mass/charge/Parent: "
                                       << track->GetTrackID() << '/' << part->GetParticleName() << '/'
                                       << part->GetPDGEncoding() << '/' << part->GetPDGMass() / GeV << " GeV/"
                                       << part->GetPDGCharge() << '/' << track->GetParentID()
                                       << " Position: " << track->GetPosition() << ' '
                                       << " R/phi: " << track->GetPosition().perp() << '/' << track->GetPosition().phi()
                                       << "   4vec " << track->GetMomentum();
    if (track->GetMomentum().mag() < 0.001) {
      LogDebug("SimG4CoreCustomPhysics") << "RHStopTracer:: track has stopped, so making StopPoint";
      mStopPoints.push_back(StopPoint(track->GetDefinition()->GetParticleName(),
                                      track->GetPosition().x(),
                                      track->GetPosition().y(),
                                      track->GetPosition().z(),
                                      track->GetGlobalTime(),
                                      track->GetDefinition()->GetPDGEncoding(),
                                      track->GetDefinition()->GetPDGMass() / GeV,
                                      track->GetDefinition()->GetPDGCharge()));
    }
  }
}

Member Data Documentation

maxPdgId

int RHStopTracer::maxPdgId

private

Definition at line 49 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().

minPdgId

int RHStopTracer::minPdgId

private

Definition at line 48 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().

mStopPoints

std::vector<StopPoint> RHStopTracer::mStopPoints

private

Definition at line 53 of file RHStopTracer.h.

Referenced by produce(), and update().

mStopRegular

bool RHStopTracer::mStopRegular

private

Definition at line 46 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().

mTraceEnergy

double RHStopTracer::mTraceEnergy

private

Definition at line 47 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().

mTraceParticleName

std::string RHStopTracer::mTraceParticleName

private

Definition at line 51 of file RHStopTracer.h.

Referenced by RHStopTracer().

otherPdgId

int RHStopTracer::otherPdgId

private

Definition at line 50 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().

rePartName

std::regex RHStopTracer::rePartName

private

Definition at line 52 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().