GENIEGenerator
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Private Member Functions | Private Attributes | List of all members
genie::ReinSehgalSPPPXSec Class Reference

Computes the differential cross section for an exclusive 1-pion reaction through resonance neutrinoproduction according to the Rein-Sehgal model. More...

#include <ReinSehgalSPPPXSec.h>

Inheritance diagram for genie::ReinSehgalSPPPXSec:
Inheritance graph
[legend]
Collaboration diagram for genie::ReinSehgalSPPPXSec:
Collaboration graph
[legend]

Public Member Functions

 ReinSehgalSPPPXSec ()
 
 ReinSehgalSPPPXSec (string config)
 
virtual ~ReinSehgalSPPPXSec ()
 
double XSec (const Interaction *i, KinePhaseSpace_t k) const
 Compute the cross section for the input interaction. More...
 
double Integral (const Interaction *i) const
 
bool ValidProcess (const Interaction *i) const
 Can this cross section algorithm handle the input process? More...
 
void Configure (const Registry &config)
 
void Configure (string config)
 
- Public Member Functions inherited from genie::XSecAlgorithmI
virtual ~XSecAlgorithmI ()
 
virtual bool ValidKinematics (const Interaction *i) const
 Is the input kinematical point a physically allowed one? More...
 
- Public Member Functions inherited from genie::Algorithm
virtual ~Algorithm ()
 
virtual void FindConfig (void)
 
virtual const RegistryGetConfig (void) const
 
RegistryGetOwnedConfig (void)
 
virtual const AlgIdId (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 AlgorithmSubAlg (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)
 
double XSecNRES (const Interaction *i, KinePhaseSpace_t k) const
 
double XSec1RES (const Interaction *i, KinePhaseSpace_t k) const
 

Private Attributes

BaryonResList fResList
 
const XSecAlgorithmIfSingleResXSecModel
 
const XSecIntegratorIfXSecIntegrator
 

Additional Inherited Members

- Static Public Member Functions inherited from genie::Algorithm
static string BuildParamVectKey (const std::string &comm_name, unsigned int i)
 
static string BuildParamVectSizeKey (const std::string &comm_name)
 
static string BuildParamMatKey (const std::string &comm_name, unsigned int i, unsigned int j)
 
static string BuildParamMatRowSizeKey (const std::string &comm_name)
 
static string BuildParamMatColSizeKey (const std::string &comm_name)
 
- Protected Member Functions inherited from genie::XSecAlgorithmI
 XSecAlgorithmI ()
 
 XSecAlgorithmI (string name)
 
 XSecAlgorithmI (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)
 
RegistryExtractLocalConfig (const Registry &in) const
 
RegistryExtractLowerConfig (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 >
int GetParamVect (const std::string &comm_name, std::vector< T > &v, bool is_top_call=true) const
 Handle to load vectors of parameters. More...
 
int GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
 
template<class T >
int GetParamMat (const std::string &comm_name, TMatrixT< T > &mat, bool is_top_call=true) const
 Handle to load matrix of parameters. More...
 
template<class T >
int GetParamMatSym (const std::string &comm_name, TMatrixTSym< T > &mat, bool is_top_call=true) const
 
int GetParamMatKeys (const std::string &comm_name, std::vector< RgKey > &k, 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 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...
 
AlgMapfOwnedSubAlgMp
 local pool for owned sub-algs (taken out of the factory pool) More...
 

Detailed Description

Computes the differential cross section for an exclusive 1-pion reaction through resonance neutrinoproduction according to the Rein-Sehgal model.

      The cross section is computed for an input list of resonances
      as the sum of the Rein-Sehgal single resonance cross sections
      weighted:

      \li With the value of their Breit-Wigner distributions at the given
          W,Q^2 (The code for BW weighting is included in the single
          resonance cross section algorithm. The user needs to make sure
          that he does not run the single resonance cross section code with
          a configuration that inhibits weighting).

      \li With the isospin Glebsch-Gordon coefficient determining the
          contribution of each resonance to the exclusive final state.

      \li With the BR for the produced resonance to decay into the given
          exclusive final state.

      In this algorithm we follow the non-coherent approach: we sum
      the weighted resonance production cross sections rather than the
      resonance production amplitudes.

      Is a concrete implementation of the XSecAlgorithmI interface.
References:
D.Rein and L.M.Sehgal, Neutrino Excitation of Baryon Resonances and Single Pion Production, Ann.Phys.133, 79 (1981)
Author
Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool
Created:
November 22, 2004
License:
Copyright (c) 2003-2024, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 55 of file ReinSehgalSPPPXSec.h.

Constructor & Destructor Documentation

ReinSehgalSPPPXSec::ReinSehgalSPPPXSec ( )

Definition at line 27 of file ReinSehgalSPPPXSec.cxx.

27  :
28 XSecAlgorithmI("genie::ReinSehgalSPPPXSec")
29 {
30 
31 }
ReinSehgalSPPPXSec::ReinSehgalSPPPXSec ( string  config)

Definition at line 33 of file ReinSehgalSPPPXSec.cxx.

33  :
34 XSecAlgorithmI("genie::ReinSehgalSPPPXSec", config)
35 {
36 
37 }
ReinSehgalSPPPXSec::~ReinSehgalSPPPXSec ( )
virtual

Definition at line 39 of file ReinSehgalSPPPXSec.cxx.

40 {
41 
42 }

Member Function Documentation

void ReinSehgalSPPPXSec::Configure ( const Registry config)
virtual

Configure the algorithm with an external registry The registry is merged with the top level registry if it is owned, Otherwise a copy of it is added with the highest priority

Reimplemented from genie::Algorithm.

Definition at line 160 of file ReinSehgalSPPPXSec.cxx.

References genie::Algorithm::Configure(), and LoadConfig().

161 {
162  Algorithm::Configure(config);
163  this->LoadConfig();
164 }
virtual void Configure(const Registry &config)
Definition: Algorithm.cxx:62
void ReinSehgalSPPPXSec::Configure ( string  config)
virtual

Configure the algorithm from the AlgoConfigPool based on param_set string given in input An algorithm contains a vector of registries coming from different xml configuration files, which are loaded according a very precise prioriy This methods will load a number registries in order of priority: 1) "Tunable" parameter set from CommonParametes. This is loaded with the highest prioriry and it is designed to be used for tuning procedure Usage not expected from the user. 2) For every string defined in "CommonParame" the corresponding parameter set will be loaded from CommonParameter.xml 3) parameter set specified by the config string and defined in the xml file of the algorithm 4) if config is not "Default" also the Default parameter set from the same xml file will be loaded Effectively this avoids the repetion of a parameter when it is not changed in the requested configuration

