An abstract class. It avoids implementing the HadronizationModelI interface, leaving it for the concrete subclasses (KNO, Pythia,...). It propagates some common methods to the concrete subclasses. More...

#include <HadronizationModelBase.h>

Inheritance diagram for genie::HadronizationModelBase:

Public Member Functions
virtual void	Initialize (void) const =0

virtual TClonesArray *	Hadronize (const Interaction *) const =0

virtual double	Weight (void) const =0

virtual PDGCodeList *	SelectParticles (const Interaction *) const =0

virtual TH1D *	MultiplicityProb (const Interaction , Option_t opt="") const =0

Public Member Functions inherited from genie::HadronizationModelI
virtual	~HadronizationModelI ()

Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()

virtual void	Configure (const Registry &config)

virtual void	Configure (string config)

virtual void	FindConfig (void)

virtual const Registry &	GetConfig (void) const

Registry *	GetOwnedConfig (void)

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

	HadronizationModelBase (string name)

	HadronizationModelBase (string name, string config)

virtual	~HadronizationModelBase ()

double	Wmin (void) const
	Various utility methods common to hadronization models. More...

double	MaxMult (const Interaction *i) const

void	ApplyRijk (const Interaction i, bool norm, TH1D mp) const

TH1D *	CreateMultProbHist (double maxmult) const

Protected Member Functions inherited from genie::HadronizationModelI
	HadronizationModelI ()

	HadronizationModelI (string name)

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

Registry *	ExtractLocalConfig (const Registry &in) const

Registry *	ExtractLowerConfig (const Registry &in, const string &alg_key) const
	Split an incoming configuration Registry into a block valid for the sub-algo identified by alg_key. More...

template<class T >
bool	GetParam (const RgKey &name, T &p, bool is_top_call=true) const

template<class T >
bool	GetParamDef (const RgKey &name, T &p, const T &def) const

template<class T >
bool	GetParamVect (const std::string &comm_name, std::vector< T > &v, unsigned int max, bool is_top_call=true) const

int	AddTopRegistry (Registry *rp, bool owns=true)
	add registry with top priority, also update ownership More...

int	AddLowRegistry (Registry *rp, bool owns=true)
	add registry with lowest priority, also update ownership More...

int	MergeTopRegistry (const Registry &r)

int	AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
	Add registries with top priority, also udated Ownerships. More...

Protected Attributes
double	fWcut
	configuration data common to all hadronizers More...

double	fRvpCCm2
	neugen's Rijk: vp, CC, multiplicity = 2 More...

double	fRvpCCm3
	neugen's Rijk: vp, CC, multiplicity = 3 More...

double	fRvpNCm2
	neugen's Rijk: vp, NC, multiplicity = 2 More...

double	fRvpNCm3
	neugen's Rijk: vp, NC, multiplicity = 3 More...

double	fRvnCCm2
	neugen's Rijk: vn, CC, multiplicity = 2 More...

double	fRvnCCm3
	neugen's Rijk: vn, CC, multiplicity = 3 More...

double	fRvnNCm2
	neugen's Rijk: vn, NC, multiplicity = 2 More...

double	fRvnNCm3
	neugen's Rijk: vn, NC, multiplicity = 3 More...

double	fRvbpCCm2
	neugen's Rijk: vbp, CC, multiplicity = 2 More...

double	fRvbpCCm3
	neugen's Rijk: vbp, CC, multiplicity = 3 More...

double	fRvbpNCm2
	neugen's Rijk: vbp, NC, multiplicity = 2 More...

double	fRvbpNCm3
	neugen's Rijk: vbp, NC, multiplicity = 3 More...

double	fRvbnCCm2
	neugen's Rijk: vbn, CC, multiplicity = 2 More...

double	fRvbnCCm3
	neugen's Rijk: vbn, CC, multiplicity = 3 More...

double	fRvbnNCm2
	neugen's Rijk: vbn, NC, multiplicity = 2 More...

double	fRvbnNCm3
	neugen's Rijk: vbn, NC, multiplicity = 3 More...

Protected Attributes inherited from genie::Algorithm
bool	fAllowReconfig

bool	fOwnsSubstruc
	true if it owns its substructure (sub-algs,...) More...

AlgId	fID
	algorithm name and configuration set More...

vector< Registry * >	fConfVect

vector< bool >	fOwnerships
	ownership for every registry in fConfVect 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

An abstract class. It avoids implementing the HadronizationModelI interface, leaving it for the concrete subclasses (KNO, Pythia,...). It propagates some common methods to the concrete subclasses.

Author: Costas Andreopoulos <costas.andreopoulos stfc.ac.uk> University of Liverpool & STFC Rutherford Appleton Lab

August 17, 2004

Definition at line 28 of file HadronizationModelBase.h.

Constructor & Destructor Documentation

HadronizationModelBase::HadronizationModelBase ( void )

