GENIEGenerator
Loading...
Searching...
No Matches
Public Member Functions | Private Member Functions | Private Attributes | List of all members
genie::GLRESPXSec Class Reference

Differential cross section for glashow resonance. More...

#include <GLRESPXSec.h>

Inheritance diagram for genie::GLRESPXSec:
[legend]
Collaboration diagram for genie::GLRESPXSec:
[legend]

Public Member Functions

 GLRESPXSec ()
 GLRESPXSec (string config)
virtual ~GLRESPXSec ()
double XSec (const Interaction *i, KinePhaseSpace_t k) const
 Compute the cross section for the input interaction.
double Integral (const Interaction *i) const
bool ValidProcess (const Interaction *i) const
 Can this cross section algorithm handle the input process?
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?
Public Member Functions inherited from genie::Algorithm
virtual ~Algorithm ()
virtual void FindConfig (void)
virtual const Registry & GetConfig (void) const
Registry * GetOwnedConfig (void)
virtual const AlgId & Id (void) const
 Get algorithm ID.
virtual AlgStatus_t GetStatus (void) const
 Get algorithm status.
virtual bool AllowReconfig (void) const
virtual AlgCmp_t Compare (const Algorithm *alg) const
 Compare with input algorithm.
virtual void SetId (const AlgId &id)
 Set algorithm ID.
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.

Private Member Functions

void LoadConfig (void)

Private Attributes

const XSecIntegratorI * fXSecIntegrator
 diff. xsec integrator
double fWmin
 Minimum value of W.
Born * born

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)
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.
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.
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.
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
int AddLowRegistry (Registry *rp, bool owns=true)
 add registry with lowest priority, also update ownership
int MergeTopRegistry (const Registry &r)
int AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
 Add registries with top priority, also udated Ownerships.
Protected Attributes inherited from genie::Algorithm
bool fAllowReconfig
bool fOwnsSubstruc
 true if it owns its substructure (sub-algs,...)
AlgId fID
 algorithm name and configuration set
vector< Registry * > fConfVect
vector< bool > fOwnerships
 ownership for every registry in fConfVect
AlgStatus_t fStatus
 algorithm execution status
AlgMap * fOwnedSubAlgMp
 local pool for owned sub-algs (taken out of the factory pool)

Detailed Description

Differential cross section for glashow resonance.

Author
Alfonso Garcia <aagarciasoto \at km3net.de> IFIC & Harvard University
References:\n Phys. Rev. D 100, 091301 (2019)
Created:\n Dec 8, 2021
License:\n Copyright (c) 2003-2025, The GENIE Collaboration
For the full text of the license visit http://copyright.genie-mc.org or see $GENIE/LICENSE

Definition at line 31 of file GLRESPXSec.h.

Constructor & Destructor Documentation

◆ GLRESPXSec() [1/2]

GLRESPXSec::GLRESPXSec ( )

Definition at line 27 of file GLRESPXSec.cxx.

27 :
28XSecAlgorithmI("genie::GLRESPXSec")
29{
30 born = new Born();
31
32}

References born, and genie::XSecAlgorithmI::XSecAlgorithmI().

◆ GLRESPXSec() [2/2]

GLRESPXSec::GLRESPXSec ( string config)

Definition at line 34 of file GLRESPXSec.cxx.

34 :
35XSecAlgorithmI("genie::GLRESPXSec", config)
36{
37
38}

References genie::XSecAlgorithmI::XSecAlgorithmI().

◆ ~GLRESPXSec()

GLRESPXSec::~GLRESPXSec ( )
virtual

Definition at line 40 of file GLRESPXSec.cxx.

41{
42
43}

Member Function Documentation

◆ Configure() [1/2]

void GLRESPXSec::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 137 of file GLRESPXSec.cxx.

138{
139
140 Algorithm::Configure(config);
141 this->LoadConfig();
142}
genie::Algorithm::Configure
virtual void Configure(const Registry &config)
Definition Algorithm.cxx:62
genie::GLRESPXSec::LoadConfig
void LoadConfig(void)
Definition GLRESPXSec.cxx:150

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

◆ Configure() [2/2]

void GLRESPXSec::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 144 of file GLRESPXSec.cxx.

145{
146 Algorithm::Configure(config);
147 this->LoadConfig();
148}

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

◆ Integral()

double GLRESPXSec::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 114 of file GLRESPXSec.cxx.

115{
116 double xsec = fXSecIntegrator->Integrate(this,interaction);
117
118 return xsec;
119}
genie::GLRESPXSec::fXSecIntegrator
const XSecIntegratorI * fXSecIntegrator
diff. xsec integrator
Definition GLRESPXSec.h:51

