GauginosPhysicsTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// class header
#include "GauginosPhysicsTool.h"
// package headers
#include "GMSBNeutralino.hh"
#include "GMSBGravitino.hh"
// Geant4 physics lists
#include "G4ProcessManager.hh"
#include "G4ParticleTable.hh"
#include "G4hIonisation.hh"
#include "G4hMultipleScattering.hh"
#include "G4ParticleDefinition.hh"
#include "G4Transportation.hh"
#include "G4MuIonisation.hh"
// CLHEP headers
#include "CLHEP/Units/SystemOfUnits.h"
 
//-----------------------------------------------------------------------------
// Implementation file for class : GauginosPhysicsTool
//
// 15-05-2015 : Edoardo Farina
//-----------------------------------------------------------------------------
 
//=============================================================================
// Standard constructor, initializes variables
//=============================================================================
GauginosPhysicsTool::GauginosPhysicsTool( const std::string& type,
                                          const std::string& nam,const IInterface* parent )
  : base_class ( type, nam , parent )
{
  m_physicsOptionType = G4AtlasPhysicsOption::Type::BSMPhysics;
 
  declareProperty("GravitinoMass",m_GravitinoMass=0.108E-04*CLHEP::GeV,"Gravitino Mass");
  declareProperty("GravitinoWidth",m_GravitinoWidth=0.*CLHEP::GeV,"Gravitino Width");
  declareProperty("GravitinoCharge",m_GravitinoCharge=0,"Gravitino charge");
  declareProperty("GravitinoPDGCode",m_GravitinoPDGCode=1000039,"Gravitino PDG CODE");
  declareProperty("GravitinoLifetime",m_GravitinoLifetime=-1 ,"Gravitino Lifetime");
  declareProperty("GravitinoStable",m_GravitinoStable=true ,"Gravitino Stable");
  declareProperty("GravitinoShortlived",m_GravitinoShortlived=false ,"Gravitino Shortlived");
 
  declareProperty("NeutralinoMass",m_NeutralinoMass=118.848*CLHEP::GeV,"Neutralino Mass");
  declareProperty("NeutralinoWidth",m_NeutralinoWidth=0.*CLHEP::GeV,"Neutralino Width");
  declareProperty("NeutralinoCharge",m_NeutralinoCharge=0,"Neutralino charge");
  declareProperty("NeutralinoPDGCode",m_NeutralinoPDGCode=1000022,"Neutralino PDG CODE");
  declareProperty("NeutralinoStable",m_NeutralinoStable=true,"Neutralino Stable");
  declareProperty("NeutralinoLifetime",m_NeutralinoLifetime=-1,"Neutralino Lifetime");
  declareProperty("NeutralinoShortlived",m_NeutralinoShortlived=false,"Neutralino Shortlived");
 
}
 
//=============================================================================
// Destructor
//=============================================================================
 
GauginosPhysicsTool::~GauginosPhysicsTool()
{
}
 
//=============================================================================
// Initialize
//=============================================================================
StatusCode GauginosPhysicsTool::initialize( )
{
  ATH_MSG_DEBUG("GauginosPhysicsTool initialize(  )");
 
  this->SetPhysicsName(name());
  return StatusCode::SUCCESS;
}
 
 
GauginosPhysicsTool* GauginosPhysicsTool::GetPhysicsOption()
{
  return this;
}
 
void GauginosPhysicsTool::ConstructParticle()
{
  ATH_MSG_DEBUG("Create particle of Gauginos" );
 
  GMSBNeutralino::Definition(m_NeutralinoMass, m_NeutralinoWidth, m_NeutralinoCharge, m_NeutralinoPDGCode, m_NeutralinoStable, m_NeutralinoLifetime, m_NeutralinoShortlived );
 
  GMSBGravitino::Definition(m_GravitinoMass, m_GravitinoWidth, m_GravitinoCharge, m_GravitinoPDGCode, m_GravitinoStable, m_GravitinoLifetime, m_NeutralinoShortlived);
 
}
 
void GauginosPhysicsTool::ConstructProcess()
{
  ATH_MSG_DEBUG(" Construct Process for the Gauginos being run");
}