Message: Re: The confused of Cherenkov photons Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: The confused of Cherenkov photons 

Forum: Particles
Re: None The confused of Cherenkov photons (yanxin0221)
Re: Feedback Re: The confused of Cherenkov photons (Gumplinger Peter)
Re: None Re: The confused of Cherenkov photons (yanxin0221)
Re: Disagree Re: The confused of Cherenkov photons (Gumplinger Peter)
Re: Question Re: The confused of Cherenkov photons (yanxin0221)
Date: 16 Sep, 2013
From: Gumplinger Peter <Gumplinger Peter>

> I have solved the reflected Cherenkov photons according to your suggestion.But

> The new problem is that the light spectrum is wrong.Can you help me find errors

> in my code.

What spectrum is wrong? The Cherenkov spectrum or your scintillation spectrum?

You scintillation spectrum is a combination of a flat fast component:

  G4double ScintilFast[nEntries] =
            { 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
              1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00,
              1.00, 1.00, 1.00, 1.00 };

and this slow component:

  G4double ScintilSlow[nEntries] =
            { 0.01, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00,
              7.00, 8.00, 9.00, 8.00, 7.00, 6.00, 4.00,
              3.00, 2.00, 1.00, 0.01, 1.00, 2.00, 3.00,
              4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 8.00,
              7.00, 6.00, 5.00, 4.00 };

in the energy range of:

  G4double PhotonEnergy[nEntries] =
            { 2.034*eV, 2.068*eV, 2.103*eV, 2.139*eV,
              2.177*eV, 2.216*eV, 2.256*eV, 2.298*eV,
              2.341*eV, 2.386*eV, 2.433*eV, 2.481*eV,
              2.532*eV, 2.585*eV, 2.640*eV, 2.697*eV,
              2.757*eV, 2.820*eV, 2.885*eV, 2.954*eV,
              3.026*eV, 3.102*eV, 3.181*eV, 3.265*eV,
              3.353*eV, 3.446*eV, 3.545*eV, 3.649*eV,
              3.760*eV, 3.877*eV, 4.002*eV, 4.136*eV };

The exact combination you requested is:

  myMPT1->AddConstProperty("YIELDRATIO",0.8);

My guess is that you don't whan a flat slow scintillation spectrum at all. If you don't, then don't specify any and don't set the "YIELDRATIO". If you only want Cherenkov light, don't add the G4Scintillation process. And if you only want scintillation light, don't add the Cherenkov process in your PhysicsList.

With this my tutoring is exhausted.

 Add Message Add Message
to: "Re: The confused of Cherenkov photons"

 Subscribe Subscribe

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


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]