Message: 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 pi minus stoping energy deposit at 150 MeV/c  

Keywords: pion stopping energy deposit
Forum: Physics List
Date: 04 Jan, 2004
From: Ran Shneor <Ran Shneor>

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

 

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

1 None: Re: pi minus stoping energy deposit at 150 MeV/c   (Vladimir IVANTCHENKO - 04 Jan, 2004)
 Add Message Add Message
to: "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 ]