samplePDFDUNEBeamNDGar Class Reference

Base class for handling beam ND GAR samples. More...

#include <samplePDFDUNE/samplePDFDUNEBeamNDGar.h>

Inheritance diagram for samplePDFDUNEBeamNDGar:

Collaboration diagram for samplePDFDUNEBeamNDGar:

Public Types
enum	KinematicTypes {   kTrueNeutrinoEnergy , kRecoNeutrinoEnergy , kTrueXPos , kTrueYPos ,   kTrueZPos , kTrueRad , kNMuonsRecoOverTruth , kRecoLepEnergy ,   kTrueLepEnergy , kRecoXPos , kRecoYPos , kRecoZPos ,   kRecoRad , kLepPT , kLepPZ , kPionMultiplicity ,   kNRecoParticles , kInFDV , kTrueMinusRecoEnergyRatio , kTrueMinusRecoEnergy ,   kNTrueMuons , kNRecoMuons }
	Enum to identify kinematics. More...

Public Member Functions
	samplePDFDUNEBeamNDGar (std::string mc_version, covarianceXsec *xsec_cov)
	Constructor.
	~samplePDFDUNEBeamNDGar ()
	destructor

Protected Member Functions
void	Init ()
	Initialises object.
int	setupExperimentMC (int iSample)
	Function to setup MC from file.
void	setupFDMC (int iSample)
	Tells FD base which variables to point to/be set to.
void	SetupWeightPointers ()
	Sets up pointers weights for each event (oscillation/xsec/etc.)
void	SetupSplines ()
	Sets up splines.
double	CalcXsecWeightFunc (int iSample, int iEvent)
	NOT IMPLEMENTED: Dunder method to calculate xsec weights.
void	applyShifts (int iSample, int iEvent)
	NOT IMPLEMENTED: Apply kinematic shifts.
const double *	GetPointerToKinematicParameter (KinematicTypes KinPar, int iSample, int iEvent)
	Returns pointer to kinemtatic parameter for event in Structs DUNE.
const double *	GetPointerToKinematicParameter (double KinematicVariable, int iSample, int iEvent)
	Returns pointer to kinemtatic parameter for event in Structs DUNE.
const double *	GetPointerToKinematicParameter (std::string KinematicParameter, int iSample, int iEvent)
	Returns pointer to kinemtatic parameter for event in Structs DUNE.
double	ReturnKinematicParameter (double KinematicVariable, int iSample, int iEvent)
	Returns pointer to kinemtatic parameter for event in Structs DUNE.
double	ReturnKinematicParameter (std::string KinematicParameter, int iSample, int iEvent)
	Returns pointer to kinemtatic parameter for event in Structs DUNE.
std::vector< double >	ReturnKinematicParameterBinning (std::string KinematicParameter)
	Gets binning for a given parameter.
std::vector< double >	ReturnKinematicParameterBinning (KinematicTypes KinematicParameter)
	Gets binning for a given parameter.
int	ReturnKinematicParameterFromString (std::string KinematicParameterStr)
	Get kinematic parameter ID from string name.
std::string	ReturnStringFromKinematicParameter (int KinematicParameter)
	Gets name of kinematic parmaeter.

Protected Attributes
std::vector< struct dunemc_base >	dunendgarmcSamples
	Array filled with MC samples for each oscillation channel.
TFile *	_sampleFile
	File containing sample objects.
TTree *	_data
	TTree containing sample Data.
double	pot
	Value of POT used for sample.
TString	_nutype
int	_mode
double	_ev
double	_erec
double	_erec_nue
double	_elep_reco
double	_LepNuAngle
int	_reco_numu
int	_reco_nue
double	_BeRPA_cvwgt = 1
int	_isCC
int	_nuPDGunosc
int	_nuPDG
int	_run
int	_isND
int	_isFHC
double	_vtx_x
double	_vtx_y
double	_vtx_z
double	_LepTheta
double	_Q2
bool	iscalo_reco
float	muonscore_threshold
caf::StandardRecord *	sr = new caf::StandardRecord()

Detailed Description

Base class for handling beam ND GAR samples.

Definition at line 24 of file samplePDFDUNEBeamNDGar.h.

Member Enumeration Documentation

KinematicTypes

enum samplePDFDUNEBeamNDGar::KinematicTypes

Enum to identify kinematics.

Enumerator
kTrueNeutrinoEnergy
kRecoNeutrinoEnergy
kTrueXPos
kTrueYPos
kTrueZPos
kTrueRad
kNMuonsRecoOverTruth
kRecoLepEnergy
kTrueLepEnergy
kRecoXPos
kRecoYPos
kRecoZPos
kRecoRad
kLepPT
kLepPZ
kPionMultiplicity
kNRecoParticles
kInFDV
kTrueMinusRecoEnergyRatio
kTrueMinusRecoEnergy
kNTrueMuons
kNRecoMuons

