#include <DatabasePDG.h>

Public Member Functions
int	CheckImpossibleDecays (bool dump=kFALSE)

	DatabasePDG ()

void	DumpData (bool dumpAll=kFALSE)

char *	GetDecayFilename ()

double	GetMaximumMass ()

double	GetMaximumWidth ()

double	GetMinimumMass ()

double	GetMinimumWidth ()

int	GetNAllowedChannels (ParticlePDG *particle, double motherMass)

int	GetNParticles (bool all=kFALSE)

char *	GetParticleFilename ()

ParticlePDG *	GetPDGParticle (char *name)

ParticlePDG *	GetPDGParticle (int pdg)

ParticlePDG *	GetPDGParticleByIndex (int index)

bool	GetPDGParticleStatus (char *name)

bool	GetPDGParticleStatus (int pdg)

bool	GetPDGParticleStatusByIndex (int index)

bool	GetUseCharmParticles ()

bool	IsChannelAllowed (DecayChannel *channel, double motherMass)

bool	LoadData ()

void	SetDecayFilename (char *filename)

void	SetMassRange (double min, double max)

void	SetMaximumMass (double value)

void	SetMaximumWidth (double value)

void	SetMinimumMass (double value)

void	SetMinimumWidth (double value)

void	SetParticleFilename (char *filename)

void	SetUseCharmParticles (bool flag)

void	SetWidthRange (double min, double max)

void	UseThisListOfParticles (char *filename, bool exclusive=kTRUE)

	~DatabasePDG ()

Private Member Functions
bool	LoadDecays ()

bool	LoadParticles ()

void	SortParticles ()

Private Attributes
char	fDecayFilename [256]

double	fMaximumMass

double	fMaximumWidth

double	fMinimumMass

double	fMinimumWidth

int	fNParticles

char	fParticleFilename [256]

ParticlePDG *	fParticles [kMaxParticles]

bool	fStatus [kMaxParticles]

bool	fUseCharmParticles

Detailed Description

Definition at line 34 of file DatabasePDG.h.

Constructor & Destructor Documentation

◆ DatabasePDG()

DatabasePDG::DatabasePDG ( )

Definition at line 28 of file DatabasePDG.cc.

                          {
   fNParticles = 0;
  
   std::string file1 = "GeneratorInterface/Hydjet2Interface/data/particles.data";
   edm::FileInPath f1(file1);
   particlesDATAstr = ((f1.fullPath()).c_str());
  
   std::string file2 = "GeneratorInterface/Hydjet2Interface/data/tabledecay.txt";
   edm::FileInPath f2(file2);
   tableDECAYstr = ((f2.fullPath()).c_str());
  
   strcpy(fParticleFilename, particlesDATAstr);
   strcpy(fDecayFilename, tableDECAYstr);
   for (int i = 0; i < kMaxParticles; i++) {
     fParticles[i] = new ParticlePDG();
     fStatus[i] = kFALSE;
   }
   fUseCharmParticles = kTRUE;
   fMinimumWidth = 0.;
   fMaximumWidth = 10.;
   fMinimumMass = 0.;
   fMaximumMass = 200.;
 }

References DeadROC_duringRun::f1, DeadROC_duringRun::f2, timingPdfMaker::file1, edmOneToOneComparison::file2, mps_fire::i, kMaxParticles, particlesDATAstr, AlCaHLTBitMon_QueryRunRegistry::string, and tableDECAYstr.

◆ ~DatabasePDG()

DatabasePDG::~DatabasePDG ( )

Definition at line 52 of file DatabasePDG.cc.

                           {
   for (int i = 0; i < kMaxParticles; i++)
     delete fParticles[i];
 }

References mps_fire::i, and kMaxParticles.

Member Function Documentation

◆ CheckImpossibleDecays()

int DatabasePDG::CheckImpossibleDecays ( bool dump = kFALSE )

