KalFitCylinder Class Reference

Cylinder is an Element whose shape is a cylinder. More...

#include <KalFitCylinder.h>

Inheritance diagram for KalFitCylinder:

Public Member Functions
	KalFitCylinder (const KalFitMaterial *material, double radius, double thick, double length, double z0)
	Constructor.
virtual double	intersect (const KalFitTrack &track, HepPoint3D &x) const
	Find intersection with Helix.
virtual double	intersect (const KalFitTrack &track, HepPoint3D &x, const HepPoint3D &point) const
bool	isInside (const HepPoint3D &x) const
	Check if the position x is inside the current cylinder.
bool	isInside2 (const HepPoint3D &x) const
virtual double	radius (void) const
	Extract radius of the cylinder.
Public Member Functions inherited from KalFitElement
	KalFitElement (const KalFitMaterial *material)
	Constructor.
virtual	~KalFitElement ()
	Destructor.
virtual double	intersect (const KalFitTrack &track, HepPoint3D &x) const =0
	Find intersection with Helix.
virtual double	intersect (const KalFitTrack &track, HepPoint3D &x, const HepPoint3D &point) const =0
const KalFitMaterial &	material (void) const
	Return Material of which Element is made.
virtual double	radius (void) const =0
	Return distance in r-phi plane.
virtual void	updateTrack (KalFitTrack &track, int index)
	Update track according the material properties of the current element.
virtual void	updateTrack_rphi (KalFitTrack &track, int index)
virtual void	updateTrack_alreadyfound (KalFitTrack &track, int index)
virtual void	asso_rphi (Lpav &circ, KalFitTrack &track)
virtual void	asso_rphi (KalFitTrack &track)
bool	operator< (const KalFitElement &a) const
	Operator.

Protected Attributes
double	ro_
double	ri_
double	zf_
double	zb_
Protected Attributes inherited from KalFitElement
const KalFitMaterial *	material_
	Material of this element.

Additional Inherited Members
Static Public Member Functions inherited from KalFitElement
static int	muls (void)
static int	loss (void)
static void	muls (int i)
static void	loss (int i)
Static Protected Attributes inherited from KalFitElement
static int	muls_
	Flags.
static int	loss_

Detailed Description

Cylinder is an Element whose shape is a cylinder.

Definition at line 21 of file KalFitCylinder.h.

Constructor & Destructor Documentation

KalFitCylinder()

KalFitCylinder::KalFitCylinder ( const KalFitMaterial *  material, double  radius, double  thick, double  length, double  z0  )

inline

Constructor.

Definition at line 30 of file KalFitCylinder.h.

    : KalFitElement(material)
  {
    ri_ = radius;
    ro_ = radius + thick;
    zf_ = z0 + length * .5;
    zb_ = z0 - length * .5;
  }

Member Function Documentation

intersect() [1/2]

double KalFitCylinder::intersect ( const KalFitTrack &  track, HepPoint3D &  x  ) const

virtual

Find intersection with Helix.

Implements KalFitElement.

Definition at line 23 of file KalFitCylinder.cxx.

{
  double dPhi[4];
  dPhi[0] = track.intersect_cylinder(ro_);
  if(dPhi[0] == 0) return -1;
  dPhi[1] = track.intersect_cylinder(ri_);
  if(dPhi[1] == 0) return -1;
  dPhi[2] = track.intersect_xy_plane(zf_);
  dPhi[3] = track.intersect_xy_plane(zb_);
 
  int n[2];
  int j = 0;
  for(int i = 0; i < 4 && j < 2; i++){
    HepPoint3D xx = track.x(dPhi[i]);
    if(isInside(xx)) n[j++] = i;
  }
  if(j < 2) return -1;
 
  x = track.x((dPhi[n[0]] + dPhi[n[1]]) * .5);
 
  double tanl = track.tanl();
  //cout<<"KalFitCylinder: track radius"<<track.radius()<<" dphi0 "
  //    <<dPhi[n[0]]<<" dphi1 "<<dPhi[n[1]]<<" tanl "<<tanl<<endl;
  return fabs(track.radius() * (dPhi[n[0]] - dPhi[n[1]])
             * sqrt(1 + tanl * tanl));
}