Definition at line 36 of file samplePDFDUNEBeamNDGar.h.

{kTrueNeutrinoEnergy, kRecoNeutrinoEnergy, kTrueXPos, kTrueYPos, kTrueZPos, kTrueRad, kNMuonsRecoOverTruth, kRecoLepEnergy, kTrueLepEnergy, kRecoXPos, kRecoYPos, kRecoZPos, kRecoRad, kLepPT, kLepPZ, kPionMultiplicity, kNRecoParticles, kInFDV, kTrueMinusRecoEnergyRatio, kTrueMinusRecoEnergy, kNTrueMuons, kNRecoMuons};

Constructor & Destructor Documentation

samplePDFDUNEBeamNDGar()

samplePDFDUNEBeamNDGar::samplePDFDUNEBeamNDGar	(	std::string	mc_version,
		covarianceXsec *	xsec_cov )

Constructor.

Parameters

mc_version	Configuration file
xsec_cov	cross-section covariance matrix

Definition at line 11 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                              : samplePDFFDBase(mc_version_, XsecCov_) {
  Initialise();
}

~samplePDFDUNEBeamNDGar()

samplePDFDUNEBeamNDGar::~samplePDFDUNEBeamNDGar

(

)

destructor

Definition at line 15 of file samplePDFDUNEBeamNDGar.cpp.

                                                {
}

Member Function Documentation

applyShifts()

void samplePDFDUNEBeamNDGar::applyShifts ( int iSample, int iEvent )

inlineprotected

NOT IMPLEMENTED: Apply kinematic shifts.

Parameters

iSample	Sample Number
iEvent	Event number

Definition at line 65 of file samplePDFDUNEBeamNDGar.h.

{}

CalcXsecWeightFunc()

double samplePDFDUNEBeamNDGar::CalcXsecWeightFunc ( int iSample, int iEvent )

inlineprotected

NOT IMPLEMENTED: Dunder method to calculate xsec weights.

Parameters

iSample	sample ID
iEvent	Event number

Definition at line 60 of file samplePDFDUNEBeamNDGar.h.

{return 1.;}

GetPointerToKinematicParameter() [1/3]

const double * samplePDFDUNEBeamNDGar::GetPointerToKinematicParameter ( double KinematicVariable, int iSample, int iEvent )

protected

Returns pointer to kinemtatic parameter for event in Structs DUNE.

Parameters

KinematicVariable	Kinematic parameter ID as double (gets cast -> int)
iSample	Sample ID
iEvent	Event ID

Returns

Value of kinetmatic parameter corresponding for a given event

Definition at line 309 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                                                      {
  KinematicTypes KinPar = (KinematicTypes) std::round(KinematicVariable);
  return GetPointerToKinematicParameter(KinPar,iSample,iEvent);
}

GetPointerToKinematicParameter() [2/3]

const double * samplePDFDUNEBeamNDGar::GetPointerToKinematicParameter ( KinematicTypes KinPar, int iSample, int iEvent )

protected

Returns pointer to kinemtatic parameter for event in Structs DUNE.

Parameters

KinPar	Kinematic parameter enum val
iSample	Sample ID
iEvent	Event ID

Returns

Pointer to KinPar for a given event

