Geant4-11
Public Member Functions | Private Attributes
G4RTPrimaryGeneratorAction Class Reference

#include <G4RTPrimaryGeneratorAction.hh>

Inheritance diagram for G4RTPrimaryGeneratorAction:
G4VUserPrimaryGeneratorAction

Public Member Functions

 G4RTPrimaryGeneratorAction ()
 
virtual void GeneratePrimaries (G4Event *anEvent)
 
void SetUp ()
 
virtual ~G4RTPrimaryGeneratorAction ()
 

Private Attributes

G4bool distortionOn
 
G4ThreeVector eyeDirection
 
G4ThreeVector eyePosition
 
G4double headAngle
 
G4int nColumn
 
G4int nRow
 
G4ParticleDefinitionparticle_definition
 
G4double particle_energy
 
G4ThreeVector particle_polarization
 
G4double particle_time
 
G4VPhysicalVolumepWorld
 
G4double stepAngle
 
G4ThreeVector up
 
G4double viewSpan
 
G4double viewSpanX
 
G4double viewSpanY
 
EInside whereisit
 

Detailed Description

Definition at line 40 of file G4RTPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

◆ G4RTPrimaryGeneratorAction()

G4RTPrimaryGeneratorAction::G4RTPrimaryGeneratorAction ( )

Definition at line 41 of file G4RTPrimaryGeneratorAction.cc.

42{
46 particle_time = 0.0;
48
49 pWorld = 0;
51
52 nRow = 0;
53 nColumn = 0;
54
57 up = G4ThreeVector(0,1,0);
58 headAngle = 0.0;
59 viewSpan = 0.0;
60 stepAngle = 0.0;
61 viewSpanX = 0.0;
62 viewSpanY = 0.0;
63
64 distortionOn = false;
65}
CLHEP::Hep3Vector G4ThreeVector
G4ParticleDefinition * particle_definition
@ kInside
Definition: geomdefs.hh:70
static constexpr double GeV
static const G4LorentzVector zero(0., 0., 0., 0.)

References distortionOn, eyeDirection, eyePosition, CLHEP::GeV, headAngle, kInside, nColumn, nRow, particle_definition, particle_energy, particle_polarization, particle_time, pWorld, stepAngle, up, viewSpan, viewSpanX, viewSpanY, whereisit, and anonymous_namespace{G4CascadeDeexciteBase.cc}::zero.

◆ ~G4RTPrimaryGeneratorAction()

G4RTPrimaryGeneratorAction::~G4RTPrimaryGeneratorAction ( )
virtual

Definition at line 67 of file G4RTPrimaryGeneratorAction.cc.

68{;}

Member Function Documentation

◆ GeneratePrimaries()

void G4RTPrimaryGeneratorAction::GeneratePrimaries ( G4Event anEvent)
virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 70 of file G4RTPrimaryGeneratorAction.cc.

