InDet::SiSpacePointsSeedMaker_LowMomentum Class Reference

#include <SiSpacePointsSeedMaker_LowMomentum.h>

Inheritance diagram for InDet::SiSpacePointsSeedMaker_LowMomentum:

Collaboration diagram for InDet::SiSpacePointsSeedMaker_LowMomentum:

Public Member Functions
Standard tool methods
	SiSpacePointsSeedMaker_LowMomentum (const std::string &, const std::string &, const IInterface *)
virtual	~SiSpacePointsSeedMaker_LowMomentum ()=default
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
Methods to initialize tool for new event or region
virtual void	newEvent (const EventContext &ctx, EventData &data, int iteration) const override
virtual void	newRegion (const EventContext &ctx, EventData &data, const std::vector< IdentifierHash > &vPixel, const std::vector< IdentifierHash > &vSCT) const override
virtual void	newRegion (const EventContext &ctx, EventData &data, const std::vector< IdentifierHash > &vPixel, const std::vector< IdentifierHash > &vSCT, const IRoiDescriptor &iRD) const override
Methods to initilize different strategies of seeds production
virtual void	find2Sp (EventData &data, const std::list< Trk::Vertex > &lv) const override
	with two space points with or without vertex constraint
virtual void	find3Sp (const EventContext &ctx, EventData &data, const std::list< Trk::Vertex > &lv) const override
	with three space points with or without vertex constraint
virtual void	find3Sp (const EventContext &ctx, EventData &data, const std::list< Trk::Vertex > &lv, const double *zVertex) const override
	with three space points with or without vertex constraint with information about min and max Z of the vertex
virtual void	findVSp (const EventContext &ctx, EventData &data, const std::list< Trk::Vertex > &lv) const override
	with variable number space points with or without vertex constraint Variable means (2,3,4,....) any number space points
Iterator through seeds pseudo collection
produced accordingly methods find
virtual const SiSpacePointsSeed *	next (const EventContext &ctx, EventData &data) const override
virtual void	writeNtuple (const SiSpacePointsSeed *seed, const Trk::Track *track, int seedType, long eventNumber) const override
virtual bool	getWriteNtupleBoolProperty () const override

Private Attributes
Data handles
SG::ReadHandleKey< SpacePointContainer >	m_spacepointsSCT {this, "SpacePointsSCTName", "SCT_SpacePoints", "SCT space points container"}
SG::ReadHandleKey< SpacePointContainer >	m_spacepointsPixel {this, "SpacePointsPixelName", "PixelSpacePoints", "Pixel space points container"}
SG::ReadHandleKey< SpacePointOverlapCollection >	m_spacepointsOverlap {this, "SpacePointsOverlapName", "OverlapSpacePoints"}
SG::ReadHandleKey< Trk::PRDtoTrackMap >	m_prdToTrackMap {this,"PRDtoTrackMap","","option PRD-to-track association"}
SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey { this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" }
SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCondObjInputKey
Properties, which will not be changed after construction
BooleanProperty	m_pixel {this, "usePixel", true}
BooleanProperty	m_sct {this, "useSCT", true}
BooleanProperty	m_useOverlap {this, "useOverlapSpCollection", false}
IntegerProperty	m_maxsize {this, "maxSize", 2000}
IntegerProperty	m_maxsizeSP {this, "maxSizeSP", 1500}
IntegerProperty	m_maxOneSize {this, "maxSeedsForSpacePoint", 5}
FloatProperty	m_r1min {this, "minRadius1", 0.}
FloatProperty	m_r1max {this, "maxRadius1", 600.}
FloatProperty	m_r2min {this, "minRadius2", 0.}
FloatProperty	m_r2max {this, "maxRadius2", 600.}
FloatProperty	m_r3min {this, "minRadius3", 0.}
FloatProperty	m_drmin {this, "mindRadius", 10.}
FloatProperty	m_drmax {this, "maxdRadius", 200.}
FloatProperty	m_zmin {this, "minZ", -250.}
FloatProperty	m_zmax {this, "maxZ", +250.}
FloatProperty	m_r_rmax {this, "radMax", 200.}
FloatProperty	m_r_rstep {this, "radStep", 2.}
FloatProperty	m_dzver {this, "maxdZver", 5.}
FloatProperty	m_dzdrver {this, "maxdZdRver", 0.02}
FloatProperty	m_diver {this, "maxdImpact", 7.}
FloatProperty	m_ptmax {this, "pTmax", 500.}
Properties, which can be updated in initialize
FloatProperty	m_etamax {this, "etaMax", 2.7}
FloatProperty	m_ptmin {this, "pTmin", 100.}
FloatProperty	m_fieldScale {this, "fieldScale", 1.}
Data members, which are updated only in initialize
bool	m_initialized {false}
int	m_outputlevel {0}
int	m_fNmax {0}
int	m_r_size {0}
int	m_rfz_b [SizeRFZ] {}
int	m_rfz_t [SizeRFZ] {}
int	m_rfz_ib [SizeRFZ][SizeI] {}
int	m_rfz_it [SizeRFZ][SizeI] {}
float	m_dzdrmin {0.}
float	m_dzdrmax {0.}
float	m_r3max {0.}
float	m_iptmin {0.}
float	m_iptmax {1.f/400.f}
float	m_sF {0.}

Print internal tool parameters and status
enum	Size { SizeRF =20 , SizeZ =11 , SizeRFZ =SizeRF*SizeZ , SizeI =9 }
	enum for array sizes More...
virtual MsgStream &	dump (EventData &data, MsgStream &out) const override

Disallow default instantiation, copy, assignment
	SiSpacePointsSeedMaker_LowMomentum ()=delete
	SiSpacePointsSeedMaker_LowMomentum (const SiSpacePointsSeedMaker_LowMomentum &)=delete
SiSpacePointsSeedMaker_LowMomentum &	operator= (const SiSpacePointsSeedMaker_LowMomentum &)=delete
MsgStream &	dumpConditions (EventData &data, MsgStream &out) const
void	buildFrameWork ()
void	buildBeamFrameWork (EventData &data) const
void	newOneSeed (EventData &data, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const float &, const float &) const
void	fillLists (EventData &data) const
void	production3Sp (const EventContext &ctx, EventData &data) const
void	production3Sp (EventData &data, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, std::vector< InDet::SiSpacePointForSeed * >::iterator *, int, int, int &, float) const
void	findNext (const EventContext &ctx, EventData &data) const
bool	isZCompatible (EventData &data, float &, float &, float &) const
bool	isUsed (const Trk::SpacePoint *, const Trk::PRDtoTrackMap &prd_to_track_map) const
void	initializeEventData (EventData &data) const
static MsgStream &	dumpEvent (EventData &data, MsgStream &out)
static SiSpacePointForSeed *	newSpacePoint (EventData &data, const Trk::SpacePoint *const &)
static void	newSeed (EventData &data, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const float &)
static void	newSeed (EventData &data, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const Trk::SpacePoint *&, const float &)
static void	fillSeeds (EventData &data)
static void	erase (EventData &data)
static void	production2Sp (EventData &data)
static bool	newVertices (EventData &data, const std::list< Trk::Vertex > &)
static void	convertToBeamFrameWork (EventData &data, const Trk::SpacePoint *const &, float *)

Detailed Description

Definition at line 52 of file SiSpacePointsSeedMaker_LowMomentum.h.

Member Enumeration Documentation

Size

enum InDet::SiSpacePointsSeedMaker_LowMomentum::Size

private

enum for array sizes

Enumerator
SizeRF
SizeZ
SizeRFZ
SizeI

Definition at line 125 of file SiSpacePointsSeedMaker_LowMomentum.h.

              {SizeRF=20,
               SizeZ=11,
               SizeRFZ=SizeRF*SizeZ,
               SizeI=9};

Constructor & Destructor Documentation

SiSpacePointsSeedMaker_LowMomentum() [1/3]

InDet::SiSpacePointsSeedMaker_LowMomentum::SiSpacePointsSeedMaker_LowMomentum	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p )

