MedicalBeam Class Reference

#include <MedicalBeam.hh>

Inheritance diagram for MedicalBeam:

Public Types
enum	FieldShape { SQUARE =0 , CIRCLE }

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

Protected Member Functions
G4ThreeVector	GenerateBeamDirection () const

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

Detailed Description

Definition at line 45 of file MedicalBeam.hh.

Member Enumeration Documentation

FieldShape

enum MedicalBeam::FieldShape

Enumerator
SQUARE
CIRCLE

Definition at line 47 of file MedicalBeam.hh.

{ SQUARE=0, CIRCLE };

Constructor & Destructor Documentation

MedicalBeam()

MedicalBeam::MedicalBeam

(

)

Definition at line 48 of file MedicalBeam.cc.

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

References cm, and fieldXY.

~MedicalBeam()

MedicalBeam::~MedicalBeam

(

)

Definition at line 62 of file MedicalBeam.cc.

{
}

Member Function Documentation

GenerateBeamDirection()

G4ThreeVector MedicalBeam::GenerateBeamDirection ( ) const

protected

Definition at line 69 of file MedicalBeam.cc.

{
  // 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 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.;
  }
 
  G4double dcos{0.0};
  G4double dsin{0.0};
  G4double dphi{0.0};
 
  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;
  }
 
  G4double z= dcos;
 
  return G4ThreeVector(x,y,z);
}

References DBL_MAX, G4INCL::DeJongSpin::shoot(), sqr(), SQUARE, and twopi.

GeneratePrimaries()

void MedicalBeam::GeneratePrimaries ( G4Event *  anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 114 of file MedicalBeam.cc.

{
  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);
}

References ns, sqr(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

GetFieldR()

G4double MedicalBeam::GetFieldR ( ) const

inline

Definition at line 139 of file MedicalBeam.hh.

{ return fieldR; }

References fieldR.

Referenced by BOOST_PYTHON_MODULE().

GetFieldShape()

MedicalBeam::FieldShape MedicalBeam::GetFieldShape ( ) const

inline

Definition at line 118 of file MedicalBeam.hh.

{ return fieldShape; }

References fieldShape.

Referenced by BOOST_PYTHON_MODULE().

GetFieldX()

G4double MedicalBeam::GetFieldX ( ) const

inline

Definition at line 130 of file MedicalBeam.hh.

{ return fieldXY[0]; }

References fieldXY.

Referenced by BOOST_PYTHON_MODULE().

GetFieldY()

G4double MedicalBeam::GetFieldY ( ) const

inline

Definition at line 133 of file MedicalBeam.hh.

{ return fieldXY[1]; }

References fieldXY.

Referenced by BOOST_PYTHON_MODULE().

GetKineticE()

G4double MedicalBeam::GetKineticE ( ) const

inline

Definition at line 106 of file MedicalBeam.hh.

{ return kineticE; }

References kineticE.

Referenced by BOOST_PYTHON_MODULE().

GetParticleDefinition()

const G4ParticleDefinition * MedicalBeam::GetParticleDefinition ( ) const

inline

Definition at line 100 of file MedicalBeam.hh.

{ return particle; }

References particle.

GetSourcePosition()

G4ThreeVector MedicalBeam::GetSourcePosition ( ) const

inline

Definition at line 112 of file MedicalBeam.hh.

{ return sourcePosition; }

References sourcePosition.

Referenced by BOOST_PYTHON_MODULE().

GetSSD()

G4double MedicalBeam::GetSSD ( ) const

inline

Definition at line 124 of file MedicalBeam.hh.

{ return SSD; }

References SSD.

Referenced by BOOST_PYTHON_MODULE().

SetFieldR()

void MedicalBeam::SetFieldR ( G4double  r )

inline

Definition at line 136 of file MedicalBeam.hh.

{ fieldR= r; }

References fieldR.

Referenced by BOOST_PYTHON_MODULE().

SetFieldShape()

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

inline

Definition at line 115 of file MedicalBeam.hh.

{ fieldShape= shape; }

References fieldShape.

Referenced by BOOST_PYTHON_MODULE().

SetFieldXY()

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

inline

Definition at line 127 of file MedicalBeam.hh.

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

References fieldXY.

SetKineticE()

void MedicalBeam::SetKineticE ( G4double  e )

inline

Definition at line 103 of file MedicalBeam.hh.

{ kineticE= e; }

References kineticE.

Referenced by BOOST_PYTHON_MODULE().

SetParticleDefinition()

void MedicalBeam::SetParticleDefinition ( G4ParticleDefinition *  pd )

inline

Definition at line 97 of file MedicalBeam.hh.

{ particle= pd; }

References particle.

SetSourcePosition()

void MedicalBeam::SetSourcePosition ( const G4ThreeVector &  pos )

inline

Definition at line 109 of file MedicalBeam.hh.

{ sourcePosition= pos; }

References pos, and sourcePosition.

Referenced by BOOST_PYTHON_MODULE().

SetSSD()

void MedicalBeam::SetSSD ( G4double  ssd )

inline

Definition at line 121 of file MedicalBeam.hh.

{ SSD= ssd; }

References SSD.

Referenced by BOOST_PYTHON_MODULE().

Field Documentation

fieldR

G4double MedicalBeam::fieldR

protected

Definition at line 57 of file MedicalBeam.hh.

Referenced by GetFieldR(), and SetFieldR().

fieldShape

FieldShape MedicalBeam::fieldShape

protected

Definition at line 55 of file MedicalBeam.hh.

Referenced by GetFieldShape(), and SetFieldShape().

fieldXY

G4double MedicalBeam::fieldXY[2]

protected

Definition at line 56 of file MedicalBeam.hh.

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

kineticE

G4double MedicalBeam::kineticE

protected

Definition at line 51 of file MedicalBeam.hh.

Referenced by GetKineticE(), mcscore.MCParticle::printout(), and SetKineticE().

particle

G4ParticleDefinition* MedicalBeam::particle

protected

Definition at line 50 of file MedicalBeam.hh.

Referenced by emcalc_gui.MainWindow::cb_select_particle(), emcalc_gui.MainWindow::cb_show_plot(), GetParticleDefinition(), and SetParticleDefinition().

sourcePosition

G4ThreeVector MedicalBeam::sourcePosition

protected

Definition at line 52 of file MedicalBeam.hh.

Referenced by GetSourcePosition(), and SetSourcePosition().

SSD

G4double MedicalBeam::SSD

protected

Definition at line 54 of file MedicalBeam.hh.

Referenced by GetSSD(), and SetSSD().

---

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

• environments/g4py/tests/g4pytest/MedicalBeam/MedicalBeam.hh
• environments/g4py/tests/g4pytest/MedicalBeam/MedicalBeam.cc