Definition at line 252 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                                                               {
  double* KinematicValue;
 
 switch(KinematicParameter) {
 case kTrueNeutrinoEnergy:
   KinematicValue = &dunendgarmcSamples[iSample].rw_etru[iEvent]; 
   break;
 case kRecoNeutrinoEnergy:
   KinematicValue = &dunendgarmcSamples[iSample].rw_erec[iEvent];
   break;
 case kTrueXPos:
   KinematicValue = &dunendgarmcSamples[iSample].rw_vtx_x[iEvent];
   break;
 case kTrueYPos:
   KinematicValue = &dunendgarmcSamples[iSample].rw_vtx_y[iEvent];
   break;
 case kTrueZPos:
   KinematicValue = &dunendgarmcSamples[iSample].rw_vtx_z[iEvent];
   break;
 case kTrueRad:
   KinematicValue = &dunendgarmcSamples[iSample].rw_rad[iEvent];
   break;
 case kNMuonsRecoOverTruth:
   KinematicValue = &dunendgarmcSamples[iSample].nmuonsratio[iEvent];
   break;
 case kRecoLepEnergy:
   KinematicValue = &dunendgarmcSamples[iSample].rw_elep_reco[iEvent];
   break;
 case kTrueLepEnergy:
   KinematicValue = &dunendgarmcSamples[iSample].rw_elep_true[iEvent];
   break;
 case kRecoXPos:
   KinematicValue = &dunendgarmcSamples[iSample].rw_reco_vtx_x[iEvent];
   break;
 case kRecoYPos:
   KinematicValue = &dunendgarmcSamples[iSample].rw_reco_vtx_y[iEvent];
   break;
 case kRecoZPos:
   KinematicValue = &dunendgarmcSamples[iSample].rw_reco_vtx_z[iEvent];
   break;
 case kRecoRad:
   KinematicValue = &dunendgarmcSamples[iSample].rw_reco_rad[iEvent];
   break;
 case kLepPT:
   KinematicValue = &dunendgarmcSamples[iSample].rw_lep_pT[iEvent];
   break;
 case kLepPZ:
   KinematicValue = &dunendgarmcSamples[iSample].rw_lep_pZ[iEvent];
   break;
 default:
   MACH3LOG_ERROR("Did not recognise Kinematic Parameter type...");
   throw MaCh3Exception(__FILE__, __LINE__);
 }
 
 return KinematicValue;
}

GetPointerToKinematicParameter() [3/3]

const double * samplePDFDUNEBeamNDGar::GetPointerToKinematicParameter ( std::string KinematicParameter, int iSample, int iEvent )

protected

Returns pointer to kinemtatic parameter for event in Structs DUNE.

Parameters

KinematicParameter	Kinematic parameter name as string (gets cast -> int)
iSample	Sample ID
iEvent	Event ID

Returns

Value of kinematic parameter corresponding for a given event

Definition at line 314 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                                                          {
  KinematicTypes KinPar = static_cast<KinematicTypes>(ReturnKinematicParameterFromString(KinematicParameter));
  return GetPointerToKinematicParameter(KinPar,iSample,iEvent);
}

Init()

void samplePDFDUNEBeamNDGar::Init ( )

protected

Initialises object.

Definition at line 18 of file samplePDFDUNEBeamNDGar.cpp.

                                  {
  dunendgarmcSamples.resize(nSamples,dunemc_base());
 
  iscalo_reco = SampleManager->raw()["SampleBools"]["iscalo_reco"].as<bool>(); //NK determine what reco used
}

ReturnKinematicParameter() [1/2]

double samplePDFDUNEBeamNDGar::ReturnKinematicParameter ( double KinematicVariable, int iSample, int iEvent )

protected

Returns pointer to kinemtatic parameter for event in Structs DUNE.

Parameters

KinematicVariable	Kinematic parameter ID as double (gets cast -> int)
iSample	Sample ID
iEvent	Event ID

Returns

Value of kinetmatic parameter corresponding for a given event

Definition at line 319 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                                         {
  return *GetPointerToKinematicParameter(KinematicVariable, iSample, iEvent);
}

ReturnKinematicParameter() [2/2]

double samplePDFDUNEBeamNDGar::ReturnKinematicParameter ( std::string KinematicParameter, int iSample, int iEvent )

protected

Returns pointer to kinemtatic parameter for event in Structs DUNE.

Parameters

KinematicParameter	Kinematic parameter name as string (gets cast -> int)
iSample	Sample ID
iEvent	Event ID

Returns

Value of kinematic parameter corresponding for a given event

Definition at line 323 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                                             {
  return *GetPointerToKinematicParameter(KinematicParameter, iSample, iEvent);
}

ReturnKinematicParameterBinning() [1/2]

std::vector< double > samplePDFDUNEBeamNDGar::ReturnKinematicParameterBinning ( KinematicTypes KinematicParameter )

protected

Gets binning for a given parameter.

Parameters

KinPar

Parameter ID

Returns

Vector containing parameter bins

