L1TrackPar.h

Go to the documentation of this file.

#ifndef L1TrackPar_h
#define L1TrackPar_h 1
 
#include "CbmL1Def.h"
 
class L1TrackPar{
 
  public:
 
    fvec
      x,y,tx,ty,qp,z,
      C00,
      C10, C11,
      C20, C21, C22,
      C30, C31, C32, C33,
      C40, C41, C42, C43, C44,
      chi2, NDF                ;
 
    L1TrackPar():
        x(0),y(0),tx(0),ty(0),qp(0),z(0),
    C00(0),
    C10(0), C11(0),
    C20(0), C21(0), C22(0),
    C30(0), C31(0), C32(0), C33(0),
    C40(0), C41(0), C42(0), C43(0), C44(0),
    chi2(0), NDF(0)
    {};
    L1TrackPar(double *T, double *C):
      x(T[0]),
      y(T[1]),
      tx(T[2]),
      ty(T[3]),
      qp(T[4]),
      z(T[5]),
 
      C00(C[0]),
      C10(C[1]),
      C11(C[2]),
      C20(C[3]),
      C21(C[4]),
      C22(C[5]),
      C30(C[6]),
      C31(C[7]),
      C32(C[8]),
      C33(C[9]),
      C40(C[10]),
      C41(C[11]),
      C42(C[12]),
      C43(C[13]),
      C44(C[14]),
      chi2(0), NDF(0)
      {};
      
    void SetOneEntry( const int i0, const L1TrackPar &T1, const int i1 );
  
    void Print(int i = -1);
// void L1Extrapolate
// ( 
// //  L1TrackPar &T, // input track parameters (x,y,tx,ty,Q/p) and cov.matrix
//  fvec        z_out  , // extrapolate to this z position
//  fvec       qp0    , // use Q/p linearisation at this value
//  L1FieldRegion &F
//  );
 
    
} _fvecalignment;
 
 // =============================================================================================
 
inline void L1TrackPar::Print(int i)
{
  std::cout.setf(std::ios::scientific,std::ios::floatfield);
  if (i == -1){
    std::cout << "T = " << std::endl;
    std::cout << x  << std::endl;
    std::cout << y  << std::endl;
    std::cout << tx << std::endl;
    std::cout << ty << std::endl;
    std::cout << qp << std::endl;
    std::cout << z  << std::endl;
  }
  else{
    std::cout << "T = ";
    std::cout << x[i]  << " ";
    std::cout << y[i]  << " ";
    std::cout << tx[i] << " ";
    std::cout << ty[i] << " ";
    std::cout << qp[i] << " ";
    std::cout << z[i]  << std::endl;
    std::cout << "C = ";
    std::cout << C00[i]  << " ";
    std::cout << C11[i]  << " ";
    std::cout << C22[i]  << " ";
    std::cout << C33[i]  << " ";
    std::cout << C44[i]  << std::endl;
  }
}
 
inline void L1TrackPar::SetOneEntry( const int i0, const L1TrackPar &T1, const int i1 )
{
  x[i0] = T1.x[i1];
  y[i0] = T1.y[i1];
  tx[i0] = T1.tx[i1];
  ty[i0] = T1.ty[i1];
  qp[i0] = T1.qp[i1];
  z[i0] = T1.z[i1];
  C00[i0] = T1.C00[i1];
  C10[i0] = T1.C10[i1];
  C11[i0] = T1.C11[i1];
  C20[i0] = T1.C20[i1];
  C21[i0] = T1.C21[i1];
  C22[i0] = T1.C22[i1];
  C30[i0] = T1.C30[i1];
  C31[i0] = T1.C31[i1];
  C32[i0] = T1.C32[i1];
  C33[i0] = T1.C33[i1];
  C40[i0] = T1.C40[i1];
  C41[i0] = T1.C41[i1];
  C42[i0] = T1.C42[i1];
  C43[i0] = T1.C43[i1];
  C44[i0] = T1.C44[i1];
   
  chi2[i0] = T1.chi2[i1];
  NDF[i0] = T1.NDF[i1];
} // SetOneEntry
 
#endif