MedicalBeam Class Reference

#include "MedicalBeam.hh"

Inheritance diagram for MedicalBeam:

Public Types
enum	FieldShape { SQUARE =0, CIRCLE }

Public Member Functions
	MedicalBeam ()
	~MedicalBeam ()
void	SetParticleDefinition (G4ParticleDefinition *pd)
const G4ParticleDefinition *	GetParticleDefinition () const
void	SetKineticE (G4double e)
G4double	GetKineticE () const
void	SetSourcePosition (const G4ThreeVector &pos)
G4ThreeVector	GetSourcePosition () const
void	SetFieldShape (FieldShape shape)
FieldShape	GetFieldShape () const
void	SetSSD (G4double ssd)
G4double	GetSSD () const
void	SetFieldXY (G4double fx, G4double fy)
G4double	GetFieldX () const
G4double	GetFieldY () const
void	SetFieldR (G4double r)
G4double	GetFieldR () const
virtual void	GeneratePrimaries (G4Event *anEvent)

Protected Member Functions
G4ThreeVector	GenerateBeamDirection () const

Protected Attributes
G4ParticleDefinition *	particle
G4double	kineticE
G4ThreeVector	sourcePosition
G4double	SSD
FieldShape	fieldShape
G4double	fieldXY [2]
G4double	fieldR

Detailed Description

Definition at line 45 of file MedicalBeam.hh.

Member Enumeration Documentation

enum MedicalBeam::FieldShape

Enumerator
SQUARE
CIRCLE

Definition at line 47 of file MedicalBeam.hh.

{ SQUARE=0, CIRCLE };

Constructor & Destructor Documentation

MedicalBeam::MedicalBeam

(

)

Definition at line 48 of file MedicalBeam.cc.

References fieldXY.

  : particle(0), 
    kineticE(1.*MeV),
    sourcePosition(G4ThreeVector()),
    SSD(1.*m), 
    fieldShape(MedicalBeam::SQUARE),
    fieldR(10.*cm)
{
  fieldXY[0]= fieldXY[1]= 10.*cm;
}

MedicalBeam::~MedicalBeam

(

)

Definition at line 62 of file MedicalBeam.cc.

{
}

Member Function Documentation

G4ThreeVector MedicalBeam::GenerateBeamDirection ( ) const

protected

Definition at line 69 of file MedicalBeam.cc.

References util::abs(), fieldR, fieldShape, fieldXY, sqr(), SQUARE, SSD, x, y, and z.

Referenced by GeneratePrimaries().

{
  // uniform distribution in a limitted solid angle
  G4double dr;
  if(fieldShape == MedicalBeam::SQUARE) {
    dr= std::sqrt(sqr(fieldXY[0]/2.)+sqr(fieldXY[1]/2.));
  } else {
    dr= fieldR;
  }

  G4double sin0= dr/SSD;
  G4double cos0= std::sqrt(1.-sqr(sin0));

  G4double dcos, dsin, dphi, z;

  G4double x= DBL_MAX;
  G4double y= DBL_MAX;

  G4double xmax, ymax;
  if(fieldShape == MedicalBeam::SQUARE) {
    xmax= fieldXY[0]/2./SSD;
    ymax= fieldXY[1]/2./SSD;
  } else {
    xmax= ymax= DBL_MAX-1.;
  }

  while(! (std::abs(x)< xmax && std::abs(y)< ymax) ) {
    dcos= RandFlat::shoot(cos0, 1.);
    dsin= std::sqrt(1.-sqr(dcos));
    dphi= RandFlat::shoot(0., twopi);

    x= std::cos(dphi)*dsin*dcos;
    y= std::sin(dphi)*dsin*dcos;
  }
  z= dcos;

  return G4ThreeVector(x,y,z);
}

void MedicalBeam::GeneratePrimaries ( G4Event *  anEvent )

virtual

Definition at line 111 of file MedicalBeam.cc.

References GenerateBeamDirection(), kineticE, particle, sourcePosition, and sqr().

{
  if(particle==0) return;

  // create a new vertex
  G4PrimaryVertex* vertex= new G4PrimaryVertex(sourcePosition, 0.*ns);

  // momentum
  G4double mass= particle-> GetPDGMass();
  G4double p= std::sqrt(sqr(mass+kineticE)-sqr(mass));
  G4ThreeVector pmon= p*GenerateBeamDirection();
  G4PrimaryParticle* primary= new G4PrimaryParticle(particle, 
                                                    pmon.x(), 
                                                    pmon.y(), 
                                                    pmon.z());
  // set primary to vertex
  vertex-> SetPrimary(primary);

  // set vertex to event
  anEvent-> AddPrimaryVertex(vertex);
}

G4double MedicalBeam::GetFieldR ( ) const

inline

Definition at line 139 of file MedicalBeam.hh.

References fieldR.

Referenced by BOOST_PYTHON_MODULE().

{ return fieldR; }

