Public Member Functions | Private Member Functions | Private Attributes | Friends | List of all members
genie::HNIntranuke2014 Class Reference

#include <HNIntranuke2014.h>

Inheritance diagram for genie::HNIntranuke2014:
genie::Intranuke2014 genie::EventRecordVisitorI genie::Algorithm

Public Member Functions

 HNIntranuke2014 ()
 
 HNIntranuke2014 (string config)
 
 ~HNIntranuke2014 ()
 
void ProcessEventRecord (GHepRecord *event_rec) const
 
- Public Member Functions inherited from genie::Intranuke2014
 Intranuke2014 ()
 
 Intranuke2014 (string name)
 
 Intranuke2014 (string name, string config)
 
 ~Intranuke2014 ()
 
void Configure (const Registry &config)
 Configure the algorithm. More...
 
void Configure (string param_set)
 Configure the algorithm. More...
 
- Public Member Functions inherited from genie::EventRecordVisitorI
virtual ~EventRecordVisitorI ()
 
- Public Member Functions inherited from genie::Algorithm
virtual ~Algorithm ()
 
virtual void FindConfig (void)
 Lookup configuration from the config pool. More...
 
virtual const Registry & GetConfig (void) const
 Get configuration registry. More...
 
Registry * GetOwnedConfig (void)
 Get a writeable version of an owned configuration Registry. More...
 
virtual const AlgId & Id (void) const
 Get algorithm ID. More...
 
virtual AlgStatus_t GetStatus (void) const
 Get algorithm status. More...
 
virtual bool AllowReconfig (void) const
 
virtual AlgCmp_t Compare (const Algorithm *alg) const
 Compare with input algorithm. More...
 
virtual void SetId (const AlgId &id)
 Set algorithm ID. More...
 
virtual void SetId (string name, string config)
 
const Algorithm * SubAlg (const RgKey &registry_key) const
 
void AdoptConfig (void)
 
void AdoptSubstructure (void)
 
virtual void Print (ostream &stream) const
 Print algorithm info. More...
 

Private Member Functions

void LoadConfig (void)
 
void SimulateHadronicFinalState (GHepRecord *ev, GHepParticle *p) const
 
INukeFateHN_t HadronFateHN (const GHepParticle *p) const
 
double FateWeight (int pdgc, INukeFateHN_t fate) const
 
void ElasHN (GHepRecord *ev, GHepParticle *p, INukeFateHN_t fate) const
 
void AbsorbHN (GHepRecord *ev, GHepParticle *p, INukeFateHN_t fate) const
 
void InelasticHN (GHepRecord *ev, GHepParticle *p) const
 
void GammaInelasticHN (GHepRecord *ev, GHepParticle *p, INukeFateHN_t fate) const
 
bool HandleCompoundNucleusHN (GHepRecord *ev, GHepParticle *p) const
 
int HandleCompoundNucleus (GHepRecord *ev, GHepParticle *p, int mom) const
 

Private Attributes

int nuclA
 value of A for the target nucleus in hA mode More...
 
double fNucQEFac
 

Friends

class IntranukeTester
 

Additional Inherited Members

- Protected Member Functions inherited from genie::Intranuke2014
void TransportHadrons (GHepRecord *ev) const
 
void GenerateVertex (GHepRecord *ev) const
 
bool NeedsRescattering (const GHepParticle *p) const
 
bool CanRescatter (const GHepParticle *p) const
 
bool IsInNucleus (const GHepParticle *p) const
 
void SetTrackingRadius (const GHepParticle *p) const
 
double GenerateStep (GHepRecord *ev, GHepParticle *p) const
 
- Protected Member Functions inherited from genie::EventRecordVisitorI
 EventRecordVisitorI ()
 
 EventRecordVisitorI (string name)
 
 EventRecordVisitorI (string name, string config)
 
- Protected Member Functions inherited from genie::Algorithm
 Algorithm ()
 
 Algorithm (string name)
 
 Algorithm (string name, string config)
 
void Initialize (void)
 
void DeleteConfig (void)
 
void DeleteSubstructure (void)
 
- Protected Attributes inherited from genie::Intranuke2014
double fTrackingRadius
 tracking radius for the nucleus in the current event More...
 
TGenPhaseSpace fGenPhaseSpace
 a phase space generator More...
 
INukeHadroData2014 * fHadroData2014
 a collection of h+N,h+A data & calculations More...
 
AlgFactory * fAlgf
 algorithm factory instance More...
 
const NuclearModelI * fNuclmodel
 nuclear model used to generate fermi momentum More...
 
int fRemnA
 remnant nucleus A More...
 
int fRemnZ
 remnant nucleus Z More...
 
TLorentzVector fRemnP4
 P4 of remnant system. More...
 
GEvGenMode_t fGMode
 event generation mode (lepton+A, hadron+A, ...) More...
 
double fR0
 effective nuclear size param More...
 
double fNR
 param multiplying the nuclear radius, determining how far to track hadrons beyond the "nuclear boundary" More...
 
double fNucRmvE
 binding energy to subtract from cascade nucleons More...
 
double fDelRPion
 factor by which Pion Compton wavelength gets multiplied to become nuclear size enhancement More...
 
double fDelRNucleon
 factor by which Nucleon Compton wavelength gets multiplied to become nuclear size enhancement More...
 
double fHadStep
 step size for intranuclear hadron transport More...
 
double fNucAbsFac
 absorption xsec correction factor (hN Mode) More...
 
double fNucCEXFac
 charge exchange xsec correction factor (hN Mode) More...
 
double fEPreEq
 threshold for pre-equilibrium reaction More...
 
double fFermiFac
 testing parameter to modify fermi momentum More...
 
double fFreeStep
 produced particle free stem, in fm More...
 
double fFermiMomentum
 whether or not particle collision is pauli blocked More...
 
bool fDoFermi
 whether or not to do fermi mom. More...
 
bool fDoMassDiff
 whether or not to do mass diff. mode More...
 
bool fDoCompoundNucleus
 whether or not to do compound nucleus considerations More...
 
- Protected Attributes inherited from genie::Algorithm
bool fAllowReconfig
 
bool fOwnsConfig
 true if it owns its config. registry More...
 
bool fOwnsSubstruc
 true if it owns its substructure (sub-algs,...) More...
 
AlgId fID
 algorithm name and configuration set More...
 
Registry * fConfig
 config. (either owned or pointing to config pool) More...
 
AlgStatus_t fStatus
 algorithm execution status More...
 
AlgMap * fOwnedSubAlgMp
 local pool for owned sub-algs (taken out of the factory pool) More...
 

Detailed Description

Definition at line 51 of file HNIntranuke2014.h.

Constructor & Destructor Documentation

HNIntranuke2014::HNIntranuke2014 ( )

Definition at line 87 of file HNIntranuke2014.cxx.

87  :
88 Intranuke2014("genie::HNIntranuke2014")
89 {
90 
91 }
HNIntranuke2014::HNIntranuke2014 ( string  config)

Definition at line 93 of file HNIntranuke2014.cxx.

93  :
94 Intranuke2014("genie::HNIntranuke2014",config)
95 {
96 
97 }
generate_config.config
config
Definition: generate_config.py:3
HNIntranuke2014::~HNIntranuke2014 ( )

Definition at line 99 of file HNIntranuke2014.cxx.

100 {
101 
102 }

Member Function Documentation

void HNIntranuke2014::AbsorbHN ( GHepRecord *  ev,
GHepParticle *  p,
INukeFateHN_t  fate 
) const
private

Definition at line 335 of file HNIntranuke2014.cxx.