Definition at line 27 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

  : base_class(t, n, p)
{
}

~SiSpacePointsSeedMaker_LowMomentum()

virtual InDet::SiSpacePointsSeedMaker_LowMomentum::~SiSpacePointsSeedMaker_LowMomentum ( )

virtualdefault

SiSpacePointsSeedMaker_LowMomentum() [2/3]

InDet::SiSpacePointsSeedMaker_LowMomentum::SiSpacePointsSeedMaker_LowMomentum ( )

privatedelete

SiSpacePointsSeedMaker_LowMomentum() [3/3]

InDet::SiSpacePointsSeedMaker_LowMomentum::SiSpacePointsSeedMaker_LowMomentum ( const SiSpacePointsSeedMaker_LowMomentum & )

privatedelete

Member Function Documentation

buildBeamFrameWork()

void InDet::SiSpacePointsSeedMaker_LowMomentum::buildBeamFrameWork ( EventData & data ) const

private

Definition at line 660 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{ 
  SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey };
 
  const Amg::Vector3D &cb =     beamSpotHandle->beamPos();
  double     tx = std::tan(beamSpotHandle->beamTilt(0));
  double     ty = std::tan(beamSpotHandle->beamTilt(1));
 
  double ph   = std::atan2(ty,tx);
  double th   = std::acos(1.f/std::sqrt(1.f+tx*tx+ty*ty));
  double sint = std::sin(th);
  double cost = std::cos(th);
  double sinp = std::sin(ph);
  double cosp = std::cos(ph);
  
  data.xbeam[0] =  static_cast<float>(cb.x());
  data.xbeam[1] =  static_cast<float>(cost*cosp*cosp+sinp*sinp);
  data.xbeam[2] =  static_cast<float>(cost*sinp*cosp-sinp*cosp);
  data.xbeam[3] = -static_cast<float>(sint*cosp               );
  
  data.ybeam[0] =  static_cast<float>(cb.y());
  data.ybeam[1] =  static_cast<float>(cost*cosp*sinp-sinp*cosp);
  data.ybeam[2] =  static_cast<float>(cost*sinp*sinp+cosp*cosp);
  data.ybeam[3] = -static_cast<float>(sint*sinp               );
  
  data.zbeam[0] =  static_cast<float>(cb.z());
  data.zbeam[1] =  static_cast<float>(sint*cosp);
  data.zbeam[2] =  static_cast<float>(sint*sinp);
  data.zbeam[3] =  static_cast<float>(cost);
}

buildFrameWork()

void InDet::SiSpacePointsSeedMaker_LowMomentum::buildFrameWork ( )

private

Definition at line 568 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  m_ptmin = std::max( std::abs(m_ptmin), float(50.*m_fieldScale));
  m_iptmax  = 1.f/std::abs(m_ptmax);
  m_iptmin  = 1.f/std::abs(m_ptmin);
  m_etamax    = std::abs(m_etamax);
  m_dzdrmax   = 1.f/std::tan(2.f*std::atan(std::exp(-m_etamax)));
  m_dzdrmin   =-m_dzdrmax;
  m_r3max     = m_r_rmax;
 
  // Build radius sorted containers
  //
  m_r_size = static_cast<int>((m_r_rmax+.1f)/m_r_rstep);
 
  // Build radius-azimuthal sorted containers
  //
  constexpr float pi2 = 2.*M_PI;
  const int   NFmax = SizeRF;
  const float sFmax   = static_cast<float>(NFmax )/pi2;
  const float sFmin   = 100./60.;
  m_sF = m_ptmin/m_fieldScale /60.f;
  if (m_sF > sFmax ) m_sF = sFmax;
  else if (m_sF < sFmin) m_sF = sFmin;
  m_fNmax = static_cast<int>(pi2*m_sF);
  if (m_fNmax >=NFmax) m_fNmax = NFmax-1;
 
  // Build maps for radius-azimuthal-Z sorted collections 
  //
  for (int f=0; f<=m_fNmax; ++f) {
    int fb = f-1; if (fb<0      ) fb=m_fNmax;
    int ft = f+1; if (ft>m_fNmax) ft=0;
    
    // For each azimuthal region loop through all Z regions
    //
    for (int z=0; z<SizeZ; ++z) {
      int a        = f *SizeZ+z;
      int b        = fb*SizeZ+z;
      int c        = ft*SizeZ+z;
      m_rfz_b [a]    = 3; m_rfz_t [a]    = 3;
      m_rfz_ib[a][0] = a; m_rfz_it[a][0] = a;
      m_rfz_ib[a][1] = b; m_rfz_it[a][1] = b;
      m_rfz_ib[a][2] = c; m_rfz_it[a][2] = c;
      if (z==5) {
        m_rfz_t [a]    = 9;
        m_rfz_it[a][3] = a+1;
        m_rfz_it[a][4] = b+1;
        m_rfz_it[a][5] = c+1;
        m_rfz_it[a][6] = a-1;
        m_rfz_it[a][7] = b-1;
        m_rfz_it[a][8] = c-1;
      } else if (z> 5) {
        m_rfz_b [a]    = 6;
        m_rfz_ib[a][3] = a-1;
        m_rfz_ib[a][4] = b-1;
        m_rfz_ib[a][5] = c-1;
        if (z<10) {
          m_rfz_t [a]    = 6;
          m_rfz_it[a][3] = a+1;
          m_rfz_it[a][4] = b+1;
          m_rfz_it[a][5] = c+1;
        }
      } else {
        m_rfz_b [a]    = 6;
        m_rfz_ib[a][3] = a+1;
        m_rfz_ib[a][4] = b+1;
        m_rfz_ib[a][5] = c+1;
        if (z>0) {
          m_rfz_t [a]    = 6;
          m_rfz_it[a][3] = a-1;
          m_rfz_it[a][4] = b-1;
          m_rfz_it[a][5] = c-1;
        }
      }
      if (z==3) {
        m_rfz_b[a]      = 9;
        m_rfz_ib[a][6] = a+2;
        m_rfz_ib[a][7] = b+2;
        m_rfz_ib[a][8] = c+2;
      } else if (z==7) {
        m_rfz_b[a]      = 9;
        m_rfz_ib[a][6] = a-2;
        m_rfz_ib[a][7] = b-2;
        m_rfz_ib[a][8] = c-2;
      }
    }
  }
}

convertToBeamFrameWork()

void InDet::SiSpacePointsSeedMaker_LowMomentum::convertToBeamFrameWork ( EventData & data, const Trk::SpacePoint *const & sp, float * r )

staticprivate

Definition at line 695 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  r[0] = static_cast<float>(sp->globalPosition().x())-data.xbeam[0];
  r[1] = static_cast<float>(sp->globalPosition().y())-data.ybeam[0];
  r[2] = static_cast<float>(sp->globalPosition().z())-data.zbeam[0];
}

dump()

MsgStream & InDet::SiSpacePointsSeedMaker_LowMomentum::dump ( EventData & data, MsgStream & out ) const

overridevirtual

Definition at line 369 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  if (data.nprint) return dumpEvent(data, out);
  return dumpConditions(data, out);
}

dumpConditions()

MsgStream & InDet::SiSpacePointsSeedMaker_LowMomentum::dumpConditions ( EventData & data, MsgStream & out ) const

private

