118 if (!aMaterialPropertiesTable)
127 G4int NumPhotons = 1;
131 G4double MeanNumberOfPhotons = aMaterialPropertiesTable->
136 if (NumPhotons <= 0) {
160 WLSTime = aMaterialPropertiesTable->
169 G4int NumberOfPhotons = NumPhotons;
171 for (
G4int i = 0; i < NumPhotons; i++) {
177 for (
G4int j = 1; j <= 100; j++) {
182 sampledEnergy = WLSIntegral->
GetEnergy(CIIvalue);
185 G4cout <<
"sampledEnergy = " << sampledEnergy <<
G4endl;
189 if (sampledEnergy <= primaryEnergy)
break;
194 if (sampledEnergy > primaryEnergy) {
196 G4cout <<
" *** One less WLS photon will be returned ***" <<
G4endl;
199 if (NumberOfPhotons == 0) {
201 G4cout <<
" *** No WLS photon can be sampled for this primary ***"
212 G4double sint = std::sqrt((1.-cost)*(1.+cost));
237 sinp = std::sin(phi);
238 cosp = std::cos(phi);
240 photonPolarization = cosp * photonPolarization + sinp * perp;
242 photonPolarization = photonPolarization.
unit();
250 (photonPolarization.
x(),
251 photonPolarization.
y(),
252 photonPolarization.
z());
266 new G4Track(aWLSPhoton,aSecondaryTime,aSecondaryPosition);
277 G4cout <<
"\n Exiting from G4OpWLS::DoIt -- NumberOfSecondaries = "
305 for (
G4int i=0 ; i < numOfMaterials; i++)
313 G4Material* aMaterial = (*theMaterialTable)[i];
318 if (aMaterialPropertiesTable) {
328 G4double currentIN = (*theWLSVector)[0];
330 if (currentIN >= 0.0) {
338 aPhysicsOrderedFreeVector->
339 InsertValues(currentPM , currentCII);
354 currentPM = theWLSVector->
Energy(j);
355 currentIN = (*theWLSVector)[j];
357 currentCII = 0.5 * (prevIN + currentIN);
359 currentCII = prevCII +
360 (currentPM - prevPM) * currentCII;
362 aPhysicsOrderedFreeVector->
363 InsertValues(currentPM, currentCII);
366 prevCII = currentCII;
399 if ( aMaterialPropertyTable ) {
400 AttenuationLengthVector = aMaterialPropertyTable->
402 if ( AttenuationLengthVector ){
403 AttenuationLength = AttenuationLengthVector->
404 Value(thePhotonEnergy);
414 return AttenuationLength;
425 else if (name ==
"exponential")
435 "generator does not exist");
G4double Energy(size_t index) const
G4VWLSTimeGeneratorProfile * WLSTimeGeneratorProfile
G4VParticleChange * PostStepDoIt(const G4Track &aTrack, const G4Step &aStep)
void SetPolarization(G4double polX, G4double polY, G4double polZ)
G4int GetNumberOfSecondaries() const
static size_t GetNumberOfMaterials()
static G4MaterialTable * GetMaterialTable()
void SetKineticEnergy(G4double aEnergy)
G4OpWLS(const G4String &processName="OpWLS", G4ProcessType type=fOptical)
G4double GetEnergy(G4double aValue)
const G4TouchableHandle & GetTouchableHandle() const
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
void insertAt(size_t, G4PhysicsVector *)
void AddSecondary(G4Track *aSecondary)
G4ParticleChange aParticleChange
void BuildPhysicsTable(const G4ParticleDefinition &aParticleType)
void UseTimeProfile(const G4String name)
G4bool ConstPropertyExists(const char *key) const
virtual void Initialize(const G4Track &)
const G4ThreeVector & GetPosition() const
static constexpr double twopi
static G4OpticalPhoton * OpticalPhoton()
void SetTouchableHandle(const G4TouchableHandle &apValue)
const G4String & GetProcessName() const
G4double GetGlobalTime() const
G4StepPoint * GetPostStepPoint() const
G4Material * GetMaterial() const
Hep3Vector cross(const Hep3Vector &) const
std::vector< G4Material * > G4MaterialTable
G4PhysicsTable * theIntegralTable
G4double GetMeanFreePath(const G4Track &aTrack, G4double, G4ForceCondition *)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
virtual G4double GenerateTime(const G4double time_constant)=0
void SetParentID(const G4int aValue)
void SetProcessSubType(G4int)
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
G4MaterialPropertyVector * GetProperty(const char *key, G4bool warning=false)
G4long G4Poisson(G4double mean)
G4double GetTotalEnergy() const
void ProposeTrackStatus(G4TrackStatus status)
const G4DynamicParticle * GetDynamicParticle() const
G4MaterialPropertiesTable * GetMaterialPropertiesTable() const
size_t GetVectorLength() const
void SetNumberOfSecondaries(G4int totSecondaries)