Message: Regarding Proton and Triton Scattering Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Regarding Proton and Triton Scattering 

Forum: Hadronic Processes
Date: 30 Jul, 2007
From: Sarah Xu <Sarah Xu>

Hi All,

I'm trying to see how the protons and tritions produced in the 3He(n, p)3H are dissipating thier energy inside the helium3 chamber. In order to do that, I first tried to simply shoot a triton/proton from the surface into the 3He chamber with energies in hundred keV range. The following is the code for my physics list:

#include "NcdPhysicsList.hh"
#include "G4ParticleTypes.hh"
#include "G4ProcessManager.hh"
#include "G4ParticleDefinition.hh"
#include "G4ParticleTable.hh"

NcdPhysicsList::NcdPhysicsList()
{}
NcdPhysicsList::~NcdPhysicsList()
{;}
//////////////////////////////////////////
//Construct Particles
void NcdPhysicsList::ConstructParticle()
{
  ConstructHadrons();
  ConstructBosons();
  ConstructBaryons();
  ConstructLeptons();
  ConstructMesons();
}
void NcdPhysicsList::ConstructHadrons(){ 
...

//////////////////////////////////////////// //Construct process

#include "G4HadronElasticProcess.hh"
#include "G4TritonInelasticProcess.hh"
#include "G4ProtonInelasticProcess.hh"

#include "G4LETritonInelastic.hh"
#include "G4LEProtonInelastic.hh"
#include "G4LElastic.hh"
void NcdPhysicsList::ConstructProcess()
{
   AddTransportation();
   ConstructHaProc();
}

void NcdPhysicsList::ConstructHaProc()
{
  G4ParticleDefinition* triton = G4Triton::TritonDefinition();
  G4ProcessManager* pManager = triton->GetProcessManager();
  //
  G4HadronElasticProcess* hadElastProc = new G4HadronElasticProcess();
  G4LElastic* elasticModel = new G4LElastic();
  hadElastProc->RegisterMe(elasticModel);
  pManager->AddDiscreteProcess(hadElastProc);
  //
      G4TritonInelasticProcess* tinelProc = new G4TritonInelasticProcess();
      G4LETritonInelastic* TriModel = new G4LETritonInelastic();
      tinelProc->RegisterMe(TriModel);
      pManager->AddDiscreteProcess(tinelProc); 
//
  G4ParticleDefinition* proton = G4Proton::ProtonDefinition();
  G4ProcessManager* pmanager = proton->GetProcessManager();
  pmanager->AddDiscreteProcess(hadElastProc);
  //
   G4ProtonInelasticProcess* proProc = new G4ProtonInelasticProcess();
   G4LEProtonInelastic* proModel = new G4LEProtonInelastic();
   proProc->RegisterMe(proModel);
   pmanager->AddDiscreteProcess(proProc);
}

void NcdPhysicsList::SetCuts()
{
  SetCutsWithDefault();
}

And after I tried to run the process, only transportation process is involved, for instance:

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

Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0      0 fm  -2.54 cm      0 fm    573 keV     0 eV      0 fm      0 fm   Helium_Tub_physical    initStep
    1  -2.76 mm  -2.52 cm   1.94 mm    573 keV     0 eV   3.37 mm   3.37 mm   Helium_Tub_physical  Transportation
    2    -20 cm  -1.45 cm   14.1 cm    573 keV     0 eV   24.1 cm   24.5 cm   OutOfWorld  Transportation
Track (trackID 1, parentID 0) is processed with stopping code 2

I'm wondering if anyone could shine some light on me. Any help is highly appreciated.

Best Sarah

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

1 None: Re: Regarding Proton and Triton Scattering   (Vladimir IVANTCHENKO - 30 Jul, 2007)
2 None: Re: Regarding Proton and Triton Scattering   (Sarah Xu - 30 Jul, 2007)
1 None: Re: Regarding Proton and Triton Scattering   (Vladimir IVANTCHENKO - 31 Jul, 2007)
 Add Message Add Message
to: "Regarding Proton and Triton Scattering"

 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 ]