NuTransition Class Reference

#include <NuTransition.h>

List of all members.

Public Member Functions
	NuTransition (const NuXMLConfig *xmlConfig)
	~NuTransition ()
void	Fill (NuMCEvent &mc)
void	Reweight (NuEvent &nu)
NuEvent *	GetTransitionEvent (const NuEvent &nu)
Bool_t	IsSetToDoTransitions () const
Private Member Functions
Float_t	Baseline () const
Float_t	OscillationWeight (const Float_t energy, const Int_t inu, const NuXMLConfig *lxmlConfig=0) const
void	DoSystematicShifts (NuEvent &event)
Double_t	GetXSecGraph (NuEvent &nu, bool opp_charge=false) const
Double_t	GetXSecGraph (NuMCEvent &nu, bool opp_charge=false) const
Double_t	GetXSecNeugen (NuEvent &nu, bool opp_charge=false) const
Private Attributes
TH1D *	numu_flux
TH1D *	nubar_flux
TH1D *	flux_weight
TGraph *	fNuMuXGraph
TGraph *	fNuBarXGraph
NuXMLConfig *	fxmlConfig
NuSystematic *	fSyst
bool	applyTranSyst
bool	firstRunTrans
bool	doNothing

---

Detailed Description

Definition at line 18 of file NuTransition.h.

---

Constructor & Destructor Documentation

NuTransition::NuTransition

(

const NuXMLConfig *

xmlConfig

)

Definition at line 46 of file NuTransition.cxx.

References applyTranSyst, NuXMLConfig::BinningScheme(), NuXMLConfig::DM2Bar(), NuXMLConfig::DM2Nu(), doNothing, firstRunTrans, flux_weight, NuSystematic::FluxSyst(), fNuBarXGraph, fNuMuXGraph, fSyst, fxmlConfig, Msg::kError, Msg::kInfo, MAXMSG, MSG, nubar_flux, numu_flux, NuUtilities::RecoBins(), NuXMLConfig::SN2Bar(), NuXMLConfig::SN2Nu(), NuXMLConfig::TransitionProb(), and NuUtilities::TrueBins().

