genie::Intranuke2015 Class Reference

#include <Intranuke2015.h>

Inheritance diagram for genie::Intranuke2015:

Public Member Functions
	Intranuke2015 ()
	Intranuke2015 (string name)
	Intranuke2015 (string name, string config)
	~Intranuke2015 ()
virtual void	ProcessEventRecord (GHepRecord *event_rec) const
void	Configure (const Registry &config)
	Configure the algorithm. More...
void	Configure (string param_set)
	Configure the algorithm. More...
virtual string	GetINukeMode () const
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...

Protected Member Functions
virtual void	LoadConfig (void)=0
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
virtual void	SimulateHadronicFinalState (GHepRecord *ev, GHepParticle *p) const =0
virtual int	HandleCompoundNucleus (GHepRecord *ev, GHepParticle *p, int mom) const =0
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
double	fTrackingRadius
	tracking radius for the nucleus in the current event More...
TGenPhaseSpace	fGenPhaseSpace
	a phase space generator More...
INukeHadroData2015 *	fHadroData2015
	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...
bool	fUseOset
	Oset model for low energy pion in hN. More...
bool	fAltOset
	NuWro's table-based implementation (not recommended) More...
bool	fXsecNNCorr
	use nuclear medium correction for NN cross section 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...

Friends
class	IntranukeTester

Detailed Description

Definition at line 54 of file Intranuke2015.h.

Constructor & Destructor Documentation

Intranuke2015::Intranuke2015

(

)

Definition at line 93 of file Intranuke2015.cxx.

                             :
EventRecordVisitorI()
{

}

Intranuke2015::Intranuke2015

(

string

name

)

Definition at line 99 of file Intranuke2015.cxx.

                                        :
EventRecordVisitorI(name)
{

}

Intranuke2015::Intranuke2015	(	string	name,
		string	config
	)

Definition at line 105 of file Intranuke2015.cxx.

                                                       :
EventRecordVisitorI(name, config)
{

}

Intranuke2015::~Intranuke2015

(

)

Definition at line 111 of file Intranuke2015.cxx.

{

}

Member Function Documentation

bool Intranuke2015::CanRescatter ( const GHepParticle *  p ) const

protected

Definition at line 233 of file Intranuke2015.cxx.

{
// checks whether a particle that needs to be rescattered, can in fact be 
// rescattered by this cascade MC

  assert(p);
  return  ( p->Pdg() == kPdgPiP     || 
            p->Pdg() == kPdgPiM     || 
            p->Pdg() == kPdgPi0     ||
            p->Pdg() == kPdgProton  ||
            p->Pdg() == kPdgNeutron ||
            //      p->Pdg() == kPdgGamma   ||
            p->Pdg() == kPdgKP      //||
            //      p->Pdg() == kPdgKM
          );
}

void Intranuke2015::Configure ( const Registry &  config )

virtual

Configure the algorithm.

Reimplemented from genie::Algorithm.

Definition at line 414 of file Intranuke2015.cxx.

{
  Algorithm::Configure(config);
  this->LoadConfig();
}

void Intranuke2015::Configure ( string  config )

virtual

Configure the algorithm.

Reimplemented from genie::Algorithm.

Definition at line 420 of file Intranuke2015.cxx.

{
  Algorithm::Configure(param_set);
  this->LoadConfig();
}

double Intranuke2015::GenerateStep ( GHepRecord *  ev, GHepParticle *  p  ) const

protected

Definition at line 394 of file Intranuke2015.cxx.

{
// Generate a step (in fermis) for particle p in the input event.
// Computes the mean free path L and generate an 'interaction' distance d 
// from an exp(-d/L) distribution

  RandomGen * rnd = RandomGen::Instance();

  string fINukeMode = this->GetINukeMode();
  
  double L = utils::intranuke2015::MeanFreePath(p->Pdg(), *p->X4(), *p->P4(), fRemnA,
                                                fRemnZ, fDelRPion, fDelRNucleon, fUseOset, fAltOset, fXsecNNCorr, fINukeMode);
  double d = -1.*L * TMath::Log(rnd->RndFsi().Rndm());

  LOG("Intranuke2015", pDEBUG)
            << "Mean free path = " << L << " fm / "
                              << "Generated path length = " << d << " fm";
  return d;
}