Definition at line 350 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                               {
  std::vector<double> binningVector;
  switch(KinPar){
  case kTrueNeutrinoEnergy:
    for(double ibins =0; ibins<10*10; ibins++){
      double binval = ibins/10;
      binningVector.push_back(binval);
    }
    break;
  case kRecoNeutrinoEnergy:
    for(double ibins =0; ibins<10*10; ibins++){
      double binval = ibins/10;
      binningVector.push_back(binval);
    } 
    break;
  case kRecoXPos:
  case kTrueXPos:
    for(double ibins =0; ibins<259*2; ibins++){
      binningVector.push_back(ibins-259);
    }
    break;
  case kRecoYPos:
  case kTrueYPos:
    for(double ibins =0; ibins<277*2; ibins++){
      binningVector.push_back(ibins-277-150);
    }
    break;
  case kRecoZPos:
  case kTrueZPos:
    for(double ibins =0; ibins<277*2; ibins++){
      binningVector.push_back(ibins-277+1486);
    }
    break;
  case kPionMultiplicity:
    for(double ibins =0; ibins<10; ibins++){
      binningVector.push_back(ibins);
    }
    break;
  case kNRecoParticles:
    for(double ibins =0; ibins<50; ibins++){
      binningVector.push_back(ibins);
    }
    break; 
  case kInFDV:
    for(double ibins =0; ibins<3; ibins++){
      binningVector.push_back(ibins);
    }
    break;
  case kNMuonsRecoOverTruth:
  case kTrueMinusRecoEnergyRatio:
    for(double ibins =0; ibins<20*10; ibins++){
      binningVector.push_back(-10+(double)(ibins)/10);
    }
    break;
  case kTrueMinusRecoEnergy:
    for(double ibins =0; ibins<20*10; ibins++){
      binningVector.push_back(-10+(double)(ibins)/10);
    }
    break;
  case kNTrueMuons:
  case kNRecoMuons:
    for(double ibins =0; ibins<10; ibins++){
      binningVector.push_back(ibins);
    }
    break;
  case kRecoLepEnergy:
    for(double ibins =0; ibins<10*10; ibins++){
      binningVector.push_back((double)(ibins)/10);
    } 
    break;
  case kTrueLepEnergy:
    for(double ibins =0; ibins<10*10; ibins++){
      binningVector.push_back((double)(ibins)/10);
    } 
    break;
  case kTrueRad:
  case kRecoRad:
    for(double ibins =0; ibins<300; ibins++){
      binningVector.push_back(ibins);
    }
    break;
  case kLepPT:
  case kLepPZ:
    for(double ibins =0; ibins<10*10; ibins++){
      binningVector.push_back((double)(ibins)/10);
    }
    break;
  default:
    for(double ibins =0; ibins<10*100; ibins++){
      binningVector.push_back(ibins/100);
    }
    break;
  }
  
  return binningVector;
}

ReturnKinematicParameterBinning() [2/2]

std::vector< double > samplePDFDUNEBeamNDGar::ReturnKinematicParameterBinning ( std::string KinematicParameter )

protected

Gets binning for a given parameter.

Parameters

KinematicParameterStr

Parameter name

Returns

Vector containing parameter bins

Definition at line 345 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                                         {
  KinematicTypes KinPar = static_cast<KinematicTypes>(ReturnKinematicParameterFromString(KinematicParameterStr));
  return ReturnKinematicParameterBinning(KinPar);
}

ReturnKinematicParameterFromString()

int samplePDFDUNEBeamNDGar::ReturnKinematicParameterFromString ( std::string KinematicParameterStr )

protected

Get kinematic parameter ID from string name.

Parameters

KinematicStr

Returns

Parameter ID

Definition at line 447 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                              {
  return -1;
}

ReturnStringFromKinematicParameter()

std::string samplePDFDUNEBeamNDGar::ReturnStringFromKinematicParameter ( int KinematicParameter )

protected

Gets name of kinematic parmaeter.

Parameters

KinPar

Parameter ID

Returns

Name of parameter

Definition at line 451 of file samplePDFDUNEBeamNDGar.cpp.

                                                                                           {
  return "";
}

setupExperimentMC()

int samplePDFDUNEBeamNDGar::setupExperimentMC ( int iSample )

protected

Function to setup MC from file.

Parameters

iSample

sample ID

Returns

Total number of events

