#include <L1Field.h>

Public Member Functions
	L1FieldRegion ()

	L1FieldRegion (float reg[10])

CAFieldValue	Get (const float_v z)

void	Set (const CAFieldValue &B0, const float_v B0z, const CAFieldValue &B1, const float_v B1z, const CAFieldValue &B2, const float_v B2z)

void	Set (const CAFieldValue &B0, const float_v B0z, const CAFieldValue &B1, const float_v B1z)

void	Shift (float_v z)

void	SetOneEntry (const int i0, const L1FieldRegion &f1, const int i1)

void	SetOneEntry (const L1FieldRegion &f1, const int i1)

void	GetOneEntry (float reg[10], const int iVec)

Public Attributes
float_v	cx0

float_v	cx1

float_v	cx2

float_v	cy0

float_v	cy1

float_v	cy2

float_v	cz0

float_v	cz1

float_v	cz2

float_v	z0

Friends
ostream &	operator<< (ostream &out, L1FieldRegion &B)

Detailed Description

Definition at line 115 of file L1Field.h.

Constructor & Destructor Documentation

L1FieldRegion::L1FieldRegion ( )

inline

Definition at line 118 of file L1Field.h.

                    :
       cx0(Vc::Zero), cx1(Vc::Zero), cx2(Vc::Zero),
       cy0(Vc::Zero), cy1(Vc::Zero), cy2(Vc::Zero),
       cz0(Vc::Zero), cz1(Vc::Zero), cz2(Vc::Zero),
       z0(Vc::Zero)
     {}

L1FieldRegion::L1FieldRegion ( float reg[10] )

inline

Definition at line 125 of file L1Field.h.

                                 :
       cx0(reg[0]), cx1(reg[1]), cx2(reg[2]),
       cy0(reg[3]), cy1(reg[4]), cy2(reg[5]),
       cz0(reg[6]), cz1(reg[7]), cz2(reg[8]),
       z0(reg[9])
     {}

Member Function Documentation

CAFieldValue L1FieldRegion::Get ( const float_v z )

inline

Definition at line 137 of file L1Field.h.

References cx0, cx1, cx2, cy0, cy1, cy2, cz0, cz1, cz2, dz, and z0.

                                      {
       float_v dz = (z-z0);
       float_v dz2 = dz*dz;
       CAFieldValue B;
       B.x = cx0 + cx1*dz + cx2*dz2;
       B.y = cy0 + cy1*dz + cy2*dz2;
       B.z = cz0 + cz1*dz + cz2*dz2;
       return B;
     }

void L1FieldRegion::GetOneEntry	(	float	reg[10],
		const int	iVec
	)

inline

Definition at line 239 of file L1Field.h.

References cx0, cx1, cx2, cy0, cy1, cy2, cz0, cz1, cz2, and z0.

     {
       reg[0] = cx0[iVec];
       reg[1] = cx1[iVec];
       reg[2] = cx2[iVec];
       reg[3] = cy0[iVec];
       reg[4] = cy1[iVec];
       reg[5] = cy2[iVec];
       reg[6] = cz0[iVec];
       reg[7] = cz1[iVec];
       reg[8] = cz2[iVec];
       reg[9] = z0[iVec];
     }

void L1FieldRegion::Set	(	const CAFieldValue &	B0,
		const float_v	B0z,
		const CAFieldValue &	B1,
		const float_v	B1z,
		const CAFieldValue &	B2,
		const float_v	B2z
	)

inline

Definition at line 147 of file L1Field.h.

References cx0, cx1, cx2, cy0, cy1, cy2, cz0, cz1, cz2, f, rcp(), and z0.

     {
       z0 = B0z;
       float_v dz1 = B1z-B0z, dz2 = B2z-B0z;
       
       float_v deti = dz1*dz2*(dz2-dz1);
       float_m mask = abs(deti) > float_v(1.e-8f);
       
       float_v det = rcp(float_v(dz1*dz2*(dz2-dz1)));
       float_v w21 = -dz2*det;
       float_v w22 = dz1*det;
       float_v w11 = -dz2*w21;
       float_v w12 = -dz1*w22;
             
       float_v dB1 = B1.x - B0.x;
       float_v dB2 = B2.x - B0.x;
       cx0(mask) = B0.x;
       cx1(mask) = dB1*w11 + dB2*w12 ;
       cx2(mask) = dB1*w21 + dB2*w22 ;
       
       dB1 = B1.y - B0.y;
       dB2 = B2.y - B0.y;
       cy0(mask) = B0.y;
       cy1(mask) = dB1*w11 + dB2*w12 ;
       cy2(mask) = dB1*w21 + dB2*w22  ;
 
       dB1 = B1.z - B0.z;
       dB2 = B2.z - B0.z;
       cz0(mask) = B0.z;
       cz1(mask) = dB1*w11 + dB2*w12 ;
       cz2(mask) = dB1*w21 + dB2*w22;
     }

