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

#include <SiSpacePointsSeedMaker_HeavyIon.h>

Inheritance diagram for InDet::SiSpacePointsSeedMaker_HeavyIon:

Collaboration diagram for InDet::SiSpacePointsSeedMaker_HeavyIon:

Public Member Functions
Standard tool methods
	SiSpacePointsSeedMaker_HeavyIon (const std::string &, const std::string &, const IInterface *)
virtual	~SiSpacePointsSeedMaker_HeavyIon ()=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
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::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"}
SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCondObjInputKey {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"}
SG::ReadHandleKey< Trk::PRDtoTrackMap >	m_prdToTrackMap {this,"PRDtoTrackMap","","option PRD-to-track association"}
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", 20000}
IntegerProperty	m_maxsizeSP {this, "maxSizeSP", 4000}
IntegerProperty	m_maxOneSize {this, "maxSeedsForSpacePoint", 6}
UnsignedIntegerProperty	m_maxNumberVertices {this, "maxNumberVertices", 1}
FloatProperty	m_r1minv {this, "minVRadius1", 0.}
FloatProperty	m_r1maxv {this, "maxVRadius1", 60.}
FloatProperty	m_r2minv {this, "minVRadius2", 70.}
FloatProperty	m_r2maxv {this, "maxVRadius2", 200.}
FloatProperty	m_r3minv {this, "minVRadius3", 0.}
FloatProperty	m_drmin_b {this, "mindRadius", 10.}
FloatProperty	m_drmin_t {this, "mindRadiusTop", 10.}
FloatProperty	m_drmax {this, "maxdRadius", 270.}
FloatProperty	m_zmin {this, "minZ", -250.}
FloatProperty	m_zmax {this, "maxZ", +250.}
FloatProperty	m_r_rmin {this, "radMin", 43.}
FloatProperty	m_r_rmax {this, "radMax", 600.}
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", 10.}
FloatProperty	m_diverpps {this, "maxdImpactPPS", 1.7}
FloatProperty	m_diversss {this, "maxdImpactSSS", 1000.}
BooleanProperty	m_useVertexPosition {this, "useVertexPosition", false}
BooleanProperty	m_writeNtuple {this, "WriteNtuple", false, "Flag to write Validation Ntuples"}
Properties, which can be updated in initialize
FloatProperty	m_etamax {this, "etaMax", 2.7}
FloatProperty	m_ptmin {this, "pTmin", 500.}
FloatProperty	m_ptmax {this, "pTmax", std::numeric_limits<float>::max()}
FloatProperty	m_fieldScale {this, "fieldScale", 1.}
Data members, which are not updated at all
float	m_drminv {20.}
bool	m_initialized {false}
int	m_outputlevel {0}
int	m_r_size {0}
int	m_fNmax {0}
int	m_fvNmax {0}
int	m_rfz_b [SizeRFZ] {}
int	m_rfz_t [SizeRFZ] {}
int	m_rfz_ib [SizeRFZ][SizeI] {}
int	m_rfz_it [SizeRFZ][SizeI] {}
int	m_rfzv_n [SizeRFZV] {}
int	m_rfzv_i [SizeRFZV][SizeIV] {}
float	m_dzdrmin {0.}
float	m_dzdrmax {0.}
float	m_ipt {0.}
float	m_iptmax {0.}
float	m_ipt2 {0.}
float	m_COF {0.}
float	m_sF {0.}
float	m_sFv {0.}
ServiceHandle< ITHistSvc >	m_thistSvc {this, "THistSvc", "THistSvc", "Histogramming svc"}
TTree *	m_outputTree {nullptr}
std::mutex	m_mutex
std::string	m_treeName {""}
std::string	m_treeFolder {"/valNtuples/"}
float m_d0	ATLAS_THREAD_SAFE = 0
float m_z0	ATLAS_THREAD_SAFE = 0
float m_pt	ATLAS_THREAD_SAFE = 0
float m_eta	ATLAS_THREAD_SAFE = 0
double m_x1	ATLAS_THREAD_SAFE = 0
double m_x2	ATLAS_THREAD_SAFE = 0
double m_x3	ATLAS_THREAD_SAFE = 0
double m_y1	ATLAS_THREAD_SAFE = 0
double m_y2	ATLAS_THREAD_SAFE = 0
double m_y3	ATLAS_THREAD_SAFE = 0
double m_z1	ATLAS_THREAD_SAFE = 0
double m_z2	ATLAS_THREAD_SAFE = 0
double m_z3	ATLAS_THREAD_SAFE = 0
double m_r1	ATLAS_THREAD_SAFE = 0
double m_r2	ATLAS_THREAD_SAFE = 0
double m_r3	ATLAS_THREAD_SAFE = 0
float m_quality	ATLAS_THREAD_SAFE = 0
int m_type	ATLAS_THREAD_SAFE = 0
double m_dzdr_t	ATLAS_THREAD_SAFE = 0
double m_dzdr_b	ATLAS_THREAD_SAFE = 0
double m_dr_b	ATLAS_THREAD_SAFE = 0
double m_dr_t	ATLAS_THREAD_SAFE = 0
bool m_givesTrack	ATLAS_THREAD_SAFE = 0
float m_trackPt	ATLAS_THREAD_SAFE = 0
float m_trackEta	ATLAS_THREAD_SAFE = 0
long m_eventNumber	ATLAS_THREAD_SAFE = 0

Print internal tool parameters and status
enum	Size { SizeRF =53 , SizeZ =11 , SizeRFZ =SizeRFSizeZ , SizeI =9 , SizeRFV =100 , SizeZV =3 , SizeRFZV =SizeRFVSizeZV , SizeIV =6 }
	enum for array sizes More...
virtual MsgStream &	dump (EventData &data, MsgStream &out) const override

Disallow default instantiation, copy, assignment
	SiSpacePointsSeedMaker_HeavyIon ()=delete
	SiSpacePointsSeedMaker_HeavyIon (const SiSpacePointsSeedMaker_HeavyIon &)=delete
SiSpacePointsSeedMaker_HeavyIon &	operator= (const SiSpacePointsSeedMaker_HeavyIon &)=delete
MsgStream &	dumpConditions (EventData &data, MsgStream &out) const
void	buildFrameWork ()
void	buildBeamFrameWork (const EventContext &ctx, 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	production2Sp (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	production3SpNoVertex (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
bool	newVertices (EventData &data, const std::list< Trk::Vertex > &) const
void	findNext (EventData &data) const
bool	isZCompatible (EventData &data, float &, float &, float &) const
void	initializeEventData (EventData &data) const
bool	isUsed (const Trk::SpacePoint *sp, const Trk::PRDtoTrackMap &prd_to_track_map) 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 float	dZVertexMin (EventData &data, float &)
static void	convertToBeamFrameWork (EventData &data, const Trk::SpacePoint *const &, std::array< float, 3 > &)

Detailed Description

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

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

Definition at line 54 of file SiSpacePointsSeedMaker_HeavyIon.h.

Member Enumeration Documentation

◆ Size

enum InDet::SiSpacePointsSeedMaker_HeavyIon::Size

private

enum for array sizes

Enumerator
SizeRF
SizeZ
SizeRFZ
SizeI
SizeRFV
SizeZV
SizeRFZV
SizeIV

Definition at line 127 of file SiSpacePointsSeedMaker_HeavyIon.h.

              {SizeRF=53,
               SizeZ=11,
               SizeRFZ=SizeRF*SizeZ,
               SizeI=9,
               SizeRFV=100,
               SizeZV=3,
               SizeRFZV=SizeRFV*SizeZV,
               SizeIV=6};

Constructor & Destructor Documentation

◆ SiSpacePointsSeedMaker_HeavyIon() [1/3]

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

Definition at line 30 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

  : base_class(t, n, p)
{
}

◆ ~SiSpacePointsSeedMaker_HeavyIon()

virtual InDet::SiSpacePointsSeedMaker_HeavyIon::~SiSpacePointsSeedMaker_HeavyIon ( )

virtualdefault

◆ SiSpacePointsSeedMaker_HeavyIon() [2/3]

InDet::SiSpacePointsSeedMaker_HeavyIon::SiSpacePointsSeedMaker_HeavyIon ( )

privatedelete

◆ SiSpacePointsSeedMaker_HeavyIon() [3/3]

InDet::SiSpacePointsSeedMaker_HeavyIon::SiSpacePointsSeedMaker_HeavyIon ( const SiSpacePointsSeedMaker_HeavyIon & )

privatedelete

Member Function Documentation

◆ buildBeamFrameWork()

void InDet::SiSpacePointsSeedMaker_HeavyIon::buildBeamFrameWork	(	const EventContext &	ctx,
		EventData &	data ) const

private

Definition at line 807 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle{m_beamSpotKey, ctx};
 
  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_HeavyIon::buildFrameWork ( )

private

Definition at line 670 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  m_ptmin = std::max( std::abs(m_ptmin), float(100.*m_fieldScale));
  m_etamax    = std::abs(m_etamax)                 ;
  m_dzdrmax   = 1.f/std::tan(2.f*std::atan(exp(-m_etamax)));
  m_dzdrmin   =-m_dzdrmax                      ;
  m_COF       =  134*.05*9.                    ;
  m_ipt       = 1.f/std::abs(.9f*m_ptmin)            ;
  m_iptmax  = m_ptmax == std::numeric_limits<float>::max()
    ? 0.f
    : 1.f/std::abs(1.1f*m_ptmax);
  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.;
  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 radius-azimuthal-Z sorted containers for Z-vertices
  //
  const int   NFtmax = SizeRFV;
  const float sFvmax = static_cast<float>(NFtmax)/pi2;
  m_sFv = m_ptmin/m_fieldScale /120.f;
  if (m_sFv > sFvmax) m_sFv = sFvmax; 
  m_fvNmax = static_cast<int>(pi2*m_sFv);
  if (m_fvNmax>=NFtmax) m_fvNmax = NFtmax-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;
      }
    }
  }
 
  // Build maps for radius-azimuthal-Z sorted collections for Z
  //
  for (int f=0; f<=m_fvNmax; ++f) {
    int fb = f-1; if (fb<0       ) fb=m_fvNmax;
    int ft = f+1; if (ft>m_fvNmax) ft=0;
    
    // For each azimuthal region loop through central Z regions
    //
    for (int z=0; z<SizeZV; ++z) {      
      int a  = f *SizeZV+z;
      int b  = fb*SizeZV+z;
      int c  = ft*SizeZV+z;
      m_rfzv_n[a]    = 3;
      m_rfzv_i[a][0] = a;
      m_rfzv_i[a][1] = b;
      m_rfzv_i[a][2] = c;
      if (z>1) {
    m_rfzv_n[a]    = 6;
    m_rfzv_i[a][3] = a-1;
    m_rfzv_i[a][4] = b-1;
    m_rfzv_i[a][5] = c-1;
      } else if (z<1) {
    m_rfzv_n[a]    = 6;
    m_rfzv_i[a][3] = a+1;
    m_rfzv_i[a][4] = b+1;
    m_rfzv_i[a][5] = c+1;
      }
    }
  }
}