337 {
338  // handles pi+d->2p, pi-d->nn, pi0 d->pn absorbtion, all using pi+d values
339 
340  int pdgc = p->Pdg();
341 
342 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
343  LOG("HNIntranuke2014", pDEBUG)
344  << "AbsorbHN() is invoked for a : " << p->Name()
345  << " whose fate is : " << INukeHadroFates::AsString(fate);
346 #endif
347 
348  // check fate
349  if(fate!=kIHNFtAbs)
350  {
351  LOG("HNIntranuke2014", pWARN)
352  << "AbsorbHN() cannot handle fate: " << INukeHadroFates::AsString(fate);
353  return;
354  }
355 
356  // random number generator
357  RandomGen * rnd = RandomGen::Instance();
358 
359  // Notes on the kinematics
360  // -- Simple variables are used for efficiency
361  // -- Variables are numbered according to particle
362  // -- -- #1 -> incoming particle
363  // -- -- #2 -> target (here, 2_1 and 2_2 for individual particles)
364  // -- -- #3 -> scattered incoming (Particle tracked in hA mode)
365  // -- -- #4 -> other scattered particle
366  // -- Suffix "L" is for lab frame, suffix "CM" is for center of mass frame
367  // -- Subscript "z" is for parallel component, "t" is for transverse
368 
369  int pcode, t1code, t2code, scode, s2code; // particles
370  double M1, M2_1, M2_2, M3, M4; // rest energies, in GeV
371  double E1L, P1L, E2L, P2L, E3L, P3L, E4L, P4L;
372  double P1zL, P2zL;
373  double beta, gm; // speed and gamma for CM frame in lab
374  double Et, E2CM;
375  double C3CM, S3CM; // cos and sin of scattering angle
376  double Theta1, Theta2, theta5;
377  double PHI3; // transverse scattering angle
378  double E1CM, E3CM, E4CM, P3CM;
379  double P3zL, P3tL, P4zL, P4tL;
380  double E2_1L, E2_2L;
381  TVector3 tP2_1L, tP2_2L, tP1L, tP2L, tPtot, P1zCM, P2zCM;
382  TVector3 tP3L, tP4L;
383  TVector3 bDir, tTrans, tbeta, tVect;
384 
385  // Library instance for reference
386  PDGLibrary * pLib = PDGLibrary::Instance();
387 
388  // Handle fermi target
389  Target target(ev->TargetNucleus()->Pdg());
390 
391  // Target should be a deuteron, but for now
392  // handling it as seperate nucleons
393  if(pdgc==211) // pi-plus
394  {
395  pcode = 211;
396  t1code = 2212; // proton
397  t2code = 2112; // neutron
398  scode = 2212;
399  s2code = 2212;
400  }
401  else if(pdgc==-211) // pi-minus
402  {
403  pcode = -211;
404  t1code = 2212;
405  t2code = 2112;
406  scode = 2112;
407  s2code = 2112;
408  }
409  else if(pdgc==111) // pi-zero
410  {
411  pcode = 111;
412  t1code = 2212;
413  t2code = 2112;
414  scode = 2212;
415  s2code = 2112;
416  }
417  else
418  {
419  LOG("HNIntranuke2014", pWARN)
420  << "AbsorbHN() cannot handle probe: " << pdgc;
421  return;
422  }
423 
424  // assign proper masses
425  M1 = pLib->Find(pcode) ->Mass();
426  M2_1 = pLib->Find(t1code)->Mass();
427  M2_2 = pLib->Find(t2code)->Mass();
428  M3 = pLib->Find(scode) ->Mass();
429  M4 = pLib->Find(s2code)->Mass();
430 
431  // handle fermi momentum
432  if(fDoFermi)
433  {
434  target.SetHitNucPdg(t1code);
435  fNuclmodel->GenerateNucleon(target);
436  tP2_1L=fFermiFac * fNuclmodel->Momentum3();
437  E2_1L = TMath::Sqrt(tP2_1L.Mag2() + M2_1*M2_1);
438 
439  target.SetHitNucPdg(t2code);
440  fNuclmodel->GenerateNucleon(target);
441  tP2_2L=fFermiFac * fNuclmodel->Momentum3();
442  E2_2L = TMath::Sqrt(tP2_2L.Mag2() + M2_2*M2_2);
443  }
444  else
445  {
446  tP2_1L.SetXYZ(0.0, 0.0, 0.0);
447  E2_1L = M2_1;
448 
449  tP2_2L.SetXYZ(0.0, 0.0, 0.0);
450  E2_2L = M2_2;
451  }
452 
453  E2L = E2_1L + E2_2L;
454 
455  // adjust p to reflect scattering
456  // get random scattering angle
457  C3CM = fHadroData2014->IntBounce(p,t1code,scode,fate);
458  if (C3CM<-1.)
459  {
460  p->SetStatus(kIStStableFinalState);
461  ev->AddParticle(*p);
462  return;
463  }
464  S3CM = TMath::Sqrt(1.0 - C3CM*C3CM);
465 
466  // Get lab energy and momenta
467  E1L = p->E();
468  if(E1L<0.001) E1L=0.001;
469  P1L = TMath::Sqrt(E1L*E1L - M1*M1);
470  tP1L = p->P4()->Vect();
471  tP2L = tP2_1L + tP2_2L;
472  P2L = tP2L.Mag();
473  tPtot = tP1L + tP2L;
474 
475  // get unit vectors and angles needed for later
476  bDir = tPtot.Unit();
477  Theta1 = tP1L.Angle(bDir);
478  Theta2 = tP2L.Angle(bDir);
479 
480  // get parallel and transverse components
481  P1zL = P1L*TMath::Cos(Theta1);
482  P2zL = P2L*TMath::Cos(Theta2);
483  tVect.SetXYZ(1,0,0);
484  if(TMath::Abs((tVect - bDir).Mag())<.01) tVect.SetXYZ(0,1,0);
485  theta5 = tVect.Angle(bDir);
486  tTrans = (tVect - TMath::Cos(theta5)*bDir).Unit();
487 
488  // calculate beta and gamma
489  tbeta = tPtot * (1.0 / (E1L + E2L));
490  beta = tbeta.Mag();
491  gm = 1.0 / TMath::Sqrt(1.0 - beta*beta);
492 
493  // boost to CM frame to get scattered particle momenta
494  E1CM = gm*E1L - gm*beta*P1zL;
495  P1zCM = gm*P1zL*bDir - gm*tbeta*E1L;
496  E2CM = gm*E2L - gm*beta*P2zL;
497  P2zCM = gm*P2zL*bDir - gm*tbeta*E2L;
498  Et = E1CM + E2CM;
499  E3CM = (Et*Et + (M3*M3) - (M4*M4)) / (2.0*Et);
500  E4CM = Et - E3CM;
501  P3CM = TMath::Sqrt(E3CM*E3CM - M3*M3);
502 
503  // boost back to lab
504  P3zL = gm*beta*E3CM + gm*P3CM*C3CM;
505  P3tL = P3CM*S3CM;
506  P4zL = gm*beta*E4CM + gm*P3CM*(-C3CM);
507  P4tL = P3CM*(-S3CM);
508 
509  P3L = TMath::Sqrt(P3zL*P3zL + P3tL*P3tL);
510  P4L = TMath::Sqrt(P4zL*P4zL + P4tL*P4tL);
511 
512  // check for too small values
513  // may introduce error, so warn if it occurs
514  if(!(TMath::Finite(P3L))||P3L<.001)
515  {
516  LOG("HNIntranuke2014",pINFO)
517  << "Particle 3 " << M3 << " momentum small or non-finite: " << P3L
518  << "\n" << "--> Assigning .001 as new momentum";
519  P3tL = 0;
520  P3zL = .001;
521  P3L = .001;
522  E3L = TMath::Sqrt(P3L*P3L + M3*M3);
523  }
524 
525  if(!(TMath::Finite(P4L))||P4L<.001)
526  {
527  LOG("HNIntranuke2014",pINFO)
528  << "Particle 4 " << M4 << " momentum small or non-finite: " << P4L
529  << "\n" << "--> Assigning .001 as new momentum";
530  P4tL = 0;
531  P4zL = .001;
532  P4L = .001;
533  E4L = TMath::Sqrt(P4L*P4L + M4*M4);
534  }
535 
536  // pauli blocking
537  if(P3L < fFermiMomentum || P4L < fFermiMomentum)
538  {
539 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
540  LOG("HNIntranuke2014",pINFO) << "AbsorbHN failed: Pauli blocking";
541 #endif
542  p->SetStatus(kIStHadronInTheNucleus);
543 
544  //utils::intranuke2014::StepParticle(p,fFreeStep,fTrackingRadius); //not needed
545 
546  ev->AddParticle(*p);
547  return;
548  }
549 
550  // handle remnant nucleus updates
551  fRemnZ--;
552  fRemnA -=2;
553  fRemnP4 -= TLorentzVector(tP2_1L,E2_1L);
554  fRemnP4 -= TLorentzVector(tP2_2L,E2_2L);
555 
556  // get random phi angle, distributed uniformally in 360 deg
557  PHI3 = 2 * kPi * rnd->RndFsi().Rndm();
558 
559  tP3L = P3zL*bDir + P3tL*tTrans;
560  tP4L = P4zL*bDir + P4tL*tTrans;
561 
562  tP3L.Rotate(PHI3,bDir); // randomize transverse components
563  tP4L.Rotate(PHI3,bDir);
564 
565  E3L = TMath::Sqrt(P3L*P3L + M3*M3);
566  E4L = TMath::Sqrt(P4L*P4L + M4*M4);
567 
568  // create t particle w/ appropriate momenta, code, and status
569  // set target's mom to be the mom of the hadron that was cloned
570  GHepParticle * t = new GHepParticle(*p);
571  t->SetFirstMother(p->FirstMother());
572  t->SetLastMother(p->LastMother());
573 
574  TLorentzVector t4P4L(tP4L,E4L);
575  t->SetPdgCode(s2code);
576  t->SetMomentum(t4P4L);
577  t->SetStatus(kIStHadronInTheNucleus);
578 
579  // adjust p to reflect scattering
580  TLorentzVector t4P3L(tP3L,E3L);
581  p->SetPdgCode(scode);
582  p->SetMomentum(t4P3L);
583  p->SetStatus(kIStHadronInTheNucleus);
584 
585 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
586  LOG("HNIntranuke2014", pDEBUG)
587  << "|p3| = " << (P3L) << ", E3 = " << (E3L);
588  LOG("HNIntranuke2014", pDEBUG)
589  << "|p4| = " << (P4L) << ", E4 = " << (E4L);
590 #endif
591 
592  ev->AddParticle(*p);
593  ev->AddParticle(*t);
594 
595  delete t; // delete particle clone
596 }
genie::GHepParticle::SetFirstMother
void SetFirstMother(int m)
Definition: GHepParticle.h:131
genie::RandomGen::RndFsi
TRandom3 & RndFsi(void) const
rnd number generator used by intranuclear cascade monte carlos
Definition: RandomGen.h:60
genie::GHepParticle::E
double E(void) const
Get energy.
Definition: GHepParticle.h:90
genie::RandomGen::Instance
static RandomGen * Instance()
Access instance.
Definition: RandomGen.cxx:79
generate_CCQE_events.t
t
Definition: generate_CCQE_events.py:68
genie::Intranuke2014::fRemnZ
int fRemnZ
remnant nucleus Z
Definition: Intranuke2014.h:97
genie::Intranuke2014::fHadroData2014
INukeHadroData2014 * fHadroData2014
a collection of h+N,h+A data & calculations
Definition: Intranuke2014.h:93
genie::GHepParticle::SetMomentum
void SetMomentum(const TLorentzVector &p4)
Definition: GHepParticle.cxx:254
genie::Intranuke2014::fFermiMomentum
double fFermiMomentum
whether or not particle collision is pauli blocked
Definition: Intranuke2014.h:113
genie::RandomGen
A singleton holding random number generator classes. All random number generation in GENIE should tak...
Definition: RandomGen.h:30
generate_CCQE_events.target
target
Definition: generate_CCQE_events.py:44
genie::GHepParticle::LastMother
int LastMother(void) const
Definition: GHepParticle.h:68
genie::GHepParticle::Pdg
int Pdg(void) const
Definition: GHepParticle.h:64
genie::GHepParticle::FirstMother
int FirstMother(void) const
Definition: GHepParticle.h:67
genie::GHepParticle::Name
string Name(void) const
Name that corresponds to the PDG code.
Definition: GHepParticle.cxx:153
genie::Intranuke2014::fRemnA
int fRemnA
remnant nucleus A
Definition: Intranuke2014.h:96
genie::Intranuke2014::fNuclmodel
const NuclearModelI * fNuclmodel
nuclear model used to generate fermi momentum
Definition: Intranuke2014.h:95
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::Target
A Neutrino Interaction Target. Is a transparent encapsulation of quite different physical systems suc...
Definition: Target.h:41
genie::GHepRecord::TargetNucleus
virtual GHepParticle * TargetNucleus(void) const
Definition: GHepRecord.cxx:323
pINFO
#define pINFO
Definition: Messenger.h:53
pWARN
#define pWARN
Definition: Messenger.h:51
genie::Intranuke2014::fFermiFac
double fFermiFac
testing parameter to modify fermi momentum
Definition: Intranuke2014.h:111
genie::PDGLibrary::Instance
static PDGLibrary * Instance(void)
Definition: PDGLibrary.cxx:45
genie::GHepParticle::SetLastMother
void SetLastMother(int m)
Definition: GHepParticle.h:132
genie::PDGLibrary
Singleton class to load & serve a TDatabasePDG.
Definition: PDGLibrary.h:27
genie::GHepParticle::SetStatus
void SetStatus(GHepStatus_t s)
Definition: GHepParticle.h:125
genie::NuclearModelI::Momentum3
virtual TVector3 Momentum3(void) const
Definition: NuclearModelI.h:59
genie::INukeHadroFates::AsString
static string AsString(INukeFateHN_t fate)
Definition: INukeHadroFates.h:68
genie::PDGLibrary::Find
TParticlePDG * Find(int pdgc)
Definition: PDGLibrary.cxx:63
genie::NuclearModelI::GenerateNucleon
virtual bool GenerateNucleon(const Target &) const =0
genie::GHepRecord::AddParticle
virtual void AddParticle(const GHepParticle &p)
Definition: GHepRecord.cxx:526
genie::GHepParticle::P4
TLorentzVector * P4(void) const
Definition: GHepParticle.h:79
genie::constants::kPi
static const double kPi
Definition: Constants.h:38
genie::GHepParticle
STDHEP-like event record entry that can fit a particle or a nucleus.
Definition: GHepParticle.h:40
genie::Intranuke2014::fRemnP4
TLorentzVector fRemnP4
P4 of remnant system.
Definition: Intranuke2014.h:98
genie::INukeHadroData2014::IntBounce
double IntBounce(const GHepParticle *p, int target, int s1, INukeFateHN_t fate)
Definition: INukeHadroData2014.cxx:1517
pDEBUG
#define pDEBUG
Definition: Messenger.h:54
genie::GHepParticle::SetPdgCode
void SetPdgCode(int c)
Definition: GHepParticle.cxx:248
genie::Intranuke2014::fDoFermi
bool fDoFermi
whether or not to do fermi mom.
Definition: Intranuke2014.h:114
void HNIntranuke2014::ElasHN ( GHepRecord *  ev,
GHepParticle *  p,
INukeFateHN_t  fate 
) const
private

