Re: [Fwd: Neutron capture in Gd-doped liquid scintillator.]

Dear Ivana
many thanks, we will try this. I will let you know about the results. Best regards, Dmitry

Ivana Hrivnacova пишет:
> Dear Dmitry,
>
> In VMC, there is no automatic scoring like the
> screen output in Fluka which would work generally
> (also for G4).
> You have to combine the methods for access to secondaries to get
> the info about the residual nuclei.
> The following 3 methods will give you an info about all secondaries
> at each step of your beam particle:
> virtual Int_t NSecondaries() const=0;
> virtual void GetSecondary(Int_t isec, Int_t& particleId,
> TLorentzVector& position,
> TLorentzVector& momentum) =0;
> virtual TMCProcess ProdProcess(Int_t isec) const =0;
>
> Or you can account the info when tracking the secondary
> particles (nuclei) using the methods:
> // Return the PDG of the particle transported
> virtual Int_t TrackPid() const =0;
>
> // Return the charge of the track currently transported
> virtual Double_t TrackCharge() const =0;
> ...
> virtual void TrackPosition(Double_t &x, Double_t &y, Double_t
> &z) const =0;
> virtual void TrackMomentum(Double_t &px, Double_t &py, Double_t
> &pz, Double_t &etot) const =0;
> ...
>
> The common defined in fluka_vmc/source/Ffheavy.h
> is interfaced in TFluka via GetSecondary() method.
>
> I suggest you to try to modify one of the VMC examples
> for your case with a very simple geometry first and
> try to see if you are able to combine the existing VMC methods
> to get all info about neclei that you need.
>
> Best regards,
>
> Ivana
>
>
>
> Dmitry Naumov wrote:
>> Dear Ivana,
>>
>> Concerning FLUKA there is a common block RESNUC (see its description
>> in "fluka2008.3/flukapro/(RESNUC)" file ) and the screen output which
>> can be obtained as explained here
>> http://www.fluka.org/fluka.php?id=man_onl&sub=97#RES_put
>>
>> What we would like to access is the following. Having a beam
>> particle to know:
>> - target particle
>> - list of secondaries including nuclei
>> - process occured
>>
>> As an example we want to simulate neutron propagation through Gd
>> dopped liquid scintillator with help of both G4 and FLUKA (thus -
>> VMC). We seek for a Gd capture of neutrons, produced gammas, nuclei
>> remnants. If this info is already available it would be just great!
>> Perhaps, this is present (fully, or partially in
>> fluka_vmc/source/Ffheavy.h, but I do not know about a general VMC
>> access to this info).
>>
>> Kind regards, Dmitry
>>
>>
>> Ivana Hrivnacova пишет:
>>> Dear Dmitry,
>>>
>>> There is now a reply to my question in the forum;
>>> I include it also below in this mail.
>>> In principle, this procedure could be done via VMC,
>>> as we have interfaced the necessary functions,
>>> like the creator process for a secondary etc.
>>> Could you explain in more detail what exactly
>>> do you need and which are the access functions
>>> in Fluka that you want to use and which you do
>>> not find in VMC?
>>>
>>> Thank you,
>>>
>>> Ivana
>>>
>>> *** Discussion title: Hadronic Processes
>>>
>>> There is no easy way to do this, but you could try:
>>>
>>> 1) Loop over tracks after the capture and test that they are created by
>>> the capture process. Do this by calling track->GetCreatorProcess() and
>>> checking the process name.
>>>
>>> 2) Then get the particle by doing track->GetDynamicParticle().
>>>
>>> 3) Once you have the dynamic particle you can check that the track is a
>>> nucleus, and also look at things like mass, charge, momentum, etc.
>>>
>>>
>>> Dmitry Naumov wrote:
>>>> Dear colleagues,
>>>>
>>>> I am not an expert in Fluka, but it seems that Fluka has these
>>>> access functions, what is left is to organize vmc like interface to
>>>> them. If G4 does not have these functions then perhaps vmc
>>>> functions should return something like 'false' value.
>>>> Best regards, Dmitry
>>>>
>>>> Ivana Hrivnacova пишет:
>>>>> Hi,
>>>>>
>>>>> I don't know. I have now posted the question in G4 hypernews:
>>>>> http://hypernews.slac.stanford.edu/HyperNews/geant4/get/hadronprocess/973.html
>>>>>
>>>>> And what about Fluka?
>>>>>
>>>>> Best regards,
>>>>>
>>>>> Ivana
>>>>>
>>>>> Carminati Federico wrote:
>>>>>> Hello Iliya,
>>>>>> I frankly ignore whether this information is available in G4.
>>>>>> Ivana, do you know? Best regards,
>>>>>> Federico Carminati
>>>>>> CERN-PH 1211 Geneva 23
>>>>>> Switzerland
>>>>>> Tel: +41 22 76 7495:
>>>>>
>>>>>> Fax: +41 22 76 68505
>>>>>> Mobile: +41 76 487 4843
>>>>>>
>>>>>>
>>>>>>
>>>>>> On 24 Dec 2008, at 16:05, Illya I. Shapoval wrote:
>>>>>>
>>>>>>> Dear VMC team! I need an effective technic for the neutron
>>>>>>> capture and final nucleus state monitoring in my set of tasks
>>>>>>> too. It would be realy interesting to hear an answer on Maxim's
>>>>>>> question.
>>>>>>>
>>>>>>> With respect,
>>>>>>> Illya Shapoval
>>>>>>> ----- Original Message -----
>>>>>>> *From:* Dmitry Naumov <mailto:naumov_at_nusun.jinr.ru>
>>>>>>> *To:* vmc_at_lxbuild091.cern.ch <mailto:vmc_at_lxbuild091.cern.ch>
>>>>>>> *Sent:* Wednesday, December 24, 2008 10:57 AM
>>>>>>> *Subject:* [VMC] [Fwd: Neutron capture in Gd-doped liquid
>>>>>>> scintillator.]
>>>>>>>
>>>>>>>
>>>>>>> Dear VMC team it seems that this message was either lost and
>>>>>>> did
>>>>>>> not arrived to you or left without an asnwer. As we need an
>>>>>>> answer
>>>>>>> could anyone please address it once? Thanks! Dmitry
>>>>>>>
>>>>>>>
>>>>>>> -------- Исходное сообщение --------
>>>>>>> Тема: Neutron capture in Gd-doped liquid scintillator.
>>>>>>> Дата: Wed, 26 Nov 2008 17:17:56 +0300
>>>>>>> От: Maxim Gonchar <gmaxfl_at_gmail.com>
>>>>>>> Отвечать: gmaxfl_at_gmail.com
>>>>>>> Кому: vmc_at_root.cern.ch
>>>>>>> Копия: Dmitry Naumov <naumov_at_nusun.jinr.ru>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Hello,
>>>>>>>
>>>>>>> I want to simulate neutron propagation through the Gd-doped
>>>>>>> liquid scintillator. Is there any way to get information
>>>>>>> about the nucleus that captures neutron?
>>>>>>>
>>>>>>> Thank you in advance.
>>>>>>>
>>>>>>> regards, Maxim Gonchar
>>>>>>> mailto:gmaxfl_at_gmail.com
>>>>>>
>>>>>
>>>
>
Received on Wed Jan 14 2009 - 07:29:33 CET