G4RTPrimaryGeneratorAction Class Reference

#include <G4RTPrimaryGeneratorAction.hh>

Inheritance diagram for G4RTPrimaryGeneratorAction:

Collaboration diagram for G4RTPrimaryGeneratorAction:

Public Member Functions
	G4RTPrimaryGeneratorAction ()
virtual	~G4RTPrimaryGeneratorAction ()
virtual void	GeneratePrimaries (G4Event *anEvent)
void	SetUp ()
Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()
virtual	~G4VUserPrimaryGeneratorAction ()

Private Attributes
G4ParticleDefinition *	particle_definition
G4double	particle_energy
G4double	particle_time
G4ThreeVector	particle_polarization
G4VPhysicalVolume *	pWorld
EInside	whereisit
G4int	nColumn
G4int	nRow
G4ThreeVector	eyePosition
G4ThreeVector	eyeDirection
G4ThreeVector	up
G4double	headAngle
G4double	viewSpan
G4double	stepAngle
G4double	viewSpanX
G4double	viewSpanY
G4bool	distortionOn

Detailed Description

Definition at line 41 of file G4RTPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

G4RTPrimaryGeneratorAction()

G4RTPrimaryGeneratorAction::G4RTPrimaryGeneratorAction

(

)

Definition at line 42 of file G4RTPrimaryGeneratorAction.cc.

{
  G4ThreeVector zero;
  particle_definition = 0;
  particle_energy = 1.0*CLHEP::GeV;
  particle_time = 0.0;
  particle_polarization = zero;

  pWorld = 0;
  whereisit = kInside;

  nRow = 0;
  nColumn = 0;
  
  eyePosition = zero;
  eyeDirection = zero;
  up = G4ThreeVector(0,1,0);
  headAngle = 0.0;
  viewSpan = 0.0;
  stepAngle = 0.0;
  viewSpanX = 0.0;
  viewSpanY = 0.0;

  distortionOn = false;
}

~G4RTPrimaryGeneratorAction()

G4RTPrimaryGeneratorAction::~G4RTPrimaryGeneratorAction ( )

virtual

Definition at line 68 of file G4RTPrimaryGeneratorAction.cc.

{;}

Member Function Documentation

GeneratePrimaries()

