Class for track candidates generation using space points information for standard Atlas geometry. More...

#include <SiSpacePointsSeedMaker_BeamGas.h>

Inheritance diagram for InDet::SiSpacePointsSeedMaker_BeamGas:

Collaboration diagram for InDet::SiSpacePointsSeedMaker_BeamGas:

Public Member Functions
Standard tool methods
	SiSpacePointsSeedMaker_BeamGas (const std::string &, const std::string &, const IInterface *)
virtual	~SiSpacePointsSeedMaker_BeamGas ()=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
virtual void	find2Sp (EventData &data, const std::list< Trk::Vertex > &lv) const override
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"}
Read handle for conditions object to get the field cache
SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCondObjInputKey
Properties, which will not be changed after construction
BooleanProperty	m_pixel {this, "usePixel", true}
BooleanProperty	m_sct {this, "useSCT""useSCT", true}
BooleanProperty	m_useOverlap {this, "useOverlapSpCollection", true}
IntegerProperty	m_maxsize {this, "maxSize", 20000}
IntegerProperty	m_maxsizeSP {this, "maxSizeSP", 1500}
IntegerProperty	m_maxOneSize {this, "maxSeedsForSpacePoint", 5}
FloatProperty	m_drmax {this, "maxdRadius", 270.}
FloatProperty	m_zmin {this, "minZ", -5000.}
FloatProperty	m_zmax {this, "maxZ", +5000.}
FloatProperty	m_r_rmax {this, "radMax", 600.}
FloatProperty	m_r_rstep {this, "radStep", 2.}
FloatProperty	m_r1min {this, "minRadius1", 0.}
FloatProperty	m_r1max {this, "maxRadius1", 450.}
FloatProperty	m_r2min {this, "minRadius2", 0.}
FloatProperty	m_r2max {this, "maxRadius2", 600.}
FloatProperty	m_r3min {this, "minRadius3", 0.}
FloatProperty	m_r3max {this, "maxRadius3", 600.}
FloatProperty	m_drmin {this, "mindRadius", 10.}
FloatProperty	m_diver {this, "maxdImpact", 10.}
FloatProperty	m_diverpps {this, "maxdImpactPPS", 1.7}
Properties, which can be updated in buildFrameWork method which is called in initialize method
FloatProperty	m_ptmin {this, "pTmin", 1000.}
FloatProperty	m_etamax {this, "etaMax", 5.3}
FloatProperty	m_fieldScale {this, "fieldScale", 1.}
Data members, which are updated in only initialize
bool	m_initialized {false}
int	m_outputlevel {0}
int	m_r_size {0}
int	m_fNmax {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_dzdrmax {0.}
float	m_dzdrmin {0.}
float	m_COF {0.}
float	m_ipt {0.}
float	m_ipt2 {0.}
float	m_sF {0.}

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

Disallow default instantiation, copy, assignment
	SiSpacePointsSeedMaker_BeamGas ()=delete
	SiSpacePointsSeedMaker_BeamGas (const SiSpacePointsSeedMaker_BeamGas &)=delete
SiSpacePointsSeedMaker_BeamGas &	operator= (const SiSpacePointsSeedMaker_BeamGas &)=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 (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 &) const
void	findNext (EventData &data) const
bool	isZCompatible (float &) const
bool	isUsed (const Trk::SpacePoint *sp, 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 void	convertToBeamFrameWork (EventData &data, const Trk::SpacePoint const &sp, float r)

Detailed Description

Class for track candidates generation using space points information for standard Atlas geometry.

In AthenaMT, event dependent cache inside SiSpacePointsSeedMaker_BeamGas is not preferred. SiSpacePointsSeedMakerEventData = EventData class holds event dependent data for SiSpacePointsSeedMaker_BeamGas. Its object is instantiated in SiSPSeededTrackFinder::execute.

Definition at line 54 of file SiSpacePointsSeedMaker_BeamGas.h.

Member Enumeration Documentation

◆ Size

enum InDet::SiSpacePointsSeedMaker_BeamGas::Size

private

enum for array sizes

Enumerator
SizeRF
SizeZ
SizeRFZ
SizeI

Definition at line 127 of file SiSpacePointsSeedMaker_BeamGas.h.

              {SizeRF=53,
               SizeZ=11,
               SizeRFZ=SizeRF*SizeZ,
               SizeI=12};

Constructor & Destructor Documentation

◆ SiSpacePointsSeedMaker_BeamGas() [1/3]

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

Definition at line 27 of file SiSpacePointsSeedMaker_BeamGas.cxx.

  : base_class(t, n, p)
{
}

◆ ~SiSpacePointsSeedMaker_BeamGas()

virtual InDet::SiSpacePointsSeedMaker_BeamGas::~SiSpacePointsSeedMaker_BeamGas ( )

virtualdefault

◆ SiSpacePointsSeedMaker_BeamGas() [2/3]

InDet::SiSpacePointsSeedMaker_BeamGas::SiSpacePointsSeedMaker_BeamGas ( )

privatedelete

◆ SiSpacePointsSeedMaker_BeamGas() [3/3]

InDet::SiSpacePointsSeedMaker_BeamGas::SiSpacePointsSeedMaker_BeamGas ( const SiSpacePointsSeedMaker_BeamGas & )

privatedelete

Member Function Documentation

◆ buildBeamFrameWork()

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

private

Definition at line 680 of file SiSpacePointsSeedMaker_BeamGas.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   = atan2(ty,tx);
  double th   = acos(1./sqrt(1.+tx*tx+ty*ty));
  double sint = sin(th);
  double cost = cos(th);
  double sinp = sin(ph);
  double cosp = 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_BeamGas::buildFrameWork ( )

private

Definition at line 589 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  m_ptmin = std::max( std::abs(m_ptmin), float(300.*m_fieldScale));
  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_COF       =  134*.05*9.;
  m_ipt       = 1.f/std::abs(.9f*m_ptmin);
  m_ipt2      = m_ipt*m_ipt;
 
  // Build radius sorted containers
  //
  m_r_size = static_cast<int>((m_r_rmax+.1)/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_b [a]    = 9;
    m_rfz_it[a][3] = a+1;  m_rfz_ib[a][3] = a+1;
    m_rfz_it[a][4] = b+1;  m_rfz_ib[a][4] = b+1;
    m_rfz_it[a][5] = c+1;  m_rfz_ib[a][5] = c+1;
    m_rfz_it[a][6] = a-1;  m_rfz_ib[a][6] = a-1;
    m_rfz_it[a][7] = b-1;  m_rfz_ib[a][7] = b-1;
    m_rfz_it[a][8] = c-1;  m_rfz_ib[a][8] = c-1;
      } else if (z> 5) {
 
    m_rfz_b [a]    = 6;   m_rfz_t [a]    = 6;
    m_rfz_ib[a][3] = a-1;  m_rfz_it[a][3] = a-1;
    m_rfz_ib[a][4] = b-1;  m_rfz_it[a][4] = b-1;
    m_rfz_ib[a][5] = c-1;  m_rfz_it[a][5] = c-1;
 
    if (z<10) {
 
      m_rfz_t [a]    = 9;  m_rfz_b [a]    = 9;
      m_rfz_it[a][6] = a+1; m_rfz_ib[a][6] = a+1;
      m_rfz_it[a][7] = b+1; m_rfz_ib[a][7] = b+1;
      m_rfz_it[a][8] = c+1; m_rfz_ib[a][8] = c+1;
    }
      } else {
 
    m_rfz_b [a]    = 6;   m_rfz_t [a]    = 6;
    m_rfz_ib[a][3] = a+1;   m_rfz_it[a][3] = a+1;
    m_rfz_ib[a][4] = b+1;   m_rfz_it[a][4] = b+1;
    m_rfz_ib[a][5] = c+1;   m_rfz_it[a][5] = c+1;
 
    if (z>0) {
 
      m_rfz_t [a]    = 9; m_rfz_b [a]    = 9;
      m_rfz_it[a][6] = a-1; m_rfz_ib[a][6] = a-1;
      m_rfz_it[a][7] = b-1; m_rfz_ib[a][7] = b-1;
      m_rfz_it[a][8] = c-1; m_rfz_ib[a][8] = c-1;
    }
      }
    }
  }
}

