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 26 of file ShipBellField.h.

Constructor & Destructor Documentation

ShipBellField() [1/3]

ShipBellField::ShipBellField

(

)

Default constructor

Definition at line 29 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.  )

explicit

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 37 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 )

explicit

Constructor from ShipFieldPar

Definition at line 48 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 64 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 75 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 104 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 125 of file ShipBellField.cxx.

{ return 0.; }

IncludeTarget()

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

Definition at line 68 of file ShipBellField.cxx.

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

Print()

void ShipBellField::Print ( )

virtual

Screen output

Definition at line 129 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 62 of file ShipBellField.h.

fInclTarget

Bool_t ShipBellField::fInclTarget

private

Definition at line 64 of file ShipBellField.h.

fMiddle

Double_t ShipBellField::fMiddle

private

Definition at line 61 of file ShipBellField.h.

fOrient

Int_t ShipBellField::fOrient

private

Definition at line 63 of file ShipBellField.h.

fPeak

Double_t ShipBellField::fPeak

private

Field parameters

Definition at line 60 of file ShipBellField.h.

targetL

Double_t ShipBellField::targetL

private

Definition at line 67 of file ShipBellField.h.

targetXY

Double_t ShipBellField::targetXY

private

Definition at line 65 of file ShipBellField.h.

targetZ0

Double_t ShipBellField::targetZ0

private

Definition at line 66 of file ShipBellField.h.

---

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

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