References fXSecIntegrator.

◆ LoadConfig()

void GLRESPXSec::LoadConfig ( void )
private

Definition at line 150 of file GLRESPXSec.cxx.

151{
152
153 //-- load the differential cross section integrator
154 fXSecIntegrator = dynamic_cast<const XSecIntegratorI *> (this->SubAlg("XSec-Integrator"));
155 assert(fXSecIntegrator);
156
157 GetParam( "Xsec-Wmin", fWmin ) ;
158
159
160}
genie::Algorithm::GetParam
bool GetParam(const RgKey &name, T &p, bool is_top_call=true) const
genie::Algorithm::SubAlg
const Algorithm * SubAlg(const RgKey &registry_key) const
Definition Algorithm.cxx:345
genie::GLRESPXSec::fWmin
double fWmin
Minimum value of W.
Definition GLRESPXSec.h:53

References fWmin, fXSecIntegrator, genie::Algorithm::GetParam(), and genie::Algorithm::SubAlg().

Referenced by Configure(), and Configure().

◆ ValidProcess()

bool GLRESPXSec::ValidProcess ( const Interaction * i) const
virtual

Can this cross section algorithm handle the input process?

Implements genie::XSecAlgorithmI.

Definition at line 121 of file GLRESPXSec.cxx.

122{
123 if(interaction->TestBit(kISkipProcessChk)) return true;
124
125 const ProcessInfo & proc_info = interaction->ProcInfo();
126 if(!proc_info.IsGlashowResonance()) return false;
127 if(!proc_info.IsWeakCC()) return false;
128
129 const InitialState & init_state = interaction -> InitState();
130 if(!pdg::IsAntiNuE(init_state.ProbePdg())) return false;
131
132 if(init_state.Tgt().HitNucIsSet()) return false;
133
134 return true;
135}
genie::InitialState::Tgt
const Target & Tgt(void) const
Definition InitialState.h:66
genie::InitialState::ProbePdg
int ProbePdg(void) const
Definition InitialState.h:64
genie::ProcessInfo::IsWeakCC
bool IsWeakCC(void) const
Definition ProcessInfo.cxx:203
genie::ProcessInfo::IsGlashowResonance
bool IsGlashowResonance(void) const
Definition ProcessInfo.cxx:156
genie::Target::HitNucIsSet
bool HitNucIsSet(void) const
Definition Target.cxx:283
genie::pdg::IsAntiNuE
bool IsAntiNuE(int pdgc)
Definition PDGUtils.cxx:173
genie::kISkipProcessChk
const UInt_t kISkipProcessChk
if set, skip process validity checks
Definition Interaction.h:47

References genie::Target::HitNucIsSet(), genie::pdg::IsAntiNuE(), genie::ProcessInfo::IsGlashowResonance(), genie::ProcessInfo::IsWeakCC(), genie::kISkipProcessChk, genie::InitialState::ProbePdg(), genie::Interaction::ProcInfo(), and genie::InitialState::Tgt().

Referenced by XSec().

◆ XSec()

double GLRESPXSec::XSec ( const Interaction * i,
KinePhaseSpace_t k ) const
virtual

Compute the cross section for the input interaction.

Implements genie::XSecAlgorithmI.

Definition at line 45 of file GLRESPXSec.cxx.

