Public Member Functions | Static Public Attributes | Private Member Functions | Private Attributes | List of all members
genie::rew::GReWeightNuXSecDIS Class Reference

Reweighting GENIE DIS neutrino-nucleus cross sections. More...

#include <GReWeightNuXSecDIS.h>

Inheritance diagram for genie::rew::GReWeightNuXSecDIS:
genie::rew::GReWeightI

Public Member Functions

 GReWeightNuXSecDIS ()
 
 ~GReWeightNuXSecDIS ()
 
bool IsHandled (GSyst_t syst)
 does the current weight calculator handle the input nuisance param? More...
 
void SetSystematic (GSyst_t syst, double val)
 update the value for the specified nuisance param More...
 
void Reset (void)
 set all nuisance parameters to default values More...
 
void Reconfigure (void)
 propagate updated nuisance parameter values to actual MC, etc More...
 
double CalcWeight (const EventRecord &event)
 calculate a weight for the input event using the current nuisance param values More...
 
void SetMode (int m)
 
void RewNue (bool tf)
 
void RewNuebar (bool tf)
 
void RewNumu (bool tf)
 
void RewNumubar (bool tf)
 
void RewCC (bool tf)
 
void RewNC (bool tf)
 
void SetAhtBYPath (string p)
 
void SetBhtBYPath (string p)
 
void SetCV1uBYPath (string p)
 
void SetCV2uBYPath (string p)
 
void SetWminCut (double W)
 
void SetQ2minCut (double Q2)
 
- Public Member Functions inherited from genie::rew::GReWeightI
virtual ~GReWeightI ()
 

Static Public Attributes

static const int kModeABCV12u = 0
 
static const int kModeABCV12uShape = 1
 

Private Member Functions

void Init (void)
 
double CalcWeightABCV12u (const genie::EventRecord &event)
 rew. Aht,Bht,CV1u,CV2u More...
 
double CalcWeightABCV12uShape (const genie::EventRecord &event)
 rew. AhtShape,BhtShape,CV1uShape,CV2uShape More...
 

Private Attributes

XSecAlgorithmIfXSecModelDef
 default model More...
 
XSecAlgorithmIfXSecModel
 tweaked model More...
 
RegistryfXSecModelConfig
 config in tweaked model More...
 
bool fRewNue
 reweight nu_e? More...
 
bool fRewNuebar
 reweight nu_e_bar? More...
 
bool fRewNumu
 reweight nu_mu? More...
 
bool fRewNumubar
 reweight nu_mu_bar? More...
 
bool fRewCC
 reweight CC? More...
 
bool fRewNC
 reweight NC? More...
 
int fMode
 0: Aht,Bht,CV1u,CV2u, 1:AhtShape,BhtShape,CV1uShape,CV2uShape More...
 
double fWmin
 W_{min} cut. Reweight only events with W > W_{min}. More...
 
double fQ2min
 Q2_{min} cut. Reweight only events with Q2 > Q2_{min}. More...
 
double fAhtBYTwkDial
 tweak dial for BY parameter: Aht More...
 
double fBhtBYTwkDial
 tweak dial for BY parameter: Bht More...
 
double fCV1uBYTwkDial
 tweak dial for BY parameter: CV1u More...
 
double fCV2uBYTwkDial
 tweak dial for BY parameter: CV2u More...
 
double fAhtBYDef
 
double fBhtBYDef
 
double fCV1uBYDef
 
double fCV2uBYDef
 
double fAhtBYCur
 
double fBhtBYCur
 
double fCV1uBYCur
 
double fCV2uBYCur
 
string fAhtBYPath
 
string fBhtBYPath
 
string fCV1uBYPath
 
string fCV2uBYPath
 

Additional Inherited Members

- Protected Member Functions inherited from genie::rew::GReWeightI
 GReWeightI ()
 

Detailed Description

Reweighting GENIE DIS neutrino-nucleus cross sections.

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

Jim Dobson <J.Dobson07 imperial.ac.uk> Imperial College London

Aug 1, 2009

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

Definition at line 39 of file GReWeightNuXSecDIS.h.

Constructor & Destructor Documentation