void Intranuke2015::GenerateVertex ( GHepRecord *  ev ) const

protected

Definition at line 153 of file Intranuke2015.cxx.

{
// Sets a vertex in the nucleus periphery
// Called onlt in hadron/photon-nucleus interactions.

  GHepParticle * nucltgt = evrec->TargetNucleus();
  assert(nucltgt);

  RandomGen * rnd = RandomGen::Instance();
  TVector3 vtx(999999.,999999.,999999.);

  // *** For h+A events (test mode): 
  // Assume a hadron beam with uniform intensity across an area, 
  // so we need to choose events uniformly within that area.
  double x=999999., y=999999., epsilon = 0.001;
  double R2  = TMath::Power(fTrackingRadius,2.);
  double rp2 = TMath::Power(x,2.) + TMath::Power(y,2.);
  while(rp2 > R2-epsilon) {
      x = (fTrackingRadius-epsilon) * rnd->RndFsi().Rndm();
      y = -fTrackingRadius + 2*fTrackingRadius * rnd->RndFsi().Rndm();
      y -= ((y>0) ? epsilon : -epsilon);
      rp2 = TMath::Power(x,2.) + TMath::Power(y,2.);
  }
  vtx.SetXYZ(x,y, -1.*TMath::Sqrt(TMath::Max(0.,R2-rp2)) + epsilon);

  // get the actual unit vector along the incoming hadron direction
  TVector3 direction = evrec->Probe()->P4()->Vect().Unit();

  // rotate the vtx position
  vtx.RotateUz(direction);
  
  LOG("Intranuke2015", pNOTICE) 
     << "Generated vtx @ R = " << vtx.Mag() << " fm / " 
     << print::Vec3AsString(&vtx);

  TObjArrayIter piter(evrec);
  GHepParticle * p = 0;
  while( (p = (GHepParticle *) piter.Next()) )
  {
    if(pdg::IsPseudoParticle(p->Pdg())) continue;
    if(pdg::IsIon           (p->Pdg())) continue;

    p->SetPosition(vtx.x(), vtx.y(), vtx.z(), 0.);
  }
}

virtual string genie::Intranuke2015::GetINukeMode ( ) const

inlinevirtual

Reimplemented in genie::HAIntranuke2015, and genie::HNIntranuke2015.

Definition at line 72 of file Intranuke2015.h.

{return "XX2015";};

virtual int genie::Intranuke2015::HandleCompoundNucleus ( GHepRecord *  ev, GHepParticle *  p, int  mom  ) const

protectedpure virtual

Implemented in genie::HAIntranuke2015, and genie::HNIntranuke2015.

bool Intranuke2015::IsInNucleus ( const GHepParticle *  p ) const

protected

Definition at line 250 of file Intranuke2015.cxx.

{
// check whether the input particle is still within the nucleus
//
  return (p->X4()->Vect().Mag() < fTrackingRadius + fHadStep);
}

virtual void genie::Intranuke2015::LoadConfig ( void  )

protectedpure virtual

Implemented in genie::HAIntranuke2015, and genie::HNIntranuke2015.

bool Intranuke2015::NeedsRescattering ( const GHepParticle *  p ) const

protected

Definition at line 214 of file Intranuke2015.cxx.

{
// checks whether the particle should be rescattered

  assert(p);

  if(fGMode == kGMdHadronNucleus ||
     fGMode == kGMdPhotonNucleus) {
    // hadron/photon-nucleus scattering propagate the incoming particle
    return (
      (p->Status() == kIStInitialState || p->Status() == kIStHadronInTheNucleus)
       && !pdg::IsIon(p->Pdg()));
  }
  else {
   // attempt to rescatter anything marked as 'hadron in the nucleus'
   return (p->Status() == kIStHadronInTheNucleus);
  }
}

void Intranuke2015::ProcessEventRecord ( GHepRecord *  event_rec ) const

virtual

Implements genie::EventRecordVisitorI.

Reimplemented in genie::HAIntranuke2015, and genie::HNIntranuke2015.

