Classes | Functions
genie::utils::res Namespace Reference

Baryon Resonance utilities. More...

Classes

struct  CacheBWNorm
 cached breit-wigner normalization factor More...
 

Functions

const char * AsString (Resonance_t res)
 resonance id -> string More...
 
Resonance_t FromString (const char *res)
 string -> resonance id More...
 
int PdgCode (Resonance_t res, int Q)
 (resonance id, charge) -> PDG code More...
 
Resonance_t FromPdgCode (int pdgc)
 PDG code -> resonance id. More...
 
bool IsBaryonResonance (int pdgc)
 is input a baryon resonance? More...
 
bool IsDelta (Resonance_t res)
 is it a Delta resonance? More...
 
bool IsN (Resonance_t res)
 is it an N resonance? More...
 
double Mass (Resonance_t res)
 resonance mass (GeV) More...
 
double Width (Resonance_t res)
 resonance width (GeV) More...
 
double BWNorm (Resonance_t res, double N0ResMaxNWidths=6, double N2ResMaxNWidths=2, double GnResMaxNWidths=4)
 breit-wigner normalization factor More...
 
int OrbitalAngularMom (Resonance_t res)
 orbital angular momentum More...
 
int ResonanceIndex (Resonance_t res)
 resonance idx, quark model / SU(6) More...
 

Detailed Description

Baryon Resonance utilities.

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

November 25, 2004

Copyright (c) 2003-2019, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org or see $GENIE/LICENSE

Function Documentation

const char * genie::utils::res::AsString ( Resonance_t  res)

resonance id -> string

Definition at line 49 of file BaryonResUtils.cxx.

References genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kS11_1535, genie::kS11_1650, and genie::kS31_1620.

Referenced by genie::ReinSehgalRESXSecWithCache::CacheBranchName(), genie::ReinSehgalRESXSecWithCacheFast::CacheBranchName(), genie::ReinSehgalRESXSecWithCache::CacheResExcitationXSec(), genie::ReinSehgalRESXSecWithCacheFast::CacheResExcitationXSec(), novarwgt::internal::GetGenieKnobName(), ana::GetGenieKnobSyst(), caf::CAFMaker::InitializeOutfile(), genie::ReinSehgalRESXSec::Integrate(), genie::ReinSehgalRESXSecFast::Integrate(), genie::ReinSehgalSPPXSec::Integrate(), main(), genie::BaryonResList::Print(), genie::XclsTag::Print(), rwgt::MakeGENIEReweightTable::produce(), genie::BaryonResList::ResonanceName(), SaveGraphsToRootFile(), genie::RSPPResonanceSelector::SelectResonance(), genie::ReinSehgalRESPXSec::XSec(), genie::BSKLNBaseRESPXSec2014::XSec(), genie::ReinSehgalSPPPXSec::XSec(), and genie::ReinSehgalSPPPXSec::XSecNRES().

50 {
51  switch(res) {
52  case kP33_1232 : return "P33(1232)" ; break;
53  case kS11_1535 : return "S11(1535)" ; break;
54  case kD13_1520 : return "D13(1520)" ; break;
55  case kS11_1650 : return "S11(1650)" ; break;
56  case kD13_1700 : return "D13(1700)" ; break;
57  case kD15_1675 : return "D15(1675)" ; break;
58  case kS31_1620 : return "S31(1620)" ; break;
59  case kD33_1700 : return "D33(1700)" ; break;
60  case kP11_1440 : return "P11(1440)" ; break;
61  case kP33_1600 : return "P33(1600)" ; break;
62  case kP13_1720 : return "P13(1720)" ; break;
63  case kF15_1680 : return "F15(1680)" ; break;
64  case kP31_1910 : return "P31(1910)" ; break;
65  case kP33_1920 : return "P33(1920)" ; break;
66  case kF35_1905 : return "F35(1905)" ; break;
67  case kF37_1950 : return "F37(1950)" ; break;
68  case kP11_1710 : return "P11(1710)" ; break;
69  case kF17_1970 : return "F17(1970)" ; break;
70  default: break;
71  }
72  return "unknown resonance!";
73 }
double genie::utils::res::BWNorm ( Resonance_t  res,
double  N0ResMaxNWidths = 6,
double  N2ResMaxNWidths = 2,
double  GnResMaxNWidths = 4 
)

breit-wigner normalization factor

Definition at line 487 of file BaryonResUtils.cxx.

References genie::utils::bwfunc::BreitWignerL(), genie::utils::res::CacheBWNorm::cache, MECModelEnuComparisons::i, genie::kNoResonance, Mass(), OrbitalAngularMom(), ResonanceIndex(), sum, W, Width(), Wmax, and Wmin.

Referenced by main(), genie::ReinSehgalRESPXSec::XSec(), and genie::BSKLNBaseRESPXSec2014::XSec().

488 {
489  static genie::utils::res::CacheBWNorm cbwn;
490  if (res==kNoResonance) return -1;
491  if (cbwn.cache[res]!=0) return cbwn.cache[res];
492 
493 
494  // Get baryon resonance parameters
495  int IR = utils::res::ResonanceIndex (res);
496  int LR = utils::res::OrbitalAngularMom (res);
497  double MR = utils::res::Mass (res);
498  double WR = utils::res::Width (res);
499 
500  // imported part of code from src/contrib/misc/bwnorm.C
501 
502  double NW = GnResMaxNWidths;
503  if(IR==2) NW = N2ResMaxNWidths;
504  if(IR==0) NW = N0ResMaxNWidths;
505 
506  double Wmin = 0.001;
507  double Wmax = MR + NW*WR;
508  int N = 1000* TMath::Nint( (Wmax-Wmin)/WR );
509  if(N%2==0) N++;
510 
511  double dW = (Wmax-Wmin)/(N-1);
512 
513  double sum = 0.5 * (genie::utils::bwfunc::BreitWignerL(Wmin,LR,MR,WR,1.0) + genie::utils::bwfunc::BreitWignerL(Wmax,LR,MR,WR,1.0));
514 
515  for(int i=1; i<N-1; i++) {
516  double W = Wmin + i*dW;
517  sum += ( genie::utils::bwfunc::BreitWignerL(W,LR,MR,WR,1.0) * (i%2+1) );
518  }
519  sum *= (2.*dW/3.);
520  cbwn.cache[res]=sum;
521  return cbwn.cache[res];
522 }
genie::utils::res::CacheBWNorm
cached breit-wigner normalization factor
Definition: BaryonResUtils.h:51
Wmax
double Wmax[kNWBins]
Definition: piN_asymmetric_decay.C:40
genie::utils::res::Mass
double Mass(Resonance_t res)
resonance mass (GeV)
Definition: BaryonResUtils.cxx:435
genie::utils::res::Width
double Width(Resonance_t res)
resonance width (GeV)
Definition: BaryonResUtils.cxx:461
genie::utils::bwfunc::BreitWignerL
double BreitWignerL(double W, int L, double mass, double width0, double norm)
Definition: BWFunc.cxx:107
genie::utils::res::OrbitalAngularMom
int OrbitalAngularMom(Resonance_t res)
orbital angular momentum
Definition: BaryonResUtils.cxx:524
Wmin
double Wmin[kNWBins]
Definition: piN_asymmetric_decay.C:39
sum
Double_t sum
Definition: plot.C:31
W
#define W(x)
genie::utils::res::CacheBWNorm::cache
std::map< genie::EResonance, double > cache
Definition: BaryonResUtils.h:73
genie::utils::res::ResonanceIndex
int ResonanceIndex(Resonance_t res)
resonance idx, quark model / SU(6)
Definition: BaryonResUtils.cxx:555
Resonance_t genie::utils::res::FromPdgCode ( int  pdgc)

