Message: HadronicProcessStore - how to get Isotopic Cross Sections? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None HadronicProcessStore - how to get Isotopic Cross Sections? 

Forum: Hadronic Processes
Date: 29 May, 2014
From: Moritz Kütt <Moritz Kütt>

Hi,

I'm trying to extract specific isotopic cross sections from the G4NDL set, mainly for plotting.

When I use the G4HadronicProcessStore::GetInelasticCrossSectionPerIsotope(...) method it always returns 0, even if the PerAtom method returns a value. As far as I understand both methods, the latter should be the weighted sum of all natural occuring isotope cross sections?

I tried this by adding a column to the Hadr00 example.

    if(fVerbose>0) G4cout << std::setw(12) << xs/barn;  
    xs = store->GetInelasticCrossSectionPerAtom(particle,e,elm,mat);
    xtot += xs;
    if(fVerbose>0) G4cout << " " << std::setw(12) << xs/barn;  
    xs = store->GetInelasticCrossSectionPerIsotope(particle,e,13,27);

It always returns 0 (see attached file). Similarly for other elements/isotopes.

What am I'm doing wrong?

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/05/29/07.06-36822-forum.txt

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Idea: Re: HadronicProcessStore - how to get Isotopic Cross Sections?   (Vladimir Ivanchenko - 02 Jun, 2014)
(_ More: Re: HadronicProcessStore - how to get Isotopic Cross Sections?   (Gumplinger Peter - 19 Jan, 2016)
 Add Message Add Message
to: "HadronicProcessStore - how to get Isotopic Cross Sections?"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]