Dear Alfredo,
Thanks for your quick reply.
Q: could you please explain a bit better what do you mean by "%"?
Is it the fraction of K+ at the end of your setup for an incoming beam of
K0, K0L, K0s or what?
A: The yield rate(%) is defined by
N_K+
K+ Yield(%) = -------- × 100(%)
N_Tot
N_K+ : event numbers of K+ produced by the inelastic process;
N_Tot: total event numbers simulated.
Eg : I simulate total 10000 events of K0, I find that there are
3516 events will produce K+ by inelastic interaction during K0
propagate through the high density material. Then,
K+ yield rate = 3516/10000 * 100%= 35.16%
How to count the N_K+, the code is in the ENTRY USDRAW of mgdraw.f
as attachment. Here I also show them below:
* ------------------------ Kaon+ count
NKAONP = 0
OPEN(UNIT = 80, FILE = "usdraw.dat", ACCESS = 'APPEND')
DO 1000 ip = 1, NP
IF (ICODE.EQ.101 .AND. KPART(ip).EQ.15) THEN
NKAONP = NKAONP + 1
WRITE (80,*) NCASE, NKAONP
END IF
1000 CONTINUE
*-----------------------------------------
PS: I also print out the track infos trying to understand the difference,
1. In fluka, the interaction process like followed:
K0 : K0 will scattering off the nucleus by elastic or inelastic, or decay.
No K0L or K0S shown up.
K0L: One condition is that K0L will be transformed into K0(K0bar) by elastic
process, then K0(K0bar) will go through by elastic, inelastic and decay.
The other condition is that K0L will directly scattering with the nucleus by
elastic or inelastic, or decay.
K0S: Just like K0L.
2. In Geant4, the interaction process like followed:
K0 : It will be directly transformed into K0L or K0S by 50% probability,
then K0L or K0S will transport through the materail by elastic or inelastic or decay.
No K0(K0bar) shown up.
K0L: It will transport the material by elastic, inelastic process and decay,
No K0(K0bar) shown up.
K0S: Just like K0L.
3. For Geant4, it is clear that the K+ yield rate from K0 is equal to the half
of the summation yield rate from K0L and K0S, just as shown in the attached JPG file.
For FLUKA, why the K+ yield rate is not simply equal to the half of the summation
yield rate from K0L and K0S ?
Looking forward to your suggestions.
Thanks very much.
Best wishes,
Hang
> -----原始邮件-----
> 发件人: alfredo.ferrari_at_mi.infn.it
> 发送时间: 2021-12-10 17:05:41 (星期五)
> 收件人: "\"胡航\"" <huhang_at_ihep.ac.cn>
> 抄送: "fluka-discuss_at_fluka.org" <fluka-discuss_at_fluka.org>
> 主题: Re: [fluka-discuss]: Fluka vs Geant4: KaonPlus production rate
>
> Dear Hang
>
> could you please explain a bit better what do you mean by "%"?
>
> Is it the fraction of K+ at the end of your setup for an incoming beam of
> K0, K0L, K0s or what?
>
> There could surely be differences between FLUKA and Geant4 wrt the
> interaction models, however in order to assess if this is the reason for
> the differences you see some extra infos are required.
>
> Thanks
> Alfredo
>
> > Dear experts,
> >
> >
> >
> >
> > I am trying to simulate KaonZero/KaonZeroLong/KaonZeroShort which
> > propagate through a high density material
> >
> > to get the yield rate of KaonPlus. The detailed setup as shown in the
> > input file as attachment.
> >
> > I use the USDRAW entry in mgdraw.f to count the KaonPlus produced by
> > inelastic process.
> >
> > Also I use the Geant4 to the same simualtion. The comparsion of the
> > KaonPlus production rate as the attached JPG file.
> >
> >
> >
> > Here is my question: Why there is difference in the KaonPlus yield
> > rate? Does due to the cross section or handling in the
> >
> > interaction process?
> >
> >
> >
> > Looking forward to your reply.
> >
> >
> >
> >
> >
> >
> > Best regards,
> >
> >
> >
> >
> > Hang
>
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Sat Dec 11 2021 - 04:34:56 CET