#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 (int pdg)

ParticlePDG *	GetPDGParticle (char *name)

ParticlePDG *	GetPDGParticleByIndex (int index)

bool	GetPDGParticleStatus (int pdg)

bool	GetPDGParticleStatus (char *name)

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 ( )

Definition at line 29 of file DatabasePDG.cc.

References connectstrParser::f1, connectstrParser::f2, edmOneToOneComparison::file1, edmOneToOneComparison::file2, edm::FileInPath::fullPath(), mps_fire::i, kMaxParticles, particlesDATAstr, AlCaHLTBitMon_QueryRunRegistry::string, and tableDECAYstr.

                          {
   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.;
 }

DatabasePDG::~DatabasePDG ( )

Definition at line 53 of file DatabasePDG.cc.

References mps_fire::i, and kMaxParticles.

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

Member Function Documentation

int DatabasePDG::CheckImpossibleDecays ( bool dump = kFALSE )

Definition at line 372 of file DatabasePDG.cc.

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

Referenced by GetMaximumMass().

                                                 {
   // 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;
 }

void DatabasePDG::DumpData ( bool dumpAll = kFALSE )

Definition at line 318 of file DatabasePDG.cc.

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

Referenced by GetMaximumMass().

                                        {
   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;
   }
 }

char* DatabasePDG::GetDecayFilename ( )

inline

Definition at line 82 of file DatabasePDG.h.

References Vispa.Plugins.EdmBrowser.EdmDataAccessor::all(), fDecayFilename, GetNParticles(), GetPDGParticle(), GetPDGParticleByIndex(), GetPDGParticleStatus(), GetPDGParticleStatusByIndex(), diffTreeTool::index, and dataset::name.

82 {return fDecayFilename;}

DatabasePDG::fDecayFilename

char fDecayFilename[256]

Definition: DatabasePDG.h:40

double DatabasePDG::GetMaximumMass ( )

inline

Definition at line 94 of file DatabasePDG.h.

References CheckImpossibleDecays(), FrontierConditions_GlobalTag_cff::dump, DumpData(), fMaximumMass, GetNAllowedChannels(), and IsChannelAllowed().

94 {return fMaximumMass;};

DatabasePDG::fMaximumMass

double fMaximumMass

Definition: DatabasePDG.h:45

double DatabasePDG::GetMaximumWidth ( )

inline

Definition at line 92 of file DatabasePDG.h.

References fMaximumWidth.

92 {return fMaximumWidth;};

DatabasePDG::fMaximumWidth

double fMaximumWidth

Definition: DatabasePDG.h:43

double DatabasePDG::GetMinimumMass ( )

inline

Definition at line 93 of file DatabasePDG.h.

References fMinimumMass.

93 {return fMinimumMass;};

DatabasePDG::fMinimumMass

double fMinimumMass

Definition: DatabasePDG.h:44

double DatabasePDG::GetMinimumWidth ( )

inline

Definition at line 91 of file DatabasePDG.h.

References fMinimumWidth.

91 {return fMinimumWidth;};

DatabasePDG::fMinimumWidth

double fMinimumWidth

Definition: DatabasePDG.h:42

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

Definition at line 650 of file DatabasePDG.cc.

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

Referenced by Decay(), and GetMaximumMass().

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

int DatabasePDG::GetNParticles ( bool all = kFALSE )

Definition at line 574 of file DatabasePDG.cc.

References mps_fire::i.

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

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

char* DatabasePDG::GetParticleFilename ( )

inline

Definition at line 81 of file DatabasePDG.h.

References fParticleFilename.

81 {return fParticleFilename;}

DatabasePDG::fParticleFilename

char fParticleFilename[256]

Definition: DatabasePDG.h:39

ParticlePDG * DatabasePDG::GetPDGParticle ( int pdg )

Definition at line 222 of file DatabasePDG.cc.

References mps_fire::i.

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

                                                 {
   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 0x0;
   }
   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 0x0;
 }