Definition at line 56 of file samplePDFDUNEBeamNDGar.cpp.

                                                         {
 
  dunemc_base *duneobj = &(dunendgarmcSamples[iSample]);
  int nutype = sample_nutype[iSample];
  int oscnutype = sample_oscnutype[iSample];
  bool signal = sample_signal[iSample];
  
  MACH3LOG_INFO("-------------------------------------------------------------------");
  MACH3LOG_INFO("Input File: {}", mc_files.at(iSample));
  
  _sampleFile = TFile::Open(mc_files.at(iSample).c_str(), "READ");
  _data = (TTree*)_sampleFile->Get("cafTree");
 
  if(_data){
    MACH3LOG_INFO("Found \"caf\" tree in {}", mc_files[iSample]);
    MACH3LOG_INFO("With number of entries: {}", _data->GetEntries());
  }
  else{
    MACH3LOG_ERROR("Could not find \"caf\" tree in {}", mc_files[iSample]);
    throw MaCh3Exception(__FILE__, __LINE__);
  }
  
  _data->SetBranchStatus("*", 1);
  _data->SetBranchAddress("rec", &sr);
 
  duneobj->norm_s = 1.0;
  duneobj->pot_s = (pot)/1e21;
 
  duneobj->nEvents = _data->GetEntries();
  duneobj->nutype = nutype;
  duneobj->oscnutype = oscnutype;
  duneobj->signal = signal;
 
  // allocate memory for dunendgarmc variables
  duneobj->rw_yrec = new double[duneobj->nEvents];
  duneobj->rw_elep_reco = new double[duneobj->nEvents];
  duneobj->rw_etru = new double[duneobj->nEvents];
  duneobj->rw_erec = new double[duneobj->nEvents];
  duneobj->flux_w = new double[duneobj->nEvents];
  duneobj->rw_isCC = new int[duneobj->nEvents];
  duneobj->rw_nuPDGunosc = new int[duneobj->nEvents];
  duneobj->rw_nuPDG = new int[duneobj->nEvents];
  duneobj->rw_berpaacvwgt = new double[duneobj->nEvents]; 
 
  duneobj->mode = new double[duneobj->nEvents];
 
  duneobj->nproton = new int[duneobj->nEvents];
  duneobj->nneutron = new int[duneobj->nEvents];
  duneobj->npip = new int[duneobj->nEvents];
  duneobj->npim = new int[duneobj->nEvents];
  duneobj->npi0 = new int[duneobj->nEvents];
 
  duneobj->nrecomuon = new int[duneobj->nEvents];
  duneobj->ntruemuon = new int[duneobj->nEvents];
  duneobj->nmuonsratio = new double[duneobj->nEvents];
  duneobj->ntruemuonprim = new int[duneobj->nEvents];
 
  duneobj->nrecoparticles = new int[duneobj->nEvents];
  duneobj->in_fdv = new bool[duneobj->nEvents];
  duneobj->rw_elep_true = new double[duneobj->nEvents];
 
  duneobj->rw_vtx_x = new double[duneobj->nEvents];
  duneobj->rw_vtx_y = new double[duneobj->nEvents];
  duneobj->rw_vtx_z = new double[duneobj->nEvents];
  duneobj->rw_rad = new double[duneobj->nEvents];
 
  duneobj->rw_lep_pT = new double[duneobj->nEvents];
  duneobj->rw_lep_pZ = new double[duneobj->nEvents];
 
  duneobj->rw_reco_vtx_x = new double[duneobj->nEvents];
  duneobj->rw_reco_vtx_y = new double[duneobj->nEvents];
  duneobj->rw_reco_vtx_z = new double[duneobj->nEvents];
  duneobj->rw_reco_rad = new double[duneobj->nEvents];
 
  duneobj->Target = new int[duneobj->nEvents];
 
  int num_no_ixns =0;
  int num_no_recparticles = 0;
  int num_in_fdv = 0;
  int num_in_fdv_noreco = 0;
  int num_notin_fdv =0;
  int num_nanenergy =0;
  int num_nanparticles =0;
 
  //FILL DUNE STRUCT
  for (int i = 0; i < (duneobj->nEvents); ++i) { // Loop through tree
    _data->GetEntry(i);
    double radius = pow((pow((sr->mc.nu[0].vtx.y+150),2) + pow((sr->mc.nu[0].vtx.z-1486),2)),0.5);
    if(std::abs(sr->mc.nu[0].vtx.x)<=209.0 &&  radius<=227.02){
      num_in_fdv++;
      duneobj->in_fdv[i] = 1;
    } else{
      num_notin_fdv++;
      duneobj->in_fdv[i] = 0;
    }
    
    if(sr->common.ixn.ngsft == 0){
      //duneobj->rw_erec[i] = (double)(0);
      float erec_total =0;
      duneobj->rw_elep_reco[i] = double(0);
      duneobj->rw_yrec[i] = (double)(0);
      num_no_ixns++;
      duneobj->nrecoparticles[i] = (int)(0);
    } else{
      duneobj->nrecoparticles[i] = (int)(0);
      float erec_total =0;
      float elep_reco =0;
      float muonscore = muonscore_threshold;
      int nixns = (int)(sr->common.ixn.ngsft);
      for(int i_ixn =0; i_ixn<nixns; i_ixn++) {
    int nrecoparticles = (int)(sr->common.ixn.gsft[i_ixn].part.ngsft);
    duneobj->nrecoparticles[i] += (int)(sr->common.ixn.gsft[i_ixn].part.ngsft);
    int nanparticles = 0;
    if(nrecoparticles ==0){
      double radius = pow((pow((sr->mc.nu[0].vtx.y+150.),2) + pow((sr->mc.nu[0].vtx.z-1486.),2)),0.5);
      if(std::abs(sr->mc.nu[0].vtx.x)<=209.0 || radius<=227.02) {
        num_in_fdv_noreco++;
      }   
      num_no_recparticles++;}
    for(int i_part =0; i_part<nrecoparticles; i_part++) {
      float erec_part = (float)(sr->common.ixn.gsft[i_ixn].part.gsft[i_part].E);
      if(std::isnan(erec_part)){nanparticles++;}
      erec_total+=erec_part;
      if((float)(sr->common.ixn.gsft[i_ixn].part.gsft[i_part].score.gsft_pid.muon_score>muonscore)){
        if(erec_part>elep_reco){
          elep_reco = erec_part;
          duneobj->rw_reco_vtx_x[i] = (double)((float)(sr->common.ixn.gsft[i_ixn].part.gsft[i_part].start.x));
          duneobj->rw_reco_vtx_y[i] = (double)((float)(sr->common.ixn.gsft[i_ixn].part.gsft[i_part].start.y));
          duneobj->rw_reco_vtx_z[i] = (double)((float)(sr->common.ixn.gsft[i_ixn].part.gsft[i_part].start.z));
          duneobj->rw_lep_pT[i] = (double)(pow(pow((float)(sr->common.ixn.gsft[i_ixn].part.gsft[i_part].p.x), 2) + pow((float)(sr->common.ixn.gsft[i_ixn].part.gsft[i_part].p.y), 2), 0.5));
          duneobj->rw_lep_pZ[i] = (double)(sr->common.ixn.gsft[i_ixn].part.gsft[i_part].p.z);
        }
        duneobj->nrecomuon[i]++; 
      }
    }
    num_nanparticles = num_nanparticles + (nanparticles/nrecoparticles);
      } //ADD PRIMARY LEPTON ENERGY ELEP_RECO
      if(std::isnan(erec_total)){num_nanenergy++; erec_total = (float)(sr->common.ixn.gsft[0].Enu.lep_calo);}
      if(iscalo_reco){duneobj->rw_erec[i]=(double)(sr->common.ixn.gsft[0].Enu.lep_calo);}
      else{duneobj->rw_erec[i]=(double)(erec_total);}
      duneobj->rw_elep_reco[i] = (double)(elep_reco);
    }
    
    if(duneobj->rw_erec[i] != 0){duneobj->rw_yrec[i] = (double)(((duneobj->rw_erec[i])-(duneobj->rw_elep_reco[i]))/(duneobj->rw_erec[i]));}
    else{duneobj->rw_yrec[i] = (double)(0);}
    duneobj->rw_etru[i] = (double)(sr->mc.nu[0].E); // in GeV
    duneobj->rw_isCC[i] = (int)(sr->mc.nu[0].iscc);
    duneobj->rw_nuPDGunosc[i] = sr->mc.nu[0].pdgorig;
    duneobj->rw_nuPDG[i] = sr->mc.nu[0].pdg;
    duneobj->rw_berpaacvwgt[i] = _BeRPA_cvwgt;
    
    int ntrueparticles = (int)(sr->mc.nu[0].nprim);
    for(int i_truepart =0; i_truepart<ntrueparticles; i_truepart++){
      if(std::abs(sr->mc.nu[0].prim[i_truepart].pdg) == 13){
    duneobj->ntruemuon[i]++;
    duneobj->ntruemuonprim[i]++;
      }
    }
    int ntruesecparticles = (int)(sr->mc.nu[0].nsec);
    for(int i_truepart =0; i_truepart<ntruesecparticles; i_truepart++){
      if(std::abs(sr->mc.nu[0].sec[i_truepart].pdg) == 13){
    duneobj->ntruemuon[i]++;
      }
    }
    
    duneobj->nproton[i] = sr->mc.nu[0].nproton;
    duneobj->nneutron[i] = sr->mc.nu[0].nneutron;
    duneobj->npip[i] = sr->mc.nu[0].npip;
    duneobj->npim[i] = sr->mc.nu[0].npim;
    duneobj->npi0[i] = sr->mc.nu[0].npi0;
    
    duneobj->nmuonsratio[i] = (double)(duneobj->nrecomuon[i])/(double)(duneobj->ntruemuonprim[i]);
    duneobj->rw_vtx_x[i] = (double)(sr->mc.nu[0].vtx.x);
    duneobj->rw_vtx_y[i] = (double)(sr->mc.nu[0].vtx.y);
    duneobj->rw_vtx_z[i] = (double)(sr->mc.nu[0].vtx.z);
    
    duneobj->rw_rad[i] = (double)(pow((pow((duneobj->rw_vtx_y[i]+150),2) + pow((duneobj->rw_vtx_z[i]-1486),2)),0.5)); 
    duneobj->rw_reco_rad[i] = (double)(pow(pow((duneobj->rw_reco_vtx_y[i]+150),2) + pow((duneobj->rw_reco_vtx_z[i]-1486), 2), 0.5));
    duneobj->rw_elep_true[i] = (double)(sr->mc.nu[0].prim[0].p.E);
    
    //Assume everything is on Argon for now....
    duneobj->Target[i] = 40;
    
    _mode = sr->mc.nu[0].mode;
    _isCC = (int)(sr->mc.nu[0].iscc);
    
    int mode= TMath::Abs(_mode);       
    duneobj->mode[i]=(double)GENIEMode_ToMaCh3Mode(mode, _isCC);
    
    duneobj->flux_w[i] = 1.0;
  }
 
  _sampleFile->Close();
  return duneobj->nEvents;
}

