#include <G4VFieldModel.hh>

Inheritance diagram for G4VFieldModel:

Public Types
enum	Representation { fullArrow , lightArrow }

Public Member Functions
	G4VFieldModel (const G4String &typeOfField, const G4String &symbol="", const G4VisExtent &extentForField=G4VisExtent(), const std::vector< G4PhysicalVolumesSearchScene::Findings > &pvFindings=std::vector< G4PhysicalVolumesSearchScene::Findings >(), G4int nDataPointsPerHalfScene=10, Representation representation=Representation::fullArrow, G4int arrow3DLineSegmentsPerCircle=6)

virtual	~G4VFieldModel ()

virtual void	DescribeYourselfTo (G4VGraphicsScene &sceneHandler)

Public Member Functions inherited from G4VModel
	G4VModel (const G4Transform3D &modelTransformation=G4Transform3D(), const G4ModelingParameters *=0)

virtual	~G4VModel ()

virtual void	DescribeYourselfTo (G4VGraphicsScene &)=0

const G4ModelingParameters *	GetModelingParameters () const

const G4String &	GetType () const

virtual G4String	GetCurrentDescription () const

virtual G4String	GetCurrentTag () const

const G4VisExtent &	GetExtent () const

const G4VisExtent &	GetTransformedExtent () const

const G4String &	GetGlobalDescription () const

const G4String &	GetGlobalTag () const

const G4Transform3D &	GetTransformation () const

void	SetModelingParameters (const G4ModelingParameters *)

void	SetExtent (const G4VisExtent &)

void	SetType (const G4String &)

void	SetGlobalDescription (const G4String &)

void	SetGlobalTag (const G4String &)

void	SetTransformation (const G4Transform3D &)

virtual G4bool	Validate (G4bool warn=true)

Protected Member Functions
virtual void	GetFieldAtLocation (const G4Field *field, const G4Point3D &position, G4double time, G4Point3D &result) const =0

Additional Inherited Members
Protected Attributes inherited from G4VModel
G4String	fType

G4String	fGlobalTag

G4String	fGlobalDescription

G4VisExtent	fExtent

G4Transform3D	fTransform

const G4ModelingParameters *	fpMP

Detailed Description

Definition at line 46 of file G4VFieldModel.hh.

Member Enumeration Documentation

◆ Representation

enum G4VFieldModel::Representation

Enumerator
fullArrow
lightArrow

Definition at line 50 of file G4VFieldModel.hh.

50{fullArrow, lightArrow};

G4VFieldModel::lightArrow

@ lightArrow

Definition: G4VFieldModel.hh:50

G4VFieldModel::fullArrow

@ fullArrow

Definition: G4VFieldModel.hh:50

Constructor & Destructor Documentation

◆ G4VFieldModel()

G4VFieldModel::G4VFieldModel	(	const G4String &	typeOfField,
		const G4String &	symbol = `""`,
		const G4VisExtent &	extentForField = `G4VisExtent()`,
		const std::vector< G4PhysicalVolumesSearchScene::Findings > &	pvFindings = `std::vector<G4PhysicalVolumesSearchScene::Findings>()`,
		G4int	nDataPointsPerHalfScene = `10`,
		Representation	representation = `Representation::fullArrow`,
		G4int	arrow3DLineSegmentsPerCircle = `6`
	)

Definition at line 60 of file G4VFieldModel.cc.

: fExtentForField(extentForField)
, fPVFindings(pvFindings)
, fNDataPointsPerMaxHalfScene(nDataPointsPerMaxHalfScene)
, fRepresentation(representation)
, fArrow3DLineSegmentsPerCircle(arrow3DLineSegmentsPerCircle)
, fTypeOfField(typeOfField)
, fArrowPrefix(symbol)
{
  fType = "G4"+typeOfField+"FieldModel";
  fGlobalTag = fType;
 
  std::ostringstream oss;
  oss << ':' << fNDataPointsPerMaxHalfScene
  << ':' << fArrow3DLineSegmentsPerCircle;
  if (fExtentForField == G4VisExtent::GetNullExtent()) {
    oss << " whole scene";
  } else {
    oss
    << ':' << fExtentForField.GetXmin()
    << ':' << fExtentForField.GetXmax()
    << ':' << fExtentForField.GetYmin()
    << ':' << fExtentForField.GetYmax()
    << ':' << fExtentForField.GetZmin()
    << ':' << fExtentForField.GetZmax();
  }
  for (const auto& findings: fPVFindings) {
    oss
    << ',' << findings.fpFoundPV->GetName()
    << ':' << findings.fFoundPVCopyNo;
  }
  if (fRepresentation == Representation::fullArrow) {
    oss << " full arrow";
  } else if (fRepresentation == Representation::lightArrow) {
    oss << " light arrow";
  }
 
  fGlobalDescription = fType + oss.str();
}

◆ ~G4VFieldModel()

G4VFieldModel::~G4VFieldModel ( )

virtual