GReWeightNuXSecDIS::GReWeightNuXSecDIS ( )

Definition at line 56 of file GReWeightNuXSecDIS.cxx.

57 {
58  this->Init();
59 }
GReWeightNuXSecDIS::~GReWeightNuXSecDIS ( )

Definition at line 61 of file GReWeightNuXSecDIS.cxx.

62 {
63 #ifdef _G_REWEIGHT_DIS_DEBUG_
64  fTestFile->cd();
65  fTestNtp ->Write();
66  fTestFile->Close();
67  delete fTestFile;
68 #endif
69 }

Member Function Documentation

double GReWeightNuXSecDIS::CalcWeight ( const EventRecord event)
virtual

calculate a weight for the input event using the current nuisance param values

Implements genie::rew::GReWeightI.

Definition at line 179 of file GReWeightNuXSecDIS.cxx.

180 {
181  bool tweaked =
182  (TMath::Abs(fAhtBYTwkDial) > controls::kASmallNum) ||
183  (TMath::Abs(fBhtBYTwkDial) > controls::kASmallNum) ||
184  (TMath::Abs(fCV1uBYTwkDial) > controls::kASmallNum) ||
185  (TMath::Abs(fCV2uBYTwkDial) > controls::kASmallNum);
186  if(!tweaked) return 1.0;
187 
188  Interaction * interaction = event.Summary();
189 
190  bool is_dis = interaction->ProcInfo().IsDeepInelastic();
191  if(!is_dis) return 1.;
192 
193  bool charm = event.Summary()->ExclTag().IsCharmEvent(); // skip DIS charm
194  if(charm) return 1.;
195 
196  bool is_cc = interaction->ProcInfo().IsWeakCC();
197  bool is_nc = interaction->ProcInfo().IsWeakNC();
198  if(is_cc && !fRewCC) return 1.;
199  if(is_nc && !fRewNC) return 1.;
200 
201  int nupdg = interaction->InitState().ProbePdg();
202  if(nupdg==kPdgNuMu && !fRewNumu ) return 1.;
203  if(nupdg==kPdgAntiNuMu && !fRewNumubar) return 1.;
204  if(nupdg==kPdgNuE && !fRewNue ) return 1.;
205  if(nupdg==kPdgAntiNuE && !fRewNuebar ) return 1.;
206 
207  bool selected = true;
208  double W = interaction->Kine().W (selected);
209  double Q2 = interaction->Kine().Q2(selected);
210  bool passes_kine_cuts = (W>=fWmin && Q2>=fQ2min);
211  if(!passes_kine_cuts) return 1.;
212 
213  //
214  // calculate weight
215  //
216 
217  double wght = 1.;
218 
219  if(fMode == kModeABCV12u) {
220  wght = this->CalcWeightABCV12u(event);
221  }
222  else
223  if(fMode == kModeABCV12uShape) {
224  wght = this->CalcWeightABCV12uShape(event);
225  }
226 
227 #ifdef _G_REWEIGHT_DIS_DEBUG_
228  double E = interaction->InitState().ProbeE(kRfHitNucRest);
229  double x = interaction->Kine().x(true);
230  double y = interaction->Kine().y(true);
231  int ccnc = (is_cc) ? 1 : 0;
232  int nuc = interaction->InitState().Tgt().HitNucPdg();
233  int qrk = interaction->InitState().Tgt().HitQrkPdg();
234  int sea = (interaction->InitState().Tgt().HitSeaQrk()) ? 1 : 0;
235  fTestNtp->Fill(E,x,y,nupdg,nuc,qrk,sea,ccnc,wght);
236 #endif
237 
238  return wght;
239 }
double W(bool selected=false) const
Definition: Kinematics.cxx:167
bool HitSeaQrk(void) const
Definition: Target.cxx:306
bool IsWeakCC(void) const
const int kPdgNuE
Definition: PDGCodes.h:25
double fBhtBYTwkDial
tweak dial for BY parameter: Bht
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:642
int HitNucPdg(void) const
Definition: Target.cxx:311
int HitQrkPdg(void) const
Definition: Target.cxx:249
const int kPdgAntiNuE
Definition: PDGCodes.h:26
const int kPdgNuMu
Definition: PDGCodes.h:27
double x(bool selected=false) const
Definition: Kinematics.cxx:109
double y(bool selected=false) const
Definition: Kinematics.cxx:122
Summary information for an interaction.
Definition: Interaction.h:53
double y
double fCV1uBYTwkDial
tweak dial for BY parameter: CV1u
bool IsWeakNC(void) const
int fMode
0: Aht,Bht,CV1u,CV2u, 1:AhtShape,BhtShape,CV1uShape,CV2uShape
bool fRewNuebar
reweight nu_e_bar?
const Kinematics & Kine(void) const
Definition: Interaction.h:68
int ProbePdg(void) const
Definition: InitialState.h:54
static const double kASmallNum
Definition: Controls.h:40
const int kPdgAntiNuMu
Definition: PDGCodes.h:28
double fCV2uBYTwkDial
tweak dial for BY parameter: CV2u
double fWmin
W_{min} cut. Reweight only events with W > W_{min}.
double fAhtBYTwkDial
tweak dial for BY parameter: Aht
bool IsDeepInelastic(void) const
Definition: ProcessInfo.cxx:98
double fQ2min
Q2_{min} cut. Reweight only events with Q2 > Q2_{min}.
double CalcWeightABCV12u(const genie::EventRecord &event)
rew. Aht,Bht,CV1u,CV2u
const InitialState & InitState(void) const
Definition: Interaction.h:66
const ProcessInfo & ProcInfo(void) const
Definition: Interaction.h:67
double Q2(bool selected=false) const
Definition: Kinematics.cxx:135
const Target & Tgt(void) const
Definition: InitialState.h:56
double CalcWeightABCV12uShape(const genie::EventRecord &event)
rew. AhtShape,BhtShape,CV1uShape,CV2uShape
double ProbeE(RefFrame_t rf) const
bool fRewNumubar
reweight nu_mu_bar?
double GReWeightNuXSecDIS::CalcWeightABCV12u ( const genie::EventRecord event)
private

