InDet::SiSpacePointsSeedMaker_HeavyIon Class Reference

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
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_drmin {this, "mindRadius", 10.}
FloatProperty	m_drmax {this, "maxdRadius", 270.}
FloatProperty	m_zmin {this, "minZ", -250.}
FloatProperty	m_zmax {this, "maxZ", +250.}
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.}
Properties, which can be updated in initialize
FloatProperty	m_etamax {this, "etaMax", 2.7}
FloatProperty	m_ptmin {this, "pTmin", 500.}
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_ipt2 {0.}
float	m_COF {0.}
float	m_sF {0.}
float	m_sFv {0.}

Print internal tool parameters and status
enum	Size {   SizeRF =53 , SizeZ =11 , SizeRFZ =SizeRF*SizeZ , SizeI =9 ,   SizeRFV =100 , SizeZV =3 , SizeRFZV =SizeRFV*SizeZV , 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
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 &, float *)

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 52 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 125 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 26 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 718 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 584 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_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, float * r )

staticprivate

Definition at line 753 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 390 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 402 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<<"| |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     | "
     <<std::setw(12)<<std::setprecision(5)<<m_drmin
     <<"                              |"<<endmsg;
  out<<"| max space points dR     | "
     <<std::setw(12)<<std::setprecision(5)<<m_drmax
     <<"                              |"<<endmsg;
  out<<"| max dZ    impact        | "
     <<std::setw(12)<<std::setprecision(5)<<m_dzver 
     <<"                              |"<<endmsg;
  out<<"| max dZ/dR impact        | "
     <<std::setw(12)<<std::setprecision(5)<<m_dzdrver 
     <<"                              |"<<endmsg;
  out<<"| max       impact        | "
     <<std::setw(12)<<std::setprecision(5)<<m_diver
     <<"                              |"<<endmsg;
  out<<"| 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 522 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 1421 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 825 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 765 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 1502 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 73 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 280 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 316 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 347 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  find3Sp(ctx, data, lv);
}

findNext()

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

private

Definition at line 547 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 358 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 1533 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

                                                                           {
    return false;
}

initialize()

StatusCode InDet::SiSpacePointsSeedMaker_HeavyIon::initialize ( )

overridevirtual

Definition at line 36 of file SiSpacePointsSeedMaker_HeavyIon.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());
 
  // 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_HeavyIon::initializeEventData ( EventData & data ) const

private

Definition at line 1518 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
}

isZCompatible()

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

private

Definition at line 1407 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  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

Definition at line 82 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) {
 
    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 < 43. || 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) {
 
    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;
      InDet::SiSpacePointForSeed* sps = newSpacePoint(data, sp);
      int ir = static_cast<int>(sps->radius()*irstep);
          if (ir>irmax) ir = irmax;
      data.r_Sorted[ir].push_back(sps);
          ++data.r_map[ir];
      if (data.r_map[ir]==1) data.r_index[data.nr++] = ir;
      ++data.ns;
    }
      }
    }
  }
  fillLists(data);
}

newOneSeed()

void InDet::SiSpacePointsSeedMaker_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 1358 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 169 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 267 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 1459 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 1480 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 1435 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

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

newVertices()

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

private

Definition at line 564 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 1395 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 857 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

{
  if (data.nsazv<2) return;
 
  std::vector<InDet::SiSpacePointForSeed*>::iterator r0,r0e,r,re;
  int nseed = 0;
 
  // Loop thorugh 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              ;
        if (std::abs(B*data.K) > m_ipt*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 955 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 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();
      } 
      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 1006 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();
 
    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();
    float dR = R-Rb;
    if (dR > m_drmax) {rb[i]=r; continue;}
    if (dR < m_drmin) 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 (++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[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 (++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 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.        ;
 
    //
    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 || 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);
      }
    }
    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 1178 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();
 
    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();
    float dR = R-Rb;
    if (dR > m_drmax) {
          rb[i]=r;
          continue;
        }
    if (dR < m_drmin) 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 (++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[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 (++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;
    } else if (Im > m_diversss) {
          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_HeavyIon::writeNtuple ( const SiSpacePointsSeed * seed, const Trk::Track * track, int seedType, long eventNumber ) const

overridevirtual

Definition at line 1530 of file SiSpacePointsSeedMaker_HeavyIon.cxx.

                                                                                                                {
}

Member Data Documentation

m_beamSpotKey

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

private

Definition at line 143 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

m_COF

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

private

Definition at line 204 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0.};

m_diver

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

private

Definition at line 170 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxdImpact", 10.};

m_diverpps

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

private

Definition at line 171 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxdImpactPPS", 1.7};

m_diversss

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

private

Definition at line 172 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxdImpactSSS", 1000.};

m_drmax

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

private

Definition at line 163 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxdRadius", 270.};

m_drmin

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

private

Definition at line 162 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "mindRadius", 10.};

m_drminv

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

private

Definition at line 184 of file SiSpacePointsSeedMaker_HeavyIon.h.

{20.};

m_dzdrmax

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

private

Definition at line 201 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0.};

m_dzdrmin

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

private

Definition at line 200 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0.};

