test_resdec_piangle_rew.C File Reference

Go to the source code of this file.

Functions
	for (int i=0;i< ndial;i++)

Variables
double	p12iso = 0.50
double	p32rs = 0.75
double	p12rs = 0.25
const int	ndial = 11
const double	dial_step = 2./(ndial-1)
double	epsilon = 1E-6
const int	ncosthetabins =1000
TH1D *	h [ndial]

Function Documentation

for

(

)

Definition at line 16 of file test_resdec_piangle_rew.C.

References dial_step, MECModelEnuComparisons::i, calib::j, and ndial.

                             {
    double dial = -1. + i*dial_step;
    h[i] = new TH1D("a","",ncosthetabins,-1.,1.);
    h[i]->SetDirectory(0);
    if( TMath::Abs(dial)    < epsilon )  { h[i]->SetLineWidth(3); }
    if( TMath::Abs(dial-1.) < epsilon )  { h[i]->SetLineWidth(3); h[i]->SetLineStyle(kDashed); }
    for(int j=1; j<=h[i]->GetNbinsX(); j++) {
      double costheta = h[i]->GetBinCenter(j);
      double P2       = 0.5 * (3.*costheta*costheta - 1.);
      double Wiso     = 1 - p32iso * P2 + p12iso * P2; // = 1.0
      double Wrs      = 1 - p32rs  * P2 + p12rs  * P2;
      double Wdef     = Wiso;
      double Wtwk     = dial*Wrs + (1-dial)*Wiso;
      double wght = 1.;
      if(Wdef>0. && Wtwk>0.) {
        wght = Wtwk/Wdef;
      }
      h[i]->SetBinContent(j,wght);
    }//j
  }//i

Variable Documentation

const double dial_step = 2./(ndial-1)

Definition at line 8 of file test_resdec_piangle_rew.C.

Referenced by for().

double epsilon = 1E-6

Definition at line 10 of file test_resdec_piangle_rew.C.

Referenced by CheckChargeConservation(), CheckEnergyMomentumConservation(), genie::utils::math::CholeskyDecomposition(), DDTOnMonPlotMaker(), beamlinereco::WCHitFinderAlg::expandCluster(), ana::covmx::CovarianceMatrix::ForcePosDef(), genie::Intranuke::GenerateVertex(), genie::Intranuke2018::GenerateVertex(), GetBinByEdge(), isCloseEnough(), osetUtils::isEqual(), genie::BardinIMDRadCorPXSec::Li2(), ana::LikelihoodCovMxExperiment::LikelihoodCovMxExperiment(), ana::FluxDecomp::MakeWeightsNumuFromKaon(), ana::BENDecomp::MakeWeightsNumuFromKaon(), evg::MCTruthToDk2NuHackItr::pretty_darn_close(), evgen::CosmicsGen::ProjectMuonsToDetectorBigBox(), genie::SmithMonizUtils::QEL_EnuMin_SM(), beamlinereco::WCHitFinderAlg::run_DBSCAN(), genie::DISHadronicSystemGenerator::SimulateFormationZone(), genie::utils::intranuke::StepParticle(), genie::utils::intranuke2018::StepParticle(), genie::RosenbluthPXSec::XSec(), and ana::covmx::CovarianceMatrix::~CovarianceMatrix().

TH1D* h[ndial]

Definition at line 14 of file test_resdec_piangle_rew.C.

const int ncosthetabins =1000

Definition at line 12 of file test_resdec_piangle_rew.C.

const int ndial = 11

Definition at line 7 of file test_resdec_piangle_rew.C.

Referenced by for().

double p12iso = 0.50

Definition at line 3 of file test_resdec_piangle_rew.C.

double p12rs = 0.25

Definition at line 5 of file test_resdec_piangle_rew.C.

double p32rs = 0.75

Definition at line 4 of file test_resdec_piangle_rew.C.