Singleton class to load & serve a TDatabasePDG. More...

#include <PDGLibrary.h>

Collaboration diagram for genie::PDGLibrary:

Classes
struct	Cleaner

Public Member Functions
TDatabasePDG *	DBase (void)

TParticlePDG *	Find (int pdgc, bool must_exist=true)

void	ReloadDBase (void)

void	AddDarkMatter (double mass, double med_ratio)

Static Public Member Functions
static PDGLibrary *	Instance (void)

Private Member Functions
	PDGLibrary ()

	PDGLibrary (const PDGLibrary &config_pool)

virtual	~PDGLibrary ()

bool	LoadDBase (void)

bool	AddDarkSector ()

bool	AddHNL (void)

Private Attributes
TDatabasePDG *	fDatabasePDG

Static Private Attributes
static PDGLibrary *	fInstance = 0

Friends
struct	Cleaner

Detailed Description

Singleton class to load & serve a TDatabasePDG.

Author: Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool

Changes required to implement the GENIE Boosted Dark Matter module were installed by Josh Berger (Univ. of Wisconsin)

Changes required to implement the GENIE Dark Neutrino module were installed by Iker de Icaza (Univ. of Sussex)

Changes required to implement the GENIE BeamHNL module were installed by John Plows (Univ. of Oxford)

Created:: May 06, 2004

License:: Copyright (c) 2003-2024, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 35 of file PDGLibrary.h.

Constructor & Destructor Documentation

PDGLibrary::PDGLibrary ( )

private

Definition at line 40 of file PDGLibrary.cxx.

References AddDarkSector(), AddHNL(), fInstance, LoadDBase(), LOG, pERROR, pFATAL, and pINFO.

Referenced by Instance().

 {
   if( ! LoadDBase() ) LOG("PDG", pERROR) << "Could not load PDG data";
 
 #ifdef __GENIE_DARK_NEUTRINO_ENABLED__
   LOG("PDG", pINFO) << "Loading Dark sector Info";
   if ( ! AddDarkSector() ) { 
     LOG("PDG", pFATAL) << "Could not load Dark Neutrino data";
     exit(78);
   }
 #endif // __GENIE_DARK_NEUTRINO_ENABLED__
 
 #ifdef __GENIE_HEAVY_NEUTRAL_LEPTON_ENABLED__
   LOG("PDG", pINFO) << "Loading Heavy Neutral Lepton data";
   if( ! AddHNL() ){
     LOG("PDG", pFATAL) << "Could not load Heavy Neutral Lepton data";
     exit(78);
   }
 #endif // #ifdef __GENIE_HEAVY_NEUTRAL_LEPTON_ENABLED__
   
   fInstance =  0;
 }

genie::PDGLibrary::PDGLibrary ( const PDGLibrary & config_pool )

private

PDGLibrary::~PDGLibrary ( )

privatevirtual

Definition at line 63 of file PDGLibrary.cxx.

References fInstance.

 {
   fInstance = 0;
 }

Member Function Documentation

void PDGLibrary::AddDarkMatter	(	double	mass,
		double	med_ratio
	)

Definition at line 153 of file PDGLibrary.cxx.

References fDatabasePDG, genie::kPdgDarkMatter, and genie::kPdgMediator.

Referenced by main().

 {
 // Add dark matter particle to PDG database
 
   double med_mass = mass*med_ratio;
   TParticlePDG * dm_particle = fDatabasePDG->GetParticle(kPdgDarkMatter);
   TParticlePDG * med_particle = fDatabasePDG->GetParticle(kPdgMediator);
   if (!dm_particle) {
     // Name Title Mass Stable Width Charge Class PDG
     fDatabasePDG->AddParticle("chi_dm","chi_dm",mass,true,0.,0,"DarkMatter",kPdgDarkMatter);
   }
   else {
     assert(dm_particle->Mass() == mass);
   }
   if (!med_particle) {
     // Name Title Mass Stable Width Charge Class PDG
     fDatabasePDG->AddParticle("Z_prime","Z_prime",med_mass,true,0.,0,"DarkMatter",kPdgMediator);
   }
   else {
     assert(med_particle->Mass() == med_mass);
   }
 }