Definition at line 381 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  int n = 42-m_spacepointsPixel.key().size();
  std::string s2; for (int i=0; i<n; ++i) s2.append(" "); s2.append("|");
  n     = 42-m_spacepointsSCT.key().size();
  std::string s3; for (int i=0; i<n; ++i) s3.append(" "); s3.append("|");
  n     = 42-m_spacepointsOverlap.key().size();
  std::string s4; for (int i=0; i<n; ++i) s4.append(" "); s4.append("|");
  n     = 42-m_beamSpotKey.key().size();
  std::string s5; for (int i=0; i<n; ++i) s5.append(" "); s5.append("|");
 
  out<<"|---------------------------------------------------------------------|"
     <<endmsg;
  out<<"| Pixel    space points   | "<<m_spacepointsPixel.key() <<s2
     <<endmsg;
  out<<"| SCT      space points   | "<<m_spacepointsSCT.key()<<s3
     <<endmsg;
  out<<"| Overlap  space points   | "<<m_spacepointsOverlap.key()<<s4
     <<endmsg;
  out<<"| BeamConditionsService   | "<<m_beamSpotKey.key()<<s5
     <<endmsg;
  out<<"| usePixel                | "
     <<std::setw(12)<<m_pixel 
     <<"                              |"<<endmsg;
  out<<"| useSCT                  | "
     <<std::setw(12)<<m_sct 
     <<"                              |"<<endmsg;
  out<<"| Use PRD-to-track assoc.?| "
     <<std::setw(12)<< (!m_prdToTrackMap.key().empty() ? "yes" : "no ")
     <<"                              |"<<endmsg;
  out<<"| maxSize                 | "
     <<std::setw(12)<<m_maxsize 
     <<"                              |"<<endmsg;
  out<<"| maxSizeSP               | "
     <<std::setw(12)<<m_maxsizeSP
     <<"                              |"<<endmsg;
  out<<"| pTmin  (mev)            | "
     <<std::setw(12)<<std::setprecision(5)<<m_ptmin
     <<"                              |"<<endmsg;
  out<<"| pTmax  (mev)            | "
     <<std::setw(12)<<std::setprecision(5)<<m_ptmax
     <<"                              |"<<endmsg;
  out<<"| |eta|               <=  | " 
     <<std::setw(12)<<std::setprecision(5)<<m_etamax
     <<"                              |"<<endmsg;
  out<<"| max radius SP           | "
     <<std::setw(12)<<std::setprecision(5)<<m_r_rmax 
     <<"                              |"<<endmsg;
  out<<"| radius step             | "
     <<std::setw(12)<<std::setprecision(5)<<m_r_rstep
     <<"                              |"<<endmsg;
  out<<"| min Z-vertex position   | "
     <<std::setw(12)<<std::setprecision(5)<<m_zmin
     <<"                              |"<<endmsg;
  out<<"| max Z-vertex position   | "
     <<std::setw(12)<<std::setprecision(5)<<m_zmax
     <<"                              |"<<endmsg;
  out<<"| min radius first  SP(3) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r1min
     <<"                              |"<<endmsg;
  out<<"| min radius second SP(3) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r2min
     <<"                              |"<<endmsg;
  out<<"| min radius last   SP(3) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r3min
     <<"                              |"<<endmsg;
  out<<"| max radius first  SP(3) | "
     <<std::setw(12)<<std::setprecision(4)<<m_r1max
     <<"                              |"<<endmsg;
  out<<"| max radius second SP(3) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r2max
     <<"                              |"<<endmsg;
  out<<"| max radius last   SP(3) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r3max
     <<"                              |"<<endmsg;
  out<<"| min space points dR     | "
     <<std::setw(12)<<std::setprecision(5)<<m_drmin
     <<"                              |"<<endmsg;
  out<<"| max space points dR     | "
     <<std::setw(12)<<std::setprecision(5)<<m_drmax
     <<"                              |"<<endmsg;
  out<<"| max dZ    impact        | "
     <<std::setw(12)<<std::setprecision(5)<<m_dzver 
     <<"                              |"<<endmsg;
  out<<"| max dZ/dR impact        | "
     <<std::setw(12)<<std::setprecision(5)<<m_dzdrver 
     <<"                              |"<<endmsg;
  out<<"| max       impact        | "
     <<std::setw(12)<<std::setprecision(5)<<m_diver
     <<"                              |"<<endmsg;
  out<<"|---------------------------------------------------------------------|"
     <<endmsg;
  out<<"| Beam X center           | "
     <<std::setw(12)<<std::setprecision(5)<<data.xbeam[0]
     <<"                              |"<<endmsg;
  out<<"| Beam Y center           | "
     <<std::setw(12)<<std::setprecision(5)<<data.ybeam[0]
     <<"                              |"<<endmsg;
  out<<"| Beam Z center           | "
     <<std::setw(12)<<std::setprecision(5)<<data.zbeam[0]
     <<"                              |"<<endmsg;
  out<<"| Beam X-axis direction   | "
     <<std::setw(12)<<std::setprecision(5)<<data.xbeam[1]
     <<std::setw(12)<<std::setprecision(5)<<data.xbeam[2]
     <<std::setw(12)<<std::setprecision(5)<<data.xbeam[3]
     <<"      |"<<endmsg;
  out<<"| Beam Y-axis direction   | "
     <<std::setw(12)<<std::setprecision(5)<<data.ybeam[1]
     <<std::setw(12)<<std::setprecision(5)<<data.ybeam[2]
     <<std::setw(12)<<std::setprecision(5)<<data.ybeam[3]
     <<"      |"<<endmsg;
  out<<"| Beam Z-axis direction   | "
     <<std::setw(12)<<std::setprecision(5)<<data.zbeam[1]
     <<std::setw(12)<<std::setprecision(5)<<data.zbeam[2]
     <<std::setw(12)<<std::setprecision(5)<<data.zbeam[3]
     <<"      |"<<endmsg;
  out<<"|---------------------------------------------------------------------|"
     <<endmsg;
  return out;
}

dumpEvent()

MsgStream & InDet::SiSpacePointsSeedMaker_LowMomentum::dumpEvent ( EventData & data, MsgStream & out )

staticprivate

Definition at line 506 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  out<<"|---------------------------------------------------------------------|"
     <<endmsg;
  out<<"| data.ns                    | "
     <<std::setw(12)<<data.ns
     <<"                              |"<<endmsg;
  out<<"| data.nsaz                  | "
     <<std::setw(12)<<data.nsaz
     <<"                              |"<<endmsg;
  out<<"| seeds                   | "
     <<std::setw(12)<<data.l_seeds.size()
     <<"                              |"<<endmsg;
  out<<"|---------------------------------------------------------------------|"
     <<endmsg;
  return out;
}

erase()

void InDet::SiSpacePointsSeedMaker_LowMomentum::erase ( EventData & data )

staticprivate

Definition at line 752 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  for (int i=0; i!=data.nr; ++i) {
    int n = data.r_index[i];
    data.r_map[n] = 0;
    data.r_Sorted[n].clear();
  }
 
  for (int i=0; i!=data.nrfz; ++i) {
    int n = data.rfz_index[i];
    data.rfz_map[n] = 0;
    data.rfz_Sorted[n].clear();
  }
 
  data.state = 0;
  data.ns    = 0;
  data.nsaz  = 0;
  data.nr    = 0;
  data.nrfz  = 0;
}

fillLists()

void InDet::SiSpacePointsSeedMaker_LowMomentum::fillLists ( EventData & data ) const

private

