SimG4Core/CustomPhysics/src/CMSAPrime.cc - master - githubmirror/cmssw

SimG4Core/CustomPhysics/src/CMSAPrime.cc
80 строк3 KB
Vladimir
code-format
06 июл 2025, 15:06
06 июл 2025, 15:064a85dd5
#include "SimG4Core/CustomPhysics/interface/CMSAPrime.h"
#include <CLHEP/Units/SystemOfUnits.h>

CMSAPrime* CMSAPrime::theAPrime = nullptr;

CMSAPrime::CMSAPrime(const G4String& aName,
                     G4double mass,
                     G4double width,
                     G4double charge,
                     G4int iSpin,
                     G4int iParity,
                     G4int iConjugation,
                     G4int iIsospin,
                     G4int iIsospin3,
                     G4int gParity,
                     const G4String& pType,
                     G4int lepton,
                     G4int baryon,
                     G4int encoding,
                     G4bool stable,
                     G4double lifetime,
                     G4DecayTable* decaytable)
    : G4ParticleDefinition(aName,
                           mass,
                           width,
                           charge,
                           iSpin,
                           iParity,
                           iConjugation,
                           iIsospin,
                           iIsospin3,
                           gParity,
                           pType,
                           lepton,
                           baryon,
                           encoding,
                           stable,
                           lifetime,
                           decaytable) {}

CMSAPrime* CMSAPrime::APrime(double apmass) {
  if (!theAPrime) {
    const G4String& name = "A^1";
    G4double mass = apmass * CLHEP::MeV;
    G4double width = 0.;
    G4double charge = 0;
    G4int iSpin = 0;
    G4int iParity = 0;
    G4int iConjugation = 0;
    G4int iIsospin = 0;
    G4int iIsospin3 = 0;
    G4int gParity = 0;
    const G4String& pType = "APrime";
    G4int lepton = 0;
    G4int baryon = 0;
    G4int encoding = 9994;
    G4bool stable = true;
    G4double lifetime = -1;
    G4DecayTable* decaytable = nullptr;

    theAPrime = new CMSAPrime(name,
                              mass,
                              width,
                              charge,
                              iSpin,
                              iParity,
                              iConjugation,
                              iIsospin,
                              iIsospin3,
                              gParity,
                              pType,
                              lepton,
                              baryon,
                              encoding,
                              stable,
                              lifetime,
                              decaytable);
  }
  return theAPrime;
}