rew. Aht,Bht,CV1u,CV2u

Definition at line 241 of file GReWeightNuXSecDIS.cxx.

242 {
243  Interaction * interaction = event.Summary();
244 
245  interaction->KinePtr()->UseSelectedKinematics();
246 
247  double old_xsec = event.DiffXSec();
248  double old_weight = event.Weight();
249  double twk_xsec = fXSecModel->XSec(interaction, kPSxyfE);
250  double weight = old_weight * (twk_xsec/old_xsec);
251 
252  interaction->KinePtr()->ClearRunningValues();
253 
254  return weight;
255 }
Kinematics * KinePtr(void) const
Definition: Interaction.h:73
virtual double XSec(const Interaction *i, KinePhaseSpace_t k=kPSfE) const =0
Compute the cross section for the input interaction.
XSecAlgorithmI * fXSecModel
tweaked model
Summary information for an interaction.
Definition: Interaction.h:53
void UseSelectedKinematics(void)
Definition: Kinematics.cxx:369
weight
Definition: test.py:293
void ClearRunningValues(void)
Definition: Kinematics.cxx:357
double GReWeightNuXSecDIS::CalcWeightABCV12uShape ( const genie::EventRecord event)
private

rew. AhtShape,BhtShape,CV1uShape,CV2uShape

Definition at line 257 of file GReWeightNuXSecDIS.cxx.

