TRTPhysicsTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#include "TRTPhysicsTool.h"
#include "TRTTransitionRadiation.h"
 
#include "G4ProcessManager.hh"
#include "G4ParticleTable.hh"
#include "G4VProcess.hh"
#include "G4Version.hh"
#include "G4AutoDelete.hh"
 
#include "globals.hh"
#include <iostream>
#include <memory>
 
//-----------------------------------------------------------------------------
// Implementation file for class : TRTPhysicsTool
//
// 18-05-2015 Edoardo Farina
//-----------------------------------------------------------------------------
 
#if G4VERSION_NUMBER > 1029
#define PARTICLEITERATOR (this->GetParticleIterator())
#elif G4VERSION_NUMBER > 1009
#define PARTICLEITERATOR aParticleIterator
#else
#define PARTICLEITERATOR theParticleIterator
#endif
 
 
//=============================================================================
// Standard constructor, initializes variables
//=============================================================================
TRTPhysicsTool::TRTPhysicsTool(const std::string& type, const std::string& name,
                               const IInterface* parent)
    : base_class(type, name, parent) {
  m_physicsOptionType = G4AtlasPhysicsOption::Type::GlobalProcesses;
 
  declareProperty("XMLFile", m_xmlFile="TRgeomodelgeometry.xml");
}
 
//=============================================================================
// Initialize
//=============================================================================
StatusCode TRTPhysicsTool::initialize( )
{
  ATH_MSG_DEBUG("TRTPhysicsTool initialize()");
  return StatusCode::SUCCESS;
}
 
//=============================================================================
// Return the physics constructor
//=============================================================================
auto TRTPhysicsTool::GetPhysicsOption() -> UPPhysicsConstructor {
  return std::make_unique<TRTPhysicsTool::PhysicsConstructor>(name(), this->msgLevel(), m_xmlFile);
}
 
//=============================================================================
// Particle construction; not implemented
//=============================================================================
void TRTPhysicsTool::PhysicsConstructor::ConstructParticle() {}
 
//=============================================================================
// Physics process construction
//=============================================================================
void TRTPhysicsTool::PhysicsConstructor::ConstructProcess() {
  ATH_MSG_DEBUG("TRTPhysicsTool::ConstructProcess() - start");
 
  // Use the Geant4 garbage collection mechanism to clean this up.
  // It's just a convenient way to clean up in multi-threading.
  auto trProc = new TRTTransitionRadiation("XTR", m_xmlFile) ;
  G4AutoDelete::Register(trProc);
 
  PARTICLEITERATOR->reset();
  while( (*PARTICLEITERATOR)() ) {
    G4ParticleDefinition* particle = PARTICLEITERATOR->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    if ( particle->GetPDGCharge() != 0.0 && particle->GetPDGMass() != 0.0 ) {
      trProc->SetVerboseLevel(1);
      ATH_MSG_DEBUG("TRT Process Added to "<<particle->GetParticleName());
      pmanager->AddDiscreteProcess(trProc);
    }
  }
  ATH_MSG_DEBUG("TRTPhysicsTool::ConstructProcess() - end");
}