CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Static Private Attributes
CustomParticleFactory Class Reference

#include <CustomParticleFactory.h>

Public Member Functions

 CustomParticleFactory ()
 
const std::vector< G4ParticleDefinition * > & getCustomParticles ()
 
void loadCustomParticles (const std::string &filePath)
 
 ~CustomParticleFactory ()=default
 

Private Member Functions

void addCustomParticle (int pdgCode, double mass, const std::string &name)
 
G4DecayTable * getAntiDecayTable (int pdgId, G4DecayTable *theDecayTable)
 
G4DecayTable * getDecayTable (std::ifstream *configFile, int pdgId)
 
void getMassTable (std::ifstream *configFile)
 
std::string ToLower (std::string str)
 

Static Private Attributes

static bool loaded = false
 
static std::vector< G4ParticleDefinition * > m_particles
 

Detailed Description

Definition at line 12 of file CustomParticleFactory.h.

Constructor & Destructor Documentation

◆ CustomParticleFactory()

CustomParticleFactory::CustomParticleFactory ( )
explicit

Definition at line 29 of file CustomParticleFactory.cc.

29 {}

◆ ~CustomParticleFactory()

CustomParticleFactory::~CustomParticleFactory ( )
default

Member Function Documentation

◆ addCustomParticle()

void CustomParticleFactory::addCustomParticle ( int  pdgCode,
double  mass,
const std::string &  name 
)
private

Definition at line 92 of file CustomParticleFactory.cc.

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, CMSAntiSIMP::Definition(), CMSSIMP::Definition(), createfilelist::int, m_particles, EgHLTOffHistBins_cfi::mass, Skims_PA_cff::name, CustomPDGParser::s_charge(), CustomPDGParser::s_isDphoton(), CustomPDGParser::s_isgluinoHadron(), CustomPDGParser::s_isMesonino(), CustomPDGParser::s_isSbaryon(), CustomPDGParser::s_issbottomHadron(), CustomPDGParser::s_isSIMP(), CustomPDGParser::s_isSLepton(), CustomPDGParser::s_isstopHadron(), CustomPDGParser::s_spin(), CustomParticle::SetCloud(), CustomParticle::SetSpectator(), Validation_hcalonly_cfi::sign, TopDecayID::stable, str, AlCaHLTBitMon_QueryRunRegistry::string, and ApeEstimator_cff::width.

Referenced by getMassTable().