Definition at line 598 of file HNIntranuke2014.cxx.

600 {
601  // scatters particle within nucleus
602 
603 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
604  LOG("HNIntranuke2014", pDEBUG)
605  << "ElasHN() is invoked for a : " << p->Name()
606  << " whose fate is : " << INukeHadroFates::AsString(fate);
607 #endif
608 
609  if(fate!=kIHNFtCEx && fate!=kIHNFtElas)
610  {
611  LOG("HNIntranuke2014", pWARN)
612  << "ElasHN() cannot handle fate: " << INukeHadroFates::AsString(fate);
613  return;
614  }
615 
616  // Random number generator
617  RandomGen * rnd = RandomGen::Instance();
618 
619  // vars for incoming particle, target, and scattered pdg codes
620  int pcode = p->Pdg();
621  int tcode, scode, s2code;
622  double ppcnt = (double) fRemnZ / (double) fRemnA; // % of protons
623 
624  // Select a target randomly, weighted to #
625  // -- Unless, of course, the fate is CEx,
626  // -- in which case the target may be deterministic
627  // Also assign scattered particle code
628  if(fate==kIHNFtCEx)
629  {
630  if(pcode==kPdgPiP) {tcode = kPdgNeutron; scode = kPdgPi0; s2code = kPdgProton;}
631  else if(pcode==kPdgPiM) {tcode = kPdgProton; scode = kPdgPi0; s2code = kPdgNeutron;}
632  else
633  {
634  // for pi0
635  tcode = (rnd->RndFsi().Rndm()<=ppcnt)?(kPdgProton) :(kPdgNeutron);
636  scode = (tcode == kPdgProton) ?(kPdgPiP) :(kPdgPiM);
637  s2code = (tcode == kPdgProton) ?(kPdgNeutron):(kPdgProton);
638  }
639  }
640  else
641  {
642  tcode = (rnd->RndFsi().Rndm()<=ppcnt)?(kPdgProton):(kPdgNeutron);
643  scode = pcode;
644  s2code = tcode;
645  }
646 
647  // get random scattering angle
648  double C3CM = fHadroData2014->IntBounce(p,tcode,scode,fate);
649  if (C3CM<-1.)
650  {
651  p->SetStatus(kIStStableFinalState);
652  ev->AddParticle(*p);
653  return;
654  }
655 
656  // create scattered particle
657  GHepParticle * t = new GHepParticle(*p);
658  t->SetPdgCode(tcode);
659  double Mt = t->Mass();
660  //t->SetMomentum(TLorentzVector(0,0,0,Mt));
661 
662  // handle fermi momentum
663  if(fDoFermi)
664  {
665  // Handle fermi target
666  Target target(ev->TargetNucleus()->Pdg());
667 
668  target.SetHitNucPdg(tcode);
669  fNuclmodel->GenerateNucleon(target);
670  TVector3 tP3L = fFermiFac * fNuclmodel->Momentum3();
671  double tE = TMath::Sqrt(tP3L.Mag2() + Mt*Mt);
672  t->SetMomentum(TLorentzVector(tP3L,tE));
673  }
674  else
675  {
676  t->SetMomentum(TLorentzVector(0,0,0,Mt));
677  }
678 
679  bool pass = utils::intranuke2014::TwoBodyCollision(ev,pcode,tcode,scode,s2code,C3CM,
680  p,t,fRemnA,fRemnZ,fRemnP4,kIMdHN);
681 
682 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
683  LOG("HNIntranuke2014",pDEBUG)
684  << "|p3| = " << P3L << ", E3 = " << E3L;
685  LOG("HNIntranuke2014",pDEBUG)
686  << "|p4| = " << P4L << ", E4 = " << E4L;
687 #endif
688 
689  if (pass==true)
690  {
691 
692  //utils::intranuke2014::StepParticle(p,fFreeStep,fTrackingRadius); //not needed
693  //utils::intranuke2014::StepParticle(t,fFreeStep,fTrackingRadius);
694 
695  ev->AddParticle(*p);
696  ev->AddParticle(*t);
697  } else
698  {
699 
700  delete t; //fixes a troublesome memory leak
701 
702  LOG("HNIntranuke2014", pINFO) << "Elastic in hN failed calling TwoBodyCollision";
703  exceptions::INukeException exception;
704  exception.SetReason("hN scattering kinematics through TwoBodyCollision failed");
705  throw exception;
706  }
707 
708  delete t;
709 
710 }
genie::exceptions::INukeException
An exception thrown by SimulateHadronState for kinematics problems. TwoBodyCollision/Kinematics used ...
Definition: INukeException.h:38
genie::RandomGen::RndFsi
TRandom3 & RndFsi(void) const
rnd number generator used by intranuclear cascade monte carlos
Definition: RandomGen.h:60
genie::utils::intranuke2014::TwoBodyCollision
bool TwoBodyCollision(GHepRecord *ev, int pcode, int tcode, int scode, int s2code, double C3CM, GHepParticle *p, GHepParticle *t, int &RemnA, int &RemnZ, TLorentzVector &RemnP4, EINukeMode mode=kIMdHA)
Intranuke utility functions.
Definition: INukeUtils2014.cxx:628
genie::RandomGen::Instance
static RandomGen * Instance()
Access instance.
Definition: RandomGen.cxx:79
generate_CCQE_events.t
t
Definition: generate_CCQE_events.py:68
genie::Intranuke2014::fRemnZ
int fRemnZ
remnant nucleus Z
Definition: Intranuke2014.h:97
genie::Intranuke2014::fHadroData2014
INukeHadroData2014 * fHadroData2014
a collection of h+N,h+A data & calculations
Definition: Intranuke2014.h:93
genie::GHepParticle::SetMomentum
void SetMomentum(const TLorentzVector &p4)
Definition: GHepParticle.cxx:254
genie::GHepParticle::Mass
double Mass(void) const
Mass that corresponds to the PDG code.
Definition: GHepParticle.cxx:161
genie::RandomGen
A singleton holding random number generator classes. All random number generation in GENIE should tak...
Definition: RandomGen.h:30
generate_CCQE_events.target
target
Definition: generate_CCQE_events.py:44
genie::GHepParticle::Pdg
int Pdg(void) const
Definition: GHepParticle.h:64
genie::GHepParticle::Name
string Name(void) const
Name that corresponds to the PDG code.
Definition: GHepParticle.cxx:153
genie::Intranuke2014::fRemnA
int fRemnA
remnant nucleus A
Definition: Intranuke2014.h:96
pass
const HLTPathStatus pass
Definition: EventSelExc_t.cpp:90
genie::Intranuke2014::fNuclmodel
const NuclearModelI * fNuclmodel
nuclear model used to generate fermi momentum
Definition: Intranuke2014.h:95
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::Target
A Neutrino Interaction Target. Is a transparent encapsulation of quite different physical systems suc...
Definition: Target.h:41
genie::kPdgPiP
const int kPdgPiP
Definition: PDGCodes.h:132
genie::kPdgPi0
const int kPdgPi0
Definition: PDGCodes.h:134
genie::GHepRecord::TargetNucleus
virtual GHepParticle * TargetNucleus(void) const
Definition: GHepRecord.cxx:323
pINFO
#define pINFO
Definition: Messenger.h:53
pWARN
#define pWARN
Definition: Messenger.h:51
genie::Intranuke2014::fFermiFac
double fFermiFac
testing parameter to modify fermi momentum
Definition: Intranuke2014.h:111
genie::GHepParticle::SetStatus
void SetStatus(GHepStatus_t s)
Definition: GHepParticle.h:125
genie::NuclearModelI::Momentum3
virtual TVector3 Momentum3(void) const
Definition: NuclearModelI.h:59
genie::INukeHadroFates::AsString
static string AsString(INukeFateHN_t fate)
Definition: INukeHadroFates.h:68
genie::kPdgPiM
const int kPdgPiM
Definition: PDGCodes.h:133
genie::NuclearModelI::GenerateNucleon
virtual bool GenerateNucleon(const Target &) const =0
genie::GHepRecord::AddParticle
virtual void AddParticle(const GHepParticle &p)
Definition: GHepRecord.cxx:526
genie::kPdgProton
const int kPdgProton
Definition: PDGCodes.h:62
genie::kPdgNeutron
const int kPdgNeutron
Definition: PDGCodes.h:64
genie::GHepParticle
STDHEP-like event record entry that can fit a particle or a nucleus.
Definition: GHepParticle.h:40
genie::Intranuke2014::fRemnP4
TLorentzVector fRemnP4
P4 of remnant system.
Definition: Intranuke2014.h:98
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
genie::INukeHadroData2014::IntBounce
double IntBounce(const GHepParticle *p, int target, int s1, INukeFateHN_t fate)
Definition: INukeHadroData2014.cxx:1517
pDEBUG
#define pDEBUG
Definition: Messenger.h:54
genie::GHepParticle::SetPdgCode
void SetPdgCode(int c)
Definition: GHepParticle.cxx:248
genie::Intranuke2014::fDoFermi
bool fDoFermi
whether or not to do fermi mom.
Definition: Intranuke2014.h:114
double HNIntranuke2014::FateWeight ( int  pdgc,
INukeFateHN_t  fate 
) const
private