protected

Definition at line 32 of file HadronizationModelBase.cxx.

                                                    :
 HadronizationModelI()
 {
 
 }

HadronizationModelBase::HadronizationModelBase ( string name )

protected

Definition at line 38 of file HadronizationModelBase.cxx.

                                                           :
 HadronizationModelI(name)
 {
 
 }

HadronizationModelBase::HadronizationModelBase	(	string	name,
		string	config
	)

protected

Definition at line 44 of file HadronizationModelBase.cxx.

                                                                          :
 HadronizationModelI(name, config)
 {
 
 }

HadronizationModelBase::~HadronizationModelBase ( )

protectedvirtual

Definition at line 50 of file HadronizationModelBase.cxx.

 {
 
 }

Member Function Documentation

void HadronizationModelBase::ApplyRijk	(	const Interaction *	i,
		bool	norm,
		TH1D *	mp
	)		const

protected

Definition at line 80 of file HadronizationModelBase.cxx.

 {
 // Apply the NEUGEN multiplicity probability scaling factors
 //
   if(!mp) return;
 
   const InitialState & init_state = interaction->InitState();
   int probe_pdg = init_state.ProbePdg();
   int nuc_pdg   = init_state.Tgt().HitNucPdg();
 
   const ProcessInfo & proc_info = interaction->ProcInfo();
   bool is_CC = proc_info.IsWeakCC();
   bool is_NC = proc_info.IsWeakNC();
   bool is_EM = proc_info.IsEM();
   // EDIT
   bool is_dm = proc_info.IsDarkMatter();
 
   //
   // get the R2, R3 factors
   //
 
   double R2=1., R3=1.;
 
   // weak CC or NC case
   // EDIT
   if(is_CC || is_NC || is_dm) {
      bool is_nu    = pdg::IsNeutrino     (probe_pdg); 
      bool is_nubar = pdg::IsAntiNeutrino (probe_pdg);
      bool is_p     = pdg::IsProton       (nuc_pdg);
      bool is_n     = pdg::IsNeutron      (nuc_pdg);
      bool is_dmi   = pdg::IsDarkMatter   (probe_pdg);  // EDIT
      if((is_nu && is_p) || (is_dmi && is_p))  {
          R2 = (is_CC) ? fRvpCCm2 : fRvpNCm2;
          R3 = (is_CC) ? fRvpCCm3 : fRvpNCm3;
      } else 
       if((is_nu && is_n) || (is_dmi && is_n)) {
          R2 = (is_CC) ? fRvnCCm2 : fRvnNCm2;
          R3 = (is_CC) ? fRvnCCm3 : fRvnNCm3;
       } else 
       if(is_nubar && is_p)  {
          R2 = (is_CC) ? fRvbpCCm2 :   fRvbpNCm2;
          R3 = (is_CC) ? fRvbpCCm3 :   fRvbpNCm3;
       } else 
       if(is_nubar && is_n) {
          R2 = (is_CC) ? fRvbnCCm2 : fRvbnNCm2;
          R3 = (is_CC) ? fRvbnCCm3 : fRvbnNCm3;
       } else {
          LOG("BaseHad", pERROR) 
             << "Invalid initial state: " << init_state;
      }
   }//cc||nc?
 
   // EM case (apply the NC tuning factors)
 
   if(is_EM) {
      bool is_l     = pdg::IsNegChargedLepton (probe_pdg); 
      bool is_lbar  = pdg::IsPosChargedLepton (probe_pdg);
      bool is_p     = pdg::IsProton           (nuc_pdg);
      bool is_n     = pdg::IsNeutron          (nuc_pdg);
      if(is_l && is_p)  {
          R2 = fRvpNCm2;
          R3 = fRvpNCm3;
       } else 
       if(is_l && is_n) {
          R2 = fRvnNCm2;
          R3 = fRvnNCm3;
       } else 
       if(is_lbar && is_p)  {
          R2 = fRvbpNCm2;
          R3 = fRvbpNCm3;
       } else 
       if(is_lbar && is_n) {
          R2 = fRvbnNCm2;
          R3 = fRvbnNCm3;
       } else {
          LOG("BaseHad", pERROR) 
             << "Invalid initial state: " << init_state;
      }
   }//em?
 
   //
   // Apply to the multiplicity probability distribution
   //
 
   int nbins = mp->GetNbinsX();
   for(int i = 1; i <= nbins; i++) {
      int n = TMath::Nint( mp->GetBinCenter(i) ); 
 
      double R=1;
      if      (n==2) R=R2;
      else if (n==3) R=R3;
 
      if(n==2 || n==3) {
         double P   = mp->GetBinContent(i);
         double Psc = R*P;
         LOG("BaseHad", pDEBUG) 
           << "n=" << n << "/ Scaling factor R = " 
                               << R << "/ P " << P << " --> " << Psc;
         mp->SetBinContent(i, Psc);
      }
      if(n>3) break;
   }
 
   // renormalize the histogram?
   if(norm) {
      double histo_norm = mp->Integral("width");
      if(histo_norm>0) mp->Scale(1.0/histo_norm);
   }
 }