◆ convertToBeamFrameWork()

void InDet::SiSpacePointsSeedMaker_HeavyIon::convertToBeamFrameWork	(	EventData &	data,
		const Trk::SpacePoint *const &	sp,
		std::array< float, 3 > &	r )

staticprivate

Definition at line 842 of file SiSpacePointsSeedMaker_HeavyIon.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_HeavyIon::dump	(	EventData &	data,
		MsgStream &	out ) const

overridevirtual

Definition at line 470 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

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

◆ dumpConditions()

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

private

Definition at line 482 of file SiSpacePointsSeedMaker_HeavyIon.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<<"| 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(2) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r1minv
     <<"                              |"<<endmsg;
  out<<"| min radius second SP(2) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r2minv
     <<"                              |"<<endmsg;
  out<<"| max radius first  SP(2) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r1maxv
     <<"                              |"<<endmsg;
  out<<"| max radius second SP(2) | "
     <<std::setw(12)<<std::setprecision(5)<<m_r2maxv
     <<"                              |"<<endmsg;
  out<<"| min space points dR (b)  | "
     <<std::setw(12)<<std::setprecision(5)<<m_drmin_b
     <<"                              |"<<endmsg;
  out<<"| min space points dR (t)     | "
     <<std::setw(12)<<std::setprecision(5)<<m_drmin_t
     <<"                              |"<<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<<"| max       impact pps    | "
     <<std::setw(12)<<std::setprecision(5)<<m_diverpps
     <<"                              |"<<endmsg;
  out<<"| max       impact sss    | "
     <<std::setw(12)<<std::setprecision(5)<<m_diversss
     <<"                              |"<<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_HeavyIon::dumpEvent	(	EventData &	data,
		MsgStream &	out )

