Message: Re: pi minus stoping energy deposit at 150 MeV/c Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: pi minus stoping energy deposit at 150 MeV/c 

Forum: Physics List
Re: None pi minus stoping energy deposit at 150 MeV/c (Ran Shneor)
Date: 04 Jan, 2004
From: Vladimir Ivanchenko <Vladimir Ivanchenko>

On Sun, 4 Jan 2004, Ran Shneor wrote:

> *** Discussion title: Physics List
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/phys-list/42"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> Hi
> 
> I am trying to simulate the energy deposit by pi- hitting an infinit
>  scintillator at 150 MeV/c. I expect to get a total energy deposit 
> of ~200 MeV which is the result of the kinetic energy and the 
> pi-  mass. (wich in theory should be all deposited in the scintillator 
> matirial). At the moment I am getting an energy deposit which is equal to the 
> kinetic energy (~150 MeV).
> 
> this is the physics list I am using -
> 
> #ifndef BBPL_VAR
> #define BBPL_VAR 1
> 
> #include "BB6PhysicsList.hh"
> #include "globals.hh"
> #include "G4ProcessManager.hh"
> #include "G4ParticleTypes.hh"
> #include "G4MultipleScattering.hh"
> 
> #include "G4hIonisation.hh"
> 
> #include "G4LCapture.hh"
> #include "G4HadronElasticProcess.hh"
> #include "G4NeutronInelasticProcess.hh"
> 
> #include "G4ProtonInelasticProcess.hh"
> #include "G4PionMinusInelasticProcess.hh"
> #include "G4PiMinusAbsorptionAtRest.hh"
> #include "G4PiMinusStopC.hh"
> #include "G4PionMinusAbsorptionAtRest.hh"
> 
> BB2PhysicsList::BB2PhysicsList() {;}
> 
> BB2PhysicsList::~BB2PhysicsList() {;}
> 
> void BB2PhysicsList::ConstructParticle() {
> 
>   G4Neutron::NeutronDefinition();
>   G4Proton::ProtonDefinition();
>   G4PionMinus::PionMinusDefinition();
>   G4PionPlus::PionPlusDefinition();
>   G4PionZero::PionZeroDefinition();
> }
> 
> void BB2PhysicsList::ConstructProcess()
> {
>   // Define transportation process
>   AddTransportation();
> 
>   // Taken from example 2 of G4
>   G4ParticleTable* partable = G4ParticleTable::GetParticleTable();
>     G4ParticleDefinition* particle = partable->FindParticle("pi-");
> 
>   if (particle==NULL) { G4cout << "\n NULL !!" ; }
> 
>   G4String partname = particle->GetParticleName();
> 
>   G4PionMinusAbsorptionAtRest pmar;
>   G4PiMinusStopC pmsc;
>   G4PiMinusAbsorptionAtRest pmar1;
>   G4PionMinusInelasticProcess pmip;
> 
> // G4ProcessManager * pmanager =
> G4Neutron::Neutron()->GetProcessManager(); // G4ProcessManager *
> pmanager = G4Proton::Proton()->GetProcessManager();
> 
> G4ProcessManager * pmanager =
> G4PionMinus::PionMinus()->GetProcessManager();
> 
>  pmanager->AddDiscreteProcess((G4VProcess *) &pmar);
>  pmanager->AddDiscreteProcess((G4VProcess *) &pmar1);
>  pmanager->AddDiscreteProcess((G4VProcess *) &pmsc);
>  pmanager->AddDiscreteProcess((G4VProcess *) &pmip);
> 
> }
> 
> void BB2PhysicsList::SetCuts() //a faire ### {
> 
>   G4int temp = GetVerboseLevel();
> 
>   SetCutsWithDefault();
>   SetCutValue(1.*mm,"pi-");
> 
> SetVerboseLevel(temp); }
> 
> #endif
> 
> 100000 x thanks
> 
> Ran Shneor
> 
>  
> 

First of all, one should not include all possible processes to simulate
pi- stopping at rest - they are alternative. It is difficult to predit how 
such PhysicsList will work.

The best is to use one of PhysicsLists from the web page of G4 Hadronics.

VI

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

 Add Message Add Message
to: "Re: pi minus stoping energy deposit at 150 MeV/c"

 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 ]