PDG code -> resonance id.

Definition at line 98 of file BaryonResUtils.cxx.

References genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF35_1905, genie::kF37_1950, genie::kNoResonance, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1920, genie::kPdgD13m1520_N0, genie::kPdgD13m1520_NP, genie::kPdgD13m1700_N0, genie::kPdgD13m1700_NP, genie::kPdgD15m1675_N0, genie::kPdgD15m1675_NP, genie::kPdgD33m1700_Delta0, genie::kPdgD33m1700_DeltaM, genie::kPdgD33m1700_DeltaP, genie::kPdgD33m1700_DeltaPP, genie::kPdgF15m1680_N0, genie::kPdgF15m1680_NP, genie::kPdgF35m1905_Delta0, genie::kPdgF35m1905_DeltaM, genie::kPdgF35m1905_DeltaP, genie::kPdgF35m1905_DeltaPP, genie::kPdgF37m1950_Delta0, genie::kPdgF37m1950_DeltaM, genie::kPdgF37m1950_DeltaP, genie::kPdgF37m1950_DeltaPP, genie::kPdgP11m1440_N0, genie::kPdgP11m1440_NP, genie::kPdgP11m1710_N0, genie::kPdgP11m1710_NP, genie::kPdgP13m1720_N0, genie::kPdgP13m1720_NP, genie::kPdgP31m1910_Delta0, genie::kPdgP31m1910_DeltaM, genie::kPdgP31m1910_DeltaP, genie::kPdgP31m1910_DeltaPP, genie::kPdgP33m1232_Delta0, genie::kPdgP33m1232_DeltaM, genie::kPdgP33m1232_DeltaP, genie::kPdgP33m1232_DeltaPP, genie::kPdgP33m1920_Delta0, genie::kPdgP33m1920_DeltaM, genie::kPdgP33m1920_DeltaP, genie::kPdgP33m1920_DeltaPP, genie::kPdgS11m1535_N0, genie::kPdgS11m1535_NP, genie::kPdgS11m1650_N0, genie::kPdgS11m1650_NP, genie::kPdgS31m1620_Delta0, genie::kPdgS31m1620_DeltaM, genie::kPdgS31m1620_DeltaP, genie::kPdgS31m1620_DeltaPP, genie::kS11_1535, genie::kS11_1650, and genie::kS31_1620.

Referenced by genie::BaryonResonanceDecayer::EvolveDeltaDecayWidth().

99 {
100  // Delta(1232) code from the standard PDG table.
101  // Higher resonance codes come from MINOS extensions to PDG tables.
102 
103  switch(pdgc) {
104 
105  case (kPdgP33m1232_DeltaM ) : /* Delta- */
106  case (kPdgP33m1232_Delta0 ) : /* Delta0 */
107  case (kPdgP33m1232_DeltaP ) : /* Delta+ */
108  case (kPdgP33m1232_DeltaPP) : /* Delta++ */
109  return kP33_1232; break;
110 
111  case (kPdgS11m1535_N0) : /* N0 */
112  case (kPdgS11m1535_NP) : /* N+ */
113  return kS11_1535; break;
114 
115  case (kPdgD13m1520_N0) : /* N0 */
116  case (kPdgD13m1520_NP) : /* N+ */
117  return kD13_1520; break;
118 
119  case (kPdgS11m1650_N0) : /* N0 */
120  case (kPdgS11m1650_NP) : /* N+ */
121  return kS11_1650; break;
122 
123  case (kPdgD13m1700_N0) : /* N0 */
124  case (kPdgD13m1700_NP) : /* N+ */
125  return kD13_1700; break;
126 
127  case (kPdgD15m1675_N0) : /* N0 */
128  case (kPdgD15m1675_NP) : /* N+ */
129  return kD15_1675; break;
130 
131  case (kPdgS31m1620_DeltaM ) : /* Delta- */
132  case (kPdgS31m1620_Delta0 ) : /* Delta0 */
133  case (kPdgS31m1620_DeltaP ) : /* Delta+ */
134  case (kPdgS31m1620_DeltaPP) : /* Delta++ */
135  return kS31_1620; break;
136 
137  case (kPdgD33m1700_DeltaM ) : /* Delta- */
138  case (kPdgD33m1700_Delta0 ) : /* Delta0 */
139  case (kPdgD33m1700_DeltaP ) : /* Delta+ */
140  case (kPdgD33m1700_DeltaPP) : /* Delta++ */
141  return kD33_1700; break;
142 
143  case (kPdgP11m1440_N0) : /* N0 */
144  case (kPdgP11m1440_NP) : /* N+ */
145  return kP11_1440; break;
146 
147  case (kPdgP13m1720_N0) : /* N0 */
148  case (kPdgP13m1720_NP) : /* N+ */
149  return kP13_1720; break;
150 
151  case (kPdgF15m1680_N0) : /* N0 */
152  case (kPdgF15m1680_NP) : /* N+ */
153  return kF15_1680; break;
154 
155  case (kPdgP31m1910_DeltaM ) : /* Delta- */
156  case (kPdgP31m1910_Delta0 ) : /* Delta0 */
157  case (kPdgP31m1910_DeltaP ) : /* Delta+ */
158  case (kPdgP31m1910_DeltaPP) : /* Delta++ */
159  return kP31_1910; break;
160 
161  case (kPdgP33m1920_DeltaM ) : /* Delta- */
162  case (kPdgP33m1920_Delta0 ) : /* Delta0 */
163  case (kPdgP33m1920_DeltaP ) : /* Delta+ */
164  case (kPdgP33m1920_DeltaPP) : /* Delta++ */
165  return kP33_1920; break;
166 
167  case (kPdgF35m1905_DeltaM ) : /* Delta- */
168  case (kPdgF35m1905_Delta0 ) : /* Delta0 */
169  case (kPdgF35m1905_DeltaP ) : /* Delta+ */
170  case (kPdgF35m1905_DeltaPP) : /* Delta++ */
171  return kF35_1905; break;
172 
173  case (kPdgF37m1950_DeltaM ) : /* Delta- */
174  case (kPdgF37m1950_Delta0 ) : /* Delta0 */
175  case (kPdgF37m1950_DeltaP ) : /* Delta+ */
176  case (kPdgF37m1950_DeltaPP) : /* Delta++ */
177  return kF37_1950; break;
178 
179  case (kPdgP11m1710_N0) : /* N0 */
180  case (kPdgP11m1710_NP) : /* N+ */
181  return kP11_1710; break;
182  }
183 
184  return kNoResonance;
185 }
genie::kPdgP31m1910_DeltaP
const int kPdgP31m1910_DeltaP
Definition: PDGCodes.h:118
genie::kPdgS31m1620_DeltaM
const int kPdgS31m1620_DeltaM
Definition: PDGCodes.h:102
genie::kPdgP33m1232_DeltaPP
const int kPdgP33m1232_DeltaPP
Definition: PDGCodes.h:91
genie::kPdgF35m1905_DeltaM
const int kPdgF35m1905_DeltaM
Definition: PDGCodes.h:124
genie::kPdgD15m1675_N0
const int kPdgD15m1675_N0
Definition: PDGCodes.h:100
genie::kPdgF35m1905_Delta0
const int kPdgF35m1905_Delta0
Definition: PDGCodes.h:125
genie::kPdgF15m1680_NP
const int kPdgF15m1680_NP
Definition: PDGCodes.h:115
genie::kPdgD33m1700_Delta0
const int kPdgD33m1700_Delta0
Definition: PDGCodes.h:107
genie::kPdgD13m1700_NP
const int kPdgD13m1700_NP
Definition: PDGCodes.h:99
genie::kPdgS11m1650_N0
const int kPdgS11m1650_N0
Definition: PDGCodes.h:96
genie::kPdgF37m1950_DeltaM
const int kPdgF37m1950_DeltaM
Definition: PDGCodes.h:128
genie::kPdgF37m1950_DeltaP
const int kPdgF37m1950_DeltaP
Definition: PDGCodes.h:130
genie::kPdgS11m1535_N0
const int kPdgS11m1535_N0
Definition: PDGCodes.h:92
genie::kPdgP33m1232_DeltaP
const int kPdgP33m1232_DeltaP
Definition: PDGCodes.h:90
genie::kPdgP33m1232_DeltaM
const int kPdgP33m1232_DeltaM
Definition: PDGCodes.h:88
genie::kPdgP31m1910_DeltaM
const int kPdgP31m1910_DeltaM
Definition: PDGCodes.h:116
genie::kPdgF37m1950_Delta0
const int kPdgF37m1950_Delta0
Definition: PDGCodes.h:129
genie::kPdgD33m1700_DeltaPP
const int kPdgD33m1700_DeltaPP
Definition: PDGCodes.h:109
genie::kPdgF35m1905_DeltaP
const int kPdgF35m1905_DeltaP
Definition: PDGCodes.h:126
genie::kPdgD15m1675_NP
const int kPdgD15m1675_NP
Definition: PDGCodes.h:101
genie::kPdgF35m1905_DeltaPP
const int kPdgF35m1905_DeltaPP
Definition: PDGCodes.h:127
genie::kPdgS31m1620_DeltaPP
const int kPdgS31m1620_DeltaPP
Definition: PDGCodes.h:105
genie::kPdgP31m1910_DeltaPP
const int kPdgP31m1910_DeltaPP
Definition: PDGCodes.h:119
genie::kPdgP33m1920_DeltaPP
const int kPdgP33m1920_DeltaPP
Definition: PDGCodes.h:123
genie::kPdgP31m1910_Delta0
const int kPdgP31m1910_Delta0
Definition: PDGCodes.h:117
genie::kPdgP33m1920_DeltaP
const int kPdgP33m1920_DeltaP
Definition: PDGCodes.h:122
genie::kPdgD33m1700_DeltaM
const int kPdgD33m1700_DeltaM
Definition: PDGCodes.h:106
genie::kPdgP33m1920_DeltaM
const int kPdgP33m1920_DeltaM
Definition: PDGCodes.h:120
genie::kPdgF37m1950_DeltaPP
const int kPdgF37m1950_DeltaPP
Definition: PDGCodes.h:131
genie::kPdgD13m1520_NP
const int kPdgD13m1520_NP
Definition: PDGCodes.h:95
genie::kPdgF15m1680_N0
const int kPdgF15m1680_N0
Definition: PDGCodes.h:114
genie::kPdgP11m1710_N0
const int kPdgP11m1710_N0
Definition: PDGCodes.h:132
genie::kPdgS11m1650_NP
const int kPdgS11m1650_NP
Definition: PDGCodes.h:97
genie::kPdgP33m1232_Delta0
const int kPdgP33m1232_Delta0
Definition: PDGCodes.h:89
genie::kPdgP13m1720_N0
const int kPdgP13m1720_N0
Definition: PDGCodes.h:112
genie::kPdgP11m1440_N0
const int kPdgP11m1440_N0
Definition: PDGCodes.h:110
genie::kPdgP33m1920_Delta0
const int kPdgP33m1920_Delta0
Definition: PDGCodes.h:121
genie::kPdgD13m1520_N0
const int kPdgD13m1520_N0
Definition: PDGCodes.h:94
genie::kPdgP11m1440_NP
const int kPdgP11m1440_NP
Definition: PDGCodes.h:111
genie::kPdgS31m1620_Delta0
const int kPdgS31m1620_Delta0
Definition: PDGCodes.h:103
genie::kPdgP11m1710_NP
const int kPdgP11m1710_NP
Definition: PDGCodes.h:133
genie::kPdgS11m1535_NP
const int kPdgS11m1535_NP
Definition: PDGCodes.h:93
genie::kPdgS31m1620_DeltaP
const int kPdgS31m1620_DeltaP
Definition: PDGCodes.h:104
genie::kPdgP13m1720_NP
const int kPdgP13m1720_NP
Definition: PDGCodes.h:113
genie::kPdgD13m1700_N0
const int kPdgD13m1700_N0
Definition: PDGCodes.h:98
genie::kPdgD33m1700_DeltaP
const int kPdgD33m1700_DeltaP
Definition: PDGCodes.h:108
Resonance_t genie::utils::res::FromString ( const char *  res)