TH1D * HadronizationModelBase::CreateMultProbHist ( double maxmult ) const

protected

Definition at line 68 of file HadronizationModelBase.cxx.

 {
   double minmult = 2;
   int    nbins   = TMath::Nint(maxmult-minmult+1);
 
   TH1D * mult_prob = new TH1D("mult_prob", 
       "hadronic multiplicity distribution", nbins, minmult-0.5, maxmult+0.5);
   mult_prob->SetDirectory(0);
 
   return mult_prob;
 }

virtual TClonesArray* genie::HadronizationModelBase::Hadronize ( const Interaction * ) const

pure virtual

Implements genie::HadronizationModelI.

Implemented in genie::KNOHadronization, and genie::PythiaHadronization.

virtual void genie::HadronizationModelBase::Initialize ( void ) const

pure virtual

Don't implement the HadronizationModelI interface Leave it for the concrete implementations (KNO, Pythia,...)

Implements genie::HadronizationModelI.

Implemented in genie::KNOHadronization, and genie::PythiaHadronization.

double HadronizationModelBase::MaxMult ( const Interaction * i ) const

protected

Definition at line 60 of file HadronizationModelBase.cxx.

 {
   double W = interaction->Kine().W();
 
   double maxmult = TMath::Floor(1 + (W-kNeutronMass)/kPionMass);
   return maxmult;
 }

virtual TH1D* genie::HadronizationModelBase::MultiplicityProb	(	const Interaction *	,
		Option_t *	opt = `""`
	)		const

pure virtual

Implements genie::HadronizationModelI.

Implemented in genie::KNOHadronization, and genie::PythiaHadronization.

virtual PDGCodeList* genie::HadronizationModelBase::SelectParticles ( const Interaction * ) const

pure virtual

Implements genie::HadronizationModelI.

Implemented in genie::KNOHadronization, and genie::PythiaHadronization.

virtual double genie::HadronizationModelBase::Weight ( void ) const

pure virtual

Implements genie::HadronizationModelI.

Implemented in genie::KNOHadronization, and genie::PythiaHadronization.

double HadronizationModelBase::Wmin ( void ) const

protected

Various utility methods common to hadronization models.

Definition at line 55 of file HadronizationModelBase.cxx.

 {
   return (kNucleonMass+kPionMass);
 }

Member Data Documentation

double genie::HadronizationModelBase::fRvbnCCm2

protected

neugen's Rijk: vbn, CC, multiplicity = 2

Definition at line 68 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvbnCCm3

protected

neugen's Rijk: vbn, CC, multiplicity = 3

Definition at line 69 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvbnNCm2

protected

neugen's Rijk: vbn, NC, multiplicity = 2

Definition at line 70 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvbnNCm3

protected

neugen's Rijk: vbn, NC, multiplicity = 3

Definition at line 71 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvbpCCm2

protected

neugen's Rijk: vbp, CC, multiplicity = 2

Definition at line 64 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvbpCCm3

protected

neugen's Rijk: vbp, CC, multiplicity = 3

Definition at line 65 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvbpNCm2

protected

neugen's Rijk: vbp, NC, multiplicity = 2

Definition at line 66 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvbpNCm3

protected

neugen's Rijk: vbp, NC, multiplicity = 3

Definition at line 67 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvnCCm2

protected

neugen's Rijk: vn, CC, multiplicity = 2

Definition at line 60 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvnCCm3

protected

neugen's Rijk: vn, CC, multiplicity = 3

Definition at line 61 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvnNCm2

protected

neugen's Rijk: vn, NC, multiplicity = 2

Definition at line 62 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvnNCm3

protected

neugen's Rijk: vn, NC, multiplicity = 3

Definition at line 63 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvpCCm2

protected

neugen's Rijk: vp, CC, multiplicity = 2

Definition at line 56 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvpCCm3

protected

neugen's Rijk: vp, CC, multiplicity = 3

Definition at line 57 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvpNCm2

protected

neugen's Rijk: vp, NC, multiplicity = 2

Definition at line 58 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fRvpNCm3

protected

neugen's Rijk: vp, NC, multiplicity = 3

Definition at line 59 of file HadronizationModelBase.h.

double genie::HadronizationModelBase::fWcut

protected

configuration data common to all hadronizers

neugen's Rijk applied for W<Wcut (see DIS/RES join scheme)

Definition at line 55 of file HadronizationModelBase.h.

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

Generator/src/Physics/Hadronization/HadronizationModelBase.h
Generator/src/Physics/Hadronization/HadronizationModelBase.cxx

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation