F04ElementField Class Reference

#include <F04ElementField.hh>

Inheritance diagram for F04ElementField:

Public Member Functions
	F04ElementField (const G4ThreeVector, G4LogicalVolume *)
	Constructor. More...
void	Construct ()
	the actual implementation constructs the F04ElementField More...
virtual	~F04ElementField ()
	Destructor. More...
void	SetMaxStep (G4double stp)
	SetMaxStep(G4double) sets the max. step size. More...
G4double	GetMaxStep ()
	GetMaxStep() returns the max. step size. More...
void	SetColor (G4String c)
	SetColor(G4String) sets the color. More...
G4String	GetColor ()
	GetColor() returns the color. More...
void	SetGlobalPoint (const G4double point[4])
bool	IsInBoundingBox (const G4double point[4]) const
virtual void	AddFieldValue (const G4double point[4], G4double field[6]) const =0
virtual G4double	GetLength ()=0
virtual G4double	GetWidth ()=0
virtual G4double	GetHeight ()=0

Static Public Member Functions
static G4VisAttributes *	GetVisAttribute (G4String color)
	GetVisAttribute() returns the appropriate G4VisAttributes. More...

Protected Attributes
G4LogicalVolume *	fVolume
G4AffineTransform	fGlobal2local

Detailed Description

Definition at line 53 of file F04ElementField.hh.

Constructor & Destructor Documentation

F04ElementField::F04ElementField	(	const G4ThreeVector	c,
		G4LogicalVolume *	lv
	)

Constructor.

Definition at line 44 of file F04ElementField.cc.

References F04GlobalField::AddElementField(), test::c, DBL_MAX, fVolume, F04GlobalField::GetObject(), GetVisAttribute(), python.hepunit::m, python.hepunit::MeV, python.hepunit::mm, python.hepunit::ms, G4UserLimits::SetMaxAllowedStep(), G4LogicalVolume::SetUserLimits(), G4UserLimits::SetUserMaxTime(), G4UserLimits::SetUserMaxTrackLength(), G4UserLimits::SetUserMinEkine(), and G4LogicalVolume::SetVisAttributes().

{
  fCenter = c;

  fMinX = fMinY = fMinZ = -DBL_MAX;
  fMaxX = fMaxY = fMaxZ =  DBL_MAX;

  F04GlobalField::GetObject()->AddElementField(this);

  fColor = "1,1,1";

  fUserLimits = new G4UserLimits();

  fVolume = lv;
  fVolume->SetVisAttributes(GetVisAttribute(fColor));

  fMaxStep = 1*mm;

  fUserLimits->SetMaxAllowedStep(fMaxStep);

  fUserLimits->SetUserMaxTrackLength(500.*m);
  fUserLimits->SetUserMaxTime(10*ms);
  fUserLimits->SetUserMinEkine(0.1*MeV);
//  fUserLimits->SetUserMinRange(1*mm);

  fVolume->SetUserLimits(fUserLimits);
}

virtual F04ElementField::~F04ElementField ( )

inlinevirtual

Destructor.

Definition at line 69 of file F04ElementField.hh.

{}

Member Function Documentation

virtual void F04ElementField::AddFieldValue ( const G4double  point[4], G4double  field[6]  ) const

pure virtual

AddFieldValue() will add the field value for this element to field[]. Implementations must be sure to verify that point[] is within the field region, and do nothing if not. point[] is in global coordinates and geant4 units; x,y,z,t. field[] is in geant4 units; Bx,By,Bz,Ex,Ey,Ez. For efficiency, the caller may (but need not) call IsInBoundingBox(point), and only call this function if that returns true.

Implemented in F04SimpleSolenoid, and F04FocusSolenoid.

Referenced by F04GlobalField::GetFieldValue().

void F04ElementField::Construct

(

void

)

the actual implementation constructs the F04ElementField

Definition at line 74 of file F04ElementField.cc.

