Dear Paola,
I double checked my results, and indeed there is a bug in my analysis
program. I was double counting the energy loss...
Thank you very much for your time!
Shirley
On Fri, Jun 8, 2018 at 8:27 AM, Paola Sala <paola.sala_at_mi.infn.it> wrote:
> Dear Shirley
> I changed your mgdraw output to formatted, just to read it easily, and I
> do not see the problem. Dumping the quantities
> NCASE, JTRACK DTQUEN(1) DTRACK(1), where dtquen is quenched and
> dtrack is unquenched, I always obtain dtquen<dtrack as expected:
> 1 3 5.26920596E-09 5.31032152E-09
> 1 3 1.79176452E-03 2.06275354E-03
> 1 3 2.61539157E-04 2.96100159E-04
> and so on.
> When I sum all energies for each event, I get quenched around 0.7 GeV,
> unquenched around 0.8 GeV.
> Maybe you would like to send also your analysis program?
> Paola
>
> > Dear FLUKA users,
> >
> > I'm trying to enable quenching in materials and record continuous energy
> > loss in userdump MGDRAW routine.
> >
> > I inject 1 GeV electron in liquid argon. Without quenching, if I count
> > all
> > continuous energy loss through DTRACK, I get ~ 0.9 GeV energy, which
> seems
> > reasonable. But when I include quenching through UDQUENCH, somehow the
> > total continuous energy loss through either the un-quenched(?) DTRACK or
> > DTQUEN for every primary event is ~ 1.5 GeV > 1 GeV.
> >
> > Could anyone see where I made mistakes? I attached my input files.
> >
> > Thank you very much!
> >
> > Shirley
> >
>
>
> Paola Sala
> INFN Milano
> tel. Milano +39-0250317374
> tel. CERN +41-227679148
>
>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Sat Jun 09 2018 - 01:20:27 CEST