#include <G4QInelasticCHIPS_HPBuilder.hh>

Public Member Functions
	G4QInelasticCHIPS_HPBuilder (G4int verbose)

virtual	~G4QInelasticCHIPS_HPBuilder ()

void	Build ()

Protected Attributes
G4ParticleTable *	theParticleTable

G4ParticleTable::G4PTblDicIterator *	theParticleIterator

Detailed Description

Definition at line 72 of file G4QInelasticCHIPS_HPBuilder.hh.

Constructor & Destructor Documentation

◆ G4QInelasticCHIPS_HPBuilder()

G4QInelasticCHIPS_HPBuilder::G4QInelasticCHIPS_HPBuilder ( G4int verbose )

Definition at line 51 of file G4QInelasticCHIPS_HPBuilder.cc.

                                                                 :
    verbose(ver)
    , wasActivated(false)
    , inelastic(0)
    , theInProcessMixer(0)
    , theNgProcessMixer(0)
    , theFiProcessMixer(0) 
    , theNeutronInelastic(0)
    , theNeutronFission(0)
    , theNeutronCapture(0)
    , theCHIPSInelastic(0)
    , theCHIPSNGamma(0)
    , theHPNeutron(0)
{
  // pointer to the particle table
  theParticleTable = G4ParticleTable::GetParticleTable();
  theParticleIterator = theParticleTable->GetIterator();
}

◆ ~G4QInelasticCHIPS_HPBuilder()

G4QInelasticCHIPS_HPBuilder::~G4QInelasticCHIPS_HPBuilder ( )

virtual

Definition at line 70 of file G4QInelasticCHIPS_HPBuilder.cc.

{
  if(wasActivated)
  {
    delete inelastic;
    delete theCHIPSInelastic;
    delete theCHIPSNGamma;
    //delete theCHIPSFission;
    delete theNeutronInelastic;
    delete theNeutronCapture;
    delete theNeutronFission;
    delete theHPNeutron;
    delete theInProcessMixer;
    delete theNgProcessMixer;
    delete theFiProcessMixer;
  }
}

Member Function Documentation

◆ Build()

void G4QInelasticCHIPS_HPBuilder::Build ( )

Definition at line 88 of file G4QInelasticCHIPS_HPBuilder.cc.

