MonopolePhysicsTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// class header
#include "MonopolePhysicsTool.h"
 
#include <memory>
// package headers
#include "CustomMonopole.h"
#include "CustomMonopoleFactory.h"
#include "G4mplAtlasTransportation.h"
#include "G4mplAtlasIonisation.hh"
// Geant4 headers
#include "G4Version.hh"
#include "G4Transportation.hh"
#include "G4MuIonisation.hh"
#include "G4ProcessManager.hh"
#include "G4hIonisation.hh"
 
//-----------------------------------------------------------------------------
// Implementation file for class : MonopolePhysicsTool
//
// 14-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
//=============================================================================
MonopolePhysicsTool::MonopolePhysicsTool( const std::string& type,
                                          const std::string& nam,const IInterface* parent )
  : base_class ( type, nam , parent )
{
  m_physicsOptionType = G4AtlasPhysicsOption::Type::BSMPhysics;
}
 
//=============================================================================
// Destructor
//=============================================================================
 
MonopolePhysicsTool::~MonopolePhysicsTool()
{
}
 
//=============================================================================
// Initialize
//=============================================================================
StatusCode MonopolePhysicsTool::initialize( )
{
  ATH_MSG_DEBUG("MonopolePhysicsTool initialize( )");
  return StatusCode::SUCCESS;
}
 
auto MonopolePhysicsTool::GetPhysicsOption() -> UPPhysicsConstructor {
  return std::make_unique<MonopolePhysicsTool::PhysicsConstructor>(name(),
                                                                   msgLevel());
}
 
void MonopolePhysicsTool::PhysicsConstructor::ConstructParticle() {
  ATH_MSG_DEBUG(" ConstructParticle for the Monopole being run");
  CustomMonopoleFactory::instance();
}
void MonopolePhysicsTool::PhysicsConstructor::ConstructProcess() {
 
  PARTICLEITERATOR->reset();
 
  while((*PARTICLEITERATOR)())
    {
 
      CustomMonopole* particle = dynamic_cast<CustomMonopole*>(PARTICLEITERATOR->value());
      if(CustomMonopoleFactory::instance().isCustomMonopole(particle))
        {
          ATH_MSG_DEBUG( particle->GetParticleName() << " is Custom"  );
 
          G4ProcessManager* pmanager = particle->GetProcessManager();
 
          G4double magnCharge = particle->MagneticCharge();
 
 
          pmanager->RemoveProcess(0);
          pmanager->AddProcess(new G4mplAtlasTransportation(particle), -1, 0, 0);
 
          if (magnCharge != 0.0){
            pmanager->AddProcess(new G4mplAtlasIonisation(magnCharge, G4String("mplAtlasIonisation")), -1, 1, 1);
 
          }
 
 
          if(particle->GetPDGCharge() != 0.0 && magnCharge == 0.0) { // don't apply process for dyons, electric ionisation is applied in G4mplAtlasIonisationWithDeltaModel.cxx 
            pmanager->AddProcess(new G4hIonisation(),  -1, 2, 2);
          }
          pmanager->DumpInfo();
 
 
        }
    }
 
  ATH_MSG_DEBUG("ConstructProcess for Monopole finished");
}