GENIEGenerator
Loading...
Searching...
No Matches
Functions
gtestDecay.cxx File Reference
#include <ostream>
#include <iomanip>
#include <RVersion.h>
#include <TClonesArray.h>
#include <TParticlePDG.h>
#include <TIterator.h>
#include "Framework/Algorithm/Algorithm.h"
#include "Framework/Algorithm/AlgFactory.h"
#include "Framework/Conventions/Units.h"
#include "Framework/Conventions/Constants.h"
#include "Physics/Decay/Decayer.h"
#include "Physics/Decay/PythiaDecayer.h"
#include "Framework/GHEP/GHepStatus.h"
#include "Framework/GHEP/GHepParticle.h"
#include "Framework/GHEP/GHepRecord.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGLibrary.h"
Include dependency graph for gtestDecay.cxx:

Go to the source code of this file.

Functions

ostream & operator<< (ostream &stream, const TClonesArray *particle_list)
ostream & operator<< (ostream &stream, const GHepParticle *particle)
void TestPythiaTauDecays (void)
void Decay (const Decayer *decayer, int pdgc, double E, int ndecays)
int main (int, char **)

Function Documentation

◆ Decay()

void Decay ( const Decayer * decayer,
int pdgc,
double E,
int ndecays )

Definition at line 179 of file gtestDecay.cxx.

180{
181 DecayerInputs_t dinp;
182
183 TLorentzVector p4;
184 p4.SetE(E);
185 p4.SetTheta(0.);
186 p4.SetPhi(0.);
187
188 dinp.PdgCode = pdgc;
189 dinp.P4 = &p4;
190
191 PDGLibrary * pdglib = PDGLibrary::Instance();
192 TParticlePDG * pp = pdglib->Find(pdgc);
193 if(!pp) return;
194
195 LOG("test",pINFO)
196 << "Decaying a " << pp->GetName() << " with E = " << p4.Energy();
197
198 // Perform the decay a few times & print-out the decay products
199 for(int idec = 0; idec < ndecays; idec++) {
200
201 // Decay
202 TClonesArray * particle_list = decayer->Decay(dinp);
203
204 if(!particle_list) {
205 LOG("test",pWARN)
206 << "\n ** Decay nu.: " << idec << " ==> NULL particle list";
207 continue;
208 }
209
210 // Print-out
211 LOG("test",pINFO)
212 << "\n ** Decay nu.: " << idec
213 << " (weight = " << decayer->Weight() << ") : \n "
214 << particle_list;
215
216 // Clean-up
217 particle_list->Delete();
218 delete particle_list;
219 }//ndecays
220
221}
pINFO
#define pINFO
Definition Messenger.h:62
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE,...
Definition Messenger.h:96
pWARN
#define pWARN
Definition Messenger.h:60
genie::PDGLibrary
Singleton class to load & serve a TDatabasePDG.
Definition PDGLibrary.h:36
genie::PDGLibrary::Instance
static PDGLibrary * Instance(void)
Definition PDGLibrary.cxx:68
genie::PDGLibrary::Find
TParticlePDG * Find(int pdgc, bool must_exist=true)
Definition PDGLibrary.cxx:86

References genie::PDGLibrary::Find(), genie::PDGLibrary::Instance(), LOG, pINFO, and pWARN.

Referenced by TestPythiaTauDecays().

◆ main()

int main ( int ,
char **  )

Definition at line 56 of file gtestDecay.cxx.

57{
58 TestPythiaTauDecays();
59
60 return 0;
61}
TestPythiaTauDecays
void TestPythiaTauDecays(void)
Definition gtestDecay.cxx:63

References TestPythiaTauDecays().

◆ operator<<() [1/2]

ostream & operator<< ( ostream & stream,
const GHepParticle * particle )

Definition at line 250 of file gtestDecay.cxx.

