Message: Re: tot (p+12C) non elatic XS is very law Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: tot (p+12C) non elatic XS is very law 

Forum: Hadronic Processes
Re: None tot (p+12C) non elatic XS is very law (Jeyasingam Jeyasugiththan)
Re: None Re: tot (p+12C) non elatic XS is very law (michel maire)
Date: 18 Jan, 2013
From: Jeyasingam Jeyasugiththan <Jeyasingam Jeyasugiththan>

Hi michel,

Thank you very much for your mail with details. I already deactivated the EM and elastic process etc using /process/inactivate and only proton inelastic is active.

      ____________________________________________________________________

                            CROSS SECTION MEASUREMENTS                    
      ____________________________________________________________________

	 The run consists of 100000 proton of 7 MeV through 20 m   of Carbon (density: 2.267 g/cm3 )

	 Process calls frequency --->	ProtonInelastic= 100000

	 fTotalcount	100000

	 MeanFreePath:	1.2883 m   +- 1.2892 m  		massic: 292.06 g/cm2 
	 CrossSection:	0.0077621 cm^-1 			massic: 342.4 um2/mg
	 CrossSection:	68.218 mb 

	 Verification : crossSections from G4HadronicProcessStore:
	 ProtonInelastic: Cross Section = 67.989 mb   massic XS = 341.25 um2/mg

	 total massic XS = 341.25 um2/mg

with regards Sugi

On Fri, 18 Jan 2013 14:26:03 GMT, michel maire wrote:

> On Thu, 17 Jan 2013 17:20:29 GMT, Jeyasingam Jeyasugiththan wrote:
> 
> > I use Hadr03 example to calculate (proton+Carbon) total non-elastic reaction cross section.  An event is killed at the first interaction of the incident Proton. Target is 12Carbon (infinite length to be sure that interaction
> >  occurs at any events ).  First 
> >  mean free path is calculated and then XS/volume is calculated by (1/meanfreepath). Finally XS/Atom is calculated using following equation
> > 
> >  XS/Atom =(XS/Volume)*(M/(density*Avogadro number))
> > 
> > My Physics list
> > 
> >     AddPhysicsList("standard_opt3");
> >     hadronPhys.push_back( new HadronPhysicsQGSP_BIC());
> >     hadronPhys.push_back( new G4EmExtraPhysics());
> >     hadronPhys.push_back( new G4HadronElasticPhysics());
> >     hadronPhys.push_back( new G4StoppingPhysics());
> >     hadronPhys.push_back( new G4IonBinaryCascadePhysics());
> >     hadronPhys.push_back( new G4NeutronTrackingCut());
> > 
> >  XS/Atom values in the low energy are very lower than the actual values (ple see the attachment). 
> > 
> > The XS/atom is strongly depend on density of target material. How can I
> > solve it? I would greatly appreciate your helps.
> > 
> >    Attachment:
> >       http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/01/17/09.18-75950-XS.pdf
> > 
> > Thanks in advance
> > 
> > with regards Sugi
> > 
> 
>  As explained in Readme, the electromagnetic processes must not be registered in Hadr03 context: no energy loss or multiple scattering before hadronic interaction.
>  Indeed, only QGSP_BIC (and ElasticPhysics) are relevant here.
>  For the rest, I let experts in hadronic physics judge the coherence of your physics list.
>  I attach results of 'default' Hadr03, with G4HadronInelasticQBBC physics constructor.
>       http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/01/18/06.22-52853-sugi.mac.txt
>       http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2013/01/18/06.22-28669-sugi.log.txt
> 

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

1 None: Re: tot (p+12C) non elatic XS is very law   (michel maire - 18 Jan, 2013)
(_ Question: Re: tot (p+12C) non elatic XS is very law   (Jeyasingam Jeyasugiththan - 22 Jan, 2013)
(_ None: Re: tot (p+12C) non elatic XS is very law   (michel maire - 22 Jan, 2013)
 Add Message Add Message
to: "Re: tot (p+12C) non elatic XS is very law"

 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 ]