Message: Sciintillation Emission Spectrum table Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Sciintillation Emission Spectrum table 

Forum: Processes Involving Optical Photons
Date: 21 Jan, 2010
From: Ron Mathis <Ron Mathis>

Hello, I am trying to insert an Emission Spectrum table for the scintillator I am using. The resulting spectrum of the optical photons at the photocathode consists of disconnected sections of points, many of which have the wrong slope or a negative slope when it should be positive etc. The number of disconnected sections matches the number of entries in the table. If the number of entries is changed the number of disconnected sections changes accordingly. It is behaving as if the interpolation between points is not working correctly. The following is the code I am using.

  G4Material* EJ260 = new G4Material("EJ260",density= 1.023*g/cm3, nel=2);
  EJ260->AddElement(H, natoms=10);
  EJ260->AddElement(C, natoms= 9); 

 const G4int NUMej = 26;
  G4double ej_PhotonEnergy[NUMej] = {2.07*eV, 2.088*eV, 2.111*eV, 2.133*eV,
				     2.156*eV, 2.18*eV, 2.204*eV, 2.229*eV,
				     2.254*eV, 2.28*eV, 2.307*eV, 2.334*eV,
				     2.362*eV, 2.39*eV, 2.42*eV, 2.449*eV,
				     2.48*eV, 2.511*eV, 2.544*eV, 2.577*eV,
				     2.61*eV, 2.645*eV, 2.681*eV, 2.718*eV,
				     2.756*eV, 4.5*eV};
  G4double scint_FAST[NUMej] = {0., 0.01, 0.02, 0.03, 0.04, 0.05, 0.06,
				0.085, 0.12, 0.16, 0.22, 0.28, 0.32, 
				0.36, 0.44, 0.6, 0.83, 1.0, 0.9, 0.62, 
				0.36, 0.12, 0.01, 0., 0., 0.};
  G4double ej_RIND[NUMej] = {1.58, 1.58, 1.58, 1.58, 1.58, 1.58, 1.58, 
			     1.58, 1.58, 1.58, 1.58, 1.58, 1.58, 1.58,
			     1.58, 1.58, 1.58, 1.58, 1.58, 1.58, 1.58,
			     1.58, 1.58, 1.58, 1.58, 1.58};  // refractive index
  G4double ej_ABSL[NUMej] = {2.*m, 2.*m, 2.*m, 2.*m, 2.*m, 2.*m, 2.*m, 
			     2.*m, 2.*m, 2.*m, 2.*m, 2.*m, 2.*m, 2.*m,
			     2.*m, 2.*m, 2.*m, 2.*m, 2.*m, 2.*m, 2.*m,
			     2.*m, 2.*m, 2.*m, 2.*m, 2.*m};
			     // atten length (Absorption Length)
  G4MaterialPropertiesTable* EJ260_mt = new G4MaterialPropertiesTable();
  EJ260_mt->AddProperty(
	    "FASTCOMPONENT", ej_PhotonEnergy, scint_FAST, NUMej);
  EJ260_mt->AddProperty("RINDEX", ej_PhotonEnergy, ej_RIND, NUMej);
  EJ260_mt->AddProperty("ABSLENGTH", ej_PhotonEnergy, ej_ABSL, NUMej);

  EJ260_mt->AddConstProperty("SCINTILLATIONYIELD", 9200./MeV);
  EJ260_mt->AddConstProperty("RESOLUTIONSCALE", 2.0);//Broaden distribution
  EJ260_mt->AddConstProperty("FASTTIMECONSTANT", 9.2*ns);
  EJ260_mt->AddConstProperty("YIELDRATIO", 1.0); //No slow component
  EJ260->SetMaterialPropertiesTable(EJ260_mt);

Any suggestions? Thank you for your help Ron Mathis

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

1 Feedback: Re: Sciintillation Emission Spectrum table   (Benton Pahlka - 21 Jan, 2010)
2 Feedback: Re: Sciintillation Emission Spectrum table   (Gumplinger Peter - 21 Jan, 2010)
1 More: Re: Sciintillation Emission Spectrum table   (Ron Mathis - 22 Jan, 2010)
(_ Warning: Re: Sciintillation Emission Spectrum table   (Gumplinger Peter - 22 Jan, 2010)
(_ Feedback: Re: Sciintillation Emission Spectrum table   (Ron Mathis - 25 Jan, 2010)
 Add Message Add Message
to: "Sciintillation Emission Spectrum table"

 Subscribe Subscribe

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