Message: Scoring double differential cross section neutron production in 800 MeV p + 2cm Pb Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Scoring double differential cross section neutron production in 800 MeV p + 2cm Pb 

Forum: Biasing and Scoring
Date: 10 Oct, 2011
From: Sylt.Zhang <Sylt.Zhang>

Dear Geant4 experts,
   I want to score the  double differential cross section  of the neutron production at the given angle( 0.8GeV proton + 2cm Pb(target) at 30 degree emission angle), and I think it can be converted from the energy spectrum dN/dE.
So I using following code:

 const G4TrackVector* secondary = fpSteppingManager->GetSecondary();
  for (size_t lp=0; lp<(*secondary).size(); lp++) {
    G4String  particlename = (*secondary)[lp]->GetDefinition()->GetParticleName();
    if (particlename=="neutron"){
      G4double theta = (*secondary)[lp]->GetMomentumDirection().theta()/degree;
      if(29.5 <= theta && theta <=30.5){
	G4double energy = (*secondary)[lp]->GetKineticEnergy();
	myHitsData.CS=energy;  //filling the Tree

The result didn't agree with the Fluka simulation and experimental data,but the shape is much more like.Are there any scaling factor ?
  the attachment is my results, and I scale the Geant4 0.01. I don't know why gean4 more larger ?

 Looking forward to receive your help!

Sylt.Zhang

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/10/10/23.54-93890-0.8GeVp2cmPb_30.png

 Add Message Add Message
to: "Scoring double differential cross section neutron production in 800 MeV p + 2cm Pb"

 Subscribe Subscribe

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