SiSpacePointForSeed.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#include <cmath>
 
 
 
#include "SiSPSeededTrackFinderData/SiSpacePointForSeed.h"
 
#include "InDetPrepRawData/SiCluster.h"
#include "InDetReadoutGeometry/SiDetectorElement.h"
#include "TrkSpacePoint/SpacePoint.h"
#include "TrkSurfaces/Surface.h"
 
namespace InDet {
 
  SiSpacePointForSeed::SiSpacePointForSeed
  (const Trk::SpacePoint* sp, std::span<float const, 3> r) 
  {
    set(sp,r); m_param = 0.;
  }
 
  SiSpacePointForSeed::SiSpacePointForSeed
  (const Trk::SpacePoint* sp, std::span<float const, 3> r, std::span<float const, 4> sc) 
  {
    set(sp,r,sc); m_param = 0.;
  }
 

  // Set
 
  void SiSpacePointForSeed::set
  (const Trk::SpacePoint* sp,std::span<float const, 3> r)
  {
    spacepoint = sp  ;
    m_x        = r[0];
    m_y        = r[1];
    m_z        = r[2];
    m_r        =std::sqrt(m_x*m_x+m_y*m_y);
    m_q        = 100000.;
 
    const InDet::SiCluster*           c  = static_cast<const InDet::SiCluster*>(sp->clusterList().first);
    const InDetDD::SiDetectorElement* de = c ->detectorElement();
 
    if( de->isPixel() ) {
      
      const Amg::MatrixX& v =  c->localCovariance();
      float f22 = float(v(1,1) );
      float wid = float(c->width().z());
      float cov = wid*wid*.08333f; if(cov < f22) cov = f22;
      if(de->isBarrel()) {m_covz = 9.f*cov; m_covr = .06f;}
      else               {m_covr = 9.f*cov; m_covz = .06f;}
      m_sn = nullptr;
    }
    else                {
 
      const Amg::MatrixX& v = sp->localCovariance();
      float f22 = float(v(1,1));
      if(de->isBarrel()) {m_covz = 8.f*f22; m_covr = .1f;} 
      else               {m_covr = 8.f*f22; m_covz = .1f;} 
      m_sn =  &sp->clusterList().second->detectorElement()->surface();
    }
    m_su = &sp->clusterList().first->detectorElement()->surface();
  } 

  // Set with error correction 
  // sc[0] - barrel pixels error correction
  // sc[1] - endcap pixels 
  // sc[2] - barrel sct
  // sc[3] - endcap sct 
 
  void SiSpacePointForSeed::set
  (const Trk::SpacePoint* sp,std::span<float const, 3> r,std::span<float const, 4> sc)
  {
    spacepoint = sp  ;
    m_x        = r[0];
    m_y        = r[1];
    m_z        = r[2];
    m_r        =std::sqrt(m_x*m_x+m_y*m_y);
    m_q        = 100000.;
 
    const InDet::SiCluster*           c  = static_cast<const InDet::SiCluster*>(sp->clusterList().first);
    const InDetDD::SiDetectorElement* de = c ->detectorElement();
 
    if( de->isPixel() ) {
      
      const Amg::MatrixX& v =  c->localCovariance();
      float f22 = float(v(1,1));
      float wid = float(c->width().z());
      float cov = wid*wid*.08333f; if(cov < f22) cov = f22;
      if(de->isBarrel()) {m_covz = 9.f*cov*sc[0]; m_covr = .06f;}
      else               {m_covr = 9.f*cov*sc[1]; m_covz = .06f;}
      m_sn = nullptr;
    }
    else                {
 
      const Amg::MatrixX& v = sp->localCovariance();
      float f22 = float(v(1,1));
      if(de->isBarrel()) {m_covz = 8.f*f22*sc[2]; m_covr = .1f;} 
      else               {m_covr = 8.f*f22*sc[3]; m_covz = .1f;} 
      m_sn =  &sp->clusterList().second->detectorElement()->surface();
    }
    m_su = &sp->clusterList().first->detectorElement()->surface();
  }
} // end of name space