staticprivate

Definition at line 608 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

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

◆ dZVertexMin()

float InDet::SiSpacePointsSeedMaker_HeavyIon::dZVertexMin	(	EventData &	data,
		float &	Z )

staticprivate

Definition at line 1560 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  float dZm = std::numeric_limits<float>::max();
  for (const float& v: data.l_vertex) {
    float dZ = std::abs(v-Z);
    if (dZ<dZm) dZm = dZ;
  }
  return dZm;
}

◆ erase()

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

staticprivate

Definition at line 914 of file SiSpacePointsSeedMaker_HeavyIon.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();
  }
 
  for (int i=0; i<data.nrfzv; ++i) {
    int n = data.rfzv_index[i];
    data.rfzv_map[n] = 0;
    data.rfzv_Sorted[n].clear();
  }
  data.state = 0;
  data.ns    = 0;
  data.nsaz  = 0;
  data.nsazv = 0;
  data.nr    = 0;
  data.nrfz  = 0;
  data.nrfzv = 0;
}

◆ fillLists()

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

private

Definition at line 854 of file SiSpacePointsSeedMaker_HeavyIon.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;
 
      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;
 
      if ((*r)->spacepoint->clusterList().second == 0 && z>=3 && z<=7) { 
    z<=4 ? z=0 : z>=6 ? z=2 : z=1;
    
    // Azimuthal angle and Z-coordinate sort for fast vertex search
    //
    f = static_cast<int>(F*m_sFv);
        if (f < 0) f += m_fvNmax;
        else if (f> m_fvNmax) f -= m_fvNmax;
    n = f*SizeZV+z;
        ++data.nsazv;
    data.rfzv_Sorted[n].push_back(*r);
        if (!data.rfzv_map[n]++) data.rfzv_index[data.nrfzv++] = n;
      }
      r = data.r_Sorted[i].erase(r);
    }
    data.r_map[i] = 0;
  }
  data.nr    = 0;
  data.state = 0;
}

◆ fillSeeds()

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

staticprivate

Definition at line 1641 of file SiSpacePointsSeedMaker_HeavyIon.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_HeavyIon::finalize ( )

overridevirtual

Definition at line 117 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  return AlgTool::finalize();
}

◆ find2Sp()

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

overridevirtual

with two space points with or without vertex constraint

Definition at line 360 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  data.izvertex = not lv.empty();
 
  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.fvNmin  = 0   ;
    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_HeavyIon::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 396 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  data.izvertex = not lv.empty();
 
  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.fvNmin  = 0               ;
    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_HeavyIon::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 427 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  find3Sp(ctx, data, lv);
}

◆ findNext()

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

private

Definition at line 633 of file SiSpacePointsSeedMaker_HeavyIon.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_HeavyIon::findVSp	(	const EventContext &	ctx,
		EventData &	data,
		const std::list< Trk::Vertex > &	lv ) const

overridevirtual

Definition at line 438 of file SiSpacePointsSeedMaker_HeavyIon.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.fvNmin  = 0               ;
    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_HeavyIon::getWriteNtupleBoolProperty ( ) const

overridevirtual

Definition at line 1716 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

                                                                           {
  return m_writeNtuple;
}

◆ initialize()

StatusCode InDet::SiSpacePointsSeedMaker_HeavyIon::initialize ( )

overridevirtual

Definition at line 40 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  StatusCode sc = AlgTool::initialize();
  ATH_CHECK( m_prdToTrackMap.initialize( !m_prdToTrackMap.key().empty()));
  
  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());
 
  // 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;
 
  if (m_writeNtuple) {
    ATH_CHECK( m_thistSvc.retrieve() );
 
    m_treeName = (std::string("SeedTree_")+name());
    std::replace( m_treeName.begin(), m_treeName.end(), '.', '_' );
 
    m_outputTree = new TTree( m_treeName.c_str() , "SeedMakerValTool"); 
 
    m_outputTree->Branch("eventNumber",    &m_eventNumber); 
    m_outputTree->Branch("d0",             &m_d0);
    m_outputTree->Branch("z0",             &m_z0);
    m_outputTree->Branch("pt",             &m_pt);
    m_outputTree->Branch("eta",            &m_eta);
    m_outputTree->Branch("x1",             &m_x1);
    m_outputTree->Branch("x2",             &m_x2);
    m_outputTree->Branch("x3",             &m_x3);
    m_outputTree->Branch("y1",             &m_y1);
    m_outputTree->Branch("y2",             &m_y2);
    m_outputTree->Branch("y3",             &m_y3);
    m_outputTree->Branch("z1",             &m_z1);
    m_outputTree->Branch("z2",             &m_z2);
    m_outputTree->Branch("z3",             &m_z3);
    m_outputTree->Branch("r1",             &m_r1);
    m_outputTree->Branch("r2",             &m_r2);
    m_outputTree->Branch("r3",             &m_r3);
    m_outputTree->Branch("quality",        &m_quality);
    m_outputTree->Branch("seedType",       &m_type);
    m_outputTree->Branch("givesTrack",     &m_givesTrack);
    m_outputTree->Branch("dzdr_b",         &m_dzdr_b);
    m_outputTree->Branch("dzdr_t",         &m_dzdr_t);
    m_outputTree->Branch("dr_b",           &m_dr_b);
    m_outputTree->Branch("dr_t",           &m_dr_t);
    m_outputTree->Branch("track_pt",       &m_trackPt);
    m_outputTree->Branch("track_eta",      &m_trackEta);
 
    TString fullTreeName = m_treeFolder + m_treeName;
    ATH_CHECK(  m_thistSvc->regTree( fullTreeName.Data(), m_outputTree )  );
  }
  
  return sc;
}