string -> resonance id

Definition at line 75 of file BaryonResUtils.cxx.

References genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kNoResonance, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kS11_1535, genie::kS11_1650, and genie::kS31_1620.

Referenced by genie::BaryonResList::DecodeFromNameList().

76 {
77  if ( strcmp( res,"P33(1232)" ) == 0 ) return kP33_1232;
78  else if( strcmp( res,"S11(1535)" ) == 0 ) return kS11_1535;
79  else if( strcmp( res,"D13(1520)" ) == 0 ) return kD13_1520;
80  else if( strcmp( res,"S11(1650)" ) == 0 ) return kS11_1650;
81  else if( strcmp( res,"D13(1700)" ) == 0 ) return kD13_1700;
82  else if( strcmp( res,"D15(1675)" ) == 0 ) return kD15_1675;
83  else if( strcmp( res,"S31(1620)" ) == 0 ) return kS31_1620;
84  else if( strcmp( res,"D33(1700)" ) == 0 ) return kD33_1700;
85  else if( strcmp( res,"P11(1440)" ) == 0 ) return kP11_1440;
86  else if( strcmp( res,"P33(1600)" ) == 0 ) return kP33_1600;
87  else if( strcmp( res,"P13(1720)" ) == 0 ) return kP13_1720;
88  else if( strcmp( res,"F15(1680)" ) == 0 ) return kF15_1680;
89  else if( strcmp( res,"P31(1910)" ) == 0 ) return kP31_1910;
90  else if( strcmp( res,"P33(1920)" ) == 0 ) return kP33_1920;
91  else if( strcmp( res,"F35(1905)" ) == 0 ) return kF35_1905;
92  else if( strcmp( res,"F37(1950)" ) == 0 ) return kF37_1950;
93  else if( strcmp( res,"P11(1710)" ) == 0 ) return kP11_1710;
94  else if( strcmp( res,"F17(1970)" ) == 0 ) return kF17_1970;
95  else return kNoResonance;
96 }
bool genie::utils::res::IsBaryonResonance ( int  pdgc)

is input a baryon resonance?

Definition at line 302 of file BaryonResUtils.cxx.