Reimplemented from genie::Algorithm.

Definition at line 166 of file ReinSehgalSPPPXSec.cxx.

References genie::Algorithm::Configure(), and LoadConfig().

167 {
168  Algorithm::Configure(config);
169  this->LoadConfig();
170 }
virtual void Configure(const Registry &config)
Definition: Algorithm.cxx:62
double ReinSehgalSPPPXSec::Integral ( const Interaction i) const
virtual

Integrate the model over the kinematic phase space available to the input interaction (kinematical cuts can be included)

Implements genie::XSecAlgorithmI.

Definition at line 137 of file ReinSehgalSPPPXSec.cxx.

References fXSecIntegrator, and genie::XSecIntegratorI::Integrate().

138 {
139  return fXSecIntegrator->Integrate(this,interaction);
140 }
const XSecIntegratorI * fXSecIntegrator
virtual double Integrate(const XSecAlgorithmI *model, const Interaction *interaction) const =0
void ReinSehgalSPPPXSec::LoadConfig ( void  )
private

Definition at line 172 of file ReinSehgalSPPPXSec.cxx.

References genie::BaryonResList::Clear(), genie::BaryonResList::DecodeFromNameList(), fResList, fSingleResXSecModel, fXSecIntegrator, genie::Algorithm::GetParam(), and genie::Algorithm::SubAlg().

Referenced by Configure().

173 {
174 // load the single resonance cross section algorithm specified in the config.
175 
177  dynamic_cast<const XSecAlgorithmI *> (this->SubAlg("SingleRESDiffXSecAlg"));
178  assert(fSingleResXSecModel);
179 
180  //-- Create a BaryonResList by decoding the resonance list from
181  // the XML input
182  // The list of resonances can be specified as a string with
183  // comma separated resonance names (eg "P33(1233),S11(1535),D13(1520)")
184  // The BaryonResList can also decode lists of pdg-codes or
185  // resonance-ids (Resonance_t enumerations).
186  // Support for this will be added here as well.
187 
188  fResList.Clear();
189 
190 
191  string resonances ;
192  GetParam( "ResonanceNameList", resonances ) ;
193  fResList.DecodeFromNameList(resonances);
194 
195  //-- load the differential cross section integrator
197  dynamic_cast<const XSecIntegratorI *> (this->SubAlg("XSec-Integrator"));
198  assert(fXSecIntegrator);
199 }
Cross Section Calculation Interface.
Cross Section Integrator Interface.
void DecodeFromNameList(string list, string delimiter=",")
const XSecIntegratorI * fXSecIntegrator
bool GetParam(const RgKey &name, T &p, bool is_top_call=true) const
const XSecAlgorithmI * fSingleResXSecModel
const Algorithm * SubAlg(const RgKey &registry_key) const
Definition: Algorithm.cxx:345
bool ReinSehgalSPPPXSec::ValidProcess ( const Interaction i) const
virtual

