|Message: Re: Does GEANT4 accurately model the trapping efficiency for a scintillating optical fiber?||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)|
thanks for your quick response.
> > Problem in brief: I am trying to model the trapping efficiency of a
> > scintillating optical fiber and the efficiency I measure is more than
> > 10x the manufacturer's value.
> You trap photons that are produced inside the fiber. The Cherenkov
> photons even have a specific cone distribution relative to the radiating
> particle's direction. To understand your situation it would be best to
> eliminate the Cherenkov photons (no such process) in your simulation for
> the moment.
I can certainly turn off Cherenkov. That could reduce the number of variables for now.
> > I am trying to model the scintillation and Cherenkov light for a Bicron
> > BCF-10 scintillating fiber. factory values for the refractive indices of
> > the fiber are n2 = 1.49 (cladding), and n1= 1.6 (fiber core).
> > Calculating the trapping efficiency using the equation eff =
> > 0.5(1-n2/n1) gives eff = 3.44%, which is the manufacturer stated value.
> Real fibers aren't perfectly smooth. A little roughness introduced
> changes the situation dramatically.
> > Strangely, when I calculate eff = (# photons absorbed in PMT)/(# of
> > scintillation photons + # of Cherenkov photons produced), i get eff =
> > 37%. (steppingAction used to keep track of produced and absorbed photons
> > is attached.)
> Looks ok.
> > Do I misunderstand something in the opticalPhysics process? Shouldn't
> > the number of photons trapped in the fiber and transported to the PMT be
> > equal to the value calculable based on the ratio of refractive indices
> > of the cladding and the fiber? I understand that value only describes
> > the meridial rays, but the skew rays only add a few percent.
> I would think that almost all of your rays are skew rays. Why do you
> think the opposite?
I agree the rays are mostly skew. I'm just confused as to why I'm detecting so many of them. Shouldn't the detection efficiency I'm calculating be closer to 1/2 the manufacturer stated trapping efficiency?