◆ initializeEventData()

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

private

Definition at line 1657 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

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

◆ isUsed()

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

inlineprivate

Definition at line 318 of file SiSpacePointsSeedMaker_HeavyIon.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_HeavyIon::isZCompatible	(	EventData &	data,
		float &	Zv,
		float &	R,
		float &	T ) const

private

Definition at line 1541 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  if (m_useVertexPosition) {
    if (Zv < data.zCollisionMinimum || Zv > data.zCollisionMaximum) return false;
  }
  else {
    if (Zv < m_zmin || Zv > m_zmax) return false;
  }
  if (!data.izvertex) 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)*sqrt(1.+T*T);
}

◆ newEvent()

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

overridevirtual

read the prd to track map, in case we want to use it.

usual rejection of SP used in previous track finding passes if we run with the PRT to track map + check of the max radius

SP creation, entry into list of the data object

Definition at line 126 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  data.trigger = false;
  if (!m_pixel && !m_sct) return;
  erase(data);
  buildBeamFrameWork(ctx, 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_HeavyIon: 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./m_r_rstep;
  int   irmax  = m_r_size-1  ;
  const Trk::PRDtoTrackMap *prd_to_track_map_cptr = nullptr;
  if (not m_prdToTrackMap.key().empty()) {
    SG::ReadHandle<Trk::PRDtoTrackMap> 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<m_r_rmin || r>=m_r_rmax) continue;
      if (prd_to_track_map_cptr and 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<m_r_rmin || r>=m_r_rmax) continue;
      if (prd_to_track_map_cptr and 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;
    }
      }
    }
  }
 
  if (m_useOverlap && !data.checketa) {
    SG::ReadHandle<SpacePointOverlapCollection> spacepointsOverlap{m_spacepointsOverlap, ctx};
    if (spacepointsOverlap.isValid()) {
      for (const Trk::SpacePoint* sp: *spacepointsOverlap) {
    float r = sp->r();
    if (r<m_r_rmin || 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);
    if (!sps) continue;
 
    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_HeavyIon::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 1492 of file SiSpacePointsSeedMaker_HeavyIon.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_HeavyIon::newRegion	(	const EventContext &	ctx,
		EventData &	data,
		const std::vector< IdentifierHash > &	vPixel,
		const std::vector< IdentifierHash > &	vSCT ) const

overridevirtual

Definition at line 249 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  if (not data.initialized) initializeEventData(data);
 
  data.trigger = false;
  if (!m_pixel && !m_sct) return;
  erase(data);
 
  buildBeamFrameWork(ctx, 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_HeavyIon: 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./m_r_rstep;
  int irmax = m_r_size-1;
 
  // 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;
      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;
      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_HeavyIon::newRegion	(	const EventContext &	ctx,
		EventData &	data,
		const std::vector< IdentifierHash > &	vPixel,
		const std::vector< IdentifierHash > &	vSCT,
		const IRoiDescriptor &	iRD ) const

overridevirtual

Definition at line 347 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

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

◆ newSeed() [1/2]

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

staticprivate

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

staticprivate

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

staticprivate

Definition at line 1574 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  InDet::SiSpacePointForSeed* sps = nullptr;
 
  std::array<float, 3> r{0, 0, 0};
  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_HeavyIon::newVertices	(	EventData &	data,
		const std::list< Trk::Vertex > &	lV ) const

private

Definition at line 650 of file SiSpacePointsSeedMaker_HeavyIon.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()));
    if (data.l_vertex.size() >= m_maxNumberVertices) break;
  }
  return false;
}

◆ next()

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

overridevirtual

Definition at line 1529 of file SiSpacePointsSeedMaker_HeavyIon.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_HeavyIon & InDet::SiSpacePointsSeedMaker_HeavyIon::operator= ( const SiSpacePointsSeedMaker_HeavyIon & )

privatedelete

◆ production2Sp()

void InDet::SiSpacePointsSeedMaker_HeavyIon::production2Sp ( EventData & data ) const

private

