CDCPathBasicVarSet Class Reference

Var set to store basic quantities related to CDC CKF. More...

#include <CDCPathBasicVarSet.h>

Inheritance diagram for CDCPathBasicVarSet:

Public Types
using	Object
	Type from which variables should be extracted.

Public Member Functions
virtual bool	extract (const BaseCDCPathFilter::Object *path) override
	Generate and assign the variables from the object.
std::vector< Named< Float_t * > >	getNamedVariables (const std::string &prefix) override
	Getter for the named references to the individual variables Base implementation returns empty vector.
std::vector< Named< Float_t * > >	getNamedVariables ()
	Getter for the named references to the individual variables.
MayBePtr< Float_t >	find (const std::string &varName) override
	Pointer to the variable with the given name.
virtual bool	extract (const Object *obj)
	Main method that extracts the variable values from the complex object.
bool	extract (const Object &obj)
	Method for extraction from an object instead of a pointer.
std::map< std::string, Float_t >	getNamedValues (const std::string &prefix) const
	Getter for a map of names to float values.
std::map< std::string, Float_t >	getNamedValues () const
	Getter for a map of names to float values.
void	initialize () override
	Receive and dispatch signal before the start of the event processing.
void	beginRun () override
	Receive and dispatch signal for the beginning of a new run.
void	beginEvent () override
	Receive and dispatch signal for the start of a new event.
void	endRun () override
	Receive and dispatch signal for the end of the run.
void	terminate () override
	Receive and dispatch Signal for termination of the event processing.
virtual void	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix)
	Expose the set of parameters of the filter to the module parameter list.

Protected Member Functions
Float_t	get () const
	Getter for the value of the ith variable. Static version.
Float_t &	var ()
	Reference getter for the value of the ith variable. Static version.
AssignFinite< Float_t >	finitevar ()
	Reference getter for the value of the ith variable. Transforms non-finite values to finite value.
void	addProcessingSignalListener (ProcessingSignalListener *psl)
	Register a processing signal listener to be notified.
int	getNProcessingSignalListener ()
	Get the number of currently registered listeners.

Static Protected Member Functions
static constexpr int	named (const char *name)
	Getter for the index from the name.

Private Types
using	Super
	Type of the super class.

Private Attributes
StoreObjPtr< EventMetaData >	m_eventMetaData
	use this to get the eventNumber
FixedSizeNamedFloatTuple< CDCPathBasicVarNames >	m_variables
	Memory for nVars floating point values.
std::vector< ProcessingSignalListener * >	m_subordinaryProcessingSignalListeners
	References to subordinary signal processing listener contained in this findlet.
bool	m_initialized
	Flag to keep track whether initialization happened before.
bool	m_terminated
	Flag to keep track whether termination happened before.
std::string	m_initializedAs
	Name of the type during initialisation.

Static Private Attributes
static const size_t	nVars
	Number of floating point values represented by this class.

Detailed Description

Var set to store basic quantities related to CDC CKF.

Definition at line 104 of file CDCPathBasicVarSet.h.

Member Typedef Documentation

Object

using Object

inherited

Type from which variables should be extracted.

Definition at line 44 of file VarSet.h.

Super

using Super

privateinherited

Type of the super class.

Definition at line 40 of file VarSet.h.

Member Function Documentation

addProcessingSignalListener()

void addProcessingSignalListener ( ProcessingSignalListener * psl )

protectedinherited

Register a processing signal listener to be notified.

Definition at line 53 of file CompositeProcessingSignalListener.cc.

{
  m_subordinaryProcessingSignalListeners.push_back(psl);
}

beginEvent()

void beginEvent ( )

overrideinherited

Receive and dispatch signal for the start of a new event.

Definition at line 36 of file CompositeProcessingSignalListener.cc.

{
  Super::beginEvent();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->beginEvent();
  }
}

beginRun()

void beginRun ( )

overrideinherited

Receive and dispatch signal for the beginning of a new run.