Definition at line 707 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  constexpr float pi2 = 2.*M_PI;
  
  for (int i=0; i<m_r_size;  ++i) {
 
    if (!data.r_map[i]) continue;
 
    for (InDet::SiSpacePointForSeed* r  : data.r_Sorted[i]) {
      
      // Azimuthal angle sort
      //
      float F = r->phi();
      if (F<0.) F+=pi2;
 
      int   f = static_cast<int>(F*m_sF);
      if (f < 0) f = m_fNmax;
      else if (f > m_fNmax) f = 0;
 
      int z;
      float Z = r->z();
 
      // Azimuthal angle and Z-coordinate sort
      //
      if (Z>0.) {
        Z< 250.?z=5:Z< 450.?z=6:Z< 925.?z=7:Z< 1400.?z=8:Z< 2500.?z=9:z=10;
      } else {
        Z>-250.?z=5:Z>-450.?z=4:Z>-925.?z=3:Z>-1400.?z=2:Z>-2500.?z=1:z= 0;
      }
      int n = f*SizeZ+z;
      ++data.nsaz;
      data.rfz_Sorted[n].push_back(r);
      if (!data.rfz_map[n]++) data.rfz_index[data.nrfz++] = n;
    }
    data.r_Sorted[i].clear();
    data.r_map[i] = 0;
  }
  data.nr    = 0;
  data.state = 0;
}

fillSeeds()

void InDet::SiSpacePointsSeedMaker_LowMomentum::fillSeeds ( EventData & data )

staticprivate

Definition at line 1158 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  std::multimap<float,InDet::SiSpacePointsSeed*>::iterator 
    l  = data.mapOneSeeds.begin(),
    le = data.mapOneSeeds.end  ();
 
  for (; l!=le; ++l) {
    if (data.i_seede!=data.l_seeds.end()) {
      InDet::SiSpacePointsSeed* s = &(*data.i_seede++);
      *s = *(*l).second;
    } else {
      data.l_seeds.emplace_back(*(*l).second);
      data.i_seede = data.l_seeds.end();
    }
  }
}

finalize()

StatusCode InDet::SiSpacePointsSeedMaker_LowMomentum::finalize ( )

overridevirtual

Definition at line 77 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  return AlgTool::finalize();
}

find2Sp()

void InDet::SiSpacePointsSeedMaker_LowMomentum::find2Sp ( EventData & data, const std::list< Trk::Vertex > & lv ) const

overridevirtual

with two space points with or without vertex constraint

Definition at line 268 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  int mode = 0;
  if (lv.begin()!=lv.end()) mode = 1;
  bool newv = newVertices(data, lv);
  
  if (newv || !data.state || data.nspoint!=2 || data.mode!=mode || data.nlist) {
    data.i_seede = data.l_seeds.begin();
    data.state   = 1;
    data.nspoint = 2;
    data.nlist   = 0;
    data.mode    = mode;
    data.endlist = true;
    data.fNmin   = 0;
    data.zMin    = 0;
    production2Sp(data);
  }
  data.i_seed = data.l_seeds.begin();
  
  if (m_outputlevel<=0) {
    data.nprint=1;
    dump(data, msg(MSG::DEBUG));
  }
}

find3Sp() [1/2]

void InDet::SiSpacePointsSeedMaker_LowMomentum::find3Sp ( const EventContext & ctx, EventData & data, const std::list< Trk::Vertex > & lv ) const

overridevirtual

with three space points with or without vertex constraint

Definition at line 300 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  int mode = 2;
  if (lv.begin()!=lv.end()) mode = 3;
  bool newv = newVertices(data, lv);
 
  if (newv || !data.state || data.nspoint!=3 || data.mode!=mode || data.nlist) {
    data.i_seede = data.l_seeds.begin();
    data.state   = 1;
    data.nspoint = 3;
    data.nlist   = 0;
    data.mode    = mode;
    data.endlist = true;
    data.fNmin   = 0;
    data.zMin    = 0;
    production3Sp(ctx, data);
  }
  data.i_seed = data.l_seeds.begin();
 
  if (m_outputlevel<=0) {
    data.nprint=1;
    dump(data, msg(MSG::DEBUG));
  }
}

find3Sp() [2/2]

void InDet::SiSpacePointsSeedMaker_LowMomentum::find3Sp ( const EventContext & ctx, EventData & data, const std::list< Trk::Vertex > & lv, const double * zVertex ) const

overridevirtual

with three space points with or without vertex constraint with information about min and max Z of the vertex

Definition at line 327 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  find3Sp(ctx, data, lv);
}

findNext()

void InDet::SiSpacePointsSeedMaker_LowMomentum::findNext ( const EventContext & ctx, EventData & data ) const

private

Definition at line 528 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (data.endlist) return;
  
  data.i_seede = data.l_seeds.begin();
  if (data.mode==0 || data.mode==1) {
    production2Sp(data);
  } else if (data.mode==2 || data.mode==3) {
    production3Sp(ctx, data);
  } else if (data.mode==5 || data.mode==6) {
    production3Sp(ctx, data);
  }
 
  data.i_seed = data.l_seeds.begin();
  ++data.nlist;
}                       

findVSp()

void InDet::SiSpacePointsSeedMaker_LowMomentum::findVSp ( const EventContext & ctx, EventData & data, const std::list< Trk::Vertex > & lv ) const

overridevirtual

with variable number space points with or without vertex constraint Variable means (2,3,4,....) any number space points

Definition at line 338 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  int mode = 5;
  if (lv.begin()!=lv.end()) mode = 6;
  bool newv = newVertices(data, lv);
  
  if (newv || !data.state || data.nspoint!=4 || data.mode!=mode || data.nlist) {
    data.i_seede = data.l_seeds.begin();
    data.state   = 1;
    data.nspoint = 4;
    data.nlist   = 0;
    data.mode    = mode;
    data.endlist = true;
    data.fNmin   = 0;
    data.zMin    = 0;
    production3Sp(ctx, data);
  }
  data.i_seed = data.l_seeds.begin();
 
  if (m_outputlevel<=0) {
    data.nprint=1;
    dump(data, msg(MSG::DEBUG));
  }
}

getWriteNtupleBoolProperty()

bool InDet::SiSpacePointsSeedMaker_LowMomentum::getWriteNtupleBoolProperty ( ) const

overridevirtual

Definition at line 1190 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

                                                                              {
    return false;
}

initialize()

StatusCode InDet::SiSpacePointsSeedMaker_LowMomentum::initialize ( )

overridevirtual

Definition at line 37 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  StatusCode sc = AlgTool::initialize();
 
  ATH_CHECK(m_spacepointsPixel.initialize(m_pixel));
  ATH_CHECK(m_spacepointsSCT.initialize(m_sct));
  ATH_CHECK(m_spacepointsOverlap.initialize(m_useOverlap));
 
  // Get beam geometry
  //
  ATH_CHECK(m_beamSpotKey.initialize());
 
  ATH_CHECK( m_fieldCondObjInputKey.initialize());
 
  // PRD-to-track association (optional)
  ATH_CHECK( m_prdToTrackMap.initialize( !m_prdToTrackMap.key().empty()));
 
  // Build framework
  //
  buildFrameWork();
 
  // Get output print level
  //
  m_outputlevel = msg().level()-MSG::DEBUG;
  if (m_outputlevel<=0) {
    EventData data;
    initializeEventData(data);
    data.nprint=0;
    dump(data, msg(MSG::DEBUG));
  }
 
  m_initialized = true;
 
  return sc;
}

initializeEventData()

void InDet::SiSpacePointsSeedMaker_LowMomentum::initializeEventData ( EventData & data ) const

private

Definition at line 1175 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

                                                                                       {
  data.initialize(EventData::ToolType::LowMomentum,
                  m_maxsizeSP,
                  m_maxOneSize,
                  0, // maxsize not used
                  m_r_size,
                  0, // sizeRF not used
                  SizeRFZ,
                  0, // sizeRFZV not used
                  false); // checkEta not used
}

isUsed()

bool InDet::SiSpacePointsSeedMaker_LowMomentum::isUsed ( const Trk::SpacePoint * sp, const Trk::PRDtoTrackMap & prd_to_track_map ) const

inlineprivate

