Message: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD 

Forum: Processes Involving Optical Photons
Re: Question G4Exception-START: ELECTRONSCINTILLATIONYIELD (Sergio Lo Meo)
Re: Feedback Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD (Gumplinger Peter)
Re: Question Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD (Sergio Lo Meo)
Date: 12 May, 2016
From: Gumplinger Peter <Gumplinger Peter>

Dear Sergio,

Okay, it escaped me yesterday, the scintillation yields must be specified as non-constant property with AddProperty (and not with AddConstProperty). The array is the scintillation light yield as function of incident particle energy!

To account for potential nonlinearity and scintillation photon density along the track, light (L) is produced according to:

L_currentStep = L(PreStepKE) - L(PreStepKE - EDep)

What you specify is the number of scintillation photons that an incident particle of energy E, produces in the scintillator if it were to stop in the scintillator!

> In my DetectorConstruction I set:
> 
>  G4double SCY_Energy[nentries] = {1.0*MeV, 20.0*MeV, 40.0*MeV, 60.0*MeV, 80.0*MeV, 100.0*MeV, 120.0*MeV, 140.0*MeV,
>                                   160.0*MeV,180.0*MeV,200.0*MeV,300.0*MeV,400.0*MeV, 500.0*MeV, 600.0*MeV,
>                                   700.0*MeV, 800.0*MeV, 900.0*MeV, 1000.0*MeV};

So, these are the incident particle kinetic energies.

> G4double SCYELECTRON[nentries]={500.0/MeV, 1750./MeV,3000./MeV,4200./MeV,6000./MeV, 7200./MeV, 9000./MeV, 10000./MeV,
>                                 10000./MeV, 10000./MeV,10000.0/MeV, 10000./MeV, 10000./MeV, 10000./MeV,10000.0/MeV,
>                                 10000./MeV, 10000./MeV, 10000./MeV,10000.0/MeV};

What is wrong is that contrary to SCINTILLATIONYIELD, where you give the number of photons per unit energy deposited, here you give the number of photons produced by the incident particle at that energy (e.g. don't divide by energy).

Perhaps, you can change your code and define multiple G4Scintillation processes, each attached to a different particle (type) and you specify for each process the SCINTILLATIONYIELD per deposited energy.

Best regards, Peter

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Question: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Sergio Lo Meo - 13 May, 2016)
(_ Feedback: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Gumplinger Peter - 13 May, 2016)
(_ Ok: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Sergio Lo Meo - 22 May, 2016)
(_ Question: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Nicholas DePorzio - 09 Aug, 2016)
(_ Feedback: Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD   (Gumplinger Peter - 11 Aug, 2016)
 Add Message Add Message
to: "Re: G4Exception-START: ELECTRONSCINTILLATIONYIELD"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews