KalFitElement.cxx

Go to the documentation of this file.

//-----------------------------------------------------------------------
// File from KalFit module
//
// Filename : KalFitElement.cc
//------------------------------------------------------------------------
// Description :
// Element is a base class which represents an element of detector
//------------------------------------------------------------------------
// Modif :
//------------------------------------------------------------------------
#include "CLHEP/Geometry/Point3D.h"
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Point3D<double> HepPoint3D;
#endif
// #include "TrackUtil/Helix.h"
#include "KalFitAlg/helix/Helix.h"
 
#include "KalFitAlg/KalFitElement.h"
#include "KalFitAlg/KalFitMaterial.h"
#include "KalFitAlg/KalFitTrack.h"
 
int KalFitElement::loss_;
int KalFitElement::muls_;
 
void KalFitElement::updateTrack( KalFitTrack& track, int index ) {
 
  HepPoint3D x;
  double     path = intersect( track, x );
 
  // cout<<"KalFitElement: path= "<<path<<" intersect x "<<x<<endl;
 
  if ( path > 0 )
  {
    // move pivot
 
    // std::cout<<"KalFitElement: track helix1= "<<track.a()<<std::endl;
    // std::cout<<" KalFitTrack::numf_: "<<KalFitTrack::numf_<<std::endl;
    // std::cout<<" KalFitTrack::muls_: "<<muls_<<std::endl;
    // std::cout<<" KalFitTrack::loss_: "<<loss_<<std::endl;
 
    track.pivot_numf( x );
 
    // std::cout<<"KalFitElement: track helix2= "<<track.a()<<std::endl;
 
    // multiple scattering and energy loss
    int index_element( index );
    if ( index_element == 0 ) index_element = 1;
    if ( muls_ ) track.ms( path, *material_, index_element );
    if ( loss_ ) track.eloss( path, *material_, index_element );
 
    // cout<<"KalFitElement: track helix3= "<<track.a()<<endl;
  }
  // cout<<"KalfitElement: track helix2= "<<track.a()<<endl;
}
 
void KalFitElement::updateTrack_rphi( KalFitTrack& track, int index ) {
  HepPoint3D x;
  double     path = intersect( track, x );
  if ( path > 0 )
  {
    // move pivot
    track.pivot_numf( x );
 
    // multiple scattering and energy loss
    if ( muls_ ) track.ms( path, *material_, index );
    if ( loss_ ) track.eloss( path, *material_, index );
  }
}
 
void KalFitElement::updateTrack_alreadyfound( KalFitTrack& track, int index ) {
  HepPoint3D x;
  double     path = intersect( track, x );
  if ( path > 0 )
  {
    // move pivot
    track.pivot_numf( x );
 
    // multiple scattering and energy loss
    if ( muls_ ) track.ms( path, *material_, index );
    if ( loss_ ) track.eloss( path, *material_, index );
  }
}
 
void KalFitElement::asso_rphi( Lpav& circ, KalFitTrack& track ) {}
 
void KalFitElement::asso_rphi( KalFitTrack& track ) {}
 
void KalFitElement::muls( int i ) { muls_ = i; }
 
void KalFitElement::loss( int i ) { loss_ = i; }
 
int KalFitElement::muls( void ) { return muls_; }
 
int KalFitElement::loss( void ) { return loss_; }