intersect() [2/2]

double KalFitCylinder::intersect ( const KalFitTrack &  track, HepPoint3D &  x, const HepPoint3D &  point  ) const

virtual

Implements KalFitElement.

Definition at line 52 of file KalFitCylinder.cxx.

{
 
  const double ro = sqrt(point.x()*point.x()+point.y()*point.y());
 
  //std::cout<<" ro: "<<ro<<std::endl;
  
  double dPhi[4];
  dPhi[0] = track.intersect_cylinder(ro);
  if(dPhi[0] == 0) return -1;
  dPhi[1] = track.intersect_cylinder(ri_);
  if(dPhi[1] == 0) return -1;
  dPhi[2] = track.intersect_xy_plane(zf_);
  dPhi[3] = track.intersect_xy_plane(zb_);
 
  //for(int ii=0; ii<4; ii++)
    //std::cout<<"dPhi["<<ii<<"]"<<dPhi[ii]<<std::endl;
  
  int n[2];
  int j = 0;
  for(int i = 0; i < 4 && j < 2; i++){
    HepPoint3D xx = track.x(dPhi[i]);
    if(isInside(xx)) n[j++] = i;
  }
 
  if(j < 2) return -1;
 
  x = track.x((dPhi[n[0]] + dPhi[n[1]]) * .5);
  
  double tanl = track.tanl();
 
  return fabs(track.radius() * (dPhi[n[0]] - dPhi[n[1]])
      * sqrt(1 + tanl * tanl));
 }

isInside()

bool KalFitCylinder::isInside

(

const HepPoint3D &

x

)

const

Check if the position x is inside the current cylinder.

Definition at line 91 of file KalFitCylinder.cxx.

{
  double r = x.perp();
  double z = x.z();
  //std::cout<<"r: "<<r<<" z: "<<z<<" ri: "<<ri_<<" ro: "<<ro_<<" zb_: "<<zb_<<"zf: "<<zf_<<std::endl;
  
  return (r >= ri_ - FLT_EPSILON &&
          r <= ro_ + FLT_EPSILON &&
          z >= zb_ - FLT_EPSILON &&
          z <= zf_ + FLT_EPSILON);
}

Referenced by intersect().

isInside2()

bool KalFitCylinder::isInside2

(

const HepPoint3D &

x

)

const

Definition at line 104 of file KalFitCylinder.cxx.

{
  double r = x.perp();
  double z = x.z();   
  //std::cout<<"r: "<<r<<" z: "<<z<<" ri: "<<ri_<<" ro: "<<ro_<<" zb_: "<<zb_<<"zf: "<<zf_<<std::endl;
 
  return (r <= ro_ + FLT_EPSILON &&
          z >= zb_ - FLT_EPSILON &&
          z <= zf_ + FLT_EPSILON);
}

radius()

virtual double KalFitCylinder::radius ( void  ) const

inlinevirtual

Extract radius of the cylinder.

Implements KalFitElement.

Definition at line 53 of file KalFitCylinder.h.

{ return ro_; }

Referenced by KalFitCylinder().

Member Data Documentation

ri_

double KalFitCylinder::ri_

protected

Definition at line 24 of file KalFitCylinder.h.

Referenced by intersect(), isInside(), and KalFitCylinder().

ro_

double KalFitCylinder::ro_

protected

Definition at line 23 of file KalFitCylinder.h.

Referenced by intersect(), isInside(), isInside2(), KalFitCylinder(), and radius().

zb_

double KalFitCylinder::zb_

protected

Definition at line 26 of file KalFitCylinder.h.

Referenced by intersect(), isInside(), isInside2(), and KalFitCylinder().

zf_

double KalFitCylinder::zf_

protected

Definition at line 25 of file KalFitCylinder.h.

Referenced by intersect(), isInside(), isInside2(), and KalFitCylinder().

---

The documentation for this class was generated from the following files:

• KalFitCylinder.h
• KalFitCylinder.cxx