Definition at line 311 of file HNIntranuke2014.cxx.

312 {
313  // turn fates off if the remnant nucleus does not have the number of p,n
314  // required
315 
316  int np = fRemnZ;
317  int nn = fRemnA - fRemnZ;
318 
319  if (np < 1 && nn < 1)
320  {
321  LOG("HNIntranuke2014", pERROR) << "** Nothing left in nucleus!!! **";
322  return 0;
323  }
324  else
325  {
326  if (fate == kIHNFtCEx && pdgc==kPdgPiP ) { return (nn>=1) ? 1. : 0.; }
327  if (fate == kIHNFtCEx && pdgc==kPdgPiM ) { return (np>=1) ? 1. : 0.; }
328  if (fate == kIHNFtAbs) { return ((nn>=1) && (np>=1)) ? 1. : 0.; }
329  if (fate == kIHNFtCmp ) { return ((pdgc==kPdgProton||pdgc==kPdgNeutron)&&fDoCompoundNucleus&&fRemnA>5) ? 1. : 0.; }
330 
331  }
332  return 1.;
333 }
pERROR
#define pERROR
Definition: Messenger.h:50
genie::Intranuke2014::fRemnZ
int fRemnZ
remnant nucleus Z
Definition: Intranuke2014.h:97
genie::Intranuke2014::fRemnA
int fRemnA
remnant nucleus A
Definition: Intranuke2014.h:96
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::kPdgPiP
const int kPdgPiP
Definition: PDGCodes.h:132
genie::Intranuke2014::fDoCompoundNucleus
bool fDoCompoundNucleus
whether or not to do compound nucleus considerations
Definition: Intranuke2014.h:116
genie::kPdgPiM
const int kPdgPiM
Definition: PDGCodes.h:133
genie::kPdgProton
const int kPdgProton
Definition: PDGCodes.h:62
genie::kPdgNeutron
const int kPdgNeutron
Definition: PDGCodes.h:64
void HNIntranuke2014::GammaInelasticHN ( GHepRecord *  ev,
GHepParticle *  p,
INukeFateHN_t  fate 
) const
private

Definition at line 753 of file HNIntranuke2014.cxx.

754 {
755  // This function handles pion photoproduction reactions
756 
757 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
758  LOG("HNIntranuke2014", pDEBUG)
759  << "GammaInelasticHN() is invoked for a : " << p->Name()
760  << " whose fate is : " << INukeHadroFates::AsString(fate);
761 #endif
762 
763  if(fate!=kIHNFtInelas && p->Pdg()!=kPdgGamma)
764  {
765  LOG("HNIntranuke2014", pWARN)
766  << "GammaInelasticHN() cannot handle fate: " << INukeHadroFates::AsString(fate);
767  return;
768  }
769 
770  // random number generator
771  RandomGen * rnd = RandomGen::Instance();
772 
773  // vars for incoming particle, target, and scattered reaction products
774  double ppcnt = (double) fRemnZ / (double) fRemnA; // % of protons
775  int pcode = p->Pdg();
776  int tcode = (rnd->RndFsi().Rndm()<=ppcnt)?(kPdgProton):(kPdgNeutron);
777  int scode, s2code;
778  double ke = p->KinE() / units::MeV;
779 
780  LOG("HNIntranuke2014", pNOTICE)
781  << "Particle code: " << pcode << ", target: " << tcode;
782 
783 
784  if (rnd->RndFsi().Rndm() * (fHadroData2014 -> XSecGamp_fs() -> Evaluate(ke) +
785  fHadroData2014 -> XSecGamn_fs() -> Evaluate(ke) )
786  <= fHadroData2014 -> XSecGamp_fs() -> Evaluate(ke) ) { scode = kPdgProton; }
787  else { scode = kPdgNeutron; }
788 
789  //scode=fHadroData2014->AngleAndProduct(p,tcode,C3CM,fate);
790  //double C3CM = 0.0; // cos of scattering angle
791  double C3CM = fHadroData2014->IntBounce(p,tcode,scode,fate);
792 
793  if ((tcode == kPdgProton ) && (scode==kPdgProton )) {s2code=kPdgPi0;}
794  else if ((tcode == kPdgProton ) && (scode==kPdgNeutron)) {s2code=kPdgPiP;}
795  else if ((tcode == kPdgNeutron) && (scode==kPdgProton )) {s2code=kPdgPiM;}
796  else if ((tcode == kPdgNeutron) && (scode==kPdgNeutron)) {s2code=kPdgPi0;}
797  else {
798  LOG("HNIntranuke2014", pERROR)
799  << "Error: could not determine particle final states";
800  ev->AddParticle(*p);
801  return;
802  }
803 
804  LOG("HNIntranuke2014", pNOTICE)
805  << "GammaInelastic fate: " << INukeHadroFates::AsString(fate);
806  LOG("HNIntranuke2014", pNOTICE)
807  << " final state: " << scode << " and " << s2code;
808  LOG("HNIntranuke2014", pNOTICE)
809  << " scattering angle: " << C3CM;
810 
811  GHepParticle * t = new GHepParticle(*p);
812  t->SetPdgCode(tcode);
813  double Mt = t->Mass();
814 
815  // handle fermi momentum
816  if(fDoFermi)
817  {
818  // Handle fermi target
819  Target target(ev->TargetNucleus()->Pdg());
820 
821  target.SetHitNucPdg(tcode);
822  fNuclmodel->GenerateNucleon(target);
823  TVector3 tP3L = fFermiFac * fNuclmodel->Momentum3();
824  double tE = TMath::Sqrt(tP3L.Mag2() + Mt*Mt);
825  t->SetMomentum(TLorentzVector(tP3L,tE));
826  }
827  else
828  {
829  t->SetMomentum(TLorentzVector(0,0,0,Mt));
830  }
831 
832  bool pass = utils::intranuke2014::TwoBodyCollision(ev,pcode,tcode,scode,s2code,C3CM,
833  p,t,fRemnA,fRemnZ,fRemnP4,kIMdHN);
834 
835 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
836  LOG("HNIntranuke2014",pDEBUG)
837  << "|p3| = " << P3L << ", E3 = " << E3L;
838  LOG("HNIntranuke2014",pDEBUG)
839  << "|p4| = " << P4L << ", E4 = " << E4L;
840 #endif
841 
842  if (pass==true)
843  {
844  //p->SetStatus(kIStStableFinalState);
845  //t->SetStatus(kIStStableFinalState);
846  ev->AddParticle(*p);
847  ev->AddParticle(*t);
848  } else
849  {
850  ev->AddParticle(*p);
851  }
852 
853  delete t;
854 
855 }
genie::RandomGen::RndFsi
TRandom3 & RndFsi(void) const
rnd number generator used by intranuclear cascade monte carlos
Definition: RandomGen.h:60
pERROR
#define pERROR
Definition: Messenger.h:50
genie::utils::intranuke2014::TwoBodyCollision
bool TwoBodyCollision(GHepRecord *ev, int pcode, int tcode, int scode, int s2code, double C3CM, GHepParticle *p, GHepParticle *t, int &RemnA, int &RemnZ, TLorentzVector &RemnP4, EINukeMode mode=kIMdHA)
Intranuke utility functions.
Definition: INukeUtils2014.cxx:628
genie::RandomGen::Instance
static RandomGen * Instance()
Access instance.
Definition: RandomGen.cxx:79
generate_CCQE_events.t
t
Definition: generate_CCQE_events.py:68
genie::units::MeV
static const double MeV
Definition: Units.h:130
genie::Intranuke2014::fRemnZ
int fRemnZ
remnant nucleus Z
Definition: Intranuke2014.h:97
genie::Intranuke2014::fHadroData2014
INukeHadroData2014 * fHadroData2014
a collection of h+N,h+A data & calculations
Definition: Intranuke2014.h:93
genie::GHepParticle::SetMomentum
void SetMomentum(const TLorentzVector &p4)
Definition: GHepParticle.cxx:254
genie::GHepParticle::Mass
double Mass(void) const
Mass that corresponds to the PDG code.
Definition: GHepParticle.cxx:161
genie::RandomGen
A singleton holding random number generator classes. All random number generation in GENIE should tak...
Definition: RandomGen.h:30
generate_CCQE_events.target
target
Definition: generate_CCQE_events.py:44
genie::GHepParticle::Pdg
int Pdg(void) const
Definition: GHepParticle.h:64
genie::GHepParticle::Name
string Name(void) const
Name that corresponds to the PDG code.
Definition: GHepParticle.cxx:153
genie::Intranuke2014::fRemnA
int fRemnA
remnant nucleus A
Definition: Intranuke2014.h:96
pass
const HLTPathStatus pass
Definition: EventSelExc_t.cpp:90
genie::Intranuke2014::fNuclmodel
const NuclearModelI * fNuclmodel
nuclear model used to generate fermi momentum
Definition: Intranuke2014.h:95
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::kPdgGamma
const int kPdgGamma
Definition: PDGCodes.h:163
genie::Target
A Neutrino Interaction Target. Is a transparent encapsulation of quite different physical systems suc...
Definition: Target.h:41
genie::kPdgPiP
const int kPdgPiP
Definition: PDGCodes.h:132
genie::kPdgPi0
const int kPdgPi0
Definition: PDGCodes.h:134
genie::GHepRecord::TargetNucleus
virtual GHepParticle * TargetNucleus(void) const
Definition: GHepRecord.cxx:323
pWARN
#define pWARN
Definition: Messenger.h:51
genie::Intranuke2014::fFermiFac
double fFermiFac
testing parameter to modify fermi momentum
Definition: Intranuke2014.h:111
genie::GHepParticle::KinE
double KinE(bool mass_from_pdg=false) const
Get kinetic energy.
Definition: GHepParticle.cxx:177
genie::NuclearModelI::Momentum3
virtual TVector3 Momentum3(void) const
Definition: NuclearModelI.h:59
genie::INukeHadroFates::AsString
static string AsString(INukeFateHN_t fate)
Definition: INukeHadroFates.h:68
genie::kPdgPiM
const int kPdgPiM
Definition: PDGCodes.h:133
genie::NuclearModelI::GenerateNucleon
virtual bool GenerateNucleon(const Target &) const =0
genie::GHepRecord::AddParticle
virtual void AddParticle(const GHepParticle &p)
Definition: GHepRecord.cxx:526
genie::kPdgProton
const int kPdgProton
Definition: PDGCodes.h:62
pNOTICE
#define pNOTICE
Definition: Messenger.h:52
genie::kPdgNeutron
const int kPdgNeutron
Definition: PDGCodes.h:64
genie::GHepParticle
STDHEP-like event record entry that can fit a particle or a nucleus.
Definition: GHepParticle.h:40
genie::Intranuke2014::fRemnP4
TLorentzVector fRemnP4
P4 of remnant system.
Definition: Intranuke2014.h:98
genie::INukeHadroData2014::IntBounce
double IntBounce(const GHepParticle *p, int target, int s1, INukeFateHN_t fate)
Definition: INukeHadroData2014.cxx:1517
pDEBUG
#define pDEBUG
Definition: Messenger.h:54
genie::GHepParticle::SetPdgCode
void SetPdgCode(int c)
Definition: GHepParticle.cxx:248
genie::Intranuke2014::fDoFermi
bool fDoFermi
whether or not to do fermi mom.
Definition: Intranuke2014.h:114
INukeFateHN_t HNIntranuke2014::HadronFateHN ( const GHepParticle *  p) const
private