Definition at line 946 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  if (data.nsazv<2) return;
 
  std::vector<InDet::SiSpacePointForSeed*>::iterator r0,r0e,r,re;
  int nseed = 0;
 
  // Loop through all azimuthal regions
  //
  for (int f=data.fvNmin; f<=m_fvNmax; ++f) {
 
    // For each azimuthal region loop through Z regions
    //
    int z = 0;
    if (!data.endlist) z = data.zMin;
    for (; z<SizeZV; ++z) {
      int a  = f*SizeZV+z;
      if (!data.rfzv_map[a]) continue;
      r0  = data.rfzv_Sorted[a].begin();
      r0e = data.rfzv_Sorted[a].end  ();
 
      if (!data.endlist) {
        r0 = data.rMin;
        data.endlist = true;
      }
 
      // Loop through trigger space points
      //
      for (; r0!=r0e; ++r0) {
    float X  = (*r0)->x();
    float Y  = (*r0)->y();
    float R  = (*r0)->radius();
    if (R<m_r2minv) continue;
        if (R>m_r2maxv) break;
    float Z  = (*r0)->z();
    float ax = X/R;
    float ay = Y/R;
 
    // Bottom links production
    //
    int NB = m_rfzv_n[a];
    for (int i=0; i<NB; ++i) {
      int an = m_rfzv_i[a][i];
      if (!data.rfzv_map[an]) continue;
 
      r  =  data.rfzv_Sorted[an].begin();
      re =  data.rfzv_Sorted[an].end  ();
      
      for (; r!=re; ++r) {
        float Rb =(*r)->radius();
        if (Rb<m_r1minv) continue;
            if (Rb>m_r1maxv) break;
        float dR = R-Rb;
        if (dR<m_drminv) break;
            if (dR>m_drmax) continue;
        float dZ = Z-(*r)->z();
        float Tz = dZ/dR;
            if (Tz<m_dzdrmin || Tz>m_dzdrmax) continue;
        float Zo = Z-R*Tz;
 
        // Comparison with vertices Z coordinates
        //
        if (!isZCompatible(data, Zo, Rb, Tz)) continue;
 
        // Momentum cut
        //
        float dx =(*r)->x()-X;
        float dy =(*r)->y()-Y;
        float x  = dx*ax+dy*ay          ;
        float y  =-dx*ay+dy*ax          ;
        float xy = x*x+y*y              ; if (xy == 0.) continue;
        float r2 = 1./xy                ;
        float Ut = x*r2                 ;
        float Vt = y*r2                 ;
        float UR = Ut*R+1.              ; if (UR == 0.) continue;
        float A  = Vt*R/UR              ;
        float B  = Vt-A*Ut              ;
        // min pt cut
        if (std::abs(B*data.K) > m_ipt*std::sqrt(1.f+A*A)) continue;
        // max pt cut
        if (std::abs(B * data.K) < m_iptmax * std::sqrt(1.f + A * A)) continue;
            ++nseed;
        newSeed(data, (*r)->spacepoint, (*r0)->spacepoint,Zo);
      }
    }
    if (nseed < m_maxsize) continue;
    data.endlist=false;
        data.rMin = (++r0);
        data.fvNmin=f;
        data.zMin=z;
    return;
      }
    }
  }
  data.endlist = true;
}

◆ production3Sp() [1/2]

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

private

Definition at line 1047 of file SiSpacePointsSeedMaker_HeavyIon.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 through 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();
      } 
      if (data.izvertex) {
    production3Sp(data, rb, rbe, rt, rte, NB, NT, nseed);
      } else {
    production3SpNoVertex(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_HeavyIon::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 1098 of file SiSpacePointsSeedMaker_HeavyIon.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_r2minv) continue;
    if (R > m_r2maxv) break;
 
    const Trk::SpacePoint* SP0 = (*r0)->spacepoint;
    if (SP0->clusterList().second) break;
 
    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_r1minv) {
      rb[i]=r;
      continue;
    }
    // if (Rb > m_r1maxv) break;
    float dR = R-Rb;
    if (dR > m_drmax) {rb[i]=r; continue;}
    if (dR < m_drmin_b) break;
    if ((*r)->sur()==sur0) continue;
 
    float Tz = (Z-(*r)->z())/dR;
 
    if (Tz < m_dzdrmin || Tz > m_dzdrmax) continue;
    
    // Comparison with vertices Z coordinates
    //
    float Zo = Z-R*Tz;
        if (!isZCompatible(data, Zo, Rb, Tz)) continue;
    data.SP[Nb] = (*r);
    if (m_writeNtuple) data.SP[Nb]->setDZDR(Tz);
        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();
    if (Rt < m_r3minv) {
      rt[i]=r;
      continue;
    }
    float dR = Rt-R;
        if (dR<m_drmin_t) {
          rt[i]=r;
          continue;
        }
    if (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
    //
    float Zo = Z-R*Tz;
        if (!isZCompatible(data, Zo, R ,Tz)) continue;
    data.SP[Nt] = (*r);
    if (m_writeNtuple) data.SP[Nt]->setDZDR(Tz);
        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 ipt2  = m_ipt2    ;
    float ipt2max = m_iptmax * m_iptmax;
    float K     = data.K       ;
    float K2    = K*K       ;
    float COF   = m_COF     ;
    float ipt2K = ipt2/K2   ;
    float ipt2Kmax = ipt2max/K2;
    float ipt2C = ipt2*COF  ;
    float COFK  = COF*K2    ;
    covr0      *= 2.        ;
    covz0      *= 2.        ;
 
    //
    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   ;
      float dZ    = dZVertexMin(data, Zob);
      float Iz    = (dZ*dZ)/Tzb2 ;
 
      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 || B2 < ipt2Kmax*S2 || dT*S2 > B2*CSA) continue;
    float Im  = std::abs((A-B*R)*R)                  ;
 
    if ( Im > imc ) continue;
    Im = Im*Im+Iz;
    newOneSeed(data, SPb, SP0, data.SP[t]->spacepoint, Zob, Im);
 
    if (m_writeNtuple) {
      float theta = std::atan(1.f / Ts);
      data.SP[t]->setEta(-std::log(std::tan(0.5f * theta)));
      data.SP[t]->setPt(std::sqrt(S2 / B2) / (1000.f * data.K)); 
    }
      }
    }
    nseed += data.mapOneSeeds.size();
    fillSeeds(data);
    if (nseed>=m_maxsize) {
      data.endlist=false;
      ++r0;
      data.rMin = r0;
      return;
    } 
  }
}

◆ production3SpNoVertex()