251{
252 stream
253 << std::scientific << setprecision(6);
254
255 stream
256 << setfill(' ') << setw(10) << p->Name()
257 << setfill(' ') << setw(10) << p->Pdg()
258 << setfill(' ') << setw(10) << p->Status()
259 << setfill(' ') << setw(15) << p->Energy()
260 << setfill(' ') << setw(15) << p->Px()
261 << setfill(' ') << setw(15) << p->Py()
262 << setfill(' ') << setw(15) << p->Pz()
263 << setfill(' ') << setw(15) << p->Vt() /(units::mm)
264 << setfill(' ') << setw(15) << p->Vx() /(units::mm)
265 << setfill(' ') << setw(15) << p->Vy() /(units::mm)
266 << setfill(' ') << setw(15) << p->Vz() /(units::mm)
267 << endl;
268
269 return stream;
270}
genie::units::mm
static constexpr double mm
Definition Units.h:65

References genie::GHepParticle::Energy(), genie::units::mm, genie::GHepParticle::Name(), genie::GHepParticle::Pdg(), genie::GHepParticle::Px(), genie::GHepParticle::Py(), genie::GHepParticle::Pz(), genie::GHepParticle::Status(), genie::GHepParticle::Vt(), genie::GHepParticle::Vx(), genie::GHepParticle::Vy(), and genie::GHepParticle::Vz().

◆ operator<<() [2/2]

ostream & operator<< ( ostream & stream,
const TClonesArray * particle_list )

Definition at line 223 of file gtestDecay.cxx.

224{
225 GHepParticle * p = 0;
226 TObjArrayIter particle_iter(particle_list);
227
228
229 stream
230 << setfill(' ') << setw(10) << "name "
231 << setfill(' ') << setw(10) << "PDG"
232 << setfill(' ') << setw(10) << "Status"
233 << setfill(' ') << setw(15) << "E (GeV)"
234 << setfill(' ') << setw(15) << "Px (GeV/c)"
235 << setfill(' ') << setw(15) << "Py (GeV/c)"
236 << setfill(' ') << setw(15) << "Pz (GeV/c)"
237 << setfill(' ') << setw(15) << "t (mm/c)"
238 << setfill(' ') << setw(15) << "x (mm)"
239 << setfill(' ') << setw(15) << "y (mm)"
240 << setfill(' ') << setw(15) << "z (mm)"
241 << endl;
242
243 while( (p = (GHepParticle *) particle_iter.Next()) ) stream << p;
244
245 stream << setfill('-') << setw(100) << "|";
246
247 return stream;
248}
genie::GHepParticle
STDHEP-like event record entry that can fit a particle or a nucleus.
Definition GHepParticle.h:39

◆ TestPythiaTauDecays()

void TestPythiaTauDecays ( void )

Definition at line 63 of file gtestDecay.cxx.