Definition at line 385 of file DatabasePDG.cc.

                                                 {
   // Check the database for impossible decays
   int nImpossibleDecays = 0;
   for (int currPart = 0; currPart < fNParticles; currPart++) {
     if (!fStatus[currPart])
       continue;
     int allChannels = fParticles[currPart]->GetNDecayChannels();
     int allowedChannels = GetNAllowedChannels(fParticles[currPart], fParticles[currPart]->GetMass());
     if (dump) {
       cout << "Particle " << fParticles[currPart]->GetPDG() << " has " << allChannels
            << " decay channels specified in the database" << endl;
       cout << " Allowed channels assuming table mass = " << allowedChannels << endl;
     }
     if (dump && allChannels > 0 && allowedChannels == 0) {
       cout << "**********************************************************************" << endl;
       cout << "       All channels for this particles are not allowed" << endl;
       cout << "**********************************************************************" << endl;
     }
     if (dump && fParticles[currPart]->GetWidth() > 0. && allChannels == 0) {
       cout << "**********************************************************************" << endl;
       cout << "    Particle has finite width but no decay channels specified" << endl;
       cout << "**********************************************************************" << endl;
     }
     for (int currChannel = 0; currChannel < fParticles[currPart]->GetNDecayChannels(); currChannel++) {
       double motherMass = fParticles[currPart]->GetMass();
       double daughtersSumMass = 0.;
       for (int currDaughter = 0; currDaughter < fParticles[currPart]->GetDecayChannel(currChannel)->GetNDaughters();
            currDaughter++) {
         ParticlePDG *daughter =
             GetPDGParticle(fParticles[currPart]->GetDecayChannel(currChannel)->GetDaughterPDG(currDaughter));
         daughtersSumMass += daughter->GetMass();
       }
       if (daughtersSumMass >= motherMass) {
         nImpossibleDecays++;
         if (dump) {
           cout << "Imposible decay for particle " << fParticles[currPart]->GetPDG() << endl;
           cout << "  Channel: " << fParticles[currPart]->GetPDG() << " --> ";
           for (int currDaughter = 0; currDaughter < fParticles[currPart]->GetDecayChannel(currChannel)->GetNDaughters();
                currDaughter++) {
             ParticlePDG *daughter =
                 GetPDGParticle(fParticles[currPart]->GetDecayChannel(currChannel)->GetDaughterPDG(currDaughter));
             cout << daughter->GetPDG() << " ";
           }
           cout << endl;
           cout << "  Mother particle mass = " << motherMass << endl;
           cout << "  Daughters sum mass   = " << daughtersSumMass << endl;
         }
       }
     }
   }
   return nImpossibleDecays;
 }

References gather_cfg::cout, FrontierConditions_GlobalTag_cff::dump, ParticlePDG::GetMass(), and ParticlePDG::GetPDG().

◆ DumpData()

void DatabasePDG::DumpData ( bool dumpAll = kFALSE )

Definition at line 320 of file DatabasePDG.cc.

                                        {
   cout << "***********************************************************************************************************"
        << endl;
   cout << "Dumping all the information contained in the database..." << endl;
   int nDecays = 0;
   int nGoodStatusDecays = 0;
   int nGoodStatusParticles = 0;
   for (int currPart = 0; currPart < fNParticles; currPart++) {
     nGoodStatusParticles += (fStatus[currPart] ? 1 : 0);
     nGoodStatusDecays += (fStatus[currPart] ? fParticles[currPart]->GetNDecayChannels() : 0);
     nDecays += fParticles[currPart]->GetNDecayChannels();
     if (!(dumpAll || (!dumpAll && fStatus[currPart])))
       continue;
     cout << "###### Particle: " << fParticles[currPart]->GetName() << " with PDG code "
          << fParticles[currPart]->GetPDG() << endl;
     cout << "   status          = " << fStatus[currPart] << endl;
     cout << "   mass            = " << fParticles[currPart]->GetMass() << " GeV" << endl;
     cout << "   width           = " << fParticles[currPart]->GetWidth() << " GeV" << endl;
     cout << "   2*spin          = " << int(2. * fParticles[currPart]->GetSpin()) << endl;
     cout << "   2*isospin       = " << int(2. * fParticles[currPart]->GetIsospin()) << endl;
     cout << "   2*isospin3      = " << int(2. * fParticles[currPart]->GetIsospinZ()) << endl;
     cout << "   u,d quarks      = " << int(fParticles[currPart]->GetLightQNumber()) << endl;
     cout << "   s quarks        = " << int(fParticles[currPart]->GetStrangeQNumber()) << endl;
     cout << "   c quarks        = " << int(fParticles[currPart]->GetCharmQNumber()) << endl;
     cout << "   anti u,d quarks = " << int(fParticles[currPart]->GetLightAQNumber()) << endl;
     cout << "   anti s quarks   = " << int(fParticles[currPart]->GetStrangeAQNumber()) << endl;
     cout << "   anti c quarks   = " << int(fParticles[currPart]->GetCharmAQNumber()) << endl;
     cout << "   baryon number   = " << int(fParticles[currPart]->GetBaryonNumber()) << endl;
     cout << "   strangeness     = " << int(fParticles[currPart]->GetStrangeness()) << endl;
     cout << "   charmness       = " << int(fParticles[currPart]->GetCharmness()) << endl;
     cout << "   electric charge = " << int(fParticles[currPart]->GetElectricCharge()) << endl;
     cout << "   full branching  = " << fParticles[currPart]->GetFullBranching() << endl;
     cout << "   decay modes     = " << fParticles[currPart]->GetNDecayChannels() << endl;
     for (int currChannel = 0; currChannel < fParticles[currPart]->GetNDecayChannels(); currChannel++) {
       cout << "   channel " << currChannel + 1 << " with branching "
            << fParticles[currPart]->GetDecayChannel(currChannel)->GetBranching() << endl;
       cout << "   daughters PDG codes: ";
       double daughtersMass = 0.0;
       for (int currDaughter = 0; currDaughter < fParticles[currPart]->GetDecayChannel(currChannel)->GetNDaughters();
            currDaughter++) {
         cout << fParticles[currPart]->GetDecayChannel(currChannel)->GetDaughterPDG(currDaughter) << "\t";
         ParticlePDG *daughter =
             GetPDGParticle(fParticles[currPart]->GetDecayChannel(currChannel)->GetDaughterPDG(currDaughter));
         daughtersMass += daughter->GetMass();
       }
       cout << endl;
       cout << "   daughters sum mass = " << daughtersMass << endl;
     }
   }
   if (dumpAll) {
     cout << "Finished dumping information for " << fNParticles << " particles with " << nDecays
          << " decay channels in total." << endl;
     cout << "**********************************************************************************************************"
             "***"
          << endl;
   } else {
     cout << "Finished dumping information for " << nGoodStatusParticles << "(" << fNParticles << ")"
          << " particles with " << nGoodStatusDecays << "(" << nDecays << ")"
          << " decay channels in total." << endl;
     cout << "**********************************************************************************************************"
             "***"
          << endl;
   }
 }