Definition at line 58 of file G4VFieldModel.cc.

58{;}

Member Function Documentation

◆ DescribeYourselfTo()

void G4VFieldModel::DescribeYourselfTo ( G4VGraphicsScene & sceneHandler )

virtual

Implements G4VModel.

Definition at line 109 of file G4VFieldModel.cc.

                                                                     {
//  G4cout << "G4VFieldModel::DescribeYourselfTo" << G4endl;
 
  G4TransportationManager* tMgr =
  G4TransportationManager::GetTransportationManager();
  assert(tMgr);
  G4Navigator* navigator = tMgr->GetNavigatorForTracking();
  assert(navigator);
 
  G4FieldManager* globalFieldMgr = tMgr->GetFieldManager();
  const G4Field* globalField = 0;
  const G4String intro = "G4VFieldModel::DescribeYourselfTo: ";
  if (globalFieldMgr) {
    if (globalFieldMgr->DoesFieldExist()) {
      globalField = globalFieldMgr->GetDetectorField();
      if (!globalField) {
        static G4bool warned = false;
        if (!warned) {
          G4cout << intro << "Null global field pointer." << G4endl;
          warned = true;
        }
      }
    }
  } else {
    static G4bool warned = false;
    if (!warned) {
      G4cout << intro << "No global field manager." << G4endl;
      warned = true;
    }
  }
 
  G4VisExtent sceneExtent = sceneHandler.GetExtent();
  const G4double& xMin = sceneExtent.GetXmin();
  const G4double& yMin = sceneExtent.GetYmin();
  const G4double& zMin = sceneExtent.GetZmin();
  const G4double& xMax = sceneExtent.GetXmax();
  const G4double& yMax = sceneExtent.GetYmax();
  const G4double& zMax = sceneExtent.GetZmax();
  const G4double xHalfScene = 0.5 * (xMax - xMin);
  const G4double yHalfScene = 0.5 * (yMax - yMin);
  const G4double zHalfScene = 0.5 * (zMax - zMin);
  const G4double xSceneCentre = 0.5 * (xMax + xMin);
  const G4double ySceneCentre = 0.5 * (yMax + yMin);
  const G4double zSceneCentre = 0.5 * (zMax + zMin);
  const G4double maxHalfScene =
  std::max(xHalfScene,std::max(yHalfScene,zHalfScene));
  if (maxHalfScene <= 0.) {
    G4cout << "Scene extent non-positive." << G4endl;
    return;
  }
 
  // Constants
  const G4double interval = maxHalfScene / fNDataPointsPerMaxHalfScene;
  const G4int nDataPointsPerXHalfScene = G4int(xHalfScene / interval);
  const G4int nDataPointsPerYHalfScene = G4int(yHalfScene / interval);
  const G4int nDataPointsPerZHalfScene = G4int(zHalfScene / interval);
  const G4int nXSamples = 2 * nDataPointsPerXHalfScene + 1;
  const G4int nYSamples = 2 * nDataPointsPerYHalfScene + 1;
  const G4int nZSamples = 2 * nDataPointsPerZHalfScene + 1;
  const G4int nSamples = nXSamples * nYSamples * nZSamples;
  const G4double arrowLengthMax = 0.8 * interval;
 
  // Working vectors for field values, etc.
  std::vector<G4Point3D> Field(nSamples);          // Initialises to (0,0,0)
  std::vector<G4Point3D> xyz(nSamples);            // Initialises to (0,0,0)
  G4double FieldMagnitudeMax = -std::numeric_limits<G4double>::max();
 
  // Get field values and ascertain maximum field.
  for (G4int i = 0; i < nXSamples; i++) {
    G4double x = xSceneCentre + (i - nDataPointsPerXHalfScene) * interval;
 
    for (G4int j = 0; j < nYSamples; j++) {
      G4double y = ySceneCentre + (j - nDataPointsPerYHalfScene) * interval;
 
      for (G4int k = 0; k < nZSamples; k++) {
        G4double z = zSceneCentre + (k - nDataPointsPerZHalfScene) * interval;
 
        // Calculate indices into working vectors
        const G4int ijk = i * nYSamples * nZSamples + j * nZSamples + k;
    xyz[ijk].set(x,y,z);
 
        G4ThreeVector pos(x,y,z);
 
        // Check if point is in extent for field
        if (fExtentForField != G4VisExtent::GetNullExtent()) {
          const auto& ext = fExtentForField;  // Alias
          if (x < ext.GetXmin() || x > ext.GetXmax() ||
              y < ext.GetYmin() || y > ext.GetYmax() ||
              z < ext.GetZmin() || z > ext.GetZmax())
            continue;
        }
 
        // Check if point is in findings
        if (!fPVFindings.empty()) {
          G4bool isInPV = false;
          for (const auto& findings: fPVFindings) {
            G4VPhysicalVolume* pv = findings.fpFoundPV;
            G4int copyNo = findings.fFoundPVCopyNo;
            G4VSolid* solid = pv->GetLogicalVolume()->GetSolid();
            G4PVParameterised* pvParam = dynamic_cast<G4PVParameterised*>(pv);
            if (pvParam) {
              auto* param = pvParam->GetParameterisation();
              solid = param->ComputeSolid(copyNo,pvParam);
              solid->ComputeDimensions(param,copyNo,pvParam);
            }
            // Transform point to local coordinate system
            const auto& transform = findings.fFoundObjectTransformation;
            auto rotation = transform.getRotation();
            auto translation = transform.getTranslation();
            G4ThreeVector lPos = pos; lPos -= translation; lPos.transform(rotation.invert());
            if (solid->Inside(lPos)==kInside) {
              isInPV = true;
              break;
            }
          }
          if (!isInPV) continue;
        }
        // Point is in findings - or there were no findings
 
        // Find volume and field at this location.
        const G4VPhysicalVolume* pPV =
        navigator->LocateGlobalPointAndSetup(pos,0,false,true);
        const G4Field* field = globalField;
        if (pPV) {
          // Get logical volume.
          const G4LogicalVolume* pLV = pPV->GetLogicalVolume();
          if (pLV) {
            // Value for Region, if any, overrides
            G4Region* pRegion = pLV->GetRegion();
            if (pRegion) {
              G4FieldManager* pRegionFieldMgr = pRegion->GetFieldManager();
              if (pRegionFieldMgr) {
                field = pRegionFieldMgr->GetDetectorField();
                // G4cout << "Region with field" << G4endl;
              }
            }
            // 'Local' value from logical volume, if any, overrides
            G4FieldManager* pLVFieldMgr = pLV->GetFieldManager();
            if (pLVFieldMgr) {
              field = pLVFieldMgr->GetDetectorField();
              // G4cout << "Logical volume with field" << G4endl;
            }
          }
        }
 
    G4double time = 0.; // FIXME:  Can we get event time in some way?
 
    // Subclasses will have implemented this for their own field
    GetFieldAtLocation(field, xyz[ijk], time, Field[ijk]);
 
    G4double mag = Field[ijk].mag();
    if (mag > FieldMagnitudeMax) FieldMagnitudeMax = mag;
      } // for (k, z
    }   // for (j, y
  } // for (i, x
 
  if (FieldMagnitudeMax <= 0.) {
    G4cout << "No " << fTypeOfField << " field in this extent." << G4endl;
    return;
  }
 
  for (G4int i = 0; i < nSamples; i++) {
    const G4double Fmag = Field[i].mag();
    const G4double f = Fmag / FieldMagnitudeMax;
    if (f <= 0.) continue;  // Skip zero field locations
 
    G4double red = 0., green = 0., blue = 0., alpha = 1.;
    if (f < 0.5) {  // Linear colour scale: 0->0.5->1 is red->green->blue.
      green = 2. * f;
      red = 2. * (0.5 - f);
    } else {
      blue = 2. * (f - 0.5);
      green = 2. * (1.0 - f);
    }
    const G4Colour arrowColour(red,green,blue,alpha);
 
    // Very small arrows are difficult to see. Better to draw a line.
    G4bool drawAsLine = false;
    switch (fRepresentation) {
      case Representation::fullArrow:
        if (f < 0.1) {
          drawAsLine = true;
        }
        break;
      case Representation::lightArrow:
        drawAsLine = true;
        break;
      default:
        break;
    }
 
    // Head of arrow depends on field direction and strength...
    G4double arrowLength = arrowLengthMax * f;
    // ...but limit the length so it's visible.
    if (f < 0.01) arrowLength = arrowLengthMax * 0.01;
    const G4Point3D head = xyz[i] + arrowLength*Field[i]/Fmag;
 
    if (drawAsLine) {
      G4Polyline FArrowLite;
      G4VisAttributes va(arrowColour);
      va.SetLineWidth(2.);
      FArrowLite.SetVisAttributes(va);
      FArrowLite.push_back(xyz[i]);
      FArrowLite.push_back(head);
      sceneHandler.BeginPrimitives();
      sceneHandler.AddPrimitive(FArrowLite);
      sceneHandler.EndPrimitives();
    } else {
      G4ArrowModel FArrow(xyz[i].x(), xyz[i].y(), xyz[i].z(),
                          head.x(), head.y(), head.z(),
                          arrowLength/5, arrowColour,
                          fArrowPrefix+"Field",
                          fArrow3DLineSegmentsPerCircle);
      FArrow.DescribeYourselfTo(sceneHandler);
    }
  } // for (i, nSamples
}

◆ GetFieldAtLocation()

virtual void G4VFieldModel::GetFieldAtLocation	(	const G4Field *	field,
		const G4Point3D &	position,
		G4double	time,
		G4Point3D &	result
	)		const

protectedpure virtual

Implemented in G4ElectricFieldModel, and G4MagneticFieldModel.

Referenced by DescribeYourselfTo().

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

Public Types

Public Member Functions

Protected Member Functions

Additional Inherited Members