Message: antiproton inelastic scattering cross section Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question antiproton inelastic scattering cross section 

Forum: Hadronic Processes
Date: 24 Aug, 2009
From: Tsuguo Aramaki <Tsuguo Aramaki>

Hello. 
I've been trying to modify the antiproton inelastic scattering cross section (especially for low energy) based on the paper,  Duperray et al. (2005). I simply added some lines at line 1566 in geant4/source/processes/hadronic/cross_sections/src/G4HadronCrossSections.cc and recompiled it, but it didn't work...  
*****
<line 1561>
// Apply correction factors
      sigin = crin*(sigin + sigel)*std::pow(a, alph);
      sigel = crel*0.36*sigel*std::pow(a, 1.17);
      sigel = sigel*partel[ipart1];
      sigin = sigin*partin[ipart1];
<added lines here, line 1566>
      if (ipart == 15) // antiproton
      {
        G4double m0 = aParticle->GetMass()/GeV;	
	G4double T = std::sqrt(m0*m0 + p*p) - m0;
	sigin = 661.0*(1.0+0.0036*std::pow(T,-0.774)-0.902*std::pow(T, 0.0151))*std::pow(a,0.667);	  
      }
.....
*****
I actually don't know if this is the correct way, so could anyone please give me some advice?
Thank you!
Tsuguo

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

1 Idea: Re: antiproton inelastic scattering cross section   (Vladimir Ivanchenko - 03 Sep, 2009)
(_ None: Re: antiproton inelastic scattering cross section   (Tsuguo Aramaki - 03 Sep, 2009)
(_ Question: Re: antiproton inelastic scattering cross section   (Huaimin Liu - 26 Sep, 2009)
 Add Message Add Message
to: "antiproton inelastic scattering cross section"

 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 ]