InDet::SiZvertexMaker_xk Class Reference

Class for primary vertices z-coordinates generation (xKalman algorithm). More...

#include <SiZvertexMaker_xk.h>

Inheritance diagram for InDet::SiZvertexMaker_xk:

Collaboration diagram for InDet::SiZvertexMaker_xk:

Public Member Functions
Standard tool methods
	SiZvertexMaker_xk (const std::string &, const std::string &, const IInterface *)
virtual	~SiZvertexMaker_xk ()=default
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
Methods to initialize tool for new event or region
virtual std::list< Trk::Vertex >	newEvent (const EventContext &ctx, SiSpacePointsSeedMakerEventData &data) const override
virtual std::list< Trk::Vertex >	newRegion (const EventContext &ctx, SiSpacePointsSeedMakerEventData &data, const std::vector< IdentifierHash > &, const std::vector< IdentifierHash > &) const override
virtual std::list< Trk::Vertex >	newRegion (const EventContext &ctx, SiSpacePointsSeedMakerEventData &data, const std::vector< IdentifierHash > &, const std::vector< IdentifierHash > &, const IRoiDescriptor &) const override
Print internal tool parameters and status
virtual MsgStream &	dump (MsgStream &out) const override

Protected Attributes
Tool handles
ToolHandle< InDet::ISiSpacePointsSeedMaker >	m_seedsgenerator {this, "SeedMakerTool", "InDet::SiSpacePointsSeedMaker_ATLxk"}

Properties
IntegerProperty	m_nspoint {this, "SeedSize", 2}
IntegerProperty	m_histsize {this, "HistSize", 500}
IntegerProperty	m_mincontent {this, "minContent", 20}
DoubleProperty	m_zmin {this, "Zmin", -250.}
DoubleProperty	m_zmax {this, "Zmax", +250.}
DoubleProperty	m_ratio {this, "minRatio", 0.25}
std::list< Trk::Vertex >	production (const EventContext &ctx, SiSpacePointsSeedMakerEventData &data) const
MsgStream &	dumpConditions (MsgStream &out) const

Detailed Description

Class for primary vertices z-coordinates generation (xKalman algorithm).

In AthenaMT, event dependent cache inside SiZvertexMaker_xk is not preferred. SiSpacePointsSeedMakerEventData class holds event dependent data for SiZvertexMaker_xk and ISiSpacePointsSeedMaker. Its object is instantiated in SiSPSeededTrackFinder::execute.

Definition at line 49 of file SiZvertexMaker_xk.h.

Constructor & Destructor Documentation

SiZvertexMaker_xk()

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

Definition at line 27 of file SiZvertexMaker_xk.cxx.

  : base_class(t, n, p)
{
}

~SiZvertexMaker_xk()

virtual InDet::SiZvertexMaker_xk::~SiZvertexMaker_xk ( )

virtualdefault

Member Function Documentation

dump()

MsgStream & InDet::SiZvertexMaker_xk::dump ( MsgStream & out ) const

overridevirtual

Definition at line 220 of file SiZvertexMaker_xk.cxx.

{
  out << std::endl;
  return dumpConditions(out);
}

dumpConditions()

MsgStream & InDet::SiZvertexMaker_xk::dumpConditions ( MsgStream & out ) const

protected

Definition at line 230 of file SiZvertexMaker_xk.cxx.

{
  int n = 62-m_seedsgenerator.type().size();
  std::string s1;
  for (int i=0; i<n; ++i) s1.append(" ");
  s1.append("|");
 
  out<<"|----------------------------------------------------------------------"
     <<"-------------------|"
     <<std::endl;
  out<<"| SP seeds maker tool     | "
     <<m_seedsgenerator.type()<<s1<<std::endl;
  out<<"| Histogram size          | "
     <<std::setw(12)<<m_histsize
     <<"                                                  |"<<std::endl;
  out<<"| SP seed length (2 or 3) | "
     <<std::setw(12)<<m_nspoint
     <<"                                                  |"<<std::endl;
  out<<"| min content of hist.bin | "
     <<std::setw(12)<<m_mincontent
     <<"                                                  |"<<std::endl;
  out<<"| min vertex Z-coordinate | "
     <<std::setw(12)<<std::setprecision(4)<<m_zmin
     <<"                                                  |"<<std::endl;
  out<<"| max vertex Z-coordinate | "
     <<std::setw(12)<<std::setprecision(4)<<m_zmax
     <<"                                                  |"<<std::endl;
  out<<"| min. ratio to max bin   | "
     <<std::setw(12)<<std::setprecision(4)<<m_ratio
     <<"                                                  |"<<std::endl;
  out<<"|----------------------------------------------------------------------"
     <<"-------------------|"
     <<std::endl;
 
  return out;
}

