Message: different results for different Hadron Inelastic Physics Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None different results for different Hadron Inelastic Physics 

Forum: Hadronic Processes
Date: 31 Aug, 2016
From: riazi <riazi>

Dear Sir/Madam
I run the example Hadr06 with 4 different Hadron Inelastic Physics using Geant4.10.02. I changed the target material to Li7 by activation of command /testhadr/det/setIsotopeMat Li7 3 7 0.535 g/cm3 in run1.mac file and I changed the source to proton with following two commands in run1.mac file 
 /gun/particle proton
/gun/energy 5 MeV 

I observed that the generated particles were different for four physics. for example there is no neutron for G4HadronPhysicsQGSP_BIC_HP physics.

The results for each of them is reported below.

1) for Physic List G4HadronPhysicsFTFP_BERT_HP the results is: The run is 1000000 proton of 5 MeV through 5 mm of Li7 (density: 535 mg/cm3)

 Process calls frequency :
     RadioactiveDecay=   3597        Transportation=    290                 compt=      4
                eIoni=      8                 hIoni= 999950            hadElastic=   3586
              ionIoni=   3679       protonInelastic=     53 

 List of generated particles:
            Be7:       4  Emean =  599.54 keV	( 399.77 keV --> 748.55 keV)
            Li7:    3593  Emean =  457.04 keV	( 3.0891 eV  --> 2.1779 MeV)
          alpha:      82  Emean =  6.0271 MeV	( 1.6711 MeV --> 13.096 MeV)
             e-:       8  Emean =  6.7369 keV	( 58.885 eV  --> 19.923 keV)
          gamma:     137  Emean =  2.5616 MeV	( 68.095 keV --> 15.946 MeV)
        neutron:       4  Emean =  944.43 keV	( 138.38 keV --> 1.8739 MeV)
           nu_e:       4  Emean =  861.78 keV	( 861.78 keV --> 861.78 keV)
         proton:       7  Emean =  1.9699 MeV	( 794.39 keV --> 2.7427 MeV)

 Mean energy deposit per event = 5.0003 MeV;  rms = 67.366 keV
 Mean energy flow per event    = 358.11 eV ;  rms = 63.56 keV

 List of particles leaving the absorber :
          gamma:     137  Emean =  2.5612 MeV	( 65.057 keV --> 15.946 MeV) 	Eflow/event = 350.88 eV 
        neutron:       4  Emean =  944.43 keV	( 138.38 keV --> 1.8739 MeV) 	Eflow/event = 3.7777 eV 
           nu_e:       4  Emean =  861.78 keV	( 861.78 keV --> 861.78 keV) 	Eflow/event = 3.4471 eV 

2) For Physics list G4HadronPhysicsQGSP_BIC_HP the results is:

The run is 1000000 proton of 5 MeV through 5 mm of Li7 (density: 535 mg/cm3)

 Process calls frequency :
     RadioactiveDecay=   3586                 hIoni= 999943            hadElastic=   3590
              ionIoni=   3690       protonInelastic=     53 

 List of generated particles:
            Li7:    3586  Emean =  457.23 keV	( 3.0891 eV  --> 2.1779 MeV)
          alpha:     104  Emean =  10.303 MeV	( 7.5011 MeV --> 14.746 MeV)

 Mean energy deposit per event = 5.0009 MeV;  rms = 125.09 keV
 Mean energy flow per event    = 0 meV;  rms = 0 meV

3) For Physics list G4HadronInelasticQBBC the results is: The run is 1000000 proton of 5 MeV through 5 mm of Li7 (density: 535 mg/cm3)

 Process calls frequency :
     RadioactiveDecay=   3586                 hIoni= 999943            hadElastic=   3590
              ionIoni=   3690       protonInelastic=     53 

 List of generated particles:
            Li7:    3586  Emean =  457.23 keV	( 3.0891 eV  --> 2.1779 MeV)
          alpha:     104  Emean =  10.303 MeV	( 7.5011 MeV --> 14.746 MeV)

 Mean energy deposit per event = 5.0009 MeV;  rms = 125.09 keV
 Mean energy flow per event    = 0 meV;  rms = 0 meV

4) For Physics list G4HadronPhysicsINCLXX the results is:
 List of generated particles:
            Be7:       1  Emean =   677.3 keV	( 677.3 keV --> 677.3 keV)
            Li7:    3589  Emean =  458.73 keV	( 3.0891 eV  --> 2.1779 MeV)
   Li7[477.612]:       1  Emean =  11.296 eV 	( 11.296 eV  --> 11.296 eV )
          alpha:     102  Emean =   10.28 MeV	( 7.4558 MeV --> 14.441 MeV)
             e-:       1  Emean =  60.227 eV 	( 60.227 eV  --> 60.227 eV )
          gamma:       3  Emean =  889.33 keV	( 477.59 keV --> 1.4796 MeV)
        neutron:       1  Emean =  310.71 keV	( 310.71 keV --> 310.71 keV)
           nu_e:       1  Emean =  384.21 keV	( 384.21 keV --> 384.21 keV)
         proton:       1  Emean =  2.8502 MeV	( 2.8502 MeV --> 2.8502 MeV)

 Mean energy deposit per event = 5.0009 MeV;  rms = 123.91 keV
 Mean energy flow per event    = 3.3629 eV ;  rms = 2.395 keV

 List of particles leaving the absorber :
          gamma:       3  Emean =  889.33 keV	( 477.59 keV --> 1.4796 MeV) 	Eflow/event = 2.668 eV 
        neutron:       1  Emean =  310.71 keV	( 310.71 keV --> 310.71 keV) 	Eflow/event = 310.71 meV
           nu_e:       1  Emean =  384.21 keV	( 384.21 keV --> 384.21 keV) 	Eflow/event = 384.21 meV

With Best Regards, Mobarak

 Add Message Add Message
to: "different results for different Hadron Inelastic Physics"

 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 ]