◆ convertToBeamFrameWork()

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

staticprivate

Definition at line 715 of file SiSpacePointsSeedMaker_BeamGas.cxx.

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

◆ dump()

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

overridevirtual

Definition at line 421 of file SiSpacePointsSeedMaker_BeamGas.cxx.

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

◆ dumpConditions()

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

private

Definition at line 433 of file SiSpacePointsSeedMaker_BeamGas.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<<"| 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<<"| |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       impact        | "
     <<std::setw(12)<<std::setprecision(5)<<m_diver
     <<"                              |"<<endmsg;
  out<<"| max       impact pps    | "
     <<std::setw(12)<<std::setprecision(5)<<m_diverpps
     <<"                              |"<<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_BeamGas::dumpEvent	(	EventData &	data,
		MsgStream &	out )

staticprivate

Definition at line 550 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  out<<"|---------------------------------------------------------------------|"
     <<endmsg;
  out<<"| ns                    | "
     <<std::setw(12)<<data.ns
     <<"                              |"<<endmsg;
  out<<"| 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_BeamGas::erase ( EventData & data )

staticprivate

Definition at line 779 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  for (int i=0; i!=data.nr;    ++i) {
    int n = data.r_index[i];
    data.r_map[n] = 0;
    data.r_Sorted[n].erase(data.r_Sorted[n].begin(), data.r_Sorted[n].end());
  }
 
  for (int i=0; i!=data.nrf;   ++i) {
    int n = data.rf_index[i];
    data.rf_map[n] = 0;
    data.rf_Sorted[n].erase(data.rf_Sorted[n].begin(), data.rf_Sorted[n].end());
  }
 
  for (int i=0; i!=data.nrfz;  ++i) {
    int n = data.rfz_index[i];
    data.rfz_map[n] = 0;
    data.rfz_Sorted[n].erase(data.rfz_Sorted[n].begin(), data.rfz_Sorted[n].end());
  }
 
  data.state = 0;
  data.ns    = 0;
  data.nsaz  = 0;
  data.nr    = 0;
  data.nrf   = 0;
  data.nrfz  = 0;
}