finalize()

StatusCode InDet::SiZvertexMaker_xk::finalize ( )

overridevirtual

Definition at line 52 of file SiZvertexMaker_xk.cxx.

{
  return AlgTool::finalize();
}

initialize()

StatusCode InDet::SiZvertexMaker_xk::initialize ( )

overridevirtual

Definition at line 37 of file SiZvertexMaker_xk.cxx.

{
  // Get tool for seed generator
  //
  ATH_CHECK(m_seedsgenerator.retrieve());
 
  dump(msg(MSG::DEBUG));
 
  return StatusCode::SUCCESS;
}

newEvent()

std::list< Trk::Vertex > InDet::SiZvertexMaker_xk::newEvent ( const EventContext & ctx, SiSpacePointsSeedMakerEventData & data ) const

overridevirtual

Definition at line 61 of file SiZvertexMaker_xk.cxx.

{
  m_seedsgenerator->newEvent(ctx, data);
  return production(ctx, data);
}

newRegion() [1/2]

std::list< Trk::Vertex > InDet::SiZvertexMaker_xk::newRegion ( const EventContext & ctx, SiSpacePointsSeedMakerEventData & data, const std::vector< IdentifierHash > & vPixel, const std::vector< IdentifierHash > & vSCT ) const

overridevirtual

Definition at line 71 of file SiZvertexMaker_xk.cxx.

{
  m_seedsgenerator->newRegion(ctx, data, vPixel, vSCT);
  return production(ctx, data);
}

newRegion() [2/2]

std::list< Trk::Vertex > InDet::SiZvertexMaker_xk::newRegion ( const EventContext & ctx, SiSpacePointsSeedMakerEventData & data, const std::vector< IdentifierHash > & vPixel, const std::vector< IdentifierHash > & vSCT, const IRoiDescriptor & PhEt ) const

overridevirtual

Definition at line 83 of file SiZvertexMaker_xk.cxx.

{
  m_seedsgenerator->newRegion(ctx, data, vPixel, vSCT, PhEt);
  return production(ctx, data);
}

production()

std::list< Trk::Vertex > InDet::SiZvertexMaker_xk::production ( const EventContext & ctx, SiSpacePointsSeedMakerEventData & data ) const

protected

Definition at line 96 of file SiZvertexMaker_xk.cxx.

