TG4ExtDecayerPhysics.cxx

Go to the documentation of this file.

//------------------------------------------------
// The Geant4 Virtual Monte Carlo package
// Copyright (C) 2007 - 2014 Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see geant4_vmc/LICENSE.
// Contact: root-vmc@cern.ch
//-------------------------------------------------
 
 
#include "TG4ExtDecayerPhysics.h"
#include "TG4ExtDecayer.h"
#include "TG4Globals.h"
 
#include <TVirtualMC.h>
#include <TVirtualMCDecayer.h>
 
#include <G4Decay.hh>
#include <G4ParticleDefinition.hh>
#include <G4ProcessManager.hh>
 
//_____________________________________________________________________________
TG4ExtDecayerPhysics::TG4ExtDecayerPhysics(const G4String& name)
  : TG4VPhysicsConstructor(name),
    fMessenger(this),
    fDecayProcess(0),
    fSelection(),
    fSkipNeutrino(false)
{
}
 
//_____________________________________________________________________________
TG4ExtDecayerPhysics::TG4ExtDecayerPhysics(
  G4int theVerboseLevel, const G4String& name)
  : TG4VPhysicsConstructor(name, theVerboseLevel),
    fMessenger(this),
    fDecayProcess(0),
    fSelection(),
    fSkipNeutrino(false)
{
}
 
//_____________________________________________________________________________
TG4ExtDecayerPhysics::~TG4ExtDecayerPhysics()
{
 
  delete fDecayProcess;
}
 
//
// protected methods
//
 
//_____________________________________________________________________________
void TG4ExtDecayerPhysics::ConstructParticle()
{
}
 
//_____________________________________________________________________________
void TG4ExtDecayerPhysics::ConstructProcess()
{
 
  // Check if VMC decayer is defined
  TVirtualMCDecayer* mcDecayer = gMC->GetDecayer();
  if (!mcDecayer) {
    // TG4Globals::Warning(
    //  "TG4ExtDecayerPhysics", "ConstructProcess",
    //  "No VMC external decayer defined.");
    return;
  }
 
  // Create Geant4 external decayer
  TG4ExtDecayer* tg4Decayer = new TG4ExtDecayer(mcDecayer);
  tg4Decayer->VerboseLevel(VerboseLevel());
  tg4Decayer->SetSkipNeutrino(fSkipNeutrino);
  // The tg4Decayer is deleted in G4Decay destructor
  // But we may have a problem if there are more than one
  // instances of G4Decay process
 
  auto aParticleIterator = GetParticleIterator();
  aParticleIterator->reset();
  while ((*aParticleIterator)()) {
    G4ParticleDefinition* particle = aParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
 
    // skip particles which do not have process manager
    if (!pmanager) continue;
 
    if (fSelection.find(particle->GetParticleName()) != std::string::npos) {
 
      if (VerboseLevel() > 1) {
        G4cout << "Switching off Geant4 decay table for: "
               << particle->GetParticleName() << G4endl;
      }
 
      // Unset the decay table for particles in a selection;
      // for the particles in selection, the external decayer
      // will have a priority over Geant4 decay table
      particle->SetDecayTable(0);
    }
 
    if (VerboseLevel() > 1) {
      G4cout << "Setting ext decayer for: "
             << aParticleIterator->value()->GetParticleName() << G4endl;
    }
 
    G4ProcessVector* processVector = pmanager->GetProcessList();
    for (size_t i = 0; i < processVector->length(); i++) {
 
      G4Decay* decay = dynamic_cast<G4Decay*>((*processVector)[i]);
      if (decay) decay->SetExtDecayer(tg4Decayer);
    }
  }
 
  if (VerboseLevel() > 0) {
    G4cout << "### "
           << "External decayer physics constructed." << G4endl;
  }
}