setupFDMC()

void samplePDFDUNEBeamNDGar::setupFDMC ( int iSample )

protected

Tells FD base which variables to point to/be set to.

Parameters

iSample

Sample ID

Definition at line 327 of file samplePDFDUNEBeamNDGar.cpp.

                                                  {
  dunemc_base *duneobj = &(dunendgarmcSamples[iSample]);
  fdmc_base *fdobj = &(MCSamples[iSample]);
  
  fdobj->nutype = duneobj->nutype;
  fdobj->oscnutype = duneobj->oscnutype;
  fdobj->signal = duneobj->signal;
  fdobj->SampleDetID = SampleDetID;
  
  for(int iEvent = 0 ;iEvent < fdobj->nEvents ; ++iEvent){
    fdobj->rw_etru[iEvent] = &(duneobj->rw_etru[iEvent]);
    fdobj->mode[iEvent] = &(duneobj->mode[iEvent]);
    fdobj->Target[iEvent] = &(duneobj->Target[iEvent]);
    
  }
  
}

SetupSplines()

void samplePDFDUNEBeamNDGar::SetupSplines ( )

protected

Sets up splines.

Todo

move all of the spline setup into core

Definition at line 24 of file samplePDFDUNEBeamNDGar.cpp.

                                          {
  if(XsecCov->GetNumParamsFromDetID(SampleDetID, kSpline) > 0){
    MACH3LOG_INFO("Found {} splines for this sample so I will create a spline object", XsecCov->GetNumParamsFromDetID(SampleDetID, kSpline));
    splinesDUNE* DUNESplines = new splinesDUNE(XsecCov);
    splineFile = (splineFDBase*)DUNESplines;
    InitialiseSplineObject();
  }
  else{
    MACH3LOG_INFO("Found {} splines for this sample so I will not load or evaluate splines", XsecCov->GetNumParamsFromDetID(SampleDetID, kSpline));
    splineFile = nullptr;
  }
 
  return;
}

