xAOD::TrackingDetails Namespace Reference

Typedefs
using	GenVecFourMom_t = ROOT::Math::LorentzVector<ROOT::Math::PxPyPzM4D<double> >
	Base 4 Momentum type for TrackParticle.

Enumerations
enum	covMatrixIndex {   d0_index =0 , z0_index =1 , phi_index =2 , th_index =3 ,   qp_index =4 }

Functions
float	charge (float qOverP)
GenVecFourMom_t	genvecP4 (float qOverP, float thetaT, float phiT, double m)
bool	definingParametersCovMatrixOffDiagCompr (std::span< const float > covMatrixOffDiag)
xAOD::ParametersCovMatrix_t	definingParametersCovMatrix (std::span< const float > covMatrixDiag, std::span< const float > covMatrixOffDiag, bool &valid)

Variables
constexpr std::size_t	COVMATRIX_OFFDIAG_VEC_COMPR_SIZE = 6
constexpr std::array< std::pair< covMatrixIndex, covMatrixIndex >, COVMATRIX_OFFDIAG_VEC_COMPR_SIZE >	covMatrixComprIndexPairs

Typedef Documentation

GenVecFourMom_t

using xAOD::TrackingDetails::GenVecFourMom_t = ROOT::Math::LorentzVector<ROOT::Math::PxPyPzM4D<double> >

Base 4 Momentum type for TrackParticle.

Definition at line 38 of file TrackingDetails.h.

Enumeration Type Documentation

covMatrixIndex

enum xAOD::TrackingDetails::covMatrixIndex

Enumerator
d0_index
z0_index
phi_index
th_index
qp_index

Definition at line 40 of file TrackingDetails.h.

{d0_index=0, z0_index=1, phi_index=2, th_index=3, qp_index=4};

Function Documentation

charge()

float xAOD::TrackingDetails::charge ( float qOverP )

inlinenodiscard

Definition at line 47 of file TrackingDetails.h.

                                                     {
      return (qOverP > 0) ? 1 : ((qOverP < 0) ? -1 : 0);
    }

definingParametersCovMatrix()

xAOD::ParametersCovMatrix_t xAOD::TrackingDetails::definingParametersCovMatrix ( std::span< const float > covMatrixDiag, std::span< const float > covMatrixOffDiag, bool & valid )

inlinenodiscard

Definition at line 74 of file TrackingDetails.h.

                                                                                                                                                                             {
 
      // Set up the result matrix.
      xAOD::ParametersCovMatrix_t cov;
      cov.setZero();
      valid = true;
 
      // Set the diagonal elements of the matrix.
      if( ( static_cast< int >( covMatrixDiag.size() ) == cov.rows() ) ) {
 
        // Access the "raw" variable.
        // Set the diagonal elements using the raw variable.
        for( int i = 0; i < cov.rows(); ++i ) {
          cov( i, i ) = covMatrixDiag[ i ];
        }
      } else {
        valid = false;
        // If the variable is not available/set, set the matrix to identity.
        cov.setIdentity();
      }
 
      bool offDiagCompr = definingParametersCovMatrixOffDiagCompr(covMatrixOffDiag);
 
      // Set the off-diagonal elements of the matrix.
      if(!offDiagCompr){
 
        if( ( static_cast< int >( covMatrixOffDiag.size() ) == ( ( ( cov.rows() - 1 ) * cov.rows() ) / 2 ) ) ) {
 
          // Set the off-diagonal elements using the raw variable.
          std::size_t vecIndex = 0;
          for( int i = 1; i < cov.rows(); ++i ) {
            for( int j = 0; j < i; ++j, ++vecIndex ) {
              float offDiagCoeff = cov(i,i)>0 && cov(j,j)>0 ? covMatrixOffDiag[vecIndex]*sqrt(cov(i,i)*cov(j,j)) : 0;
              cov.fillSymmetric( i, j, offDiagCoeff );
            }
          }
        }
 
        else valid = false;
 
      }
 
      else{ //Compressed case
 
        if( ( static_cast< int >( covMatrixOffDiag.size() ) == COVMATRIX_OFFDIAG_VEC_COMPR_SIZE ) ) {
          // Set the off-diagonal elements using the raw variable.
 
          const auto& vecPairIndex = covMatrixComprIndexPairs;
 
          for(unsigned int k=0; k<COVMATRIX_OFFDIAG_VEC_COMPR_SIZE; ++k){
            std::pair<covMatrixIndex,covMatrixIndex> pairIndex = vecPairIndex[k];
            covMatrixIndex i = pairIndex.first;
            covMatrixIndex j = pairIndex.second;
            float offDiagCoeff = cov(i,i)>0 && cov(j,j)>0 ? covMatrixOffDiag[k]*sqrt(cov(i,i)*cov(j,j)) : 0;
            cov.fillSymmetric( i, j, offDiagCoeff );
          }
 
        }
 
        else valid = false;
 
      }
 
 
      // Return the filled matrix.
      return cov;
 
    }

definingParametersCovMatrixOffDiagCompr()

bool xAOD::TrackingDetails::definingParametersCovMatrixOffDiagCompr ( std::span< const float > covMatrixOffDiag )

inlinenodiscard

Definition at line 65 of file TrackingDetails.h.

                                                                                                               {
 
      bool flag = false;
      flag = (static_cast< int >(covMatrixOffDiag.size())==COVMATRIX_OFFDIAG_VEC_COMPR_SIZE);
      return flag;
    }

genvecP4()

GenVecFourMom_t xAOD::TrackingDetails::genvecP4 ( float qOverP, float thetaT, float phiT, double m )

inlinenodiscard

Definition at line 52 of file TrackingDetails.h.

                                                                                                    {
      using namespace std;
      float p = 10.e6; // 10 TeV (default value for very high pt muons, with qOverP==0)
      if (fabs(qOverP)>0.) p = 1/fabs(qOverP);
      float sinTheta= sin(thetaT);
      float px = p*sinTheta*cos(phiT);
      float py = p*sinTheta*sin(phiT);
      float pz = p*cos(thetaT);
      return GenVecFourMom_t(px, py, pz, m);
    }

Variable Documentation

COVMATRIX_OFFDIAG_VEC_COMPR_SIZE

std::size_t xAOD::TrackingDetails::COVMATRIX_OFFDIAG_VEC_COMPR_SIZE = 6

constexpr

Definition at line 41 of file TrackingDetails.h.

covMatrixComprIndexPairs

std::array< std::pair<covMatrixIndex,covMatrixIndex>, COVMATRIX_OFFDIAG_VEC_COMPR_SIZE > xAOD::TrackingDetails::covMatrixComprIndexPairs

constexpr

Initial value:

{{
      {d0_index,phi_index}, {z0_index,th_index}, {d0_index,qp_index},
      {z0_index,qp_index}, {phi_index,qp_index}, {th_index,qp_index} }}

Definition at line 42 of file TrackingDetails.h.

                                                                                                                            {{
      {d0_index,phi_index}, {z0_index,th_index}, {d0_index,qp_index},
      {z0_index,qp_index}, {phi_index,qp_index}, {th_index,qp_index} }};