◆ fillLists()

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

private

Definition at line 730 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  constexpr float pi2 = 2.*M_PI;
 
  std::vector<InDet::SiSpacePointForSeed*>::iterator r;
  
  for (int i=0; i!= m_r_size;  ++i) {
 
    if (!data.r_map[i]) continue;
    r = data.r_Sorted[i].begin();
 
    while (r!=data.r_Sorted[i].end()) {
      
      // 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;
 
      data.rf_Sorted[f].push_back(*r);
      if (!data.rf_map[f]++) data.rf_index[data.nrf++] = f;
 
      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].erase(r++);
    }
    data.r_map[i] = 0;
  }
  data.nr    = 0;
  data.state = 0;
}

◆ fillSeeds()

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

staticprivate

Definition at line 1187 of file SiSpacePointsSeedMaker_BeamGas.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_BeamGas::finalize ( )

overridevirtual

Definition at line 77 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  return AlgTool::finalize();
}

◆ find2Sp()

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

overridevirtual

Definition at line 324 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  int mode = 0;
  if (lv.begin()!=lv.end()) mode = 1;
  
  if (!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_BeamGas::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 355 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  int mode = 2;
  if (lv.begin()!=lv.end()) mode = 3;
 
  if (!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(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_BeamGas::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 380 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  find3Sp(ctx, data, lv);
}

◆ findNext()

void InDet::SiSpacePointsSeedMaker_BeamGas::findNext ( EventData & data ) const

private

Definition at line 572 of file SiSpacePointsSeedMaker_BeamGas.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(data);
  else if (data.mode==5 || data.mode==6) production3Sp(data);
 
  data.i_seed = data.l_seeds.begin();
  ++data.nlist;
}                       

◆ findVSp()

void InDet::SiSpacePointsSeedMaker_BeamGas::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 391 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  int mode = 5;
  if (lv.begin()!=lv.end()) mode = 6;
  
  if (!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(data);
  }
  data.i_seed = data.l_seeds.begin();
 
  if (m_outputlevel<=0) {
    data.nprint=1;
    dump(data, msg(MSG::DEBUG));
  }
}