Definition at line 255 of file SiSpacePointsSeedMaker_LowMomentum.h.

  {
    const Trk::PrepRawData* d = sp->clusterList().first;
    if (!d || !prd_to_track_map.isUsed(*d)) return false;
 
    d = sp->clusterList().second;
    if (!d || prd_to_track_map.isUsed(*d)) return true;
 
    return false;
  }

isZCompatible()

bool InDet::SiSpacePointsSeedMaker_LowMomentum::isZCompatible ( EventData & data, float & Zv, float & R, float & T ) const

private

Definition at line 1072 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (Zv < m_zmin || Zv > m_zmax) return false;
 
  if (data.l_vertex.empty()) return true;
 
  float dZmin = std::numeric_limits<float>::max();
  for (const float& v : data.l_vertex) {
    float dZ = std::abs(v-Zv);
    if (dZ<dZmin) dZmin=dZ;
  }
  return dZmin < (m_dzver+m_dzdrver*R)*std::sqrt(1.f+T*T);
}

newEvent()

void InDet::SiSpacePointsSeedMaker_LowMomentum::newEvent ( const EventContext & ctx, EventData & data, int iteration ) const

overridevirtual

Definition at line 86 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  data.trigger = false;
  if (!m_pixel && !m_sct) return;
  erase(data);
  data.i_spforseed = data.l_spforseed.begin();
  buildBeamFrameWork(data);
 
  float irstep = 1./m_r_rstep;
  int   irmax  = m_r_size-1;
 
  SG::ReadHandle<Trk::PRDtoTrackMap>  prd_to_track_map;
  const Trk::PRDtoTrackMap *prd_to_track_map_cptr = nullptr;
  if (!m_prdToTrackMap.key().empty()) {
    prd_to_track_map=SG::ReadHandle<Trk::PRDtoTrackMap>(m_prdToTrackMap, ctx);
    if (!prd_to_track_map.isValid()) {
      ATH_MSG_ERROR("Failed to read PRD to track association map: " << m_prdToTrackMap.key());
    }
    prd_to_track_map_cptr = prd_to_track_map.cptr();
  }
 
  // Get pixels space points containers from store gate 
  //
  if (m_pixel) {
 
    SG::ReadHandle<SpacePointContainer> spacepointsPixel{m_spacepointsPixel, ctx};
    if (spacepointsPixel.isValid()) {
 
      for (const SpacePointCollection* spc: *spacepointsPixel) {
        for (const Trk::SpacePoint* sp: *spc) {
 
          float r = sp->r();
          if (r<0. || r>=m_r_rmax) continue;
          if (prd_to_track_map_cptr && isUsed(sp,*prd_to_track_map_cptr)) continue;
 
          InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
 
          int ir = static_cast<int>(sps->radius()*irstep);
          if (ir>irmax) ir = irmax;
          data.r_Sorted[ir].push_back(sps);
          ++data.r_map[ir];
          if (data.r_map[ir]==1) data.r_index[data.nr++] = ir;
          ++data.ns;
        }
      }
    }
  }
 
  // Get sct space points containers from store gate 
  //
  if (m_sct) {
 
    SG::ReadHandle<SpacePointContainer> spacepointsSCT{m_spacepointsSCT, ctx};
    if (spacepointsSCT.isValid()) {
 
      for (const SpacePointCollection* spc: *spacepointsSCT) {
        for (const Trk::SpacePoint* sp: *spc) {
 
          float r = sp->r();
          if (r<0. || r>=m_r_rmax) continue;
          if (prd_to_track_map_cptr && isUsed(sp,*prd_to_track_map_cptr)) continue;
 
          InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
 
          int ir = static_cast<int>(sps->radius()*irstep);
          if (ir>irmax) ir = irmax;
          data.r_Sorted[ir].push_back(sps);
          ++data.r_map[ir];
          if (data.r_map[ir]==1) data.r_index[data.nr++] = ir;
          ++data.ns;
        }
      }
    }
  }
  fillLists(data);
}

newOneSeed()

void InDet::SiSpacePointsSeedMaker_LowMomentum::newOneSeed ( EventData & data, const Trk::SpacePoint *& p1, const Trk::SpacePoint *& p2, const Trk::SpacePoint *& p3, const float & z, const float & q ) const

private

Definition at line 1025 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (data.nOneSeeds < m_maxOneSize) {
    data.OneSeeds[data.nOneSeeds].erase();
    data.OneSeeds[data.nOneSeeds].add(p1);
    data.OneSeeds[data.nOneSeeds].add(p2);
    data.OneSeeds[data.nOneSeeds].add(p3);
    data.OneSeeds[data.nOneSeeds].setZVertex(static_cast<double>(z));
    data.mapOneSeeds.insert(std::make_pair(q, &(data.OneSeeds[data.nOneSeeds])));
    ++data.nOneSeeds;
  } else {
    std::multimap<float,InDet::SiSpacePointsSeed*>::reverse_iterator 
      l = data.mapOneSeeds.rbegin();
    if ((*l).first <= q) return;
    
    InDet::SiSpacePointsSeed* s = (*l).second;
    s->erase     (  );
    s->add       (p1);
    s->add       (p2);
    s->add       (p3);
    s->setZVertex(static_cast<double>(z));
    std::multimap<float,InDet::SiSpacePointsSeed*>::iterator 
      i = data.mapOneSeeds.insert(std::make_pair(q,s));
    for (++i; i!=data.mapOneSeeds.end(); ++i) {
      if ((*i).second==s) {
        data.mapOneSeeds.erase(i);
        return;
      }
    }
  }
}

newRegion() [1/2]

void InDet::SiSpacePointsSeedMaker_LowMomentum::newRegion ( const EventContext & ctx, EventData & data, const std::vector< IdentifierHash > & vPixel, const std::vector< IdentifierHash > & vSCT ) const

overridevirtual

Definition at line 169 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

                                                                                      {
  if (not data.initialized) initializeEventData(data);
 
  data.trigger = false;
  if (!m_pixel && !m_sct) return;
  erase(data);
  data.i_spforseed = data.l_spforseed.begin();
  buildBeamFrameWork(data);
 
  int   irmax  = m_r_size-1;
  float irstep = 1.f/m_r_rstep;
 
  SG::ReadHandle<Trk::PRDtoTrackMap>  prd_to_track_map;
  const Trk::PRDtoTrackMap  *prd_to_track_map_cptr = nullptr;
  if (!m_prdToTrackMap.key().empty()) {
    prd_to_track_map=SG::ReadHandle<Trk::PRDtoTrackMap>(m_prdToTrackMap, ctx);
    if (!prd_to_track_map.isValid()) {
      ATH_MSG_ERROR("Failed to read PRD to track association map: " << m_prdToTrackMap.key());
    }
    prd_to_track_map_cptr = prd_to_track_map.cptr();
  }
 
  // Get pixels space points containers from store gate 
  //
  if (m_pixel && !vPixel.empty()) {
 
    SG::ReadHandle<SpacePointContainer> spacepointsPixel{m_spacepointsPixel, ctx};
    if (spacepointsPixel.isValid()) {
      // Loop through all trigger collections
      //
      for (const IdentifierHash& l: vPixel) {
        const auto *w = spacepointsPixel->indexFindPtr(l);
        if (w==nullptr) continue;
        for (const Trk::SpacePoint* sp: *w) {
          float r = sp->r();
          if (r<0. || r>=m_r_rmax) continue;
          if (prd_to_track_map_cptr && isUsed(sp,*prd_to_track_map_cptr)) continue;
 
          InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
 
          int ir = static_cast<int>(sps->radius()*irstep);
          if (ir>irmax) ir = irmax;
          data.r_Sorted[ir].push_back(sps);
          ++data.r_map[ir];
          if (data.r_map[ir]==1) data.r_index[data.nr++] = ir;
          ++data.ns;
        }
      }
    }
  }
 
  // Get sct space points containers from store gate 
  //
  if (m_sct && !vSCT.empty()) {
 
    SG::ReadHandle<SpacePointContainer> spacepointsSCT{m_spacepointsSCT, ctx};
    if (spacepointsSCT.isValid()) {
      // Loop through all trigger collections
      //
      for (const IdentifierHash& l: vSCT) {
        const auto *w = spacepointsSCT->indexFindPtr(l);
        if (w==nullptr) continue;
        for (const Trk::SpacePoint* sp: *w) {
          float r = sp->r();
          if (r<0. || r>=m_r_rmax) continue;
          if (prd_to_track_map_cptr && isUsed(sp,*prd_to_track_map_cptr)) continue;
          InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
          int ir = static_cast<int>(sps->radius()*irstep);
          if (ir>irmax) ir = irmax;
          data.r_Sorted[ir].push_back(sps);
          ++data.r_map[ir];
          if (data.r_map[ir]==1) data.r_index[data.nr++] = ir;
          ++data.ns;
        }
      }
    }
  }
  fillLists(data);
}