{
    // Get our local xmlConfig copy, if it is safe to do so
    if (!xmlConfig) {
        // If the XML file does not exist, we can't do anything
        doNothing = true;
        
        MAXMSG("NuTransition",Msg::kInfo,5)
        << "No NuXMLConfig object supplied: "
        << "performing no systematic shifts, oscillations, or transitions."
        << endl;
        return;
    }
    
    fxmlConfig = (NuXMLConfig*)xmlConfig->Clone();
    firstRunTrans = true;

    MSG("NuTransition",Msg::kInfo) << "Applying parameters: " 
    << "dm2nu = " << xmlConfig->DM2Nu()   << ", sn2nu = " << xmlConfig->SN2Nu()
    << ", dm2bar = " << xmlConfig->DM2Bar() << ", sn2bar = " << xmlConfig->SN2Bar()
    << ", alpha = " << xmlConfig->TransitionProb() << endl;
    
    doNothing  = false;
    
    if (doNothing) return;
    
    fSyst = new NuSystematic(*fxmlConfig);
    
    // Get the binning scheme set up
    const NuUtilities cuts;
  
    //Binning scheme
    NuBinningScheme::NuBinningScheme_t binningScheme =
        static_cast<NuBinningScheme::NuBinningScheme_t>(xmlConfig->BinningScheme());
  
    std::vector<Double_t> vReco = cuts.RecoBins(binningScheme);
    std::vector<Double_t> vTrue = cuts.TrueBins(binningScheme);
    
    Int_t numRecoBins = vReco.size() - 1;
    Int_t numTrueBins = vTrue.size() - 1;
    
    Double_t *recoBinsArray = &(vReco[0]);
    Double_t *trueBinsArray = &(vTrue[0]);
    
    // Creat the reweighting histograms
    numu_flux   = new TH1D("numu_flux","numu_flux",numTrueBins,trueBinsArray);
    nubar_flux  = new TH1D("nubar_flux","nubar_flux",numRecoBins,recoBinsArray);
    flux_weight = new TH1D("flux_weight","flux_weight",numRecoBins,recoBinsArray);
    
    // Determine how the systematics are applied
    if (NuSystematic::FluxSyst(*xmlConfig)) { 
        MSG("NuTransition", Msg::kInfo) << "Not applying " << xmlConfig->Name() 
        << " to transitions b/c it's a flux systematic." << endl;
        applyTranSyst = false;
    }
    else {
        applyTranSyst = true;
    }
    
    // Get the XSec graphs from file
    TDirectory *start = gDirectory;
    
    MSG("NuTransition",Msg::kInfo) << "Get cross-section graphs from file." << endl;
    
    const char * fxsecfilename = getenv("xf");
    
    MSG("NuTransition",Msg::kInfo) << "Opening xsec file: " << fxsecfilename << endl;
    TFile *xsecfile = new TFile(fxsecfilename,"READ");
  
    // Check to see if this failed
    if (!xsecfile->IsOpen()) {
      MSG("NuTransition", Msg::kError) << "\n\n"
        << "*****************************************************************************\n"
        << "Failed to open cross section file.\n"
        << "Currently the cross section file is found through the 'xf' environment variable.\n"
        << "Perhaps you need to issue a command such as:\n   "
        << "   export xf=\"$SRT_PRIVATE_CONTEXT/NtupleUtils/data/xsec_minos_modbyrs4_v3_5_0_mk.root\"\n"
        << "*****************************************************************************\n\n"
        << endl;
      throw runtime_error("Failed to read crosssection file");
    }
    
    // NuMu
    TH1F *fhNuMuCCCrossSections = (TH1F*) xsecfile->Get("h_numu_cc_tot");
    Float_t *x = new Float_t[fhNuMuCCCrossSections->GetNbinsX()];
    Float_t *y = new Float_t[fhNuMuCCCrossSections->GetNbinsX()];
    for(int i=0;i<fhNuMuCCCrossSections->GetNbinsX();i++) {
        x[i] = fhNuMuCCCrossSections->GetBinCenter(i+1);
        y[i] = fhNuMuCCCrossSections->GetBinContent(i+1);
    }
    fNuMuXGraph = new TGraph(fhNuMuCCCrossSections->GetNbinsX(),x,y);
    if (x) {delete[] x; x = 0;}
    if (y) {delete[] y; y = 0;}
    
    // NuMuBar
    TH1F *fhNuBarCCCrossSections = (TH1F*) xsecfile->Get("h_numubar_cc_tot");
    x = new Float_t[fhNuBarCCCrossSections->GetNbinsX()];
    y = new Float_t[fhNuBarCCCrossSections->GetNbinsX()];
    for(int i=0;i<fhNuBarCCCrossSections->GetNbinsX();i++) {
        x[i] = fhNuBarCCCrossSections->GetBinCenter(i+1);
        y[i] = fhNuBarCCCrossSections->GetBinContent(i+1);
    }
    fNuBarXGraph = new TGraph(fhNuBarCCCrossSections->GetNbinsX(),x,y);
    if (x) {delete[] x; x = 0;}
    if (y) {delete[] y; y = 0;}
    
    xsecfile->Close();
    if (xsecfile){delete xsecfile; xsecfile = 0;}
    
    gDirectory = start;

}

NuTransition::~NuTransition

(

)

Definition at line 160 of file NuTransition.cxx.

References flux_weight, nubar_flux, and numu_flux.

{
    if (numu_flux) delete numu_flux; numu_flux = 0;
    if (nubar_flux) delete nubar_flux; nubar_flux = 0;
    if (flux_weight) delete flux_weight; flux_weight = 0;
}

---

Member Function Documentation

Float_t NuTransition::Baseline

(

)

const [inline, private]

Definition at line 39 of file NuTransition.h.

{return 735.0;};

void NuTransition::DoSystematicShifts

(

NuEvent &

event

)

[private]

Definition at line 331 of file NuTransition.cxx.

References fSyst, NuXMLConfig::FullTitle(), fxmlConfig, Msg::kInfo, MAXMSG, and NuSystematic::Shift().

Referenced by Reweight().

{
    MAXMSG("NuTransition",Msg::kInfo,5)
    << "Performing systematic shift "
    << fxmlConfig->FullTitle()
    << endl;
    fSyst->Shift(event);
}

void NuTransition::Fill

(

NuMCEvent &

mc

)

Definition at line 168 of file NuTransition.cxx.

References doNothing, GetXSecGraph(), NuMCEvent::iaction, NuMCEvent::inu, Msg::kInfo, MAXMSG, NuMCEvent::neuEnMC, nubar_flux, numu_flux, and NuMCEvent::rw.

