GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » GEANT3 dEdx for low energy protons
GEANT3 dEdx for low energy protons [message #6306] Tue, 08 April 2008 18:02 Go to next message
Sebastian Neubert is currently offline  Sebastian Neubert
Messages: 282
Registered: March 2006
Location: Munich
first-grade participant

From: *e18.physik.tu-muenchen.de
Hi!

We observe a remarkable behaviour with GEANT3 when we look at the energy loss distributions for 200MeV protons.

From BetheBloch we would expect something like dE/dx=22keV/cm in the TPC gas.

Please have a look to the following plots, where you see the dE/dx (in keV/cm) from MC-Points plotted against the dx (in cm) for different step limitations:

dx<2.5cm
index.php?t=getfile&id=4340&private=0

dx<3cm
index.php?t=getfile&id=4341&private=0

dx<4cm
index.php?t=getfile&id=4342&private=0

One observes that as soon as a step is smaller than 3cm the dEdx calculation gives a wrong result, which is too low by a factor of 3 to 4 (6keV/cm instead of 22keV/cm).

It is remarkable, that this behaviour is observed for any short step, regardless of if it was triggered by the step limiter or by anything else. You see this in the third plot, where the points with too low dEdx still persist!

If we shut off the straggeling and just use the energy loss tables of GEANT (LOSS=4 option in SetCuts.C) we get the exepcted value for the dE/dx.

How should we deal with this?

Cheers!
Viola and Sebastian.
  • Attachment: dEdx2.5.png
    (Size: 14.16KB, Downloaded 775 times)
  • Attachment: dEdx3.png
    (Size: 14.42KB, Downloaded 710 times)
  • Attachment: dEdx4.png
    (Size: 14.68KB, Downloaded 760 times)


Sebastian Neubert
Technische Universität München
Department Physik E18
sneubert@e18.physik.tu-muenchen.de
tel: +49-8928912592

[Updated on: Tue, 08 April 2008 18:03]

Report message to a moderator

Re: GEANT3 dEdx for low energy protons [message #6454 is a reply to message #6306] Wed, 16 April 2008 14:23 Go to previous messageGo to next message
Sebastian Neubert is currently offline  Sebastian Neubert
Messages: 282
Registered: March 2006
Location: Munich
first-grade participant

From: *e18.physik.tu-muenchen.de
Hi!

We have found out what was going wrong with the energy loss.

The option "STRA" was turned to 1 in g3config. As I read from the forum contributions from november Stefano suggested this, because this setting is used in Hades.

I have put some debug output into the geant3 code and observe that for STRA=1 only the G3STREN function gives a contribution to the energy loss with the behaviour that I have shown lately.

When you switch to the Urban model (STRA=0) it works.

What should we do?

Sebastian.


Sebastian Neubert
Technische Universität München
Department Physik E18
sneubert@e18.physik.tu-muenchen.de
tel: +49-8928912592
Re: GEANT3 dEdx for low energy protons [message #6577 is a reply to message #6454] Wed, 23 April 2008 10:00 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *physik.uni-giessen.de
Could you please try to change in SetCuts.C the variable cut1 from 1MeV to 1keV?
In the MVD case it helped a bit. I wonder if it works even for TPC.
Re: GEANT3 dEdx for low energy protons [message #6578 is a reply to message #6577] Wed, 23 April 2008 10:21 Go to previous message
Sebastian Neubert is currently offline  Sebastian Neubert
Messages: 282
Registered: March 2006
Location: Munich
first-grade participant

From: *e18.physik.tu-muenchen.de
Hi Stefano!

We have played around with the cut variable a lot. It does not change the picture. The only thing that helps is to use the Urban model which according to the GEANT3 manual is also the recommended thing to do.

Sebastian.


Sebastian Neubert
Technische Universität München
Department Physik E18
sneubert@e18.physik.tu-muenchen.de
tel: +49-8928912592
Previous Topic: libPGen pgenerators
Next Topic: Compilation problems with DpmEvtGen
Goto Forum:
  


Current Time: Mon Nov 25 15:37:05 CET 2024

Total time taken to generate the page: 0.01120 seconds