Definition at line 200 of file HNIntranuke2014.cxx.

201 {
202 // Select a hadron fate in HN mode
203 //
204  RandomGen * rnd = RandomGen::Instance();
205 
206  // get pdgc code & kinetic energy in MeV
207  int pdgc = p->Pdg();
208  double ke = p->KinE() / units::MeV;
209 
210  LOG("HNIntranuke2014", pINFO)
211  << "Selecting hN fate for " << p->Name() << " with KE = " << ke << " MeV";
212 
213  // try to generate a hadron fate
214  unsigned int iter = 0;
215  while(iter++ < kRjMaxIterations) {
216 
217  // handle pions
218  //
219  if (pdgc==kPdgPiP || pdgc==kPdgPiM || pdgc==kPdgPi0) {
220 
221  double frac_cex = this->FateWeight(pdgc, kIHNFtCEx)
222  * fHadroData2014->Frac(pdgc, kIHNFtCEx, ke, fRemnA, fRemnZ);
223  double frac_elas = this->FateWeight(pdgc, kIHNFtElas)
224  * fHadroData2014->Frac(pdgc, kIHNFtElas, ke, fRemnA, fRemnZ);
225  double frac_inel = this->FateWeight(pdgc, kIHNFtInelas)
226  * fHadroData2014->Frac(pdgc, kIHNFtInelas, ke, fRemnA, fRemnZ);
227  double frac_abs = this->FateWeight(pdgc, kIHNFtAbs)
228  * fHadroData2014->Frac(pdgc, kIHNFtAbs, ke, fRemnA, fRemnZ);
229  frac_cex *= fNucCEXFac; // scaling factors
230  frac_abs *= fNucAbsFac;
231  frac_elas *= fNucQEFac;
232  if(pdgc==kPdgPi0) frac_abs*= 0.665; //isospin factor
233 
234  LOG("HNIntranuke2014", pINFO)
235  << "\n frac{" << INukeHadroFates::AsString(kIHNFtCEx) << "} = " << frac_cex
236  << "\n frac{" << INukeHadroFates::AsString(kIHNFtElas) << "} = " << frac_elas
237  << "\n frac{" << INukeHadroFates::AsString(kIHNFtInelas) << "} = " << frac_inel
238  << "\n frac{" << INukeHadroFates::AsString(kIHNFtAbs) << "} = " << frac_abs;
239 
240  // compute total fraction (can be <1 if fates have been switched off)
241  double tf = frac_cex +
242  frac_elas +
243  frac_inel +
244  frac_abs;
245 
246  double r = tf * rnd->RndFsi().Rndm();
247 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
248  LOG("HNIntranuke2014", pDEBUG) << "r = " << r << " (max = " << tf << ")";
249 #endif
250 
251  double cf=0; // current fraction
252 
253  if(r < (cf += frac_cex )) return kIHNFtCEx; //cex
254  if(r < (cf += frac_elas )) return kIHNFtElas; //elas
255  if(r < (cf += frac_inel )) return kIHNFtInelas; //inelas
256  if(r < (cf += frac_abs )) return kIHNFtAbs; //abs
257 
258  LOG("HNIntranuke2014", pWARN)
259  << "No selection after going through all fates! "
260  << "Total fraction = " << tf << " (r = " << r << ")";
261  ////////////////////////////
262  return kIHNFtUndefined;
263  }
264 
265  // handle nucleons
266  else if (pdgc==kPdgProton || pdgc==kPdgNeutron) {
267 
268  double frac_elas = this->FateWeight(pdgc, kIHNFtElas)
269  * fHadroData2014->Frac(pdgc, kIHNFtElas, ke, fRemnA, fRemnZ);
270  double frac_inel = this->FateWeight(pdgc, kIHNFtInelas)
271  * fHadroData2014->Frac(pdgc, kIHNFtInelas, ke, fRemnA, fRemnZ);
272  double frac_cmp = this->FateWeight(pdgc, kIHNFtCmp)
273  * fHadroData2014->Frac(pdgc, kIHNFtCmp, ke, fRemnA , fRemnZ);
274 
275  LOG("HNIntranuke2014", pINFO)
276  << "\n frac{" << INukeHadroFates::AsString(kIHNFtElas) << "} = " << frac_elas
277  << "\n frac{" << INukeHadroFates::AsString(kIHNFtInelas) << "} = " << frac_inel;
278 
279  // compute total fraction (can be <1 if fates have been switched off)
280  double tf = frac_elas + frac_inel + frac_cmp;
281 
282  double r = tf * rnd->RndFsi().Rndm();
283 
284 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
285  LOG("HNIntranuke2014", pDEBUG) << "r = " << r << " (max = " << tf << ")";
286 #endif
287 
288  double cf=0; // current fraction
289  if(r < (cf += frac_elas )) return kIHNFtElas; // elas
290  if(r < (cf += frac_inel )) return kIHNFtInelas; // inelas
291  if(r < (cf += frac_cmp )) return kIHNFtCmp; // cmp
292 
293  LOG("HNIntranuke2014", pWARN)
294  << "No selection after going through all fates! "
295  << "Total fraction = " << tf << " (r = " << r << ")";
296  //////////////////////////
297  return kIHNFtUndefined;
298  }
299 
300  // handle gamma -- does not currently consider the elastic case
301  else if (pdgc==kPdgGamma) return kIHNFtInelas;
302 
303  // handle kaon -- currently only elastic case
304  else if (pdgc==kPdgKP) return kIHNFtElas;
305 
306  }//iterations
307 
308  return kIHNFtUndefined;
309 }
genie::RandomGen::RndFsi
TRandom3 & RndFsi(void) const
rnd number generator used by intranuclear cascade monte carlos
Definition: RandomGen.h:60
genie::RandomGen::Instance
static RandomGen * Instance()
Access instance.
Definition: RandomGen.cxx:79
genie::units::MeV
static const double MeV
Definition: Units.h:130
genie::Intranuke2014::fRemnZ
int fRemnZ
remnant nucleus Z
Definition: Intranuke2014.h:97
genie::Intranuke2014::fHadroData2014
INukeHadroData2014 * fHadroData2014
a collection of h+N,h+A data & calculations
Definition: Intranuke2014.h:93
tf
Definition: tf_graph.h:23
genie::RandomGen
A singleton holding random number generator classes. All random number generation in GENIE should tak...
Definition: RandomGen.h:30
genie::GHepParticle::Pdg
int Pdg(void) const
Definition: GHepParticle.h:64
genie::GHepParticle::Name
string Name(void) const
Name that corresponds to the PDG code.
Definition: GHepParticle.cxx:153
genie::Intranuke2014::fRemnA
int fRemnA
remnant nucleus A
Definition: Intranuke2014.h:96
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::HNIntranuke2014::FateWeight
double FateWeight(int pdgc, INukeFateHN_t fate) const
Definition: HNIntranuke2014.cxx:311
genie::kPdgGamma
const int kPdgGamma
Definition: PDGCodes.h:163
genie::INukeHadroData2014::Frac
double Frac(int hpdgc, INukeFateHN_t fate, double ke, int targA=0, int targZ=0) const
Definition: INukeHadroData2014.cxx:1486
genie::kPdgKP
const int kPdgKP
Definition: PDGCodes.h:146
genie::kPdgPiP
const int kPdgPiP
Definition: PDGCodes.h:132
genie::kPdgPi0
const int kPdgPi0
Definition: PDGCodes.h:134
pINFO
#define pINFO
Definition: Messenger.h:53
pWARN
#define pWARN
Definition: Messenger.h:51
genie::Intranuke2014::fNucCEXFac
double fNucCEXFac
charge exchange xsec correction factor (hN Mode)
Definition: Intranuke2014.h:109
genie::GHepParticle::KinE
double KinE(bool mass_from_pdg=false) const
Get kinetic energy.
Definition: GHepParticle.cxx:177
genie::INukeHadroFates::AsString
static string AsString(INukeFateHN_t fate)
Definition: INukeHadroFates.h:68
genie::kPdgPiM
const int kPdgPiM
Definition: PDGCodes.h:133
genie::kPdgProton
const int kPdgProton
Definition: PDGCodes.h:62
genie::controls::kRjMaxIterations
static const unsigned int kRjMaxIterations
Definition: Controls.h:27
genie::kPdgNeutron
const int kPdgNeutron
Definition: PDGCodes.h:64
genie::Intranuke2014::fNucAbsFac
double fNucAbsFac
absorption xsec correction factor (hN Mode)
Definition: Intranuke2014.h:108
plot_model.r
r
Definition: plot_model.py:48
pDEBUG
#define pDEBUG
Definition: Messenger.h:54
int HNIntranuke2014::HandleCompoundNucleus ( GHepRecord *  ev,
GHepParticle *  p,
int  mom 
) const
privatevirtual

