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 >>>>> >>>> >>
-- ============================================================== e-mail: Ivana.Hrivnacova_at_cern.ch address: Institut Physique Nucleaire, 91406 Orsay, France phone: (33) 01 69 15 65 94 ==============================================================Received on Tue Jan 13 2009 - 17:30:56 CET
This archive was generated by hypermail 2.2.0 : Wed Jan 14 2009 - 11:25:02 CET