References fGlobal2local, fVolume, GetHeight(), GetLength(), G4TransportationManager::GetTransportationManager(), GetWidth(), G4Navigator::GetWorldVolume(), G4AffineTransform::Inverse(), local, G4Navigator::LocateGlobalPointAndSetup(), SetGlobalPoint(), G4Navigator::SetWorldVolume(), G4AffineTransform::TransformPoint(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

{
  G4Navigator* theNavigator =
                    G4TransportationManager::GetTransportationManager()->
                                                 GetNavigatorForTracking();

  if (!fNavigator) {
     fNavigator = new G4Navigator();
     if ( theNavigator->GetWorldVolume() )
               fNavigator->SetWorldVolume(theNavigator->GetWorldVolume());
   }

  fNavigator->LocateGlobalPointAndSetup(fCenter,0,false);

  G4TouchableHistoryHandle touchable = fNavigator->
                                         CreateTouchableHistoryHandle();

  G4int depth = touchable->GetHistoryDepth();
  for (G4int i = 0; i<depth; ++i) {
      if(touchable->GetVolume()->GetLogicalVolume() == fVolume)break;
      touchable->MoveUpHistory();
  }

  // set fGlobal2local transform
  fGlobal2local = touchable->GetHistory()->GetTopTransform();

  // set global bounding box
  G4double local[4], global[4];

  G4ThreeVector globalPosition;
  local[3] = 0.0;
  for (int i=0; i<2; ++i) {
      local[0] = (i==0 ? -1.0 : 1.0) * GetWidth()/2.;
      for (int j=0; j<2; ++j) {
          local[1] = (j==0 ? -1.0 : 1.0) * GetHeight()/2.;
          for (int k=0; k<2; ++k) {
              local[2] = (k==0 ? -1.0 : 1.0) * GetLength()/2.;
              G4ThreeVector localPosition(local[0],local[1],local[2]);
              globalPosition =
                         fGlobal2local.Inverse().TransformPoint(localPosition);
              global[0] = globalPosition.x();
              global[1] = globalPosition.y();
              global[2] = globalPosition.z();
              SetGlobalPoint(global);
           }
      }
  }
}

G4String F04ElementField::GetColor ( )

inline

GetColor() returns the color.

Definition at line 90 of file F04ElementField.hh.

{ return fColor; }

virtual G4double F04ElementField::GetHeight ( )

pure virtual

Implemented in F04SimpleSolenoid.

Referenced by Construct().

virtual G4double F04ElementField::GetLength ( )

pure virtual

Implemented in F04SimpleSolenoid.

Referenced by Construct().

G4double F04ElementField::GetMaxStep ( )

inline

GetMaxStep() returns the max. step size.

Definition at line 80 of file F04ElementField.hh.

{ return fMaxStep; }

G4VisAttributes * F04ElementField::GetVisAttribute ( G4String  color )

static

GetVisAttribute() returns the appropriate G4VisAttributes.

Definition at line 125 of file F04ElementField.cc.

References blue, G4cout, G4endl, G4VisAttributes::Invisible, red, and G4VisAttributes::SetDaughtersInvisible().

Referenced by F04ElementField(), and SetColor().

{
   G4VisAttributes* p = NULL;
   if(color.size() > 0 &&
     (isdigit(color.c_str()[0]) || color.c_str()[0] == '.')) {
        G4double red=0.0, green=0.0, blue=0.0;
        if (sscanf(color.c_str(),"%lf,%lf,%lf",&red,&green,&blue) == 3) {
           p = new G4VisAttributes(true,G4Color(red,green,blue));
        } else {
           G4cout << " Invalid color " << color << G4endl;
        }
   }

   if (!p) p = new G4VisAttributes(G4VisAttributes::Invisible);
   p->SetDaughtersInvisible(false);

   return p;
}

virtual G4double F04ElementField::GetWidth ( )

pure virtual

Implemented in F04SimpleSolenoid.

Referenced by Construct().

bool F04ElementField::IsInBoundingBox ( const G4double  point[4] ) const

inline

IsInBoundingBox() returns true if the point is within the global bounding box - global coordinates.

Definition at line 113 of file F04ElementField.hh.

Referenced by F04GlobalField::GetFieldValue().

    {
      if(point[2] < fMinZ || point[2] > fMaxZ) return false;
      if(point[0] < fMinX || point[0] > fMaxX) return false;
      if(point[1] < fMinY || point[1] > fMaxY) return false;
      return true;
    }

void F04ElementField::SetColor ( G4String  c )

inline

SetColor(G4String) sets the color.

Definition at line 83 of file F04ElementField.hh.

References test::c, fVolume, GetVisAttribute(), and G4LogicalVolume::SetVisAttributes().

Referenced by F04GlobalField::ConstructField().

    {
      fColor = c;
      fVolume->SetVisAttributes(GetVisAttribute(fColor));
    }

void F04ElementField::SetGlobalPoint ( const G4double  point[4] )

inline

SetGlobalPoint() ensures that the point is within the global bounding box of this ElementField's global coordinates. Normally called 8 times for the corners of the local bounding box, after a local->global coordinate transform. If never called, the global bounding box is infinite. BEWARE: if called only once, the bounding box is just a point.

Definition at line 101 of file F04ElementField.hh.

References DBL_MAX.

Referenced by Construct().

    {
      if(fMinX == -DBL_MAX || fMinX > point[0]) fMinX = point[0];
      if(fMinY == -DBL_MAX || fMinY > point[1]) fMinY = point[1];
      if(fMinZ == -DBL_MAX || fMinZ > point[2]) fMinZ = point[2];
      if(fMaxX ==  DBL_MAX || fMaxX < point[0]) fMaxX = point[0];
      if(fMaxY ==  DBL_MAX || fMaxY < point[1]) fMaxY = point[1];
      if(fMaxZ ==  DBL_MAX || fMaxZ < point[2]) fMaxZ = point[2];
    }

void F04ElementField::SetMaxStep ( G4double  stp )

inline

SetMaxStep(G4double) sets the max. step size.

Definition at line 72 of file F04ElementField.hh.

References fVolume, G4UserLimits::SetMaxAllowedStep(), and G4LogicalVolume::SetUserLimits().

Referenced by F04GlobalField::ConstructField().

    {
      fMaxStep = stp;
      fUserLimits->SetMaxAllowedStep(fMaxStep);
      fVolume->SetUserLimits(fUserLimits);
    }

Field Documentation

G4AffineTransform F04ElementField::fGlobal2local

protected

Definition at line 140 of file F04ElementField.hh.

Referenced by F04FocusSolenoid::AddFieldValue(), F04SimpleSolenoid::AddFieldValue(), and Construct().

G4LogicalVolume* F04ElementField::fVolume

protected

Definition at line 138 of file F04ElementField.hh.

Referenced by Construct(), F04ElementField(), F04SimpleSolenoid::F04SimpleSolenoid(), SetColor(), and SetMaxStep().

---

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

• F04ElementField.hh
• F04ElementField.cc