Message: Re: neutron_HP_Physics and NIST_Material Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: neutron_HP_Physics and NIST_Material 

Forum: Hadronic Processes
Re: None neutron_HP_Physics and NIST_Material (Sylvia Studeny)
Re: None Re: neutron_HP_Physics and NIST_Material (Vladimir Ivantchenko )
Re: None Re: neutron_HP_Physics and NIST_Material (Sylvia Studeny)
Re: None Re: neutron_HP_Physics and NIST_Material (Vladimir Ivantchenko )
Re: None Re: neutron_HP_Physics and NIST_Material (Sylvia Studeny)
Date: 08 May, 2006
From: Vladimir Ivantchenko <vnivanch@mail.cern.ch>

On Mon, 8 May 2006, Sylvia Studeny wrote:

> *** Discussion title: Hadronic Processes
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/hadronprocess/491/1/1/1/1"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> Hi,
> 
> I changed the hadron physic list (concerning n) to:
> 
>  // neutron Physics
>    pManager = G4Neutron::Neutron()->GetProcessManager();
> 
>    // add process
>    G4HadronElasticProcess* thenElasticProcess 
>                          = new G4HadronElasticProcess();
>    thenElasticProcess->AddDataSet(new G4NeutronHPElasticData); 
> 
>    G4LElastic* thenElasticModel = new G4LElastic();
>   // thenElasticModel->SetMinEnergy(20.*MeV);
>  //  G4NeutronHPElastic* theNeutronHPElasticModel= new G4NeutronHPElastic();
>  //  theNeutronHPElasticModel ->SetMaxEnergy(20.*MeV);
>    thenElasticProcess->RegisterMe(thenElasticModel);
>  //  thenElasticProcess->RegisterMe(theNeutronHPElasticModel);
> 
>    pManager->AddDiscreteProcess(thenElasticProcess);
> 
>    G4NeutronInelasticProcess* theNeutronInelasticProcess 
>                          = new G4NeutronInelasticProcess();
>    theNeutronInelasticProcess->AddDataSet(new G4NeutronHPInelasticData);  
> 
>   // G4NeutronHPInelastic* theNeutronHPInelasticModel= new G4NeutronHPInelastic();
>   // theNeutronHPInelasticModel ->SetMaxEnergy(20.*MeV);
>    G4BinaryCascade* theNeutronBCModel = new G4BinaryCascade();
>  //  theNeutronBCModel->SetMinEnergy(20.*MeV);
>    theNeutronBCModel->SetMaxEnergy(6.*GeV);
>    G4LENeutronInelastic* theNeutronLEPModel = new G4LENeutronInelastic();
>    theNeutronLEPModel->SetMinEnergy(4.*GeV);
>    G4HENeutronInelastic* theNeutronHEPModel = new G4HENeutronInelastic();
>  //  theNeutronInelasticProcess->RegisterMe(theNeutronHPInelasticModel);
>    theNeutronInelasticProcess->RegisterMe(theNeutronBCModel);
>    theNeutronInelasticProcess->RegisterMe(theNeutronLEPModel);
>    theNeutronInelasticProcess->RegisterMe(theNeutronHEPModel);
>    pManager->AddDiscreteProcess(theNeutronInelasticProcess);
> 
>    G4HadronFissionProcess* thenFission
>                          = new G4HadronFissionProcess();
>    thenFission->AddDataSet(new G4NeutronHPFissionData);  
> 
>    G4LFission* thenFissionModel = new G4LFission();
>  //  thenFissionModel->SetMinEnergy(20.*MeV);
> //   G4NeutronHPFission* theNeutronHPFissionModel= new G4NeutronHPFission();
>  //  theNeutronHPFissionModel ->SetMaxEnergy(20.*MeV);
>    thenFission->RegisterMe(thenFissionModel);
> //   thenFission->RegisterMe(theNeutronHPFissionModel);
>    pManager->AddDiscreteProcess(thenFission);
> 
>    G4HadronCaptureProcess* thenCapture
>                          = new G4HadronCaptureProcess();
>    thenCapture->AddDataSet(new G4NeutronHPCaptureData);
> 
>    G4LCapture* thenCaptureModel = new G4LCapture();
>  //  thenCaptureModel->SetMinEnergy(20.*MeV);
>  //  G4NeutronHPCapture* theNeutronHPCaptureModel= new G4NeutronHPCapture();
>  //  theNeutronHPCaptureModel ->SetMaxEnergy(20.*MeV);
>    thenCapture->RegisterMe(thenCaptureModel);
>  //  thenCapture->RegisterMe(theNeutronHPCaptureModel);
>    pManager->AddDiscreteProcess(thenCapture);
> 
> There seems to be the same problem! Maybe there is something wrong with
> my cross section data? (I did set the environment variable
> NeutronHPCrossSections=/home/studeny/geant4/data/G4NDL3.7 )
> 
> *********************************************************************************************************
> * 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     1.03     -300     0.72         1        0        0         0       World initStep
>     1     1.03       -5     0.72         1        0      295       295 PhantomSens Transportation
>     2     1.03       -5     0.72         0        0 9.15e-19       295 PhantomSens NeutronInelastic
>     :----- List of 2ndaries - #SpawnInStep=  3(Rest= 0,Along= 0,Post= 3), #SpawnTotal=  3 ---------------
>     :      1.03        -5      0.72     0.722            neutron
>     :      1.03        -5      0.72     0.246              gamma
>     :      1.03        -5      0.72    0.0312          Cl35[0.0]
>     :----------------------------------------------------------------- EndOf2ndaries Info ---------------
> 
> When I change the material from G4_SKIN_ICRP to G4_WATER, there are also
> other reaktions, like :
> 
> *********************************************************************************************************
> * 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    -1.62     -300     1.93         1        0        0         0       World initStep
>     1    -1.62       -5     1.93         1        0      295       295 PhantomSens Transportation
>     2    -1.62       -4     1.93         1        0        1       296 PhantomSens Transportation
>     3    -1.62       -3     1.93         1        0        1       297 PhantomSens Transportation
>     4    -1.62       -2     1.93         1        0        1       298 PhantomSens Transportation
>     5    -1.62       -1     1.93         1        0        1       299 PhantomSens Transportation
>     6    -1.62        0     1.93         1        0        1       300 PhantomSens Transportation
>     7    -1.62        1     1.93         1        0        1       301 PhantomSens Transportation
>     8    -1.62        2     1.93         1        0        1       302 PhantomSens Transportation
>     9    -1.62        3     1.93         1        0        1       303 PhantomSens Transportation
>    10    -1.62        4     1.93         1        0        1       304 PhantomSens Transportation
>    11    -1.62        5     1.93         1        0        1       305       World Transportation
>    12    -1.62      500     1.93         1        0      495       800  OutOfWorld Transportation
> >>> Event 0
> 
> what can I do? Best regards Sylvia
> 

This G4_SCIN_ICRP has 13 components including even Fe. In water we have
only H and O.  Mass fractions of more heavy elements are small. However,
if there is a problem in one of inelastic cross section, than interaction
may happen very frequently.

What you can try is to increase verbosity level and to see mfp for water 
and skin - this will demonstrate how different are total cross sections.

VI

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

1 None: Re: neutron_HP_Physics and NIST_Material   (Sylvia Studeny - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Sylvia Studeny - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Vladimir Ivantchenko - 08 May, 2006)
(_ None: Re: neutron_HP_Physics and NIST_Material   (Sylvia Studeny - 09 May, 2006)
 Add Message Add Message
to: "Re: neutron_HP_Physics and NIST_Material"

 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 ]