Reweighter of thin target nC interactions. More...

#include <ThinTargetnCPionReweighter.h>

Inheritance diagram for NeutrinoFluxReweight::ThinTargetnCPionReweighter:

Public Member Functions
	ThinTargetnCPionReweighter (int iuniv, const ParameterTable &cv_pars, const ParameterTable &univ_pars)

virtual	~ThinTargetnCPionReweighter ()

virtual bool	canReweight (const InteractionData &aa)
	can the particular instance of this class reweight this interaction? More...

virtual double	calculateWeight (const InteractionData &aa)
	calculate a weight for this interaction given the central value parameters and the parameters for this universe. The weight is something like: f(cv)/f(MC) * f(univ)/f(cv) where cv in this case corresponds to the best value of the parameter, given the data. If univ_pars=cv_pars then we are calculating a central value weight More...

Public Member Functions inherited from NeutrinoFluxReweight::IInteractionReweighting
virtual	~IInteractionReweighting ()

Private Attributes
int	iUniv

const ParameterTable &	cvPars

const ParameterTable &	univPars

ThinTargetpCPionReweighter *	tt_pCPionRew

Detailed Description

Reweighter of thin target nC interactions.

Definition at line 20 of file ThinTargetnCPionReweighter.h.

Constructor & Destructor Documentation

NeutrinoFluxReweight::ThinTargetnCPionReweighter::ThinTargetnCPionReweighter	(	int	iuniv,
		const ParameterTable &	cv_pars,
		const ParameterTable &	univ_pars
	)

Definition at line 8 of file ThinTargetnCPionReweighter.cpp.

                                                                                                                                  :iUniv(iuniv),cvPars(cv_pars),univPars(univ_pars){
     
   }

NeutrinoFluxReweight::ThinTargetnCPionReweighter::~ThinTargetnCPionReweighter ( )

virtual

Definition at line 12 of file ThinTargetnCPionReweighter.cpp.

                                                           {
     
   }

Member Function Documentation

double NeutrinoFluxReweight::ThinTargetnCPionReweighter::calculateWeight ( const InteractionData & inter_data )

virtual

calculate a weight for this interaction given the central value parameters and the parameters for this universe. The weight is something like: f(cv)/f(MC) * f(univ)/f(cv) where cv in this case corresponds to the best value of the parameter, given the data. If univ_pars=cv_pars then we are calculating a central value weight

Implements NeutrinoFluxReweight::IInteractionReweighting.

Definition at line 44 of file ThinTargetnCPionReweighter.cpp.

                                                                              {
     
     MakeReweight*  makerew =  MakeReweight::getInstance();
     if(iUniv==-1)tt_pCPionRew = (makerew->cv_rw)->THINTARGET_PC_PION_Universe;
     else tt_pCPionRew = (makerew->vec_rws[iUniv])->THINTARGET_PC_PION_Universe;
 
     double inc_mom[3]  = {aa.Inc_P4[0], aa.Inc_P4[1], aa.Inc_P4[2]};
     double prod_mom[3] = {aa.Prod_P4[0],aa.Prod_P4[1],aa.Prod_P4[2]};
     double vtx_int[3]  = {aa.Vtx[0],aa.Vtx[1],aa.Vtx[2]};
     InteractionData intData(aa.gen, inc_mom,2212,prod_mom,-1*aa.Prod_pdg,aa.Vol,aa.Proc,vtx_int);
 
     double wgt = tt_pCPionRew->calculateWeight(intData);
 
     if(wgt<0)std::cout<<"TTNCPI check wgt(<0) "<<iUniv<<" "<<aa.Inc_P<<" "<<aa.xF<<" "<<aa.Pt<<" "<<aa.Prod_pdg<<std::endl;
     return wgt;
       
   }

bool NeutrinoFluxReweight::ThinTargetnCPionReweighter::canReweight ( const InteractionData & aa )

virtual

can the particular instance of this class reweight this interaction?

Implements NeutrinoFluxReweight::IInteractionReweighting.

Definition at line 15 of file ThinTargetnCPionReweighter.cpp.

                                                                        {
 
     MakeReweight*  makerew =  MakeReweight::getInstance();
     if(iUniv==-1)tt_pCPionRew = (makerew->cv_rw)->THINTARGET_PC_PION_Universe;
     else tt_pCPionRew = (makerew->vec_rws[iUniv])->THINTARGET_PC_PION_Universe;
     
     //checking:
     std::string mode(getenv("MODE"));
     if(aa.Inc_pdg != 2112)return false;
     if(aa.Inc_P < 12.0)return false;
     //volume check: 
     bool is_wrong_volume = aa.Vol != "TGT1" && aa.Vol != "BudalMonitor" && aa.Vol != "Budal_HFVS" && aa.Vol != "Budal_VFHS";
     if( (mode=="REF") || (mode=="OPT") ){
      // is_wrong_volume = aa.Vol != "TargetFinHorizontal" && aa.Vol != "TargetNoSplitSegment";
       is_wrong_volume = aa.Vol != "TargetFinHorizontal" && aa.Vol != "TargetNoSplitSegment" && aa.Vol!="tCoreLog";     
     }
     if(is_wrong_volume)return false;
     //
     if(aa.Prod_pdg != 211 && aa.Prod_pdg != -211)return false;
     
     double inc_mom[3]  = {aa.Inc_P4[0], aa.Inc_P4[1], aa.Inc_P4[2]};
     double prod_mom[3] = {aa.Prod_P4[0],aa.Prod_P4[1],aa.Prod_P4[2]};
     double vtx_int[3]  = {aa.Vtx[0],aa.Vtx[1],aa.Vtx[2]};
     
     InteractionData intData(aa.gen, inc_mom,2212,prod_mom,-1*aa.Prod_pdg,aa.Vol,aa.Proc,vtx_int);
     return tt_pCPionRew->canReweight(intData);
     
   }