SetupWeightPointers()

void samplePDFDUNEBeamNDGar::SetupWeightPointers ( )

protected

Sets up pointers weights for each event (oscillation/xsec/etc.)

Definition at line 41 of file samplePDFDUNEBeamNDGar.cpp.

                                                 {
  for (int i = 0; i < (int)dunendgarmcSamples.size(); ++i) {
    for (int j = 0; j < dunendgarmcSamples[i].nEvents; ++j) {
      MCSamples[i].ntotal_weight_pointers[j] = 6;
      MCSamples[i].total_weight_pointers[j] = new const double*[MCSamples[i].ntotal_weight_pointers[j]];
      MCSamples[i].total_weight_pointers[j][0] = &(dunendgarmcSamples[i].pot_s);
      MCSamples[i].total_weight_pointers[j][1] = &(dunendgarmcSamples[i].norm_s);
      MCSamples[i].total_weight_pointers[j][2] = MCSamples[i].osc_w_pointer[j];
      MCSamples[i].total_weight_pointers[j][3] = &(dunendgarmcSamples[i].rw_berpaacvwgt[j]);
      MCSamples[i].total_weight_pointers[j][4] = &(dunendgarmcSamples[i].flux_w[j]);
      MCSamples[i].total_weight_pointers[j][5] = &(MCSamples[i].xsec_w[j]);
    }
  }
}

Member Data Documentation

_BeRPA_cvwgt

double samplePDFDUNEBeamNDGar::_BeRPA_cvwgt = 1

protected

Definition at line 146 of file samplePDFDUNEBeamNDGar.h.

_data

TTree* samplePDFDUNEBeamNDGar::_data

protected