71{
72 // Note: We don't use G4ParticleGun here, as instantiating a G4ParticleGun
73 // object causes creation of UI commands and corresponding UI messenger
74 // that interfare with normal G4ParticleGun UI commands.
75
76 // evId = iRow * nColumn + iColumn
77 G4int evId = anEvent->GetEventID();
78 G4int iRow = evId / nColumn;
79 G4int iColumn = evId % nColumn;
80 G4double angleX = -(viewSpanX/2. - G4double(iColumn)*stepAngle);
81 G4double angleY = viewSpanY/2. - G4double(iRow)*stepAngle;
82 G4ThreeVector rayDirection;
83 if(distortionOn)
84 { rayDirection = G4ThreeVector(-std::tan(angleX)/std::cos(angleY),std::tan(angleY)/std::cos(angleX),1.0); }
85 else
86 { rayDirection = G4ThreeVector(-std::tan(angleX),std::tan(angleY),1.0); }
87 G4double cp = std::cos(eyeDirection.phi());
88 G4double sp = std::sqrt(1.-cp*cp);
89 G4double ct = std::cos(eyeDirection.theta());
90 G4double st = std::sqrt(1.-ct*ct);
91 G4double gam = std::atan2(ct*cp*up.x()+ct*sp*up.y()-st*up.z(), -sp*up.x()+cp*up.y());
92 rayDirection.rotateZ(-gam);
93 rayDirection.rotateZ(headAngle);
94 rayDirection.rotateUz(eyeDirection);
95
96 G4ThreeVector rayPosition(eyePosition);
97 if (whereisit != kInside) {
98 // Eye position is outside the world, so move it inside.
99 G4double outsideDistance = pWorld->GetLogicalVolume()->GetSolid()->
100 DistanceToIn(rayPosition,rayDirection);
101 if(outsideDistance != kInfinity)
102 { rayPosition = rayPosition + (outsideDistance+0.001)*rayDirection; }
103 else
104 {
105 // Ray does not intercept world at all.
106 // Return without primary particle.
107 return;
108 }
109 }
110
111 // create a new vertex
112 G4PrimaryVertex* vertex = new G4PrimaryVertex(rayPosition,particle_time);
113
114 // create new primaries and set them to the vertex
118 particle->SetMass( mass );
119 particle->SetMomentumDirection( rayDirection.unit() );
123 vertex->SetPrimary( particle );
124
125 anEvent->AddPrimaryVertex( vertex );
126}
static const G4double cp
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
double z() const
Hep3Vector unit() const
double phi() const
double theta() const
double x() const
Hep3Vector & rotateZ(double)
Definition: ThreeVector.cc:107
double y() const
Hep3Vector & rotateUz(const Hep3Vector &)
Definition: ThreeVector.cc:33
G4int GetEventID() const
Definition: G4Event.hh:118
void AddPrimaryVertex(G4PrimaryVertex *aPrimaryVertex)
Definition: G4Event.hh:121
G4VSolid * GetSolid() const
void SetPolarization(const G4ThreeVector &pol)
void SetKineticEnergy(G4double eKin)
void SetMomentumDirection(const G4ThreeVector &p)
void SetMass(G4double mas)
void SetPrimary(G4PrimaryParticle *pp)
G4LogicalVolume * GetLogicalVolume() const
static const G4double kInfinity
Definition: geomdefs.hh:41

References G4Event::AddPrimaryVertex(), cp, distortionOn, eyeDirection, eyePosition, G4InuclParticleNames::gam, G4Event::GetEventID(), G4VPhysicalVolume::GetLogicalVolume(), G4ParticleDefinition::GetPDGMass(), G4LogicalVolume::GetSolid(), headAngle, kInfinity, kInside, nColumn, particle_definition, particle_energy, particle_polarization, particle_time, CLHEP::Hep3Vector::phi(), pWorld, CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::rotateZ(), G4PrimaryParticle::SetKineticEnergy(), G4PrimaryParticle::SetMass(), G4PrimaryParticle::SetMomentumDirection(), G4PrimaryParticle::SetPolarization(), G4PrimaryVertex::SetPrimary(), G4InuclParticleNames::sp, stepAngle, CLHEP::Hep3Vector::theta(), CLHEP::Hep3Vector::unit(), up, viewSpanX, viewSpanY, whereisit, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

◆ SetUp()

void G4RTPrimaryGeneratorAction::SetUp ( )

Definition at line 128 of file G4RTPrimaryGeneratorAction.cc.

129{
131 particle_definition = particleTable->FindParticle("geantino");
133 {
134 G4String msg;
135 msg = " G4RayTracer uses geantino to trace the ray, but your physics list does not\n";
136 msg += "define G4Geantino. Please add G4Geantino in your physics list.";
137 G4Exception("G4RTPrimaryGeneratorAction::SetUp","VisRayTracer00101",FatalException,msg);
138 }
139
141 nRow = rt->nRow;
142 nColumn = rt->nColumn;
145 viewSpan = rt->viewSpan;
146 stepAngle = viewSpan/100.;
150
152 GetNavigatorForTracking()->GetWorldVolume();
154}
@ FatalException
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
static G4ParticleTable * GetParticleTable()
static G4TheMTRayTracer * theInstance
G4ThreeVector eyeDirection
G4ThreeVector eyePosition
static G4TransportationManager * GetTransportationManager()
virtual EInside Inside(const G4ThreeVector &p) const =0