References genie::kPdgD13m1520_N0, genie::kPdgD13m1520_NP, genie::kPdgD13m1700_N0, genie::kPdgD13m1700_NP, genie::kPdgD15m1675_N0, genie::kPdgD15m1675_NP, genie::kPdgD33m1700_Delta0, genie::kPdgD33m1700_DeltaM, genie::kPdgD33m1700_DeltaP, genie::kPdgD33m1700_DeltaPP, genie::kPdgF15m1680_N0, genie::kPdgF15m1680_NP, genie::kPdgF35m1905_Delta0, genie::kPdgF35m1905_DeltaM, genie::kPdgF35m1905_DeltaP, genie::kPdgF35m1905_DeltaPP, genie::kPdgF37m1950_Delta0, genie::kPdgF37m1950_DeltaM, genie::kPdgF37m1950_DeltaP, genie::kPdgF37m1950_DeltaPP, genie::kPdgP11m1440_N0, genie::kPdgP11m1440_NP, genie::kPdgP11m1710_N0, genie::kPdgP11m1710_NP, genie::kPdgP13m1720_N0, genie::kPdgP13m1720_NP, genie::kPdgP31m1910_Delta0, genie::kPdgP31m1910_DeltaM, genie::kPdgP31m1910_DeltaP, genie::kPdgP31m1910_DeltaPP, genie::kPdgP33m1232_Delta0, genie::kPdgP33m1232_DeltaM, genie::kPdgP33m1232_DeltaP, genie::kPdgP33m1232_DeltaPP, genie::kPdgP33m1920_Delta0, genie::kPdgP33m1920_DeltaM, genie::kPdgP33m1920_DeltaP, genie::kPdgP33m1920_DeltaPP, genie::kPdgS11m1535_N0, genie::kPdgS11m1535_NP, genie::kPdgS11m1650_N0, genie::kPdgS11m1650_NP, genie::kPdgS31m1620_Delta0, genie::kPdgS31m1620_DeltaM, genie::kPdgS31m1620_DeltaP, and genie::kPdgS31m1620_DeltaPP.

Referenced by genie::pdg::IsBaryonResonance(), genie::PythiaDecayer::IsHandled(), genie::BaryonResonanceDecayer::IsHandled(), and genie::Decayer::IsUnstable().

303 {
304  // Delta(1232) code from the standard PDG table.
305  // Higher resonance codes come from MINOS extensions to PDG tables.
306 
307  switch(pdgc) {
308 
309  /* ------ P33(1232) ------*/
310  case (kPdgP33m1232_DeltaM ) : /* Delta- */
311  case (kPdgP33m1232_Delta0 ) : /* Delta0 */
312  case (kPdgP33m1232_DeltaP ) : /* Delta+ */
313  case (kPdgP33m1232_DeltaPP) : /* Delta++ */
314 
315  /* ------ S11(1535) ------*/
316  case (kPdgS11m1535_N0) : /* N0 */
317  case (kPdgS11m1535_NP) : /* N+ */
318 
319  /* ------ D13(1520) ------*/
320  case (kPdgD13m1520_N0) : /* N0 */
321  case (kPdgD13m1520_NP) : /* N+ */
322 
323  /* ------ S11(1650) ------*/
324  case (kPdgS11m1650_N0) : /* N0 */
325  case (kPdgS11m1650_NP) : /* N+ */
326 
327  /* ------ D13(1700) ------*/
328  case (kPdgD13m1700_N0) : /* N0 */
329  case (kPdgD13m1700_NP) : /* N+ */
330 
331  /* ------ D15(1675) ------*/
332  case (kPdgD15m1675_N0) : /* N0 */
333  case (kPdgD15m1675_NP) : /* N+ */
334 
335  /* ------ S31(1620) ------*/
336  case (kPdgS31m1620_DeltaM ) : /* Delta- */
337  case (kPdgS31m1620_Delta0 ) : /* Delta0 */
338  case (kPdgS31m1620_DeltaP ) : /* Delta+ */
339  case (kPdgS31m1620_DeltaPP) : /* Delta++ */
340 
341  /* ------ D33(1700) ------*/
342  case (kPdgD33m1700_DeltaM ) : /* Delta- */
343  case (kPdgD33m1700_Delta0 ) : /* Delta0 */
344  case (kPdgD33m1700_DeltaP ) : /* Delta+ */
345  case (kPdgD33m1700_DeltaPP) : /* Delta++ */
346 
347  /* ------ P11(1440) ------*/
348  case (kPdgP11m1440_N0) : /* N0 */
349  case (kPdgP11m1440_NP) : /* N+ */
350 
351  /* ------ P33(1600) ------*/
352  // are you?
353 
354  /* ------ P13(1720) ------*/
355  case (kPdgP13m1720_N0) : /* N0 */
356  case (kPdgP13m1720_NP) : /* N+ */
357 
358  /* ------ F15(1680) ------*/
359  case (kPdgF15m1680_N0) : /* N0 */
360  case (kPdgF15m1680_NP) : /* N+ */
361 
362  /* ------ P31(1910) ------*/
363  case (kPdgP31m1910_DeltaM ) : /* Delta- */
364  case (kPdgP31m1910_Delta0 ) : /* Delta0 */
365  case (kPdgP31m1910_DeltaP ) : /* Delta+ */
366  case (kPdgP31m1910_DeltaPP) : /* Delta++ */
367 
368  /* ------ P33(1920) ------*/
369  case (kPdgP33m1920_DeltaM ) : /* Delta- */
370  case (kPdgP33m1920_Delta0 ) : /* Delta0 */
371  case (kPdgP33m1920_DeltaP ) : /* Delta+ */
372  case (kPdgP33m1920_DeltaPP) : /* Delta++ */
373 
374  /* ------ F35(1905) ------*/
375  case (kPdgF35m1905_DeltaM ) : /* Delta- */
376  case (kPdgF35m1905_Delta0 ) : /* Delta0 */
377  case (kPdgF35m1905_DeltaP ) : /* Delta+ */
378  case (kPdgF35m1905_DeltaPP) : /* Delta++ */
379 
380  /* ------ F37(1950) ------*/
381  case (kPdgF37m1950_DeltaM ) : /* Delta- */
382  case (kPdgF37m1950_Delta0 ) : /* Delta0 */
383  case (kPdgF37m1950_DeltaP ) : /* Delta+ */
384  case (kPdgF37m1950_DeltaPP) : /* Delta++ */
385 
386  /* ------ P11(1710) ------*/
387  case (kPdgP11m1710_N0) : /* N0 */
388  case (kPdgP11m1710_NP) : /* N+ */
389 
390  /* ------ F17(1970) ------*/
391  // are you?
392 
393  return true;
394  }
395 
396  return false;
397 }
genie::kPdgP31m1910_DeltaP
const int kPdgP31m1910_DeltaP
Definition: PDGCodes.h:118
genie::kPdgS31m1620_DeltaM
const int kPdgS31m1620_DeltaM
Definition: PDGCodes.h:102
genie::kPdgP33m1232_DeltaPP
const int kPdgP33m1232_DeltaPP
Definition: PDGCodes.h:91
genie::kPdgF35m1905_DeltaM
const int kPdgF35m1905_DeltaM
Definition: PDGCodes.h:124
genie::kPdgD15m1675_N0
const int kPdgD15m1675_N0
Definition: PDGCodes.h:100
genie::kPdgF35m1905_Delta0
const int kPdgF35m1905_Delta0
Definition: PDGCodes.h:125
genie::kPdgF15m1680_NP
const int kPdgF15m1680_NP
Definition: PDGCodes.h:115
genie::kPdgD33m1700_Delta0
const int kPdgD33m1700_Delta0
Definition: PDGCodes.h:107
genie::kPdgD13m1700_NP
const int kPdgD13m1700_NP
Definition: PDGCodes.h:99
genie::kPdgS11m1650_N0
const int kPdgS11m1650_N0
Definition: PDGCodes.h:96
genie::kPdgF37m1950_DeltaM
const int kPdgF37m1950_DeltaM
Definition: PDGCodes.h:128
genie::kPdgF37m1950_DeltaP
const int kPdgF37m1950_DeltaP
Definition: PDGCodes.h:130
genie::kPdgS11m1535_N0
const int kPdgS11m1535_N0
Definition: PDGCodes.h:92
genie::kPdgP33m1232_DeltaP
const int kPdgP33m1232_DeltaP
Definition: PDGCodes.h:90
genie::kPdgP33m1232_DeltaM
const int kPdgP33m1232_DeltaM
Definition: PDGCodes.h:88
genie::kPdgP31m1910_DeltaM
const int kPdgP31m1910_DeltaM
Definition: PDGCodes.h:116
genie::kPdgF37m1950_Delta0
const int kPdgF37m1950_Delta0
Definition: PDGCodes.h:129
genie::kPdgD33m1700_DeltaPP
const int kPdgD33m1700_DeltaPP
Definition: PDGCodes.h:109
genie::kPdgF35m1905_DeltaP
const int kPdgF35m1905_DeltaP
Definition: PDGCodes.h:126
genie::kPdgD15m1675_NP
const int kPdgD15m1675_NP
Definition: PDGCodes.h:101
genie::kPdgF35m1905_DeltaPP
const int kPdgF35m1905_DeltaPP
Definition: PDGCodes.h:127
genie::kPdgS31m1620_DeltaPP
const int kPdgS31m1620_DeltaPP
Definition: PDGCodes.h:105
genie::kPdgP31m1910_DeltaPP
const int kPdgP31m1910_DeltaPP
Definition: PDGCodes.h:119
genie::kPdgP33m1920_DeltaPP
const int kPdgP33m1920_DeltaPP
Definition: PDGCodes.h:123
genie::kPdgP31m1910_Delta0
const int kPdgP31m1910_Delta0
Definition: PDGCodes.h:117
genie::kPdgP33m1920_DeltaP
const int kPdgP33m1920_DeltaP
Definition: PDGCodes.h:122
genie::kPdgD33m1700_DeltaM
const int kPdgD33m1700_DeltaM
Definition: PDGCodes.h:106
genie::kPdgP33m1920_DeltaM
const int kPdgP33m1920_DeltaM
Definition: PDGCodes.h:120
genie::kPdgF37m1950_DeltaPP
const int kPdgF37m1950_DeltaPP
Definition: PDGCodes.h:131
genie::kPdgD13m1520_NP
const int kPdgD13m1520_NP
Definition: PDGCodes.h:95
genie::kPdgF15m1680_N0
const int kPdgF15m1680_N0
Definition: PDGCodes.h:114
genie::kPdgP11m1710_N0
const int kPdgP11m1710_N0
Definition: PDGCodes.h:132
genie::kPdgS11m1650_NP
const int kPdgS11m1650_NP
Definition: PDGCodes.h:97
genie::kPdgP33m1232_Delta0
const int kPdgP33m1232_Delta0
Definition: PDGCodes.h:89
genie::kPdgP13m1720_N0
const int kPdgP13m1720_N0
Definition: PDGCodes.h:112
genie::kPdgP11m1440_N0
const int kPdgP11m1440_N0
Definition: PDGCodes.h:110
genie::kPdgP33m1920_Delta0
const int kPdgP33m1920_Delta0
Definition: PDGCodes.h:121
genie::kPdgD13m1520_N0
const int kPdgD13m1520_N0
Definition: PDGCodes.h:94
genie::kPdgP11m1440_NP
const int kPdgP11m1440_NP
Definition: PDGCodes.h:111
genie::kPdgS31m1620_Delta0
const int kPdgS31m1620_Delta0
Definition: PDGCodes.h:103
genie::kPdgP11m1710_NP
const int kPdgP11m1710_NP
Definition: PDGCodes.h:133
genie::kPdgS11m1535_NP
const int kPdgS11m1535_NP
Definition: PDGCodes.h:93
genie::kPdgS31m1620_DeltaP
const int kPdgS31m1620_DeltaP
Definition: PDGCodes.h:104
genie::kPdgP13m1720_NP
const int kPdgP13m1720_NP
Definition: PDGCodes.h:113
genie::kPdgD13m1700_N0
const int kPdgD13m1700_N0
Definition: PDGCodes.h:98
genie::kPdgD33m1700_DeltaP
const int kPdgD33m1700_DeltaP
Definition: PDGCodes.h:108
bool genie::utils::res::IsDelta ( Resonance_t  res)