{
    if (doNothing) return;
    
    // CC-only
    if (mc.iaction==1) { 
        double fw = 1;
        
        // Correct for Cross section
        fw *= this->GetXSecGraph(mc);
        
        if (mc.inu == 14) {
            numu_flux->Fill(mc.neuEnMC, mc.rw/fw);
            MAXMSG("NuTransition",Msg::kInfo,5) << "===> Fill numu, xsec=" << fw << endl;
        }
        else if (mc.inu == -14) {
            nubar_flux->Fill(mc.neuEnMC, mc.rw/fw);
            MAXMSG("NuTransition",Msg::kInfo,5) << "===> Fill nubar, xsec=" << fw << endl;
        }
        else{
            MAXMSG("NuTransition",Msg::kInfo,5) << "===> Fill nothing, ntyp=" << mc.inu << endl;
        }
    }
}

NuEvent* NuTransition::GetTransitionEvent

(

const NuEvent &

nu

)

Double_t NuTransition::GetXSecGraph	(	NuMCEvent &	nu,
		bool	opp_charge = false
	)			const [private]

Definition at line 387 of file NuTransition.cxx.

References fNuBarXGraph, fNuMuXGraph, NuMCEvent::inu, NuMCEvent::neuEnMC, and NueConvention::numu.

{
    Double_t nuBarXSec = fNuBarXGraph->Eval(nu.neuEnMC,0,"");
    Double_t nuMuXSec =  fNuMuXGraph->Eval( nu.neuEnMC,0,"");
    
    bool numu = (nu.inu > 0);
    if (opp_charge) numu = !numu;
       
    if (numu) return nuMuXSec;
    else      return nuBarXSec;
}

Double_t NuTransition::GetXSecGraph	(	NuEvent &	nu,
		bool	opp_charge = false
	)			const [private]

Definition at line 374 of file NuTransition.cxx.

References fNuBarXGraph, fNuMuXGraph, NuEvent::inu, NuEvent::neuEnMC, and NueConvention::numu.

Referenced by Fill().

{
    Double_t nuBarXSec = fNuBarXGraph->Eval(nu.neuEnMC,0,"");
    Double_t nuMuXSec =  fNuMuXGraph->Eval( nu.neuEnMC,0,"");
    
    bool numu = (nu.inu > 0);
    if (opp_charge) numu = !numu;
    
    if (numu) return nuMuXSec;
    else      return nuBarXSec;
}

Double_t NuTransition::GetXSecNeugen	(	NuEvent &	nu,
		bool	opp_charge = false
	)			const [private]

Definition at line 401 of file NuTransition.cxx.

References e_undefined_init_state, e_vA, e_vbA, e_vbn, e_vbN, e_vbp, e_vn, e_vN, e_vp, NuEvent::hadronicFinalStateMC, NuEvent::iaction, init(), NuEvent::initialStateMC, NuEvent::inu, NuEvent::iresonance, Msg::kError, Msg::kInfo, Msg::kWarning, MSG, NuEvent::neuEnMC, NuEvent::neuPxMC, NuEvent::neuPyMC, NuEvent::neuPzMC, NuEvent::nucleusMC, neugen_wrapper::offshell_diff_xsec(), NuEvent::q2MC, neugen_config::set_best_parameters(), NuEvent::tgtEnMC, NuEvent::tgtPxMC, NuEvent::tgtPyMC, NuEvent::tgtPzMC, NuEvent::w2MC, NuEvent::xMC, and NuEvent::yMC.