47{
48
49 if(! this -> ValidProcess (interaction) ) return 0.;
50
51 const InitialState & init_state = interaction -> InitState();
52 const Kinematics & kinematics = interaction -> Kine();
53 const XclsTag & xclstag = interaction -> ExclTag();
54
55 int loutpdg = xclstag.FinalLeptonPdg();
56
57 double mlout = interaction->FSPrimLepton()->Mass(); //mass of charged lepton
58 double mlin = kElectronMass; //mass of incoming charged lepton
59
60 double Enuin = init_state.ProbeE(kRfLab);
61 double s = born->GetS(mlin,Enuin);
62
63 double n1 = kinematics.GetKV(kKVn1);
64 double n2 = kinematics.GetKV(kKVn2);
65 double t = born->GetT(mlin,mlout,s,n1);
66 if (t>0) return 0.;
67
68 //nlo correction
69 double zeta = born->GetReAlpha()/kPi*(2.*TMath::Log(TMath::Sqrt(-t)/kElectronMass)-1.);
70 double omx = TMath::Power(n2, 1./zeta );
71 double pdf_soft = TMath::Exp(zeta*(3./4.-TMath::EulerGamma()))/TMath::Gamma(1.+zeta) + omx*(omx-2.)/2./n2;
72 if ( omx<0. || omx>1. ) return 0.;
73 double s_r = s*(1. - omx);
74 double t_r = t*(1. - omx);
75
76 //http://users.jyu.fi/~tulappi/fysh300sl11/l2.pdf [slide 22]
77 //remember we always define nuout as p4
78 double Enuout = (mlin*mlin-t_r)/2./mlin;
79 if ( !born->IsInPhaseSpace(mlin,mlout,Enuin,Enuout) ) return 0.;
80
81 double xsec = kPi/4./(s-mlin*mlin) * pdf_soft ;
82
83 if ( pdg::IsPion(loutpdg) ) {
84 if ( TMath::Sqrt(s_r)<fWmin ) return 0.; // The W limit is because hadronization might have issues at low W (as in PYTHIA6).
85 xsec *= 64.41/10.63;
86 }
87
88 double ME = 0.;
89 if ( loutpdg == kPdgElectron ) ME = born->PXSecCCRNC(s_r,t_r,mlin,mlout);
90 else ME = born->PXSecCCR (s_r,t_r,mlin,mlout);
91 xsec *= TMath::Max(0.,ME);
92
93 //----- If requested return the free electron xsec even for nuclear target
94 if( interaction->TestBit(kIAssumeFreeElectron) ) return xsec;
95
96 //----- Scale for the number of scattering centers at the target
97 int Ne = init_state.Tgt().Z(); // num of scattering centers
98 xsec *= Ne;
99
100 if(kps!=kPSn1n2fE) {
101 LOG("GLRESPXSec", pWARN)
102 << "Doesn't support transformation from "
103 << KinePhaseSpace::AsString(kPSn1n2fE) << " to "
104 << KinePhaseSpace::AsString(kps);
105 xsec = 0;
106 }
107
108 LOG("GLRESPXSec", pINFO) << "dxsec/dn1dn2 (E= " << Enuin << ", n1= " << n1 << ", n2=" << n2 << ") = " << xsec;
109
110 return xsec;
111
112}
pINFO
#define pINFO
Definition Messenger.h:62
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE,...
Definition Messenger.h:96
pWARN
#define pWARN
Definition Messenger.h:60
genie::GLRESPXSec::ValidProcess
bool ValidProcess(const Interaction *i) const
Can this cross section algorithm handle the input process?
Definition GLRESPXSec.cxx:121
genie::InitialState::ProbeE
double ProbeE(RefFrame_t rf) const
Definition InitialState.cxx:384
genie::KinePhaseSpace::AsString
static string AsString(KinePhaseSpace_t kps)
Definition KinePhaseSpace.h:85
genie::Kinematics::GetKV
double GetKV(KineVar_t kv) const
Definition Kinematics.cxx:323
genie::Target::Z
int Z(void) const
Definition Target.h:68
genie::XclsTag::FinalLeptonPdg
int FinalLeptonPdg(void) const
Definition XclsTag.h:74
genie::pdg::IsPion
bool IsPion(int pdgc)
Definition PDGUtils.cxx:326
genie::kIAssumeFreeElectron
const UInt_t kIAssumeFreeElectron
Definition Interaction.h:50
genie::kKVn2
@ kKVn2
Definition KineVar.h:62
genie::kKVn1
@ kKVn1
Definition KineVar.h:61
genie::kRfLab
@ kRfLab
Definition RefFrame.h:26
genie::kPdgElectron
const int kPdgElectron
Definition PDGCodes.h:35

References genie::KinePhaseSpace::AsString(), born, genie::XclsTag::FinalLeptonPdg(), genie::Interaction::FSPrimLepton(), fWmin, genie::Kinematics::GetKV(), genie::pdg::IsPion(), genie::constants::kElectronMass, genie::kIAssumeFreeElectron, genie::kKVn1, genie::kKVn2, genie::kPdgElectron, genie::constants::kPi, genie::kPSn1n2fE, genie::kRfLab, LOG, pINFO, genie::InitialState::ProbeE(), pWARN, genie::InitialState::Tgt(), ValidProcess(), and genie::Target::Z().

Member Data Documentation

◆ born

Born* genie::GLRESPXSec::born
private

Definition at line 55 of file GLRESPXSec.h.

Referenced by GLRESPXSec(), and XSec().

◆ fWmin

double genie::GLRESPXSec::fWmin
private

Minimum value of W.

Definition at line 53 of file GLRESPXSec.h.

Referenced by LoadConfig(), and XSec().

◆ fXSecIntegrator

const XSecIntegratorI* genie::GLRESPXSec::fXSecIntegrator
private

diff. xsec integrator

Definition at line 51 of file GLRESPXSec.h.

Referenced by Integral(), and LoadConfig().

The documentation for this class was generated from the following files:
Generated on for GENIEGenerator by doxygen 1.14.0