bool PDGLibrary::AddDarkSector ( )

private

Definition at line 193 of file PDGLibrary.cxx.

References genie::AlgConfigPool::CommonList(), fDatabasePDG, genie::Registry::GetDouble(), genie::AlgConfigPool::Instance(), genie::kPdgAntiDarkNeutrino, genie::kPdgDarkNeutrino, genie::kPdgDNuMediator, LOG, and pERROR.

Referenced by PDGLibrary().

 {
   // Add dark neutrino particles to PDG database
 
   const Registry * reg = AlgConfigPool::Instance()->CommonList("Dark", "Masses");
   if(!reg) {
     LOG("PDG", pERROR) << "The Dark Sector masses not available.";
     return false;
   }
   TParticlePDG * dnu_particle = fDatabasePDG->GetParticle(kPdgDarkNeutrino);
   TParticlePDG * anti_dnu_particle = fDatabasePDG->GetParticle(kPdgAntiDarkNeutrino);
   TParticlePDG * med_particle = fDatabasePDG->GetParticle(kPdgDNuMediator);
   if (!dnu_particle) {
     // Name Title Mass Stable Width Charge Class PDG
     fDatabasePDG->AddParticle("nu_D","#nu_{D}",reg->GetDouble("Dark-NeutrinoMass"),
                               true,0.,0,"DarkNeutrino",kPdgDarkNeutrino);
   }
   if (!anti_dnu_particle) {
     // Name Title Mass Stable Width Charge Class PDG
     fDatabasePDG->AddParticle("nu_D_bar","#bar{#nu}_{D}",reg->GetDouble("Dark-NeutrinoMass"),
                               true,0.,0,"DarkNeutrino",kPdgAntiDarkNeutrino);
   }
   if (!med_particle) {
     // Name Title Mass Stable Width Charge Class PDG
     fDatabasePDG->AddParticle("Z_D","Z_{D}",reg->GetDouble("Dark-MediatorMass"),
                               true,0.,0,"DarkNeutrino",kPdgDNuMediator);
   }
   return true;
 }

bool PDGLibrary::AddHNL ( void )

private

Definition at line 176 of file PDGLibrary.cxx.

References genie::AlgConfigPool::CommonList(), fDatabasePDG, genie::Registry::GetDouble(), genie::AlgConfigPool::Instance(), genie::kPdgHNL, LOG, and pERROR.

Referenced by PDGLibrary().

 {
   // Add HNL to PDG database
   const Registry * reg = AlgConfigPool::Instance()->CommonList("HNL", "ParameterSpace");
   if (!reg) {
     LOG("PDG", pERROR) << "Cannot find HNL ParameterSpace param_set";
     return false;
   }
   TParticlePDG * hnl = fDatabasePDG->GetParticle(kPdgHNL);
   if (!hnl) {
     // Name Title Mass Stable Width Charge Class PDG
     fDatabasePDG->AddParticle("HNL","HNL",reg->GetDouble("HNL-Mass"),true,0.,0,"HNL",kPdgHNL);
     fDatabasePDG->AddParticle("HNLBar","HNLBar",reg->GetDouble("HNL-Mass"),true,0.,0,"HNL",-1*kPdgHNL);
   }
   return true;
 }

TDatabasePDG * PDGLibrary::DBase ( void )

Definition at line 81 of file PDGLibrary.cxx.

References fDatabasePDG.

Referenced by genie::evtlib::EventLibraryInterface::LoadRecords(), and genie::evtlib::EvtLibPXSec::LoadXSecs().

 {
   return fDatabasePDG;
 }

TParticlePDG * PDGLibrary::Find	(	int	pdgc,
		bool	must_exist = `true`
	)

Definition at line 86 of file PDGLibrary.cxx.

References fDatabasePDG, LOG, and pERROR.

 {
 
   TParticlePDG * ret = fDatabasePDG->GetParticle(pdgc);
   if(ret) return ret;
 
   if ( must_exist ) {
     LOG("PDG", pERROR) << "Requested missing particle with PDG: " << pdgc ;
   }
 
   return ret ;
 }