newRegion() [2/2]

void InDet::SiSpacePointsSeedMaker_LowMomentum::newRegion ( const EventContext & ctx, EventData & data, const std::vector< IdentifierHash > & vPixel, const std::vector< IdentifierHash > & vSCT, const IRoiDescriptor & iRD ) const

overridevirtual

Definition at line 255 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  newRegion(ctx, data, vPixel, vSCT);
}

newSeed() [1/2]

void InDet::SiSpacePointsSeedMaker_LowMomentum::newSeed ( EventData & data, const Trk::SpacePoint *& p1, const Trk::SpacePoint *& p2, const float & z )

staticprivate

Definition at line 1115 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (data.i_seede!=data.l_seeds.end()) {
    InDet::SiSpacePointsSeed* s = &(*data.i_seede++);
    s->erase();
    s->add(p1);
    s->add(p2);
    s->setZVertex(static_cast<double>(z));
  } else {
    data.l_seeds.emplace_back(p1, p2, z);
    data.i_seede = data.l_seeds.end();
  }
}

newSeed() [2/2]

void InDet::SiSpacePointsSeedMaker_LowMomentum::newSeed ( EventData & data, const Trk::SpacePoint *& p1, const Trk::SpacePoint *& p2, const Trk::SpacePoint *& p3, const float & z )

staticprivate

Definition at line 1136 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (data.i_seede!=data.l_seeds.end()) {
    InDet::SiSpacePointsSeed* s = &(*data.i_seede++);
    s->erase();
    s->add(p1);
    s->add(p2);
    s->add(p3);
    s->setZVertex(static_cast<double>(z));
  } else {
    data.l_seeds.emplace_back(p1, p2, p3, z);
    data.i_seede = data.l_seeds.end();
  }
}

newSpacePoint()

InDet::SiSpacePointForSeed * InDet::SiSpacePointsSeedMaker_LowMomentum::newSpacePoint ( EventData & data, const Trk::SpacePoint *const & sp )

staticprivate

Definition at line 1091 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  InDet::SiSpacePointForSeed* sps = nullptr;
 
  float r[3];
  convertToBeamFrameWork(data, sp, r);
 
  if (data.i_spforseed!=data.l_spforseed.end()) {
    sps = &(*data.i_spforseed++);
    sps->set(sp, r);
  } else {
    data.l_spforseed.emplace_back(sp, r);
    sps = &(data.l_spforseed.back());
    data.i_spforseed = data.l_spforseed.end();
  }
      
  return sps;
}

newVertices()

bool InDet::SiSpacePointsSeedMaker_LowMomentum::newVertices ( EventData & data, const std::list< Trk::Vertex > & lV )

staticprivate

Definition at line 549 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  unsigned int s1 = data.l_vertex.size();
  unsigned int s2 = lV.size();
 
  if (s1==0 && s2==0) return false;
 
  data.l_vertex.clear();
  
  for (const Trk::Vertex& v : lV) {
    data.l_vertex.insert(static_cast<float>(v.position().z()));
  }
  return false;
}

next()

const InDet::SiSpacePointsSeed * InDet::SiSpacePointsSeedMaker_LowMomentum::next ( const EventContext & ctx, EventData & data ) const

overridevirtual

Definition at line 1060 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  if (data.i_seed==data.i_seede) {
    findNext(ctx, data);
    //cppcheck-suppress identicalInnerCondition
    if (data.i_seed==data.i_seede) return nullptr;
  } 
  return &(*data.i_seed++);
}

operator=()

SiSpacePointsSeedMaker_LowMomentum & InDet::SiSpacePointsSeedMaker_LowMomentum::operator= ( const SiSpacePointsSeedMaker_LowMomentum & )

privatedelete

production2Sp()

void InDet::SiSpacePointsSeedMaker_LowMomentum::production2Sp ( EventData & data )

staticprivate

Definition at line 777 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  data.endlist = true;
}

production3Sp() [1/2]

void InDet::SiSpacePointsSeedMaker_LowMomentum::production3Sp ( const EventContext & ctx, EventData & data ) const

private

Definition at line 786 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
  if (data.nsaz<3) return;
 
  float K = 0.;
  double f[3], gP[3] ={10.,10.,0.};
 
  MagField::AtlasFieldCache    fieldCache;
 
  // Get field cache object
  SG::ReadCondHandle<AtlasFieldCacheCondObj> readHandle{m_fieldCondObjInputKey, ctx};
  const AtlasFieldCacheCondObj* fieldCondObj{*readHandle};
  if (fieldCondObj == nullptr) {
    ATH_MSG_ERROR("SiSpacePointsSeedMaker_LowMomentum: Failed to retrieve AtlasFieldCacheCondObj with key " << m_fieldCondObjInputKey.key());
    return;
  }
  fieldCondObj->getInitializedCache (fieldCache);
 
  if (fieldCache.solenoidOn()) {
    fieldCache.getFieldZR(gP, f);
 
    K = 2./(300.*f[2]);
  } else {
    K = 2./(300.* 5. );
  }
 
  const int   ZI[SizeZ]= {5,6,7,8,9,10,4,3,2,1,0};
  std::vector<InDet::SiSpacePointForSeed*>::iterator rt[9],rte[9],rb[9],rbe[9];
  int nseed = 0;
 
  // Loop thorugh all azimuthal regions
  //
  for (int f=data.fNmin; f<=m_fNmax; ++f) {
    // For each azimuthal region loop through all Z regions
    //
    int z = 0;
    if (!data.endlist) z = data.zMin;
 
    for (; z<SizeZ; ++z) {
      int a  = f *SizeZ+ZI[z];
      if (!data.rfz_map[a]) continue;
      int NB = 0, NT = 0;
      for (int i=0; i!=m_rfz_b[a]; ++i) {
        int an =  m_rfz_ib[a][i];
        if (!data.rfz_map[an]) continue;
        rb [NB] = data.rfz_Sorted[an].begin();
        rbe[NB++] = data.rfz_Sorted[an].end();
      } 
      for (int i=0; i!=m_rfz_t[a]; ++i) {
        int an =  m_rfz_it[a][i];
        if (!data.rfz_map[an]) continue;
        rt [NT] = data.rfz_Sorted[an].begin();
        rte[NT++] = data.rfz_Sorted[an].end();
      } 
      production3Sp(data, rb, rbe, rt, rte, NB, NT, nseed, K);
      if (!data.endlist) {
        data.fNmin=f;
        data.zMin = z;
        return;
      }
    }
  }
  data.endlist = true;
}

production3Sp() [2/2]

