15 #include "CLHEP/Random/RandGaussQ.h"
17 #include "G4ParticleDefinition.hh"
18 #include "G4ParticleGun.hh"
22 G4int nPhotonsGun = 1;
23 G4int nPhotonsBeam = 1;
24 G4double Energy = 1.15 * eV;
54 G4ParticleTable * theParticleTable = G4ParticleTable::GetParticleTable();
55 G4ParticleDefinition * theOpticalPhoton = theParticleTable->FindParticle(
"opticalphoton");
59 theParticleGun->SetParticlePosition(G4ThreeVector(-500.0 * cm, 0.0 * cm, 0.0 * cm));
60 theParticleGun->SetParticleMomentumDirection(G4ThreeVector(5.0, 3.0, 0.0));
80 unsigned int seed = rng->
mySeed();
87 const G4int nLaserRings = 2;
88 const G4int nLaserBeams = 8;
91 G4double LaserPositionZ = -2057.5 * mm;
94 G4double LaserRingRadius[nLaserRings] = {564.0 * mm, 840.0 * mm};
97 G4double LaserPhi0 = 0.392699082;
100 G4double LaserRingSigmaX[nLaserRings] = {0.5 * mm, 0.5 * mm};
101 G4double LaserRingSigmaY[nLaserRings] = {0.5 * mm, 0.5 * mm};
104 G4ParticleTable * theParticleTable = G4ParticleTable::GetParticleTable();
105 G4ParticleDefinition * theOpticalPhoton= theParticleTable->FindParticle(
"opticalphoton");
108 for (
int theRing = 0; theRing < nLaserRings; theRing++)
111 for (
int theBeam = 0; theBeam < nLaserBeams; theBeam++)
115 G4double LaserPositionPhi = LaserPhi0 + G4double(theBeam * G4double(G4double(2 *
M_PI) / nLaserBeams));
118 G4double LaserPositionX =
cos(LaserPositionPhi) * LaserRingRadius[theRing];
119 G4double LaserPositionY =
sin(LaserPositionPhi) * LaserRingRadius[theRing];
122 for (
int theParticle = 0; theParticle <
thenParticle; theParticle++)
125 CLHEP::RandGaussQ aGaussObjX( *
theDRand48Engine, LaserPositionX, LaserRingSigmaX[theRing] );
126 CLHEP::RandGaussQ aGaussObjY( *
theDRand48Engine, LaserPositionY, LaserRingSigmaY[theRing] );
128 G4double theXPosition = aGaussObjX.fire();
129 G4double theYPosition = aGaussObjY.fire();
130 G4double theZPosition = LaserPositionZ;
135 theParticleGun->SetParticlePosition(G4ThreeVector(theXPosition, theYPosition, theZPosition));
139 for (
int theDirection = 0; theDirection < 2; theDirection++)
142 if (theDirection == 0)
144 theParticleGun->SetParticleMomentumDirection(G4ThreeVector(0.0, 0.0, 1.0));
151 if (theDirection == 1)
153 theParticleGun->SetParticleMomentumDirection(G4ThreeVector(0.0, 0.0, -1.0));
173 if (
theParticleGun->GetParticleDefinition()->GetParticleName() !=
"opticalphoton" )
175 edm::LogWarning(
"SimLaserAlignment:LaserBeamsTEC2") <<
"<LaserBeamsTEC2::setOptPhotonPolar()>: WARNING! The ParticleGun is not an optical photon";
180 G4ThreeVector normal(1.0, 0.0, 0.0);
181 G4ThreeVector kphoton =
theParticleGun->GetParticleMomentumDirection();
182 G4ThreeVector product = normal.cross(kphoton);
183 G4double modul2 = product * product;
185 G4ThreeVector e_perpendicular(0.0, 0.0, 1.0);
187 if ( modul2 > 0.0 ) { e_perpendicular = (1.0 /
sqrt(modul2)) * product; }
189 G4ThreeVector e_parallel = e_perpendicular.cross(kphoton);
191 G4ThreeVector polar =
cos(Angle) * e_parallel +
sin(Angle) * e_perpendicular;
void setOptPhotonPolar(G4double Angle)
set the polarisation of the photons
G4ParticleGun * theParticleGun
Sin< T >::type sin(const T &t)
LaserBeamsTEC2()
default constructor
~LaserBeamsTEC2()
destructor
void GeneratePrimaries(G4Event *myEvent)
shoot optical photons into the detector at the beginning of an event
CLHEP::DRand48Engine * theDRand48Engine
Cos< T >::type cos(const T &t)
virtual uint32_t mySeed() const =0
Exists for backward compatibility.