GSI Forum: Monte Carlo Engines

Home » PANDA » PandaRoot » Monte Carlo Engines » NIEL

Show: Today's Messages :: Polls :: Message Navigator

NIEL [message #22449]

Wed, 11 April 2018 06:31

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

Report message to a moderator

Re: NIEL [message #22468 is a reply to message #22449]

Thu, 12 April 2018 13:29

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

Report message to a moderator

Re: NIEL [message #22469 is a reply to message #22468]

Thu, 12 April 2018 19:50

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

Report message to a moderator

Re: NIEL [message #22473 is a reply to message #22469]

Fri, 13 April 2018 03:52

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

Report message to a moderator

Re: NIEL [message #22501 is a reply to message #22449]

Tue, 17 April 2018 16:55

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

Report message to a moderator