References gather_cfg::cout, ParticlePDG::GetMass(), and createfilelist::int.

◆ GetDecayFilename()

char* DatabasePDG::GetDecayFilename ( )

inline

Definition at line 82 of file DatabasePDG.h.

82 { return fDecayFilename; }

References fDecayFilename.

◆ GetMaximumMass()

double DatabasePDG::GetMaximumMass ( )

inline

Definition at line 94 of file DatabasePDG.h.

94 { return fMaximumMass; };

References fMaximumMass.

◆ GetMaximumWidth()

double DatabasePDG::GetMaximumWidth ( )

inline

Definition at line 92 of file DatabasePDG.h.

92 { return fMaximumWidth; };

References fMaximumWidth.

◆ GetMinimumMass()

double DatabasePDG::GetMinimumMass ( )

inline

Definition at line 93 of file DatabasePDG.h.

93 { return fMinimumMass; };

References fMinimumMass.

◆ GetMinimumWidth()

double DatabasePDG::GetMinimumWidth ( )

inline

Definition at line 91 of file DatabasePDG.h.

91 { return fMinimumWidth; };

References fMinimumWidth.

◆ GetNAllowedChannels()

int DatabasePDG::GetNAllowedChannels	(	ParticlePDG *	particle,
		double	motherMass
	)

Definition at line 657 of file DatabasePDG.cc.

                                                                              {
   int nAllowedChannels = 0;
   for (int i = 0; i < particle->GetNDecayChannels(); i++)
     nAllowedChannels += (IsChannelAllowed(particle->GetDecayChannel(i), motherMass) ? 1 : 0);
   return nAllowedChannels;
 }

References ParticlePDG::GetDecayChannel(), ParticlePDG::GetNDecayChannels(), and mps_fire::i.

Referenced by Decay().

◆ GetNParticles()

int DatabasePDG::GetNParticles ( bool all = kFALSE )

Definition at line 582 of file DatabasePDG.cc.

                                        {
   if (all)
     return fNParticles;
  
   int nGoodStatus = 0;
   for (int i = 0; i < fNParticles; i++)
     if (fStatus[i])
       nGoodStatus++;
   return nGoodStatus;
 }

References python.cmstools::all(), and mps_fire::i.

Referenced by GrandCanonical::BaryonDensity(), GrandCanonical::CharmDensity(), GrandCanonical::ElectroDensity(), GrandCanonical::EnergyDensity(), gen::Hydjet2Hadronizer::initializeForInternalPartons(), GrandCanonical::StrangeDensity(), and NAStrangeDensity::StrangenessDensity().

◆ GetParticleFilename()

char* DatabasePDG::GetParticleFilename ( )

inline

Definition at line 81 of file DatabasePDG.h.

81 { return fParticleFilename; }

References fParticleFilename.

◆ GetPDGParticle() [1/2]

ParticlePDG * DatabasePDG::GetPDGParticle ( char * name )

Definition at line 268 of file DatabasePDG.cc.

                                                    {
   int nFindings = 0;
   int firstTimeIndex = 0;
   for (int i = 0; i < fNParticles; i++) {
     if (!strcmp(name, fParticles[i]->GetName())) {
       if (nFindings == 0)
         firstTimeIndex = i;
       nFindings++;
     }
   }
   if (nFindings == 1)
     return fParticles[firstTimeIndex];
   if (nFindings == 0) {
     edm::LogWarning("DatabasePDG") << "The particle required with name (" << name
                                    << ") was not found in the database!!";
     return nullptr;
   }
   if (nFindings >= 2) {
     edm::LogWarning("DatabasePDG") << "The particle required with name (" << name << ") was found with " << nFindings
                                    << " entries in the database. Check it out !!" << endl
                                    << "Returning the first instance found";
     return fParticles[firstTimeIndex];
   }
   return nullptr;
 }