Definition at line 116 of file Intranuke2015.cxx.

{
  // Do not continue if there is no nuclear target
  GHepParticle * nucltgt = evrec->TargetNucleus();
  if (!nucltgt) {
    LOG("HNIntranuke2015", pINFO) << "No nuclear target found - INTRANUKE exits";
    return;
  }

  // Decide tracking radius for the current nucleus (few * R0 * A^1/3)
  this->SetTrackingRadius(nucltgt);

  // Understand what the event generation mode is (hadron/photon-nucleus,
  // lepton-nucleus, nucleon decay) from the input event.
  // The determined mode has an effect on INTRANUKE behaviour (how to lookup
  // the residual nucleus, whether to set an intranuclear vtx etc) but it
  // does not affect the INTRANUKE physics.
  fGMode = evrec->EventGenerationMode();

  // For lepton-nucleus scattering and for nucleon decay intranuclear vtx 
  // position (in the target nucleus coord system) is set elsewhere.
  // This method only takes effect in hadron/photon-nucleus interactions.
  // In this special mode, an interaction vertex is set at the periphery 
  // of the target nucleus.
  if(fGMode == kGMdHadronNucleus ||
     fGMode == kGMdPhotonNucleus)
  {
    this->GenerateVertex(evrec);  
  }

  // Now transport all hadrons outside the tracking radius.
  // Stepping part is common for both HA and HN.
  // Once it has been estabished that an interaction takes place then
  // HA and HN specific code takes over in order to simulate the final state.
  this->TransportHadrons(evrec);
}

void Intranuke2015::SetTrackingRadius ( const GHepParticle *  p ) const

protected

Definition at line 199 of file Intranuke2015.cxx.

{
  assert(p && pdg::IsIon(p->Pdg()));
  double A = p->A();
  fTrackingRadius = fR0 * TMath::Power(A, 1./3.);

  // multiply that by some input factor so that hadrons are tracked
  // beyond the nuclear 'boundary' since the nuclear density distribution
  // is not zero there
  fTrackingRadius *= fNR; 

  LOG("Intranuke2015", pNOTICE) 
      << "Setting tracking radius to R = " << fTrackingRadius;
}

virtual void genie::Intranuke2015::SimulateHadronicFinalState ( GHepRecord *  ev, GHepParticle *  p  ) const

protectedpure virtual

Implemented in genie::HNIntranuke2015, and genie::HAIntranuke2015.

void Intranuke2015::TransportHadrons ( GHepRecord *  ev ) const

protected

Definition at line 257 of file Intranuke2015.cxx.