{
  if(wasActivated) return;
  wasActivated = true;
  theParticleIterator->reset();
  inelastic = new G4QInelastic();
  while( (*theParticleIterator)() )
  {
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4String pname = particle->GetParticleName();
    if(pname == "kaon-" || pname == "kaon+" || pname == "kaon0S"  ||  pname == "kaon0L" ||
       //pname == "pi-" || pname == "pi+"   || pname == "neutron" ||  pname == "proton" ||
       pname == "pi-"   || pname == "pi+"   ||  pname == "proton" ||
       pname == "lambda"       || pname == "sigma+"       || pname == "sigma0"       ||
       pname == "sigma-"       || pname == "xi0" || pname == "xi-" || pname == "omega-" ||
       pname == "anti_proton"  || pname == "anti_neutron" || pname == "anti_lambda"  ||
       pname == "anti_sigma+"  || pname == "anti_sigma0"  || pname == "anti_sigma-"  ||
       pname == "anti_xi0"     || pname == "anti_xi-"     || pname == "anti_omega-"  )
    {
      if(verbose>1)
        G4cout<< "__G4QInelCHIPS_HPBuilder: "<< pname <<" is defined here"<<G4endl;
      G4ProcessManager* pmanager = particle->GetProcessManager();
      pmanager->AddDiscreteProcess(inelastic);
      if(verbose>1)G4cout<<"###>G4QInelasticCHIPS_HPBuilder: "<<inelastic->GetProcessName()
                         <<" is added for "<<pname<<G4endl;
    }
    else if(pname == "neutron")
    {
      if(verbose>1)
        G4cout<< "__G4QInelCHIPS_HPBuilder: "<< pname <<" is defined here"<<G4endl;
      G4ProcessManager* pmanager = particle->GetProcessManager();
      // The model definition for neutrons (needed for HP implementation)
#ifdef debug
    G4cout<<"G4QInelasticCHIPS_HPBuilder::Build: before NeutronBuild"<<G4endl;
#endif
      theCHIPSInelastic   = new G4QInelastic();
      theCHIPSNGamma      = new G4QNGamma();
      //theCHIPSFission   = new G4QFission();
      theNeutronInelastic = new G4NeutronInelasticProcess();
      theNeutronCapture   = new G4HadronCaptureProcess();
      theNeutronFission   = new G4HadronFissionProcess();
      theInProcessMixer   = new G4QDiscProcessMixer("Mixed NeutronInelastic", particle);
      theNgProcessMixer   = new G4QDiscProcessMixer("Mixed NGamma", particle);
      theFiProcessMixer   = new G4QDiscProcessMixer("Mixed NFission", particle);
#ifdef debug
    G4cout<<"G4QInelasticCHIPS_HPBuilder::Build: before Build HP processes"<<G4endl;
#endif
      theHPNeutron = new G4NeutronHPBuilder;
      theHPNeutron->Build(theNeutronInelastic);
      theHPNeutron->Build(theNeutronCapture);
      theHPNeutron->Build(theNeutronFission);
#ifdef debug
    G4cout<<"G4QInelasticCHIPS_HPBuilder::Build: before QIn="<<theCHIPSInelastic<<G4endl;
#endif
      theInProcessMixer->AddDiscreteProcess(theCHIPSInelastic, 1.E8*megaelectronvolt);
#ifdef debug
    G4cout<<"G4QInelasticCHIPS_HPBuilder::Build: befr HPI="<<theNeutronInelastic<<G4endl;
#endif
      theInProcessMixer->AddDiscreteProcess(theNeutronInelastic, 19.9*megaelectronvolt);
      
#ifdef debug
    G4cout<<"G4QInelasticCHIPS_HPBuilder::Build: before QNG="<<theCHIPSNGamma<<G4endl;
#endif
      theNgProcessMixer->AddDiscreteProcess(theCHIPSNGamma, 1.E8*megaelectronvolt);
#ifdef debug
    G4cout<<"G4QInelasticCHIPS_HPBuilder::Build: before HPC="<<theNeutronCapture<<G4endl;
#endif
      theNgProcessMixer->AddDiscreteProcess(theNeutronCapture, 19.9*megaelectronvolt);
      
      //theFiProcessMixer->AddDiscreteProcess(theCHIPSFission, 1.E8*megaelectronvolt);
      //theFiProcessMixer->AddDiscreteProcess(theNeutronFission, 19.9*megaelectronvolt);
      
#ifdef debug
    G4cout<<"G4QInelasticCHIPS_HPBuilder::Build: before ProcessAdd"<<G4endl;
#endif
      pmanager->AddDiscreteProcess(theInProcessMixer); // Mix CHIPS+HP for neutronInelastic
      if(verbose>1)
        G4cout<<"###>G4QInelasticCHIPS_HPBuilder: "<<theCHIPSInelastic->GetProcessName()
              <<" is added for "<<pname<<G4endl;
      pmanager->AddDiscreteProcess(theNgProcessMixer);  // Mix CHIPS+HP for (n,gamma)
      if(verbose>1)
        G4cout<<"###>G4QInelasticCHIPS_HPBuilder: "<<theCHIPSNGamma->GetProcessName()
              <<" is added for "<<pname<<G4endl;
      pmanager->AddDiscreteProcess(theNeutronFission); // Only HP for fission
      //pmanager->AddDiscreteProcess(theFiProcessMixer); // Mix CHIPS+HP for fission
      if(verbose>1)
        G4cout<<"###>G4QInelasticCHIPS_HPBuilder: "<<theNeutronFission->GetProcessName()
              <<" is added for "<<pname<<G4endl;
    }
  }
}