References mps_fire::i, and Skims_PA_cff::name.

◆ GetPDGParticle() [2/2]

ParticlePDG * DatabasePDG::GetPDGParticle ( int pdg )

Definition at line 218 of file DatabasePDG.cc.

                                                 {
   int nFindings = 0;
   int firstTimeIndex = 0;
   for (int i = 0; i < fNParticles; i++) {
     if (pdg == fParticles[i]->GetPDG()) {
       if (nFindings == 0)
         firstTimeIndex = i;
       nFindings++;
     }
   }
   if (nFindings == 1)
     return fParticles[firstTimeIndex];
   if (nFindings == 0) {
     edm::LogWarning("DatabasePDG") << "The particle required with PDG: " << pdg << " was not found in the database!!";
     return nullptr;
   }
   if (nFindings >= 2) {
     edm::LogWarning("DatabasePDG") << "The particle required with PDG: " << pdg << " was found with " << nFindings
                                    << " entries in the database. Check it out !!" << endl
                                    << "Returning the first instance found";
     return fParticles[firstTimeIndex];
   }
   return nullptr;
 }

References mps_fire::i.

Referenced by Decay(), gen::Hydjet2Hadronizer::generatePartonsAndHadronize(), and gen::Hydjet2Hadronizer::initializeForInternalPartons().

◆ GetPDGParticleByIndex()

ParticlePDG * DatabasePDG::GetPDGParticleByIndex ( int index )

Definition at line 198 of file DatabasePDG.cc.

                                                          {
   if (index < 0 || index > fNParticles) {
     edm::LogWarning("DatabasePDG") << "Particle index is negative or too big !!" << endl
                                    << " It must be inside this range: (0, " << fNParticles - 1 << ")" << endl
                                    << " Returning null pointer!!";
     return nullptr;
   }
   return fParticles[index];
 }

Referenced by GrandCanonical::BaryonDensity(), GrandCanonical::CharmDensity(), GrandCanonical::ElectroDensity(), GrandCanonical::EnergyDensity(), gen::Hydjet2Hadronizer::initializeForInternalPartons(), GrandCanonical::StrangeDensity(), and NAStrangeDensity::StrangenessDensity().

◆ GetPDGParticleStatus() [1/2]

bool DatabasePDG::GetPDGParticleStatus ( char * name )

Definition at line 294 of file DatabasePDG.cc.

                                                  {
   int nFindings = 0;
   int firstTimeIndex = 0;
   for (int i = 0; i < fNParticles; i++) {
     if (!strcmp(name, fParticles[i]->GetName())) {
       if (nFindings == 0)
         firstTimeIndex = i;
       nFindings++;
     }
   }
   if (nFindings == 1)
     return fStatus[firstTimeIndex];
   if (nFindings == 0) {
     edm::LogWarning("DatabasePDG") << "The particle required with name (" << name
                                    << ") was not found in the database!!";
     return kFALSE;
   }
   if (nFindings >= 2) {
     edm::LogWarning("DatabasePDG") << "The particle status required for name (" << name << ") was found with "
                                    << nFindings << " entries in the database. Check it out !!" << endl
                                    << "Returning the first instance found";
     return fStatus[firstTimeIndex];
   }
   return kFALSE;
 }

References mps_fire::i, and Skims_PA_cff::name.

◆ GetPDGParticleStatus() [2/2]

bool DatabasePDG::GetPDGParticleStatus ( int pdg )

Definition at line 243 of file DatabasePDG.cc.

                                               {
   int nFindings = 0;
   int firstTimeIndex = 0;
   for (int i = 0; i < fNParticles; i++) {
     if (pdg == fParticles[i]->GetPDG()) {
       if (nFindings == 0)
         firstTimeIndex = i;
       nFindings++;
     }
   }
   if (nFindings == 1)
     return fStatus[firstTimeIndex];
   if (nFindings == 0) {
     edm::LogWarning("DatabasePDG") << "The particle required with PDG: " << pdg << " was not found in the database!!";
     return kFALSE;
   }
   if (nFindings >= 2) {
     edm::LogWarning("DatabasePDG") << "The particle status required for PDG: " << pdg << " was found with " << nFindings
                                    << " entries in the database. Check it out !!" << endl
                                    << "Returning the status of first instance found";
     return fStatus[firstTimeIndex];
   }
   return kFALSE;
 }

References mps_fire::i.

◆ GetPDGParticleStatusByIndex()