PDGLibrary * PDGLibrary::Instance ( void )

static

Definition at line 68 of file PDGLibrary.cxx.

References genie::PDGLibrary::Cleaner::DummyMethodAndSilentCompiler(), fInstance, LOG, PDGLibrary(), and pINFO.

 {
   if(fInstance == 0) {
     LOG("PDG", pINFO) << "PDGLibrary late initialization";
 
     static PDGLibrary::Cleaner cleaner;
     cleaner.DummyMethodAndSilentCompiler();
 
     fInstance = new PDGLibrary;
   }
   return fInstance;
 }

bool PDGLibrary::LoadDBase ( void )

private

Definition at line 100 of file PDGLibrary.cxx.

References genie::AlgConfigPool::CommonList(), fDatabasePDG, genie::utils::xml::GetString(), genie::AlgConfigPool::Instance(), LOG, pERROR, and pINFO.

Referenced by PDGLibrary(), and ReloadDBase().

 {
   fDatabasePDG = TDatabasePDG::Instance();
 
   // loading PDG data from $GENIE/config/
   const char* altpdgtable = gSystem->Getenv("GENIE_PDG_TABLE");
   if ( altpdgtable ) {
     if ( ! (gSystem->AccessPathName(altpdgtable) ) ) {
         LOG("PDG", pINFO) << "Load PDG data from $GENIE_PDG_TABLE: "
                           << altpdgtable;
         fDatabasePDG->ReadPDGTable( altpdgtable );
         return true;
     }
   }
 
   if ( gSystem->Getenv("GENIE") ) {
     string base_dir = string( gSystem->Getenv("GENIE") );
     base_dir += string("/data/evgen/catalogues/pdg/") ; 
 
     string file_name = "genie_pdg_table.txt" ; 
     const Registry * reg = AlgConfigPool::Instance()->CommonList("Param", "PDG");
     if( reg ) {
       file_name = reg -> GetString("PDG-TableName") ;
       LOG("PDG", pINFO) << "Found file name specification: " << file_name ;
 
     }
     
     string path = base_dir + file_name ;
 
     if ( ! (gSystem->AccessPathName(path.c_str()) ) ) {
         LOG("PDG", pINFO) << "Load PDG data from: " << path;
         fDatabasePDG->ReadPDGTable( path.c_str() );
         return true;
     }
   }
 
   // no PDG data in $GENIE/config/ - Try $ROOTSYS/etc/
 
   if(gSystem->Getenv("ROOTSYS")) {
     string base_dir  = string( gSystem->Getenv("ROOTSYS") );
     string path = base_dir  + string("/etc/pdg_table.txt");
 
     if ( !(gSystem->AccessPathName(path.c_str())) ) {
         LOG("PDG", pINFO) << "Load PDG data from: " << path;
         fDatabasePDG->ReadPDGTable( path.c_str() );
         return true;
      }
   }
 
   LOG("PDG", pERROR) << " *** The PDG extensions will not be loaded!! ***";
   return false;
 };

void PDGLibrary::ReloadDBase ( void )

Definition at line 224 of file PDGLibrary.cxx.

References fDatabasePDG, LoadDBase(), LOG, and pERROR.

Referenced by main().

 {
   if(fDatabasePDG) {
     delete fDatabasePDG;
   }
 
   if( ! LoadDBase() ) LOG("PDG", pERROR) << "Could not load PDG data";
 }

Friends And Related Function Documentation

friend struct Cleaner

friend

Definition at line 72 of file PDGLibrary.h.

Member Data Documentation

TDatabasePDG* genie::PDGLibrary::fDatabasePDG

private

Definition at line 61 of file PDGLibrary.h.

Referenced by AddDarkMatter(), AddDarkSector(), AddHNL(), DBase(), Find(), LoadDBase(), and ReloadDBase().

PDGLibrary * PDGLibrary::fInstance = 0

staticprivate

Definition at line 60 of file PDGLibrary.h.

Referenced by Instance(), PDGLibrary(), genie::PDGLibrary::Cleaner::~Cleaner(), and ~PDGLibrary().

The documentation for this class was generated from the following files:

Classes

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation