ShipBellField Class Reference

#include <ShipBellField.h>

Inheritance diagram for ShipBellField:

Collaboration diagram for ShipBellField:

Public Member Functions
	ShipBellField ()
	ShipBellField (const char *name, Double_t Bpeak, Double_t Zmiddle, Int_t fOrient=1, Double_t Btube=500.)
	ShipBellField (ShipFieldPar *fieldPar)
virtual	~ShipBellField ()
virtual Double_t	GetBx (Double_t x, Double_t y, Double_t z)
virtual Double_t	GetBy (Double_t x, Double_t y, Double_t z)
virtual Double_t	GetBz (Double_t x, Double_t y, Double_t z)
void	IncludeTarget (Double_t xy, Double_t z, Double_t l)
virtual void	Print ()

Private Member Functions
	ClassDef (ShipBellField, 2)

Private Attributes
Double_t	fPeak
Double_t	fMiddle
Double_t	fBtube
Int_t	fOrient
Bool_t	fInclTarget
Double_t	targetXY
Double_t	targetZ0
Double_t	targetL

Detailed Description

Definition at line 27 of file ShipBellField.h.

Constructor & Destructor Documentation

ShipBellField() [1/3]

ShipBellField::ShipBellField

(

)

Default constructor

Definition at line 25 of file ShipBellField.cxx.

  : FairField(),
    fPeak(0.),
    fMiddle(0.),
    fBtube(0.)
{
  fType = 1;
  fInclTarget=kFALSE;
}

ShipBellField() [2/3]

ShipBellField::ShipBellField	(	const char *	name,
		Double_t	Bpeak,
		Double_t	Zmiddle,
		Int_t	fOrient = 1,
		Double_t	Btube = 500.
	)

Standard constructor

Parameters

name	Object name
Bpeak	peak field..
Zmiddle	middle of the magnet (global coordinates)
Btube	largest radius of the tube ellips (inside)

Definition at line 39 of file ShipBellField.cxx.

  : FairField(name),
    fPeak(Peak),
    fMiddle(Middle),
    fBtube(Btube)
{
  fType=1;
  fInclTarget=kFALSE;
  fOrient = orientation;
  fBtube = Btube;
}

ShipBellField() [3/3]

ShipBellField::ShipBellField

(

ShipFieldPar *

fieldPar

)

Constructor from ShipFieldPar

Definition at line 55 of file ShipBellField.cxx.

  : FairField(),
    fPeak(0.),
    fMiddle(0.),
    fBtube(0.),
    fInclTarget(kFALSE)
{
  if ( ! fieldPar ) {
    cerr << "-W- ShipBellField::ShipBellField: empty parameter container!"
     << endl;
    fType=0;
  }
  else {
    fPeak   = fieldPar->GetPeak();
    fMiddle = fieldPar->GetMiddle();
    fBtube  = fieldPar->GetBtube();
    fType   = fieldPar->GetType();
  }
}

~ShipBellField()

ShipBellField::~ShipBellField ( )

virtual

Destructor

Definition at line 79 of file ShipBellField.cxx.

{ }

Member Function Documentation

ClassDef()

ShipBellField::ClassDef ( ShipBellField  , 2    )

private

GetBx()

Double_t ShipBellField::GetBx ( Double_t  x, Double_t  y, Double_t  z  )

virtual

Get components of field at a given point

Parameters

x,y,z

Point coordinates [cm]

Definition at line 90 of file ShipBellField.cxx.

                                                                {
  if (fOrient==1){ return 0.;}
  else {
    Double_t zlocal=fabs(z-fMiddle)/100.; //zlocal: convert cm->m
   //old Bell field from Wilfried Flegel
   //Double_t bx= fPeak/(1.+pow(fabs(zlocal)/2.1,6.));
   //new field based on simulation of Davide Tommasini (22/1/2015)
    
    Double_t bx=0.;
    //field in box 20 cm larger than inner tube.
    if ( (fabs(x)<2.7*m) && (fabs(y)<fBtube+0.2*m))  {
      if (zlocal<3.8) {
        bx=0.14361*exp( -0.5 * pow((zlocal-0.45479E-01)/2.5046,2.));
      }else if (zlocal<11.9) {
        bx=0.19532-0.61512E-01*zlocal+0.68447E-02*pow(zlocal,2.)-0.25672E-03*pow(zlocal,3.);
      }
      bx=((fPeak/tesla)/0.14361)*bx*tesla;
    }
    // cout << "Bell GetBX " << x << ", " << y << ", " << z << ", Bx= " << bx << endl;
   return bx;
  }
}

GetBy()

Double_t ShipBellField::GetBy ( Double_t  x, Double_t  y, Double_t  z  )

virtual

Definition at line 117 of file ShipBellField.cxx.

                                                                {
  Double_t by = 0.;
  if (fInclTarget && z < targetZ0+targetL && z > targetZ0){
 // check if in target area
    if (fabs(x)<targetXY && fabs(y)<targetXY) {by = 1*tesla;}
    if (fabs(x)>2*targetXY && fabs(x)<3*targetXY && fabs(y)<targetXY) {by = -1*tesla;}
  }
  else if (fOrient==1){
   Double_t zlocal=(z-fMiddle)/100.;
   by= fPeak/(1.+pow(fabs(zlocal)/2.1,6.));
   //cout << "Bell GetBY " << z << ", By= " << by << endl;
  }
  return by;
 
}

GetBz()

Double_t ShipBellField::GetBz ( Double_t  x, Double_t  y, Double_t  z  )

virtual

Definition at line 137 of file ShipBellField.cxx.

                                                                {
  return 0.;
}

IncludeTarget()

void ShipBellField::IncludeTarget	(	Double_t	xy,
		Double_t	z,
		Double_t	l
	)

Definition at line 83 of file ShipBellField.cxx.

                                                                     {
  fInclTarget = kTRUE;
  targetXY = xy;
  targetZ0 = z;
  targetL = l;
}

Print()

void ShipBellField::Print ( )

virtual

Screen output

Definition at line 145 of file ShipBellField.cxx.

                          {
  cout << "======================================================" << endl;
  cout << "----  " << fTitle << " : " << fName << endl;
  cout << "----" << endl;
  cout << "----  Field type    : constant" << endl;
  cout << "----" << endl;
  cout << "----  Field regions : " << endl;
  cout.precision(4);
  cout << "======================================================" << endl;
}

Member Data Documentation

fBtube

Double_t ShipBellField::fBtube

private

Definition at line 71 of file ShipBellField.h.

fInclTarget

Bool_t ShipBellField::fInclTarget

private

Definition at line 73 of file ShipBellField.h.

fMiddle

Double_t ShipBellField::fMiddle

private

Definition at line 70 of file ShipBellField.h.

fOrient

Int_t ShipBellField::fOrient

private

Definition at line 72 of file ShipBellField.h.

fPeak

Double_t ShipBellField::fPeak

private

Field parameters

Definition at line 69 of file ShipBellField.h.

targetL

Double_t ShipBellField::targetL

private

Definition at line 76 of file ShipBellField.h.

targetXY

Double_t ShipBellField::targetXY

private

Definition at line 74 of file ShipBellField.h.

targetZ0

Double_t ShipBellField::targetZ0

private

Definition at line 75 of file ShipBellField.h.

---

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

• field/ShipBellField.h
• field/ShipBellField.cxx