#include <cassert>
#include <cstdlib>
#include <cctype>
#include <string>
#include <vector>
#include <sstream>
#include <map>
#include <TRotation.h>
#include <TH1D.h>
#include <TH3D.h>
#include <TTree.h>
#include <TLorentzVector.h>
#include "Framework/Algorithm/Algorithm.h"
#include "Framework/Algorithm/AlgFactory.h"
#include "Framework/EventGen/XSecAlgorithmI.h"
#include "Framework/Conventions/Constants.h"
#include "Framework/Conventions/Units.h"
#include "Framework/EventGen/GFluxI.h"
#include "Framework/EventGen/GMCJDriver.h"
#include "Framework/EventGen/GMCJMonitor.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGLibrary.h"
#include "Framework/Utils/XSecSplineList.h"
#include "Framework/Utils/StringUtils.h"
#include "Framework/Utils/SystemUtils.h"
#include "Framework/Utils/PrintUtils.h"
#include "Framework/Utils/UnitUtils.h"
#include "Framework/Utils/KineUtils.h"
#include "Framework/Utils/AppInit.h"
#include "Framework/Utils/RunOpt.h"
#include "Framework/Utils/CmdLnArgParser.h"
Go to the source code of this file.
Functions | |
void | GetCommandLineArgs (int argc, char **argv) |
void | PrintSyntax (void) |
GFluxI * | GetFlux (void) |
void | GenerateUpNu (GFluxI *flux_driver) |
TH3D * | BuildEmuEnuCosThetaPdf (int nu_code) |
TH1D * | GetEmuPdf (double Enu, double costheta, const TH3D *pdf3d) |
TVector3 | GetDetectorVertex (double CosTheta, double Enu) |
double | GetCrossSection (int nu_code, double Enu, double Emu) |
double | ProbabilityEmu (int nu_code, double Enu, double Emu) |
int | main (int argc, char **argv) |
Variables | |
Long_t | gOptRunNu |
string | gOptFluxSim |
map< int, string > | gOptFluxFiles |
int | gOptNev = -1 |
double | gOptDetectorSide |
long int | gOptRanSeed |
string | gOptInpXSecFile |
const double | a = 2e+6 |
const double | e = 500e+9 |
double | kDefOptDetectorSide = 1e+5 |
TH3D * BuildEmuEnuCosThetaPdf | ( | int | nu_code | ) |
Definition at line 444 of file gUpMuFluxGen.cxx.
References stan::math::exp(), h3, MECModelEnuComparisons::i, calib::j, test_ParserArtEvents::log, and ProbabilityEmu().
Referenced by main().
Definition at line 416 of file gUpMuFluxGen.cxx.
References genie::GFluxI::GenerateNext(), sanity_check_grl::keep, genie::kPdgAntiNuMu, genie::kPdgNuMu, LOG, genie::GFluxI::Momentum(), make_associated_cosmic_defs::p4, genie::GFluxI::PdgCode(), and pINFO.
Referenced by main().
Definition at line 643 of file gUpMuFluxGen.cxx.
References genie::CmdLnArgParser::ArgAsDouble(), genie::CmdLnArgParser::ArgAsInt(), genie::CmdLnArgParser::ArgAsLong(), genie::CmdLnArgParser::ArgAsString(), exit(), shutoffs::filename, genie::PDGLibrary::Find(), flux, genie::gAbortingInErr, gOptDetectorSide, gOptFluxFiles, gOptFluxSim, gOptInpXSecFile, gOptNev, gOptRanSeed, gOptRunNu, plot_validation_datamc::help, MECModelEnuComparisons::i, genie::PDGLibrary::Instance(), genie::RunOpt::Instance(), kDefOptDetectorSide, LOG, genie::CmdLnArgParser::OptionExists(), plot_validation_datamc::parser, pDEBUG, pFATAL, pINFO, pNOTICE, genie::utils::print::PrintFramedMesg(), PrintSyntax(), genie::RunOpt::ReadFromCommandLine(), and genie::utils::str::Split().
Referenced by main().
double GetCrossSection | ( | int | nu_code, |
double | Enu, | ||
double | Emu | ||
) |
Definition at line 538 of file gUpMuFluxGen.cxx.
References genie::AlgFactory::AdoptAlgorithm(), genie::units::cm2, genie::Interaction::DISCC(), dx, dy, MECModelEnuComparisons::i, genie::AlgFactory::Instance(), calib::j, genie::Interaction::KinePtr(), genie::constants::kNucleonMass, genie::kPdgNeutron, genie::kPdgProton, genie::kPdgTgtFreeN, genie::kPdgTgtFreeP, genie::kPSxyfE, genie::utils::kinematics::Q2(), genie::Kinematics::SetQ2(), genie::Kinematics::SetW(), genie::Kinematics::Setx(), genie::Kinematics::Sety(), W, submit_syst::x, genie::XSecAlgorithmI::XSec(), genie::utils::kinematics::XYtoWQ2(), submit_syst::y, and ymax.
Referenced by main(), and ProbabilityEmu().
TVector3 GetDetectorVertex | ( | double | CosTheta, |
double | Enu | ||
) |
Definition at line 288 of file gUpMuFluxGen.cxx.
References gOptDetectorSide, genie::RandomGen::Instance(), kPi, submit_hadd::l, make_associated_cosmic_defs::p3, Munits::rad, generate_hists::rnd, genie::RandomGen::RndFlux(), ana::Sqrt(), submit_syst::x, submit_syst::y, and test::z.
Referenced by main().
TH1D * GetEmuPdf | ( | double | Enu, |
double | costheta, | ||
const TH3D * | pdf3d | ||
) |
Definition at line 595 of file gUpMuFluxGen.cxx.
References genie::flux::GAtmoFlux::AddFluxFile(), exit(), shutoffs::filename, genie::gAbortingInErr, genie::flux::GAtmoFlux::GenerateWeighted(), gOptFluxFiles, gOptFluxSim, genie::flux::GAtmoFlux::LoadFluxData(), LOG, pFATAL, and genie::flux::GAtmoFlux::SetRadii().
Referenced by main().
Definition at line 172 of file gUpMuFluxGen.cxx.
References BuildEmuEnuCosThetaPdf(), genie::RunOpt::BuildTune(), exit(), GenerateUpNu(), GetCommandLineArgs(), GetCrossSection(), GetDetectorVertex(), GetEmuPdf(), GetFlux(), gOptInpXSecFile, gOptNev, gOptRanSeed, gOptRunNu, genie::RunOpt::Instance(), genie::kPdgAntiNuMu, genie::kPdgNuMu, LOG, genie::utils::app_init::MesgThresholds(), genie::GFluxI::Momentum(), outf, genie::GFluxI::PdgCode(), pFATAL, pNOTICE, genie::utils::app_init::RandGen(), genie::GFluxI::Weight(), Write(), and genie::utils::app_init::XSecTable().
Definition at line 843 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs().
double ProbabilityEmu | ( | int | nu_code, |
double | Enu, | ||
double | Emu | ||
) |
Definition at line 435 of file gUpMuFluxGen.cxx.
References a, e, GetCrossSection(), and genie::constants::kNA.
Referenced by BuildEmuEnuCosThetaPdf().
const double a = 2e+6 |
Definition at line 164 of file gUpMuFluxGen.cxx.
Referenced by ana::_HistAxis< T >::_HistAxis(), genie::Spline::Add(), evdb::View2D::AddArc(), ana::KrigeKernel::AddPoint(), ana::AlmostEqual(), slicer::S4DParamCalc::analyze(), bpfit::BPFdEdxHistoMaker::analyze(), slicer::SlicerAna::analyze(), bpfit::BPFTmvaTrainer::analyze(), air::AirKalmanAna::analyze(), novaddt::HotMapMaker::AppendCold(), novaddt::HotMapMaker::AppendHot(), novaddt::HotMapMaker::AppendRate(), novaddt::utils::DetectorUtils::AttenuationInFiber(), genie::KNOHadronization::AverageChMult(), gibuu::GiBUURegen::Basis(), stan::math::beta_rng(), calib::BetheBlochAnalytic::BetheBlochAnalytic(), genie::utils::nuclear::BindEnergy(), by_charge(), genie::QPMDISStrucFuncBase::Calculate(), beamlinereco::WCTrackAlg::calculateTheMomentum(), earms::GridSearch::ChooseBestProngDir(), genie::utils::kinematics::CohQ2Lim(), ana::CombineExposures(), boost::python::map_indexing_suite< Container, NoProxy, final_map_derived_policies< Container, NoProxy > >::compare_index(), genie::RSHelicityAmplModelCC::Compute(), genie::RSHelicityAmplModelNCp::Compute(), genie::RSHelicityAmplModelNCn::Compute(), genie::alvarezruso::AREikonalSolution::Con(), CLHEP::MixMaxRng::convert1double(), genie::geometry::PlaneParam::ConvertMaster2Top(), genie::mueloss::BetheBlochModel::dE_dx(), fhicl::detail::decode(), fhicl::detail::decode_tuple(), calib::BetheBlochAnalytic::Delta(), genie::Spline::Divide(), genie::mueloss::gsl::BezrukovBugaevIntegrand::DoEval(), genie::mueloss::gsl::PetrukhinShestakovIntegrand::DoEval(), genie::utils::gsl::dXSec_Log_Wrapper::DoEval(), stan::mcmc::softabs_metric< Model, BaseRNG >::dphi_dq(), stan::mcmc::softabs_metric< Model, BaseRNG >::dtau_dq(), earms::ElasticArms::E(), earmsnt(), dt::Cand::EstimateStraightLine(), osc::_PMNS< T >::EvalEqn11(), osc::_PMNS< T >::EvalEqn2(), stan::math::fabs(), fiducial_mass(), art::View< T >::fill(), fhicl::detail::find_an_any(), util::find_median(), ana::GenieMultiverseSpectra::FindBandBoundaries(), lem::FindMatches(), cmf::SpectrumPredictionMaker::FindOscillationPoints(), fit_thresh_corrs(), fit_thresh_corrs_fb(), stan::math::fmax(), stan::math::fmin(), g4lbneCrossChecks(), g4minervaCrossChecks(), GammaFit(), stan::math::gaussian_dlm_obs_lpdf(), genie::BBA03ELFormFactorsModel::Gen(), stan::lang::generate_register_mpi(), ana::get3FlavorAna2020SmallXsecSysts(), fhicl::ParameterSet::get_one_(), ana::getAna2017SmallXsecSysts(), ana::getAna2018SmallXsecSysts(), genie::Spline::GetAsTSpline(), numusand::NumuSandFxs::getAveTrackdEdxLast4Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast6Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast8Cells(), trackinfo::TrackInfoFxs::getAveTrackdEdxLastCells(), HepGeom::BasicVector3D< float >::getEta(), genie::alvarezruso::AlvarezRusoCOHPiPDXSec::GetLeptonMass(), evg::MCTruthToDk2NuHackItr::GetMCTruthPtr(), ana::GenieMultiverseSpectra::GetSigmaHistograms(), LSTME::Png3dVars::getSortedIndices(), SliceLID::Png3dVars::getSortedIndices(), stan::math::grad_inc_beta(), stan::math::grad_reg_lower_inc_gamma(), cheat::BackTracker::HitCollectionEfficiency(), stan::math::hypergeometric_log(), stan::math::hypergeometric_lpmf(), stan::math::inc_beta_dda(), genie::utils::kinematics::InelYLim_X(), genie::utils::kinematics::electromagnetic::InelYLim_X(), daqdataformats::VERSION_NAMESPACE::RawConfigurationBlock::init(), daqdataformats::VERSION_NAMESPACE::RawSummaryDroppedMicroblock::init(), daqdataformats::VERSION_NAMESPACE::RawTriggerTime::init(), rb::Track::InterpolateDir(), rb::Track::InterpolateXY(), stan::math::inv_Phi(), fhicl::ParameterSet::is_key_to_atom(), isCloseEnough(), genie::geometry::PlaneParam::IsValid(), hough::IslandsTH2::LabelIslands(), ana::LabelsAndBins::LabelsAndBins(), earms::ElasticArms::LeadingArm(), novaddt::NNFastMMTrigger::LinFit(), zcl::FmmTrackerValidation::LinFit(), novaddt::UpMuTrigger::LinFit(), zcl::FastMMStudy::LinFit(), novaddt::UpMuAna::LinFit(), upmuana::UpMuAnalysis::LinFit(), upmuana::UpMuRecoAna::LinFit(), novaddt::UpMuTrigger::LinFitLLR(), novaddt::UpMuAna::LinFitLLR(), upmuana::LlrUpmu::LLR(), upmuana::UpMuAnalysis::LLR(), upmuana::UpMuRecoAna::LLR(), genie::NievesQELCCPXSec::LmunuAnumu(), nerd::NERDEval::load_anchors(), g4n::LoadNeutronFates(), genie::MECHadronTensor::LoadTensorTables(), stan::math::log_modified_bessel_first_kind(), make_cm_pullterm_pdf(), make_xsec_tuning_hists_mp5(), make_xsec_wgts_2018_hists(), nerd::NERDProng::MakeCluster(), novaddt::HotMapMaker::MakeMask(), MakeTestStandPlots_AllRuns(), MakeTestStandPlots_NewRuns(), stan::math::matrix_exp_2x2(), MAX(), nl::TimeSeriesData::MeanValue(), dt::DiscreteTracker::MergeParasiteTracks(), ana::PredictionExtendToPeripheral::MergePeripheral(), min(), evd::RawDrawingOptions::ModuleList(), evd::RecoDrawingOptions::ModuleList(), hough::MultiHough2P::MultiMap(), genie::Spline::Multiply(), ana::TargetCount::NNuclei(), genie::geometry::PlaneParam::Normalize(), ana::Nue2018ThreeBinAxis(), Nue2018ThreeBinAxis_Reduced(), ana::Nue2018ThreeBinAxis_Reduced(), nus17_fiducial_mass(), OnMonPlotMaker(), ana::_Cut< T >::operator!(), art::operator!(), art::operator!=(), art::operator&&(), ana::_Cut< T >::operator&&(), order_cuts::operator()(), order_chns::operator()(), ana::_Var< T >::operator*(), stan::math::operator*(), HepGeom::BasicVector3D< T >::operator*(), HepGeom::BasicVector3D< float >::operator*=(), ana::_Var< T >::operator+(), stan::math::operator+(), stan::math::operator++(), ana::_Var< T >::operator-(), stan::math::operator-(), stan::math::operator--(), ana::_Var< T >::operator/(), stan::math::operator/(), HepGeom::BasicVector3D< T >::operator/(), HepGeom::BasicVector3D< float >::operator/=(), art::operator<(), lem::LiteHit::operator<(), ana::operator<(), nova::database::Column::operator<(), nutools::dbi::Column::operator<(), nova::dbi::Column::operator<(), om::OStream::operator<<(), art::operator<<(), nova::dbi::operator<<(), nova::database::operator<<(), genie::flux::operator<<(), nova::database::Column::operator<=(), nutools::dbi::Column::operator<=(), nova::dbi::Column::operator<=(), art::operator==(), nova::database::Column::operator==(), art::operator>(), ana::operator>(), nova::database::Column::operator>(), nutools::dbi::Column::operator>(), nova::dbi::Column::operator>(), nova::database::Column::operator>=(), nutools::dbi::Column::operator>=(), nova::dbi::Column::operator>=(), ana::_Cut< T >::operator||(), art::operator||(), pie_arc(), genie::geometry::PlaneParam::PlaneParam(), plot_3NDvsFD_FHC(), plot_3NDvsFD_RHC(), PlotPionPlots(), sim::PhotonSignal::PoissonLambda(), g4n::NeutronSubstitutionProcess::PostStepDoIt(), pots(), genie::P33PaschosLalakulichPXSec::PPiStar(), evg::MCTruthToDk2NuHackItr::pretty_darn_close(), genie::geometry::PlaneParam::Print(), ProbabilityEmu(), genie::FGMBodekRitchie::ProbDistro(), genie::LocalFGM::ProbDistro(), cvntf::CVNCosmicTFFilt::produce(), cvn::CVNCosmicMapper::produce(), earms::ElasticArmsHS::produce(), bpfit::BPFEnergyEstimator::produce(), bpfit::BPFPIdMaker::produce(), genie::SmithMonizUtils::Q2lim1_SM(), genie::SmithMonizUtils::Q2lim2_SM(), ana::RefineSeeds(), HepGeom::Reflect3D::Reflect3D(), genie::NievesQELCCPXSec::relLindhardIm(), hough::MultiHough2P::ReweightHits(), genie::alvarezruso::integrationtools::RG202D(), genie::alvarezruso::integrationtools::RG482D(), art::right(), HepGeom::Rotate3D::Rotate3D(), stan::model::rvalue(), stan::mcmc::softabs_metric< Model, BaseRNG >::sample_p(), genie::BYStrucFunc::ScalingVar(), stan::math::matrix_exp_action_handler::set_approximation_parameter(), setA(), earms::ElasticArms::SetArm(), ana::covmx::Sample::SetAxis(), SetGraphErrors(), ana::FitSinSq2Theta13::SetValue(), daqdataformats::VERSION_NAMESPACE::RawMicroBlockHeader::setVersion(), HepGeom::BasicVector3D< float >::setX(), HepGeom::BasicVector3D< float >::setY(), HepGeom::BasicVector3D< float >::setZ(), genie::alvarezruso::integrationtools::SG20R(), genie::alvarezruso::integrationtools::SG48R(), show_x_cut(), slid::ShowerLID::ShowerLID(), hough::Hough2P::SmoothMap(), hough::MultiHough2P::SmoothMap(), solve(), ana::MCMCSamples::SortedLLs(), evd::RawDataDrawer::StickFigure(), earms::ElasticArms::SumVi(), art::ValidHandle< T >::swap(), TEST(), test_vec(), calib::ThresholdCorrMap::ThresholdCorrAt(), art::detail::throw_if_invalid(), lem::dec::Forest::TrainSingle(), rb::Track::Trajectory(), ncid::NCNNKeras::TransMomFraction(), earms::ElasticArms::UpdateDaa(), earms::ElasticArms::UpdateMia(), earms::ElasticArms::UpdateVia(), uptime(), osc::VacuumHamiltonian(), ana::Var2D(), ana::Var3D(), genie::geometry::PlaneParam::Vd(), genie::geometry::PlaneParam::Vn(), genie::SmithMonizUtils::vQES_SM_lim(), fhicl::detail::ParameterWalker< C >::walk_over(), genie::BSKLNBaseRESPXSec2014::XSec(), Y11(), Y12(), Y22(), rb::Track::ZipWith(), and calib::ThresholdCorrMap::~ThresholdCorrMap().
const double e = 500e+9 |
Definition at line 165 of file gUpMuFluxGen.cxx.
Referenced by ProbabilityEmu().
double gOptDetectorSide |
Definition at line 159 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs(), and GetDetectorVertex().
Definition at line 157 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs(), and GetFlux().
string gOptFluxSim |
Definition at line 156 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs(), and GetFlux().
string gOptInpXSecFile |
Definition at line 161 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs(), and main().
int gOptNev = -1 |
Definition at line 158 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs(), and main().
long int gOptRanSeed |
Definition at line 160 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs(), and main().
Long_t gOptRunNu |
Definition at line 155 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs(), and main().
double kDefOptDetectorSide = 1e+5 |
Definition at line 169 of file gUpMuFluxGen.cxx.
Referenced by GetCommandLineArgs().