Definition at line 33 of file CompositeProcessingSignalListener.cc.

{
  Super::beginRun();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->beginRun();
  }
}

endRun()

void endRun ( )

overrideinherited

Receive and dispatch signal for the end of the run.

Definition at line 39 of file CompositeProcessingSignalListener.cc.

{
  for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
    psl->endRun();
  }
  Super::endRun();
}

exposeParameters()

virtual void exposeParameters ( ModuleParamList * moduleParamList, const std::string & prefix )

inlinevirtualinherited

Expose the set of parameters of the filter to the module parameter list.

Parameters

moduleParamList	Module parameter list to add the new parameters to
prefix	Prefix (or sometimes postfix) to be added to all nely added parameters.

Definition at line 47 of file CompositeProcessingSignalListener.h.

      {};

extract() [1/3]

bool extract ( const BaseCDCPathFilter::Object * path )

overridevirtual

Generate and assign the variables from the object.

Definition at line 21 of file CDCPathBasicVarSet.cc.

{
  // check if hit belongs to same seed
  const auto& seed = path->front();
  const auto* seedRecoTrack = seed.getSeed();
 
  double arcLength_total = 0.;
  std::vector<double> arcLengths;
  arcLengths.reserve(path->size() - 1);
  std::vector<double> hitDistances;
  hitDistances.reserve(path->size() - 1);
  std::vector<double> flipPos(8, 0);
 
  unsigned int chargeFlip = 0;
  int lastCharge = seedRecoTrack->getChargeSeed();
 
  for (auto const& state : *path) {
    if (state.isSeed()) {
      continue;
    }
 
    arcLengths.push_back(state.getArcLength() - arcLength_total);
    arcLength_total = state.getArcLength();
 
    hitDistances.push_back(state.getHitDistance());
 
 
    // check how many times genfit changed the charge and where it occurred
    genfit::MeasuredStateOnPlane trackState = state.getTrackState();
    int stateCharge = trackState.getCharge();
    if (lastCharge != stateCharge) {
      if (chargeFlip < flipPos.size()) {
        flipPos[chargeFlip] = state.getWireHit()->getWire().getICLayer();
      }
 
      chargeFlip += 1;
      lastCharge = stateCharge;
    }
  }
 
  // general stuff
  var<named("eventNumber")>() = m_eventMetaData->getEvent();
 
  // track properties
  var<named("totalHits")>() = path->size() - 1;
  var<named("chargeFlip")>() = chargeFlip;
 
  bool reachedEnd = false;
  if (path->size() > 1) {
    if (path->back().getWireHit()->getWire().getICLayer() == 0) {
      reachedEnd = true;
    }
  }
  var<named("reachedEnd")>() = reachedEnd;
 
  var<named("flipPos0")>() = flipPos[0];
  var<named("flipPos1")>() = flipPos[1];
  var<named("flipPos2")>() = flipPos[2];
  var<named("flipPos3")>() = flipPos[3];
 
  // seed properties
  ROOT::Math::XYZVector seedPos = seedRecoTrack->getPositionSeed();
  ROOT::Math::XYZVector seedMom = seedRecoTrack->getMomentumSeed();
  var<named("seed_r")>() = seedPos.Rho();
  var<named("seed_z")>() = seedPos.Z();
  var<named("seed_x")>() = seedPos.X();
  var<named("seed_y")>() = seedPos.Y();
  var<named("seed_p")>() = seedMom.R();
  var<named("seed_theta")>() = seedMom.Theta() * 180. / M_PI;
  var<named("seed_pt")>() = seedMom.Rho();
  var<named("seed_pz")>() = seedMom.Z();
  var<named("seed_px")>() = seedMom.X();
  var<named("seed_py")>() = seedMom.Y();
  var<named("seed_charge")>() = seedRecoTrack->getChargeSeed();
 
  var<named("totalHitsSeedTrack")>() = seedRecoTrack->getNumberOfCDCHits();
 
  // get ICLayer assigned to seed (only really defined for ECL seeds)
  const auto& wireTopology = CDC::CDCWireTopology::getInstance();
  const auto& wires = wireTopology.getWires();
  const float maxForwardZ = wires.back().getForwardZ();     // 157.615
  const float maxBackwardZ = wires.back().getBackwardZ();   // -72.0916
 
  int seedICLayer = -1;
  const float seedPosZ = seedPos.z();
  if (seedPosZ < maxForwardZ && seedPosZ > maxBackwardZ) {
    seedICLayer = 56;
  } else {
    // do straight extrapolation of seed momentum to CDC outer walls
    ROOT::Math::XYZVector seedMomZOne(seedMom * (1. / seedMom.Z()));
 
    // find closest iCLayer
    float minDist = 99999;
    for (const auto& wire : wires) {
      const float maxZ = seedPosZ > 0 ? wire.getForwardZ() : wire.getBackwardZ();
 
      const auto distance = wire.getDistance(ROOT::Math::XYZVector(seedPos - seedMomZOne * (seedPosZ - maxZ)));
      if (distance < minDist) {
        minDist = distance;
        seedICLayer = wire.getICLayer();
      }
    }
  }
  var<named("seedICLayer")>() = seedICLayer;
 
 
  // track representation
  ROOT::Math::XYZVector trackMom(0, 0, 0);
  int trackCharge = 0;
  float firstChi2 = 0;
  float lastChi2 = 0;
  float firstICLayer = 0;
  float lastICLayer = 0;
  if (path->size() > 1) {
    genfit::MeasuredStateOnPlane trackState = path->back().getTrackState();
 
    trackMom = trackState.getMom();
    trackCharge = trackState.getCharge();
 
    firstChi2 = path->at(1).getChi2();
    lastChi2 = path->back().getChi2();
 
    firstICLayer = path->at(1).getWireHit()->getWire().getICLayer();
    lastICLayer = path->back().getWireHit()->getWire().getICLayer();
  }
  var<named("track_theta")>() = trackMom.Theta() * 180. / M_PI;
  var<named("track_p")>() = trackMom.R();
  var<named("track_pt")>() = trackMom.Rho();
  var<named("track_pz")>() = trackMom.Z();
  var<named("track_px")>() = trackMom.X();
  var<named("track_py")>() = trackMom.Y();
  var<named("track_charge")>() = trackCharge;
 
  var<named("firstChi2")>() = firstChi2;
  var<named("lastChi2")>() = lastChi2;
 
  var<named("firstICLayer")>() = firstICLayer;
  var<named("lastICLayer")>() = lastICLayer;
 
  if (path->size() > 3) {
    int lastN = path->size() - 1;
    var<named("ICLayerLast0")>() = path->at(lastN).getWireHit()->getWire().getICLayer();
    var<named("ICLayerLast1")>() = path->at(lastN - 1).getWireHit()->getWire().getICLayer();
    var<named("ICLayerLast2")>() = path->at(lastN - 2).getWireHit()->getWire().getICLayer();
    var<named("IWireLast0")>() = path->at(lastN).getWireHit()->getWire().getIWire();
    var<named("IWireLast1")>() = path->at(lastN - 1).getWireHit()->getWire().getIWire();
    var<named("IWireLast2")>() = path->at(lastN - 2).getWireHit()->getWire().getIWire();
  } else {
    var<named("ICLayerLast0")>() = -1;
    var<named("ICLayerLast1")>() = -1;
    var<named("ICLayerLast2")>() = -1;
    var<named("IWireLast0")>() = -1;
    var<named("IWireLast1")>() = -1;
    var<named("IWireLast2")>() = -1;
  }
 
  if (path->size() > 10) {
    var<named("hitDistance0")>() = hitDistances[0];
    var<named("hitDistance1")>() = hitDistances[1];
    var<named("hitDistance2")>() = hitDistances[2];
    var<named("hitDistance3")>() = hitDistances[3];
    var<named("hitDistance4")>() = hitDistances[4];
    var<named("hitDistance5")>() = hitDistances[5];
    var<named("hitDistance6")>() = hitDistances[6];
    var<named("hitDistance7")>() = hitDistances[7];
    var<named("hitDistance8")>() = hitDistances[8];
    var<named("hitDistance9")>() = hitDistances[9];
    var<named("arcLength0")>() = arcLengths[0];
    var<named("arcLength1")>() = arcLengths[1];
    var<named("arcLength2")>() = arcLengths[2];
    var<named("arcLength3")>() = arcLengths[3];
    var<named("arcLength4")>() = arcLengths[4];
    var<named("arcLength5")>() = arcLengths[5];
    var<named("arcLength6")>() = arcLengths[6];
    var<named("arcLength7")>() = arcLengths[7];
    var<named("arcLength8")>() = arcLengths[8];
    var<named("arcLength9")>() = arcLengths[9];
  } else {
    var<named("hitDistance0")>() = -1;
    var<named("hitDistance1")>() = -1;
    var<named("hitDistance2")>() = -1;
    var<named("hitDistance3")>() = -1;
    var<named("hitDistance4")>() = -1;
    var<named("hitDistance5")>() = -1;
    var<named("hitDistance6")>() = -1;
    var<named("hitDistance7")>() = -1;
    var<named("hitDistance8")>() = -1;
    var<named("hitDistance9")>() = -1;
    var<named("arcLength0")>() = -1;
    var<named("arcLength1")>() = -1;
    var<named("arcLength2")>() = -1;
    var<named("arcLength3")>() = -1;
    var<named("arcLength4")>() = -1;
    var<named("arcLength5")>() = -1;
    var<named("arcLength6")>() = -1;
    var<named("arcLength7")>() = -1;
    var<named("arcLength8")>() = -1;
    var<named("arcLength9")>() = -1;
  }
 
  return true;
}