void L1FieldRegion::Set	(	const CAFieldValue &	B0,
		const float_v	B0z,
		const CAFieldValue &	B1,
		const float_v	B1z
	)

inline

Definition at line 182 of file L1Field.h.

References cx0, cx1, cx2, cy0, cy1, cy2, cz0, cz1, cz2, rcp(), and z0.

     {
       z0 = B0z[0];
       float_v dzi = rcp(float_v( B1z - B0z));
       cx0 = B0.x;
       cy0 = B0.y;
       cz0 = B0.z;
       cx1 = ( B1.x - B0.x )*dzi;
       cy1 = ( B1.y - B0.y )*dzi;
       cz1 = ( B1.z - B0.z )*dzi;
       cx2 = cy2 = cz2 = 0;
     }

void L1FieldRegion::SetOneEntry	(	const int	i0,
		const L1FieldRegion &	f1,
		const int	i1
	)

inline

Definition at line 211 of file L1Field.h.

References cx0, cx1, cx2, cy0, cy1, cy2, cz0, cz1, cz2, and z0.

Referenced by KFParticleFinder::CombineTrackPart(), KFParticleFinder::Find2DaughterDecay(), and KFParticleFinder::FindTrackV0Decay().

     {
       cx0[i0] = f1.cx0[i1];
       cx1[i0] = f1.cx1[i1];
       cx2[i0] = f1.cx2[i1];
       cy0[i0] = f1.cy0[i1];
       cy1[i0] = f1.cy1[i1];
       cy2[i0] = f1.cy2[i1];
       cz0[i0] = f1.cz0[i1];
       cz1[i0] = f1.cz1[i1];
       cz2[i0] = f1.cz2[i1];
       z0[i0] = f1.z0[i1];
     }

void L1FieldRegion::SetOneEntry	(	const L1FieldRegion &	f1,
		const int	i1
	)

inline

Definition at line 225 of file L1Field.h.

References cx0, cx1, cx2, cy0, cy1, cy2, cz0, cz1, cz2, and z0.

     {
       cx0 = f1.cx0[i1];
       cx1 = f1.cx1[i1];
       cx2 = f1.cx2[i1];
       cy0 = f1.cy0[i1];
       cy1 = f1.cy1[i1];
       cy2 = f1.cy2[i1];
       cz0 = f1.cz0[i1];
       cz1 = f1.cz1[i1];
       cz2 = f1.cz2[i1];
       z0 = f1.z0[i1];
     }

void L1FieldRegion::Shift ( float_v z )

inline

Definition at line 196 of file L1Field.h.

References cx0, cx1, cx2, cy0, cy1, cy2, cz0, cz1, cz2, dz, and z0.

     {
       float_v dz = z-z0;
       float_v cx2dz = cx2*dz;
       float_v cy2dz = cy2*dz;
       float_v cz2dz = cz2*dz;
       z0 = float_v(z[0]);
       cx0+= ( cx1 + cx2dz )*dz;
       cy0+= ( cy1 + cy2dz )*dz;
       cz0+= ( cz1 + cz2dz )*dz;
       cx1+= cx2dz + cx2dz;
       cy1+= cy2dz + cy2dz;
       cz1+= cz2dz + cz2dz;
     }

Friends And Related Function Documentation

ostream& operator<<	(	ostream &	out,
		L1FieldRegion &	B
	)

friend

Definition at line 253 of file L1Field.h.

                                                             {
     return out << " FieldRegion " << endl <<
       B.cx0 << endl << B.cx1 << endl << B.cx2 << endl <<
       B.cy0 << endl << B.cy1 << endl << B.cy2 << endl <<
       B.cz0 << endl << B.cz1 << endl << B.cz2 << endl <<
       B.z0 << endl;
   };