Message: Help with neutron process. Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Help with neutron process. 

Forum: Hadronic Processes
Date: 27 Jun, 2012
From: Dat Vu <Dat Vu>

Hi everybody. I'm simulating the process: n + B -> Li + alpha particle + gamma ray. Here is my physics list:

       if (particleName == "neutron")
{
  // elastic scattering
  G4HadronElasticProcess* theNeutronElasticProcess =
    new G4HadronElasticProcess;
  G4LElastic* theElasticModel1 = new G4LElastic;
  theNeutronElasticProcess->RegisterMe(theElasticModel1);
  theNeutronElasticProcess->RegisterMe(new G4NeutronHPElastic());
  theNeutronElasticProcess->AddDataSet(new G4NeutronHPElasticData());
  pmanager->AddDiscreteProcess(theNeutronElasticProcess);

  // inelastic scatteringG4NeutronHPInelastic
  G4NeutronInelasticProcess* theInelasticProcess =
    new G4NeutronInelasticProcess();
  G4BinaryCascade* theBinaryModel = new G4BinaryCascade();
  theInelasticProcess->RegisterMe(theBinaryModel);
  theInelasticProcess->RegisterMe(new G4HENeutronInelastic());
  theInelasticProcess->RegisterMe(new G4NeutronHPInelastic());

  theInelasticProcess->AddDataSet(new G4NeutronHPInelasticData());
  pmanager->AddDiscreteProcess(theInelasticProcess);

  // capture process
  G4LCapture* theCaptureModel = new G4LCapture;
  theCaptureProcess->RegisterMe(theCaptureModel);
  theCaptureProcess->RegisterMe(new G4NeutronHPCapture());
  theCaptureProcess->AddDataSet(new G4NeutronHPCaptureData());

  pmanager->AddDiscreteProcess(theCaptureProcess);
}

and in main(), I included:

using namespace std;

const G4int USEQGSP_BERT_HP = 1; const G4int USEMYPHYSICS = 1;

int main(int argc,char** argv) {

  CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine);

  G4RunManager * runManager = new G4RunManager;

  runManager->SetUserInitialization(new B1DetectorConstruction);
  if (USEMYPHYSICS == 1)
  {
   B1PhysicsList* physics = new B1PhysicsList;
   runManager->SetUserInitialization(physics);

  }
  if (USEQGSP_BERT_HP == 1)
  {
  G4VModularPhysicsList* physicsList = new QGSP_BERT_HP;
  physicsList->SetVerboseLevel(1);
  runManager->SetUserInitialization(physicsList);
  }
}

However, when I shoot a neutron into the detector, the neutron only scatters. And nothing else happens. What can I do so that the neutron reacts with the detector giving the alpha particle and gamma ray?

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

1 None: Re: Help with neutron process.   (michel maire - 27 Jun, 2012)
1 None: Re: Help with neutron process.   (Dat Vu - 27 Jun, 2012)
2 None: Re: Help with neutron process.   (michel maire - 29 Jun, 2012)
 Add Message Add Message
to: "Help with neutron process."

 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 ]