extract() [2/3]

bool extract ( const Object & obj )

inlineinherited

Method for extraction from an object instead of a pointer.

Definition at line 57 of file BaseVarSet.h.

      {
        return extract(&obj);
      }

extract() [3/3]

virtual bool extract ( const Object * obj )

inlinevirtualinherited

Main method that extracts the variable values from the complex object.

Returns indication whether the extraction could be completed successfully. Base implementation returns always true.

Parameters

obj	dummy object not used

Returns

: always true for the base implementation

Definition at line 51 of file BaseVarSet.h.

      {
        return true;
      }

find()

MayBePtr< Float_t > find ( const std::string & varName )

inlineoverridevirtualinherited

Pointer to the variable with the given name.

Returns nullptr if not found.

Reimplemented from BaseVarSet< CDCPathBasicVarNames::Object >.

Definition at line 64 of file VarSet.h.

      {
        return m_variables.find(varName);
      }

finitevar()

AssignFinite< Float_t > finitevar ( )

inlineprotectedinherited

Reference getter for the value of the ith variable. Transforms non-finite values to finite value.

Definition at line 130 of file VarSet.h.

      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return AssignFinite<Float_t>(m_variables[I]);
      }

get()

Float_t get ( ) const

inlineprotectedinherited

Getter for the value of the ith variable. Static version.

