RHStopTracer Class Reference

#include <RHStopTracer.h>

Inheritance diagram for RHStopTracer:

Classes
struct	StopPoint

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &)
	RHStopTracer (edm::ParameterSet const &p)
void	update (const BeginOfRun *)
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfEvent *)
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfTrack *)
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfTrack *)
	This routine will be called when the appropriate signal arrives. More...
virtual	~RHStopTracer ()
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
const G4ParticleDefinition *	mParticle
std::vector< StopPoint >	mStopPoints
bool	mStopRegular
double	mTraceEnergy
std::string	mTraceParticleName

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 16 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, mParticle, mStopRegular, mTraceEnergy, mTraceParticleName, metProducer_cfi::parameters, 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");
  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");

  mParticle = G4ParticleTable::GetParticleTable()->FindParticle(mTraceParticleName);

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

RHStopTracer::~RHStopTracer ( )

virtual

Definition at line 42 of file RHStopTracer.cc.

                            {
}

Member Function Documentation

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

virtual

Implements SimProducer.

Definition at line 99 of file RHStopTracer.cc.

References LogDebug, particleFlowDisplacedVertex_cfi::masses, eostools::move(), mStopPoints, cscdqm::h::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 *  )

virtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 45 of file RHStopTracer.cc.

References LogDebug.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

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

void RHStopTracer::update ( const BeginOfEvent *  )

virtual

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__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

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

void RHStopTracer::update ( const BeginOfTrack *  )

virtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 55 of file RHStopTracer.cc.

References GeV, LogDebug, mParticle, mStopRegular, mTraceEnergy, and HiIsolationCommonParameters_cff::track.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                      {
  const G4Track* track = (*fTrack)();
  const G4ParticleDefinition* part = track->GetDefinition();
  if(part && part == mParticle &&  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 *  )

virtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 73 of file RHStopTracer.cc.

References GeV, LogDebug, mParticle, mStopPoints, mTraceEnergy, and HiIsolationCommonParameters_cff::track.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                    {
  const G4Track* track = (*fTrack)();
  const G4ParticleDefinition* part = track->GetDefinition();
  if(part && part == mParticle &&  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) {
      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

const G4ParticleDefinition* RHStopTracer::mParticle

private

Definition at line 47 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().

std::vector<StopPoint> RHStopTracer::mStopPoints

private

Definition at line 48 of file RHStopTracer.h.

Referenced by produce(), and update().

bool RHStopTracer::mStopRegular

private

Definition at line 44 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().

double RHStopTracer::mTraceEnergy

private

Definition at line 45 of file RHStopTracer.h.

Referenced by RHStopTracer(), and update().

std::string RHStopTracer::mTraceParticleName

private

Definition at line 46 of file RHStopTracer.h.

Referenced by RHStopTracer().