◆ getWriteNtupleBoolProperty()

bool InDet::SiSpacePointsSeedMaker_BeamGas::getWriteNtupleBoolProperty ( ) const

overridevirtual

Definition at line 1219 of file SiSpacePointsSeedMaker_BeamGas.cxx.

                                                                          {
    return false;
}

◆ initialize()

StatusCode InDet::SiSpacePointsSeedMaker_BeamGas::initialize ( )

overridevirtual

Definition at line 37 of file SiSpacePointsSeedMaker_BeamGas.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_BeamGas::initializeEventData ( EventData & data ) const

private

Definition at line 1204 of file SiSpacePointsSeedMaker_BeamGas.cxx.

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

◆ isUsed()

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

inlineprivate

Definition at line 258 of file SiSpacePointsSeedMaker_BeamGas.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_BeamGas::isZCompatible ( float & Zv ) const

private

Definition at line 1111 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  return Zv > m_zmin && Zv < m_zmax;
}

◆ newEvent()

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

overridevirtual

Definition at line 86 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  if (!m_pixel && !m_sct) return;
 
  if (not data.initialized) initializeEventData(data);
 
  erase(data);
  buildBeamFrameWork(data);
 
  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_BeamGas: Failed to retrieve AtlasFieldCacheCondObj with key " << m_fieldCondObjInputKey.key());
    return;
  }
  fieldCondObj->getInitializedCache (fieldCache);
 
  if (fieldCache.solenoidOn()) {
    fieldCache.getFieldZR(gP,f);
 
    data.K = 2./(300.*f[2]);
  } else {
    data.K = 2./(300.* 5. );
  }
 
  data.i_spforseed = data.l_spforseed.begin();
 
  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) {
 
    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;
 
      int ir = static_cast<int>(r*irstep);
      InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
      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;
 
      int ir = static_cast<int>(r*irstep);
      InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
      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 overlap space points containers from store gate 
    //
    if (m_useOverlap) {
 
      SG::ReadHandle<SpacePointOverlapCollection> spacepointsOverlap{m_spacepointsOverlap, ctx};
      if (spacepointsOverlap.isValid()) {
 
        for (const Trk::SpacePoint* sp: *spacepointsOverlap) {
 
      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;
 
      int ir = static_cast<int>(r*irstep);
      InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
      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_BeamGas::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 1064 of file SiSpacePointsSeedMaker_BeamGas.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_BeamGas::newRegion	(	const EventContext &	ctx,
		EventData &	data,
		const std::vector< IdentifierHash > &	vPixel,
		const std::vector< IdentifierHash > &	vSCT ) const

overridevirtual

Definition at line 207 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  if (!m_pixel && !m_sct) return;
 
  if (not data.initialized) initializeEventData(data);
 
  erase(data);
  buildBeamFrameWork(data);
 
  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_BeamGas: Failed to retrieve AtlasFieldCacheCondObj with key " << m_fieldCondObjInputKey.key());
    return;
  }
  fieldCondObj->getInitializedCache (fieldCache);
 
  if (fieldCache.solenoidOn()) {
    fieldCache.getFieldZR(gP,f);
 
    data.K = 2./(300.*f[2]);
  } else {
    data.K = 2./(300.* 5. );
  }
 
  data.i_spforseed = data.l_spforseed.begin();
 
  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;
      int ir = static_cast<int>(r*irstep);
      InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
      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: vPixel) {
    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;
      int ir = static_cast<int>(r*irstep);
      InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
      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_BeamGas::newRegion	(	const EventContext &	ctx,
		EventData &	data,
		const std::vector< IdentifierHash > &	vPixel,
		const std::vector< IdentifierHash > &	vSCT,
		const IRoiDescriptor &	iRD ) const

overridevirtual

Definition at line 312 of file SiSpacePointsSeedMaker_BeamGas.cxx.

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

◆ newSeed() [1/2]

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

staticprivate

Definition at line 1144 of file SiSpacePointsSeedMaker_BeamGas.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_BeamGas::newSeed	(	EventData &	data,
		const Trk::SpacePoint *&	p1,
		const Trk::SpacePoint *&	p2,
		const Trk::SpacePoint *&	p3,
		const float &	z )

staticprivate

Definition at line 1165 of file SiSpacePointsSeedMaker_BeamGas.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_BeamGas::newSpacePoint	(	EventData &	data,
		const Trk::SpacePoint *const &	sp )

staticprivate

Definition at line 1120 of file SiSpacePointsSeedMaker_BeamGas.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;
}

◆ next()

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

overridevirtual

Definition at line 1099 of file SiSpacePointsSeedMaker_BeamGas.cxx.

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

◆ operator=()

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

privatedelete

◆ production2Sp()

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

staticprivate

Definition at line 811 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  data.endlist = true;
}

◆ production3Sp() [1/2]

void InDet::SiSpacePointsSeedMaker_BeamGas::production3Sp ( EventData & data ) const

private

Definition at line 820 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{ 
  if (data.nsaz<3) return;
 
  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);
      if (!data.endlist) {
        data.fNmin = f;
        data.zMin = z; return;
      }
    }
  }
  data.endlist = true;
}