64{
65 // Get the pythia decayer
66 LOG("test",pINFO)
67 << "Asking the AlgFactory for a genie::PythiaDecayer\\Default instance";
68 AlgFactory * algf = AlgFactory::Instance();
69 const Decayer * pdecayer =
70 dynamic_cast<const Decayer *> (
71 algf->GetAlgorithm("genie::PythiaDecayer","Default"));
72
73 // Decayer config print-out
74 LOG("test",pINFO) << "Algorithm name = " << pdecayer->Id().Name();
75 LOG("test",pINFO) << "Parameter set = " << pdecayer->Id().Config();
76
77 const Registry & conf_registry = pdecayer->GetConfig();
78 LOG("test", pINFO) << conf_registry;
79
80 double E = 10;
81 int ndec = 3;
82
83 // inhibit pi0 decay
84 pdecayer->InhibitDecay(kPdgPi0);
85
86 // decay tau-
87 Decay(pdecayer, kPdgTau, E, ndec);
88
89 // now inhibit all but the tau- --> nu_mu_bar + mu- + nu_tau decay channel
90 // (see $GENIE/src/contrib/misc/print_decay_channels.C)
91 // and perform some more decays
92 LOG("test",pINFO)
93 << "\n\n"
94 << " **** Inhibiting all but the `tau- --> nu_mu_bar + mu- + nu_tau' decay channel"
95 << "\n\n";
96
97 PDGLibrary * pdglib = PDGLibrary::Instance();
98
99 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(0));
100 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(2));
101 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(3));
102 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(4));
103 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(5));
104 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(6));
105 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(7));
106 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(8));
107 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(9));
108 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(10));
109 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(11));
110 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(12));
111 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(13));
112 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(14));
113 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(15));
114 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(16));
115 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(17));
116 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(18));
117 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(19));
118 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(20));
119 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(21));
120 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(22));
121 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(23));
122 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(24));
123 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(25));
124 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(26));
125 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(27));
126 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(28));
127 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(29));
128 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(30));
129 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(31));
130 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(32));
131 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(33));
132 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(34));
133 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(35));
134 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(36));
135 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(37));
136 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(38));
137 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(39));
138 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(40));
139 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(41));
140 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(42));
141 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(43));
142 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(44));
143 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(45));
144 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(46));
145 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(47));
146 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(48));
147 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(49));
148 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(50));
149 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(51));
150 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(52));
151 pdecayer->InhibitDecay(kPdgTau, pdglib->Find(kPdgTau)->DecayChannel(53));
152
153 // a few more decays
154 Decay(pdecayer, kPdgTau, E, ndec);
155
156 // restore all decay channels
157 LOG("test",pINFO)
158 << "\n\n"
159 << " **** Restoring all tau- decay channels"
160 << "\n\n";
161
162 pdecayer->UnInhibitDecay(kPdgTau);
163
164 // a few more decays
165 Decay(pdecayer, kPdgTau, E, ndec);
166
167 // now inhibit all decay channels and try to decay!
168 LOG("test",pINFO)
169 << "\n\n"
170 << " **** Inhibit all tau- decay channels"
171 << "\n\n";
172
173 pdecayer->InhibitDecay(kPdgTau);
174
175 // a few more decays
176 Decay(pdecayer, kPdgTau, E, ndec);
177}
genie::AlgFactory
The GENIE Algorithm Factory.
Definition AlgFactory.h:39
genie::AlgFactory::GetAlgorithm
const Algorithm * GetAlgorithm(const AlgId &algid)
Definition AlgFactory.cxx:75
genie::AlgFactory::Instance
static AlgFactory * Instance()
Definition AlgFactory.cxx:64
genie::AlgId::Config
string Config(void) const
Definition AlgId.h:45
genie::AlgId::Name
string Name(void) const
Definition AlgId.h:44
genie::Algorithm::GetConfig
virtual const Registry & GetConfig(void) const
Definition Algorithm.cxx:246
genie::Algorithm::Id
virtual const AlgId & Id(void) const
Get algorithm ID.
Definition Algorithm.h:98
genie::Decayer
Base class for decayer classes. Implements common configuration, allowing users to toggle on/off flag...
Definition Decayer.h:34
genie::Decayer::UnInhibitDecay
virtual void UnInhibitDecay(int pdgc, TDecayChannel *dc=0) const =0
genie::Decayer::InhibitDecay
virtual void InhibitDecay(int pdgc, TDecayChannel *dc=0) const =0
genie::Registry
A registry. Provides the container for algorithm configuration parameters.
Definition Registry.h:65
Decay
void Decay(const Decayer *decayer, int pdgc, double E, int ndecays)
Definition gtestDecay.cxx:179
genie::kPdgPi0
const int kPdgPi0
Definition PDGCodes.h:160
genie::kPdgTau
const int kPdgTau
Definition PDGCodes.h:39

References genie::AlgId::Config(), Decay(), genie::PDGLibrary::Find(), genie::AlgFactory::GetAlgorithm(), genie::Algorithm::GetConfig(), genie::Algorithm::Id(), genie::Decayer::InhibitDecay(), genie::AlgFactory::Instance(), genie::PDGLibrary::Instance(), genie::kPdgPi0, genie::kPdgTau, LOG, genie::AlgId::Name(), pINFO, and genie::Decayer::UnInhibitDecay().

Referenced by main().

Generated on for GENIEGenerator by doxygen 1.14.0