NeugenWeightCalculator Class Reference

#include <NeugenWeightCalculator.h>

Inheritance diagram for NeugenWeightCalculator:

List of all members.

Public Member Functions
	NeugenWeightCalculator (Registry *stdconfig=0)
	~NeugenWeightCalculator ()
void	Config ()
double	GetWeight (Registry *eventinfo)
double	GetWeight (MCEventInfo *event, NuParent *parent)
void	ReweightConfigReset ()
void	PrintReweightConfig (ostream &)
void	SetOptimization (bool in)
Protected Attributes
neugen_config *	fStdNeuConfig
neugen_config *	fRwtNeuConfig
neugen_wrapper *	fNuWrap
Private Attributes
bool	fIsOptimized
bool	fUpdateConfig

---

Constructor & Destructor Documentation

NeugenWeightCalculator::NeugenWeightCalculator (  Registry *  stdconfig = 0  )

Definition at line 19 of file NeugenWeightCalculator.cxx. References fIsOptimized, fNuWrap, fRwtNeuConfig, fStdNeuConfig, fUpdateConfig, neugen_config::set_best_parameters(), WeightCalculator::SetReweightConfig(), and WeightCalculator::SetStandardConfig(). : WeightCalculator(config) { fWCname = "NeugenWeightCalculator"; fStandardConfigChanged = false; fReweightConfigChanged = false; fStdNeuConfig = new neugen_config("stdNeuConfig"); fStdNeuConfig->set_best_parameters(); if(config) SetStandardConfig(config); fRwtNeuConfig = new neugen_config("rwtNeuConfig"); fRwtNeuConfig->set_best_parameters(); if(config) SetReweightConfig(config); fNuWrap = new neugen_wrapper(fStdNeuConfig); fIsOptimized = false; fUpdateConfig = false; }

NeugenWeightCalculator::~NeugenWeightCalculator (   )

Definition at line 43 of file NeugenWeightCalculator.cxx. { delete fStdNeuConfig; delete fRwtNeuConfig; delete fNuWrap; }

---

Member Function Documentation

void NeugenWeightCalculator::Config (   )  [virtual]