void G4RTPrimaryGeneratorAction::GeneratePrimaries ( G4Event *  anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 71 of file G4RTPrimaryGeneratorAction.cc.

{
  // Note: We don't use G4ParticleGun here, as instantiating a G4ParticleGun
  //  object causes creation of UI commands and corresponding UI messenger
  //  that interfare with normal G4ParticleGun UI commands.

  // evId = iRow * nColumn + iColumn
  G4int evId = anEvent->GetEventID(); 
  G4int iRow = evId / nColumn;
  G4int iColumn = evId % nColumn;
  G4double angleX = -(viewSpanX/2. - G4double(iColumn)*stepAngle);
  G4double angleY = viewSpanY/2. - G4double(iRow)*stepAngle;
  G4ThreeVector rayDirection;
  if(distortionOn)
  { rayDirection = G4ThreeVector(-std::tan(angleX)/std::cos(angleY),std::tan(angleY)/std::cos(angleX),1.0); }
  else
  { rayDirection = G4ThreeVector(-std::tan(angleX),std::tan(angleY),1.0); }
  G4double cp = std::cos(eyeDirection.phi());
  G4double sp = std::sqrt(1.-cp*cp);
  G4double ct = std::cos(eyeDirection.theta());
  G4double st = std::sqrt(1.-ct*ct);
  G4double gam = std::atan2(ct*cp*up.x()+ct*sp*up.y()-st*up.z(), -sp*up.x()+cp*up.y());
  rayDirection.rotateZ(-gam);
  rayDirection.rotateZ(headAngle);
  rayDirection.rotateUz(eyeDirection);

  G4ThreeVector rayPosition(eyePosition);
  if (whereisit != kInside) {
    // Eye position is outside the world, so move it inside.
    G4double outsideDistance = pWorld->GetLogicalVolume()->GetSolid()->
       DistanceToIn(rayPosition,rayDirection);
    if(outsideDistance != kInfinity)
    { rayPosition = rayPosition + (outsideDistance+0.001)*rayDirection; }
    else
    {
      // Ray does not intercept world at all.
      // Return without primary particle.
      return;
    }
  }
  
  // create a new vertex
  G4PrimaryVertex* vertex = new G4PrimaryVertex(rayPosition,particle_time);

  // create new primaries and set them to the vertex
  G4double mass = particle_definition->GetPDGMass();
  G4PrimaryParticle* particle = new G4PrimaryParticle(particle_definition);
  particle->SetKineticEnergy( particle_energy );
  particle->SetMass( mass );
  particle->SetMomentumDirection( rayDirection.unit() );
  particle->SetPolarization(particle_polarization.x(),
                            particle_polarization.y(),
                            particle_polarization.z());
  vertex->SetPrimary( particle );

  anEvent->AddPrimaryVertex( vertex );
}

Here is the call graph for this function:

SetUp()

void G4RTPrimaryGeneratorAction::SetUp

(

)

Definition at line 129 of file G4RTPrimaryGeneratorAction.cc.

{
  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  particle_definition = particleTable->FindParticle("geantino");
  if(!particle_definition)
  {
    G4String msg;
    msg =  " G4RayTracer uses geantino to trace the ray, but your physics list does not\n";
    msg += "define G4Geantino. Please add G4Geantino in your physics list.";
    G4Exception("G4RTPrimaryGeneratorAction::SetUp","VisRayTracer00101",FatalException,msg);
  }

  G4TheMTRayTracer* rt = G4TheMTRayTracer::theInstance;
  nRow = rt->nRow;
  nColumn = rt->nColumn;
  eyePosition = rt->eyePosition;
  eyeDirection = rt->eyeDirection;
  viewSpan = rt->viewSpan;
  stepAngle = viewSpan/100.;
  viewSpanX = stepAngle*nColumn;
  viewSpanY = stepAngle*nRow;
  distortionOn = rt->distortionOn;

  pWorld = G4TransportationManager::GetTransportationManager()->
        GetNavigatorForTracking()->GetWorldVolume();
  whereisit = pWorld->GetLogicalVolume()->GetSolid()->Inside(eyePosition);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

distortionOn

G4bool G4RTPrimaryGeneratorAction::distortionOn

private

Definition at line 72 of file G4RTPrimaryGeneratorAction.hh.

eyeDirection

G4ThreeVector G4RTPrimaryGeneratorAction::eyeDirection

private

Definition at line 64 of file G4RTPrimaryGeneratorAction.hh.

eyePosition

G4ThreeVector G4RTPrimaryGeneratorAction::eyePosition

private

Definition at line 63 of file G4RTPrimaryGeneratorAction.hh.

headAngle

G4double G4RTPrimaryGeneratorAction::headAngle

private

Definition at line 66 of file G4RTPrimaryGeneratorAction.hh.

nColumn

G4int G4RTPrimaryGeneratorAction::nColumn

private

Definition at line 60 of file G4RTPrimaryGeneratorAction.hh.

nRow

G4int G4RTPrimaryGeneratorAction::nRow

private

Definition at line 61 of file G4RTPrimaryGeneratorAction.hh.

particle_definition

G4ParticleDefinition* G4RTPrimaryGeneratorAction::particle_definition

private

Definition at line 52 of file G4RTPrimaryGeneratorAction.hh.

particle_energy

G4double G4RTPrimaryGeneratorAction::particle_energy

private

Definition at line 53 of file G4RTPrimaryGeneratorAction.hh.

particle_polarization

G4ThreeVector G4RTPrimaryGeneratorAction::particle_polarization

private

Definition at line 55 of file G4RTPrimaryGeneratorAction.hh.

particle_time

G4double G4RTPrimaryGeneratorAction::particle_time

private

Definition at line 54 of file G4RTPrimaryGeneratorAction.hh.

pWorld

G4VPhysicalVolume* G4RTPrimaryGeneratorAction::pWorld

private

Definition at line 57 of file G4RTPrimaryGeneratorAction.hh.

stepAngle

G4double G4RTPrimaryGeneratorAction::stepAngle

private

Definition at line 68 of file G4RTPrimaryGeneratorAction.hh.

up

G4ThreeVector G4RTPrimaryGeneratorAction::up

private

Definition at line 65 of file G4RTPrimaryGeneratorAction.hh.

viewSpan

G4double G4RTPrimaryGeneratorAction::viewSpan

private

Definition at line 67 of file G4RTPrimaryGeneratorAction.hh.

viewSpanX

G4double G4RTPrimaryGeneratorAction::viewSpanX

private

Definition at line 69 of file G4RTPrimaryGeneratorAction.hh.

viewSpanY

G4double G4RTPrimaryGeneratorAction::viewSpanY

private

Definition at line 70 of file G4RTPrimaryGeneratorAction.hh.

whereisit

EInside G4RTPrimaryGeneratorAction::whereisit

private

Definition at line 58 of file G4RTPrimaryGeneratorAction.hh.

---

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

• G4RTPrimaryGeneratorAction.hh
• G4RTPrimaryGeneratorAction.cc