Message: He3 not interacting Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None He3 not interacting 

Forum: Hadronic Processes
Date: 27 Oct, 2015
From: Lech Wiktor Piotrowski <Lech Wiktor Piotrowski>

Hello,

I am a Geant4 newbie. I am trying to make a simple He3 interaction in gas (modifying N01 example). I am shooting 1MeV He3 through 25 m of argon and it only transports through the gas and the world. The same is if I replace argon with Pb, so it is more or less obvious that it does not interact. I tried to follow the simple way to add physics:

ExN01PhysicsList::ExN01PhysicsList():G4VModularPhysicsList() {

	SetVerboseLevel(1);

	RegisterPhysics(new G4HadronElasticPhysics());
	RegisterPhysics(new G4IonElasticPhysics());
	RegisterPhysics(new G4EmStandardPhysics());
	RegisterPhysics(new G4IonPhysics());

    RegisterPhysics( new G4HadronInelasticQBBC(1));	
    RegisterPhysics( new G4DecayPhysics());	    
G4cerr << "constr" << G4endl;
}

void ExN01PhysicsList::ConstructParticle()
{
  // In this method, static member functions should be called
  // for all particles which you want to use.
  // This ensures that objects of these particle types will be
  // created in the program. 

  //G4Geantino::GeantinoDefinition();
  G4He3::He3Definition();
  G4GenericIon::Definition();
}

void ExN01PhysicsList::ConstructProcess()
{
  // Define transportation process

G4cerr << "regr" << G4endl;
  AddTransportation();
}

The constructor and ConstructProcess() are both called, but only transportation seems to take place. Why is the He3 not obeying the physics defined in the contructor?

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

1 None: Re: He3 not interacting   (Lech Wiktor Piotrowski - 27 Oct, 2015)
 Add Message Add Message
to: "He3 not interacting"

 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 ]