void InDet::SiSpacePointsSeedMaker_HeavyIon::production3SpNoVertex	(	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 1291 of file SiSpacePointsSeedMaker_HeavyIon.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_r2minv) continue;
    if (R > m_r2maxv) 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_r1minv) {
          rb[i]=r;
          continue;
        }
        // if (Rb > m_r1maxv) break;
    float dR = R-Rb;
    if (dR > m_drmax) {
          rb[i]=r;
          continue;
        }
    if (dR < m_drmin_b) break;
    if ((*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
    //
    float Zo = Z-R*Tz;
        if (!isZCompatible(data, Zo, Rb, Tz)) continue;
    data.SP[Nb] = (*r);
    if (m_writeNtuple) data.SP[Nb]->setDZDR(Tz);
        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();
    if (Rt < m_r3minv) {
          rt[i]=r;
          continue;
        }
    float dR = Rt-R;
        if (dR<m_drmin_t) {
          rt[i]=r;
          continue;
        }
    if (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
    //
    float Zo = Z-R*Tz;
        if (!isZCompatible(data, Zo, R, Tz)) continue;
    data.SP[Nt] = (*r);
    if (m_writeNtuple) data.SP[Nt]->setDZDR(Tz);
        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 ipt2max = m_iptmax * m_iptmax;
    float K     = data.K       ;
    float K2    = K*K       ;
    float COF   = m_COF     ;
    float ipt2K = ipt2/K2   ;
    float ipt2Kmax = ipt2max/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 || B2 < ipt2Kmax*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;
    } else if (Im > m_diversss) {
          continue;
        }
    newOneSeed(data, SPb, SP0, data.SP[t]->spacepoint, Zob, Im);
 
    if (m_writeNtuple) {
      float theta = std::atan(1.f / Ts);
      data.SP[t]->setEta(-std::log(std::tan(0.5f * theta)));
      data.SP[t]->setPt(std::sqrt(S2 / B2) / (1000.f * data.K));
        }
      }
    }
    nseed += data.mapOneSeeds.size();
    fillSeeds(data);
    if (nseed>=m_maxsize) {
      data.endlist=false;
      ++r0;
      data.rMin = r0;
      return;
    } 
  }
}

◆ writeNtuple()

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

overridevirtual

Definition at line 1669 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  if (not m_writeNtuple) return;
 
  std::lock_guard<std::mutex> lock(m_mutex);
  
  if (track != nullptr) {
    m_trackPt = track->trackParameters()->front()->pT() / 1000.f;
    m_trackEta = std::abs(track->trackParameters()->front()->eta());
  } else {
    m_trackPt = -1.;
    m_trackEta = -1.; 
  }
 
  const std::vector<const Trk::SpacePoint*>& sps = seed->spacePoints();
  
  m_d0           =   seed->d0();
  m_z0           =   seed->zVertex();
  m_eta          =   sps[0]->globalPosition().eta();
  m_x1           =   sps[0]->globalPosition().x();
  m_x2           =   sps[1]->globalPosition().x();
  m_x3           =   sps[2]->globalPosition().x();
  m_y1           =   sps[0]->globalPosition().y();
  m_y2           =   sps[1]->globalPosition().y();
  m_y3           =   sps[2]->globalPosition().y();
  m_z1           =   sps[0]->globalPosition().z();
  m_z2           =   sps[1]->globalPosition().z();
  m_z3           =   sps[2]->globalPosition().z();
  m_r1           =   sps[0]->globalPosition().perp();
  m_r2           =   sps[1]->globalPosition().perp();
  m_r3           =   sps[2]->globalPosition().perp();
  m_type         =   seedType;
  m_dzdr_b       =   ( sps[1]->globalPosition().z() -  sps[0]->globalPosition().z() )/( sps[1]->globalPosition().perp() - sps[0]->globalPosition().perp());
  m_dzdr_t       =   ( sps[2]->globalPosition().z() - sps[1]->globalPosition().z()) / (sps[2]->globalPosition().perp() - sps[1]->globalPosition().perp());
  m_dr_b         =   sps[1]->globalPosition().perp() - sps[0]->globalPosition().perp();
  m_dr_t         =   sps[2]->globalPosition().perp() - sps[1]->globalPosition().perp();
  m_pt           =   seed->pt();
  m_givesTrack   =   !(track == nullptr);
  m_eventNumber  =   eventNumber;
  
  TTree* outputTree ATLAS_THREAD_SAFE = m_outputTree;
  outputTree->Fill(); 
}

Member Data Documentation

◆ ATLAS_THREAD_SAFE [1/26]

long m_eventNumber InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 251 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [2/26]

float m_trackEta InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 250 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [3/26]

float m_trackPt InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 249 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [4/26]

bool m_givesTrack InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 248 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [5/26]

double m_dr_t InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 247 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [6/26]

double m_dr_b InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 246 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [7/26]

double m_dzdr_b InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 245 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [8/26]

double m_dzdr_t InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 244 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [9/26]

int m_type InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 243 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [10/26]

float m_quality InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 242 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [11/26]

double m_r3 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 241 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [12/26]

double m_r2 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 240 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [13/26]

double m_r1 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 239 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [14/26]

double m_z3 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 238 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [15/26]

double m_z2 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 237 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [16/26]

double m_z1 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 236 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [17/26]

double m_y3 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 235 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [18/26]

double m_y2 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 234 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [19/26]

double m_y1 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 233 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [20/26]

double m_x3 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 232 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [21/26]

double m_x2 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 231 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [22/26]

double m_x1 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 230 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [23/26]

float m_eta InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 229 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [24/26]

float m_pt InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 228 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [25/26]

float m_z0 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 227 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ ATLAS_THREAD_SAFE [26/26]

float m_d0 InDet::SiSpacePointsSeedMaker_HeavyIon::ATLAS_THREAD_SAFE = 0

mutableprivate

Definition at line 226 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ m_beamSpotKey

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

private

Definition at line 145 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_COF

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

private

Definition at line 213 of file SiSpacePointsSeedMaker_HeavyIon.h.

213{0.};

◆ m_diver

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

private

Definition at line 175 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_diverpps

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

private

Definition at line 176 of file SiSpacePointsSeedMaker_HeavyIon.h.

176{this, "maxdImpactPPS", 1.7};

◆ m_diversss

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_diversss {this, "maxdImpactSSS", 1000.}

private

Definition at line 177 of file SiSpacePointsSeedMaker_HeavyIon.h.

177{this, "maxdImpactSSS", 1000.};

◆ m_drmax

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

private

Definition at line 167 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_drmin_b

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_drmin_b {this, "mindRadius", 10.}

private

Definition at line 165 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_drmin_t

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_drmin_t {this, "mindRadiusTop", 10.}

private

Definition at line 166 of file SiSpacePointsSeedMaker_HeavyIon.h.

