GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Monte Carlo Engines » NIEL
NIEL [message #22449] Wed, 11 April 2018 06:31 Go to next message
Sergey Kononov is currently offline  Sergey Kononov
Messages: 11
Registered: September 2014
occasional visitor

From: *inp.nsk.su
Hi!
How to extract NIEL (non-ionizing energy losses) in PandaRoot? I know it is possible in Geant4.
Best regards,
Sergey


Cheers,
Sergey
Re: NIEL [message #22468 is a reply to message #22449] Thu, 12 April 2018 13:29 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *gsi.de
Dear Sergey,

We have access to the Virtual Monte-Carlo (VMC) wrapper during simulation. The class documentaion is ROOT's TVirtualMC: https://root.cern.ch/doc/master/classTVirtualMC.html
Via google I found a pdf stating simewhere that TVirtualMC cannot deliver NIEL and a direct access to Geant4 is needed.

Kind regards from GSI
Ralf
Re: NIEL [message #22469 is a reply to message #22468] Thu, 12 April 2018 19:50 Go to previous messageGo to next message
Mohammad Al-Turany is currently offline  Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany
first-grade participant
From: *dip0.t-ipconnect.de
Hi,

I forward this issue to Ivana (the main VMC author) and here is the answer:


Hi Mohammad,

It is not so straightforward.
Users can extend Geant4 VMC user action classes, eg TG4SteppingAction and access directly Geant4 classes (G4Step) which provide this information. However this mode is not in regular testing so getting this work may require some debugging.

The simplest way how to provide this information would be to add
a dedicated function for NIEL Edep in TVirtualMC interface;
the user would need to wait when this update (Root && VMC packages) is
available in tagged packages.

I have progressed with sensitive detectors, which will also require modifications in Root vmc core library, so this addition may be added with these changes. Let me know.

Please, pass my reply to your user, or to your forum, as you wish.

Cheers,

Ivana

Re: NIEL [message #22473 is a reply to message #22469] Fri, 13 April 2018 03:52 Go to previous messageGo to next message
Sergey Kononov is currently offline  Sergey Kononov
Messages: 11
Registered: September 2014
occasional visitor

From: *92.127.232.135.nsk.sibirtelecom.ru
Dear Mohammad,
Thank you very much for the answer and triggering a request to the VMC developer.
We'll wait when it will be available. We have some workaround for calculating 1 MeV neutron equvalent flux for particles but it would be good to cross-check with direct simulation.
Best regards,
Sergey


Cheers,
Sergey
Re: NIEL [message #22501 is a reply to message #22449] Tue, 17 April 2018 16:55 Go to previous messageGo to next message
Tobias Stockmanns is currently offline  Tobias Stockmanns
Messages: 489
Registered: May 2007
first-grade participant
From: *ikp.kfa-juelich.de
Dear Sergey,

sorry for my late reply.
In the MVD we use pre-stored lists which contain the conversion values to go from particle momentum and type to the corresponding neutron flux.

Have a look at detectors/mvd/MvdTools. In there are the stored lists /MvdRadDamage and the PndMvdRadDamTask which does the conversion.

I hope this helps.

Cheers,

Tobias

Re: NIEL [message #22506 is a reply to message #22501] Wed, 18 April 2018 09:27 Go to previous message
Sergey Kononov is currently offline  Sergey Kononov
Messages: 11
Registered: September 2014
occasional visitor

From: *54.237.84.nsu.ru
Dear Tobias,

Thank you for the usefull info. We'll look into the MVD software.

Cheers,
Sergey


Cheers,
Sergey
Previous Topic: Secondary neutron production threshold
Next Topic: gphysi.dat
Goto Forum:
  


Current Time: Fri Nov 01 20:22:10 CET 2024

Total time taken to generate the page: 0.00739 seconds