◆ production3Sp() [2/2]

void InDet::SiSpacePointsSeedMaker_BeamGas::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 ) const

private

Definition at line 869 of file SiSpacePointsSeedMaker_BeamGas.cxx.

{
  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;
 
    const Trk::SpacePoint* SP0 = (*r0)->spacepoint;
 
    bool pix = true;
    if (SP0->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 Tz = (Z-(*r)->z())/dR;
 
    if (Tz < m_dzdrmin || Tz > m_dzdrmax) continue;
    
    // Comparison with vertices Z coordinates
    //
    if (pix) {
      float Zo = Z-R*Tz;
          if (!isZCompatible(Zo)) 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 Tz = ((*r)->z()-Z)/dR;
 
    if (Tz < m_dzdrmin || Tz > m_dzdrmax) continue;
 
    // Comparison with vertices Z coordinates
    //
    if (pix) {
      float Zo = Z-R*Tz; if (!isZCompatible(Zo)) 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) {
 
      InDet::SiSpacePointForSeed* sp = data.SP[i];
 
      float dx  = sp->x()-X;
      float dy  = sp->y()-Y;
      float dz  = sp->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);
      float tz  = dz*dr; if (i < Nb) tz = -tz;
 
      data.Tz[i] = tz;
      data.Zo[i] = Z-R*tz;
      data.R [i] = dr;
      data.U [i] = x*r2;
      data.V [i] = y*r2;
      data.Er[i] = (covz0+sp->covz()+tz*tz*(covr0+sp->covr()))*r2;
    }
 
    float imc   = m_diver;
    float imcs  = m_diverpps;
    float ipt2  = m_ipt2;
    float K     = data.K;
    float K2    = K*K;
    float COF   = m_COF;
    float ipt2K = ipt2/K2;
    float ipt2C = ipt2*COF;
    float COFK  = COF*K2;
    covr0      *= 2.;
    covz0      *= 2.;
    
    // Three space points comparison
    //
    for (int b=0; b!=Nb; ++b) {
    
      const Trk::SpacePoint* SPb = data.SP[b]->spacepoint;
 
      float  Zob  = data.Zo[b];
      float  Tzb  = data.Tz[b];
      float  Rb2r = data.R [b]*covr0;
      float  Rb2z = data.R [b]*covz0;
      float  Erb  = data.Er[b];
      float  Vb   = data.V [b];
      float  Ub   = data.U [b];
      float  Tzb2 = (1.f+Tzb*Tzb);
      float  CSA  = Tzb2*COFK;
      float ICSA  = Tzb2*ipt2C;
 
      for (int t=Nb;  t!=Nt; ++t) {
    
    float Ts  = .5f*(Tzb+data.Tz[t]);
    float dt  =     Tzb-data.Tz[t];
    float dT  = dt*dt-Erb-data.Er[t]-data.R[t]*(Ts*Ts*Rb2r+Rb2z);
    if ( dT > ICSA) continue;
    float dU  = data.U[t]-Ub; if (dU == 0. ) continue;
    float A   = (data.V[t]-Vb)/dU;
    float S2  = 1.f+A*A;
    float B   = Vb-A*Ub;
    float B2  = B*B;
    if (B2  > ipt2K*S2 || dT*S2 > B2*CSA) continue;
    float Im  = std::abs((A-B*R)*R);
    
    if (pix) {
      if (                                             Im > imc ) continue;
      if (data.SP[t]->spacepoint->clusterList().second && Im > imcs) continue;
    }
    newOneSeed(data, SPb, SP0, data.SP[t]->spacepoint, Zob, Im);
      }
    }
    nseed += data.mapOneSeeds.size();
    fillSeeds(data);
    if (nseed>=m_maxsize) {
      data.endlist=false;
      ++r0;
      data.rMin = r0;
      return;
    } 
  }
}