void InDet::SiSpacePointsSeedMaker_LowMomentum::production3Sp ( EventData & data, std::vector< InDet::SiSpacePointForSeed * >::iterator * rb, std::vector< InDet::SiSpacePointForSeed * >::iterator * rbe, std::vector< InDet::SiSpacePointForSeed * >::iterator * rt, std::vector< InDet::SiSpacePointForSeed * >::iterator * rte, int NB, int NT, int & nseed, float K ) const

private

Definition at line 855 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

{
 
  const float COF  = 134*.07*9.;
  const float COFP = 134*.2*9.;
 
  std::vector<InDet::SiSpacePointForSeed*>::iterator r0=rb[0],r;
  if (!data.endlist) {
    r0 = data.rMin;
    data.endlist = true;
  }
 
  // Loop through all trigger space points
  //
  for (; r0!=rbe[0]; ++r0) {
    data.nOneSeeds = 0;
    data.mapOneSeeds.erase(data.mapOneSeeds.begin(), data.mapOneSeeds.end());
 
    float R = (*r0)->radius();
    if (R<m_r2min) continue;
    if (R>m_r2max) break;
 
    bool pix = true;
    if ((*r0)->spacepoint->clusterList().second) pix = false;
    const Trk::Surface* sur0 = (*r0)->sur();
    float X = (*r0)->x();
    float Y = (*r0)->y();
    float Z = (*r0)->z();
    int Nb = 0;
    
    // Bottom links production
    //
    for (int i=0; i!=NB; ++i) {
      for (r=rb[i]; r!=rbe[i]; ++r) {
        float Rb =(*r)->radius();
        if (Rb<m_r1min) {
          rb[i]=r;
          continue;
        }
        if (Rb>m_r1max) break;
 
        float dR = R-Rb;
        if (dR<m_drmin) break;
 
        if (dR>m_drmax || (*r)->sur()==sur0) continue;
 
        if ( !pix && !(*r)->spacepoint->clusterList().second) continue;
 
        float dx = X-(*r)->x();
        float dy = Y-(*r)->y();
        float dZ = Z-(*r)->z();
        data.Tz[Nb] = dZ/std::sqrt(dx*dx+dy*dy);
        if (data.Tz[Nb]<m_dzdrmin || data.Tz[Nb]>m_dzdrmax) continue;
        data.Zo[Nb] = Z-R*data.Tz[Nb];
 
        // Comparison with vertices Z coordinates
        //
        if (!isZCompatible(data, data.Zo[Nb], Rb, data.Tz[Nb])) continue;
        data.SP[Nb] = (*r);
        if (++Nb==m_maxsizeSP) goto breakb;
      }
    }
  breakb:
    if (!Nb || Nb==m_maxsizeSP) continue;
    int Nt = Nb;
    
    // Top   links production
    //
    for (int i=0; i!=NT; ++i) {
      for (r=rt[i]; r!=rte[i]; ++r) {
        float Rt =(*r)->radius();
        float dR = Rt-R;
        if (dR<m_drmin || Rt<m_r3min) {
          rt[i]=r;
          continue;
        }
        if (Rt>m_r3max || dR>m_drmax) break;
 
        if ((*r)->sur()==sur0) continue;
 
        float dx = X-(*r)->x();
        float dy = Y-(*r)->y();
        float dZ = (*r)->z()-Z;
        data.Tz[Nt]   = dZ/std::sqrt(dx*dx+dy*dy);
        if (data.Tz[Nt]<m_dzdrmin || data.Tz[Nt]>m_dzdrmax) continue;
        data.Zo[Nt]   = Z-R*data.Tz[Nt];
 
        // Comparison with vertices Z coordinates
        //
        if (!isZCompatible(data, data.Zo[Nt], Rt, data.Tz[Nt])) continue;
        data.SP[Nt] = (*r);
        if (++Nt==m_maxsizeSP) goto breakt;
      }
    }
    
  breakt:
    if (!(Nt-Nb)) continue;
 
    float covr0 = (*r0)->covr ();
    float covz0 = (*r0)->covz ();
    float ax   = X/R;
    float ay   = Y/R;
    
    for (int i=0; i!=Nt; ++i) {
      float dx = data.SP[i]->x()-X;
      float dy = data.SP[i]->y()-Y;
      float dz = data.SP[i]->z()-Z;
      float x  = dx*ax+dy*ay;
      float y  =-dx*ay+dy*ax;
      float r2 = 1.f/(x*x+y*y);
      float dr  = std::sqrt(r2);
 
      i < Nb ?  data.Tz[i] = -dz*dr :  data.Tz[i] = dz*dr;
 
      data.U [i] = x*r2;
      data.V [i] = y*r2;
      data.Er[i] = (covz0+data.SP[i]->covz()+data.Tz[i]*data.Tz[i]*(covr0+data.SP[i]->covr()))*r2;
      data.R [i] = dr;
    }
 
    // Three space points comparison
    //
    for (int b=Nb-1; b>=0; --b) {
      float SA  = 1.f+data.Tz[b]*data.Tz[b];
      for (int t=Nb;  t!=Nt; ++t) {
        float cof = COF;
        if (!data.SP[t]->spacepoint->clusterList().second) cof = COFP;
 
        float Ts = .5f*(data.Tz[b]+data.Tz[t]);
        float dT =  data.Tz[b]-data.Tz[t];
        dT        = dT*dT-data.Er[b]-data.Er[t]-2.*data.R[b]*data.R[t]*(Ts*Ts*covr0+covz0);
 
        if (dT > 0. && dT > (m_iptmin*m_iptmin)*cof*SA) continue;
        float dU = data.U[t]-data.U[b];
        if (dU == 0.) continue;
        float A  = (data.V[t]-data.V[b])/dU;
        float B  =  data.V[t]-A*data.U[t];
        float S2 = 1.f+A*A;
        float S  = std::sqrt(S2);
        float BK = std::abs(B*K);
        if (BK > m_iptmin*S || BK < m_iptmax*S) continue; // Momentum    cut
        if (dT > 0. && dT  > (BK*BK/S2)*cof*SA) continue; // Polar angle cut
 
        float Im = std::abs((A-B*R)*R);
        if (Im > m_diver) continue;
 
        newOneSeed(data, data.SP[b]->spacepoint,(*r0)->spacepoint,data.SP[t]->spacepoint,data.Zo[b],Im);
      }
    }
    nseed += data.mapOneSeeds.size();
    fillSeeds(data);
    if (nseed>=m_maxsize) {
      data.endlist=false;
      ++r0;
      data.rMin = r0;
      return;
    } 
  }
}

writeNtuple()

void InDet::SiSpacePointsSeedMaker_LowMomentum::writeNtuple ( const SiSpacePointsSeed * seed, const Trk::Track * track, int seedType, long eventNumber ) const

overridevirtual

Definition at line 1187 of file SiSpacePointsSeedMaker_LowMomentum.cxx.

                                                                                                                   {
}

Member Data Documentation

m_beamSpotKey

SG::ReadCondHandleKey<InDet::BeamSpotData> InDet::SiSpacePointsSeedMaker_LowMomentum::m_beamSpotKey { this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" }

private

Definition at line 140 of file SiSpacePointsSeedMaker_LowMomentum.h.

{ this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" };

m_diver

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_diver {this, "maxdImpact", 7.}

private

Definition at line 167 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxdImpact", 7.};

m_drmax

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_drmax {this, "maxdRadius", 200.}

private

Definition at line 160 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxdRadius", 200.};

m_drmin

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_drmin {this, "mindRadius", 10.}

private

Definition at line 159 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "mindRadius", 10.};

m_dzdrmax

float InDet::SiSpacePointsSeedMaker_LowMomentum::m_dzdrmax {0.}

private

Definition at line 189 of file SiSpacePointsSeedMaker_LowMomentum.h.

{0.};

m_dzdrmin

float InDet::SiSpacePointsSeedMaker_LowMomentum::m_dzdrmin {0.}

private

Definition at line 188 of file SiSpacePointsSeedMaker_LowMomentum.h.

{0.};

m_dzdrver

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_dzdrver {this, "maxdZdRver", 0.02}

private

Definition at line 166 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxdZdRver", 0.02};

m_dzver

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_dzver {this, "maxdZver", 5.}

private

Definition at line 165 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxdZver", 5.};

m_etamax

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_etamax {this, "etaMax", 2.7}

private

Definition at line 173 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "etaMax", 2.7};

