RkFitCylinder Class Reference

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

#include <RkFitCylinder.h>

Public Member Functions
	RkFitCylinder (const RkFitMaterial *material, double radius, double thick, double length, double z0)
	Constructor.
double	intersect (TRunge &track, HepPoint3D &x) const
	Find intersection with Helix.
double	intersect (TRunge &track, HepPoint3D &x, const HepPoint3D &point) const
void	updateTrack (TRunge &track, double y[6]) 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 Attributes
const RkFitMaterial *	material_

Protected Attributes
double	ro_
double	ri_
double	zf_
double	zb_

Detailed Description

Cylinder is an Element whose shape is a cylinder.

Definition at line 21 of file RkFitCylinder.h.

Constructor & Destructor Documentation

RkFitCylinder()

RkFitCylinder::RkFitCylinder ( const RkFitMaterial *  material, double  radius, double  thick, double  length, double  z0  )

inline

Constructor.

Definition at line 30 of file RkFitCylinder.h.

  {
    material_=material;  
    ri_ = radius;
    ro_ = radius + thick;
    zf_ = z0 + length * .5;
    zb_ = z0 - length * .5;
  }

Member Function Documentation

intersect() [1/2]

double RkFitCylinder::intersect	(	TRunge &	track,
		HepPoint3D &	x
	)		const

Find intersection with Helix.

Definition at line 21 of file RkFitCylinder.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.helix().x(dPhi[i]);
    if(isInside(xx)) n[j++] = i;
  }
  if(j < 2) return -1;
 
  x = track.helix().x((dPhi[n[0]] + dPhi[n[1]]) * .5);
 
  double tanl = track.helix().tanl();
  //cout<<"RkFitCylinder: track radius"<<track.radius()<<" dphi0 "
  //    <<dPhi[n[0]]<<" dphi1 "<<dPhi[n[1]]<<" tanl "<<tanl<<endl;
  return fabs(track.helix().radius() * (dPhi[n[0]] - dPhi[n[1]])
             * sqrt(1 + tanl * tanl));
 //   return 0;
}

Referenced by updateTrack().

intersect() [2/2]

double RkFitCylinder::intersect	(	TRunge &	track,
		HepPoint3D &	x,
		const HepPoint3D &	point
	)		const

Definition at line 51 of file RkFitCylinder.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.helix().x(dPhi[i]);
   if(isInside(xx)) n[j++] = i;
 }
 
 if(j < 2) return -1;
 
 x = track.helix().x((dPhi[n[0]] + dPhi[n[1]]) * .5);
 
 double tanl = track.helix().tanl();
 
 return fabs(track.helix().radius() * (dPhi[n[0]] - dPhi[n[1]])
     * sqrt(1 + tanl * tanl));
 }

isInside()

bool RkFitCylinder::isInside

(

const HepPoint3D &

x

)

const

Check if the position x is inside the current cylinder.

Definition at line 102 of file RkFitCylinder.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 RkFitCylinder::isInside2

(

const HepPoint3D &

x

)

const

Definition at line 115 of file RkFitCylinder.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 RkFitCylinder::radius ( void  ) const

inlinevirtual

Extract radius of the cylinder.

Definition at line 54 of file RkFitCylinder.h.

{ return ro_; }

Referenced by RkFitCylinder().

updateTrack()

void RkFitCylinder::updateTrack	(	TRunge &	track,
		double	y[6]
	)		const

Definition at line 87 of file RkFitCylinder.cxx.

                                                              {
    HepPoint3D x;
    double path = intersect(track, x);
    double mass=0.000511;
    if(path > 0){
        // move pivot
        // multiple scattering and energy loss
       // if(muls_) track.ms(path, *material_, index_element);
        track.eloss(path, material_,mass,y,1);
    //    track.Propagate(path,y);
    }
      
}

Member Data Documentation

material_

const RkFitMaterial* RkFitCylinder::material_

Definition at line 28 of file RkFitCylinder.h.

Referenced by RkFitCylinder(), and updateTrack().

ri_

double RkFitCylinder::ri_

protected

Definition at line 24 of file RkFitCylinder.h.

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

ro_

double RkFitCylinder::ro_

protected

Definition at line 23 of file RkFitCylinder.h.

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

zb_

double RkFitCylinder::zb_

protected

Definition at line 26 of file RkFitCylinder.h.

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

zf_

double RkFitCylinder::zf_

protected

Definition at line 25 of file RkFitCylinder.h.

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

---

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

• RkFitCylinder.h
• RkFitCylinder.cxx