166{this, "mindRadiusTop", 10.};

◆ m_drminv

float InDet::SiSpacePointsSeedMaker_HeavyIon::m_drminv {20.}

private

Definition at line 192 of file SiSpacePointsSeedMaker_HeavyIon.h.

192{20.};

◆ m_dzdrmax

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

private

Definition at line 209 of file SiSpacePointsSeedMaker_HeavyIon.h.

209{0.};

◆ m_dzdrmin

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

private

Definition at line 208 of file SiSpacePointsSeedMaker_HeavyIon.h.

208{0.};

◆ m_dzdrver

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

private

Definition at line 174 of file SiSpacePointsSeedMaker_HeavyIon.h.

174{this, "maxdZdRver", 0.02};

◆ m_dzver

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

private

Definition at line 173 of file SiSpacePointsSeedMaker_HeavyIon.h.

173{this, "maxdZver", 5.};

◆ m_etamax

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

private

Definition at line 184 of file SiSpacePointsSeedMaker_HeavyIon.h.

184{this, "etaMax", 2.7};

◆ m_fieldCondObjInputKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> InDet::SiSpacePointsSeedMaker_HeavyIon::m_fieldCondObjInputKey {this, "AtlasFieldCacheCondObj", "fieldCondObj", "Name of the Magnetic Field conditions object key"}

private

Definition at line 147 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_fieldScale

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

private

Definition at line 187 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_fNmax

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_fNmax {0}

private

Definition at line 200 of file SiSpacePointsSeedMaker_HeavyIon.h.

200{0};

◆ m_fvNmax

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_fvNmax {0}

private

Definition at line 201 of file SiSpacePointsSeedMaker_HeavyIon.h.

201{0};

◆ m_initialized

bool InDet::SiSpacePointsSeedMaker_HeavyIon::m_initialized {false}

private

Definition at line 197 of file SiSpacePointsSeedMaker_HeavyIon.h.

197{false};

◆ m_ipt

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

private

Definition at line 210 of file SiSpacePointsSeedMaker_HeavyIon.h.

210{0.};

◆ m_ipt2

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

private

Definition at line 212 of file SiSpacePointsSeedMaker_HeavyIon.h.

212{0.};

◆ m_iptmax

float InDet::SiSpacePointsSeedMaker_HeavyIon::m_iptmax {0.}

private

Definition at line 211 of file SiSpacePointsSeedMaker_HeavyIon.h.

211{0.};

◆ m_maxNumberVertices

UnsignedIntegerProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_maxNumberVertices {this, "maxNumberVertices", 1}

private

Definition at line 159 of file SiSpacePointsSeedMaker_HeavyIon.h.

159{this, "maxNumberVertices", 1};

◆ m_maxOneSize

IntegerProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_maxOneSize {this, "maxSeedsForSpacePoint", 6}

private

Definition at line 158 of file SiSpacePointsSeedMaker_HeavyIon.h.

158{this, "maxSeedsForSpacePoint", 6};

◆ m_maxsize

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

private

Definition at line 156 of file SiSpacePointsSeedMaker_HeavyIon.h.

156{this, "maxSize", 20000};

◆ m_maxsizeSP

IntegerProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_maxsizeSP {this, "maxSizeSP", 4000}

private

Definition at line 157 of file SiSpacePointsSeedMaker_HeavyIon.h.

157{this, "maxSizeSP", 4000};

◆ m_mutex

std::mutex InDet::SiSpacePointsSeedMaker_HeavyIon::m_mutex

mutableprivate

Definition at line 221 of file SiSpacePointsSeedMaker_HeavyIon.h.

◆ m_outputlevel

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_outputlevel {0}

private

Definition at line 198 of file SiSpacePointsSeedMaker_HeavyIon.h.

198{0};

◆ m_outputTree

TTree* InDet::SiSpacePointsSeedMaker_HeavyIon::m_outputTree {nullptr}

private

Definition at line 219 of file SiSpacePointsSeedMaker_HeavyIon.h.

219{nullptr};

◆ m_pixel

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

private

Definition at line 153 of file SiSpacePointsSeedMaker_HeavyIon.h.

153{this, "usePixel", true};

◆ m_prdToTrackMap

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

private

Definition at line 148 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_ptmax

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_ptmax {this, "pTmax", std::numeric_limits<float>::max()}

private

Definition at line 186 of file SiSpacePointsSeedMaker_HeavyIon.h.

186{this, "pTmax", std::numeric_limits<float>::max()};

◆ m_ptmin

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_ptmin {this, "pTmin", 500.}

private

Definition at line 185 of file SiSpacePointsSeedMaker_HeavyIon.h.

185{this, "pTmin", 500.};

◆ m_r1maxv

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_r1maxv {this, "maxVRadius1", 60.}

private

Definition at line 161 of file SiSpacePointsSeedMaker_HeavyIon.h.

161{this, "maxVRadius1", 60.};

◆ m_r1minv

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_r1minv {this, "minVRadius1", 0.}

private

Definition at line 160 of file SiSpacePointsSeedMaker_HeavyIon.h.

160{this, "minVRadius1", 0.};

◆ m_r2maxv

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_r2maxv {this, "maxVRadius2", 200.}

private

Definition at line 163 of file SiSpacePointsSeedMaker_HeavyIon.h.

163{this, "maxVRadius2", 200.};

◆ m_r2minv

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_r2minv {this, "minVRadius2", 70.}

private

Definition at line 162 of file SiSpacePointsSeedMaker_HeavyIon.h.

162{this, "minVRadius2", 70.};

◆ m_r3minv

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_r3minv {this, "minVRadius3", 0.}

private

Definition at line 164 of file SiSpacePointsSeedMaker_HeavyIon.h.

164{this, "minVRadius3", 0.};

◆ m_r_rmax

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

private

Definition at line 171 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_r_rmin

FloatProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_r_rmin {this, "radMin", 43.}

private

Definition at line 170 of file SiSpacePointsSeedMaker_HeavyIon.h.

170{this, "radMin", 43.};

