Message: I can't simulate B10(n,alpha)Li7 Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question I can't simulate B10(n,alpha)Li7 

Forum: Hadronic Processes
Date: 17 Sep, 2008
From: <csi505>

Hi :

Recently I have a simulation for B10(n,alpha)

Li7 reactions with Geat4.9.1. Two weeks ago, I accounted a problem , I found the energy of alpha that I simulated is wrong. The following is the

result:

Idle> /run/beamOn 10
### Run 11 start
WARNING: Scene "refresh" not found.
  /vis/scene/list to see scenes.
EventAction: Event # 0 started

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0 0.000282         0        0  0.00178       200      Target NeutronInelastic
    :----- List of 2ndaries - #SpawnInStep=  2(Rest= 0,Along= 0,Post= 2), #SpawnTotal=  2 ---------------
    :         0         0  0.000282      1.75              alpha
    :         0         0  0.000282         1           Li7[0.0]
    :----------------------------------------------------------------- EndOf2ndaries Info ---------------

*********************************************************************************************************
* G4Track Information:   Particle = Li7[0.0],   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0 0.000282         1        0        0         0      Target initStep
    1  0.00139 -0.000338 0.000272         0        1  0.00143   0.00143      Target ionIoni

*********************************************************************************************************
* G4Track Information:   Particle = alpha,   Track ID = 2,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0 0.000282      1.75        0        0         0      Target initStep
    1 -0.00342 0.000832 0.000306         0     1.75  0.00352   0.00352      Target msc

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0     -200  2.53e-08        0        0         0       World initStep
    1        0        0     -150  2.53e-08        0       50        50       Check Transportation
    2        0        0  -0.0015  2.53e-08        0      150       200      Target Transportation
    3        0        0   0.0015  2.53e-08        0    0.003       200       Check Transportation
    4        0        0      150  2.53e-08        0      150       350       World Transportation
    5        0        0      200  2.53e-08        0       50       400  OutOfWorld Transportation
RunAction: End of run actions are started

As you see, the energy of alpha is 1.75 Mev and the process which created alpha is NeutronInelastic. But the fact is 
                 1. B10 + n --> Li7(1.02Mev) + alpha(1.78Mev)   (percent:6%)
                 2. B10 + n --> Li7(0.84Mev) + alpha(1.47Mev)    (percent:94%)

So it is wrong.Here I'll give you my result of 1000000 events:

Idle> /run/beamOn 1000000
### Run 22 start
WARNING: Scene "refresh" not found.
  /vis/scene/list to see scenes.