92  {
93  if (std::abs(pdgCode) % 100 < 14 && std::abs(pdgCode) / 1000000 == 0) {
94  edm::LogError("CustomParticleFactory::addCustomParticle")
95  << "Pdg code too low " << pdgCode << " " << std::abs(pdgCode) / 1000000;
96  return;
97  }
98 
99  if (CustomPDGParser::s_isSIMP(pdgCode)) {
100  CMSSIMP *simp = CMSSIMP::Definition(mass * GeV);
101  CMSAntiSIMP *antisimp = CMSAntiSIMP::Definition(mass * GeV);
102  m_particles.push_back(simp);
103  m_particles.push_back(antisimp);
104  return;
105  }
106 
108  G4String pType = "custom";
109  G4String pSubType = "";
110  G4double spectatormass = 0.0;
111  G4ParticleDefinition *spectator = nullptr;
113  if (CustomPDGParser::s_isgluinoHadron(pdgCode)) {
114  pType = "rhadron";
115  }
116  if (CustomPDGParser::s_isSLepton(pdgCode)) {
117  pType = "sLepton";
118  }
119  if (CustomPDGParser::s_isMesonino(pdgCode)) {
120  pType = "mesonino";
121  }
122  if (CustomPDGParser::s_isSbaryon(pdgCode)) {
123  pType = "sbaryon";
124  }
125 
126  double massGeV = mass * CLHEP::GeV;
127  double width = 0.0;
128  double charge = CLHEP::eplus * CustomPDGParser::s_charge(pdgCode);
129  if (name.compare(0, 4, "~HIP") == 0) {
130  if ((name.compare(0, 7, "~HIPbar") == 0)) {
131  std::string str = name.substr(7);
132  charge = CLHEP::eplus * atoi(str.c_str()) / 3.;
133  } else {
134  std::string str = name.substr(4);
135  charge = -CLHEP::eplus * atoi(str.c_str()) / 3.;
136  }
137  }
138  if (name.compare(0, 9, "anti_~HIP") == 0) {
139  if ((name.compare(0, 12, "anti_~HIPbar") == 0)) {
140  std::string str = name.substr(12);
141  charge = -CLHEP::eplus * atoi(str.c_str()) / 3.;
142  } else {
143  std::string str = name.substr(9);
144  charge = CLHEP::eplus * atoi(str.c_str()) / 3.;
145  }
146  }
147  int spin = (int)CustomPDGParser::s_spin(pdgCode) - 1;
148  int parity = +1;
149  int conjugation = 0;
150  int isospin = 0;
151  int isospinZ = 0;
152  int gParity = 0;
153  int lepton = 0; //FIXME:
154  int baryon = 1; //FIXME:
155  bool stable = true;
156  double lifetime = -1;
157 
158  if (CustomPDGParser::s_isDphoton(pdgCode)) {
159  pType = "darkpho";
160  spin = 2;
161  parity = -1;
162  conjugation = -1;
163  isospin = 0;
164  isospinZ = 0;
165  gParity = 0;
166  lepton = 0;
167  baryon = 0;
168  stable = true;
169  lifetime = -1;
170  }
171 
172  G4DecayTable *decaytable = nullptr;
173  G4ParticleTable *theParticleTable = G4ParticleTable::GetParticleTable();
174 
175  CustomParticle *particle = new CustomParticle(name,
176  massGeV,
177  width,
178  charge,
179  spin,
180  parity,
181  conjugation,
182  isospin,
183  isospinZ,
184  gParity,
185  pType,
186  lepton,
187  baryon,
188  pdgCode,
189  stable,
190  lifetime,
191  decaytable);
192 
193  if (pType == "rhadron" && name != "~g") {
194  G4String cloudname = name + "cloud";
195  G4String cloudtype = pType + "cloud";
196  spectator = theParticleTable->FindParticle(1000021);
197  spectatormass = spectator->GetPDGMass();
198  G4double cloudmass = massGeV - spectatormass;
199  CustomParticle *tmpParticle =
200  new CustomParticle(cloudname, cloudmass, 0.0, 0, 0, +1, 0, 0, 0, 0, cloudtype, 0, +1, 0, true, -1.0, nullptr);
201  particle->SetCloud(tmpParticle);
202  particle->SetSpectator(spectator);
203 
204  edm::LogVerbatim("SimG4CoreCustomPhysics")
205  << "CustomParticleFactory: " << name << " being assigned spectator" << spectator->GetParticleName()
206  << " and cloud " << cloudname << "\n"
207  << " Masses: " << mass << " Gev, " << spectatormass / CLHEP::GeV << " GeV and "
208  << cloudmass / CLHEP::GeV << " GeV.";
209  } else if (pType == "mesonino" || pType == "sbaryon") {
210  int sign = 1;
211  if (pdgCode < 0)
212  sign = -1;
213 
214  G4String cloudname = name + "cloud";
215  G4String cloudtype = pType + "cloud";
216  if (CustomPDGParser::s_isstopHadron(pdgCode)) {
217  spectator = theParticleTable->FindParticle(1000006 * sign);
218  } else {
219  if (CustomPDGParser::s_issbottomHadron(pdgCode)) {
220  spectator = theParticleTable->FindParticle(1000005 * sign);
221  } else {
222  spectator = nullptr;
223  edm::LogError("SimG4CoreCustomPhysics") << "CustomParticleFactory: Cannot find spectator parton";
224  }
225  }
226  if (nullptr != spectator) {
227  spectatormass = spectator->GetPDGMass();
228  }
229  G4double cloudmass = massGeV - spectatormass;
230  CustomParticle *tmpParticle =
231  new CustomParticle(cloudname, cloudmass, 0.0, 0, 0, +1, 0, 0, 0, 0, cloudtype, 0, +1, 0, true, -1.0, nullptr);
232  particle->SetCloud(tmpParticle);
233  particle->SetSpectator(spectator);
234 
235  if (nullptr != spectator)
236  edm::LogVerbatim("SimG4CoreCustomPhysics")
237  << "CustomParticleFactory: " << name << " being assigned spectator" << spectator->GetParticleName()
238  << " and cloud " << cloudname << "\n"
239  << " Masses: " << mass << " Gev, " << spectatormass / CLHEP::GeV << " GeV and "
240  << cloudmass / CLHEP::GeV << " GeV.";
241  } else {
242  particle->SetCloud(nullptr);
243  particle->SetSpectator(nullptr);
244  }
245  m_particles.push_back(particle);
246 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
CustomPDGParser::s_isDphoton
static bool s_isDphoton(int pdg)
Definition: CustomPDGParser.cc:35
CustomPDGParser::s_charge
static double s_charge(int pdg)
Definition: CustomPDGParser.cc:65
TopDecayID::stable
static const int stable
Definition: TopGenEvent.h:10
CustomPDGParser::s_isSbaryon
static bool s_isSbaryon(int pdg)
Definition: CustomPDGParser.cc:50
Validation_hcalonly_cfi.sign
sign
Definition: Validation_hcalonly_cfi.py:32
CustomPDGParser::s_isgluinoHadron
static bool s_isgluinoHadron(int pdg)
Definition: CustomPDGParser.cc:5
CustomParticle::SetCloud
void SetCloud(G4ParticleDefinition *theCloud)
Definition: CustomParticle.h:43
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
CustomParticleFactory::m_particles
static std::vector< G4ParticleDefinition * > m_particles
Definition: CustomParticleFactory.h:28
CMSSIMP
Definition: CMSSIMP.h:8
CustomPDGParser::s_isMesonino
static bool s_isMesonino(int pdg)
Definition: CustomPDGParser.cc:45
CustomPDGParser::s_isSIMP
static bool s_isSIMP(int pdg)
Definition: CustomPDGParser.cc:60
CustomPDGParser::s_isSLepton
static bool s_isSLepton(int pdg)
Definition: CustomPDGParser.cc:20
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
CustomParticle::SetSpectator
void SetSpectator(G4ParticleDefinition *theSpectator)
Definition: CustomParticle.h:45
CustomPDGParser::s_spin
static double s_spin(int pdg)
Definition: CustomPDGParser.cc:140
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
EgHLTOffHistBins_cfi.mass
mass
Definition: EgHLTOffHistBins_cfi.py:34
CMSAntiSIMP::Definition
static CMSAntiSIMP * Definition(double mass)
Definition: CMSAntiSIMP.cc:12
ApeEstimator_cff.width
width
Definition: ApeEstimator_cff.py:24
CustomPDGParser::s_issbottomHadron
static bool s_issbottomHadron(int pdg)
Definition: CustomPDGParser.cc:15
CustomPDGParser::s_isstopHadron
static bool s_isstopHadron(int pdg)
Definition: CustomPDGParser.cc:10
str
#define str(s)
Definition: TestProcessor.cc:54
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
CMSSIMP::Definition
static CMSSIMP * Definition(double mass)
Definition: CMSSIMP.cc:12

◆ getAntiDecayTable()

G4DecayTable * CustomParticleFactory::getAntiDecayTable ( int  pdgId,
G4DecayTable *  theDecayTable 
)
private

Definition at line 399 of file CustomParticleFactory.cc.

References mps_fire::i, dqmiolumiharvest::j, SiStripPI::min, Skims_PA_cff::name, EgammaValidation_cff::pdgId, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by loadCustomParticles().

399  {
400  std::string name[4] = {""};
401  G4ParticleTable *theParticleTable = G4ParticleTable::GetParticleTable();
402 
403  std::string parentName = theParticleTable->FindParticle(-pdgId)->GetParticleName();
404  G4DecayTable *decaytable = new G4DecayTable();
405 
406  for (int i = 0; i < theDecayTable->entries(); ++i) {
407  G4VDecayChannel *theDecayChannel = theDecayTable->GetDecayChannel(i);
408  int nd = std::min(4, theDecayChannel->GetNumberOfDaughters());
409  for (int j = 0; j < nd; ++j) {
410  int id = theDecayChannel->GetDaughter(j)->GetAntiPDGEncoding();
411  G4ParticleDefinition *part = theParticleTable->FindParticle(id);
412  if (nullptr == part) {
413  edm::LogWarning("SimG4CoreCustomPhysics")
414  << "CustomParticleFactory: antiparticle with PDG code" << id << " not found!";
415  continue;
416  }
417  name[j] = part->GetParticleName();
418  }
419  G4PhaseSpaceDecayChannel *aDecayChannel =
420  new G4PhaseSpaceDecayChannel(parentName, theDecayChannel->GetBR(), nd, name[0], name[1], name[2], name[3]);
421  decaytable->Insert(aDecayChannel);
422  edm::LogVerbatim("SimG4CoreCustomPhysics")
423  << "CustomParticleFactory: inserted decay channel "
424  << " for pdgID= " << -pdgId << " " << parentName << " BR= " << theDecayChannel->GetBR()
425  << " Daughters: " << name[0] << " " << name[1] << " " << name[2] << " " << name[3];
426  }
427  return decaytable;
428 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
mps_fire.i
i
Definition: mps_fire.py:428
EgammaValidation_cff.pdgId
pdgId
Definition: EgammaValidation_cff.py:117
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
part
part
Definition: HCALResponse.h:20
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17

◆ getCustomParticles()

const std::vector< G4ParticleDefinition * > & CustomParticleFactory::getCustomParticles ( )

Definition at line 437 of file CustomParticleFactory.cc.

References m_particles.

437 { return m_particles; }
CustomParticleFactory::m_particles
static std::vector< G4ParticleDefinition * > m_particles
Definition: CustomParticleFactory.h:28

◆ getDecayTable()

G4DecayTable * CustomParticleFactory::getDecayTable ( std::ifstream *  configFile,
int  pdgId 
)
private

Definition at line 338 of file CustomParticleFactory.cc.

References beamvalidation::br, beamspotdip_dqm_sourceclient-live_cfg::configFile, mps_fire::i, mps_splice::line, LogDebug, Skims_PA_cff::name, EgammaValidation_cff::pdgId, AlCaHLTBitMon_QueryRunRegistry::string, and ToLower().

Referenced by loadCustomParticles().

338  {
339  double br;
340  int nDaughters;
341  int pdg[4] = {0};
342  std::string line;
343 
344  G4ParticleTable *theParticleTable = G4ParticleTable::GetParticleTable();
345 
346  std::string parentName = theParticleTable->FindParticle(pdgId)->GetParticleName();
347  G4DecayTable *decaytable = new G4DecayTable();
348 
349  while (getline(*configFile, line)) {
350  line.erase(0, line.find_first_not_of(" \t")); // remove leading whitespace
351  if (line.length() == 0)
352  continue; // skip blank lines
353  if (line.at(0) == '#' && ToLower(line).find("br") < line.npos && ToLower(line).find("nda") < line.npos)
354  continue; // skip a comment of the form: # BR NDA ID1 ID2
355  if (line.at(0) == '#' || ToLower(line).find("block") < line.npos) {
356  edm::LogInfo("SimG4CoreCustomPhysics") << "CustomParticleFactory: Finished the Decay Table ";
357  break;
358  }
359 
360  std::stringstream sstr(line);
361  sstr >> br >> nDaughters; // assume SLHA format, e.g.: 1.49435135E-01 2 -15 16 # BR(H+ -> tau+ nu_tau)
362  LogDebug("SimG4CoreCustomPhysics") << "CustomParticleFactory: Branching Ratio: " << br
363  << ", Number of Daughters: " << nDaughters;
364  if (nDaughters > 4) {
365  edm::LogError("SimG4CoreCustomPhysics")
366  << "CustomParticleFactory: Number of daughters is too large (max = 4): " << nDaughters
367  << " for pdgId: " << pdgId;
368  break;
369  }
370  if (br <= 0.0) {
371  edm::LogError("SimG4CoreCustomPhysics")
372  << "CustomParticleFactory: Branching ratio is " << br << " for pdgId: " << pdgId;
373  break;
374  }
375  std::string name[4] = {""};
376  for (int i = 0; i < nDaughters; ++i) {
377  sstr >> pdg[i];
378  LogDebug("SimG4CoreCustomPhysics") << "CustomParticleFactory: Daughter ID " << pdg[i];
379  const G4ParticleDefinition *part = theParticleTable->FindParticle(pdg[i]);
380  if (!part) {
381  edm::LogWarning("SimG4CoreCustomPhysics")
382  << "CustomParticleFactory: particle with PDG code" << pdg[i] << " not found!";
383  continue;
384  }
385  name[i] = part->GetParticleName();
386  }
388  G4PhaseSpaceDecayChannel *aDecayChannel =
389  new G4PhaseSpaceDecayChannel(parentName, br, nDaughters, name[0], name[1], name[2], name[3]);
390  decaytable->Insert(aDecayChannel);
391  edm::LogVerbatim("SimG4CoreCustomPhysics")
392  << "CustomParticleFactory: inserted decay channel "
393  << " for pdgID= " << pdgId << " " << parentName << " BR= " << br << " Daughters: " << name[0] << " "
394  << name[1] << " " << name[2] << " " << name[3];
395  }
396  return decaytable;
397 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
mps_fire.i
i
Definition: mps_fire.py:428
EgammaValidation_cff.pdgId
pdgId
Definition: EgammaValidation_cff.py:117
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
beamspotdip_dqm_sourceclient-live_cfg.configFile
configFile
Definition: beamspotdip_dqm_sourceclient-live_cfg.py:8
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
mps_splice.line
line
Definition: mps_splice.py:76
CustomParticleFactory::ToLower
std::string ToLower(std::string str)
Definition: CustomParticleFactory.cc:430
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
part
part
Definition: HCALResponse.h:20
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233

◆ getMassTable()

void CustomParticleFactory::getMassTable ( std::ifstream *  configFile)
private

Definition at line 248 of file CustomParticleFactory.cc.

References funct::abs(), addCustomParticle(), beamspotdip_dqm_sourceclient-live_cfg::configFile, mps_splice::line, EgHLTOffHistBins_cfi::mass, SiStripPI::max, Skims_PA_cff::name, EgammaValidation_cff::pdgId, CustomPDGParser::s_isgluinoHadron(), CustomPDGParser::s_isSIMP(), CustomPDGParser::s_isstopHadron(), Validation_hcalonly_cfi::sign, AlCaHLTBitMon_QueryRunRegistry::string, createJobs::tmp, and ToLower().

Referenced by loadCustomParticles().

248  {
249  int pdgId;
250  double mass;
251  std::string name, tmp;
252  std::string line;
253  G4ParticleTable *theParticleTable = G4ParticleTable::GetParticleTable();
254 
255  // This should be compatible IMO to SLHA
256  while (getline(*configFile, line)) {
257  line.erase(0, line.find_first_not_of(" \t")); // remove leading whitespace
258  if (line.length() == 0 || line.at(0) == '#')
259  continue; // skip blank lines and comments
260  if (ToLower(line).find("block") < line.npos) {
261  edm::LogInfo("SimG4CoreCustomPhysics") << "CustomParticleFactory: Finished the Mass Table ";
262  break;
263  }
264  std::stringstream sstr(line);
265  sstr >> pdgId >> mass >> tmp >> name; // Assume SLHA format, e.g.: 1000001 5.68441109E+02 # ~d_L
266 
267  mass = std::max(mass, 0.0);
268  if (theParticleTable->FindParticle(pdgId)) {
269  continue;
270  }
271 
272  edm::LogVerbatim("SimG4CoreCustomPhysics")
273  << "CustomParticleFactory: Calling addCustomParticle for pdgId: " << pdgId << ", mass " << mass << " GeV "
274  << name << ", isgluinoHadron: " << CustomPDGParser::s_isgluinoHadron(pdgId)
275  << ", isstopHadron: " << CustomPDGParser::s_isstopHadron(pdgId)
276  << ", isSIMP: " << CustomPDGParser::s_isSIMP(pdgId);
277  addCustomParticle(pdgId, mass, name);
278 
280  int pdgIdPartner = pdgId % 100;
281  G4ParticleDefinition *aParticle = theParticleTable->FindParticle(pdgIdPartner);
282  //Add antiparticles for SUSY particles only, not for rHadrons.
283  if (nullptr != aParticle) {
284  edm::LogVerbatim("SimG4CoreCustomPhysics")
285  << "CustomParticleFactory: Found partner particle for "
286  << " pdgId= " << pdgId << ", pdgIdPartner= " << pdgIdPartner << " " << aParticle->GetParticleName();
287  }
288 
289  if (nullptr != aParticle && !CustomPDGParser::s_isgluinoHadron(pdgId) && !CustomPDGParser::s_isstopHadron(pdgId) &&
290  !CustomPDGParser::s_isSIMP(pdgId) && pdgId != 1000006 && pdgId != -1000006 && pdgId != 25 && pdgId != 35 &&
291  pdgId != 36 && pdgId != 37) {
292  int sign = aParticle->GetAntiPDGEncoding() / pdgIdPartner;
293  edm::LogVerbatim("SimG4CoreCustomPhysics")
294  << "CustomParticleFactory: For " << aParticle->GetParticleName() << " pdg= " << pdgIdPartner
295  << ", GetAntiPDGEncoding() " << aParticle->GetAntiPDGEncoding() << " sign= " << sign;
296 
297  if (std::abs(sign) != 1) {
298  edm::LogVerbatim("SimG4CoreCustomPhysics")
299  << "CustomParticleFactory: sign= " << sign << " a: " << aParticle->GetAntiPDGEncoding()
300  << " b: " << pdgIdPartner;
301  aParticle->DumpTable();
302  }
303  if (sign == -1 && pdgId != 25 && pdgId != 35 && pdgId != 36 && pdgId != 37 && pdgId != 1000039) {
304  tmp = "anti_" + name;
305  if (nullptr != theParticleTable->FindParticle(-pdgId)) {
306  edm::LogVerbatim("SimG4CoreCustomPhysics")
307  << "CustomParticleFactory: Calling addCustomParticle for antiparticle with pdgId: " << -pdgId << ", mass "
308  << mass << " GeV, name " << tmp;
309  addCustomParticle(-pdgId, mass, tmp);
310  theParticleTable->FindParticle(pdgId)->SetAntiPDGEncoding(-pdgId);
311  }
312  } else
313  theParticleTable->FindParticle(pdgId)->SetAntiPDGEncoding(pdgId);
314  }
315 
316  if (pdgId == 1000039)
317  theParticleTable->FindParticle(pdgId)->SetAntiPDGEncoding(pdgId); // gravitino
318  if (pdgId == 1000024 || pdgId == 1000037 || pdgId == 37) {
319  tmp = "anti_" + name;
320  if (!theParticleTable->FindParticle(-pdgId)) {
321  edm::LogVerbatim("SimG4CoreCustomPhysics")
322  << "CustomParticleFactory: Calling addCustomParticle for antiparticle (2) with pdgId: " << -pdgId
323  << ", mass " << mass << " GeV, name " << tmp;
324  addCustomParticle(-pdgId, mass, tmp);
325  theParticleTable->FindParticle(pdgId)->SetAntiPDGEncoding(-pdgId);
326  }
327  }
328  if (pdgId == 9000006) {
329  tmp = name + "bar";
330  edm::LogVerbatim("CustomPhysics") << "Calling addCustomParticle for antiparticle with pdgId: " << -pdgId
331  << ", mass " << mass << ", name " << tmp;
332  addCustomParticle(-pdgId, mass, tmp);
333  theParticleTable->FindParticle(pdgId)->SetAntiPDGEncoding(-pdgId);
334  }
335  }
336 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
Validation_hcalonly_cfi.sign
sign
Definition: Validation_hcalonly_cfi.py:32
CustomPDGParser::s_isgluinoHadron
static bool s_isgluinoHadron(int pdg)
Definition: CustomPDGParser.cc:5
EgammaValidation_cff.pdgId
pdgId
Definition: EgammaValidation_cff.py:117
beamspotdip_dqm_sourceclient-live_cfg.configFile
configFile
Definition: beamspotdip_dqm_sourceclient-live_cfg.py:8
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
CustomPDGParser::s_isSIMP
static bool s_isSIMP(int pdg)
Definition: CustomPDGParser.cc:60
mps_splice.line
line
Definition: mps_splice.py:76
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
CustomParticleFactory::ToLower
std::string ToLower(std::string str)
Definition: CustomParticleFactory.cc:430
EgHLTOffHistBins_cfi.mass
mass
Definition: EgHLTOffHistBins_cfi.py:34
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
CustomParticleFactory::addCustomParticle
void addCustomParticle(int pdgCode, double mass, const std::string &name)
Definition: CustomParticleFactory.cc:92
CustomPDGParser::s_isstopHadron
static bool s_isstopHadron(int pdg)
Definition: CustomPDGParser.cc:10
createJobs.tmp
tmp
align.sh
Definition: createJobs.py:716
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17

◆ loadCustomParticles()

void CustomParticleFactory::loadCustomParticles ( const std::string &  filePath)

Definition at line 31 of file CustomParticleFactory.cc.

References beamspotdip_dqm_sourceclient-live_cfg::configFile, MillePedeFileConverter_cfg::e, dt4ml_dqm_sourceclient-live_cfg::filePath, getAntiDecayTable(), getDecayTable(), getMassTable(), mps_splice::line, loaded, EgammaValidation_cff::pdgId, alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, ToLower(), and ApeEstimator_cff::width.

31  {
32  if (loaded) {
33  return;
34  }
35 #ifdef G4MULTITHREADED
36  G4MUTEXLOCK(&customParticleFactoryMutex);
37  if (loaded) {
38  return;
39  }
40 #endif
41 
42  // loading once
43  loaded = true;
44  std::ifstream configFile(filePath.c_str());
45 
46  std::string line;
47  edm::LogVerbatim("SimG4CoreCustomPhysics")
48  << "CustomParticleFactory: Reading Custom Particle and G4DecayTable from \n"
49  << filePath;
50  // This should be compatible IMO to SLHA
51  G4ParticleTable *theParticleTable = G4ParticleTable::GetParticleTable();
52  while (getline(configFile, line)) {
53  line.erase(0, line.find_first_not_of(" \t")); // Remove leading whitespace.
54  if (line.length() == 0 || line.at(0) == '#') {
55  continue;
56  } // Skip blank lines and comments.
57  // The mass table begins with a line containing "BLOCK MASS"
58  if (ToLower(line).find("block") < line.npos && ToLower(line).find("mass") < line.npos) {
59  edm::LogVerbatim("SimG4CoreCustomPhysics") << "CustomParticleFactory: Retrieving mass table.";
60  getMassTable(&configFile);
61  }
62  if (line.find("DECAY") < line.npos) {
63  int pdgId;
64  double width;
65  std::string tmpString;
66  std::stringstream lineStream(line);
67  lineStream >> tmpString >> pdgId >> width;
68  // assume SLHA format, e.g.: DECAY 1000021 5.50675438E+00 # gluino decays
69  edm::LogVerbatim("SimG4CoreCustomPhysics")
70  << "CustomParticleFactory: entry to G4DecayTable: pdgID, width " << pdgId << ", " << width;
71  G4ParticleDefinition *aParticle = theParticleTable->FindParticle(pdgId);
72  if (nullptr == aParticle || width == 0.0) {
73  continue;
74  }
75  G4DecayTable *aDecayTable = getDecayTable(&configFile, pdgId);
76  aParticle->SetDecayTable(aDecayTable);
77  aParticle->SetPDGStable(false);
78  aParticle->SetPDGLifeTime(1.0 / (width * CLHEP::GeV) * 6.582122e-22 * CLHEP::MeV * CLHEP::s);
79  G4ParticleDefinition *aAntiParticle = theParticleTable->FindAntiParticle(pdgId);
80  if (nullptr != aAntiParticle && aAntiParticle->GetPDGEncoding() != pdgId) {
81  aAntiParticle->SetDecayTable(getAntiDecayTable(pdgId, aDecayTable));
82  aAntiParticle->SetPDGStable(false);
83  aAntiParticle->SetPDGLifeTime(1.0 / (width * CLHEP::GeV) * 6.582122e-22 * CLHEP::MeV * CLHEP::s);
84  }
85  }
86  }
87 #ifdef G4MULTITHREADED
88  G4MUTEXUNLOCK(&customParticleFactoryMutex);
89 #endif
90 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
EgammaValidation_cff.pdgId
pdgId
Definition: EgammaValidation_cff.py:117
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
beamspotdip_dqm_sourceclient-live_cfg.configFile
configFile
Definition: beamspotdip_dqm_sourceclient-live_cfg.py:8
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
mps_splice.line
line
Definition: mps_splice.py:76
dt4ml_dqm_sourceclient-live_cfg.filePath
filePath
CUSTOMIZE FOR ML.
Definition: dt4ml_dqm_sourceclient-live_cfg.py:46
CustomParticleFactory::ToLower
std::string ToLower(std::string str)
Definition: CustomParticleFactory.cc:430
CustomParticleFactory::getMassTable
void getMassTable(std::ifstream *configFile)
Definition: CustomParticleFactory.cc:248
CustomParticleFactory::getDecayTable
G4DecayTable * getDecayTable(std::ifstream *configFile, int pdgId)
Definition: CustomParticleFactory.cc:338
ApeEstimator_cff.width
width
Definition: ApeEstimator_cff.py:24
CustomParticleFactory::getAntiDecayTable
G4DecayTable * getAntiDecayTable(int pdgId, G4DecayTable *theDecayTable)
Definition: CustomParticleFactory.cc:399

◆ ToLower()

std::string CustomParticleFactory::ToLower ( std::string  str)
private

Definition at line 430 of file CustomParticleFactory.cc.

References mps_fire::i, and str.

Referenced by getDecayTable(), getMassTable(), and loadCustomParticles().

430  {
431  std::locale loc;
432  for (std::string::size_type i = 0; i < str.length(); ++i)
433  str.at(i) = std::tolower(str.at(i), loc);
434  return str;
435 }
mps_fire.i
i
Definition: mps_fire.py:428
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:18
str
#define str(s)
Definition: TestProcessor.cc:54

Member Data Documentation

◆ loaded

bool CustomParticleFactory::loaded = false
staticprivate

Definition at line 27 of file CustomParticleFactory.h.

Referenced by loadCustomParticles().

◆ m_particles

std::vector< G4ParticleDefinition * > CustomParticleFactory::m_particles
staticprivate

Definition at line 28 of file CustomParticleFactory.h.

Referenced by addCustomParticle(), and getCustomParticles().

Generated for CMSSW Reference Manual by doxygen 1.8.14