is it a Delta resonance?

Definition at line 399 of file BaryonResUtils.cxx.

References ana::assert(), genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kS11_1535, genie::kS11_1650, and genie::kS31_1620.

Referenced by genie::RESInteractionListGenerator::CreateInteractionList(), genie::GiBUURESFormFactor::FormFactors::F1V(), genie::GiBUURESFormFactor::FormFactors::F2V(), genie::GiBUURESFormFactor::FormFactors::FA(), genie::GiBUURESFormFactor::FormFactors::FP(), IsN(), genie::SppChannel::IsospinWeight(), genie::GiBUURESFormFactor::FormFactors::LoadTables(), SaveToRootFile(), genie::RSPPResonanceSelector::SelectResonance(), genie::ReinSehgalRESPXSec::XSec(), and genie::BSKLNBaseRESPXSec2014::XSec().

400 {
401  switch(res) {
402 
403  case kP33_1232: return true; break;
404  case kS11_1535: return false; break;
405  case kD13_1520: return false; break;
406  case kS11_1650: return false; break;
407  case kD13_1700: return false; break;
408  case kD15_1675: return false; break;
409  case kS31_1620: return true; break;
410  case kD33_1700: return true; break;
411  case kP11_1440: return false; break;
412  case kP33_1600: return true; break;
413  case kP13_1720: return false; break;
414  case kF15_1680: return false; break;
415  case kP31_1910: return true; break;
416  case kP33_1920: return true; break;
417  case kF35_1905: return true; break;
418  case kF37_1950: return true; break;
419  case kP11_1710: return false; break;
420  case kF17_1970: return false; break;
421  default:
422  // should not be here - meaningless to return anything
423  assert(false);
424  }
425  return false;
426 }
ana::assert
assert(nhit_max >=nhit_nbins)
bool genie::utils::res::IsN ( Resonance_t  res)

is it an N resonance?

Definition at line 430 of file BaryonResUtils.cxx.

References IsDelta().

Referenced by genie::GiBUURESFormFactor::FormFactors::C3A(), genie::GiBUURESFormFactor::FormFactors::C3V(), genie::GiBUURESFormFactor::FormFactors::C4A(), genie::GiBUURESFormFactor::FormFactors::C4V(), genie::GiBUURESFormFactor::FormFactors::C5A(), genie::GiBUURESFormFactor::FormFactors::C5V(), genie::GiBUURESFormFactor::FormFactors::C6A(), genie::GiBUURESFormFactor::FormFactors::C6V(), and genie::GiBUURESFormFactor::FormFactors::LoadTables().

431 {
432  return (! utils::res::IsDelta(res) );
433 }
genie::utils::res::IsDelta
bool IsDelta(Resonance_t res)
is it a Delta resonance?
Definition: BaryonResUtils.cxx:399
double genie::utils::res::Mass ( Resonance_t  res)

resonance mass (GeV)

Definition at line 435 of file BaryonResUtils.cxx.

References genie::units::GeV, genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kS11_1535, genie::kS11_1650, and genie::kS31_1620.