EventAction: Event # 0 started
EventAction: Event # 10000 started
EventAction: Event # 20000 started
EventAction: Event # 30000 started
EventAction: Event # 40000 started
EventAction: Event # 50000 started
EventAction: Event # 60000 started
EventAction: Event # 70000 started
EventAction: Event # 80000 started
EventAction: Event # 90000 started
EventAction: Event # 100000 started
EventAction: Event # 110000 started
EventAction: Event # 120000 started
EventAction: Event # 130000 started
EventAction: Event # 140000 started
EventAction: Event # 150000 started
EventAction: Event # 160000 started
EventAction: Event # 170000 started
EventAction: Event # 180000 started
EventAction: Event # 190000 started
EventAction: Event # 200000 started
EventAction: Event # 210000 started
EventAction: Event # 220000 started
EventAction: Event # 230000 started
EventAction: Event # 240000 started
EventAction: Event # 250000 started
EventAction: Event # 260000 started
EventAction: Event # 270000 started
EventAction: Event # 280000 started
EventAction: Event # 290000 started
EventAction: Event # 300000 started
EventAction: Event # 310000 started
EventAction: Event # 320000 started
EventAction: Event # 330000 started
EventAction: Event # 340000 started
EventAction: Event # 350000 started
EventAction: Event # 360000 started
EventAction: Event # 370000 started
EventAction: Event # 380000 started
EventAction: Event # 390000 started
EventAction: Event # 400000 started
EventAction: Event # 410000 started
EventAction: Event # 420000 started
EventAction: Event # 430000 started
EventAction: Event # 440000 started
EventAction: Event # 450000 started
EventAction: Event # 460000 started
EventAction: Event # 470000 started
EventAction: Event # 480000 started
EventAction: Event # 490000 started
EventAction: Event # 500000 started
EventAction: Event # 510000 started
EventAction: Event # 520000 started
EventAction: Event # 530000 started
EventAction: Event # 540000 started
EventAction: Event # 550000 started
EventAction: Event # 560000 started
EventAction: Event # 570000 started
EventAction: Event # 580000 started
EventAction: Event # 590000 started
EventAction: Event # 600000 started
EventAction: Event # 610000 started
EventAction: Event # 620000 started
EventAction: Event # 630000 started
EventAction: Event # 640000 started
EventAction: Event # 650000 started
EventAction: Event # 660000 started
EventAction: Event # 670000 started
EventAction: Event # 680000 started
EventAction: Event # 690000 started
EventAction: Event # 700000 started
EventAction: Event # 710000 started
EventAction: Event # 720000 started
EventAction: Event # 730000 started
EventAction: Event # 740000 started
EventAction: Event # 750000 started
EventAction: Event # 760000 started
EventAction: Event # 770000 started
EventAction: Event # 780000 started
EventAction: Event # 790000 started
EventAction: Event # 800000 started
EventAction: Event # 810000 started
EventAction: Event # 820000 started
EventAction: Event # 830000 started
EventAction: Event # 840000 started
EventAction: Event # 850000 started
EventAction: Event # 860000 started
EventAction: Event # 870000 started
EventAction: Event # 880000 started
EventAction: Event # 890000 started
EventAction: Event # 900000 started
EventAction: Event # 910000 started
EventAction: Event # 920000 started
EventAction: Event # 930000 started
EventAction: Event # 940000 started
EventAction: Event # 950000 started
EventAction: Event # 960000 started
EventAction: Event # 970000 started
EventAction: Event # 980000 started
EventAction: Event # 990000 started
RunAction: End of run actions are started
ERROR: G4VisCommandsViewerUpdate::SetNewValue: no current viewer.
HistoManager: End of run actions are started
========================================================
Beam particle                        neutron
Beam Energy(MeV)                     2.53e-08
Number of events                     1000000
Average energy deposit (MeV)         0.06271   RMS(MeV) 0.37
Average number of steps              1.061
Average number of gamma              0.002398
Average number of e-                 5.2e-05
Average number of e+                 0
Average number of neutrons           0
Average number of protons            0
Average number of antiprotons        0
Average number of pi+ & pi-          0
Average number of pi0                0
Average number of kaons              0
Average number of muons              0
Average number of deuterons+tritons  0
Average number of He3+alpha          0.03025
Average number of ions               0.03096
Average number of forward neutrons   0
Average number of reflected neutrons 0
Average number of leaked neutrons    0
Average number of proton leak        0
Average number of pion leak          0
========================================================

So you can see that there was few Gamma, the reaction I simulated with Geant4 is wrong , I don't know why.
 My program is modifyed from examples/extended/hadronic/Hadr01, now I'll show my Phisylist :