Implements genie::Intranuke2014.

Definition at line 857 of file HNIntranuke2014.cxx.

858 {
859 
860  // handle compound nucleus option
861  // -- Call the PreEquilibrium function
862  if( fDoCompoundNucleus && IsInNucleus(p) && pdg::IsNeutronOrProton(p->Pdg()))
863  { // random number generator
864  RandomGen * rnd = RandomGen::Instance();
865 
866  double rpreeq = rnd->RndFsi().Rndm(); // sdytman test
867 
868  if((p->KinE() < fEPreEq) )
869  {
870  if(fRemnA>5&&rpreeq<0.12)
871  {
872  GHepParticle * sp = new GHepParticle(*p);
873  sp->SetFirstMother(mom);
874  utils::intranuke2014::PreEquilibrium(ev,sp,fRemnA,fRemnZ,fRemnP4,
875  fDoFermi,fFermiFac,fNuclmodel,fNucRmvE,kIMdHN);
876  delete sp;
877  return 2;
878  }
879  else
880  {
881  // nothing left to interact with!
882  LOG("HNIntranuke2014", pNOTICE)
883  << "*** Nothing left to interact with, escaping.";
884  GHepParticle * sp = new GHepParticle(*p);
885  sp->SetFirstMother(mom);
886  sp->SetStatus(kIStStableFinalState);
887  ev->AddParticle(*sp);
888  delete sp;
889  return 1;
890  }
891  }
892  }
893  return 0;
894 }
genie::GHepParticle::SetFirstMother
void SetFirstMother(int m)
Definition: GHepParticle.h:131
genie::RandomGen::RndFsi
TRandom3 & RndFsi(void) const
rnd number generator used by intranuclear cascade monte carlos
Definition: RandomGen.h:60
genie::RandomGen::Instance
static RandomGen * Instance()
Access instance.
Definition: RandomGen.cxx:79
genie::Intranuke2014::fEPreEq
double fEPreEq
threshold for pre-equilibrium reaction
Definition: Intranuke2014.h:110
genie::Intranuke2014::fRemnZ
int fRemnZ
remnant nucleus Z
Definition: Intranuke2014.h:97
genie::Intranuke2014::IsInNucleus
bool IsInNucleus(const GHepParticle *p) const
Definition: Intranuke2014.cxx:250
genie::RandomGen
A singleton holding random number generator classes. All random number generation in GENIE should tak...
Definition: RandomGen.h:30
genie::GHepParticle::Pdg
int Pdg(void) const
Definition: GHepParticle.h:64
genie::Intranuke2014::fRemnA
int fRemnA
remnant nucleus A
Definition: Intranuke2014.h:96
genie::Intranuke2014::fNuclmodel
const NuclearModelI * fNuclmodel
nuclear model used to generate fermi momentum
Definition: Intranuke2014.h:95
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::utils::intranuke2014::PreEquilibrium
void PreEquilibrium(GHepRecord *ev, GHepParticle *p, int &RemnA, int &RemnZ, TLorentzVector &RemnP4, bool DoFermi, double FermiFac, const NuclearModelI *Nuclmodel, double NucRmvE, EINukeMode mode=kIMdHN)
Definition: INukeUtils2014.cxx:398
genie::Intranuke2014::fFermiFac
double fFermiFac
testing parameter to modify fermi momentum
Definition: Intranuke2014.h:111
genie::Intranuke2014::fNucRmvE
double fNucRmvE
binding energy to subtract from cascade nucleons
Definition: Intranuke2014.h:104
genie::Intranuke2014::fDoCompoundNucleus
bool fDoCompoundNucleus
whether or not to do compound nucleus considerations
Definition: Intranuke2014.h:116
genie::GHepParticle::KinE
double KinE(bool mass_from_pdg=false) const
Get kinetic energy.
Definition: GHepParticle.cxx:177
genie::GHepParticle::SetStatus
void SetStatus(GHepStatus_t s)
Definition: GHepParticle.h:125
genie::pdg::IsNeutronOrProton
bool IsNeutronOrProton(int pdgc)
Definition: PDGUtils.cxx:321
genie::GHepRecord::AddParticle
virtual void AddParticle(const GHepParticle &p)
Definition: GHepRecord.cxx:526
pNOTICE
#define pNOTICE
Definition: Messenger.h:52
genie::GHepParticle
STDHEP-like event record entry that can fit a particle or a nucleus.
Definition: GHepParticle.h:40
genie::Intranuke2014::fRemnP4
TLorentzVector fRemnP4
P4 of remnant system.
Definition: Intranuke2014.h:98
genie::Intranuke2014::fDoFermi
bool fDoFermi
whether or not to do fermi mom.
Definition: Intranuke2014.h:114
bool HNIntranuke2014::HandleCompoundNucleusHN ( GHepRecord *  ev,
GHepParticle *  p 
) const
private

Definition at line 896 of file HNIntranuke2014.cxx.

897 {
898  return (this->HandleCompoundNucleus(ev,p,p->FirstMother())==2);
899 }
genie::GHepParticle::FirstMother
int FirstMother(void) const
Definition: GHepParticle.h:67
genie::HNIntranuke2014::HandleCompoundNucleus
int HandleCompoundNucleus(GHepRecord *ev, GHepParticle *p, int mom) const
Definition: HNIntranuke2014.cxx:857
void HNIntranuke2014::InelasticHN ( GHepRecord *  ev,
GHepParticle *  p 
) const
private

Definition at line 712 of file HNIntranuke2014.cxx.

713 {
714  // Aaron Meyer (Jan 2010)
715  // Updated version of InelasticHN
716 
717  GHepParticle * s1 = new GHepParticle(*p);
718  GHepParticle * s2 = new GHepParticle(*p);
719  GHepParticle * s3 = new GHepParticle(*p);
720 
721  if (utils::intranuke2014::PionProduction(ev,p,s1,s2,s3,fRemnA,fRemnZ,fRemnP4,fDoFermi,fFermiFac,fFermiMomentum,fNuclmodel))
722  {
723  // set status of particles and return
724 
725  s1->SetStatus(kIStHadronInTheNucleus);
726  s2->SetStatus(kIStHadronInTheNucleus);
727  s3->SetStatus(kIStHadronInTheNucleus);
728 
729  ev->AddParticle(*s1);
730  ev->AddParticle(*s2);
731  ev->AddParticle(*s3);
732  }
733  else
734  {
735 
736  delete s1; //fixes potential memory leak
737  delete s2;
738  delete s3;
739 
740  LOG("HNIntranuke2014", pNOTICE) << "Error: could not create pion production final state";
741  exceptions::INukeException exception;
742  exception.SetReason("PionProduction in hN failed");
743  throw exception;
744  }
745 
746  delete s1;
747  delete s2;
748  delete s3;
749  return;
750 
751 }
genie::exceptions::INukeException
An exception thrown by SimulateHadronState for kinematics problems. TwoBodyCollision/Kinematics used ...
Definition: INukeException.h:38
genie::Intranuke2014::fRemnZ
int fRemnZ
remnant nucleus Z
Definition: Intranuke2014.h:97
genie::Intranuke2014::fFermiMomentum
double fFermiMomentum
whether or not particle collision is pauli blocked
Definition: Intranuke2014.h:113
genie::Intranuke2014::fRemnA
int fRemnA
remnant nucleus A
Definition: Intranuke2014.h:96
genie::Intranuke2014::fNuclmodel
const NuclearModelI * fNuclmodel
nuclear model used to generate fermi momentum
Definition: Intranuke2014.h:95
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::utils::intranuke2014::PionProduction
bool PionProduction(GHepRecord *ev, GHepParticle *p, GHepParticle *s1, GHepParticle *s2, GHepParticle *s3, int &RemnA, int &RemnZ, TLorentzVector &RemnP4, bool DoFermi, double FermiFac, double FermiMomentum, const NuclearModelI *Nuclmodel)
Definition: INukeUtils2014.cxx:1180
genie::Intranuke2014::fFermiFac
double fFermiFac
testing parameter to modify fermi momentum
Definition: Intranuke2014.h:111
genie::GHepParticle::SetStatus
void SetStatus(GHepStatus_t s)
Definition: GHepParticle.h:125
genie::GHepRecord::AddParticle
virtual void AddParticle(const GHepParticle &p)
Definition: GHepRecord.cxx:526
pNOTICE
#define pNOTICE
Definition: Messenger.h:52
genie::GHepParticle
STDHEP-like event record entry that can fit a particle or a nucleus.
Definition: GHepParticle.h:40
genie::Intranuke2014::fRemnP4
TLorentzVector fRemnP4
P4 of remnant system.
Definition: Intranuke2014.h:98
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
genie::Intranuke2014::fDoFermi
bool fDoFermi
whether or not to do fermi mom.
Definition: Intranuke2014.h:114
void HNIntranuke2014::LoadConfig ( void  )
privatevirtual