◆ writeNtuple()

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

overridevirtual

Definition at line 1216 of file SiSpacePointsSeedMaker_BeamGas.cxx.

1216 {

1217}

Member Data Documentation

◆ m_beamSpotKey

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

private

Definition at line 142 of file SiSpacePointsSeedMaker_BeamGas.h.

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

◆ m_COF

float InDet::SiSpacePointsSeedMaker_BeamGas::m_COF {0.}

private

Definition at line 193 of file SiSpacePointsSeedMaker_BeamGas.h.

193{0.};

◆ m_diver

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_diver {this, "maxdImpact", 10.}

private

Definition at line 170 of file SiSpacePointsSeedMaker_BeamGas.h.

170{this, "maxdImpact", 10.};

◆ m_diverpps

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_diverpps {this, "maxdImpactPPS", 1.7}

private

Definition at line 171 of file SiSpacePointsSeedMaker_BeamGas.h.

171{this, "maxdImpactPPS", 1.7};

◆ m_drmax

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_drmax {this, "maxdRadius", 270.}

private

Definition at line 158 of file SiSpacePointsSeedMaker_BeamGas.h.

158{this, "maxdRadius", 270.};

◆ m_drmin

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

private

Definition at line 169 of file SiSpacePointsSeedMaker_BeamGas.h.

169{this, "mindRadius", 10.};

◆ m_dzdrmax

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

private

Definition at line 191 of file SiSpacePointsSeedMaker_BeamGas.h.

191{0.};

◆ m_dzdrmin

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

private

Definition at line 192 of file SiSpacePointsSeedMaker_BeamGas.h.

192{0.};

◆ m_etamax

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_etamax {this, "etaMax", 5.3}

private

Definition at line 177 of file SiSpacePointsSeedMaker_BeamGas.h.

177{this, "etaMax", 5.3};

◆ m_fieldCondObjInputKey

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

private

Initial value:

{this, "AtlasFieldCacheCondObj", "fieldCondObj",

"Name of the Magnetic Field conditions object key"}

Definition at line 147 of file SiSpacePointsSeedMaker_BeamGas.h.

147 {this, "AtlasFieldCacheCondObj", "fieldCondObj",

148 "Name of the Magnetic Field conditions object key"};

◆ m_fieldScale

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

private

Definition at line 178 of file SiSpacePointsSeedMaker_BeamGas.h.

178{this, "fieldScale", 1.};

◆ m_fNmax

int InDet::SiSpacePointsSeedMaker_BeamGas::m_fNmax {0}

