comparison.C File Reference

#include "Riostream.h"

Go to the source code of this file.

Functions
experimentalTree	SetBranchAddress ("EdepExp",&EdepExp)
experimentalTree	SetBranchAddress ("depthExp",&depthExp)
	for (Int_t i=0;i< nentries;i++)
simulatedPeak	Scale (1/simulationNormalisationFactor)
ntupleExperimental	SetMarkerStyle (4)
simulatedPeak	SetMarkerSize (2)
ntupleExperimental	Draw ("EdepExp:depthExp")
simulatedPeak	Draw ("same")
leg	SetTextSize (0.035)
leg	SetFillColor (0)
leg	AddEntry (ntupleExperimental,"Experiment","P")
leg	AddEntry (simulatedPeak,"Simulation")
leg	Draw ()

Variables
ifstream	in
TFile *	experimentalFile = new TFile("../../../experimentalData/proton/BraggPeak/62MeVInWater.root","READ")
TTree *	experimentalTree = (TTree*)experimentalFile -> Get("Experimental62MeVInWater")
Float_t	depthExp
Float_t	EdepExp
TNtuple *	ntupleExperimental = new TNtuple("ntupleExperimental","Protons, exp. data", "depthExp:EdepExp")
Int_t	nentries = (Int_t)experimentalTree -> GetEntries()
TFile *	simulationFile = new TFile("../../../SimulationOutputs/proton/BraggPeak/protonBraggPeak.root","READ")
TH1D *	simulatedPeak = (TH1D*) simulationFile -> Get("braggPeak")
Float_t	simulationNormalisationFactor = simulatedPeak -> GetBinContent(1)
TCanvas *	c1 = new TCanvas ("c1","c1",200,10,600,400)
	leg = new TLegend(0.50,0.60,0.20,0.70)

Function Documentation

leg AddEntry	(	ntupleExperimental	,
		"Experiment"	,
		"P"
	)

leg AddEntry	(	simulatedPeak	,
		"Simulation"
	)

ntupleExperimental Draw

(

"EdepExp:depthExp"

)

simulatedPeak Draw

(

"same"

)

leg Draw

(

)

for

(

)

Definition at line 30 of file comparison.C.

References EdepExp, Fill(), GetEntry(), MECModelEnuComparisons::i, and normFactor.

   {
     experimentalTree -> GetEntry(0);
     Float_t normFactor = EdepExp;
     experimentalTree -> GetEntry(i);
     ntupleExperimental -> Fill(depthExp, EdepExp/normFactor);
     
   }

simulatedPeak Scale

(

1/

simulationNormalisationFactor

)

Referenced by CheckVertexDistribution(), DataMCPair::DrawData(), DataMCPair::DrawMCBkg(), DataMCPair::DrawMCNormSyst(), DataMCPair::DrawMCSyst(), DataMCPair::DrawWrong(), extract_hadnucleus_xsec(), ana::Hist::GetEigenStan(), GetFNRatio(), GetSpectToHist(), hA_basic_analysis(), makeHadEFrac2DPlot(), MakeOverlayCanvas(), make_reco_miniprod5_plots::MakePlots(), make_reco_plots::MakePlots(), MakeSplitCans(), NormalizeHists(), nuint09_1pi2(), nuint09_1pi3(), nuint09_1pi4(), nuint09_coh2(), nuint09_coh3(), nuint09_qel2(), nuint09_qel3(), nuint09_qel4(), nuint09_qel5(), nuint09_qel6(), NuMu2020_MakePIDPlots(), plot_ehade_quantbound(), plot_quantile_boundaries_2020(), SavePCAShifts(), ana::Scaled(), SuppressND(), and ana::BayesianMarginal::ThresholdFromTH1().

experimentalTree SetBranchAddress	(	"EdepExp"	,
		&	EdepExp
	)

experimentalTree SetBranchAddress	(	"depthExp"	,
		&	depthExp
	)

leg SetFillColor

(

0

)

h2 SetMarkerSize

(

2

)

ntupleExperimental SetMarkerStyle

(

4

)

leg SetTextSize

(

0.

035

)

Variable Documentation

TCanvas* c1 = new TCanvas ("c1","c1",200,10,600,400)

Definition at line 51 of file comparison.C.

Float_t depthExp

Definition at line 22 of file comparison.C.

Float_t EdepExp

Definition at line 22 of file comparison.C.

Referenced by for().

TFile* experimentalFile = new TFile("../../../experimentalData/proton/BraggPeak/62MeVInWater.root","READ")

Definition at line 17 of file comparison.C.

