CCalMagneticField Class Reference

#include <CCalMagneticField.hh>

Inheritance diagram for CCalMagneticField:

Collaboration diagram for CCalMagneticField:

Public Member Functions
	CCalMagneticField (const G4String &name)
	~CCalMagneticField ()
void	MagneticField (const double Point[3], double Bfield[3]) const
CLHEP::Hep3Vector	MagneticField (const CLHEP::Hep3Vector Point) const
virtual void	GetFieldValue (const double Point[3], double *Bfield) const
G4double	GetConstantFieldvalue () const
Public Member Functions inherited from G4MagneticField
	G4MagneticField ()
virtual	~G4MagneticField ()
	G4MagneticField (const G4MagneticField &r)
G4MagneticField &	operator= (const G4MagneticField &p)
G4bool	DoesFieldChangeEnergy () const
virtual void	GetFieldValue (const G4double Point[4], G4double *Bfield) const =0
Public Member Functions inherited from G4ElectroMagneticField
	G4ElectroMagneticField ()
virtual	~G4ElectroMagneticField ()
	G4ElectroMagneticField (const G4ElectroMagneticField &r)
G4ElectroMagneticField &	operator= (const G4ElectroMagneticField &p)
Public Member Functions inherited from G4Field
	G4Field (G4bool gravityOn=false)
	G4Field (const G4Field &)
virtual	~G4Field ()
G4Field &	operator= (const G4Field &p)
G4bool	IsGravityActive () const
void	SetGravityActive (G4bool OnOffFlag)
virtual G4Field *	Clone () const

Protected Member Functions
G4FieldManager *	GetGlobalFieldManager ()

Detailed Description

Definition at line 38 of file CCalMagneticField.hh.

Constructor & Destructor Documentation

CCalMagneticField::CCalMagneticField

(

const G4String &

name

)

Definition at line 42 of file CCalMagneticField.cc.

                                                             :
  fval(0), pos(0), slope(0), intercept(0) 
{
#ifdef debug
  fVerbosity = 1;
#else
  fVerbosity = 0;
#endif

  //Let's open the file
  G4cout << " ==> Opening file " << filename << " to read magnetic field..."
     << G4endl;
  G4String pathName = getenv("CCAL_GLOBALPATH");
  std::ifstream is;
  bool ok = openGeomFile(is, pathName, filename);
  
  if (ok) {
    findDO(is, G4String("FLDM"));
    is >> fval >> npts >> xoff;

    if (fVerbosity)
      G4cout << "Field value " << fval << " # points " << npts << 
    " offset in x "
         << xoff*mm << G4endl;

    if (npts > 0) {
      pos       = new G4double[npts];
      slope     = new G4double[npts];
      intercept = new G4double[npts];

      for (G4int i = 0; i < npts; i++) {
    is >> pos[i] >> slope[i] >> intercept[i];
    if (fVerbosity)
      G4cout << tab << "Position " << i << " " << pos[i] << " Slope "
         << slope[i] << " Intercept " << intercept[i] << G4endl;
      }
    }

    // Close the file
    G4cout << " ==> Closing file " << filename << G4endl;
    is.close();
  }
}

Here is the call graph for this function:

CCalMagneticField::~CCalMagneticField

(

)

Definition at line 88 of file CCalMagneticField.cc.

                                      {
  if (pos)
    delete[] pos;
  if (slope)
    delete[] slope;
  if (intercept)
    delete[] intercept;
}

Member Function Documentation

G4double CCalMagneticField::GetConstantFieldvalue ( ) const

inline

Definition at line 47 of file CCalMagneticField.hh.

{return fval;}

Here is the caller graph for this function:

void CCalMagneticField::GetFieldValue ( const double  Point[3], double *  Bfield  ) const

virtual

Definition at line 152 of file CCalMagneticField.cc.

                                                                        {
  CCalMagneticField::MagneticField(x, B);
}

Here is the call graph for this function:

G4FieldManager* CCalMagneticField::GetGlobalFieldManager ( )

protected

void CCalMagneticField::MagneticField	(	const double	Point[3],
		double	Bfield[3]
	)		const

Definition at line 100 of file CCalMagneticField.cc.

{
  G4int i=0;
  for (i=0; i<2; i++) {
    B[i]   = 0*kilogauss;
  }

  G4double m1=0;
  G4double c1=0;
  G4double xnew = x[0]/mm + xoff;
  if (npts > 0) {
    for (i=0; i<npts; i++) {
      if (xnew > pos[i]*mm) {
    m1 = slope[i];
        c1 = intercept[i];
      }
    }
  }
  G4double scor = c1 + m*xnew;
  if (scor < 0.) scor = 0.;
  if (scor > 1.) scor = 1.0;

  B[2] = scor*fval*kilogauss;
  if (fVerbosity)
    {

      G4cout << "Field at x: " << x[0]/mm << "mm (" << xnew << ") = " << 
    B[2]/tesla
         << "T (m = " << m1 << ", c = " << 
    c1 << ", scale = " << scor << ")"
         << G4endl;
    }
}

Here is the caller graph for this function:

CLHEP::Hep3Vector CCalMagneticField::MagneticField

(

const CLHEP::Hep3Vector

Point

)

const

Definition at line 136 of file CCalMagneticField.cc.

                                                 {

  G4double x[3],B[3];
  CLHEP::Hep3Vector v;
  
  x[0] = point.x();
  x[1] = point.y();
  x[2] = point.z();
  CCalMagneticField::MagneticField(x, B);
  v.setX(B[0]);   
  v.setY(B[1]);   
  v.setZ(B[2]);   
  return v;
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p02/examples/advanced/composite_calorimeter/include/CCalMagneticField.hh
• source/geant4.10.03.p02/examples/advanced/composite_calorimeter/src/CCalMagneticField.cc