Reimplemented from WeightCalculator. Definition at line 65 of file NeugenWeightCalculator.cxx. References neugen_config::coh_r0(), neugen_config::coh_rei(), neugen_config::config_name(), neugen_config::config_no(), neugen_config::dcf(), neugen_config::dis_res(), e_vbn, e_vbp, e_vn, e_vp, neugen_config::fill_named_configuration(), fRwtNeuConfig, fStdNeuConfig, fUpdateConfig, Registry::Get(), neugen_config::kno_a1(), neugen_config::kno_a2(), neugen_config::kno_a3(), neugen_config::kno_a4(), neugen_config::kno_b1(), neugen_config::kno_b2(), neugen_config::kno_b3(), neugen_config::kno_b4(), neugen_config::kno_c1(), neugen_config::kno_c2(), neugen_config::kno_c3(), neugen_config::kno_c4(), neugen_config::ma_coh(), neugen_config::ma_qe(), neugen_config::ma_res(), neugen_config::qel_eta(), neugen_config::qel_fa0(), neugen_config::res_omega(), neugen_config::res_z(), neugen_config::set_coh_r0(), neugen_config::set_coh_rei(), neugen_config::set_config_name(), neugen_config::set_config_no(), neugen_config::set_dcf(), neugen_config::set_dis_res(), neugen_config::set_kno_a1(), neugen_config::set_kno_a2(), neugen_config::set_kno_a3(), neugen_config::set_kno_a4(), neugen_config::set_kno_b1(), neugen_config::set_kno_b2(), neugen_config::set_kno_b3(), neugen_config::set_kno_b4(), neugen_config::set_kno_c1(), neugen_config::set_kno_c2(), neugen_config::set_kno_c3(), neugen_config::set_kno_c4(), neugen_config::set_ma_coh(), neugen_config::set_ma_qe(), neugen_config::set_ma_res(), neugen_config::set_nres(), neugen_config::set_pdf_group(), neugen_config::set_pdf_set(), neugen_config::set_qel_eta(), neugen_config::set_qel_fa0(), neugen_config::set_res_omega(), neugen_config::set_res_z(), neugen_config::set_wcutd(), neugen_config::set_wcutr(), neugen_config::wcutd(), and neugen_config::wcutr(). { double tempd = 0; int tempi = 0; const char *temps =0; if(fStandardConfigChanged) { bool setname = false; bool setno = false; if(fStandardConfig->Get("neugen:config_no",tempi)){ fStdNeuConfig->set_config_no(tempi); setno = true; } if(fStandardConfig->Get("neugen:config_name",temps)){ fStdNeuConfig->set_config_name(temps); setname=true; } if(setname&&!setno){ std::cout<<"You tried to set a new config name in NeugenWeightCalculator::Config()" <<" but you did not set a config no. Dire consequences will most likely result" <<std::endl; assert(0); } //if the name and number were set in the registry, //refill the standard configuration with the default values of this //named configuration //otherwise assume that the values that were put in the standardconfig //when it was created are good to start from if(setname&&setno){ fStdNeuConfig->fill_named_configuration(); } if(fStandardConfig->Get("neugen:pdf_group",tempi)) fStdNeuConfig->set_pdf_group(tempi); if(fStandardConfig->Get("neugen:pdf_set",tempi)) fStdNeuConfig->set_pdf_set(tempi); if(fStandardConfig->Get("neugen:ma_qe",tempd)) fStdNeuConfig->set_ma_qe(tempd); if(fStandardConfig->Get("neugen:ma_res",tempd)) fStdNeuConfig->set_ma_res(tempd); if(fStandardConfig->Get("neugen:ma_coh",tempd)) fStdNeuConfig->set_ma_coh(tempd); if(fStandardConfig->Get("neugen:qel_fa0",tempd)) fStdNeuConfig->set_qel_fa0(tempd); if(fStandardConfig->Get("neugen:qel_eta",tempd)) fStdNeuConfig->set_qel_eta(tempd); if(fStandardConfig->Get("neugen:res_omega",tempd)) fStdNeuConfig->set_res_omega(tempd); if(fStandardConfig->Get("neugen:res_z",tempd)) fStdNeuConfig->set_res_z(tempd); if(fStandardConfig->Get("neugen:coh_r0",tempd)) fStdNeuConfig->set_coh_r0(tempd); if(fStandardConfig->Get("neugen:coh_rei",tempd)) fStdNeuConfig->set_coh_rei(tempd); if(fStandardConfig->Get("neugen:kno_a1",tempd)) fStdNeuConfig->set_kno_a1(tempd); if(fStandardConfig->Get("neugen:kno_a2",tempd)) fStdNeuConfig->set_kno_a2(tempd); if(fStandardConfig->Get("neugen:kno_a3",tempd)) fStdNeuConfig->set_kno_a3(tempd); if(fStandardConfig->Get("neugen:kno_a4",tempd)) fStdNeuConfig->set_kno_a4(tempd); if(fStandardConfig->Get("neugen:kno_b1",tempd)) fStdNeuConfig->set_kno_b1(tempd); if(fStandardConfig->Get("neugen:kno_b2",tempd)) fStdNeuConfig->set_kno_b2(tempd); if(fStandardConfig->Get("neugen:kno_b3",tempd)) fStdNeuConfig->set_kno_b3(tempd); if(fStandardConfig->Get("neugen:kno_b4",tempd)) fStdNeuConfig->set_kno_b4(tempd); if(fStandardConfig->Get("neugen:kno_c1",tempd)) fStdNeuConfig->set_kno_c1(tempd); if(fStandardConfig->Get("neugen:kno_c2",tempd)) fStdNeuConfig->set_kno_c2(tempd); if(fStandardConfig->Get("neugen:kno_c3",tempd)) fStdNeuConfig->set_kno_c3(tempd); if(fStandardConfig->Get("neugen:kno_c4",tempd)) fStdNeuConfig->set_kno_c4(tempd); //cc if(fStandardConfig->Get("neugen:kno_r112",tempd)) fStdNeuConfig->set_dis_res(1,2,e_vp,tempd); if(fStandardConfig->Get("neugen:kno_r122",tempd)) fStdNeuConfig->set_dis_res(1,2,e_vn,tempd); if(fStandardConfig->Get("neugen:kno_r132",tempd)) fStdNeuConfig->set_dis_res(1,2,e_vbp,tempd); if(fStandardConfig->Get("neugen:kno_r142",tempd)) fStdNeuConfig->set_dis_res(1,2,e_vbn,tempd); if(fStandardConfig->Get("neugen:kno_r113",tempd)) fStdNeuConfig->set_dis_res(1,3,e_vp,tempd); if(fStandardConfig->Get("neugen:kno_r123",tempd)) fStdNeuConfig->set_dis_res(1,3,e_vn,tempd); if(fStandardConfig->Get("neugen:kno_r133",tempd)) fStdNeuConfig->set_dis_res(1,3,e_vbp,tempd); if(fStandardConfig->Get("neugen:kno_r143",tempd)) fStdNeuConfig->set_dis_res(1,3,e_vbn,tempd); //nc if(fStandardConfig->Get("neugen:kno_r212",tempd)) fStdNeuConfig->set_dis_res(2,2,e_vp,tempd); if(fStandardConfig->Get("neugen:kno_r222",tempd)) fStdNeuConfig->set_dis_res(2,2,e_vn,tempd); if(fStandardConfig->Get("neugen:kno_r232",tempd)) fStdNeuConfig->set_dis_res(2,2,e_vbp,tempd); if(fStandardConfig->Get("neugen:kno_r242",tempd)) fStdNeuConfig->set_dis_res(2,2,e_vbn,tempd); if(fStandardConfig->Get("neugen:kno_r213",tempd)) fStdNeuConfig->set_dis_res(2,3,e_vp,tempd); if(fStandardConfig->Get("neugen:kno_r223",tempd)) fStdNeuConfig->set_dis_res(2,3,e_vn,tempd); if(fStandardConfig->Get("neugen:kno_r233",tempd)) fStdNeuConfig->set_dis_res(2,3,e_vbp,tempd); if(fStandardConfig->Get("neugen:kno_r243",tempd)) fStdNeuConfig->set_dis_res(2,3,e_vbn,tempd); if(fStandardConfig->Get("neugen:wcutd",tempd)) fStdNeuConfig->set_wcutd(tempd); if(fStandardConfig->Get("neugen:wcutr",tempd)) fStdNeuConfig->set_wcutr(tempd); if(fStandardConfig->Get("neugen:nres",tempi)) fStdNeuConfig->set_nres(tempi); if(fStandardConfig->Get("neugen:dcf",tempd)) fStdNeuConfig->set_dcf(tempd); fStandardConfigChanged = false; } if(fReweightConfigChanged) { fUpdateConfig = true; //toggle using input values as absolute or scale factors: Float_t sc = 0; tempi = 0; if(fReweightConfig->Get("neugen:use_scale_factors",tempi)){ if(tempi==1) sc = 1.0; } //added by Tricia, 12-16-05 //first, set the name and number of the configuration we want to reweight to //if there is no name or number set, get them from the standard configuration bool gotname=false; bool gotno=false; if(fReweightConfig->Get("neugen:config_name",temps)){ fRwtNeuConfig->set_config_name(temps); gotname=true; } if(fReweightConfig->Get("neugen:config_no",tempi)){ fRwtNeuConfig->set_config_no(tempi); gotno=true; } if(!gotname||!gotno){ fRwtNeuConfig->set_config_name(fStdNeuConfig->config_name()); fRwtNeuConfig->set_config_no(fStdNeuConfig->config_no()); } //now, fill the neugen_config object with the default values of this //named configuration fRwtNeuConfig->fill_named_configuration(); //set values: //but, fill them based on the values that we just filled fRwtNeuConfig with //added by Tricia 12-16-05 if(fReweightConfig->Get("neugen:pdf_group",tempi)) fRwtNeuConfig->set_pdf_group(tempi); if(fReweightConfig->Get("neugen:pdf_set",tempi)) fRwtNeuConfig->set_pdf_set(tempi); if(fReweightConfig->Get("neugen:ma_qe",tempd)) fRwtNeuConfig->set_ma_qe(tempd*(1+sc*(fRwtNeuConfig->ma_qe()-1))); if(fReweightConfig->Get("neugen:ma_res",tempd)) fRwtNeuConfig->set_ma_res(tempd*(1+sc*(fRwtNeuConfig->ma_res()-1))); if(fReweightConfig->Get("neugen:ma_coh",tempd)) fRwtNeuConfig->set_ma_coh(tempd*(1+sc*(fRwtNeuConfig->ma_coh()-1))); if(fReweightConfig->Get("neugen:qel_fa0",tempd)) fRwtNeuConfig->set_qel_fa0(tempd*(1+sc*(fRwtNeuConfig->qel_fa0()-1))); if(fReweightConfig->Get("neugen:qel_eta",tempd)) fRwtNeuConfig->set_qel_eta(tempd*(1+sc*(fRwtNeuConfig->qel_eta()-1))); if(fReweightConfig->Get("neugen:res_omega",tempd)) fRwtNeuConfig->set_res_omega(tempd*(1+sc*(fRwtNeuConfig->res_omega()-1))); if(fReweightConfig->Get("neugen:res_z",tempd)) fRwtNeuConfig->set_res_z(tempd*(1+sc*(fRwtNeuConfig->res_z()-1))); if(fReweightConfig->Get("neugen:coh_r0",tempd)) fRwtNeuConfig->set_coh_r0(tempd*(1+sc*(fRwtNeuConfig->coh_r0()-1))); if(fReweightConfig->Get("neugen:coh_rei",tempd)) fRwtNeuConfig->set_coh_rei(tempd*(1+sc*(fRwtNeuConfig->coh_rei()-1))); if(fReweightConfig->Get("neugen:kno_a1",tempd)) fRwtNeuConfig->set_kno_a1(tempd*(1+sc*(fRwtNeuConfig->kno_a1()-1))); if(fReweightConfig->Get("neugen:kno_a2",tempd)) fRwtNeuConfig->set_kno_a2(tempd*(1+sc*(fRwtNeuConfig->kno_a2()-1))); if(fReweightConfig->Get("neugen:kno_a3",tempd)) fRwtNeuConfig->set_kno_a3(tempd*(1+sc*(fRwtNeuConfig->kno_a3()-1))); if(fReweightConfig->Get("neugen:kno_a4",tempd)) fRwtNeuConfig->set_kno_a4(tempd*(1+sc*(fRwtNeuConfig->kno_a4()-1))); if(fReweightConfig->Get("neugen:kno_b1",tempd)) fRwtNeuConfig->set_kno_b1(tempd*(1+sc*(fRwtNeuConfig->kno_b1()-1))); if(fReweightConfig->Get("neugen:kno_b2",tempd)) fRwtNeuConfig->set_kno_b2(tempd*(1+sc*(fRwtNeuConfig->kno_b2()-1))); if(fReweightConfig->Get("neugen:kno_b3",tempd)) fRwtNeuConfig->set_kno_b3(tempd*(1+sc*(fRwtNeuConfig->kno_b3()-1))); if(fReweightConfig->Get("neugen:kno_b4",tempd)) fRwtNeuConfig->set_kno_b4(tempd*(1+sc*(fRwtNeuConfig->kno_b4()-1))); if(fReweightConfig->Get("neugen:kno_c1",tempd)) fRwtNeuConfig->set_kno_c1(tempd*(1+sc*(fRwtNeuConfig->kno_c1()-1))); if(fReweightConfig->Get("neugen:kno_c2",tempd)) fRwtNeuConfig->set_kno_c2(tempd*(1+sc*(fRwtNeuConfig->kno_c2()-1))); if(fReweightConfig->Get("neugen:kno_c3",tempd)) fRwtNeuConfig->set_kno_c3(tempd*(1+sc*(fRwtNeuConfig->kno_c3()-1))); if(fReweightConfig->Get("neugen:kno_c4",tempd)) fRwtNeuConfig->set_kno_c4(tempd*(1+sc*(fRwtNeuConfig->kno_c4()-1))); //kno cc if(fReweightConfig->Get("neugen:kno_r112",tempd)) fRwtNeuConfig->set_dis_res(1,2,e_vp,tempd*(1+sc*(fRwtNeuConfig->dis_res(1,2,e_vp)-1))); if(fReweightConfig->Get("neugen:kno_r122",tempd)) fRwtNeuConfig->set_dis_res(1,2,e_vn,tempd*(1+sc*(fRwtNeuConfig->dis_res(1,2,e_vn)-1))); if(fReweightConfig->Get("neugen:kno_r132",tempd)) fRwtNeuConfig->set_dis_res(1,2,e_vbp,tempd*(1+sc*(fRwtNeuConfig->dis_res(1,2,e_vbp)-1))); if(fReweightConfig->Get("neugen:kno_r142",tempd)) fRwtNeuConfig->set_dis_res(1,2,e_vbn,tempd*(1+sc*(fRwtNeuConfig->dis_res(1,2,e_vbn)-1))); if(fReweightConfig->Get("neugen:kno_r113",tempd)) fRwtNeuConfig->set_dis_res(1,3,e_vp,tempd*(1+sc*(fRwtNeuConfig->dis_res(1,3,e_vp)-1))); if(fReweightConfig->Get("neugen:kno_r123",tempd)) fRwtNeuConfig->set_dis_res(1,3,e_vn,tempd*(1+sc*(fRwtNeuConfig->dis_res(1,3,e_vn)-1))); if(fReweightConfig->Get("neugen:kno_r133",tempd)) fRwtNeuConfig->set_dis_res(1,3,e_vbp,tempd*(1+sc*(fRwtNeuConfig->dis_res(1,3,e_vbp)-1))); if(fReweightConfig->Get("neugen:kno_r143",tempd)) fRwtNeuConfig->set_dis_res(1,3,e_vbn,tempd*(1+sc*(fRwtNeuConfig->dis_res(1,3,e_vbn)-1))); //kno nc if(fReweightConfig->Get("neugen:kno_r212",tempd)) fRwtNeuConfig->set_dis_res(2,2,e_vp,tempd*(1+sc*(fRwtNeuConfig->dis_res(2,2,e_vp)-1))); if(fReweightConfig->Get("neugen:kno_r222",tempd)) fRwtNeuConfig->set_dis_res(2,2,e_vn,tempd*(1+sc*(fRwtNeuConfig->dis_res(2,2,e_vn)-1))); if(fReweightConfig->Get("neugen:kno_r232",tempd)) fRwtNeuConfig->set_dis_res(2,2,e_vbp,tempd*(1+sc*(fRwtNeuConfig->dis_res(2,2,e_vbp)-1))); if(fReweightConfig->Get("neugen:kno_r242",tempd)) fRwtNeuConfig->set_dis_res(2,2,e_vbn,tempd*(1+sc*(fRwtNeuConfig->dis_res(2,2,e_vbn)-1))); if(fReweightConfig->Get("neugen:kno_r213",tempd)) fRwtNeuConfig->set_dis_res(2,3,e_vp,tempd*(1+sc*(fRwtNeuConfig->dis_res(2,3,e_vp)-1))); if(fReweightConfig->Get("neugen:kno_r223",tempd)) fRwtNeuConfig->set_dis_res(2,3,e_vn,tempd*(1+sc*(fRwtNeuConfig->dis_res(2,3,e_vn)-1))); if(fReweightConfig->Get("neugen:kno_r233",tempd)) fRwtNeuConfig->set_dis_res(2,3,e_vbp,tempd*(1+sc*(fRwtNeuConfig->dis_res(2,3,e_vbp)-1))); if(fReweightConfig->Get("neugen:kno_r243",tempd)) fRwtNeuConfig->set_dis_res(2,3,e_vbn,tempd*(1+sc*(fRwtNeuConfig->dis_res(2,3,e_vbn)-1))); //kno scales: if(fReweightConfig->Get("neugen:scale_kno_all",tempd)){ for(int c=1;c<=2;c++){ for(int m=2;m<=3;m++){ fRwtNeuConfig->set_dis_res(c,m,e_vp, fRwtNeuConfig->dis_res(c,m,e_vp)*tempd); fRwtNeuConfig->set_dis_res(c,m,e_vn, fRwtNeuConfig->dis_res(c,m,e_vn)*tempd); fRwtNeuConfig->set_dis_res(c,m,e_vbp, fRwtNeuConfig->dis_res(c,m,e_vbp)*tempd); fRwtNeuConfig->set_dis_res(c,m,e_vbn, fRwtNeuConfig->dis_res(c,m,e_vbn)*tempd); } } } if(fReweightConfig->Get("neugen:scale_kno_cc",tempd)){ for(int m=2;m<=3;m++){ fRwtNeuConfig->set_dis_res(1,m,e_vp, fRwtNeuConfig->dis_res(1,m,e_vp)*tempd); fRwtNeuConfig->set_dis_res(1,m,e_vn, fRwtNeuConfig->dis_res(1,m,e_vn)*tempd); fRwtNeuConfig->set_dis_res(1,m,e_vbp, fRwtNeuConfig->dis_res(1,m,e_vbp)*tempd); fRwtNeuConfig->set_dis_res(1,m,e_vbn, fRwtNeuConfig->dis_res(1,m,e_vbn)*tempd); } } if(fReweightConfig->Get("neugen:scale_kno_nc",tempd)){ for(int m=2;m<=3;m++){ fRwtNeuConfig->set_dis_res(2,m,e_vp, fRwtNeuConfig->dis_res(2,m,e_vp)*tempd); fRwtNeuConfig->set_dis_res(2,m,e_vn, fRwtNeuConfig->dis_res(2,m,e_vn)*tempd); fRwtNeuConfig->set_dis_res(2,m,e_vbp, fRwtNeuConfig->dis_res(2,m,e_vbp)*tempd); fRwtNeuConfig->set_dis_res(2,m,e_vbn, fRwtNeuConfig->dis_res(2,m,e_vbn)*tempd); } } if(fReweightConfig->Get("neugen:wcutd",tempd)) fRwtNeuConfig->set_wcutd(tempd*(1+sc*(fRwtNeuConfig->wcutd()-1))); if(fReweightConfig->Get("neugen:wcutr",tempd)) fRwtNeuConfig->set_wcutr(tempd*(1+sc*(fRwtNeuConfig->wcutr()-1))); if(fReweightConfig->Get("neugen:nres",tempi)) fRwtNeuConfig->set_nres(tempi); if(fReweightConfig->Get("neugen:dcf",tempd)) fRwtNeuConfig->set_dcf(tempd*(1+sc*(fRwtNeuConfig->dcf()-1))); fReweightConfigChanged = false; } }

double NeugenWeightCalculator::GetWeight (  MCEventInfo *  event, NuParent *  parent )  [virtual]

Reimplemented from WeightCalculator. Definition at line 531 of file NeugenWeightCalculator.cxx. References abs(), ccnc_enum, flavor_enum, fNuWrap, fRwtNeuConfig, fUpdateConfig, MCEventInfo::had_fs, MCEventInfo::iaction, init_state_enum, MCEventInfo::initial_state, MCEventInfo::inu, MCEventInfo::iresonance, kinematic_variable_enum, MSG, MCEventInfo::nucleus, nucleus_enum, MCEventInfo::nuE, MCEventInfo::nuPx, MCEventInfo::nuPy, MCEventInfo::nuPz, neugen_wrapper::offshell_diff_xsec(), process_enum, MCEventInfo::q2, neugen_wrapper::reweight(), MCEventInfo::stdXSec, MCEventInfo::tarE, MCEventInfo::tarPx, MCEventInfo::tarPy, MCEventInfo::tarPz, MCEventInfo::useStdXSec, MCEventInfo::w2, MCEventInfo::x, and MCEventInfo::y. { if(!event) return 1; if(!parent) { MSG("NeugenWC",Msg::kDebug) << "No NuParent object" << endl; } //event neutrino energy, momentum double nu_en = event->nuE; double nu_px = event->nuPx; double nu_py = event->nuPy; double nu_pz = event->nuPz; if(nu_en==0) return 1; //event target energy, momentum double tar_en = event->tarE; double tar_px = event->tarPx; double tar_py = event->tarPy; double tar_pz = event->tarPz; if(tar_en==0) return 1; //event kinematic variables: double kval[5] = {0,0,0,0,0}; int kid1 = -1; int kid2 = -1; kval[1] = event->q2; kval[2] = event->w2; kval[3] = event->x; kval[4] = event->y; //event type int iaction = event->iaction; int inu = event->inu; int ires = event->iresonance; int initial_state = event->initial_state; int nucleus = event->nucleus; int had_fs = event->had_fs; //translate from minossoft-speak to neugen-speak: int flavor = 4; if(abs(inu)==12) flavor = 1; else if(abs(inu)==14) flavor = 2; else if(abs(inu)==16) flavor = 3; else if(abs(inu)==1 || abs(inu)==2 || abs(inu)==3) flavor = abs(inu); else { MSG("NeugenWC",Msg::kError) << "Invalid event:inu value " << inu << endl; return 1; } int ccnc = iaction; if(ccnc==0) ccnc = 2; if(!(ccnc==1||ccnc==2)) { MSG("NeugenWC",Msg::kError) << "Invalid event:iaction value " << iaction << endl; return 1; } int process = 0; if(ires>=1001&&ires<=1004) process = ires - 1000; else if(ires>=1&&ires<=4) process = ires; else { if(ires==1005) MSG("NeugenWC",Msg::kInfo) << "Event:iresonance value is 1005 " << "- no reweighting performed" << endl; else MSG("NeugenWC",Msg::kError) << "Invalid event:iresonance (or event:process) value " << ires << endl; return 1; } //extract necessary kinematic variables: if(process==1) { //qel kid1=1; kid2=0; if(kval[1]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is QEL: q^2 kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } } else if(process==2) { //res kid1=1; kid2=2; if(kval[1]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is RES: q^2 kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } if(kval[2]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is RES: w^2 kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } //need to pass w rather than w^2 to neugen if(kval[2]>0) kval[2]=sqrt(kval[2]); else { MSG("NeugenWC",Msg::kWarning) << "w^2 kinematic variable is negative!" << " Returning weight of 1." << endl; return 1; } } else if(process==3) { //dis kid1=3; kid2=4; if(kval[3]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is DIS: x kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } if(kval[4]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is DIS: y kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } } else if(process==4) { //coh //reweighting not possible here yet return 1; } //start making neugen style objects: TLorentzVector *nu_lv = new TLorentzVector(nu_px,nu_py,nu_pz,nu_en); TLorentzVector *tar_lv = new TLorentzVector(tar_px,tar_py,tar_pz,tar_en); interaction *iact = new interaction(flavor_enum(flavor), nucleus_enum(nucleus), ccnc_enum(ccnc), init_state_enum(initial_state)); final_state *fs = new final_state(had_fs); //calculate weight double weight = 0; if(event->useStdXSec) { //calculate standard xsec and store: if(event->stdXSec==0) event->stdXSec = fNuWrap->offshell_diff_xsec(nu_lv,tar_lv, iact,process_enum(process),fs, kinematic_variable_enum(kid1),float(kval[kid1]), kinematic_variable_enum(kid2),float(kval[kid2])); if(!fIsOptimized){ weight = fNuWrap->reweight(nu_lv,tar_lv, kinematic_variable_enum(kid1),float(kval[kid1]), kinematic_variable_enum(kid2),float(kval[kid2]), iact,process_enum(process),fs, fRwtNeuConfig,event->stdXSec, fUpdateConfig); fUpdateConfig = false; }else{ //calculate reweight factor using stdXSec as denominator weight = fNuWrap->reweight(nu_lv,tar_lv, kinematic_variable_enum(kid1),float(kval[kid1]), kinematic_variable_enum(kid2),float(kval[kid2]), iact,process_enum(process),fs, fRwtNeuConfig,0,event->stdXSec); } } else { weight = fNuWrap->reweight(nu_lv,tar_lv, kinematic_variable_enum(kid1),float(kval[kid1]), kinematic_variable_enum(kid2),float(kval[kid2]), iact,process_enum(process),fs, fRwtNeuConfig); } delete nu_lv; delete tar_lv; delete iact; delete fs; return weight; }

double NeugenWeightCalculator::GetWeight (  Registry *  eventinfo  )  [virtual]

Reimplemented from WeightCalculator. Definition at line 369 of file NeugenWeightCalculator.cxx. References abs(), ccnc_enum, flavor_enum, fNuWrap, fRwtNeuConfig, Registry::Get(), init_state_enum, kinematic_variable_enum, MSG, nucleus_enum, process_enum, and neugen_wrapper::reweight(). { if(!event) return 1; //event neutrino energy, momentum double nu_en = 0; double nu_px = 0; double nu_py = 0; double nu_pz = 0; if(!event->Get("event:nu_en",nu_en)) return 1; if(!event->Get("event:nu_px",nu_px)) return 1; if(!event->Get("event:nu_py",nu_py)) return 1; if(!event->Get("event:nu_pz",nu_pz)) return 1; //event target energy, momentum double tar_en = 0; double tar_px = 0; double tar_py = 0; double tar_pz = 0; if(!event->Get("event:tar_en",tar_en)) return 1; if(!event->Get("event:tar_px",tar_px)) return 1; if(!event->Get("event:tar_py",tar_py)) return 1; if(!event->Get("event:tar_pz",tar_pz)) return 1; //event kinematic variables: double kval[5] = {0,0,0,0,0}; int kid1 = -1; int kid2 = -1; if(!event->Get("event:q2",kval[1])) kval[1] = 999999; if(!event->Get("event:w2",kval[2])) kval[2] = 999999; if(!event->Get("event:x",kval[3])) kval[3] = 999999; if(!event->Get("event:y",kval[4])) kval[4] = 999999; //event type int iaction = 0; int inu = 0; int ires = 0; int initial_state = 0; int nucleus = 0; int had_fs = 0; if(!event->Get("event:iaction",iaction)) return 1; if(!event->Get("event:inu",inu)) return 1; if(!event->Get("event:iresonance",ires)) { if(!event->Get("event:process",ires)) return 1; } if(!event->Get("event:initial_state",initial_state)) return 1; if(!event->Get("event:nucleus",nucleus)) return 1; if(!event->Get("event:had_fs",had_fs)) return 1; //translate from minossoft-speak to neugen-speak: int flavor = 4; if(abs(inu)==12) flavor = 1; else if(abs(inu)==14) flavor = 2; else if(abs(inu)==16) flavor = 3; else if(abs(inu)==1 || abs(inu)==2 || abs(inu)==3) flavor = abs(inu); else { MSG("NeugenWC",Msg::kError) << "Invalid event:inu value " << inu << endl; return 1; } int ccnc = iaction; if(ccnc==0) ccnc = 2; if(!(ccnc==1||ccnc==2)) { MSG("NeugenWC",Msg::kError) << "Invalid event:iaction value " << iaction << endl; return 1; } int process = 0; if(ires>=1001&&ires<=1004) process = ires - 1000; else if(ires>=1&&ires<=4) process = ires; else { if(ires==1005) MSG("NeugenWC",Msg::kInfo) << "Event:iresonance value is 1005 " << "- no reweighting performed" << endl; else MSG("NeugenWC",Msg::kError) << "Invalid event:iresonance (or event:process) value " << ires << endl; return 1; } //extract necessary kinematic variables: if(process==1) { //qel kid1=1; kid2=0; if(kval[1]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is QEL: q^2 kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } } else if(process==2) { //res kid1=1; kid2=2; if(kval[1]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is RES: q^2 kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } if(kval[2]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is RES: w^2 kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } //need to pass w rather than w^2 to neugen if(kval[2]>0) kval[2]=sqrt(kval[2]); else { MSG("NeugenWC",Msg::kWarning) << "w^2 kinematic variable is negative!" << " Returning weight of 1." << endl; return 1; } } else if(process==3) { //dis kid1=3; kid2=4; if(kval[3]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is DIS: x kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } if(kval[4]==999999) { MSG("NeugenWC",Msg::kWarning) << "Current event is DIS: y kinematic variable not set!" << " Returning weight of 1." << endl; return 1; } } else if(process==4) { //coh //reweighting not possible here yet return 1; } //start making neugen style objects: TLorentzVector *nu_lv = new TLorentzVector(nu_px,nu_py,nu_pz,nu_en); TLorentzVector *tar_lv = new TLorentzVector(tar_px,tar_py,tar_pz,tar_en); interaction *iact = new interaction(flavor_enum(flavor), nucleus_enum(nucleus), ccnc_enum(ccnc), init_state_enum(initial_state)); final_state *fs = new final_state(had_fs); //calculate weight double weight = 0; weight = fNuWrap->reweight(nu_lv,tar_lv, kinematic_variable_enum(kid1),float(kval[kid1]), kinematic_variable_enum(kid2),float(kval[kid2]), iact,process_enum(process),fs, fRwtNeuConfig); delete nu_lv; delete tar_lv; delete iact; delete fs; return weight; }

void NeugenWeightCalculator::PrintReweightConfig (  ostream &   )  [virtual]

Reimplemented from WeightCalculator. Definition at line 710 of file NeugenWeightCalculator.cxx. References fRwtNeuConfig, and neugen_config::print(). { if(fRwtNeuConfig) fRwtNeuConfig->print(stream); }

void NeugenWeightCalculator::ReweightConfigReset (   )  [virtual]

Reimplemented from WeightCalculator. Definition at line 51 of file NeugenWeightCalculator.cxx. References fRwtNeuConfig, neugen_config::set_best_parameters(), and WeightCalculator::SetReweightConfig(). { fRwtNeuConfig->set_best_parameters(); if(fStandardConfig) SetReweightConfig(fStandardConfig); }

void NeugenWeightCalculator::SetOptimization (  bool  in  )

Definition at line 59 of file NeugenWeightCalculator.cxx. References fIsOptimized. Referenced by SKZPWeightCalculator::Config(). { fIsOptimized = in; }

---

Member Data Documentation

bool NeugenWeightCalculator::fIsOptimized [private]

Definition at line 31 of file NeugenWeightCalculator.h. Referenced by NeugenWeightCalculator(), and SetOptimization().

neugen_wrapper* NeugenWeightCalculator::fNuWrap [protected]

Definition at line 15 of file NeugenWeightCalculator.h. Referenced by GetWeight(), and NeugenWeightCalculator().

neugen_config* NeugenWeightCalculator::fRwtNeuConfig [protected]

Definition at line 14 of file NeugenWeightCalculator.h. Referenced by Config(), GetWeight(), NeugenWeightCalculator(), PrintReweightConfig(), and ReweightConfigReset().

neugen_config* NeugenWeightCalculator::fStdNeuConfig [protected]

Definition at line 13 of file NeugenWeightCalculator.h. Referenced by Config(), and NeugenWeightCalculator().

bool NeugenWeightCalculator::fUpdateConfig [private]

Definition at line 32 of file NeugenWeightCalculator.h. Referenced by Config(), GetWeight(), and NeugenWeightCalculator().

---

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

• NeugenWeightCalculator.h
• NeugenWeightCalculator.cxx

---