{
    // Get Cross Sections //
    
    static neugen_wrapper *fNuWrap = 0;
    //static Registry std_registry;   
    static neugen_config *std_config = 0;
    
    if (!fNuWrap) {
        //std_registry.Set("neugen:config_name","MODBYRS");
        //std_registry.Set("neugen:config_no",4);
        
        //NuSystematic::SetNeugenDefaults(std_registry);
        
        std_config = new neugen_config("stdNeuConfig");
        std_config->set_best_parameters();
        //if(config) SetStandardConfig(config);
        
        fNuWrap = new neugen_wrapper(std_config);
    }
    
    //event kinematic variables:
    double kval[5] = {0,0,0,0,0};
    int kid1 = -1;
    int kid2 = -1;
    kval[1] = nu.q2MC;
    kval[2] = nu.w2MC; 
    kval[3] = nu.xMC;
    kval[4] = nu.yMC;
    
    //translate from minossoft-speak to neugen-speak:
    int flavor = 4;
    if(abs(nu.inu)==12) flavor = 1;
    else if(abs(nu.inu)==14) flavor = 2; 
    else if(abs(nu.inu)==16) flavor = 3;
    else if(abs(nu.inu)==1 || abs(nu.inu)==2 || abs(nu.inu)==3) flavor = abs(nu.inu);
    else {
        MSG("NuTransition",Msg::kError) << "Invalid event:inu value " << nu.inu << endl;
        return 0;
    }
    
    int ccnc = nu.iaction;
    if(ccnc==0) ccnc = 2;
    if(!(ccnc==1||ccnc==2)) {
        MSG("NuTransition",Msg::kError) << "Invalid event:iaction value " 
        << nu.iaction << endl;
        return 0;
    }
    
    int process = 0;
    if(nu.iresonance>=1001&&nu.iresonance<=1004) process = nu.iresonance - 1000;
    else if(nu.iresonance>=1&&nu.iresonance<=4) process = nu.iresonance;
    else {
        if(nu.iresonance==1005) MSG("NuTransition",Msg::kInfo) 
            << "Event:iresonance value is 1005 "
            << "- no reweighting performed" << endl;
        else MSG("NeugenWC",Msg::kError) 
            << "Invalid event:iresonance (or event:process) value "
            << nu.iresonance << endl;
        return 0;
    }
    
    //extract necessary kinematic variables:
    if(process==1) {   //qel
        kid1=1;
        kid2=0;
        if(kval[1]==999999) {
            MSG("NuTransition",Msg::kWarning) 
            << "Current event is QEL: q^2 kinematic variable not set!"
            << " Returning weight of 1." << endl;
            return 0;
        }    
    }
    else if(process==2) { //res
        kid1=1;
        kid2=2;
        if(kval[1]==999999) {
            MSG("NuTransition",Msg::kWarning) 
            << "Current event is RES: q^2 kinematic variable not set!"
            << " Returning weight of 1." << endl;
            return 0;
        }
        if(kval[2]==999999) {
            MSG("NuTransition",Msg::kWarning) 
            << "Current event is RES: w^2 kinematic variable not set!"
            << " Returning weight of 1." << endl;
            return 0;
        }
        //need to pass w rather than w^2 to neugen
        if(kval[2]>0) kval[2]=sqrt(kval[2]); 
        else {
            MSG("NuTransition",Msg::kWarning) 
            << "w^2 kinematic variable is negative!"
            << " Returning weight of 1." << endl;
            return 0;
        }
    }
    else if(process==3) { //dis
        kid1=3;
        kid2=4;
        if(kval[3]==999999) {
            MSG("NuTransition",Msg::kWarning) 
            << "Current event is DIS: x kinematic variable not set!"
            << " Returning weight of 1." << endl;
            return 0;
        }
        if(kval[4]==999999) {
            MSG("NuTransition",Msg::kWarning) 
            << "Current event is DIS: y kinematic variable not set!"
            << " Returning weight of 1." << endl;
            return 0;
        }
    }
    else if(process==4) { //coh
        //reweighting not possible here yet
        return 0;
    }
    
    //start making neugen style objects:
    TLorentzVector *nu_lv = new TLorentzVector(nu.neuPxMC,nu.neuPyMC,nu.neuPzMC,nu.neuEnMC);
    TLorentzVector *tar_lv = new TLorentzVector(nu.tgtPxMC,nu.tgtPyMC,nu.tgtPzMC,nu.tgtEnMC);
    
    init_state_t init = init_state_enum(nu.initialStateMC);
    
    if (opp_charge) {
        using namespace init_state;
        switch(init) {
            case e_vN:   init = e_vbN;  break;
            case e_vbN:  init = e_vN;   break;
                
            case e_vA:   init = e_vbA;  break;
            case e_vbA:  init = e_vA;   break;
                
            case e_vp:   init = e_vbp;  break;
            case e_vbp:  init = e_vp;   break;
                
            case e_vn:   init = e_vbn;  break;
            case e_vbn:  init = e_vn;   break;
                
            case e_undefined_init_state:
            default:            
                MSG("NuTransition",Msg::kWarning) << "Unknown initial state" << endl;
                return 0;
        }
        
        if (process == 1) {
            if (init == e_vbn) init = e_vbp;
            if (init == e_vp)  init = e_vn;
        }  
    }
    
    interaction *iact = new interaction(flavor_enum(flavor),
                                        nucleus_enum(nu.nucleusMC),
                                        ccnc_enum(ccnc),
                                        init); 
    
    final_state *fs = new final_state(nu.hadronicFinalStateMC);
    
    
    //calculate cross section
    Double_t sigma = 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]));
    
    
    delete nu_lv;
    delete tar_lv;
    delete iact;
    delete fs;
    
    
    return sigma;
}