258 {
259  Interaction * interaction = event.Summary();
260 
261  interaction->KinePtr()->UseSelectedKinematics();
262 
263  double old_xsec = event.DiffXSec();
264  double old_weight = event.Weight();
265  double twk_xsec = fXSecModel->XSec(interaction, kPSxyfE);
266  double weight = old_weight * (twk_xsec/old_xsec);
267 
268 //double old_integrated_xsec = event.XSec();
269  double old_integrated_xsec = fXSecModelDef -> Integral(interaction);
270  double twk_integrated_xsec = fXSecModel -> Integral(interaction);
271 
272  assert(twk_integrated_xsec > 0);
273  weight *= (old_integrated_xsec/twk_integrated_xsec);
274 
275  interaction->KinePtr()->ClearRunningValues();
276 
277  return weight;
278 }
Kinematics * KinePtr(void) const
Definition: Interaction.h:73
virtual double XSec(const Interaction *i, KinePhaseSpace_t k=kPSfE) const =0
Compute the cross section for the input interaction.
XSecAlgorithmI * fXSecModelDef
default model
XSecAlgorithmI * fXSecModel
tweaked model
Summary information for an interaction.
Definition: Interaction.h:53
void UseSelectedKinematics(void)
Definition: Kinematics.cxx:369
weight
Definition: test.py:293
void ClearRunningValues(void)
Definition: Kinematics.cxx:357
void GReWeightNuXSecDIS::Init ( void  )
private

Definition at line 280 of file GReWeightNuXSecDIS.cxx.

281 {
282  AlgId id("genie::QPMDISPXSec","Default");
283 
284  AlgFactory * algf = AlgFactory::Instance();
285 
286  Algorithm * algdef = algf->AdoptAlgorithm(id);
287  fXSecModelDef = dynamic_cast<XSecAlgorithmI*>(algdef);
289 
290  Algorithm * alg = algf->AdoptAlgorithm(id);
291  fXSecModel = dynamic_cast<XSecAlgorithmI*>(alg);
293 
295 //LOG("ReW", pNOTICE) << *fXSecModelConfig;
296 
297  this->SetMode (kModeABCV12u);
298 
299  this->RewNue (true);
300  this->RewNuebar (true);
301  this->RewNumu (true);
302  this->RewNumubar(true);
303  this->RewCC (true);
304  this->RewNC (true);
305 
306  this->SetWminCut (1.7*units::GeV);
307  this->SetQ2minCut(1.0*units::GeV*units::GeV);
308 
309  this->SetAhtBYPath ("SFAlg/A");
310  this->SetBhtBYPath ("SFAlg/B");
311  this->SetCV1uBYPath("SFAlg/Cv1U");
312  this->SetCV2uBYPath("SFAlg/Cv2U");
313 
314  fAhtBYTwkDial = 0;
315  fBhtBYTwkDial = 0;
316  fCV1uBYTwkDial = 0;
317  fCV2uBYTwkDial = 0;
322  fAhtBYCur = fAhtBYDef;
323  fBhtBYCur = fBhtBYDef;
326 
327 #ifdef _G_REWEIGHT_DIS_DEBUG_
328  fTestFile = new TFile("./dis_reweight_test.root","recreate");
329  fTestNtp = new TNtupleD("testntp","","E:x:y:nu:nuc:qrk:sea:ccnc:wght");
330 #endif
331 }
virtual const Registry & GetConfig(void) const
Get configuration registry.
Definition: Algorithm.h:63
Cross Section Calculation Interface.
double fBhtBYTwkDial
tweak dial for BY parameter: Bht
Algorithm abstract base class.
Definition: Algorithm.h:48
RgDbl GetDouble(RgKey key) const
Definition: Registry.cxx:488
XSecAlgorithmI * fXSecModelDef
default model
XSecAlgorithmI * fXSecModel
tweaked model
double fCV1uBYTwkDial
tweak dial for BY parameter: CV1u
void AdoptSubstructure(void)
Definition: Algorithm.cxx:287
Algorithm * AdoptAlgorithm(const AlgId &algid) const
Definition: AlgFactory.cxx:127
double fCV2uBYTwkDial
tweak dial for BY parameter: CV2u
double fAhtBYTwkDial
tweak dial for BY parameter: Aht
Algorithm ID (algorithm name + configuration set name)
Definition: AlgId.h:35
static AlgFactory * Instance()
Definition: AlgFactory.cxx:75
A registry. Provides the container for algorithm configuration parameters.
Definition: Registry.h:66
Registry * fXSecModelConfig
config in tweaked model
The GENIE Algorithm Factory.
Definition: AlgFactory.h:40
static const double GeV
Definition: Units.h:29
bool GReWeightNuXSecDIS::IsHandled ( GSyst_t  syst)
virtual

does the current weight calculator handle the input nuisance param?

