#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
void	registerProducts (edm::ProducerBase &iProd)

	SimProducer ()

Public Member Functions inherited from SimWatcher
	SimWatcher ()

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)

Detailed Description

Definition at line 18 of file RHStopTracer.h.

Constructor & Destructor Documentation

RHStopTracer::RHStopTracer ( edm::ParameterSet const & p )

Definition at line 20 of file RHStopTracer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), GeV, maxPdgId, minPdgId, mStopRegular, mTraceEnergy, mTraceParticleName, otherPdgId, 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 ( )

override

Definition at line 46 of file RHStopTracer.cc.

46 {

47 }

Member Function Documentation

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

overridevirtual

Implements SimProducer.

Definition at line 115 of file RHStopTracer.cc.

References electrons_cff::ids, LogDebug, particleFlowDisplacedVertex_cfi::masses, eostools::move(), mStopPoints, names, and edm::Event::put().

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

void RHStopTracer::update ( const BeginOfRun * )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 49 of file RHStopTracer.cc.

References LogDebug.

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

void RHStopTracer::update ( const BeginOfEvent * )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 54 of file RHStopTracer.cc.

References LogDebug.

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

void RHStopTracer::update ( const BeginOfTrack * )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 59 of file RHStopTracer.cc.

References funct::abs(), GeV, LogDebug, electrons_cff::matched, maxPdgId, minPdgId, mStopRegular, mTraceEnergy, otherPdgId, BPhysicsValidation_cfi::pdgid, rePartName, AlCaHLTBitMon_QueryRunRegistry::string, and HiIsolationCommonParameters_cff::track.

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

void RHStopTracer::update ( const EndOfTrack * )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 82 of file RHStopTracer.cc.

References funct::abs(), GeV, LogDebug, electrons_cff::matched, maxPdgId, minPdgId, mStopPoints, mTraceEnergy, otherPdgId, BPhysicsValidation_cfi::pdgid, rePartName, AlCaHLTBitMon_QueryRunRegistry::string, and HiIsolationCommonParameters_cff::track.

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