MedicalBeam::FieldShape MedicalBeam::GetFieldShape ( ) const

inline

Definition at line 118 of file MedicalBeam.hh.

References fieldShape.

Referenced by BOOST_PYTHON_MODULE().

{ return fieldShape; }

G4double MedicalBeam::GetFieldX ( ) const

inline

Definition at line 130 of file MedicalBeam.hh.

References fieldXY.

Referenced by BOOST_PYTHON_MODULE().

{ return fieldXY[0]; }

G4double MedicalBeam::GetFieldY ( ) const

inline

Definition at line 133 of file MedicalBeam.hh.

References fieldXY.

Referenced by BOOST_PYTHON_MODULE().

{ return fieldXY[1]; }

G4double MedicalBeam::GetKineticE ( ) const

inline

Definition at line 106 of file MedicalBeam.hh.

References kineticE.

Referenced by BOOST_PYTHON_MODULE().

{ return kineticE; }

const G4ParticleDefinition * MedicalBeam::GetParticleDefinition ( ) const

inline

Definition at line 100 of file MedicalBeam.hh.

References particle.

{ return particle; }

G4ThreeVector MedicalBeam::GetSourcePosition ( ) const

inline

Definition at line 112 of file MedicalBeam.hh.

References sourcePosition.

Referenced by BOOST_PYTHON_MODULE().

{ return sourcePosition; }

G4double MedicalBeam::GetSSD ( ) const

inline

Definition at line 124 of file MedicalBeam.hh.

References SSD.

Referenced by BOOST_PYTHON_MODULE().

{ return SSD; }

void MedicalBeam::SetFieldR ( G4double  r )

inline

Definition at line 136 of file MedicalBeam.hh.

References fieldR, and r.

Referenced by BOOST_PYTHON_MODULE().

{ fieldR= r; }

void MedicalBeam::SetFieldShape ( MedicalBeam::FieldShape  shape )

inline

Definition at line 115 of file MedicalBeam.hh.

References fieldShape.

Referenced by BOOST_PYTHON_MODULE().

{ fieldShape= shape; }

void MedicalBeam::SetFieldXY ( G4double  fx, G4double  fy  )

inline

Definition at line 127 of file MedicalBeam.hh.

References fieldXY.

{ fieldXY[0]= fx; fieldXY[1]= fy; }

void MedicalBeam::SetKineticE ( G4double  e )

inline

Definition at line 103 of file MedicalBeam.hh.

References e, and kineticE.

Referenced by BOOST_PYTHON_MODULE().

{ kineticE= e; }

void MedicalBeam::SetParticleDefinition ( G4ParticleDefinition *  pd )

inline

Definition at line 97 of file MedicalBeam.hh.

References particle.

{ particle= pd; }

void MedicalBeam::SetSourcePosition ( const G4ThreeVector &  pos )

inline

Definition at line 109 of file MedicalBeam.hh.

References sourcePosition.

Referenced by BOOST_PYTHON_MODULE().

{ sourcePosition= pos; }

void MedicalBeam::SetSSD ( G4double  ssd )

inline

Definition at line 121 of file MedicalBeam.hh.

References SSD.

Referenced by BOOST_PYTHON_MODULE().

{ SSD= ssd; }

Member Data Documentation

G4double MedicalBeam::fieldR

protected

Definition at line 57 of file MedicalBeam.hh.

Referenced by GenerateBeamDirection(), GetFieldR(), and SetFieldR().

FieldShape MedicalBeam::fieldShape

protected

Definition at line 55 of file MedicalBeam.hh.

Referenced by GenerateBeamDirection(), GetFieldShape(), and SetFieldShape().

G4double MedicalBeam::fieldXY[2]

protected

Definition at line 56 of file MedicalBeam.hh.

Referenced by GenerateBeamDirection(), GetFieldX(), GetFieldY(), MedicalBeam(), and SetFieldXY().

G4double MedicalBeam::kineticE

protected

Definition at line 51 of file MedicalBeam.hh.

Referenced by GeneratePrimaries(), GetKineticE(), and SetKineticE().

G4ParticleDefinition* MedicalBeam::particle

protected

Definition at line 50 of file MedicalBeam.hh.

Referenced by GeneratePrimaries(), GetParticleDefinition(), and SetParticleDefinition().

G4ThreeVector MedicalBeam::sourcePosition

protected

Definition at line 52 of file MedicalBeam.hh.

Referenced by GeneratePrimaries(), GetSourcePosition(), and SetSourcePosition().

G4double MedicalBeam::SSD

protected

Definition at line 54 of file MedicalBeam.hh.

Referenced by GenerateBeamDirection(), GetSSD(), and SetSSD().

---

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

• geant4.10.06.p01/environments/g4py/site-modules/primaries/MedicalBeam/MedicalBeam.hh
• geant4.10.06.p01/environments/g4py/site-modules/primaries/MedicalBeam/MedicalBeam.cc