Implements genie::rew::GReWeightI.

Definition at line 71 of file GReWeightNuXSecDIS.cxx.

72 {
73  bool handle = false;
74 
75  switch(syst) {
76  case ( kXSecTwkDial_AhtBYshape ) :
77  case ( kXSecTwkDial_BhtBYshape ) :
78  case ( kXSecTwkDial_CV1uBYshape ) :
79  case ( kXSecTwkDial_CV2uBYshape ) :
80 
81  if (fMode == kModeABCV12u ) handle = false;
82  else if (fMode == kModeABCV12uShape) handle = true;
83  break;
84 
85  case ( kXSecTwkDial_AhtBY ) :
86  case ( kXSecTwkDial_BhtBY ) :
87  case ( kXSecTwkDial_CV1uBY ) :
88  case ( kXSecTwkDial_CV2uBY ) :
89 
90  if (fMode == kModeABCV12u ) handle = true;
91  else if (fMode == kModeABCV12uShape) handle = false;
92  break;
93 
94  default:
95  handle = false;
96  break;
97  }
98 
99  return handle;
100 }
tweak the Bodek-Yang model parameter B_{ht} - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:91
int fMode
0: Aht,Bht,CV1u,CV2u, 1:AhtShape,BhtShape,CV1uShape,CV2uShape
tweak the Bodek-Yang model parameter CV2u - incl. both shape and normalization effect ...
Definition: GSyst.h:89
tweak the Bodek-Yang model parameter A_{ht} - incl. both shape and normalization effect ...
Definition: GSyst.h:86
tweak the Bodek-Yang model parameter A_{ht} - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:90
tweak the Bodek-Yang model parameter CV1u - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:92
tweak the Bodek-Yang model parameter CV1u - incl. both shape and normalization effect ...
Definition: GSyst.h:88
tweak the Bodek-Yang model parameter CV2u - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:93
tweak the Bodek-Yang model parameter B_{ht} - incl. both shape and normalization effect ...
Definition: GSyst.h:87
void GReWeightNuXSecDIS::Reconfigure ( void  )
virtual

propagate updated nuisance parameter values to actual MC, etc

Implements genie::rew::GReWeightI.

Definition at line 139 of file GReWeightNuXSecDIS.cxx.