References distortionOn, G4TheRayTracer::distortionOn, eyeDirection, G4TheRayTracer::eyeDirection, eyePosition, G4TheRayTracer::eyePosition, FatalException, G4ParticleTable::FindParticle(), G4Exception(), G4VPhysicalVolume::GetLogicalVolume(), G4ParticleTable::GetParticleTable(), G4LogicalVolume::GetSolid(), G4TransportationManager::GetTransportationManager(), G4VSolid::Inside(), nColumn, G4TheRayTracer::nColumn, nRow, G4TheRayTracer::nRow, particle_definition, pWorld, stepAngle, G4TheMTRayTracer::theInstance, viewSpan, G4TheRayTracer::viewSpan, viewSpanX, viewSpanY, and whereisit.

Referenced by G4RTWorkerInitialization::WorkerRunStart().

Field Documentation

◆ distortionOn

G4bool G4RTPrimaryGeneratorAction::distortionOn
private

◆ eyeDirection

G4ThreeVector G4RTPrimaryGeneratorAction::eyeDirection
private

◆ eyePosition

G4ThreeVector G4RTPrimaryGeneratorAction::eyePosition
private

◆ headAngle

G4double G4RTPrimaryGeneratorAction::headAngle
private

Definition at line 65 of file G4RTPrimaryGeneratorAction.hh.

Referenced by G4RTPrimaryGeneratorAction(), and GeneratePrimaries().

◆ nColumn

G4int G4RTPrimaryGeneratorAction::nColumn
private

◆ nRow

G4int G4RTPrimaryGeneratorAction::nRow
private

Definition at line 60 of file G4RTPrimaryGeneratorAction.hh.

Referenced by G4RTPrimaryGeneratorAction(), and SetUp().

◆ particle_definition

G4ParticleDefinition* G4RTPrimaryGeneratorAction::particle_definition
private

◆ particle_energy

G4double G4RTPrimaryGeneratorAction::particle_energy
private

Definition at line 52 of file G4RTPrimaryGeneratorAction.hh.

Referenced by G4RTPrimaryGeneratorAction(), and GeneratePrimaries().

◆ particle_polarization

G4ThreeVector G4RTPrimaryGeneratorAction::particle_polarization
private

Definition at line 54 of file G4RTPrimaryGeneratorAction.hh.

Referenced by G4RTPrimaryGeneratorAction(), and GeneratePrimaries().

◆ particle_time

G4double G4RTPrimaryGeneratorAction::particle_time
private

Definition at line 53 of file G4RTPrimaryGeneratorAction.hh.

Referenced by G4RTPrimaryGeneratorAction(), and GeneratePrimaries().

◆ pWorld

G4VPhysicalVolume* G4RTPrimaryGeneratorAction::pWorld
private

◆ stepAngle

G4double G4RTPrimaryGeneratorAction::stepAngle
private

◆ up

G4ThreeVector G4RTPrimaryGeneratorAction::up
private

Definition at line 64 of file G4RTPrimaryGeneratorAction.hh.

Referenced by G4RTPrimaryGeneratorAction(), and GeneratePrimaries().

◆ viewSpan

G4double G4RTPrimaryGeneratorAction::viewSpan
private

Definition at line 66 of file G4RTPrimaryGeneratorAction.hh.

Referenced by G4RTPrimaryGeneratorAction(), and SetUp().

◆ viewSpanX

G4double G4RTPrimaryGeneratorAction::viewSpanX
private

◆ viewSpanY

G4double G4RTPrimaryGeneratorAction::viewSpanY
private

◆ whereisit

EInside G4RTPrimaryGeneratorAction::whereisit
private

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