{
  std::list<Trk::Vertex> vertices;
 
  std::vector<int> histogram(m_histsize+1, 0);
  std::vector<double> z_histogram(m_histsize+1, 0.);
 
  double sZ = static_cast<double>(m_histsize)/(m_zmax-m_zmin);
 
  std::list<int>    NMAX;
  std::list<double> ZMAX;
  const Trk::SpacePoint* p0 = nullptr;
 
  std::list<Trk::Vertex> lv;  
  
  if      (m_nspoint==2) m_seedsgenerator->find2Sp(data, lv);
  else if (m_nspoint==3) m_seedsgenerator->find3Sp(ctx, data, lv);
  else                   m_seedsgenerator->findVSp(ctx, data, lv);
 
  const InDet::SiSpacePointsSeed* seed = nullptr;
  std::multimap<int,double> ver;
 
  int    nmax = 0 ;
  int    Hmax = 0 ;
  double zmax = 0.;
 
  while ((seed = m_seedsgenerator->next(ctx, data))) {
 
    std::vector<const Trk::SpacePoint*>::const_iterator
      s = seed->spacePoints().begin();
    ++s;
    if ((*s)!=p0) {
      if (Hmax!=0) {
        NMAX.push_back(nmax);
        ZMAX.push_back(zmax);
      }
      Hmax = nmax = 0;
      p0=(*s);
    }
    int n = static_cast<int>((seed->zVertex()-m_zmin)*sZ);
    if (n>=0 && n<=m_histsize) {
      if (++histogram[n]>Hmax) {
        Hmax = histogram[n];
        nmax = n;
        zmax = seed->zVertex();
      }
    } 
  } 
 
  // Clear histogram
  //
  for (int i=0; i<=m_histsize; ++i) histogram[i] = 0;
 
  // Fill histogram again using NMAX information
  //
  std::list<int>   ::iterator m, me = NMAX.end();
  std::list<double>::iterator z = ZMAX.begin();
  for (m=NMAX.begin(); m!=me; ++m) {
    ++histogram[(*m)];
    z_histogram[(*m)] += (*z++); 
  }
   
  // Max. content seach and output bank Z coordinates of vertices production
  //
  int HiZm = 0; 
  int im   = 0;
  int i    = 1;
  while (true) {
 
    if (histogram[i]>=histogram[i-1] && histogram[i]>=histogram[i+1]) {
      int H3 = histogram[i-1]+histogram[i]+histogram[i+1]; 
      if (H3>HiZm) {
        HiZm=H3;
        im=i;
      }
      if (++i==m_histsize) break;
    }
    if (++i==m_histsize) break;
  }
 
  if (im==0 || HiZm < m_mincontent) {
 
    return vertices;
  }
  double zv = (z_histogram[im]+z_histogram[im-1]+z_histogram[im+1])/static_cast<double>(HiZm);
  ver.insert(std::make_pair(-HiZm, zv));
 
  // Z-vertex coordinates 
  //
  HiZm = static_cast<int>(m_ratio*static_cast<double>(HiZm));
  if (HiZm<m_mincontent) HiZm = m_mincontent;
  i = 1;
  while (true) {
    if (histogram[i]>=histogram[i-1] && histogram[i]>=histogram[i+1] && i!=im) { 
      int H3 = histogram[i-1]+histogram[i]+histogram[i+1];
      if (H3>=HiZm) {
        zv = (z_histogram[i]+z_histogram[i-1]+z_histogram[i+1])/static_cast<double>(H3);
        ver.insert(std::make_pair(-H3, zv));
      }
      if (++i==m_histsize) break;
    }
    if (++i==m_histsize) break;
  }
 
  if (!ver.empty()) {
    std::multimap<int, double>::iterator v = ver.begin(), ve = ver.end();
    for (; v!=ve; ++v) {
      Amg::Vector3D Vp(0., 0., (*v).second);
      Trk::Vertex Ver(Vp);
      vertices.push_back(Ver);
    }
  }
 
  dump(msg(MSG::DEBUG));
 
  return vertices;
}

Member Data Documentation

m_histsize

IntegerProperty InDet::SiZvertexMaker_xk::m_histsize {this, "HistSize", 500}

protected

Definition at line 105 of file SiZvertexMaker_xk.h.

{this, "HistSize", 500};

m_mincontent

IntegerProperty InDet::SiZvertexMaker_xk::m_mincontent {this, "minContent", 20}

protected

Definition at line 106 of file SiZvertexMaker_xk.h.

{this, "minContent", 20};

m_nspoint

IntegerProperty InDet::SiZvertexMaker_xk::m_nspoint {this, "SeedSize", 2}

protected

Definition at line 104 of file SiZvertexMaker_xk.h.

{this, "SeedSize", 2};

m_ratio

DoubleProperty InDet::SiZvertexMaker_xk::m_ratio {this, "minRatio", 0.25}

protected

Definition at line 109 of file SiZvertexMaker_xk.h.

{this, "minRatio", 0.25};

m_seedsgenerator

ToolHandle<InDet::ISiSpacePointsSeedMaker> InDet::SiZvertexMaker_xk::m_seedsgenerator {this, "SeedMakerTool", "InDet::SiSpacePointsSeedMaker_ATLxk"}

protected

Definition at line 98 of file SiZvertexMaker_xk.h.

{this, "SeedMakerTool", "InDet::SiSpacePointsSeedMaker_ATLxk"};

m_zmax

DoubleProperty InDet::SiZvertexMaker_xk::m_zmax {this, "Zmax", +250.}

protected

Definition at line 108 of file SiZvertexMaker_xk.h.

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

m_zmin

DoubleProperty InDet::SiZvertexMaker_xk::m_zmin {this, "Zmin", -250.}

protected

Definition at line 107 of file SiZvertexMaker_xk.h.

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

---

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

• SiZvertexMaker_xk.h
• SiZvertexMaker_xk.cxx