140 {
142 
143  double fracerr_aht = 0.;
144  double fracerr_bht = 0.;
145  double fracerr_cv1u = 0.;
146  double fracerr_cv2u = 0.;
147 
148  if(fMode == kModeABCV12u) {
149  fracerr_aht = fracerr->OneSigmaErr (kXSecTwkDial_AhtBY );
150  fracerr_bht = fracerr->OneSigmaErr (kXSecTwkDial_BhtBY );
151  fracerr_cv1u = fracerr->OneSigmaErr (kXSecTwkDial_CV1uBY);
152  fracerr_cv2u = fracerr->OneSigmaErr (kXSecTwkDial_CV2uBY);
153  }
154  else
155  if(fMode == kModeABCV12uShape) {
156  fracerr_aht = fracerr->OneSigmaErr (kXSecTwkDial_AhtBYshape );
157  fracerr_bht = fracerr->OneSigmaErr (kXSecTwkDial_BhtBYshape );
158  fracerr_cv1u = fracerr->OneSigmaErr (kXSecTwkDial_CV1uBYshape);
159  fracerr_cv2u = fracerr->OneSigmaErr (kXSecTwkDial_CV2uBYshape);
160  }
161 
162  fAhtBYCur = fAhtBYDef * (1. + fAhtBYTwkDial * fracerr_aht );
163  fBhtBYCur = fBhtBYDef * (1. + fBhtBYTwkDial * fracerr_bht );
164  fCV1uBYCur = fCV1uBYDef * (1. + fCV1uBYTwkDial * fracerr_cv1u);
165  fCV2uBYCur = fCV2uBYDef * (1. + fCV2uBYTwkDial * fracerr_cv2u);
166 
168 
169  r.Set(fAhtBYPath, fAhtBYCur );
170  r.Set(fBhtBYPath, fBhtBYCur );
171  r.Set(fCV1uBYPath, fCV1uBYCur);
172  r.Set(fCV2uBYPath, fCV2uBYCur);
173 
174  fXSecModel->Configure(r);
175 
176 //LOG("ReW", pDEBUG) << *fXSecModel;
177 }
double fBhtBYTwkDial
tweak dial for BY parameter: Bht
tweak the Bodek-Yang model parameter B_{ht} - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:91
XSecAlgorithmI * fXSecModel
tweaked model
double OneSigmaErr(GSyst_t syst, int sign=0) const
double fCV1uBYTwkDial
tweak dial for BY parameter: CV1u
int fMode
0: Aht,Bht,CV1u,CV2u, 1:AhtShape,BhtShape,CV1uShape,CV2uShape
tweak the Bodek-Yang model parameter CV2u - incl. both shape and normalization effect ...
Definition: GSyst.h:89
virtual void Configure(const Registry &config)
Configure the algorithm.
Definition: Algorithm.cxx:70
double fCV2uBYTwkDial
tweak dial for BY parameter: CV2u
double fAhtBYTwkDial
tweak dial for BY parameter: Aht
tweak the Bodek-Yang model parameter A_{ht} - incl. both shape and normalization effect ...
Definition: GSyst.h:86
tweak the Bodek-Yang model parameter A_{ht} - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:90
A registry. Provides the container for algorithm configuration parameters.
Definition: Registry.h:66
Registry * fXSecModelConfig
config in tweaked model
tweak the Bodek-Yang model parameter CV1u - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:92
static GSystUncertainty * Instance(void)
tweak the Bodek-Yang model parameter CV1u - incl. both shape and normalization effect ...
Definition: GSyst.h:88
tweak the Bodek-Yang model parameter CV2u - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:93
void Set(RgIMapPair entry)
Definition: Registry.cxx:281
tweak the Bodek-Yang model parameter B_{ht} - incl. both shape and normalization effect ...
Definition: GSyst.h:87
void GReWeightNuXSecDIS::Reset ( void  )
virtual

set all nuisance parameters to default values

Implements genie::rew::GReWeightI.

Definition at line 129 of file GReWeightNuXSecDIS.cxx.

130 {
131  fAhtBYTwkDial = 0.;
132  fBhtBYTwkDial = 0.;
133  fCV1uBYTwkDial = 0.;
134  fCV2uBYTwkDial = 0.;
135 
136  this->Reconfigure();
137 }
double fBhtBYTwkDial
tweak dial for BY parameter: Bht
void Reconfigure(void)
propagate updated nuisance parameter values to actual MC, etc
double fCV1uBYTwkDial
tweak dial for BY parameter: CV1u
double fCV2uBYTwkDial
tweak dial for BY parameter: CV2u
double fAhtBYTwkDial
tweak dial for BY parameter: Aht
void genie::rew::GReWeightNuXSecDIS::RewCC ( bool  tf)
inline

Definition at line 61 of file GReWeightNuXSecDIS.h.

61 { fRewCC = tf; }
Definition: tf_graph.h:23
void genie::rew::GReWeightNuXSecDIS::RewNC ( bool  tf)
inline

Definition at line 62 of file GReWeightNuXSecDIS.h.

62 { fRewNC = tf; }
Definition: tf_graph.h:23
void genie::rew::GReWeightNuXSecDIS::RewNue ( bool  tf)
inline

Definition at line 57 of file GReWeightNuXSecDIS.h.

57 { fRewNue = tf; }
Definition: tf_graph.h:23
void genie::rew::GReWeightNuXSecDIS::RewNuebar ( bool  tf)
inline

Definition at line 58 of file GReWeightNuXSecDIS.h.

58 { fRewNuebar = tf; }
Definition: tf_graph.h:23
bool fRewNuebar
reweight nu_e_bar?
void genie::rew::GReWeightNuXSecDIS::RewNumu ( bool  tf)
inline

Definition at line 59 of file GReWeightNuXSecDIS.h.

59 { fRewNumu = tf; }
Definition: tf_graph.h:23
void genie::rew::GReWeightNuXSecDIS::RewNumubar ( bool  tf)
inline

