KalFitWire.cxx

Go to the documentation of this file.

//-----------------------------------------------------------------------
// File from KalFit module
//
// Filename : KalFitWire.cc
//------------------------------------------------------------------------
// Description : Description of a Mdc Wire
//------------------------------------------------------------------------
// Modif :
//------------------------------------------------------------------------
#include "CLHEP/Geometry/Point3D.h"
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Point3D<double> HepPoint3D;
#endif
#include "KalFitAlg/KalFitLayer_Mdc.h"
#include "KalFitAlg/KalFitSuper_Mdc.h"
#include "KalFitAlg/KalFitWire.h"
 
const double KalFitWire::A[NREGION] = { 8.5265E-7, 1.1368E-6, 1.2402E-6 };
const double KalFitWire::F[NLAYER]  = {
    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
    0,    47.2, 46.6, 47.0, 46,   45,   46.7, 46,   45,   44.6, 45.9, 45.9, 45.3,
    44.5, 44.0, 46.0, 45.6, 44.2, 43.6, 44.0, 43.5, 43.1, 42.9, 42.4, 42.7, 42.5,
    42.2, 41,   42,   42,   42,   41.4, 42.5, 43.0, 43.2, 43.4, 43.4 };
const double L49_2 = 220.4 * 220.4;
 
// constructor
KalFitWire::KalFitWire( const int localID, const KalFitLayer_Mdc& layer, const HepPoint3D& fwd,
                        const HepPoint3D& bck, KalFitWire* const vt, unsigned int geoID,
                        unsigned int stereo )
    : localId_( localID )
    , fwd_( fwd )
    , bck_( bck )
    , layer_( layer )
    , state_( WireHitInvalid )
    , geoID_( geoID )
    , stereo_( stereo )
    , ddl_( 0 )
    , ddr_( 0 )
    , distance_( 0 ) {
 
  int layer_ID = layer_.layerId();
  xyPosition_  = 0.5 * ( fwd_ + bck_ );
  if ( !stereo_ )
  {
    dx_ = 0;
    dy_ = 0;
    x_  = xyPosition_.x();
    y_  = xyPosition_.y();
  }
  else
  {
    dx_ = (double)fwd_.x() - (double)bck_.x();
    dy_ = (double)fwd_.y() - (double)bck_.y();
    x_  = bck_.x();
    y_  = bck_.y();
  }
  // Wire sag coeff
  Hep3Vector wire_;
  wire_ = (CLHEP::Hep3Vector)fwd - (CLHEP::Hep3Vector)bck;
  lzx_  = sqrt( wire_.z() * wire_.z() + wire_.x() * wire_.x() );
 
  // problem with layer_ and its id inside sagcoef ?!!
  //  sagcoef();
  if ( layer_ID < 3 ) A_ = A[0];
  else
  {
    if ( layer_ID < 14 ) A_ = A[1];
    else if ( layer_ID < 50 )
    {
      double f_current( F[layer_ID] );
      double l( wire_.mag() );
      A_ = A[2] * F[49] * F[49] * L49_2 / ( f_current * f_current * l * l );
    }
    else std::cout << "*** PROBLEM WIRE !!!!! " << std::endl;
  }
 
  // Neighbor :
  neighbor_[0] = (KalFitWire*)innerLeft( vt );
  neighbor_[1] = (KalFitWire*)innerRight( vt );
 
  neighbor_[2] = (KalFitWire*)left();
  neighbor_[3] = (KalFitWire*)right();
 
  neighbor_[4] = (KalFitWire*)outerLeft( vt );
  neighbor_[5] = (KalFitWire*)outerRight( vt );
}
 
KalFitWire::KalFitWire()
    : localId_( 0 )
    , fwd_( 0, 0, 0 )
    , bck_( 0, 0, 0 )
    , layer_( *(KalFitLayer_Mdc*)NULL )
    , state_( WireHitInvalid )
    , x_( 0 )
    , y_( 0 )
    , dx_( 0 )
    , dy_( 0 ) {
  // Neighbor :
  neighbor_[0] = NULL;
  neighbor_[1] = NULL;
  neighbor_[2] = NULL;
  neighbor_[3] = NULL;
  neighbor_[4] = NULL;
  neighbor_[5] = NULL;
}
 