Definition at line 85 of file VarSet.h.

      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return m_variables.get(I);
      }

getNamedValues() [1/2]

std::map< std::string, Float_t > getNamedValues ( ) const

inlineinherited

Getter for a map of names to float values.

Definition at line 96 of file BaseVarSet.h.

      {
        const std::string prefix = "";
        return this->getNamedValues(prefix);
      }

getNamedValues() [2/2]

std::map< std::string, Float_t > getNamedValues ( const std::string & prefix ) const

inlineinherited

Getter for a map of names to float values.

Parameters

prefix

Name prefix to apply to all variable names.

Definition at line 84 of file BaseVarSet.h.

      {
        std::map<std::string, Float_t> result;
        std::vector<Named<Float_t*> > namedVariables = this->getNamedVariables(prefix);
        for (const Named<Float_t*>& namedVariable : namedVariables) {
          const Float_t* variable = namedVariable;
          result[namedVariable.getName()] = *variable;
        }
        return result;
      }

getNamedVariables() [1/2]

std::vector< Named< Float_t * > > getNamedVariables ( )

inlineinherited

Getter for the named references to the individual variables.

Definition at line 74 of file BaseVarSet.h.

      {
        const std::string prefix = "";
        return this->getNamedVariables(prefix);
      }

getNamedVariables() [2/2]