{
// transport all hadrons outside the nucleus

  int inucl = -1;
  fRemnA = -1;
  fRemnZ = -1;

  //  Get 'nuclear environment' at the beginning of hadron transport 
  //  and keep track of the remnant nucleus A,Z  

  if(fGMode == kGMdHadronNucleus ||
     fGMode == kGMdPhotonNucleus)
  {
     inucl = evrec->TargetNucleusPosition();
  }
  else
  if(fGMode == kGMdLeptonNucleus || 
     fGMode == kGMdNucleonDecay) 
  {
     inucl = evrec->RemnantNucleusPosition();
  }

  LOG("Intranuke2015", pNOTICE) 
     << "Propagating hadrons within nucleus found in position = " << inucl;
  GHepParticle * nucl = evrec->Particle(inucl);
  if(!nucl) {
    LOG("Intranuke2015", pERROR) 
       << "No nucleus found in position = " << inucl;
    LOG("Intranuke2015", pERROR) 
       << *evrec;
    return;
  }
  
  fRemnA = nucl->A();
  fRemnZ = nucl->Z();

  LOG("Intranuke2015", pNOTICE)
      << "Nucleus (A,Z) = (" << fRemnA << ", " << fRemnZ << ")";

  const TLorentzVector & p4nucl = *(nucl->P4());
  fRemnP4 = p4nucl; 

  // Loop over GHEP and run intranuclear rescattering on handled particles
  TObjArrayIter piter(evrec);
  GHepParticle * p = 0;
  int icurr = -1;

  while( (p = (GHepParticle *) piter.Next()) )
  {
    icurr++;

    // Check whether the particle needs rescattering, otherwise skip it
    if( ! this->NeedsRescattering(p) ) continue;
 
    LOG("Intranuke2015", pNOTICE)
      << " >> Stepping a " << p->Name() 
                        << " with kinetic E = " << p->KinE() << " GeV";

    // Rescatter a clone, not the original particle
    GHepParticle * sp = new GHepParticle(*p); 

    // Set clone's mom to be the hadron that was cloned
    sp->SetFirstMother(icurr); 

    // Check whether the particle can be rescattered 
    if(!this->CanRescatter(sp)) {

       // if I can't rescatter it, I will just take it out of the nucleus
       LOG("Intranuke2015", pNOTICE)
              << "... Current version can't rescatter a " << sp->Name();
       sp->SetFirstMother(icurr); 
       sp->SetStatus(kIStStableFinalState);
       evrec->AddParticle(*sp);
       delete sp;
       continue; // <-- skip to next GHEP entry
    }

    // Start stepping particle out of the nucleus
    bool has_interacted = false;
    while ( this-> IsInNucleus(sp) ) 
    {
      // advance the hadron by a step
      utils::intranuke2015::StepParticle(sp, fHadStep);

      // check whether it interacts
      double d = this->GenerateStep(evrec,sp);
      has_interacted = (d<fHadStep);
      if(has_interacted) break;
    }//stepping

    //updating the position of the original particle with the position of the clone
    evrec->Particle(sp->FirstMother())->SetPosition(*(sp->X4()));
 
    if(has_interacted && fRemnA>0)  {
        // the particle interacts - simulate the hadronic interaction
      LOG("Intranuke2015", pNOTICE) 
          << "Particle has interacted at location:  " 
          << sp->X4()->Vect().Mag() << " / nucl rad= " << fTrackingRadius;
        this->SimulateHadronicFinalState(evrec,sp);
    } else if(has_interacted && fRemnA<=0) {
        // nothing left to interact with!
      LOG("Intranuke2015", pNOTICE)
          << "*** Nothing left to interact with, escaping.";
        sp->SetStatus(kIStStableFinalState);
        evrec->AddParticle(*sp);
        evrec->Particle(sp->FirstMother())->SetRescatterCode(1);
    } else {
        // the exits the nucleus without interacting - Done with it! 
        LOG("Intranuke2015", pNOTICE) 
          << "*** Hadron escaped the nucleus! Done with it.";
        sp->SetStatus(kIStStableFinalState);
        evrec->AddParticle(*sp);
        evrec->Particle(sp->FirstMother())->SetRescatterCode(1);
    }
    delete sp;

    // Current snapshot
    //LOG("Intranuke2015", pINFO) << "Current event record snapshot: " << *evrec;

  }// GHEP entries

  // Add remnant nucleus - that 'hadronic blob' has all the remaining hadronic 
  // 4p not  put explicitly into the simulated particles
  TLorentzVector v4(0.,0.,0.,0.);
  GHepParticle remnant_nucleus(
    kPdgHadronicBlob, kIStFinalStateNuclearRemnant, inucl,-1,-1,-1, fRemnP4, v4);
  evrec->AddParticle(remnant_nucleus);
  // Mark the initial remnant nucleus as an intermediate state 
  // Don't do that in the hadron/photon-nucleus scatterig mode since the initial 
  // remnant nucleus and the target nucleus coincide.
  if(fGMode != kGMdHadronNucleus &&
     fGMode != kGMdPhotonNucleus) {
     evrec->Particle(inucl)->SetStatus(kIStIntermediateState);
  }
}

Friends And Related Function Documentation

friend class IntranukeTester

friend

Definition at line 56 of file Intranuke2015.h.

Member Data Documentation

AlgFactory* genie::Intranuke2015::fAlgf

protected

algorithm factory instance

Definition at line 96 of file Intranuke2015.h.

bool genie::Intranuke2015::fAltOset

protected

NuWro's table-based implementation (not recommended)

Definition at line 120 of file Intranuke2015.h.

double genie::Intranuke2015::fDelRNucleon

