Functions | Variables
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.

31  {
33  Float_t normFactor = EdepExp;
35  ntupleExperimental -> Fill(depthExp, EdepExp/normFactor);
36 
37  }
Float_t normFactor
Float_t depthExp
Definition: comparison.C:22
Float_t EdepExp
Definition: comparison.C:22
correl_yv Fill(-(dy[iP-1][iC-1]), hyv->GetBinContent(iP, iC))
TTree * experimentalTree
Definition: comparison.C:20
TNtuple * ntupleExperimental
Definition: comparison.C:28
ntupleExperimental GetEntry(0)
simulatedPeak Scale ( 1/  simulationNormalisationFactor)
experimentalTree SetBranchAddress ( "EdepExp"  ,
EdepExp 
)
experimentalTree SetBranchAddress ( "depthExp"  ,
depthExp 
)
leg SetFillColor ( )
h2 SetMarkerSize ( )
ntupleExperimental SetMarkerStyle ( )
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.