Definition at line 60 of file GReWeightNuXSecDIS.h.

60 { fRewNumubar = tf; }
Definition: tf_graph.h:23
bool fRewNumubar
reweight nu_mu_bar?
void genie::rew::GReWeightNuXSecDIS::SetAhtBYPath ( string  p)
inline

Definition at line 63 of file GReWeightNuXSecDIS.h.

void genie::rew::GReWeightNuXSecDIS::SetBhtBYPath ( string  p)
inline

Definition at line 64 of file GReWeightNuXSecDIS.h.

void genie::rew::GReWeightNuXSecDIS::SetCV1uBYPath ( string  p)
inline

Definition at line 65 of file GReWeightNuXSecDIS.h.

void genie::rew::GReWeightNuXSecDIS::SetCV2uBYPath ( string  p)
inline

Definition at line 66 of file GReWeightNuXSecDIS.h.

void genie::rew::GReWeightNuXSecDIS::SetMode ( int  m)
inline

Definition at line 56 of file GReWeightNuXSecDIS.h.

56 { fMode = m; }
static const double m
Definition: Units.h:79
int fMode
0: Aht,Bht,CV1u,CV2u, 1:AhtShape,BhtShape,CV1uShape,CV2uShape
void genie::rew::GReWeightNuXSecDIS::SetQ2minCut ( double  Q2)
inline

Definition at line 68 of file GReWeightNuXSecDIS.h.

68 { fQ2min = Q2; }
double Q2(const Interaction *const i)
Definition: KineUtils.cxx:642
double fQ2min
Q2_{min} cut. Reweight only events with Q2 > Q2_{min}.
void GReWeightNuXSecDIS::SetSystematic ( GSyst_t  syst,
double  val 
)
virtual

update the value for the specified nuisance param

Implements genie::rew::GReWeightI.

Definition at line 102 of file GReWeightNuXSecDIS.cxx.

103 {
104  if(! this->IsHandled(syst)) return;
105 
106  switch(syst) {
107  case ( kXSecTwkDial_AhtBY ) :
108  case ( kXSecTwkDial_AhtBYshape ) :
109  fAhtBYTwkDial = twk_dial;
110  break;
111  case ( kXSecTwkDial_BhtBY ) :
112  case ( kXSecTwkDial_BhtBYshape ) :
113  fBhtBYTwkDial = twk_dial;
114  break;
115  case ( kXSecTwkDial_CV1uBY ) :
116  case ( kXSecTwkDial_CV1uBYshape ) :
117  fCV1uBYTwkDial = twk_dial;
118  break;
119  case ( kXSecTwkDial_CV2uBY ) :
120  case ( kXSecTwkDial_CV2uBYshape ) :
121  fCV2uBYTwkDial = twk_dial;
122  break;
123  default:
124  return;
125  break;
126  }
127 }
bool IsHandled(GSyst_t syst)
does the current weight calculator handle the input nuisance param?
double fBhtBYTwkDial
tweak dial for BY parameter: Bht
tweak the Bodek-Yang model parameter B_{ht} - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:91
double fCV1uBYTwkDial
tweak dial for BY parameter: CV1u
tweak the Bodek-Yang model parameter CV2u - incl. both shape and normalization effect ...
Definition: GSyst.h:89
double fCV2uBYTwkDial
tweak dial for BY parameter: CV2u
double fAhtBYTwkDial
tweak dial for BY parameter: Aht
tweak the Bodek-Yang model parameter A_{ht} - incl. both shape and normalization effect ...
Definition: GSyst.h:86
tweak the Bodek-Yang model parameter A_{ht} - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:90
tweak the Bodek-Yang model parameter CV1u - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:92
tweak the Bodek-Yang model parameter CV1u - incl. both shape and normalization effect ...
Definition: GSyst.h:88
tweak the Bodek-Yang model parameter CV2u - shape only effect to d2sigma(DIS)/dxdy ...
Definition: GSyst.h:93
tweak the Bodek-Yang model parameter B_{ht} - incl. both shape and normalization effect ...
Definition: GSyst.h:87
void genie::rew::GReWeightNuXSecDIS::SetWminCut ( double  W)
inline