Bool_t NuTransition::IsSetToDoTransitions

(

)

const [inline]

Definition at line 27 of file NuTransition.h.

{return !doNothing;};

Float_t NuTransition::OscillationWeight	(	const Float_t	energy,
		const Int_t	inu,
		const NuXMLConfig *	lxmlConfig = 0
	)			const [private]

Definition at line 343 of file NuTransition.cxx.

References NuXMLConfig::DM2Bar(), NuXMLConfig::DM2Nu(), NuOscProb::OscillationWeight(), NuXMLConfig::SN2Bar(), and NuXMLConfig::SN2Nu().

Referenced by Reweight().

{
    const NuXMLConfig *xmlConfig;
    if (lxmlConfig) xmlConfig = lxmlConfig;
    else            xmlConfig = fxmlConfig;
        
    if (14==inu){
      return NuOscProb::OscillationWeight(energy,
                                            xmlConfig->DM2Nu(),
                                            xmlConfig->SN2Nu());
    }
    else if (-14==inu){
      return NuOscProb::OscillationWeight(energy,
                                            xmlConfig->DM2Bar(),
                                            xmlConfig->SN2Bar());
    }
    if (16==inu){
      return 1-NuOscProb::OscillationWeight(energy,
                                              xmlConfig->DM2Nu(),
                                              xmlConfig->SN2Nu());
    }
    else if (-16==inu){
      return 1-NuOscProb::OscillationWeight(energy,
                                              xmlConfig->DM2Bar(),
                                              xmlConfig->SN2Bar());
    }
    else return 1; // Doesn't oscillate
}

void NuTransition::Reweight

(

NuEvent &

nu

)

Definition at line 194 of file NuTransition.cxx.

References applyTranSyst, NuEvent::detector, NuXMLConfig::DM2Bar(), NuXMLConfig::DM2Nu(), doNothing, DoSystematicShifts(), firstRunTrans, flux_weight, fxmlConfig, NuEvent::iaction, NuEvent::inu, NuEvent::inunoosc, SimFlag::kData, Detector::kFar, Msg::kInfo, MAXMSG, NuEvent::neuEnMC, nubar_flux, numu_flux, OscillationWeight(), NuEvent::rw, NuEvent::simFlag, NuXMLConfig::SN2Bar(), NuXMLConfig::SN2Nu(), NuXMLConfig::TransitionProb(), and NuUtilities::ValueAt().

Referenced by DataMCPlots::PrepareEventForPlotting(), and NuDSTAna::StdNMBAna().