Referenced by genie::RSPPHadronicSystemGenerator::AddResonanceDecayProducts(), genie::InitialState::AsString(), genie::utils::BindHitNucleon(), BWNorm(), genie::utils::ComputeFullQELPXSec(), genie::RESKinematicsGenerator::ComputeMaxXSec(), genie::SmithMonizQELCCPXSec::d2sQES_dQ2dv_SM(), genie::utils::gsl::d2XSecRESFast_dWQ2_E::d2XSecRESFast_dWQ2_E(), genie::HAIntranuke::ElasHA(), genie::HAIntranuke2018::ElasHA(), genie::BaryonResonanceDecayer::EvolveDeltaDecayWidth(), genie::NucleonDecayPrimaryVtxGenerator::GenerateFermiMomentum(), genie::MECGenerator::GenerateFermiMomentum(), genie::KNOHadronization::GenerateHadronCodes(), genie::MECGenerator::GenerateNSVInitialHadrons(), nuesand::FillNueSandbox::GetECF(), InitializeEvent(), genie::LwlynSmithQELCCPXSec::Integral(), genie::AlamSimoAtharVacasSKXSec::Integrate(), genie::FermiMover::KickHitNucleon(), main(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::SmithMonizUtils::SetInteraction(), genie::utils::intranuke::ThreeBodyKinematics(), genie::utils::intranuke2018::ThreeBodyKinematics(), genie::EmpiricalMECPXSec2015::XSec(), genie::P33PaschosLalakulichPXSec::XSec(), genie::ReinSehgalRESPXSec::XSec(), and genie::BSKLNBaseRESPXSec2014::XSec().

436 {
437  switch(res) {
438  case kP33_1232 : return 1.232 * units::GeV ; break;
439  case kS11_1535 : return 1.535 * units::GeV ; break;
440  case kD13_1520 : return 1.515 * units::GeV ; break;
441  case kS11_1650 : return 1.655 * units::GeV ; break;
442  case kD13_1700 : return 1.700 * units::GeV ; break;
443  case kD15_1675 : return 1.675 * units::GeV ; break;
444  case kS31_1620 : return 1.630 * units::GeV ; break;
445  case kD33_1700 : return 1.700 * units::GeV ; break;
446  case kP11_1440 : return 1.430 * units::GeV ; break;
447  case kP33_1600 : return 1.600 * units::GeV ; break;
448  case kP13_1720 : return 1.720 * units::GeV ; break;
449  case kF15_1680 : return 1.685 * units::GeV ; break;
450  case kP31_1910 : return 1.890 * units::GeV ; break;
451  case kP33_1920 : return 1.920 * units::GeV ; break;
452  case kF35_1905 : return 1.880 * units::GeV ; break;
453  case kF37_1950 : return 1.930 * units::GeV ; break;
454  case kP11_1710 : return 1.710 * units::GeV ; break;
455  case kF17_1970 : return 2.190 * units::GeV ; break;
456  default: break;
457  }
458  return -1;
459 }
Munits::GeV
static constexpr Double_t GeV
Definition: Munits.h:291
int genie::utils::res::OrbitalAngularMom ( Resonance_t  res)

orbital angular momentum

Definition at line 524 of file BaryonResUtils.cxx.

References exit(), genie::gAbortingInErr, genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kS11_1535, genie::kS11_1650, genie::kS31_1620, LOG, and pFATAL.

Referenced by BWNorm(), main(), genie::ReinSehgalRESPXSec::XSec(), and genie::BSKLNBaseRESPXSec2014::XSec().

525 {
526  switch(res) {
527  case kP33_1232: return 1; break;
528  case kS11_1535: return 0; break;
529  case kD13_1520: return 2; break;
530  case kS11_1650: return 0; break;
531  case kD13_1700: return 2; break;
532  case kD15_1675: return 2; break;
533  case kS31_1620: return 0; break;
534  case kD33_1700: return 2; break;
535  case kP11_1440: return 1; break;
536  case kP33_1600: return 1; break;
537  case kP13_1720: return 1; break;
538  case kF15_1680: return 3; break;
539  case kP31_1910: return 1; break;
540  case kP33_1920: return 1; break;
541  case kF35_1905: return 3; break;
542  case kF37_1950: return 3; break;
543  case kP11_1710: return 1; break;
544  case kF17_1970: return 3; break;
545  default:
546  // should not be here - meaningless to return anything
547  gAbortingInErr = true;
548  LOG("BaryonRes", pFATAL)
549  << "Unknown resonance " << res;
550  exit(1);
551  }
552  return 0;
553 }
pFATAL
#define pFATAL
Definition: Messenger.h:57
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:97
exit
exit(0)
genie::gAbortingInErr
bool gAbortingInErr
Definition: Messenger.cxx:56
int genie::utils::res::PdgCode ( Resonance_t  res,
int  Q 
)

(resonance id, charge) -> PDG code

Definition at line 187 of file BaryonResUtils.cxx.

References genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kPdgD13m1520_N0, genie::kPdgD13m1520_NP, genie::kPdgD13m1700_N0, genie::kPdgD13m1700_NP, genie::kPdgD15m1675_N0, genie::kPdgD15m1675_NP, genie::kPdgD33m1700_Delta0, genie::kPdgD33m1700_DeltaM, genie::kPdgD33m1700_DeltaP, genie::kPdgD33m1700_DeltaPP, genie::kPdgF15m1680_N0, genie::kPdgF15m1680_NP, genie::kPdgF35m1905_Delta0, genie::kPdgF35m1905_DeltaM, genie::kPdgF35m1905_DeltaP, genie::kPdgF35m1905_DeltaPP, genie::kPdgF37m1950_Delta0, genie::kPdgF37m1950_DeltaM, genie::kPdgF37m1950_DeltaP, genie::kPdgF37m1950_DeltaPP, genie::kPdgP11m1440_N0, genie::kPdgP11m1440_NP, genie::kPdgP11m1710_N0, genie::kPdgP11m1710_NP, genie::kPdgP13m1720_N0, genie::kPdgP13m1720_NP, genie::kPdgP31m1910_Delta0, genie::kPdgP31m1910_DeltaM, genie::kPdgP31m1910_DeltaP, genie::kPdgP31m1910_DeltaPP, genie::kPdgP33m1232_Delta0, genie::kPdgP33m1232_DeltaM, genie::kPdgP33m1232_DeltaP, genie::kPdgP33m1232_DeltaPP, genie::kPdgP33m1920_Delta0, genie::kPdgP33m1920_DeltaM, genie::kPdgP33m1920_DeltaP, genie::kPdgP33m1920_DeltaPP, genie::kPdgS11m1535_N0, genie::kPdgS11m1535_NP, genie::kPdgS11m1650_N0, genie::kPdgS11m1650_NP, genie::kPdgS31m1620_Delta0, genie::kPdgS31m1620_DeltaM, genie::kPdgS31m1620_DeltaP, genie::kPdgS31m1620_DeltaPP, genie::kS11_1535, genie::kS11_1650, and genie::kS31_1620.

Referenced by caf::AddParticleToVec(), genie::RSPPResonanceSelector::AddResonance(), murem::MuonRemoveAna::analyze(), slicer::SlicerAna::analyze(), fuzz::FuzzyKValidate::analyze(), bpfit::BPFTmvaTrainer::analyze(), slid::LIDTraining::analyze(), caf::FillParticleTruth(), genie::GMCJDriver::GenerateEvent1Try(), genie::GMCJDriver::GenerateFluxNeutrino(), nuesand::FillNueSandbox::GetECF(), numusand::NumuSandFxs::getMissingE(), genie::RESHadronicSystemGenerator::GetResonancePdgCode(), evtsum::EventSummary::MakeOutput(), genie::flux::GDk2NuFlux::MaxEnergy(), numusand::FillSandbox::produce(), cosrej::CosRej::produce(), simb::MCTruth::SetNeutrino(), and cheat::BackTracker::TracksToParticles().

188 {
189  // Delta(1232) code from the standard PDG table
190  // Higher resonance codes come from MINOS extensions to PDG tables
191 
192  switch(res) {
193 
194  case kP33_1232:
195  if(Q == -1) return kPdgP33m1232_DeltaM; /* Delta- */
196  if(Q == 0) return kPdgP33m1232_Delta0; /* Delta0 */
197  if(Q == 1) return kPdgP33m1232_DeltaP; /* Delta+ */
198  if(Q == 2) return kPdgP33m1232_DeltaPP; /* Delta++ */
199  break;
200 
201  case kS11_1535:
202  if(Q == 0) return kPdgS11m1535_N0; /* N0 */
203  if(Q == 1) return kPdgS11m1535_NP; /* N+ */
204  break;
205 
206  case kD13_1520:
207  if(Q == 0) return kPdgD13m1520_N0; /* N0 */
208  if(Q == 1) return kPdgD13m1520_NP; /* N+ */
209  break;
210 
211  case kS11_1650:
212  if(Q == 0) return kPdgS11m1650_N0; /* N0 */
213  if(Q == 1) return kPdgS11m1650_NP; /* N+ */
214  break;
215 
216  case kD13_1700:
217  if(Q == 0) return kPdgD13m1700_N0; /* N0 */
218  if(Q == 1) return kPdgD13m1700_NP; /* N+ */
219  break;
220 
221  case kD15_1675:
222  if(Q == 0) return kPdgD15m1675_N0; /* N0 */
223  if(Q == 1) return kPdgD15m1675_NP; /* N+ */
224  break;
225 
226  case kS31_1620:
227  if(Q == -1) return kPdgS31m1620_DeltaM; /* Delta- */
228  if(Q == 0) return kPdgS31m1620_Delta0; /* Delta0 */
229  if(Q == 1) return kPdgS31m1620_DeltaP; /* Delta+ */
230  if(Q == 2) return kPdgS31m1620_DeltaPP; /* Delta++ */
231  break;
232 
233  case kD33_1700:
234  if(Q == -1) return kPdgD33m1700_DeltaM; /* Delta- */
235  if(Q == 0) return kPdgD33m1700_Delta0; /* Delta0 */
236  if(Q == 1) return kPdgD33m1700_DeltaP; /* Delta+ */
237  if(Q == 2) return kPdgD33m1700_DeltaPP; /* Delta++ */
238  break;
239 
240  case kP11_1440:
241  if(Q == 0) return kPdgP11m1440_N0; /* N0 */
242  if(Q == 1) return kPdgP11m1440_NP; /* N+ */
243  break;
244 
245  case kP33_1600:
246  return 0;
247 
248  case kP13_1720:
249  if(Q == 0) return kPdgP13m1720_N0; /* N0 */
250  if(Q == 1) return kPdgP13m1720_NP; /* N+ */
251  break;
252 
253  case kF15_1680:
254  if(Q == 0) return kPdgF15m1680_N0; /* N0 */
255  if(Q == 1) return kPdgF15m1680_NP; /* N+ */
256  break;
257 
258  case kP31_1910:
259  if(Q == -1) return kPdgP31m1910_DeltaM; /* Delta- */
260  if(Q == 0) return kPdgP31m1910_Delta0; /* Delta0 */
261  if(Q == 1) return kPdgP31m1910_DeltaP; /* Delta+ */
262  if(Q == 2) return kPdgP31m1910_DeltaPP; /* Delta++ */
263  break;
264 
265  case kP33_1920:
266  if(Q == -1) return kPdgP33m1920_DeltaM; /* Delta- */
267  if(Q == 0) return kPdgP33m1920_Delta0; /* Delta0 */
268  if(Q == 1) return kPdgP33m1920_DeltaP; /* Delta+ */
269  if(Q == 2) return kPdgP33m1920_DeltaPP; /* Delta++ */
270  break;
271 
272  case kF35_1905:
273  if(Q == -1) return kPdgF35m1905_DeltaM; /* Delta- */
274  if(Q == 0) return kPdgF35m1905_Delta0; /* Delta0 */
275  if(Q == 1) return kPdgF35m1905_DeltaP; /* Delta+ */
276  if(Q == 2) return kPdgF35m1905_DeltaPP; /* Delta++ */
277  break;
278 
279  case kF37_1950:
280  if(Q == -1) return kPdgF37m1950_DeltaM; /* Delta- */
281  if(Q == 0) return kPdgF37m1950_Delta0; /* Delta0 */
282  if(Q == 1) return kPdgF37m1950_DeltaP; /* Delta+ */
283  if(Q == 2) return kPdgF37m1950_DeltaPP; /* Delta++ */
284  break;
285 
286  case kP11_1710:
287  if(Q == 0) return kPdgP11m1710_N0; /* N0 */
288  if(Q == 1) return kPdgP11m1710_NP; /* N+ */
289  break;
290 
291  case kF17_1970:
292  return 0;
293  break;
294 
295  default:
296  return 0;
297  }
298 
299  return 0;
300 }
genie::kPdgP31m1910_DeltaP
const int kPdgP31m1910_DeltaP
Definition: PDGCodes.h:118
genie::kPdgS31m1620_DeltaM
const int kPdgS31m1620_DeltaM
Definition: PDGCodes.h:102
genie::kPdgP33m1232_DeltaPP
const int kPdgP33m1232_DeltaPP
Definition: PDGCodes.h:91
genie::kPdgF35m1905_DeltaM
const int kPdgF35m1905_DeltaM
Definition: PDGCodes.h:124
genie::kPdgD15m1675_N0
const int kPdgD15m1675_N0
Definition: PDGCodes.h:100
genie::kPdgF35m1905_Delta0
const int kPdgF35m1905_Delta0
Definition: PDGCodes.h:125
genie::kPdgF15m1680_NP
const int kPdgF15m1680_NP
Definition: PDGCodes.h:115
genie::kPdgD33m1700_Delta0
const int kPdgD33m1700_Delta0
Definition: PDGCodes.h:107
genie::kPdgD13m1700_NP
const int kPdgD13m1700_NP
Definition: PDGCodes.h:99
genie::kPdgS11m1650_N0
const int kPdgS11m1650_N0
Definition: PDGCodes.h:96
genie::kPdgF37m1950_DeltaM
const int kPdgF37m1950_DeltaM
Definition: PDGCodes.h:128
genie::kPdgF37m1950_DeltaP
const int kPdgF37m1950_DeltaP
Definition: PDGCodes.h:130
genie::kPdgS11m1535_N0
const int kPdgS11m1535_N0
Definition: PDGCodes.h:92
genie::kPdgP33m1232_DeltaP
const int kPdgP33m1232_DeltaP
Definition: PDGCodes.h:90
genie::kPdgP33m1232_DeltaM
const int kPdgP33m1232_DeltaM
Definition: PDGCodes.h:88
genie::kPdgP31m1910_DeltaM
const int kPdgP31m1910_DeltaM
Definition: PDGCodes.h:116
genie::kPdgF37m1950_Delta0
const int kPdgF37m1950_Delta0
Definition: PDGCodes.h:129
genie::kPdgD33m1700_DeltaPP
const int kPdgD33m1700_DeltaPP
Definition: PDGCodes.h:109
genie::kPdgF35m1905_DeltaP
const int kPdgF35m1905_DeltaP
Definition: PDGCodes.h:126
genie::kPdgD15m1675_NP
const int kPdgD15m1675_NP
Definition: PDGCodes.h:101
genie::kPdgF35m1905_DeltaPP
const int kPdgF35m1905_DeltaPP
Definition: PDGCodes.h:127
genie::kPdgS31m1620_DeltaPP
const int kPdgS31m1620_DeltaPP
Definition: PDGCodes.h:105
genie::kPdgP31m1910_DeltaPP
const int kPdgP31m1910_DeltaPP
Definition: PDGCodes.h:119
genie::kPdgP33m1920_DeltaPP
const int kPdgP33m1920_DeltaPP
Definition: PDGCodes.h:123
genie::kPdgP31m1910_Delta0
const int kPdgP31m1910_Delta0
Definition: PDGCodes.h:117
genie::kPdgP33m1920_DeltaP
const int kPdgP33m1920_DeltaP
Definition: PDGCodes.h:122
genie::kPdgD33m1700_DeltaM
const int kPdgD33m1700_DeltaM
Definition: PDGCodes.h:106
genie::kPdgP33m1920_DeltaM
const int kPdgP33m1920_DeltaM
Definition: PDGCodes.h:120
genie::kPdgF37m1950_DeltaPP
const int kPdgF37m1950_DeltaPP
Definition: PDGCodes.h:131
genie::kPdgD13m1520_NP
const int kPdgD13m1520_NP
Definition: PDGCodes.h:95
genie::kPdgF15m1680_N0
const int kPdgF15m1680_N0
Definition: PDGCodes.h:114
genie::kPdgP11m1710_N0
const int kPdgP11m1710_N0
Definition: PDGCodes.h:132
genie::kPdgS11m1650_NP
const int kPdgS11m1650_NP
Definition: PDGCodes.h:97
genie::kPdgP33m1232_Delta0
const int kPdgP33m1232_Delta0
Definition: PDGCodes.h:89
genie::kPdgP13m1720_N0
const int kPdgP13m1720_N0
Definition: PDGCodes.h:112
genie::kPdgP11m1440_N0
const int kPdgP11m1440_N0
Definition: PDGCodes.h:110
genie::kPdgP33m1920_Delta0
const int kPdgP33m1920_Delta0
Definition: PDGCodes.h:121
genie::kPdgD13m1520_N0
const int kPdgD13m1520_N0
Definition: PDGCodes.h:94
genie::kPdgP11m1440_NP
const int kPdgP11m1440_NP
Definition: PDGCodes.h:111
genie::kPdgS31m1620_Delta0
const int kPdgS31m1620_Delta0
Definition: PDGCodes.h:103
genie::kPdgP11m1710_NP
const int kPdgP11m1710_NP
Definition: PDGCodes.h:133
genie::kPdgS11m1535_NP
const int kPdgS11m1535_NP
Definition: PDGCodes.h:93
genie::kPdgS31m1620_DeltaP
const int kPdgS31m1620_DeltaP
Definition: PDGCodes.h:104
genie::kPdgP13m1720_NP
const int kPdgP13m1720_NP
Definition: PDGCodes.h:113
genie::kPdgD13m1700_N0
const int kPdgD13m1700_N0
Definition: PDGCodes.h:98
genie::kPdgD33m1700_DeltaP
const int kPdgD33m1700_DeltaP
Definition: PDGCodes.h:108
int genie::utils::res::ResonanceIndex ( Resonance_t  res)

resonance idx, quark model / SU(6)

Definition at line 555 of file BaryonResUtils.cxx.

References exit(), genie::gAbortingInErr, genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kS11_1535, genie::kS11_1650, genie::kS31_1620, LOG, and pFATAL.

Referenced by BWNorm(), genie::utils::gsl::d2XSecRESFast_dWQ2_E::d2XSecRESFast_dWQ2_E(), genie::ReinSehgalRESPXSec::XSec(), and genie::BSKLNBaseRESPXSec2014::XSec().

556 {
557  switch(res) {
558  case kP33_1232: return 0; break;
559  case kS11_1535: return 1; break;
560  case kD13_1520: return 1; break;
561  case kS11_1650: return 1; break;
562  case kD13_1700: return 1; break;
563  case kD15_1675: return 1; break;
564  case kS31_1620: return 1; break;
565  case kD33_1700: return 1; break;
566  case kP11_1440: return 2; break;
567  case kP33_1600: return 2; break;
568  case kP13_1720: return 2; break;
569  case kF15_1680: return 2; break;
570  case kP31_1910: return 2; break;
571  case kP33_1920: return 2; break;
572  case kF35_1905: return 2; break;
573  case kF37_1950: return 2; break;
574  case kP11_1710: return 2; break;
575  case kF17_1970: return 2; break;
576  default:
577  // should not be here - meaningless to return anything
578  gAbortingInErr = true;
579  LOG("BaryonRes", pFATAL)
580  << "Unknown resonance " << res;
581  exit(1);
582  }
583  return 0;
584 }
pFATAL
#define pFATAL
Definition: Messenger.h:57
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:97
exit
exit(0)
genie::gAbortingInErr
bool gAbortingInErr
Definition: Messenger.cxx:56
double genie::utils::res::Width ( Resonance_t  res)

resonance width (GeV)

Definition at line 461 of file BaryonResUtils.cxx.

References genie::units::GeV, genie::kD13_1520, genie::kD13_1700, genie::kD15_1675, genie::kD33_1700, genie::kF15_1680, genie::kF17_1970, genie::kF35_1905, genie::kF37_1950, genie::kP11_1440, genie::kP11_1710, genie::kP13_1720, genie::kP31_1910, genie::kP33_1232, genie::kP33_1600, genie::kP33_1920, genie::kS11_1535, genie::kS11_1650, and genie::kS31_1620.

Referenced by BWNorm(), genie::utils::gsl::d2XSecRESFast_dWQ2_E::d2XSecRESFast_dWQ2_E(), genie::BaryonResonanceDecayer::EvolveDeltaDecayWidth(), main(), genie::P33PaschosLalakulichPXSec::XSec(), genie::ReinSehgalRESPXSec::XSec(), and genie::BSKLNBaseRESPXSec2014::XSec().

462 {
463  switch(res) {
464  case kP33_1232 : return 0.117 * units::GeV ; break;
465  case kS11_1535 : return 0.150 * units::GeV ; break;
466  case kD13_1520 : return 0.115 * units::GeV ; break;
467  case kS11_1650 : return 0.140 * units::GeV ; break;
468  case kD13_1700 : return 0.150 * units::GeV ; break;
469  case kD15_1675 : return 0.150 * units::GeV ; break;
470  case kS31_1620 : return 0.140 * units::GeV ; break;
471  case kD33_1700 : return 0.300 * units::GeV ; break;
472  case kP11_1440 : return 0.350 * units::GeV ; break;
473  case kP33_1600 : return 0.320 * units::GeV ; break;
474  case kP13_1720 : return 0.250 * units::GeV ; break;
475  case kF15_1680 : return 0.130 * units::GeV ; break;
476  case kP31_1910 : return 0.280 * units::GeV ; break;
477  case kP33_1920 : return 0.260 * units::GeV ; break;
478  case kF35_1905 : return 0.330 * units::GeV ; break;
479  case kF37_1950 : return 0.285 * units::GeV ; break;
480  case kP11_1710 : return 0.100 * units::GeV ; break;
481  case kF17_1970 : return 0.500 * units::GeV ; break;
482  default: break;
483  }
484  return -1;
485 }
Munits::GeV
static constexpr Double_t GeV
Definition: Munits.h:291
Generated on Fri Nov 27 2020 05:26:13 for NOvA by   doxygen 1.8.11