Message: Spurious warnings after adding radioactive decay Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Spurious warnings after adding radioactive decay 

Forum: Physics List
Date: 18 Jun, 2015
From: Karen Pease <Karen Pease>

I presume I did something wrong but I don't have the experience to know what. I added in radioactive decay support based on looking at the example rdecay02; my current physics list was based on FP_BERT_HP, so I couldn't just copy it directly. The code is:

PhysicsList::PhysicsList(const G4int verbose) :
  FTFP_BERT_HP(verbose),
  fEmPhysicsList(0),      
  fRaddecayList(0),       
  fParticleList(0)        
{
  G4LossTableManager::Instance();
  defaultCutValue =1.*mm;

  SetVerboseLevel(verbose);

  fParticleList = new G4DecayPhysics();
  fRaddecayList = new G4RadioactiveDecayPhysics();
  fEmPhysicsList = new G4EmStandardPhysics();
}

PhysicsList::~PhysicsList()
{
  delete fParticleList;
  delete fRaddecayList;
  delete fEmPhysicsList;
}

void PhysicsList::ConstructParticle()
{
  fParticleList->ConstructParticle();
}

void PhysicsList::ConstructProcess()
{
  G4VModularPhysicsList::ConstructProcess(); 

  fEmPhysicsList->ConstructProcess();
  fParticleList->ConstructProcess(); 
  fRaddecayList->ConstructProcess(); 

  G4ParticleDefinition* nPD = G4Neutron::Definition();
  G4ProcessVector* pvec = nPD->GetProcessManager()->GetProcessList();
  G4HadronicProcess* nElastic = 0;
  for (G4int i=0; i<pvec->size(); i++) {
    if ((*pvec)[i]->GetProcessSubType() != fHadronElastic) continue;
    nElastic = dynamic_cast<G4HadronicProcess*>((*pvec)[i]);
    break;
  }
  assert(nElastic != 0);

  G4HadronicInteraction* nElasticHP = G4HadronicInteractionRegistry::Instance()->FindModel("NeutronHPElastic"); 
  assert(nElasticHP != 0);
  nElasticHP->SetMinEnergy(4.*eV);

  nElastic->RegisterMe(new G4NeutronHPThermalScattering);
  nElastic->AddDataSet(new G4NeutronHPThermalScatteringData);        
}       

void PhysicsList::SetCuts()
{
  SetCutsWithDefault();
  if (verboseLevel>0) DumpCutValuesTable();
}

The sort of warnings that I get flooded with (thousands of lines) are along the lines of:

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : Run0111
      issued by : G4PhysicsListHelper::RegisterPorcess
Duplication of processes
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

G4PhysicsListHelper::RegisterProcess :hIoni for B+ with type/subtype =2/2 is h as same subType as hIoni for B+ It will not be added !!

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : Run0111
      issued by : G4PhysicsListHelper::RegisterPorcess
Duplication of processes
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

G4PhysicsListHelper::RegisterProcess :msc for B- with type/subtype =2/10 is ha s same subType as msc for B- It will not be added !!

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : Run0111
      issued by : G4PhysicsListHelper::RegisterPorcess
Duplication of processes
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

G4PhysicsListHelper::RegisterProcess :hIoni for B- with type/subtype =2/2 is h as same subType as hIoni for B- It will not be added !!

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : Run0111
      issued by : G4PhysicsListHelper::RegisterPorcess
Duplication of processes
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

G4PhysicsListHelper::RegisterProcess :msc for D+ with type/subtype =2/10 is ha s same subType as msc for D+ It will not be added !!

Etc. I presume that this means that RegisterProcess (Porcess?) is being called twice but I can't see why that would be happening.

Either way, radioactive decay seems to be occurring, so that's good.

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

1 Question: Re: Spurious warnings after adding radioactive decay   (Andreas Gad - 17 Aug, 2016)
(_ None: Re: Spurious warnings after adding radioactive decay   (Karen Pease - 17 Aug, 2016)
 Add Message Add Message
to: "Spurious warnings after adding radioactive decay"

 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 ]