Hi,
You first have to enable particle type-dependent scintillation by:
G4OpticalParameters::Instance()->SetScintByParticleType(true);
You need a scintillation yield for each particle of interest (example following is for e- (due to gamma ray interaction) and p (neutron recoil) for NE213 - you can also do alphas, tritons and ions, if you have the spec sheet or literature data):
MPTNE213->AddProperty(âELECTRONSCINTILLATIONYIELDâ, particleEnergy, scintYieldElectron);
MPTNE213->AddProperty(âPROTONSCINTILLATIONYIELDâ, particleEnergy, scintYieldProton);
In Geant4 you can have up to three decay constants, which are set by:
MPTNE213->AddConstProperty("SCINTILLATIONTIMECONSTANT1",3.16*ns); // Eljen EJ-301 data sheet and St Gobain for BC-501A - fast component
MPTNE213->AddConstProperty("SCINTILLATIONTIMECONSTANT2",32.3*ns); // slow component
MPTNE213->AddConstProperty("SCINTILLATIONTIMECONSTANT3",270.0*ns); // Eljen, St. Gobain both quote an additional very slow component at 270ns
The decay constant values are particle-independent, but their amplitude ratios vary by particle. The ratios are set by:
MPTNE213->AddConstProperty(âELECTRONSCINTILLATIONYIELD1â,0.87); //Relative yield of component 1 for electrons
MPTNE213->AddConstProperty(âELECTRONSCINTILLATIONYIELD2â,0.13); //Relative yield of component 2
MPTNE213->AddConstProperty(âELECTRONSCINTILLATIONYIELD3â,0.001); //Relative yield of component 3
MPTNE213->AddConstProperty(âPROTONSCINTILLATIONYIELD1â,0.80); //Relative yield of component 1 for protons (estimate from literature)
MPTNE213->AddConstProperty(âPROTONSCINTILLATIONYIELD2â,0.20); //Relative yield of component 2
MPTNE213->AddConstProperty(âPROTONSCINTILLATIONYIELD3â,0.001); //Relative yield of component 3
The Book For Application Developers has a very good description of what all needs to be done Physics Processes â Book For Application Developers 11.0 documentation.