private

Definition at line 186 of file SiSpacePointsSeedMaker_BeamGas.h.

186{0};

◆ m_initialized

bool InDet::SiSpacePointsSeedMaker_BeamGas::m_initialized {false}

private

Definition at line 183 of file SiSpacePointsSeedMaker_BeamGas.h.

183{false};

◆ m_ipt

float InDet::SiSpacePointsSeedMaker_BeamGas::m_ipt {0.}

private

Definition at line 194 of file SiSpacePointsSeedMaker_BeamGas.h.

194{0.};

◆ m_ipt2

float InDet::SiSpacePointsSeedMaker_BeamGas::m_ipt2 {0.}

private

Definition at line 195 of file SiSpacePointsSeedMaker_BeamGas.h.

195{0.};

◆ m_maxOneSize

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

private

Definition at line 157 of file SiSpacePointsSeedMaker_BeamGas.h.

157{this, "maxSeedsForSpacePoint", 5};

◆ m_maxsize

IntegerProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_maxsize {this, "maxSize", 20000}

private

Definition at line 155 of file SiSpacePointsSeedMaker_BeamGas.h.

155{this, "maxSize", 20000};

◆ m_maxsizeSP

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

private

Definition at line 156 of file SiSpacePointsSeedMaker_BeamGas.h.

156{this, "maxSizeSP", 1500};

◆ m_outputlevel

int InDet::SiSpacePointsSeedMaker_BeamGas::m_outputlevel {0}

private

Definition at line 184 of file SiSpacePointsSeedMaker_BeamGas.h.

184{0};

◆ m_pixel

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

private

Definition at line 152 of file SiSpacePointsSeedMaker_BeamGas.h.

152{this, "usePixel", true};

◆ m_prdToTrackMap

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

private

Definition at line 141 of file SiSpacePointsSeedMaker_BeamGas.h.

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

◆ m_ptmin

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_ptmin {this, "pTmin", 1000.}

private

Definition at line 176 of file SiSpacePointsSeedMaker_BeamGas.h.

176{this, "pTmin", 1000.};

◆ m_r1max

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_r1max {this, "maxRadius1", 450.}

private

Definition at line 164 of file SiSpacePointsSeedMaker_BeamGas.h.

164{this, "maxRadius1", 450.};

◆ m_r1min

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

private

Definition at line 163 of file SiSpacePointsSeedMaker_BeamGas.h.

163{this, "minRadius1", 0.};

◆ m_r2max

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

private

Definition at line 166 of file SiSpacePointsSeedMaker_BeamGas.h.

166{this, "maxRadius2", 600.};

◆ m_r2min

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

private

Definition at line 165 of file SiSpacePointsSeedMaker_BeamGas.h.

165{this, "minRadius2", 0.};

◆ m_r3max

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_r3max {this, "maxRadius3", 600.}

private

Definition at line 168 of file SiSpacePointsSeedMaker_BeamGas.h.

168{this, "maxRadius3", 600.};

◆ m_r3min

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

private

Definition at line 167 of file SiSpacePointsSeedMaker_BeamGas.h.

167{this, "minRadius3", 0.};

◆ m_r_rmax

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_r_rmax {this, "radMax", 600.}

private

Definition at line 161 of file SiSpacePointsSeedMaker_BeamGas.h.

161{this, "radMax", 600.};

◆ m_r_rstep

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

private

Definition at line 162 of file SiSpacePointsSeedMaker_BeamGas.h.

162{this, "radStep", 2.};

◆ m_r_size

int InDet::SiSpacePointsSeedMaker_BeamGas::m_r_size {0}

private

Definition at line 185 of file SiSpacePointsSeedMaker_BeamGas.h.

185{0};

◆ m_rfz_b

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

private

Definition at line 187 of file SiSpacePointsSeedMaker_BeamGas.h.

187{};

◆ m_rfz_ib

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

private

Definition at line 189 of file SiSpacePointsSeedMaker_BeamGas.h.

189{};

◆ m_rfz_it

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

private

Definition at line 190 of file SiSpacePointsSeedMaker_BeamGas.h.

190{};

◆ m_rfz_t

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

private

Definition at line 188 of file SiSpacePointsSeedMaker_BeamGas.h.

188{};

◆ m_sct

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

private

Definition at line 153 of file SiSpacePointsSeedMaker_BeamGas.h.

153{this, "useSCT""useSCT", true};

◆ m_sF

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

private

Definition at line 196 of file SiSpacePointsSeedMaker_BeamGas.h.

196{0.};

◆ m_spacepointsOverlap

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

private

Definition at line 140 of file SiSpacePointsSeedMaker_BeamGas.h.

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

◆ m_spacepointsPixel

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

private

Definition at line 139 of file SiSpacePointsSeedMaker_BeamGas.h.

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

◆ m_spacepointsSCT

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

private

Definition at line 138 of file SiSpacePointsSeedMaker_BeamGas.h.

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

◆ m_useOverlap

BooleanProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_useOverlap {this, "useOverlapSpCollection", true}

private

Definition at line 154 of file SiSpacePointsSeedMaker_BeamGas.h.

154{this, "useOverlapSpCollection", true};

◆ m_zmax

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_zmax {this, "maxZ", +5000.}

private

Definition at line 160 of file SiSpacePointsSeedMaker_BeamGas.h.

160{this, "maxZ", +5000.};

◆ m_zmin

FloatProperty InDet::SiSpacePointsSeedMaker_BeamGas::m_zmin {this, "minZ", -5000.}

private

Definition at line 159 of file SiSpacePointsSeedMaker_BeamGas.h.

159{this, "minZ", -5000.};

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

Public Member Functions

Private Attributes

Print internal tool parameters and status

Disallow default instantiation, copy, assignment

Detailed Description

Member Enumeration Documentation

◆ Size

Constructor & Destructor Documentation

◆ SiSpacePointsSeedMaker_BeamGas() [1/3]

◆ ~SiSpacePointsSeedMaker_BeamGas()

◆ SiSpacePointsSeedMaker_BeamGas() [2/3]

◆ SiSpacePointsSeedMaker_BeamGas() [3/3]

Member Function Documentation

◆ buildBeamFrameWork()

◆ buildFrameWork()

◆ convertToBeamFrameWork()

◆ dump()

◆ dumpConditions()

◆ dumpEvent()

◆ erase()

◆ fillLists()

◆ fillSeeds()

◆ finalize()

◆ find2Sp()

◆ find3Sp() [1/2]

◆ find3Sp() [2/2]

◆ findNext()

◆ findVSp()

◆ getWriteNtupleBoolProperty()

◆ initialize()

◆ initializeEventData()

◆ isUsed()

◆ isZCompatible()

◆ newEvent()

◆ newOneSeed()

◆ newRegion() [1/2]

◆ newRegion() [2/2]

◆ newSeed() [1/2]

◆ newSeed() [2/2]

◆ newSpacePoint()

◆ next()

◆ operator=()

◆ production2Sp()

◆ production3Sp() [1/2]

◆ production3Sp() [2/2]

◆ writeNtuple()

Member Data Documentation

◆ m_beamSpotKey

◆ m_COF

◆ m_diver

◆ m_diverpps

◆ m_drmax

◆ m_drmin

◆ m_dzdrmax

◆ m_dzdrmin

◆ m_etamax

◆ m_fieldCondObjInputKey

◆ m_fieldScale

◆ m_fNmax

◆ m_initialized

◆ m_ipt

◆ m_ipt2

◆ m_maxOneSize

◆ m_maxsize

◆ m_maxsizeSP

◆ m_outputlevel

◆ m_pixel

◆ m_prdToTrackMap

◆ m_ptmin

◆ m_r1max

◆ m_r1min

◆ m_r2max

◆ m_r2min

◆ m_r3max

◆ m_r3min

◆ m_r_rmax

◆ m_r_rstep

◆ m_r_size

◆ m_rfz_b

◆ m_rfz_ib