◆ m_r_rstep

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

private

Definition at line 172 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_r_size

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_r_size {0}

private

Definition at line 199 of file SiSpacePointsSeedMaker_HeavyIon.h.

199{0};

◆ m_rfz_b

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

private

Definition at line 202 of file SiSpacePointsSeedMaker_HeavyIon.h.

202{};

◆ m_rfz_ib

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

private

Definition at line 204 of file SiSpacePointsSeedMaker_HeavyIon.h.

204{};

◆ m_rfz_it

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

private

Definition at line 205 of file SiSpacePointsSeedMaker_HeavyIon.h.

205{};

◆ m_rfz_t

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

private

Definition at line 203 of file SiSpacePointsSeedMaker_HeavyIon.h.

203{};

◆ m_rfzv_i

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_rfzv_i[SizeRFZV][SizeIV] {}

private

Definition at line 207 of file SiSpacePointsSeedMaker_HeavyIon.h.

207{};

◆ m_rfzv_n

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_rfzv_n[SizeRFZV] {}

private

Definition at line 206 of file SiSpacePointsSeedMaker_HeavyIon.h.

206{};

◆ m_sct

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

private

Definition at line 154 of file SiSpacePointsSeedMaker_HeavyIon.h.

154{this, "useSCT", true};

◆ m_sF

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

private

Definition at line 214 of file SiSpacePointsSeedMaker_HeavyIon.h.

214{0.};

◆ m_sFv

float InDet::SiSpacePointsSeedMaker_HeavyIon::m_sFv {0.}

private

Definition at line 215 of file SiSpacePointsSeedMaker_HeavyIon.h.

215{0.};

◆ m_spacepointsOverlap

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

private

Definition at line 144 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_spacepointsPixel

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

private

Definition at line 143 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_spacepointsSCT

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

private

Definition at line 142 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_thistSvc

ServiceHandle<ITHistSvc> InDet::SiSpacePointsSeedMaker_HeavyIon::m_thistSvc {this, "THistSvc", "THistSvc", "Histogramming svc"}

private

Definition at line 218 of file SiSpacePointsSeedMaker_HeavyIon.h.

218{this, "THistSvc", "THistSvc", "Histogramming svc"};

◆ m_treeFolder

std::string InDet::SiSpacePointsSeedMaker_HeavyIon::m_treeFolder {"/valNtuples/"}

private

Definition at line 224 of file SiSpacePointsSeedMaker_HeavyIon.h.

224{"/valNtuples/"};

◆ m_treeName

std::string InDet::SiSpacePointsSeedMaker_HeavyIon::m_treeName {""}

private

Definition at line 223 of file SiSpacePointsSeedMaker_HeavyIon.h.

223{""};

◆ m_useOverlap

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

private

Definition at line 155 of file SiSpacePointsSeedMaker_HeavyIon.h.

155{this, "useOverlapSpCollection", false};

◆ m_useVertexPosition

BooleanProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_useVertexPosition {this, "useVertexPosition", false}

private

Definition at line 178 of file SiSpacePointsSeedMaker_HeavyIon.h.

178{this, "useVertexPosition", false};

◆ m_writeNtuple

BooleanProperty InDet::SiSpacePointsSeedMaker_HeavyIon::m_writeNtuple {this, "WriteNtuple", false, "Flag to write Validation Ntuples"}

private

Definition at line 179 of file SiSpacePointsSeedMaker_HeavyIon.h.

179{this, "WriteNtuple", false, "Flag to write Validation Ntuples"};

◆ m_zmax

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

private

Definition at line 169 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

◆ m_zmin

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

private

Definition at line 168 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

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_HeavyIon() [1/3]

◆ ~SiSpacePointsSeedMaker_HeavyIon()

◆ SiSpacePointsSeedMaker_HeavyIon() [2/3]

◆ SiSpacePointsSeedMaker_HeavyIon() [3/3]

Member Function Documentation

◆ buildBeamFrameWork()

◆ buildFrameWork()

◆ convertToBeamFrameWork()

◆ dump()

◆ dumpConditions()

◆ dumpEvent()

◆ dZVertexMin()

◆ 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()

◆ newVertices()

◆ next()

◆ operator=()

◆ production2Sp()

◆ production3Sp() [1/2]

◆ production3Sp() [2/2]

◆ production3SpNoVertex()

◆ writeNtuple()

Member Data Documentation

◆ ATLAS_THREAD_SAFE [1/26]

◆ ATLAS_THREAD_SAFE [2/26]

◆ ATLAS_THREAD_SAFE [3/26]

◆ ATLAS_THREAD_SAFE [4/26]

◆ ATLAS_THREAD_SAFE [5/26]

◆ ATLAS_THREAD_SAFE [6/26]

◆ ATLAS_THREAD_SAFE [7/26]

◆ ATLAS_THREAD_SAFE [8/26]

◆ ATLAS_THREAD_SAFE [9/26]

◆ ATLAS_THREAD_SAFE [10/26]

◆ ATLAS_THREAD_SAFE [11/26]

◆ ATLAS_THREAD_SAFE [12/26]

◆ ATLAS_THREAD_SAFE [13/26]

◆ ATLAS_THREAD_SAFE [14/26]

◆ ATLAS_THREAD_SAFE [15/26]

◆ ATLAS_THREAD_SAFE [16/26]

◆ ATLAS_THREAD_SAFE [17/26]

◆ ATLAS_THREAD_SAFE [18/26]

◆ ATLAS_THREAD_SAFE [19/26]

◆ ATLAS_THREAD_SAFE [20/26]

◆ ATLAS_THREAD_SAFE [21/26]

◆ ATLAS_THREAD_SAFE [22/26]

◆ ATLAS_THREAD_SAFE [23/26]

◆ ATLAS_THREAD_SAFE [24/26]

◆ ATLAS_THREAD_SAFE [25/26]

◆ ATLAS_THREAD_SAFE [26/26]

◆ m_beamSpotKey

◆ m_COF

◆ m_diver

◆ m_diverpps