ParticlePDG * DatabasePDG::GetPDGParticle ( char * name )

Definition at line 270 of file DatabasePDG.cc.

References mps_fire::i.

                                                    {
   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 0x0;
   }
   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 0x0;
 }

ParticlePDG * DatabasePDG::GetPDGParticleByIndex ( int index )

Definition at line 202 of file DatabasePDG.cc.

References diffTreeTool::index.

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

                                                          {
   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 0x0;
   }
   return fParticles[index];
 }

bool DatabasePDG::GetPDGParticleStatus ( int pdg )

Definition at line 246 of file DatabasePDG.cc.

References mps_fire::i.

Referenced by GetDecayFilename().

                                               {
   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;
 }

bool DatabasePDG::GetPDGParticleStatus ( char * name )

Definition at line 294 of file DatabasePDG.cc.

References mps_fire::i.

                                                  {
   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;
 }

bool DatabasePDG::GetPDGParticleStatusByIndex ( int index )

Definition at line 212 of file DatabasePDG.cc.

References diffTreeTool::index.

Referenced by GetDecayFilename().

                                                        {
   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];
 }

bool DatabasePDG::GetUseCharmParticles ( )

inline

Definition at line 90 of file DatabasePDG.h.

References fUseCharmParticles.

90 {return fUseCharmParticles;};

DatabasePDG::fUseCharmParticles

bool fUseCharmParticles

Definition: DatabasePDG.h:41

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

Definition at line 641 of file DatabasePDG.cc.

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

Referenced by GetMaximumMass().

                                                                            {
   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;
 }

bool DatabasePDG::LoadData ( )

Definition at line 66 of file DatabasePDG.cc.

Referenced by InitialState::InitialState().

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

bool DatabasePDG::LoadDecays ( )

private

Definition at line 142 of file DatabasePDG.cc.

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

                              {
   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;
 }

bool DatabasePDG::LoadParticles ( )

private

Definition at line 70 of file DatabasePDG.cc.

References EnergyCorrector::c, LogDebug, ResonanceBuilder::mass, dataset::name, lumiQueryAPI::q, alignCSCRings::s, spin(), and ApeEstimator_cff::width.

                                 {
   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;
 }

void DatabasePDG::SetDecayFilename ( char * filename )

Definition at line 62 of file DatabasePDG.cc.

                                                  {
   strcpy(fDecayFilename, filename);
 }

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

Definition at line 514 of file DatabasePDG.cc.

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

Referenced by InitialState::InitialState().

                                                      {
 
 
 
   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;
 }

void DatabasePDG::SetMaximumMass ( double value )

Definition at line 499 of file DatabasePDG.cc.

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

                                              {
   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;
 }

void DatabasePDG::SetMaximumWidth ( double value )

Definition at line 449 of file DatabasePDG.cc.

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

                                               {
   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;
 }

void DatabasePDG::SetMinimumMass ( double value )

Definition at line 484 of file DatabasePDG.cc.

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

                                              {
   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;
 }

void DatabasePDG::SetMinimumWidth ( double value )

Definition at line 434 of file DatabasePDG.cc.

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

                                               {
   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;
 }

void DatabasePDG::SetParticleFilename ( char * filename )

Definition at line 58 of file DatabasePDG.cc.

                                                     {
   strcpy(fParticleFilename, filename);
 }

void DatabasePDG::SetUseCharmParticles ( bool flag )

Definition at line 419 of file DatabasePDG.cc.

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

                                                 {
   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;
 }

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

Definition at line 464 of file DatabasePDG.cc.

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

Referenced by InitialState::InitialState().

                                                       {
   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;
 }

void DatabasePDG::SortParticles ( )

private

Definition at line 537 of file DatabasePDG.cc.

References mps_fire::i.

                                 {
 
 
   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;
 }

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

Definition at line 584 of file DatabasePDG.cc.

References runEdmFileComparison::found, mps_fire::i, kMaxParticles, and LogDebug.

                                                                        {
   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;
 }