Can this cross section algorithm handle the input process?

Implements genie::XSecAlgorithmI.

Definition at line 142 of file ReinSehgalSPPPXSec.cxx.

References genie::SppChannel::AsString(), genie::SppChannel::FromInteraction(), genie::kISkipProcessChk, genie::kSppNull, LOG, pDEBUG, and pERROR.

Referenced by XSec().

143 {
144  if(interaction->TestBit(kISkipProcessChk)) return true;
145 
146  //-- Get the requested SPP channel
147  SppChannel_t spp_channel = SppChannel::FromInteraction(interaction);
148  if( spp_channel == kSppNull ) {
149  LOG("ReinSehgalSpp", pERROR)
150  << "\n *** Insufficient SPP exclusive final state information!";
151  return false;
152  }
153 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
154  LOG("ReinSehgalSpp", pDEBUG)
155  << "Reaction: " << SppChannel::AsString(spp_channel);
156 #endif
157  return true;
158 }
static SppChannel_t FromInteraction(const Interaction *interaction)
Definition: SppChannel.h:402
#define pERROR
Definition: Messenger.h:59
enum genie::ESppChannel SppChannel_t
static string AsString(SppChannel_t channel)
Definition: SppChannel.h:76
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
const UInt_t kISkipProcessChk
if set, skip process validity checks
Definition: Interaction.h:47
#define pDEBUG
Definition: Messenger.h:63
double ReinSehgalSPPPXSec::XSec ( const Interaction i,
KinePhaseSpace_t  k 
) const
virtual

Compute the cross section for the input interaction.

Implements genie::XSecAlgorithmI.

Definition at line 44 of file ReinSehgalSPPPXSec.cxx.

References genie::utils::res::AsString(), genie::Interaction::ExclTag(), genie::BaryonResList::Find(), fResList, fSingleResXSecModel, genie::kNoResonance, LOG, pDEBUG, pWARN, genie::XclsTag::Resonance(), genie::XSecAlgorithmI::ValidKinematics(), ValidProcess(), genie::XSecAlgorithmI::XSec(), and XSecNRES().

46 {
47  if(! this -> ValidProcess (interaction) ) return 0.;
48  if(! this -> ValidKinematics (interaction) ) return 0.;
49 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
50  LOG("ReinSehgalSpp", pDEBUG)
51  << "Computing a cross section for " << *interaction;
52 #endif
53  //-- Check whether a resonance has been specified
54  // If yes, compute only the contribution of this resonance at the
55  // specified exclusive state
56 
57  Resonance_t inpres = interaction->ExclTag().Resonance();
58  if(inpres != kNoResonance) {
59 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
60  LOG("ReinSehgalSpp", pDEBUG)
61  << "Computing only the contribution from: " << utils::res::AsString(inpres);
62 #endif
63  if(!fResList.Find(inpres)) {
64  LOG("ReinSehgalSpp", pWARN)
65  << "Resonance: " << utils::res::AsString(inpres) << " was not found in my list";
66  return 0;
67  }
68  //-- Compute the contribution of this resonance
69  //-- Get the Breit-Wigner weighted xsec for exciting the resonance
70 
71  return fSingleResXSecModel->XSec(interaction,kps);
72  }
73 
74  //-- Loop over the specified list of baryon resonances and compute
75  // the cross section for the input exclusive channel
76 
77  return this->XSecNRES(interaction,kps);
78 }
enum genie::EResonance Resonance_t
bool ValidProcess(const Interaction *i) const
Can this cross section algorithm handle the input process?
virtual bool ValidKinematics(const Interaction *i) const
Is the input kinematical point a physically allowed one?
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
#define pWARN
Definition: Messenger.h:60
virtual double XSec(const Interaction *i, KinePhaseSpace_t k=kPSfE) const =0
Compute the cross section for the input interaction.
bool Find(Resonance_t res) const
double XSecNRES(const Interaction *i, KinePhaseSpace_t k) const
const char * AsString(Resonance_t res)
resonance id -&gt; string
const XSecAlgorithmI * fSingleResXSecModel
#define pDEBUG
Definition: Messenger.h:63
double genie::ReinSehgalSPPPXSec::XSec1RES ( const Interaction i,
KinePhaseSpace_t  k 
) const
private
double ReinSehgalSPPPXSec::XSecNRES ( const Interaction i,
KinePhaseSpace_t  k 
) const
private