bool DatabasePDG::GetPDGParticleStatusByIndex ( int index )

Definition at line 208 of file DatabasePDG.cc.

                                                        {
   if (index < 0 || index > fNParticles) {
     edm::LogWarning("DatabasePDG") << "Particle index is negative or too big !!" << endl
                                    << " It must be inside this range: (0, " << fNParticles - 1 << ")" << endl
                                    << " Returning null pointer!!";
     return kFALSE;
   }
   return fStatus[index];
 }

◆ GetUseCharmParticles()

bool DatabasePDG::GetUseCharmParticles ( )

inline

Definition at line 90 of file DatabasePDG.h.

90 { return fUseCharmParticles; };

References fUseCharmParticles.

◆ IsChannelAllowed()

bool DatabasePDG::IsChannelAllowed	(	DecayChannel *	channel,
		double	motherMass
	)

Definition at line 648 of file DatabasePDG.cc.

                                                                            {
   double daughtersSumMass = 0.0;
   for (int i = 0; i < channel->GetNDaughters(); i++)
     daughtersSumMass += GetPDGParticle(channel->GetDaughterPDG(i))->GetMass();
   if (daughtersSumMass <= motherMass)
     return kTRUE;
   return kFALSE;
 }

References DecayChannel::GetDaughterPDG(), DecayChannel::GetNDaughters(), and mps_fire::i.

◆ LoadData()

bool DatabasePDG::LoadData ( )

Definition at line 61 of file DatabasePDG.cc.

61 { return (LoadParticles() && LoadDecays()); }

Referenced by InitialState::InitialState().

◆ LoadDecays()

bool DatabasePDG::LoadDecays ( )

private

Definition at line 138 of file DatabasePDG.cc.

                              {
   ifstream decayFile;
   decayFile.open(fDecayFilename);
   if (!decayFile) {
     edm::LogError("DatabasePDG") << "The ASCII file containing the decays list " << fDecayFilename << " was not found";
     return kFALSE;
   }
  
   int mother_pdg, daughter_pdg[3];
   double branching;
  
   decayFile.exceptions(ios::failbit);
   while (!decayFile.eof()) {
     mother_pdg = 0;
     for (int i = 0; i < 3; i++)
       daughter_pdg[i] = 0;
     branching = -1.0;
     try {
       decayFile >> mother_pdg;
       for (int i = 0; i < 3; i++)
         decayFile >> daughter_pdg[i];
       decayFile >> branching;
     } catch (ios::failure const &problem) {
       LogDebug("DatabasePDG") << " ios:failure in decay file " << problem.what();
       break;
     }
     if ((mother_pdg != 0) && (daughter_pdg[0] != 0) && (branching >= 0)) {
       int nDaughters = 0;
       for (int i = 0; i < 3; i++)
         if (daughter_pdg[i] != 0)
           nDaughters++;
       ParticlePDG *particle = GetPDGParticle(mother_pdg);
       if (!particle) {
         LogWarning("DatabasePDG") << " Mother particle PDG (" << mother_pdg
                                   << ") not found in the particle definition list:" << mother_pdg << " >>> ";
         for (int kk = 0; kk < nDaughters; kk++)
           LogWarning("DatabasePDG") << daughter_pdg[kk] << "  ";
         return kFALSE;
       }
       for (int kk = 0; kk < nDaughters; kk++) {
         if (!GetPDGParticle(daughter_pdg[kk])) {
           LogWarning("DatabasePDG") << "Daughter particle PDG (" << daughter_pdg[kk]
                                     << ") not found in the particle definition list: " << mother_pdg << ">>> ";
           for (int kkk = 0; kkk < nDaughters; kkk++)
             LogWarning("DatabasePDG") << daughter_pdg[kkk] << "  ";
         }
       }
       DecayChannel decay(mother_pdg, branching, nDaughters, daughter_pdg);
       particle->AddChannel(decay);
     }
   }
   decayFile.close();
   int nDecayChannels = 0;
   for (int i = 0; i < fNParticles; i++) {
     nDecayChannels += fParticles[i]->GetNDecayChannels();
   }
   edm::LogInfo("DatabasePDG") << "Number of decays found in the database is " << nDecayChannels;
   return kTRUE;
 }

References ParticlePDG::AddChannel(), PA_ZEESkim_cff::decay, mps_fire::i, GetRecoTauVFromDQM_MC_cff::kk, and LogDebug.

◆ LoadParticles()

bool DatabasePDG::LoadParticles ( )

private