protected

factor by which Nucleon Compton wavelength gets multiplied to become nuclear size enhancement

Definition at line 108 of file Intranuke2015.h.

double genie::Intranuke2015::fDelRPion

protected

factor by which Pion Compton wavelength gets multiplied to become nuclear size enhancement

Definition at line 107 of file Intranuke2015.h.

bool genie::Intranuke2015::fDoCompoundNucleus

protected

whether or not to do compound nucleus considerations

Definition at line 118 of file Intranuke2015.h.

bool genie::Intranuke2015::fDoFermi

protected

whether or not to do fermi mom.

Definition at line 116 of file Intranuke2015.h.

bool genie::Intranuke2015::fDoMassDiff

protected

whether or not to do mass diff. mode

Definition at line 117 of file Intranuke2015.h.

double genie::Intranuke2015::fEPreEq

protected

threshold for pre-equilibrium reaction

Definition at line 112 of file Intranuke2015.h.

double genie::Intranuke2015::fFermiFac

protected

testing parameter to modify fermi momentum

Definition at line 113 of file Intranuke2015.h.

double genie::Intranuke2015::fFermiMomentum

protected

whether or not particle collision is pauli blocked

Definition at line 115 of file Intranuke2015.h.

double genie::Intranuke2015::fFreeStep

protected

produced particle free stem, in fm

Definition at line 114 of file Intranuke2015.h.

TGenPhaseSpace genie::Intranuke2015::fGenPhaseSpace

mutableprotected

a phase space generator

Definition at line 94 of file Intranuke2015.h.

GEvGenMode_t genie::Intranuke2015::fGMode

mutableprotected

event generation mode (lepton+A, hadron+A, ...)

Definition at line 101 of file Intranuke2015.h.

INukeHadroData2015* genie::Intranuke2015::fHadroData2015

protected

a collection of h+N,h+A data & calculations

Definition at line 95 of file Intranuke2015.h.

double genie::Intranuke2015::fHadStep

protected

step size for intranuclear hadron transport

Definition at line 109 of file Intranuke2015.h.

double genie::Intranuke2015::fNR

protected

param multiplying the nuclear radius, determining how far to track hadrons beyond the "nuclear boundary"

Definition at line 105 of file Intranuke2015.h.

double genie::Intranuke2015::fNucAbsFac

protected

absorption xsec correction factor (hN Mode)

Definition at line 110 of file Intranuke2015.h.

double genie::Intranuke2015::fNucCEXFac

protected

charge exchange xsec correction factor (hN Mode)

Definition at line 111 of file Intranuke2015.h.

const NuclearModelI* genie::Intranuke2015::fNuclmodel

protected

nuclear model used to generate fermi momentum

Definition at line 97 of file Intranuke2015.h.

double genie::Intranuke2015::fNucRmvE

protected

binding energy to subtract from cascade nucleons

Definition at line 106 of file Intranuke2015.h.

double genie::Intranuke2015::fR0

protected

effective nuclear size param

Definition at line 104 of file Intranuke2015.h.

int genie::Intranuke2015::fRemnA

mutableprotected

remnant nucleus A

Definition at line 98 of file Intranuke2015.h.

TLorentzVector genie::Intranuke2015::fRemnP4

mutableprotected

P4 of remnant system.

Definition at line 100 of file Intranuke2015.h.

int genie::Intranuke2015::fRemnZ

mutableprotected

remnant nucleus Z

Definition at line 99 of file Intranuke2015.h.

double genie::Intranuke2015::fTrackingRadius

mutableprotected

tracking radius for the nucleus in the current event

Definition at line 93 of file Intranuke2015.h.

bool genie::Intranuke2015::fUseOset

protected

Oset model for low energy pion in hN.

Definition at line 119 of file Intranuke2015.h.

bool genie::Intranuke2015::fXsecNNCorr

protected

use nuclear medium correction for NN cross section

Definition at line 121 of file Intranuke2015.h.

---

The documentation for this class was generated from the following files:

• GENIE/src/HadronTransport/Intranuke2015.h
• GENIE/src/HadronTransport/Intranuke2015.cxx