m_fieldCondObjInputKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> InDet::SiSpacePointsSeedMaker_LowMomentum::m_fieldCondObjInputKey

private

Initial value:

{this, "AtlasFieldCacheCondObj", "fieldCondObj",
                                                                           "Name of the Magnetic Field conditions object key"}

Definition at line 142 of file SiSpacePointsSeedMaker_LowMomentum.h.

                                                                       {this, "AtlasFieldCacheCondObj", "fieldCondObj",
                                                                           "Name of the Magnetic Field conditions object key"};

m_fieldScale

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_fieldScale {this, "fieldScale", 1.}

private

Definition at line 175 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "fieldScale", 1.};

m_fNmax

int InDet::SiSpacePointsSeedMaker_LowMomentum::m_fNmax {0}

private

Definition at line 182 of file SiSpacePointsSeedMaker_LowMomentum.h.

{0};

m_initialized

bool InDet::SiSpacePointsSeedMaker_LowMomentum::m_initialized {false}

private

Definition at line 180 of file SiSpacePointsSeedMaker_LowMomentum.h.

{false};

m_iptmax

float InDet::SiSpacePointsSeedMaker_LowMomentum::m_iptmax {1.f/400.f}

private

Definition at line 192 of file SiSpacePointsSeedMaker_LowMomentum.h.

{1.f/400.f};

m_iptmin

float InDet::SiSpacePointsSeedMaker_LowMomentum::m_iptmin {0.}

private

Definition at line 191 of file SiSpacePointsSeedMaker_LowMomentum.h.

{0.};

m_maxOneSize

IntegerProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_maxOneSize {this, "maxSeedsForSpacePoint", 5}

private

Definition at line 153 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxSeedsForSpacePoint", 5};

m_maxsize

IntegerProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_maxsize {this, "maxSize", 2000}

private

Definition at line 151 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxSize", 2000};

m_maxsizeSP

IntegerProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_maxsizeSP {this, "maxSizeSP", 1500}

private

Definition at line 152 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxSizeSP", 1500};

m_outputlevel

int InDet::SiSpacePointsSeedMaker_LowMomentum::m_outputlevel {0}

private

Definition at line 181 of file SiSpacePointsSeedMaker_LowMomentum.h.

{0};

m_pixel

BooleanProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_pixel {this, "usePixel", true}

private

Definition at line 148 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "usePixel", true};

m_prdToTrackMap

SG::ReadHandleKey<Trk::PRDtoTrackMap> InDet::SiSpacePointsSeedMaker_LowMomentum::m_prdToTrackMap {this,"PRDtoTrackMap","","option PRD-to-track association"}

private

Definition at line 139 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this,"PRDtoTrackMap","","option PRD-to-track association"};

m_ptmax

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_ptmax {this, "pTmax", 500.}

private

Definition at line 168 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "pTmax", 500.};

m_ptmin

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_ptmin {this, "pTmin", 100.}

private

Definition at line 174 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "pTmin", 100.};

m_r1max

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_r1max {this, "maxRadius1", 600.}

private

Definition at line 155 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxRadius1", 600.};

m_r1min

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_r1min {this, "minRadius1", 0.}

private

Definition at line 154 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "minRadius1", 0.};

m_r2max

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_r2max {this, "maxRadius2", 600.}

private

Definition at line 157 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxRadius2", 600.};

m_r2min

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_r2min {this, "minRadius2", 0.}

private

Definition at line 156 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "minRadius2", 0.};

m_r3max

float InDet::SiSpacePointsSeedMaker_LowMomentum::m_r3max {0.}

private

Definition at line 190 of file SiSpacePointsSeedMaker_LowMomentum.h.

{0.};

m_r3min

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_r3min {this, "minRadius3", 0.}

private

Definition at line 158 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "minRadius3", 0.};

m_r_rmax

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_r_rmax {this, "radMax", 200.}

private

Definition at line 163 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "radMax", 200.};

m_r_rstep

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_r_rstep {this, "radStep", 2.}

private

Definition at line 164 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "radStep", 2.};

m_r_size

int InDet::SiSpacePointsSeedMaker_LowMomentum::m_r_size {0}

private

Definition at line 183 of file SiSpacePointsSeedMaker_LowMomentum.h.

{0};

m_rfz_b

int InDet::SiSpacePointsSeedMaker_LowMomentum::m_rfz_b[SizeRFZ] {}

private

Definition at line 184 of file SiSpacePointsSeedMaker_LowMomentum.h.

{};

m_rfz_ib

int InDet::SiSpacePointsSeedMaker_LowMomentum::m_rfz_ib[SizeRFZ][SizeI] {}

private

Definition at line 186 of file SiSpacePointsSeedMaker_LowMomentum.h.

{};

m_rfz_it

int InDet::SiSpacePointsSeedMaker_LowMomentum::m_rfz_it[SizeRFZ][SizeI] {}

private

Definition at line 187 of file SiSpacePointsSeedMaker_LowMomentum.h.

{};

m_rfz_t

int InDet::SiSpacePointsSeedMaker_LowMomentum::m_rfz_t[SizeRFZ] {}

private

Definition at line 185 of file SiSpacePointsSeedMaker_LowMomentum.h.

{};

m_sct

BooleanProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_sct {this, "useSCT", true}

private

Definition at line 149 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "useSCT", true};

m_sF

float InDet::SiSpacePointsSeedMaker_LowMomentum::m_sF {0.}

private

Definition at line 193 of file SiSpacePointsSeedMaker_LowMomentum.h.

{0.};

m_spacepointsOverlap

SG::ReadHandleKey<SpacePointOverlapCollection> InDet::SiSpacePointsSeedMaker_LowMomentum::m_spacepointsOverlap {this, "SpacePointsOverlapName", "OverlapSpacePoints"}

private

Definition at line 138 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "SpacePointsOverlapName", "OverlapSpacePoints"};

m_spacepointsPixel

SG::ReadHandleKey<SpacePointContainer> InDet::SiSpacePointsSeedMaker_LowMomentum::m_spacepointsPixel {this, "SpacePointsPixelName", "PixelSpacePoints", "Pixel space points container"}

private

Definition at line 137 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "SpacePointsPixelName", "PixelSpacePoints", "Pixel space points container"};

m_spacepointsSCT

SG::ReadHandleKey<SpacePointContainer> InDet::SiSpacePointsSeedMaker_LowMomentum::m_spacepointsSCT {this, "SpacePointsSCTName", "SCT_SpacePoints", "SCT space points container"}

private

Definition at line 136 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "SpacePointsSCTName", "SCT_SpacePoints", "SCT space points container"};

m_useOverlap

BooleanProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_useOverlap {this, "useOverlapSpCollection", false}

private

Definition at line 150 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "useOverlapSpCollection", false};

m_zmax

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_zmax {this, "maxZ", +250.}

private

Definition at line 162 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "maxZ", +250.};

m_zmin

FloatProperty InDet::SiSpacePointsSeedMaker_LowMomentum::m_zmin {this, "minZ", -250.}

private

Definition at line 161 of file SiSpacePointsSeedMaker_LowMomentum.h.

{this, "minZ", -250.};

---

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

• SiSpacePointsSeedMaker_LowMomentum.h
• SiSpacePointsSeedMaker_LowMomentum.cxx