|Message: Re: Cerenkov photons||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
Can you verify that the position:
is actually on a boundary in your geometry definition.
> > EventID = 9 Step # 1 Particle: opticalphoton ParentID: 2 > Volume: CoreParam Copy Number: 0 Process: OpRayleigh Creator Process:Cerenkov > Position: (-14.5506,0.114811,-0.135845) Energy: 4.3608
This tell us that the process that defined the step is not 'Transportation' but OpRayleigh which would with ~100% probablility produce a step length that doens't exactly end on a volume boundary.
> After turning on the G4OpBoundaryProcess verbose (level 1), I get the > following output:
> Scattering Photon! > Old Momentum Direction: (0.813042,0.272653,-0.514416) > Old Polarization: (-0.376582,-0.427564,-0.821812) > New Polarization: (0.327999,0.454307,0.828264) > Polarization Change: (0.327999,0.454307,0.828264) > New Momentum Direction: (-0.879668,0.466512,0.0924714) > Momentum Change: (-0.879668,0.466512,0.0924714)
You also turned on verbosity for G4OpRayleigh because the above is written there.
> Photon at Boundary! > thePrePV: CoreParam > thePostPV: CoreParam > *** NotAtBoundary ***
This is expected. G4OpBoundaryProcess is called at the end of every step of a photon, but it only kicks in if that step is really at a boundary. It tells us it is not: e.g.
(pPostStepPoint->GetStepStatus() != fGeomBoundary)
and so it is not a surprise that both thePrePV and thePostPV are CoreParam.
What is weired is that G4OpRayleigh is being called over and over again (and so it keeps on sampling a scattering angle.
Something else is not right ... maybe unrelated to geometry.
> Currently, I have not defined/associated any optical surfaces
I have a suspision ... did you set a finite Rayleigh scattering length for PMMA?
|Inline Depth:||Outline Depth:||Add message:|