std::vector< Named< Float_t * > > getNamedVariables ( const std::string & prefix )

inlineoverridevirtualinherited

Getter for the named references to the individual variables Base implementation returns empty vector.

Reimplemented from BaseVarSet< CDCPathBasicVarNames::Object >.

Definition at line 55 of file VarSet.h.

      {
        return m_variables.getNamedVariables(prefix);
      }

getNProcessingSignalListener()

int getNProcessingSignalListener ( )

protectedinherited

Get the number of currently registered listeners.

Definition at line 56 of file CompositeProcessingSignalListener.cc.

{
  return m_subordinaryProcessingSignalListeners.size();
}

initialize()

void initialize ( )

overrideinherited

Receive and dispatch signal before the start of the event processing.

Definition at line 30 of file CompositeProcessingSignalListener.cc.

{
  Super::initialize();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->initialize();
  }
}

named()

static constexpr int named ( const char * name )

inlinestaticconstexprprotectedinherited

Getter for the index from the name.

Looks through the associated names and returns the right index if found Returns nVars (one after the last element) if not found.

Parameters

name	Name of the sought variable

Returns

Index of the name, nVars if not found.

Definition at line 78 of file VarSet.h.

      {
        return index<nVars>(AVarNames::getName, name);
      }

terminate()

void terminate ( )

overrideinherited

Receive and dispatch Signal for termination of the event processing.

Definition at line 42 of file CompositeProcessingSignalListener.cc.

{
  for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
    psl->terminate();
  }
  Super::terminate();
}

var()

Float_t & var ( )

inlineprotectedinherited

Reference getter for the value of the ith variable. Static version.

Definition at line 93 of file VarSet.h.

      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return m_variables[I];
      }

Member Data Documentation

m_eventMetaData

StoreObjPtr<EventMetaData> m_eventMetaData

private

use this to get the eventNumber

Definition at line 112 of file CDCPathBasicVarSet.h.

m_initialized

bool m_initialized

privateinherited

Flag to keep track whether initialization happened before.

Definition at line 52 of file ProcessingSignalListener.h.

m_initializedAs

std::string m_initializedAs

privateinherited

Name of the type during initialisation.

Definition at line 58 of file ProcessingSignalListener.h.

m_subordinaryProcessingSignalListeners

std::vector<ProcessingSignalListener*> m_subordinaryProcessingSignalListeners

privateinherited

References to subordinary signal processing listener contained in this findlet.

Definition at line 60 of file CompositeProcessingSignalListener.h.

m_terminated

bool m_terminated

privateinherited

Flag to keep track whether termination happened before.

Definition at line 55 of file ProcessingSignalListener.h.

m_variables

FixedSizeNamedFloatTuple<CDCPathBasicVarNames> m_variables

privateinherited

Memory for nVars floating point values.

Definition at line 138 of file VarSet.h.

nVars

const size_t nVars

staticprivateinherited

Number of floating point values represented by this class.

Definition at line 48 of file VarSet.h.

---

The documentation for this class was generated from the following files:

• tracking/ckf/cdc/filters/paths/include/CDCPathBasicVarSet.h
• tracking/ckf/cdc/filters/paths/src/CDCPathBasicVarSet.cc