TTree* experimentalTree = (TTree*)experimentalFile -> Get("Experimental62MeVInWater")

Definition at line 20 of file comparison.C.

ifstream in

Initial value:

{  
gROOT->Reset()

Definition at line 7 of file comparison.C.

Referenced by generic.views.Pages.PlotPage::AddFigure(), genie::Algorithm::AllowReconfig(), DataFile::analyzeTextFile(), EliminateBeamSpills::beginJob(), braggPeak(), checkBeam(), stan::test::performance::command(), art::detail::TwoArgInsert< cet::map_vector< T > >::concatenate(), skim::NueSkimmer::CopyMichelSlice(), skim::NueSkimmer::CopyMichelTrack(), skim::SkimmingUtils::CopyProductAndSliceAssn(), demo_trivial_xsec_analysis(), diffapdhv(), draw_decomp_plots(), draw_michel_plots(), jmshower::EID::EID(), genie::flux::GHAKKMAtmoFlux::FillFluxHisto(), genie::flux::GFLUKAAtmoFlux::FillFluxHisto(), genie::flux::GBGLRSAtmoFlux::FillFluxHisto(), om::HistoSourceBase::FindAllMatches(), ana::FixPlotName(), FixSystName(), art::flattenCollections(), fragmentAngularDistribution(), fragmentAngularDistributionGM(), fragmentAngularDistributionHistogram(), fragmentEnergy(), fragmentEnergyDistributionDifferentAngles(), fragmentYields(), fragmentYieldsPlot(), get_files(), upmuana::UpMuAnalysis::getLLR(), upmuana::UpMuRecoAna::getLLR(), genie::NtpMCTreeHeader::Init(), genie::EmpiricalMECPXSec2015::Integral(), genie::ReinSehgalRESXSec::Integrate(), genie::ReinSehgalRESXSecFast::Integrate(), InvestigatePeriods(), dbutils::gui::DAQConfigEditorDisplay::loadThresholdsAndMasks(), lsqr_fprime(), main(), makeHWwatchlist(), ana::MakeLatexCommandName(), MakeLatexCommandName(), MakeTestStandPlots_AllRuns(), MakeTestStandPlots_NewRuns(), MakeValidation(), mix::MixSimAndCosmicData::MixCellHitDetails(), mix::MixSimAndCosmicData::MixFlatDAQDataDetails(), mix::MixSimAndCosmicData::MixFLSHitsDetails(), mix::MixSimAndCosmicData::mixMCFlux(), mix::MixSimAndCosmicData::MixMCTruthDetails(), mix::MixSimAndCosmicData::MixPhotonSignalDetails(), mix::MixSimAndCosmicData::MixRawDigitDetails(), mix::MixSimAndCosmicData::MixRawTriggerDetails(), mix::MixSimAndCosmicData::MixSimParticleDetails(), mix::MixSimAndCosmicData::mixSimTrueEnergy(), ndrockmacro(), Nearline_Watchdog_Plots(), art::detail::mix_offset< cet::map_vector< P > >::offset(), OnMonAnaHistos(), OnMonAnaHistosFEB_IssueRates(), OnMonAnaHistosPIX(), art::PtrRemapper::operator()(), stan::json::parse(), genie::utils::kinematics::PhaseSpaceVolume(), plotapdhv(), PlotCSV(), genie::utils::print::PrintBanner(), PrintTable(), FormatFile::process(), read_file(), readExfor(), genie::MECHadronTensor::ReadHadTensorqzq0File(), genie::INukeHadroData2018::ReadhNFile(), genie::INukeHadroData::ReadhNFile(), readMyTextFile(), stan::test::performance::run_command(), slid::ShowerLID::SetCosTheta(), caldp::DCMSummary::SetDCMOffsetCount(), caldp::DCMSummary::SetDCMOffsetSum(), caldp::DCMSummary::SetDCMOffsetSumSquare(), caldp::DCMSummary::SetDCMUse(), slid::ShowerPID::SetDedx0(), slid::ShowerLID::SetDedx0(), jmshower::EID::SetDedx0(), slid::ShowerPID::SetDedx1(), slid::ShowerLID::SetDedx1(), jmshower::EID::SetDedx1(), slid::ShowerPID::SetDedx2(), slid::ShowerLID::SetDedx2(), jmshower::EID::SetDedx2(), slid::ShowerPID::SetDedx3(), slid::ShowerLID::SetDedx3(), jmshower::EID::SetDedx3(), slid::ShowerPID::SetDedx4(), slid::ShowerLID::SetDedx4(), jmshower::EID::SetDedx4(), slid::ShowerPID::SetDedx5(), slid::ShowerLID::SetDedx5(), jmshower::EID::SetDedx5(), slid::ShowerLID::SetDedxp0c0(), slid::ShowerLID::SetDedxp0c1(), slid::ShowerLID::SetDedxp1c0(), slid::ShowerLID::SetDedxp1c1(), slid::ShowerLID::SetDedxp2c0(), slid::ShowerLID::SetDedxp2c1(), slid::ShowerLID::SetDedxp3c0(), slid::ShowerLID::SetDedxp3c1(), slid::ShowerLID::SetDedxp4c0(), slid::ShowerLID::SetDedxp4c1(), slid::ShowerLID::SetDedxp5c0(), slid::ShowerLID::SetDedxp5c1(), slid::ShowerLID::SetEGLLL(), slid::ShowerLID::SetEGLLT(), slid::ShowerPID::SetELLL(), slid::ShowerPID::SetELLT(), slid::ShowerLID::SetEMuLLL(), slid::ShowerLID::SetEMuLLT(), slid::ShowerLID::SetENLLL(), slid::ShowerLID::SetENLLT(), slid::ShowerLID::SetEPi0LLL(), slid::ShowerLID::SetEPi0LLT(), slid::ShowerLID::SetEPiLLL(), slid::ShowerLID::SetEPiLLT(), slid::ShowerLID::SetEPLLL(), slid::ShowerLID::SetEPLLT(), slid::ShowerLID::SetGap(), slid::ShowerPID::SetGLLL(), slid::ShowerPID::SetGLLT(), slid::ShowerLID::SetHadronicEnergy(), slid::ShowerLID::SetInverseGLLL(), slid::LID::SetIsMuon(), slid::ShowerLID::SetIsMuon(), upmuana::LlrUpmu::setLLR(), slid::ShowerPID::SetMuLLL(), slid::ShowerPID::SetMuLLT(), slid::ShowerPID::SetNLLL(), slid::ShowerPID::SetNLLT(), slid::ShowerLID::SetNMIPPlanes(), slid::LID::SetNueEnergy(), slid::ShowerLID::SetNueEnergy(), slid::ShowerPID::SetPi0LLL(), slid::ShowerPID::SetPi0LLT(), slid::ShowerLID::SetPi0mass(), slid::ShowerLID::SetPi0PhotonID(), slid::ShowerPID::SetPiLLL(), slid::ShowerPID::SetPiLLT(), slid::ShowerPID::SetPLLL(), slid::ShowerPID::SetPLLT(), slid::ShowerLID::SetRadius(), slid::ShowerLID::SetShowerDepEnergy(), slid::ShowerLID::SetShowerEFrac(), slid::ShowerLID::SetShowerEnergy(), slid::LID::SetShwDepEnergy(), slid::LID::SetShwEnergy(), advi_test::SetUp(), slid::EventLID::SetValECos(), slid::ShowerLID::SetValECos(), slid::EventLID::SetValEPi0(), slid::ShowerPID::SetValEPi0(), slid::ShowerLID::SetValEPi0(), slid::ShowerPID::SetValEPi0EL(), slid::ShowerLID::SetValEPi0EL(), slid::SliceLID::SetValue(), slid::EventLID::SetValWithE(), slid::ShowerLID::SetValWithE(), slid::LID::SetVertexEnergy(), slid::ShowerLID::SetVertexEnergy(), slid::ShowerLID::SetVtxDoca(), stanc_helper(), TB_Nearline_Watchdog_Plots(), TEST(), test_exception(), TEST_F(), test_list(), test_val(), testfom(), simb::MCNeutrino::Theta(), trim(), genie::KPhaseSpace::UseInteraction(), art::detail::verifyPtrCollection(), and HEPREP::HepRepFactory::~HepRepFactory().

leg = new TLegend(0.50,0.60,0.20,0.70)

Definition at line 61 of file comparison.C.

Int_t nentries = (Int_t)experimentalTree -> GetEntries()

Definition at line 29 of file comparison.C.

TNtuple* ntupleExperimental = new TNtuple("ntupleExperimental","Protons, exp. data", "depthExp:EdepExp")

Definition at line 28 of file comparison.C.

TH1D* simulatedPeak = (TH1D*) simulationFile -> Get("braggPeak")

Definition at line 45 of file comparison.C.

TFile* simulationFile = new TFile("../../../SimulationOutputs/proton/BraggPeak/protonBraggPeak.root","READ")

Definition at line 42 of file comparison.C.

Float_t simulationNormalisationFactor = simulatedPeak -> GetBinContent(1)

Definition at line 47 of file comparison.C.