Definition at line 63 of file DatabasePDG.cc.

                                 {
   ifstream particleFile;
   particleFile.open(fParticleFilename);
   if (!particleFile) {
     edm::LogError("DatabasePDG") << "The ASCII file containing the PDG particle list " << fParticleFilename
                                  << " was not found";
     return kFALSE;
   }
  
   char name[9];
   double mass, width, spin, isospin, isospinZ, q, s, aq, as, c, ac;
   int pdg;
   int goodStatusParticles = 0;
  
   edm::LogInfo("DatabasePDG") << "Start loading particles with the following criteria:" << endl
                               << "       Use particles containing charm quarks (1:yes;0:no) : " << fUseCharmParticles
                               << endl
                               << "       Mass range                                         : (" << fMinimumMass << "; "
                               << fMaximumMass << ")" << endl
                               << "       Width range                                        : (" << fMinimumWidth
                               << "; " << fMaximumWidth << ")";
  
   particleFile.exceptions(ios::failbit);
   while (!particleFile.eof()) {
     try {
       particleFile >> name >> mass >> width >> spin >> isospin >> isospinZ >> q >> s >> aq >> as >> c >> ac >> pdg;
     } catch (ios::failure const &problem) {
       LogDebug("DatabasePDG") << " ios:failure in particle file " << problem.what();
       break;
     }
  
     fParticles[fNParticles]->SetName(name);
     fParticles[fNParticles]->SetPDG(pdg);
     fParticles[fNParticles]->SetMass(mass);
     fParticles[fNParticles]->SetWidth(width);
     fParticles[fNParticles]->SetSpin(spin);
     fParticles[fNParticles]->SetIsospin(isospin);
     fParticles[fNParticles]->SetIsospinZ(isospinZ);
     fParticles[fNParticles]->SetLightQNumber(q);
     fParticles[fNParticles]->SetStrangeQNumber(s);
     fParticles[fNParticles]->SetLightAQNumber(aq);
     fParticles[fNParticles]->SetStrangeAQNumber(as);
     fParticles[fNParticles]->SetCharmQNumber(c);
     fParticles[fNParticles]->SetCharmAQNumber(ac);
     goodStatusParticles++;
     fStatus[fNParticles] = kTRUE;
     // check if we want charmed particles
     if (!fUseCharmParticles && (c > 0 || ac > 0)) {
       fStatus[fNParticles] = kFALSE;
       goodStatusParticles--;
     }
     // check that the particle mass is inside accepted limits
     if (!(fMinimumMass <= mass && mass <= fMaximumMass)) {
       fStatus[fNParticles] = kFALSE;
       goodStatusParticles--;
     }
     // check that the particle width is inside accepted limits
     if (!(fMinimumWidth <= width && width <= fMaximumWidth)) {
       fStatus[fNParticles] = kFALSE;
       goodStatusParticles--;
     }
  
     fNParticles++;
   }
   particleFile.close();
   if (fNParticles == 0) {
     LogWarning("DatabasePDG") << " No particles were found in the file specified!!";
     return kFALSE;
   }
   SortParticles();
   edm::LogInfo("DatabasePDG") << " Particle definitions found: " << fNParticles
                               << ". Good status particles: " << goodStatusParticles;
   return kTRUE;
 }

References HltBtagPostValidation_cff::c, LogDebug, EgHLTOffHistBins_cfi::mass, Skims_PA_cff::name, dileptonTrigSettings_cff::pdg, submitPVResolutionJobs::q, alignCSCRings::s, and ApeEstimator_cff::width.

◆ SetDecayFilename()

void DatabasePDG::SetDecayFilename ( char * filename )

Definition at line 59 of file DatabasePDG.cc.

59 { strcpy(fDecayFilename, filename); }

References corrVsCorr::filename.

◆ SetMassRange()

void DatabasePDG::SetMassRange	(	double	min,
		double	max
	)

Definition at line 527 of file DatabasePDG.cc.

                                                      {
   if (fNParticles > 0) {
     fMinimumMass = min;
     fMaximumMass = max;
     for (int i = 0; i < fNParticles; i++) {
       if ((fParticles[i]->GetMass() < fMinimumMass) || (fParticles[i]->GetMass() > fMaximumMass))
         fStatus[i] = kFALSE;
     }
     SortParticles();
  
     return;
   } else {
     fMinimumMass = min;
     fMaximumMass = max;
   }
  
   return;
 }

References mps_fire::i, SiStripPI::max, and min().

Referenced by InitialState::InitialState().

◆ SetMaximumMass()

void DatabasePDG::SetMaximumMass ( double value )

Definition at line 513 of file DatabasePDG.cc.

                                              {
   if (fNParticles > 0) {
     fMaximumMass = value;
     for (int i = 0; i < fNParticles; i++) {
       if (fParticles[i]->GetMass() > fMaximumMass)
         fStatus[i] = kFALSE;
     }
     SortParticles();
     return;
   } else
     fMaximumMass = value;
   return;
 }

References mps_fire::i, and relativeConstraints::value.

◆ SetMaximumWidth()

void DatabasePDG::SetMaximumWidth ( double value )

