118 if (!aMaterialPropertiesTable)
122 aMaterialPropertiesTable->
GetProperty(
"WLSCOMPONENT");
127 G4int NumPhotons = 1;
131 G4double MeanNumberOfPhotons = aMaterialPropertiesTable->
132 GetConstProperty(
"WLSMEANNUMBERPHOTONS");
136 if (NumPhotons <= 0) {
160 WLSTime = aMaterialPropertiesTable->
161 GetConstProperty(
"WLSTIMECONSTANT");
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));
234 G4ThreeVector perp = photonMomentum.cross(photonPolarization);
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) {
321 aMaterialPropertiesTable->
GetProperty(
"WLSCOMPONENT");
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->
401 GetProperty(
"WLSABSLENGTH");
402 if ( AttenuationLengthVector ){
403 AttenuationLength = AttenuationLengthVector->
404 Value(thePhotonEnergy);
414 return AttenuationLength;
425 else if (name ==
"exponential")
435 "generator does not exist");
G4long G4Poisson(G4double mean)
G4int GetNumberOfSecondaries() const
G4MaterialPropertyVector * GetProperty(const char *key)
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
G4double GetTotalEnergy() const
const G4DynamicParticle * GetDynamicParticle() const
G4VParticleChange * PostStepDoIt(const G4Track &aTrack, const G4Step &aStep)
static G4MaterialTable * GetMaterialTable()
std::vector< G4Material * > G4MaterialTable
void SetTouchableHandle(const G4TouchableHandle &apValue)
size_t GetVectorLength() const
void UseTimeProfile(const G4String name)
G4GLOB_DLL std::ostream G4cout
virtual G4double GenerateTime(const G4double time_constant)=0
const G4ThreeVector & GetPosition() const
void BuildPhysicsTable(const G4ParticleDefinition &aParticleType)
G4VWLSTimeGeneratorProfile * WLSTimeGeneratorProfile
static const double twopi
void SetProcessSubType(G4int)
void SetPolarization(G4double polX, G4double polY, G4double polZ)
G4double Energy(size_t index) const
static size_t GetNumberOfMaterials()
const G4String & GetProcessName() const
void SetKineticEnergy(G4double aEnergy)
const G4TouchableHandle & GetTouchableHandle() const
G4double GetEnergy(G4double aValue)
G4Material * GetMaterial() const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
static G4OpticalPhoton * OpticalPhoton()
virtual void Initialize(const G4Track &)
G4MaterialPropertiesTable * GetMaterialPropertiesTable() const
void SetNumberOfSecondaries(G4int totSecondaries)
void SetParentID(const G4int aValue)
G4StepPoint * GetPostStepPoint() const
G4ParticleChange aParticleChange
G4bool ConstPropertyExists(const char *key)
G4OpWLS(const G4String &processName="OpWLS", G4ProcessType type=fOptical)
void AddSecondary(G4Track *aSecondary)
void insertAt(size_t, G4PhysicsVector *)
G4double GetGlobalTime() const
void ProposeTrackStatus(G4TrackStatus status)
G4PhysicsTable * theIntegralTable
G4ThreeVector G4ParticleMomentum
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
G4double GetMeanFreePath(const G4Track &aTrack, G4double, G4ForceCondition *)