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 )
Date: 08 May, 2006
From: Sylvia Studeny <sylvia.studeny@gsf.de>

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

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   (Vladimir Ivantchenko - 08 May, 2006)
(_ 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 ]