Definition at line 466 of file DatabasePDG.cc.

                                               {
   if (fNParticles > 0) {
     fMaximumWidth = value;
     for (int i = 0; i < fNParticles; i++) {
       if (fParticles[i]->GetWidth() > fMaximumWidth)
         fStatus[i] = kFALSE;
     }
     SortParticles();
     return;
   } else
     fMaximumWidth = value;
   return;
 }

References mps_fire::i, and relativeConstraints::value.

◆ SetMinimumMass()

void DatabasePDG::SetMinimumMass ( double value )

Definition at line 499 of file DatabasePDG.cc.

                                              {
   if (fNParticles > 0) {
     fMinimumMass = value;
     for (int i = 0; i < fNParticles; i++) {
       if (fParticles[i]->GetMass() < fMinimumMass)
         fStatus[i] = kFALSE;
     }
     SortParticles();
     return;
   } else
     fMinimumMass = value;
   return;
 }

References mps_fire::i, and relativeConstraints::value.

◆ SetMinimumWidth()

void DatabasePDG::SetMinimumWidth ( double value )

Definition at line 452 of file DatabasePDG.cc.

                                               {
   if (fNParticles > 0) {
     fMinimumWidth = value;
     for (int i = 0; i < fNParticles; i++) {
       if (fParticles[i]->GetWidth() < fMinimumWidth)
         fStatus[i] = kFALSE;
     }
     SortParticles();
     return;
   } else
     fMinimumWidth = value;
   return;
 }

References mps_fire::i, and relativeConstraints::value.

◆ SetParticleFilename()

void DatabasePDG::SetParticleFilename ( char * filename )

Definition at line 57 of file DatabasePDG.cc.

57 { strcpy(fParticleFilename, filename); }

References corrVsCorr::filename.

◆ SetUseCharmParticles()

void DatabasePDG::SetUseCharmParticles ( bool flag )

Definition at line 438 of file DatabasePDG.cc.

                                                 {
   if (fNParticles > 0) {
     fUseCharmParticles = flag;
     for (int i = 0; i < fNParticles; i++) {
       if (fParticles[i]->GetCharmQNumber() > 0 || fParticles[i]->GetCharmAQNumber())
         fStatus[i] = flag;
     }
     SortParticles();
     return;
   } else
     fUseCharmParticles = flag;
   return;
 }

References RemoveAddSevLevel::flag, and mps_fire::i.

◆ SetWidthRange()

void DatabasePDG::SetWidthRange	(	double	min,
		double	max
	)

Definition at line 480 of file DatabasePDG.cc.

                                                       {
   if (fNParticles > 0) {
     fMinimumWidth = min;
     fMaximumWidth = max;
     for (int i = 0; i < fNParticles; i++) {
       if ((fParticles[i]->GetWidth() < fMinimumWidth) || (fParticles[i]->GetWidth() > fMaximumWidth))
         fStatus[i] = kFALSE;
     }
     SortParticles();
  
     return;
   } else {
     fMinimumWidth = min;
     fMaximumWidth = max;
   }
  
   return;
 }

References mps_fire::i, SiStripPI::max, and min().

Referenced by InitialState::InitialState().

◆ SortParticles()

void DatabasePDG::SortParticles ( )

private

Definition at line 546 of file DatabasePDG.cc.

                                 {
   if (fNParticles < 2) {
     edm::LogWarning("DatabasePDG") << "No particles to sort. Load data first!!";
     return;
   }
  
   int nGoodStatus = 0;
   for (int i = 0; i < fNParticles; i++)
     if (fStatus[i])
       nGoodStatus++;
  
   if (nGoodStatus == fNParticles)  // if all particles have good status then there is nothing to do
     return;
  
   if (nGoodStatus == 0)  // no good status particles, again nothing to do
     return;
  
   int shifts = 1;
   while (shifts) {
     shifts = 0;
     for (int i = 0; i < fNParticles - 1; i++) {
       if (!fStatus[i] && fStatus[i + 1]) {  // switch if false status is imediately before a true status particle
         ParticlePDG *temporaryPointer = fParticles[i];
         fParticles[i] = fParticles[i + 1];
         fParticles[i + 1] = temporaryPointer;
         bool temporaryStatus = fStatus[i];
         fStatus[i] = fStatus[i + 1];
         fStatus[i + 1] = temporaryStatus;
         shifts++;
       }
     }
   }
  
   return;
 }

References mps_fire::i.

◆ UseThisListOfParticles()

void DatabasePDG::UseThisListOfParticles	(	char *	filename,
		bool	exclusive = `kTRUE`
	)

