Message: getting zero cross section Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question getting zero cross section 

Forum: Hadronic Processes
Date: 02 Jul, 2011
From: Summit <Summit>

Hello experts

To calculate the cross-section, I have included the following:

G4Material* material = detector->GetAbsorMaterial();

G4String particle = kinematic->GetParticleGun()->GetParticleDefinition()->GetParticleName();

G4HadronicProcessStore* store = G4HadronicProcessStore::Instance();
G4double ener, cs; 
ener=100*MeV;
const G4ParticleDefinition* particle_name = 
    G4ParticleTable::GetParticleTable()->FindParticle(particle);
G4cout<<"........................................................................\n";
G4cout<<std::setw(5) << "N" << std::setw(18)<<"Energy(MeV)"<<std::setw(18)<<"Cross section"<<"\n";
G4cout<<"........................................................................\n";
for(int i=1; i<=500; i++)
{ 
	G4cout << std::setw(5) << i << std::setw(18) << ener;  
    	cs = store->GetInelasticCrossSectionPerVolume(particle_name,ener,material);
	G4cout<<std::setw(18)<<cs<<"\n";
	ener+=100*MeV;
}

 but I am getting zero cross section as
........................................................................
    N       Energy(MeV)     Cross section
........................................................................
    1               100                 0
    2               200                 0
    3               300                 0
    4               400                 0
    5               500                 0
    6               600                 0
    7               700                 0
...
...
...

Can anyone suggest me what should I do to calculate cross section?

Thanks in advance,

Summit

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

1 Idea: Re: getting zero cross section   (Vladimir Ivanchenko - 06 Jul, 2011)
(_ None: Re: getting zero cross section   (Summit - 08 Jul, 2011)
 Add Message Add Message
to: "getting zero 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 ]