Definition at line 80 of file ReinSehgalSPPPXSec.cxx.

References genie::utils::res::AsString(), genie::SppChannel::AsString(), genie::SppChannel::BranchingRatio(), genie::Interaction::ExclTagPtr(), fResList, genie::SppChannel::FromInteraction(), fSingleResXSecModel, genie::SppChannel::IsospinWeight(), genie::kNoResonance, LOG, genie::BaryonResList::NResonances(), pDEBUG, genie::BaryonResList::ResonanceId(), genie::XclsTag::SetResonance(), and genie::XSecAlgorithmI::XSec().

Referenced by XSec().

82 {
83 // computes the 1pi cros section taking into account the contribution of all
84 // specified baryon resonances
85 
86  unsigned int nres = fResList.NResonances();
87 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
88  LOG("ReinSehgalSpp", pDEBUG)
89  << "Computing SPP cross section using " << nres << " resonances";
90 #endif
91 
92  //-- Get 1pi exclusive channel
93  SppChannel_t spp_channel = SppChannel::FromInteraction(interaction);
94 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
95  LOG("ReinSehgalSpp", pDEBUG)
96  << "SPP channel " << SppChannel::AsString(spp_channel);
97 #endif
98 
99  double xsec = 0;
100  for(unsigned int ires = 0; ires < nres; ires++) {
101 
102  //-- Get next resonance from the resonance list
103  Resonance_t res = fResList.ResonanceId(ires);
104 
105  //-- Set current resonance to interaction object
106  interaction->ExclTagPtr()->SetResonance(res);
107 
108  //-- Get the BR for the (resonance) -> (exclusive final state)
109  double br = SppChannel::BranchingRatio(/*spp_channel,*/ res);
110 
111  //-- Get the Isospin Clebsch-Gordon coefficient for the given resonance
112  // and exclusive final state
113  double igg = SppChannel::IsospinWeight(spp_channel, res);
114 
115  //-- Compute the weighted xsec
116  // (total weight = Breit-Wigner * BR * isospin Clebsch-Gordon)
117  double res_xsec_contrib = fSingleResXSecModel->XSec(interaction,kps)*br*igg;
118 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
119  LOG("ReinSehgalSpp", pDEBUG)
120  << "Contrib. from [" << utils::res::AsString(res) << "] = "
121  << "<Clebsch-Gordon = " << igg
122  << "> * <BR(->1pi) = " << br
123  << "> * <Breit-Wigner * d^nxsec/dK^n = " << rxsec
124  << "> = " << res_xsec_contrib;
125 #endif
126 
127  //-- Add contribution of this resonance to the cross section
128  xsec += res_xsec_contrib;
129  }
130 
131  //-- delete the resonance from the input interaction
132  interaction->ExclTagPtr()->SetResonance(kNoResonance);
133 
134  return xsec;
135 }
static SppChannel_t FromInteraction(const Interaction *interaction)
Definition: SppChannel.h:402
static double IsospinWeight(SppChannel_t channel, Resonance_t res)
Definition: SppChannel.h:238
enum genie::EResonance Resonance_t
unsigned int NResonances(void) const
enum genie::ESppChannel SppChannel_t
static string AsString(SppChannel_t channel)
Definition: SppChannel.h:76
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
virtual double XSec(const Interaction *i, KinePhaseSpace_t k=kPSfE) const =0
Compute the cross section for the input interaction.
const char * AsString(Resonance_t res)
resonance id -&gt; string
static double BranchingRatio(Resonance_t res)
Definition: SppChannel.h:357
const XSecAlgorithmI * fSingleResXSecModel
#define pDEBUG
Definition: Messenger.h:63
Resonance_t ResonanceId(unsigned int ires) const

Member Data Documentation

BaryonResList genie::ReinSehgalSPPPXSec::fResList
private

Definition at line 81 of file ReinSehgalSPPPXSec.h.

Referenced by LoadConfig(), XSec(), and XSecNRES().

const XSecAlgorithmI* genie::ReinSehgalSPPPXSec::fSingleResXSecModel
private

Definition at line 82 of file ReinSehgalSPPPXSec.h.

Referenced by LoadConfig(), XSec(), and XSecNRES().

const XSecIntegratorI* genie::ReinSehgalSPPPXSec::fXSecIntegrator
private

Definition at line 83 of file ReinSehgalSPPPXSec.h.

Referenced by Integral(), and LoadConfig().


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