{
    // With no xml file, do neither systematics or oscillations
    if (doNothing) return;    
 
    double wTotal = 1;
    
    // If Systematics are applied to transitioned events,
    // can simply shift the whole event
    if (applyTranSyst) {
        this->DoSystematicShifts(nu);
    }
    else {
        NuEvent * nuCopy = (NuEvent*)nu.Clone();
        this->DoSystematicShifts(*nuCopy);
        wTotal *= nuCopy->rw / nu.rw;
        delete nuCopy;
    }
    
  //bool isNuE = (nu.inu == 12 || nu.inu == -12);
    bool isNuMu = (nu.inu == 14 || nu.inu == -14);
    bool isNuTau = (nu.inu == 16 || nu.inu == -16);
  
  
    if (SimFlag::kData == nu.simFlag) {
        // don't oscillate real data
        MAXMSG("NuTransition",Msg::kInfo,1) 
        << "Not applying fake oscillations to data" << endl;
    }
    else if (Detector::kFar != nu.detector) { 
        // don't oscillate the ND data
        MAXMSG("NuTransition",Msg::kInfo,1) 
        << "Not applying fake oscillations to ND MC" << endl;
    }
    else if (1 != nu.iaction) { 
        // check for !CC (i.e. NC)
        if(isNuTau) {
            MAXMSG("NuTransition",Msg::kInfo,1) << "Removing NC events from tau file (nu.rw set to zero)" << endl;
            wTotal = 0;
        }
        MAXMSG("NuTransition",Msg::kInfo,1) << "Not oscillating NCs" << endl;
    }
    else if ( !(isNuMu || isNuTau) ) { 
        // select everything that is not NuMu(bar) or NuTau(bar) i.e. true if nue(bar)
        MAXMSG("NuTransition",Msg::kInfo,1)
        << "Not oscillating nue events (non-muon/tau (anti)neutinos)" << endl;
    }
    else if ( isNuTau && (nu.inunoosc == 12 || nu.inunoosc==-12) ) {
        // select events that started as nues and were turned in to nutaus
        MAXMSG("NuTransition",Msg::kInfo,1)
        << "Removing beam nue events from tau file (nu.rw set to zero)" << endl;
        wTotal = 0;
    }
    else if (fxmlConfig->DM2Nu() < 0.0 ||  // Do this last, so that taus and nue's are still properly
             fxmlConfig->SN2Nu() < 0.0 ||  // reweighted, even in the no-oscillation case
             fxmlConfig->DM2Bar() < 0.0 || // For example, in the no-osc case the nutaus are given a weight of 0
             fxmlConfig->SN2Bar() < 0.0){  // so have to explictly just not touch them
        MAXMSG("NuTransition",Msg::kInfo,1)
        <<"Not applying fake oscillations due to xml configuration"<<endl;
    }
    else {
        //include the oscillation weight in the general nu.rw variable
        wTotal *= this->OscillationWeight(nu.neuEnMC, nu.inu);
        MAXMSG("NuTransition",Msg::kInfo,10)
        << "Oscillating with dm2=" << fxmlConfig->DM2Nu()
        << "; weight=" << this->OscillationWeight(nu.neuEnMC, nu.inu)
        << "; energy=" << nu.neuEnMC << endl;
        
        // Only create transitioned events if there are transitions
        if (fxmlConfig->TransitionProb() > 0) {
            if (isNuTau) {
                MAXMSG("NuTransition",Msg::kInfo,1)
                << "Weighting down appeared NuTaus for transitions by " 
                << (1.0 - fxmlConfig->TransitionProb()) << endl;
                wTotal *= (1.0 - fxmlConfig->TransitionProb());
            }
            else if (nu.inu == -14) { // and this is an antineutrino
                MAXMSG("NuTransition",Msg::kInfo,10) << "Getting a transitioned event" << endl;
                
                if (firstRunTrans) {
                    flux_weight->Divide(numu_flux, nubar_flux);
                    firstRunTrans = false;
                }            
                
                // Disappearance Probability = 1 - survival probability
                Double_t transitionWeight =  1.0 - this->OscillationWeight(nu.neuEnMC,nu.inu); 
                transitionWeight *= fxmlConfig->TransitionProb();
                transitionWeight *= NuUtilities::ValueAt(flux_weight, nu.neuEnMC);
                
                // Add, not multiply, the transition weight
                wTotal += transitionWeight;
            }
        }
    }
    
    nu.rw *= wTotal;
}

---

Member Data Documentation

bool NuTransition::applyTranSyst [private]

Definition at line 49 of file NuTransition.h.

Referenced by NuTransition(), and Reweight().

bool NuTransition::doNothing [private]

Definition at line 51 of file NuTransition.h.

Referenced by Fill(), NuTransition(), and Reweight().

bool NuTransition::firstRunTrans [private]

Definition at line 50 of file NuTransition.h.

Referenced by NuTransition(), and Reweight().

TH1D* NuTransition::flux_weight [private]

Definition at line 33 of file NuTransition.h.

Referenced by NuTransition(), Reweight(), and ~NuTransition().

TGraph* NuTransition::fNuBarXGraph [private]

Definition at line 35 of file NuTransition.h.

Referenced by GetXSecGraph(), and NuTransition().

TGraph* NuTransition::fNuMuXGraph [private]

Definition at line 34 of file NuTransition.h.

Referenced by GetXSecGraph(), and NuTransition().

NuSystematic* NuTransition::fSyst [private]

Definition at line 37 of file NuTransition.h.

Referenced by DoSystematicShifts(), and NuTransition().

NuXMLConfig* NuTransition::fxmlConfig [private]

Definition at line 36 of file NuTransition.h.

Referenced by DoSystematicShifts(), NuTransition(), and Reweight().

TH1D* NuTransition::nubar_flux [private]

Definition at line 32 of file NuTransition.h.

Referenced by Fill(), NuTransition(), Reweight(), and ~NuTransition().

TH1D* NuTransition::numu_flux [private]

Definition at line 27 of file NuTransition.h.

Referenced by Fill(), NuTransition(), Reweight(), and ~NuTransition().

---

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

• NuTransition.h
• NuTransition.cxx

---