Implements genie::Intranuke2014.

Definition at line 901 of file HNIntranuke2014.cxx.

902 {
903  AlgConfigPool * confp = AlgConfigPool::Instance();
904  const Registry * gc = confp->GlobalParameterList();
905 
906  // load hadronic cross sections
907  fHadroData2014 = INukeHadroData2014::Instance();
908 
909  // fermi momentum setup
910  fAlgf = AlgFactory::Instance();
911  fNuclmodel = dynamic_cast<const NuclearModelI *>
912  (fAlgf->GetAlgorithm("genie::FGMBodekRitchie","Default"));
913 
914  // other intranuke config params
915  fR0 = fConfig->GetDoubleDef ("R0", gc->GetDouble("NUCL-R0")); // fm
916  fNR = fConfig->GetDoubleDef ("NR", gc->GetDouble("NUCL-NR"));
917  fNucRmvE = fConfig->GetDoubleDef ("NucRmvE", gc->GetDouble("INUKE-NucRemovalE")); // GeV
918  fDelRPion = fConfig->GetDoubleDef ("DelRPion", gc->GetDouble("HNINUKE-DelRPion"));
919  fDelRNucleon = fConfig->GetDoubleDef ("DelRNucleon", gc->GetDouble("HNINUKE-DelRNucleon"));
920  fHadStep = fConfig->GetDoubleDef ("HadStep", gc->GetDouble("INUKE-HadStep")); // fm
921  fNucAbsFac = fConfig->GetDoubleDef ("NucAbsFac", gc->GetDouble("INUKE-NucAbsFac"));
922  fNucQEFac = fConfig->GetDoubleDef ("NucQEFac", gc->GetDouble("INUKE-NucQEFac"));
923  fNucCEXFac = fConfig->GetDoubleDef ("NucCEXFac", gc->GetDouble("INUKE-NucCEXFac"));
924  fEPreEq = fConfig->GetDoubleDef ("EPreEq", gc->GetDouble("INUKE-Energy_Pre_Eq"));
925  fFermiFac = fConfig->GetDoubleDef ("FermiFac", gc->GetDouble("INUKE-FermiFac"));
926  fFermiMomentum = fConfig->GetDoubleDef ("FermiMomentum",gc->GetDouble("INUKE-FermiMomentum"));
927  fDoFermi = fConfig->GetBoolDef ("DoFermi", gc->GetBool("INUKE-DoFermi"));
928  fFreeStep = fConfig->GetDoubleDef ("FreeStep", gc->GetDouble("INUKE-FreeStep"));
929  fDoCompoundNucleus = fConfig->GetBoolDef ("DoCompoundNucleus", gc->GetBool("INUKE-DoCompoundNucleus"));
930 
931 
932  // report
933  LOG("HNIntranuke2014", pINFO) << "Settings for Intranuke2014 mode: " << INukeMode::AsString(kIMdHN);
934  LOG("HNIntranuke2014", pWARN) << "R0 = " << fR0 << " fermi";
935  LOG("HNIntranuke2014", pWARN) << "NR = " << fNR;
936  LOG("HNIntranuke2014", pWARN) << "DelRPion = " << fDelRPion;
937  LOG("HNIntranuke2014", pWARN) << "DelRNucleon = " << fDelRNucleon;
938  LOG("HNIntranuke2014", pWARN) << "HadStep = " << fHadStep << " fermi";
939  LOG("HNIntranuke2014", pWARN) << "NucAbsFac = " << fNucAbsFac;
940  LOG("HNIntranuke2014", pWARN) << "NucQEFac = " << fNucQEFac;
941  LOG("HNIntranuke2014", pWARN) << "NucCEXFac = " << fNucCEXFac;
942  LOG("HNIntranuke2014", pWARN) << "EPreEq = " << fEPreEq;
943  LOG("HNIntranuke2014", pWARN) << "FermiFac = " << fFermiFac;
944  LOG("HNIntranuke2014", pWARN) << "FreeStep = " << fFreeStep; // free step in fm
945  LOG("HNIntranuke2014", pWARN) << "FermiMomtm = " << fFermiMomentum;
946  LOG("HNIntranuke2014", pWARN) << "DoFermi? = " << ((fDoFermi)?(true):(false));
947  LOG("HNIntranuke2014", pWARN) << "DoCmpndNuc? = " << ((fDoCompoundNucleus)?(true):(false));
948 }
genie::INukeMode::AsString
static string AsString(INukeMode_t mode)
Definition: INukeMode.h:41
genie::Registry::GetDoubleDef
RgDbl GetDoubleDef(RgKey key, RgDbl def_opt, bool set_def=true)
Definition: Registry.cxx:549
genie::Intranuke2014::fAlgf
AlgFactory * fAlgf
algorithm factory instance
Definition: Intranuke2014.h:94
genie::Intranuke2014::fNR
double fNR
param multiplying the nuclear radius, determining how far to track hadrons beyond the "nuclear bounda...
Definition: Intranuke2014.h:103
genie::AlgConfigPool
A singleton class holding all configuration registries built while parsing all loaded XML configurati...
Definition: AlgConfigPool.h:41
genie::Intranuke2014::fEPreEq
double fEPreEq
threshold for pre-equilibrium reaction
Definition: Intranuke2014.h:110
genie::Intranuke2014::fR0
double fR0
effective nuclear size param
Definition: Intranuke2014.h:102
genie::NuclearModelI
Pure abstract base class. Defines the NuclearModelI interface to be implemented by any physics model ...
Definition: NuclearModelI.h:42
genie::Intranuke2014::fHadroData2014
INukeHadroData2014 * fHadroData2014
a collection of h+N,h+A data & calculations
Definition: Intranuke2014.h:93
genie::Registry::GetDouble
RgDbl GetDouble(RgKey key) const
Definition: Registry.cxx:488
genie::Intranuke2014::fFermiMomentum
double fFermiMomentum
whether or not particle collision is pauli blocked
Definition: Intranuke2014.h:113
genie::Intranuke2014::fDelRPion
double fDelRPion
factor by which Pion Compton wavelength gets multiplied to become nuclear size enhancement ...
Definition: Intranuke2014.h:105
genie::Intranuke2014::fFreeStep
double fFreeStep
produced particle free stem, in fm
Definition: Intranuke2014.h:112
genie::Intranuke2014::fNuclmodel
const NuclearModelI * fNuclmodel
nuclear model used to generate fermi momentum
Definition: Intranuke2014.h:95
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::AlgFactory::GetAlgorithm
const Algorithm * GetAlgorithm(const AlgId &algid)
Definition: AlgFactory.cxx:86
pINFO
#define pINFO
Definition: Messenger.h:53
genie::AlgConfigPool::GlobalParameterList
Registry * GlobalParameterList(void) const
Definition: AlgConfigPool.cxx:495
pWARN
#define pWARN
Definition: Messenger.h:51
genie::Intranuke2014::fFermiFac
double fFermiFac
testing parameter to modify fermi momentum
Definition: Intranuke2014.h:111
genie::Intranuke2014::fNucCEXFac
double fNucCEXFac
charge exchange xsec correction factor (hN Mode)
Definition: Intranuke2014.h:109
genie::Intranuke2014::fNucRmvE
double fNucRmvE
binding energy to subtract from cascade nucleons
Definition: Intranuke2014.h:104
genie::Intranuke2014::fDoCompoundNucleus
bool fDoCompoundNucleus
whether or not to do compound nucleus considerations
Definition: Intranuke2014.h:116
genie::AlgFactory::Instance
static AlgFactory * Instance()
Definition: AlgFactory.cxx:75
genie::Registry
A registry. Provides the container for algorithm configuration parameters.
Definition: Registry.h:66
genie::Registry::GetBool
RgBool GetBool(RgKey key) const
Definition: Registry.cxx:474
genie::Algorithm::fConfig
Registry * fConfig
config. (either owned or pointing to config pool)
Definition: Algorithm.h:122
genie::Intranuke2014::fDelRNucleon
double fDelRNucleon
factor by which Nucleon Compton wavelength gets multiplied to become nuclear size enhancement ...
Definition: Intranuke2014.h:106
genie::Registry::GetBoolDef
RgBool GetBoolDef(RgKey key, RgBool def_opt, bool set_def=true)
Definition: Registry.cxx:539
genie::Intranuke2014::fNucAbsFac
double fNucAbsFac
absorption xsec correction factor (hN Mode)
Definition: Intranuke2014.h:108
genie::Intranuke2014::fHadStep
double fHadStep
step size for intranuclear hadron transport
Definition: Intranuke2014.h:107
genie::AlgConfigPool::Instance
static AlgConfigPool * Instance()
Definition: AlgConfigPool.cxx:98
genie::INukeHadroData2014::Instance
static INukeHadroData2014 * Instance(void)
Definition: INukeHadroData2014.cxx:149
genie::Intranuke2014::fDoFermi
bool fDoFermi
whether or not to do fermi mom.
Definition: Intranuke2014.h:114
void HNIntranuke2014::ProcessEventRecord ( GHepRecord *  event_rec) const
virtual