// destructor
KalFitWire::~KalFitWire( void ) {}
 
// wire sag coefficient determination :
void KalFitWire::sagcoef( void ) {
 
  Hep3Vector wire_;
  wire_        = (CLHEP::Hep3Vector)fwd_ - (CLHEP::Hep3Vector)bck_;
  int layer_ID = layer_.layerId();
 
  if ( layer_ID < 3 ) { A_ = A[0]; }
  else
  {
    if ( layer_ID < 14 ) { A_ = A[1]; }
    else if ( layer_ID < 50 )
    {
      double f_current( F[layer_ID] );
      double l( wire_.mag() );
      A_ = A[2] * F[49] * F[49] * L49_2 / ( f_current * f_current * l * l );
    }
    else { std::cout << "*** PROBLEM WIRE !!!!! " << std::endl; }
  }
}
 
const KalFitWire* KalFitWire::left( void ) const {
  const KalFitWire* tmp = this;
  if ( !localId_ ) tmp += layer_.superLayer().nWire();
  return --tmp;
}
 
const KalFitWire* KalFitWire::right( void ) const {
  const KalFitWire* tmp = this;
  if ( localId_ == layer_.superLayer().localMaxId() ) { tmp -= layer_.superLayer().nWire(); }
  return ++tmp;
}
 
const KalFitWire* KalFitWire::innerLeft( KalFitWire* const vtWire ) const {
 
  if ( !layer_.localLayerId() ) return vtWire;
  const KalFitWire* tmp = this;
  if ( layer_.offset() & 1 ) { tmp -= layer_.superLayer().nWire(); }
  else
  {
    if ( localId_ ) tmp -= layer_.superLayer().nWire();
    --tmp;
  }
  return tmp;
}
 
const KalFitWire* KalFitWire::innerRight( KalFitWire* const vtWire ) const {
  if ( !layer_.localLayerId() ) return vtWire;
  const KalFitWire* tmp = this;
  if ( layer_.offset() & 1 )
  {
    if ( localId_ == layer_.superLayer().localMaxId() ) { tmp -= layer_.superLayer().nWire(); }
    tmp -= ( layer_.superLayer().nWire() - 1 );
  }
  else { tmp -= layer_.superLayer().nWire(); }
  return tmp;
}
 
const KalFitWire* KalFitWire::outerLeft( KalFitWire* const vtWire ) const {
  if ( layer_.localLayerId() == layer_.superLayer().layerMaxId() ) return vtWire;
  const KalFitWire* tmp = this;
  if ( layer_.offset() & 1 ) { tmp += layer_.superLayer().nWire(); }
  else
  {
    if ( !localId_ ) tmp += layer_.superLayer().nWire();
    tmp += ( layer_.superLayer().nWire() - 1 );
  }
  return tmp;
}
 
const KalFitWire* KalFitWire::outerRight( KalFitWire* const vtWire ) const {
  if ( layer_.localLayerId() == layer_.superLayer().layerMaxId() ) return vtWire;
  const KalFitWire* tmp = this;
  if ( layer_.offset() & 1 )
  {
    if ( localId_ ^ layer_.superLayer().localMaxId() ) { tmp += layer_.superLayer().nWire(); }
    tmp++;
  }
  else { tmp += layer_.superLayer().nWire(); }
  return tmp;
}
 
// Check left and right side to see if there is hit valid
void KalFitWire::chk_left_and_right( void ) {
  if ( ( ( **( neighbor_ + 2 ) ).state_ & WireHit ) &&
       ( ( **( neighbor_ + 3 ) ).state_ & WireHit ) )
  {
    state_ |= WireHitInvalid;
    ( **( neighbor_ + 2 ) ).state_ |= WireHitInvalid;
    ( **( neighbor_ + 3 ) ).state_ |= WireHitInvalid;
  }
}