void PhysicsList::ConstructNH()
{
   theParticleIterator->reset();
   while( (*theParticleIterator)() ){
   G4ParticleDefinition* particle = theParticleIterator->value();
   G4ProcessManager* pManager = particle->GetProcessManager();
   G4String particleName = particle->GetParticleName();
   // G4ProcessManager * pManager = G4Neutron::Neutron()->GetProcessManager();
   if(particleName=="neutron"){
//ELASTIC
G4HadronElasticProcess* thenElasticProcess 
                         = new G4HadronElasticProcess();
   G4NeutronHPElastic* thenHPElasticModel = new G4NeutronHPElastic();
   thenHPElasticModel->SetMinEnergy(4.*eV);
   thenHPElasticModel->SetMaxEnergy(20.*MeV);
   G4LElastic* thenElasticModel = new G4LElastic();
   thenElasticModel->SetMinEnergy(20.*MeV);
   G4NeutronHPElasticData* theNeutronData1 = new G4NeutronHPElasticData;
   thenElasticProcess->AddDataSet(theNeutronData1);
//INELASTIC
G4NeutronInelasticProcess* theNeutronInelasticProcess 
                         = new G4NeutronInelasticProcess(); 
   G4NeutronHPInelastic* theNeutronHPInelastic = new G4NeutronHPInelastic();
   theNeutronHPInelastic->SetMaxEnergy(20.*MeV);
   //Binary-Cascade
   /*G4BinaryCascade* theNeutronBCModel = new G4BinaryCascade();
   theNeutronBCModel->SetMinEnergy(20.*MeV);
   theNeutronBCModel->SetMaxEnergy(5.*GeV);*/
   //Bertini-Cacscade
   G4CascadeInterface* theBertiniModel = new G4CascadeInterface();
   theBertiniModel->SetMinEnergy(20.*MeV);
   theBertiniModel->SetMaxEnergy( 5.*GeV);

   G4LENeutronInelastic* theNeutronLEPModel = new G4LENeutronInelastic();
   theNeutronLEPModel->SetMinEnergy(5.*GeV);    //geht bis 20GeV
   theNeutronLEPModel->SetMaxEnergy(20.*GeV);
   G4HENeutronInelastic* theNeutronHEModel = new G4HENeutronInelastic();
   theNeutronHEModel->SetMinEnergy(20.*GeV);
   G4NeutronHPInelasticData* theNeutronData2 = new G4NeutronHPInelasticData;
   theNeutronInelasticProcess->AddDataSet(theNeutronData2);

   //theNeutronInelasticProcess->RegisterMe(theNeutronBCModel);
   theNeutronInelasticProcess->RegisterMe(theBertiniModel);
   theNeutronInelasticProcess->RegisterMe(theNeutronLEPModel);
   theNeutronInelasticProcess->RegisterMe(theNeutronHEModel);
   theNeutronInelasticProcess->RegisterMe(theNeutronHPInelastic);
   pManager->AddDiscreteProcess(theNeutronInelasticProcess);

   G4HadronFissionProcess* thenFission
                         = new G4HadronFissionProcess();
   G4LFission* thenFissionModel = new G4LFission();
   thenFission->RegisterMe(thenFissionModel);
   thenFissionModel->SetMinEnergy(20.*MeV);
   G4NeutronHPFission* thenHPFissionModel = new G4NeutronHPFission();
   thenHPFissionModel->SetMaxEnergy(20.*MeV);
   G4NeutronHPFissionData* theNeutronData4 = new G4NeutronHPFissionData;
   thenFission->AddDataSet(theNeutronData4);
   thenFission->RegisterMe(thenHPFissionModel);
   pManager->AddDiscreteProcess(thenFission);

   G4HadronCaptureProcess* thenCapture
                         = new G4HadronCaptureProcess();
   G4LCapture* thenCaptureModel = new G4LCapture();
   thenCaptureModel->SetMinEnergy(20.*MeV);
   G4NeutronHPCapture* thenHPCaptureModel = new G4NeutronHPCapture();
   thenHPCaptureModel->SetMaxEnergy(20.*MeV);
   G4NeutronHPCaptureData* theNeutronData3 = new G4NeutronHPCaptureData;
   thenCapture->AddDataSet(theNeutronData3);

   thenCapture->RegisterMe(thenCaptureModel);
   thenCapture->RegisterMe(thenHPCaptureModel);
   pManager->AddDiscreteProcess(thenCapture);
}
}
}

The energy Of n is 0.0253 Mev , the following is my run words:

/testhadr/TargetMat G4_B /testhadr/TargetRadius 3 cm /testhadr/TargetLength 3 um /testhadr/NumberDivZ 1 /testhadr/PrintModulo 10000 /testhadr/DrawTracks all /testhadr/CutsAll 1. mm /run/initialize /gun/particle neutron /gun/energy 0.0253 eV /run/beamOn 100

Two weeks passed away, I still can't sovle this problem . So I thick it may be caused by the Geant4 itself . I'm looking foward to the answer. The file I added is my program.

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

1 Idea: Re: I can't simulate B10(n,alpha)Li7   (Dundar Ucar - 17 Sep, 2008)
1 None: Re: I can't simulate B10(n,alpha)Li7   (Koi, Tatsumi - 17 Sep, 2008)
2 None: Re: I can't simulate B10(n,alpha)Li7   (csi505 - 17 Sep, 2008)
1 None: Re: I can't simulate B10(n,alpha)Li7   (Koi, Tatsumi - 18 Sep, 2008)
 Add Message Add Message
to: "I can't simulate B10(n,alpha)Li7"

 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 ]