Definition at line 67 of file GReWeightNuXSecDIS.h.

67 { fWmin = W; }
double fWmin
W_{min} cut. Reweight only events with W > W_{min}.

Member Data Documentation

double genie::rew::GReWeightNuXSecDIS::fAhtBYCur
private

Definition at line 98 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fAhtBYDef
private

Definition at line 94 of file GReWeightNuXSecDIS.h.

string genie::rew::GReWeightNuXSecDIS::fAhtBYPath
private

Definition at line 102 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fAhtBYTwkDial
private

tweak dial for BY parameter: Aht

Definition at line 90 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fBhtBYCur
private

Definition at line 99 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fBhtBYDef
private

Definition at line 95 of file GReWeightNuXSecDIS.h.

string genie::rew::GReWeightNuXSecDIS::fBhtBYPath
private

Definition at line 103 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fBhtBYTwkDial
private

tweak dial for BY parameter: Bht

Definition at line 91 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fCV1uBYCur
private

Definition at line 100 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fCV1uBYDef
private

Definition at line 96 of file GReWeightNuXSecDIS.h.

string genie::rew::GReWeightNuXSecDIS::fCV1uBYPath
private

Definition at line 104 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fCV1uBYTwkDial
private

tweak dial for BY parameter: CV1u

Definition at line 92 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fCV2uBYCur
private

Definition at line 101 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fCV2uBYDef
private

Definition at line 97 of file GReWeightNuXSecDIS.h.

string genie::rew::GReWeightNuXSecDIS::fCV2uBYPath
private

Definition at line 105 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fCV2uBYTwkDial
private

tweak dial for BY parameter: CV2u

Definition at line 93 of file GReWeightNuXSecDIS.h.

int genie::rew::GReWeightNuXSecDIS::fMode
private

0: Aht,Bht,CV1u,CV2u, 1:AhtShape,BhtShape,CV1uShape,CV2uShape

Definition at line 87 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fQ2min
private

Q2_{min} cut. Reweight only events with Q2 > Q2_{min}.

Definition at line 89 of file GReWeightNuXSecDIS.h.

bool genie::rew::GReWeightNuXSecDIS::fRewCC
private

reweight CC?

Definition at line 84 of file GReWeightNuXSecDIS.h.

bool genie::rew::GReWeightNuXSecDIS::fRewNC
private

reweight NC?

Definition at line 85 of file GReWeightNuXSecDIS.h.

bool genie::rew::GReWeightNuXSecDIS::fRewNue
private

reweight nu_e?

Definition at line 80 of file GReWeightNuXSecDIS.h.

bool genie::rew::GReWeightNuXSecDIS::fRewNuebar
private

reweight nu_e_bar?

Definition at line 81 of file GReWeightNuXSecDIS.h.

bool genie::rew::GReWeightNuXSecDIS::fRewNumu
private

reweight nu_mu?

Definition at line 82 of file GReWeightNuXSecDIS.h.

bool genie::rew::GReWeightNuXSecDIS::fRewNumubar
private

reweight nu_mu_bar?

Definition at line 83 of file GReWeightNuXSecDIS.h.

double genie::rew::GReWeightNuXSecDIS::fWmin
private

W_{min} cut. Reweight only events with W > W_{min}.

Definition at line 88 of file GReWeightNuXSecDIS.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecDIS::fXSecModel
private

tweaked model

Definition at line 77 of file GReWeightNuXSecDIS.h.

Registry* genie::rew::GReWeightNuXSecDIS::fXSecModelConfig
private

config in tweaked model

Definition at line 78 of file GReWeightNuXSecDIS.h.

XSecAlgorithmI* genie::rew::GReWeightNuXSecDIS::fXSecModelDef
private

default model

Definition at line 76 of file GReWeightNuXSecDIS.h.

const int GReWeightNuXSecDIS::kModeABCV12u = 0
static

Definition at line 42 of file GReWeightNuXSecDIS.h.

const int GReWeightNuXSecDIS::kModeABCV12uShape = 1
static

Definition at line 43 of file GReWeightNuXSecDIS.h.


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