Reimplemented from genie::Intranuke2014.

Definition at line 104 of file HNIntranuke2014.cxx.

105 {
106  LOG("HNIntranuke2014", pNOTICE)
107  << "************ Running hN2014 MODE INTRANUKE ************";
108 
109  LOG("HNIntranuke2014", pWARN)
110  << print::PrintFramedMesg(
111  "Experimental code (INTRANUKE/hN model) - Run at your own risk");
112 
113  Intranuke2014::ProcessEventRecord(evrec);
114 
115  LOG("HNIntranuke2014", pINFO) << "Done with this event";
116 }
genie::Intranuke2014::ProcessEventRecord
virtual void ProcessEventRecord(GHepRecord *event_rec) const
Definition: Intranuke2014.cxx:116
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
pINFO
#define pINFO
Definition: Messenger.h:53
pWARN
#define pWARN
Definition: Messenger.h:51
genie::utils::print::PrintFramedMesg
string PrintFramedMesg(string mesg, unsigned int nl=1, const char f='*')
Definition: PrintUtils.cxx:171
pNOTICE
#define pNOTICE
Definition: Messenger.h:52
void HNIntranuke2014::SimulateHadronicFinalState ( GHepRecord *  ev,
GHepParticle *  p 
) const
privatevirtual

Implements genie::Intranuke2014.

Definition at line 118 of file HNIntranuke2014.cxx.

119 {
120 // Simulate a hadron interaction for the input particle p in HN mode
121 //
122  if(!p || !ev)
123  {
124  LOG("HNIntranuke2014", pERROR) << "** Null input!";
125  return;
126  }
127 
128  // check particle id
129  int pdgc = p->Pdg();
130  bool is_pion = (pdgc==kPdgPiP || pdgc==kPdgPiM || pdgc==kPdgPi0);
131  // bool is_kaon = (pdgc==kPdgKP); // UNUSED - comment to quiet compiler warnings
132  bool is_baryon = (pdgc==kPdgProton || pdgc==kPdgNeutron);
133  bool is_gamma = (pdgc==kPdgGamma);
134  if(!(is_pion || is_baryon || is_gamma))
135  {
136  LOG("HNIntranuke2014", pERROR) << "** Cannot handle particle: " << p->Name();
137  return;
138  }
139  try
140  {
141  // select a fate for the input particle
142  INukeFateHN_t fate = this->HadronFateHN(p);
143 
144  // store the fate
145  ev->Particle(p->FirstMother())->SetRescatterCode((int)fate);
146 
147  if(fate == kIHNFtUndefined)
148  {
149  LOG("HNIntranuke2014", pERROR) << "** Couldn't select a fate";
150  LOG("HNIntranuke2014", pERROR) << "** Num Protons: " << fRemnZ
151  << ", Num Neutrons: "<<(fRemnA-fRemnZ);
152  LOG("HNIntranuke2014", pERROR) << "** Particle: " << "\n" << (*p);
153  //LOG("HNIntranuke2014", pERROR) << "** Event Record: " << "\n" << (*ev);
154  //p->SetStatus(kIStUndefined);
155  p->SetStatus(kIStStableFinalState);
156  ev->AddParticle(*p);
157  return;
158  }
159 
160  LOG("HNIntranuke2014", pNOTICE)
161  << "Selected " << p->Name() << " fate: " << INukeHadroFates::AsString(fate);
162 
163  // handle the reaction
164  if(fate == kIHNFtCEx || fate == kIHNFtElas)
165  {
166  this->ElasHN(ev,p,fate);
167  }
168  else if(fate == kIHNFtAbs) {this-> AbsorbHN(ev,p,fate);}
169  else if(fate == kIHNFtInelas && pdgc != kPdgGamma)
170  {
171 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
172  LOG("HNIntranuke2014", pDEBUG)
173  << "Invoking InelasticHN() for a : " << p->Name()
174  << " whose fate is : " << INukeHadroFates::AsString(fate);
175 #endif
176 
177  this-> InelasticHN(ev,p);
178  }
179  else if(fate == kIHNFtInelas && pdgc == kPdgGamma) {this-> GammaInelasticHN(ev,p,fate);}
180  else if(fate == kIHNFtNoInteraction)
181  {
182  p->SetStatus(kIStStableFinalState);
183  ev->AddParticle(*p);
184  return;
185  }
186  else if(fate == kIHNFtCmp){utils::intranuke2014::PreEquilibrium(ev,p,fRemnA,fRemnZ,fRemnP4,fDoFermi,fFermiFac,fNuclmodel,fNucRmvE,kIMdHN);}
187 
188 
189  }
190  catch(exceptions::INukeException exception)
191  {
192  this->SimulateHadronicFinalState(ev,p);
193  LOG("HNIntranuke2014", pNOTICE)
194  << "retry call to SimulateHadronicFinalState ";
195  LOG("HNIntranuke2014", pNOTICE) << exception;
196 
197  }
198 }
genie::exceptions::INukeException
An exception thrown by SimulateHadronState for kinematics problems. TwoBodyCollision/Kinematics used ...
Definition: INukeException.h:38
genie::GHepRecord::Particle
virtual GHepParticle * Particle(int position) const
Definition: GHepRecord.cxx:148
genie::HNIntranuke2014::SimulateHadronicFinalState
void SimulateHadronicFinalState(GHepRecord *ev, GHepParticle *p) const
Definition: HNIntranuke2014.cxx:118
pERROR
#define pERROR
Definition: Messenger.h:50
genie::Intranuke2014::fRemnZ
int fRemnZ
remnant nucleus Z
Definition: Intranuke2014.h:97
genie::GHepParticle::Pdg
int Pdg(void) const
Definition: GHepParticle.h:64
genie::GHepParticle::FirstMother
int FirstMother(void) const
Definition: GHepParticle.h:67
genie::GHepParticle::Name
string Name(void) const
Name that corresponds to the PDG code.
Definition: GHepParticle.cxx:153
genie::Intranuke2014::fRemnA
int fRemnA
remnant nucleus A
Definition: Intranuke2014.h:96
genie::HNIntranuke2014::InelasticHN
void InelasticHN(GHepRecord *ev, GHepParticle *p) const
Definition: HNIntranuke2014.cxx:712
genie::Intranuke2014::fNuclmodel
const NuclearModelI * fNuclmodel
nuclear model used to generate fermi momentum
Definition: Intranuke2014.h:95
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:87
genie::utils::intranuke2014::PreEquilibrium
void PreEquilibrium(GHepRecord *ev, GHepParticle *p, int &RemnA, int &RemnZ, TLorentzVector &RemnP4, bool DoFermi, double FermiFac, const NuclearModelI *Nuclmodel, double NucRmvE, EINukeMode mode=kIMdHN)
Definition: INukeUtils2014.cxx:398
genie::kPdgGamma
const int kPdgGamma
Definition: PDGCodes.h:163
genie::kPdgPiP
const int kPdgPiP
Definition: PDGCodes.h:132
genie::kPdgPi0
const int kPdgPi0
Definition: PDGCodes.h:134
genie::INukeFateHN_t
enum genie::EINukeFateHN_t INukeFateHN_t
genie::Intranuke2014::fFermiFac
double fFermiFac
testing parameter to modify fermi momentum
Definition: Intranuke2014.h:111
genie::Intranuke2014::fNucRmvE
double fNucRmvE
binding energy to subtract from cascade nucleons
Definition: Intranuke2014.h:104
genie::HNIntranuke2014::GammaInelasticHN
void GammaInelasticHN(GHepRecord *ev, GHepParticle *p, INukeFateHN_t fate) const
Definition: HNIntranuke2014.cxx:753
genie::INukeHadroFates::AsString
static string AsString(INukeFateHN_t fate)
Definition: INukeHadroFates.h:68
genie::HNIntranuke2014::AbsorbHN
void AbsorbHN(GHepRecord *ev, GHepParticle *p, INukeFateHN_t fate) const
Definition: HNIntranuke2014.cxx:335
genie::kPdgPiM
const int kPdgPiM
Definition: PDGCodes.h:133
genie::GHepRecord::AddParticle
virtual void AddParticle(const GHepParticle &p)
Definition: GHepRecord.cxx:526
genie::kPdgProton
const int kPdgProton
Definition: PDGCodes.h:62
pNOTICE
#define pNOTICE
Definition: Messenger.h:52
genie::HNIntranuke2014::ElasHN
void ElasHN(GHepRecord *ev, GHepParticle *p, INukeFateHN_t fate) const
Definition: HNIntranuke2014.cxx:598
genie::kPdgNeutron
const int kPdgNeutron
Definition: PDGCodes.h:64
genie::HNIntranuke2014::HadronFateHN
INukeFateHN_t HadronFateHN(const GHepParticle *p) const
Definition: HNIntranuke2014.cxx:200
genie::Intranuke2014::fRemnP4
TLorentzVector fRemnP4
P4 of remnant system.
Definition: Intranuke2014.h:98
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
pDEBUG
#define pDEBUG
Definition: Messenger.h:54
genie::Intranuke2014::fDoFermi
bool fDoFermi
whether or not to do fermi mom.
Definition: Intranuke2014.h:114

Friends And Related Function Documentation

friend class IntranukeTester
friend

Definition at line 53 of file HNIntranuke2014.h.

Member Data Documentation

double genie::HNIntranuke2014::fNucQEFac
private

Definition at line 80 of file HNIntranuke2014.h.

int genie::HNIntranuke2014::nuclA
mutableprivate

value of A for the target nucleus in hA mode

Definition at line 77 of file HNIntranuke2014.h.

The documentation for this class was generated from the following files:
Generated by   doxygen 1.8.11