Definition at line 593 of file DatabasePDG.cc.

                                                                        {
   if (fNParticles < 1) {
     edm::LogError("DatabasePDG") << "You must load the data before calling this function!!";
     return;
   }
  
   ifstream listFile;
   listFile.open(filename);
   if (!listFile) {
     edm::LogError("DatabasePDG") << "The ASCII file containing the PDG codes list (" << filename
                                  << ") was not found !!";
     return;
   }
  
   bool flaggedIndexes[kMaxParticles];
   for (int i = 0; i < kMaxParticles; i++)
     flaggedIndexes[i] = kFALSE;
   int pdg = 0;
   listFile.exceptions(ios::failbit);
   while (!listFile.eof()) {
     try {
       listFile >> pdg;
     } catch (ios::failure const &problem) {
       LogDebug("DatabasePDG") << "ios:failure in list file" << problem.what();
       break;
     }
     int found = 0;
     for (int i = 0; i < fNParticles; i++) {
       if (fParticles[i]->GetPDG() == pdg) {
         found++;
         flaggedIndexes[i] = kTRUE;
       }
     }
     if (!found) {
       edm::LogWarning("DatabasePDG") << "The particle with PDG code " << pdg
                                      << " was asked but not found in the database!!";
     }
     if (found > 1) {
       edm::LogWarning("DatabasePDG") << "The particle with PDG code " << pdg
                                      << " was found more than once in the database!!";
     }
   }
  
   if (exclusive) {
     for (int i = 0; i < kMaxParticles; i++)
       fStatus[i] = flaggedIndexes[i];
   } else {
     for (int i = 0; i < kMaxParticles; i++)
       fStatus[i] = (fStatus[i] && flaggedIndexes[i]);
   }
   SortParticles();
  
   return;
 }

References generator_cfi::exclusive, corrVsCorr::filename, newFWLiteAna::found, mps_fire::i, kMaxParticles, LogDebug, and dileptonTrigSettings_cff::pdg.

Member Data Documentation

◆ fDecayFilename

char DatabasePDG::fDecayFilename[256]

private

Definition at line 40 of file DatabasePDG.h.

Referenced by GetDecayFilename().

◆ fMaximumMass

double DatabasePDG::fMaximumMass

private

Definition at line 45 of file DatabasePDG.h.

Referenced by GetMaximumMass().

◆ fMaximumWidth

double DatabasePDG::fMaximumWidth

private

Definition at line 43 of file DatabasePDG.h.

Referenced by GetMaximumWidth().

◆ fMinimumMass

double DatabasePDG::fMinimumMass

private

Definition at line 44 of file DatabasePDG.h.

Referenced by GetMinimumMass().

◆ fMinimumWidth

double DatabasePDG::fMinimumWidth

private

Definition at line 42 of file DatabasePDG.h.

Referenced by GetMinimumWidth().

◆ fNParticles

int DatabasePDG::fNParticles

private

Definition at line 36 of file DatabasePDG.h.

◆ fParticleFilename

char DatabasePDG::fParticleFilename[256]

private

Definition at line 39 of file DatabasePDG.h.

Referenced by GetParticleFilename().

◆ fParticles

ParticlePDG* DatabasePDG::fParticles[kMaxParticles]

private

Definition at line 37 of file DatabasePDG.h.

◆ fStatus

bool DatabasePDG::fStatus[kMaxParticles]

private

Definition at line 38 of file DatabasePDG.h.

◆ fUseCharmParticles

bool DatabasePDG::fUseCharmParticles

private

Definition at line 41 of file DatabasePDG.h.

Referenced by GetUseCharmParticles().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DatabasePDG()

◆ ~DatabasePDG()

Member Function Documentation

◆ CheckImpossibleDecays()

◆ DumpData()

◆ GetDecayFilename()

◆ GetMaximumMass()

◆ GetMaximumWidth()

◆ GetMinimumMass()

◆ GetMinimumWidth()

◆ GetNAllowedChannels()

◆ GetNParticles()

◆ GetParticleFilename()

◆ GetPDGParticle() [1/2]

◆ GetPDGParticle() [2/2]

◆ GetPDGParticleByIndex()

◆ GetPDGParticleStatus() [1/2]

◆ GetPDGParticleStatus() [2/2]

◆ GetPDGParticleStatusByIndex()

◆ GetUseCharmParticles()

◆ IsChannelAllowed()

◆ LoadData()

◆ LoadDecays()

◆ LoadParticles()

◆ SetDecayFilename()

◆ SetMassRange()

◆ SetMaximumMass()

◆ SetMaximumWidth()

◆ SetMinimumMass()

◆ SetMinimumWidth()

◆ SetParticleFilename()

◆ SetUseCharmParticles()

◆ SetWidthRange()

◆ SortParticles()

◆ UseThisListOfParticles()

Member Data Documentation

◆ fDecayFilename

◆ fMaximumMass

◆ fMaximumWidth

◆ fMinimumMass

◆ fMinimumWidth

◆ fNParticles

◆ fParticleFilename

◆ fParticles

◆ fStatus

◆ fUseCharmParticles