m_dzdrver

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

private

Definition at line 169 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxdZdRver", 0.02};

m_dzver

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

private

Definition at line 168 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxdZver", 5.};

m_etamax

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

private

Definition at line 177 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "etaMax", 2.7};

m_fieldCondObjInputKey

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

private

Initial value:

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

Definition at line 145 of file SiSpacePointsSeedMaker_HeavyIon.h.

                                                                       {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 179 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "fieldScale", 1.};

m_fNmax

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_fNmax {0}

private

Definition at line 192 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0};

m_fvNmax

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_fvNmax {0}

private

Definition at line 193 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0};

m_initialized

bool InDet::SiSpacePointsSeedMaker_HeavyIon::m_initialized {false}

private

Definition at line 189 of file SiSpacePointsSeedMaker_HeavyIon.h.

{false};

m_ipt

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

private

Definition at line 202 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0.};

m_ipt2

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

private

Definition at line 203 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0.};

m_maxNumberVertices

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

private

Definition at line 157 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxNumberVertices", 1};

m_maxOneSize

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

private

Definition at line 156 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxSeedsForSpacePoint", 6};

m_maxsize

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

private

Definition at line 154 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxSize", 20000};

m_maxsizeSP

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

private

Definition at line 155 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxSizeSP", 4000};

m_outputlevel

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_outputlevel {0}

private

Definition at line 190 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0};

m_pixel

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

private

Definition at line 151 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "usePixel", true};

m_ptmin

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

private

Definition at line 178 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "pTmin", 500.};

m_r1maxv

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

private

Definition at line 159 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxVRadius1", 60.};

m_r1minv

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

private

Definition at line 158 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "minVRadius1", 0.};

m_r2maxv

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

private

Definition at line 161 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "maxVRadius2", 200.};

m_r2minv

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

private

Definition at line 160 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "minVRadius2", 70.};

m_r_rmax

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

private

Definition at line 166 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "radMax", 600.};

m_r_rstep

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

private

Definition at line 167 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "radStep", 2.};

m_r_size

int InDet::SiSpacePointsSeedMaker_HeavyIon::m_r_size {0}

private

Definition at line 191 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0};

m_rfz_b

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

private

Definition at line 194 of file SiSpacePointsSeedMaker_HeavyIon.h.

{};

m_rfz_ib

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

private

Definition at line 196 of file SiSpacePointsSeedMaker_HeavyIon.h.

{};

m_rfz_it

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

private

Definition at line 197 of file SiSpacePointsSeedMaker_HeavyIon.h.

{};

m_rfz_t

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

private

Definition at line 195 of file SiSpacePointsSeedMaker_HeavyIon.h.

{};

m_rfzv_i

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

private

Definition at line 199 of file SiSpacePointsSeedMaker_HeavyIon.h.

{};

m_rfzv_n

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

private

Definition at line 198 of file SiSpacePointsSeedMaker_HeavyIon.h.

{};

m_sct

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

private

Definition at line 152 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "useSCT", true};

m_sF

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

private

Definition at line 205 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0.};

m_sFv

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

private

Definition at line 206 of file SiSpacePointsSeedMaker_HeavyIon.h.

{0.};

m_spacepointsOverlap

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

private

Definition at line 142 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

m_spacepointsPixel

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

private

Definition at line 141 of file SiSpacePointsSeedMaker_HeavyIon.h.

{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 140 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

m_useOverlap

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

private

Definition at line 153 of file SiSpacePointsSeedMaker_HeavyIon.h.

{this, "useOverlapSpCollection", false};

m_zmax

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

private

Definition at line 165 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

m_zmin

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

private

Definition at line 164 of file SiSpacePointsSeedMaker_HeavyIon.h.

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

---

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

• SiSpacePointsSeedMaker_HeavyIon.h
• SiSpacePointsSeedMaker_HeavyIon.cxx