TTree containing sample Data.

Definition at line 129 of file samplePDFDUNEBeamNDGar.h.

_elep_reco

double samplePDFDUNEBeamNDGar::_elep_reco

protected

Definition at line 142 of file samplePDFDUNEBeamNDGar.h.

_erec

double samplePDFDUNEBeamNDGar::_erec

protected

Definition at line 140 of file samplePDFDUNEBeamNDGar.h.

_erec_nue

double samplePDFDUNEBeamNDGar::_erec_nue

protected

Definition at line 141 of file samplePDFDUNEBeamNDGar.h.

_ev

double samplePDFDUNEBeamNDGar::_ev

protected

Definition at line 139 of file samplePDFDUNEBeamNDGar.h.

_isCC

int samplePDFDUNEBeamNDGar::_isCC

protected

Definition at line 147 of file samplePDFDUNEBeamNDGar.h.

_isFHC

int samplePDFDUNEBeamNDGar::_isFHC

protected

Definition at line 152 of file samplePDFDUNEBeamNDGar.h.

_isND

int samplePDFDUNEBeamNDGar::_isND

protected

Definition at line 151 of file samplePDFDUNEBeamNDGar.h.

_LepNuAngle

double samplePDFDUNEBeamNDGar::_LepNuAngle

protected

Definition at line 143 of file samplePDFDUNEBeamNDGar.h.

_LepTheta

double samplePDFDUNEBeamNDGar::_LepTheta

protected

Definition at line 156 of file samplePDFDUNEBeamNDGar.h.

_mode

int samplePDFDUNEBeamNDGar::_mode

protected

Definition at line 135 of file samplePDFDUNEBeamNDGar.h.

_nuPDG

int samplePDFDUNEBeamNDGar::_nuPDG

protected

Definition at line 149 of file samplePDFDUNEBeamNDGar.h.

_nuPDGunosc

int samplePDFDUNEBeamNDGar::_nuPDGunosc

protected

Definition at line 148 of file samplePDFDUNEBeamNDGar.h.

_nutype

TString samplePDFDUNEBeamNDGar::_nutype

protected

Definition at line 134 of file samplePDFDUNEBeamNDGar.h.

_Q2

double samplePDFDUNEBeamNDGar::_Q2

protected

Definition at line 157 of file samplePDFDUNEBeamNDGar.h.

_reco_nue

int samplePDFDUNEBeamNDGar::_reco_nue

protected

Definition at line 145 of file samplePDFDUNEBeamNDGar.h.

_reco_numu

int samplePDFDUNEBeamNDGar::_reco_numu

protected

Definition at line 144 of file samplePDFDUNEBeamNDGar.h.

_run

int samplePDFDUNEBeamNDGar::_run

protected

Definition at line 150 of file samplePDFDUNEBeamNDGar.h.

_sampleFile

TFile* samplePDFDUNEBeamNDGar::_sampleFile

protected

File containing sample objects.

Definition at line 127 of file samplePDFDUNEBeamNDGar.h.

_vtx_x

double samplePDFDUNEBeamNDGar::_vtx_x

protected

Definition at line 153 of file samplePDFDUNEBeamNDGar.h.

_vtx_y

double samplePDFDUNEBeamNDGar::_vtx_y

protected

Definition at line 154 of file samplePDFDUNEBeamNDGar.h.

_vtx_z

double samplePDFDUNEBeamNDGar::_vtx_z

protected

Definition at line 155 of file samplePDFDUNEBeamNDGar.h.

dunendgarmcSamples

std::vector<struct dunemc_base> samplePDFDUNEBeamNDGar::dunendgarmcSamples

protected

Array filled with MC samples for each oscillation channel.

Definition at line 124 of file samplePDFDUNEBeamNDGar.h.

iscalo_reco

bool samplePDFDUNEBeamNDGar::iscalo_reco

protected

Definition at line 159 of file samplePDFDUNEBeamNDGar.h.

muonscore_threshold

float samplePDFDUNEBeamNDGar::muonscore_threshold

protected

Definition at line 160 of file samplePDFDUNEBeamNDGar.h.

pot

double samplePDFDUNEBeamNDGar::pot

protected

Value of POT used for sample.

Definition at line 131 of file samplePDFDUNEBeamNDGar.h.

sr

caf::StandardRecord* samplePDFDUNEBeamNDGar::sr = new caf::StandardRecord()

protected

Definition at line 162 of file samplePDFDUNEBeamNDGar.h.

---

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

• samplePDFDUNE/samplePDFDUNEBeamNDGar.h
• samplePDFDUNE/samplePDFDUNEBeamNDGar.cpp