#include <TkBfield.h>

Public Member Functions
void	getBrfz (double const __restrict__ x, double __restrict__ Brfz) const
	B out in cylindrical.
void	getBxyz (double const __restrict__ x, double __restrict__ Bxyz) const
	B out in cartesian.
	TkBfield (std::string T="3_8T")
	~TkBfield ()
Private Member Functions
void	Bcyl (double r, double z, double *__restrict__ Bw) const
Private Attributes
double	ainv
double	ap2
double	coeff
double	hb0
double	hlova
double	prm [9]

Detailed Description

B-field in Tracker volume - based on the TOSCA computation version 1103l (tuned on MTCC measured field (fall 2006))

In: x[3]: coordinates (m) Out: B[3]: Bx,By,Bz (T) (getBxyz) Out: B[3]: Br,Bf,Bz (T) (getBrfz)

Valid for r<1.15m and |z|<2.80m

Author:: V.Karimaki 080228, 080407

Definition at line 24 of file TkBfield.h.

Constructor & Destructor Documentation

TkBfield::TkBfield ( std::string T = "3_8T" )

Definition at line 13 of file TkBfield.cc.

References ainv, ap2, coeff, Exception, hb0, hlova, i, p1, p2, p3, p4, prm, and mathSSE::sqrt().

                                {
  double p1[]={4.90541,17.8768,2.02355,0.0210538,0.000321885,2.37511,0.00326725,2.07656,1.71879}; // 2.0T-2G
  double p2[]={4.41982,15.7732,3.02621,0.0197814,0.000515759,2.43385,0.00584258,2.11333,1.76079}; // 3.0T-2G
  double p3[]={4.30161,15.2586,3.51926,0.0183494,0.000606773,2.45110,0.00709986,2.12161,1.77038}; // 3.5T-2G
  double p4[]={4.24326,15.0201,3.81492,0.0178712,0.000656527,2.45818,0.00778695,2.12500,1.77436}; // 3.8T-2G
  double p5[]={4.21136,14.8824,4.01683,0.0175932,0.000695541,2.45311,0.00813447,2.11688,1.76076}; // 4.0T-2G
  prm[0]=0;
  if (fld=="2_0T") for (int i=0; i<9; i++) prm[i]=p1[i];
  if (fld=="3_0T") for (int i=0; i<9; i++) prm[i]=p2[i];
  if (fld=="3_5T") for (int i=0; i<9; i++) prm[i]=p3[i];
  if (fld=="3_8T") for (int i=0; i<9; i++) prm[i]=p4[i];
  if (fld=="4_0T") for (int i=0; i<9; i++) prm[i]=p5[i];
  //  cout<<std::endl<<"Instantiation of TkBfield with key "<<fld<<endl;
  if (!prm[0]) {
    throw cms::Exception("BadParameters") << "Undefined key - " // abort!\n";
    <<"Defined keys are: \"2_0T\" \"3_0T\" \"3_5T\" \"3_8T\" and \"4_0T\"\n";
    // exit(1);
  }
  ap2=4*prm[0]*prm[0]/(prm[1]*prm[1]);  
  hb0=0.5*prm[2]*std::sqrt(1.0+ap2);
  hlova=1/std::sqrt(ap2);
  ainv=2*hlova/prm[1];
  coeff=1/(prm[8]*prm[8]);
}

magfieldparam::TkBfield::~TkBfield ( ) [inline]

Definition at line 30 of file TkBfield.h.

{}

Member Function Documentation

void TkBfield::Bcyl	(	double	r,
		double	z,
		double *__restrict__	Bw
	)		const `[private]`

Definition at line 85 of file TkBfield.cc.

References abs, ainv, coeff, funct::exp(), hb0, hlova, prm, and v.

Referenced by getBrfz(), and getBxyz().

                                                                       {
  //  if (r<1.15&&fabs(z)<2.8) // NOTE: check omitted, is done already by the wrapper! (NA)
  z-=prm[3];                    // max Bz point is shifted in z
  double az=std::abs(z);
  double zainv=z*ainv;
  double u=hlova-zainv;
  double v=hlova+zainv;
  double fu[4],gv[4];
  ffunkti(u,fu);
  ffunkti(v,gv);
  double rat=0.5*r*ainv;
  double rat2=rat*rat;
  Bw[0]=hb0*rat*(fu[1]-gv[1]-(fu[3]-gv[3])*rat2*0.5);
  Bw[1]=0;
  Bw[2]=hb0*(fu[0]+gv[0]-(fu[2]+gv[2])*rat2);
  double corBr= prm[4]*r*z*(az-prm[5])*(az-prm[5]);
  //    double corBz=-prm[6]*exp(coeff*(az-prm[7])*(az-prm[7]));
  //    double corBz=-prm[6]/(1+coeff*(az-prm[7])*(az-prm[7]));
  //  double corBz=-prm[6]*(exp(-(z-prm[7])*(z-prm[7])/(prm[8]*prm[8]))
  //                    + exp(-(z+prm[7])*(z+prm[7])/(prm[8]*prm[8]))); // double Gaussian
  double corBz=-prm[6]*(exp(-(z-prm[7])*(z-prm[7])*coeff) +
                        exp(-(z+prm[7])*(z+prm[7])*coeff)
                        ); // double Gaussian
  Bw[0]+=corBr;
  Bw[2]+=corBz;
  //   } else {
  //     cout <<"TkBfield: The point is outside the region r<1.15m && |z|<2.8m"<<endl;
}

void TkBfield::getBrfz	(	double const *__restrict__	x,
		double *__restrict__	Brfz
	)		const

B out in cylindrical.

Definition at line 116 of file TkBfield.cc.

References Bcyl(), and mathSSE::sqrt().

                                                                                        {
  Bcyl(sqrt(x[0]*x[0]+x[1]*x[1]),x[2], Brfz);
}

void TkBfield::getBxyz	(	double const *__restrict__	x,
		double *__restrict__	Bxyz
	)		const

B out in cartesian.

Definition at line 120 of file TkBfield.cc.

References Bcyl(), alignCSCRings::r, and mathSSE::sqrt().

Referenced by OAEParametrizedMagneticField::inTeslaUnchecked().

                                                                                        {
  double Bw[3];
  double r=sqrt(x[0]*x[0]+x[1]*x[1]);
  Bcyl(r, x[2], Bw);
  double rinv=(r>0) ? 1/r:0;
  Bxyz[0]=Bw[0]*x[0]*rinv;
  Bxyz[1]=Bw[0]*x[1]*rinv;
  Bxyz[2]=Bw[2];
}