Trk::ExtrRegionPlots Class Reference

#include <TruthTrkExtrapolationPlots.h>

Inheritance diagram for Trk::ExtrRegionPlots:

Collaboration diagram for Trk::ExtrRegionPlots:

Public Member Functions
	ExtrRegionPlots (PlotBase *pParent, const std::string &sDir, std::string sRegionBegin, std::string sRegionEnd)
void	fill (const xAOD::TruthParticle &truthprt, const std::string &sRegionBegin, const std::string &sRegionEnd, float weight=1.0)
void	initialize ()
void	finalize ()
void	setDetailLevel (int iDetailLevel)
void	RegisterSubPlot (PlotBase *pPlotBase)
std::vector< HistData >	retrieveBookedHistograms ()
	Retrieve all booked histograms.
std::vector< TreeData >	retrieveBookedTrees ()
	Retrieve all booked trees.
std::vector< EfficiencyData >	retrieveBookedEfficiencies ()
	Retrieve all booked efficiency objects.
TTree *	BookTree (const std::string &name, bool prependDir=true)
	Book a TTree.
const std::string &	getDirectory ()
Methods to book monitoring histograms
Note: methods starting with capitals should be deprecated in favour of camel-cased methods
TH1D *	Book1D (const std::string &name, const std::string &labels, int nBins, float start, float end, bool prependDir=true)
	Book a TH1D histogram.
TH1D *	Book1D (const std::string &name, TH1 *refHist, const std::string &labels, bool prependDir=true)
	Book a TH1D histogram using refHist as reference for number of bins and axis range.
TH2F *	Book2D (const std::string &name, const std::string &labels, int nBinsX, float startX, float endX, int nBinsY, float startY, float endY, bool prependDir=true)
	Book a TH2F histogram.
TH2F *	Book2D (const std::string &name, TH2 *refHist, const std::string &labels, bool prependDir=true)
	Book a TH2D histogram using refHist as reference for number of bins and axis range.
TH2F *	Book2D (const std::string &name, const std::string &labels, int nBinsX, Double_t *binsX, int nBinsY, Double_t startY, Double_t endY, bool prependDir=true)
	Book a TH2F histogram with variable x axis binning.
TH3F *	Book3D (const std::string &name, const std::string &labels, int nBinsX, float startX, float endX, int nBinsY, float startY, float endY, int nBinsZ, float startZ, float endZ, bool prependDir=true)
	Book a TH3F histogram.
TH3F *	Book3D (const std::string &name, TH3 *refHist, const std::string &labels, bool prependDir=true)
	Book a TH3F histogram using refHist as reference for number of bins and axis range.
TProfile *	BookTProfile (const std::string &name, const std::string &labels, int nBinsX, float startX, float endX, float startY=-1, float endY=-1, bool prependDir=true, bool useRMS=false)
	Book a TProfile histogram.
TProfile *	BookTProfile (const std::string &name, const std::string &labels, int nBinsX, float *binsX, bool prependDir=true)
	Book a TProfile histogram with variable binning in x-axis.
TProfile *	BookTProfileRangeY (const std::string &name, const std::string &labels, int nBinsX, double *binsX, double startY, double endY, bool prependDir=true)
	Book a TProfile histogram with variable binning in x-axis and limits in y-values.
TProfile2D *	BookTProfile2D (const std::string &name, const std::string &labels, const int nBinsX, const double xlo, const double xhi, const int nBinsY, const double ylo, const double yhi, bool prependDir=true, bool useRMS=false)
	Book a TProfile 2D histogram with variable binning in x-axis and limits in y-values.
TProfile2D *	BookTProfile2D (const std::string &name, const std::string &labels, const int nBinsX, double *binsX, const int nBinsY, double *binsY, bool prependDir=true, bool useRMS=false)
	Book a TProfile 2D histogram with variable binning in x-axis and limits in y-values.
TEfficiency *	BookTEfficiency (const std::string &name, const std::string &labels, const int nBinsX, const float xlo, const float xhi, const bool prependDir=true)
	Book a (1-D) TEfficiency histogram.
TEfficiency *	BookTEfficiency (const std::string &name, const std::string &labels, const int nBinsX, const float xlo, const float xhi, const int nBinsy, const float ylo, const float yhi, const bool prependDir=true)
	Book a (2-D) TEfficiency histogram.

Public Attributes
std::string	m_sRegionBegin
std::string	m_sRegionEnd
TH1 *	dp
TH1 *	dpt
TH1 *	dpx
TH1 *	dpy
TH1 *	dpz
TH1 *	dp_extr
TH1 *	dpx_extr
TH1 *	dpy_extr
TH1 *	dpz_extr
TH2 *	dp_vs_p
TH2 *	dp_vs_phi
TH2 *	dp_vs_eta
TH2 *	dp_vs_eta_phi
TH2 *	p_vs_p
TH2 *	p_extr_vs_p_extr
TH2 *	dp_extr_vs_eta
TH1 *	dR
TH2 *	dR_vs_p
TH2 *	dR_vs_dp
TH2 *	dR_vs_eta
TH2 *	dR_vs_phi
TH2 *	dR_vs_eta_phi
TH1 *	dphi
TH1 *	dtheta
TH1 *	dAngle
TH2 *	dp_CaloEntry_IP_dp_CaloEntry_etxr

Protected Attributes
std::vector< PlotBase * >	m_vSubNodes
std::vector< HistData >	m_vBookedHistograms
std::vector< TreeData >	m_vBookedTrees
std::vector< EfficiencyData >	m_vBookedEfficiencies
std::string	m_sDirectory
int	m_iDetailLevel

Private Member Functions
void	initializePlots ()
virtual void	finalizePlots ()

Static Private Member Functions
static std::string	constructPrefix (std::string dir, bool prependDir)

Detailed Description

Definition at line 42 of file TruthTrkExtrapolationPlots.h.

Constructor & Destructor Documentation

ExtrRegionPlots()

Trk::ExtrRegionPlots::ExtrRegionPlots	(	PlotBase *	pParent,
		const std::string &	sDir,
		std::string	sRegionBegin,
		std::string	sRegionEnd )

Definition at line 121 of file TruthTrkExtrapolationPlots.cxx.

                                                         : PlotBase(pParent, sDir),
    m_sRegionBegin(std::move(sRegionBegin)),
    m_sRegionEnd(std::move(sRegionEnd)),
    dp{},dpt{},dpx{},dpy{},dpz{},dp_extr{},dpx_extr{},dpy_extr{},dpz_extr{},
    dp_vs_p{},dp_vs_phi{},dp_vs_eta{},dp_vs_eta_phi{},p_vs_p{},p_extr_vs_p_extr{},dp_extr_vs_eta{},
    dR{},dR_vs_p{},dR_vs_dp{},dR_vs_eta{},dR_vs_phi{},dR_vs_eta_phi{},dphi{},dtheta{},dAngle{},
    dp_CaloEntry_IP_dp_CaloEntry_etxr{} {
    //nop
  }

Member Function Documentation

Book1D() [1/2]

TH1D * PlotBase::Book1D ( const std::string & name, const std::string & labels, int nBins, float start, float end, bool prependDir = true )

inherited

Book a TH1D histogram.

Definition at line 94 of file PlotBase.cxx.

                                  {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  Bool_t oldstat = TH1::AddDirectoryStatus();
  TH1::AddDirectory(false);
  TH1D *hist = new TH1D((prefix + name).c_str(), labels.c_str(), nBins, start, end);
  TH1::AddDirectory(oldstat);
 
  hist->Sumw2();
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

Book1D() [2/2]

TH1D * PlotBase::Book1D ( const std::string & name, TH1 * refHist, const std::string & labels, bool prependDir = true )

inherited

Book a TH1D histogram using refHist as reference for number of bins and axis range.

Definition at line 108 of file PlotBase.cxx.

                                                                                              {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  Bool_t oldstat = TH1::AddDirectoryStatus();
  TH1::AddDirectory(false);
  TH1D *hist = new TH1D((prefix + name).c_str(), labels.c_str(), refHist->GetNbinsX(),
                        refHist->GetXaxis()->GetXbins()->GetArray());
  hist->Sumw2();
  TH1::AddDirectory(oldstat);
 
 
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

Book2D() [1/3]

TH2F * PlotBase::Book2D ( const std::string & name, const std::string & labels, int nBinsX, Double_t * binsX, int nBinsY, Double_t startY, Double_t endY, bool prependDir = true )

inherited

Book a TH2F histogram with variable x axis binning.

Definition at line 144 of file PlotBase.cxx.

                                                                  {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  Bool_t oldstat = TH2::AddDirectoryStatus();
  TH2::AddDirectory(false);
  TH2F *hist = new TH2F((prefix + name).c_str(), labels.c_str(), nBinsX, binsX, nBinsY, startY, endY);
  hist->Sumw2();
  TH2::AddDirectory(oldstat);
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

Book2D() [2/3]

TH2F * PlotBase::Book2D ( const std::string & name, const std::string & labels, int nBinsX, float startX, float endX, int nBinsY, float startY, float endY, bool prependDir = true )

inherited

Book a TH2F histogram.

Definition at line 123 of file PlotBase.cxx.

                                                            {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  Bool_t oldstat = TH2::AddDirectoryStatus();
  TH2::AddDirectory(false);
  TH2F *hist = new TH2F((prefix + name).c_str(), labels.c_str(), nBinsX, startX, endX, nBinsY, startY, endY);
  hist->Sumw2();
  TH2::AddDirectory(oldstat);
 
 
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

Book2D() [3/3]

TH2F * PlotBase::Book2D ( const std::string & name, TH2 * refHist, const std::string & labels, bool prependDir = true )

inherited

Book a TH2D histogram using refHist as reference for number of bins and axis range.

Definition at line 138 of file PlotBase.cxx.

                                                                                              {
  return Book2D(name, labels, refHist->GetNbinsX(), refHist->GetXaxis()->GetXmin(), refHist->GetXaxis()->GetXmax(),
                refHist->GetNbinsY(), refHist->GetYaxis()->GetXmin(), refHist->GetYaxis()->GetXmax(), prependDir);
}

Book3D() [1/2]

TH3F * PlotBase::Book3D ( const std::string & name, const std::string & labels, int nBinsX, float startX, float endX, int nBinsY, float startY, float endY, int nBinsZ, float startZ, float endZ, bool prependDir = true )

inherited

Book a TH3F histogram.

Definition at line 157 of file PlotBase.cxx.

                                                                                                  {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  Bool_t oldstat = TH3::AddDirectoryStatus();
  TH3::AddDirectory(false);
  TH3F *hist = new TH3F((prefix + name).c_str(),
                        labels.c_str(), nBinsX, startX, endX, nBinsY, startY, endY, nBinsZ, startZ, endZ);
  hist->Sumw2();
  TH3::AddDirectory(oldstat);
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

Book3D() [2/2]

TH3F * PlotBase::Book3D ( const std::string & name, TH3 * refHist, const std::string & labels, bool prependDir = true )

inherited

Book a TH3F histogram using refHist as reference for number of bins and axis range.

Definition at line 171 of file PlotBase.cxx.

                                                                                              {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  Bool_t oldstat = TH3::AddDirectoryStatus();
  TH3::AddDirectory(false);
  TH3F *hist = new TH3F((prefix + name).c_str(), labels.c_str(), refHist->GetNbinsX(),
                        refHist->GetXaxis()->GetXbins()->GetArray(), refHist->GetNbinsY(),
                        refHist->GetYaxis()->GetXbins()->GetArray(), refHist->GetNbinsZ(),
                        refHist->GetZaxis()->GetXbins()->GetArray());
  TH3::AddDirectory(oldstat);
 
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

BookTEfficiency() [1/2]

TEfficiency * PlotBase::BookTEfficiency ( const std::string & name, const std::string & labels, const int nBinsX, const float xlo, const float xhi, const bool prependDir = true )

inherited

Book a (1-D) TEfficiency histogram.

Definition at line 257 of file PlotBase.cxx.

                                                                                                                                                   {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  //Bool_t oldstat = TEfficiency::AddDirectoryStatus();
  TEfficiency *hist = new TEfficiency((prefix + name).c_str(), labels.c_str(), nBinsX, xlo, xhi);
  //hist->SetAutoSave(0);
  //hist->SetAtoFlush(0);
  hist->SetDirectory(nullptr);
  m_vBookedEfficiencies.emplace_back(hist, m_sDirectory);
  //TEfficiency::AddDirectory(oldstat);
  return hist;
}

BookTEfficiency() [2/2]

TEfficiency * PlotBase::BookTEfficiency ( const std::string & name, const std::string & labels, const int nBinsX, const float xlo, const float xhi, const int nBinsy, const float ylo, const float yhi, const bool prependDir = true )

inherited

Book a (2-D) TEfficiency histogram.

Definition at line 270 of file PlotBase.cxx.

                                                                                                                                                                                                       {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
 
  TEfficiency *hist = new TEfficiency((prefix + name).c_str(), labels.c_str(), nBinsX, xlo, xhi, nBinsY, ylo, yhi);
  hist->SetDirectory(nullptr);
  m_vBookedEfficiencies.emplace_back(hist, m_sDirectory);
 
  return hist;
}

BookTProfile() [1/2]

TProfile * PlotBase::BookTProfile ( const std::string & name, const std::string & labels, int nBinsX, float * binsX, bool prependDir = true )

inherited

Book a TProfile histogram with variable binning in x-axis.

Definition at line 204 of file PlotBase.cxx.

                                                                                                                {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  TProfile *hist(nullptr);
  Bool_t oldstat = TProfile::AddDirectoryStatus();
  TProfile::AddDirectory(false);
 
  hist = new TProfile((prefix + name).c_str(), labels.c_str(), nBinsX, binsX);
  TProfile::AddDirectory(oldstat);
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

BookTProfile() [2/2]

TProfile * PlotBase::BookTProfile ( const std::string & name, const std::string & labels, int nBinsX, float startX, float endX, float startY = -1, float endY = -1, bool prependDir = true, bool useRMS = false )

inherited

Book a TProfile histogram.

Definition at line 186 of file PlotBase.cxx.

                                                                               {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  TProfile *hist(nullptr);
  Bool_t oldstat = TProfile::AddDirectoryStatus();
  TProfile::AddDirectory(false);
  std::string opt = useRMS ? "S" : "";
  if ((startY == -1) and (endY == -1)) {
    hist = new TProfile((prefix + name).c_str(), labels.c_str(), nBinsX, startX, endX, opt.c_str());
  } else {
    hist = new TProfile((prefix + name).c_str(), labels.c_str(), nBinsX, startX, endX, startY, endY, opt.c_str());
  }
  TProfile::AddDirectory(oldstat);
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

BookTProfile2D() [1/2]

TProfile2D * PlotBase::BookTProfile2D ( const std::string & name, const std::string & labels, const int nBinsX, const double xlo, const double xhi, const int nBinsY, const double ylo, const double yhi, bool prependDir = true, bool useRMS = false )

inherited

Book a TProfile 2D histogram with variable binning in x-axis and limits in y-values.

Definition at line 231 of file PlotBase.cxx.

                                                       {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  Bool_t oldstat = TProfile2D::AddDirectoryStatus();
  TProfile2D::AddDirectory(false);
  std::string opt = useRMS ? "S" : "";
  TProfile2D *hist = new TProfile2D((prefix + name).c_str(), labels.c_str(), nBinsX, xlo, xhi, nBinsY, ylo, yhi, opt.c_str());
  TProfile2D::AddDirectory(oldstat);
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

BookTProfile2D() [2/2]

TProfile2D * PlotBase::BookTProfile2D ( const std::string & name, const std::string & labels, const int nBinsX, double * binsX, const int nBinsY, double * binsY, bool prependDir = true, bool useRMS = false )

inherited

Book a TProfile 2D histogram with variable binning in x-axis and limits in y-values.

Definition at line 245 of file PlotBase.cxx.

                                                                                                                                                                       {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  Bool_t oldstat = TProfile2D::AddDirectoryStatus();
  TProfile2D::AddDirectory(false);
  std::string opt = useRMS ? "S" : "";
  TProfile2D *hist = new TProfile2D((prefix + name).c_str(), labels.c_str(), nBinsX, binsX, nBinsY, binsY, opt.c_str());
  TProfile2D::AddDirectory(oldstat);
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

BookTProfileRangeY()

TProfile * PlotBase::BookTProfileRangeY ( const std::string & name, const std::string & labels, int nBinsX, double * binsX, double startY, double endY, bool prependDir = true )

inherited

Book a TProfile histogram with variable binning in x-axis and limits in y-values.

Definition at line 217 of file PlotBase.cxx.

                                                                          {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  TProfile *hist(nullptr);
  Bool_t oldstat = TProfile::AddDirectoryStatus();
  TProfile::AddDirectory(false);
 
  hist = new TProfile((prefix + name).c_str(), labels.c_str(), (Int_t) nBinsX, binsX, startY, endY);
  TProfile::AddDirectory(oldstat);
  m_vBookedHistograms.emplace_back(hist, m_sDirectory);
  return hist;
}

BookTree()

TTree * PlotBase::BookTree ( const std::string & name, bool prependDir = true )

inherited

Book a TTree.

Definition at line 281 of file PlotBase.cxx.

                                                         {
  std::string prefix = constructPrefix(m_sDirectory, prependDir);
  TTree *tree = new TTree((prefix + name).c_str(), "");
 
  tree->SetAutoSave(0);
  tree->SetAutoFlush(0);
  tree->SetDirectory(nullptr);
  m_vBookedTrees.emplace_back(tree, m_sDirectory);
  return tree;
}

constructPrefix()

std::string PlotBase::constructPrefix ( std::string dir, bool prependDir )

staticprivateinherited

Definition at line 293 of file PlotBase.cxx.

                                                        {
  if (!prependDir) {
    return "";
  }
  std::replace(dir.begin(), dir.end(), '/', '_');
  return dir;
}

fill()

void Trk::ExtrRegionPlots::fill	(	const xAOD::TruthParticle &	truthprt,
		const std::string &	sRegionBegin,
		const std::string &	sRegionEnd,
		float	weight = 1.0 )

Definition at line 227 of file TruthTrkExtrapolationPlots.cxx.

                                                                                                                             {
    SG::ConstAccessor<float> pxBeginAcc (sDetBegin + "_px");
    SG::ConstAccessor<float> pyBeginAcc (sDetBegin + "_py");
    SG::ConstAccessor<float> pzBeginAcc (sDetBegin + "_pz");
    SG::ConstAccessor<float> pxEndAcc (sDetEnd + "_px");
    SG::ConstAccessor<float> pyEndAcc (sDetEnd + "_py");
    SG::ConstAccessor<float> pzEndAcc (sDetEnd + "_pz");
    if (!pxBeginAcc.isAvailable(truthprt) ||
        !pyBeginAcc.isAvailable(truthprt) ||
        !pzBeginAcc.isAvailable(truthprt) ||
        !pxEndAcc.isAvailable(truthprt) ||
        !pyEndAcc.isAvailable(truthprt) ||
        !pzEndAcc.isAvailable(truthprt)) {
      return;
    }
 
    TVector3 vecDetBegin(pxBeginAcc(truthprt), pyBeginAcc(truthprt), pzBeginAcc(truthprt));
    TVector3 vecDetEnd(pxEndAcc(truthprt), pyEndAcc(truthprt), pzEndAcc(truthprt));
 
 
    // check if IP, since extr is not available there (sDetBegin="")
    TVector3 vecDetBegin_extr;
    if (!sDetBegin.empty()) {
      SG::ConstAccessor<float> px_extrAcc (sDetBegin + "_px_extr");
      SG::ConstAccessor<float> py_extrAcc (sDetBegin + "_py_extr");
      SG::ConstAccessor<float> pz_extrAcc (sDetBegin + "_pz_extr");
      if (!px_extrAcc.isAvailable(truthprt) ||
          !py_extrAcc.isAvailable(truthprt) ||
          !pz_extrAcc.isAvailable(truthprt)) {
        return;
      }
 
      vecDetBegin_extr.SetXYZ(px_extrAcc(truthprt), py_extrAcc(truthprt), pz_extrAcc(truthprt));
    }else {
      vecDetBegin_extr = vecDetBegin;
    }
 
    TVector3 vecDetEnd_extr;
    if (!sDetEnd.empty()) {
      SG::ConstAccessor<float> px_extrAcc (sDetEnd + "_px_extr");
      SG::ConstAccessor<float> py_extrAcc (sDetEnd + "_py_extr");
      SG::ConstAccessor<float> pz_extrAcc (sDetEnd + "_pz_extr");
      if (!px_extrAcc.isAvailable(truthprt) ||
          !py_extrAcc.isAvailable(truthprt) ||
          !pz_extrAcc.isAvailable(truthprt)) {
        return;
      }
 
      vecDetEnd_extr.SetXYZ(px_extrAcc(truthprt), py_extrAcc(truthprt), pz_extrAcc(truthprt));
    }else {
      vecDetEnd_extr = vecDetEnd;
    }
 
    dp->Fill((vecDetBegin.Mag() - vecDetEnd.Mag()) * 0.001, weight);
    dpt->Fill((vecDetBegin.Perp() - vecDetEnd.Perp()) * 0.001, weight);
    dpx->Fill((vecDetBegin.Px() - vecDetEnd.Px()) * 0.001, weight);
    dpy->Fill((vecDetBegin.Py() - vecDetEnd.Py()) * 0.001, weight);
    dpz->Fill((vecDetBegin.Pz() - vecDetEnd.Pz()) * 0.001, weight);
    dp_extr->Fill((vecDetBegin_extr.Mag() - vecDetEnd_extr.Mag()) * 0.001, weight);
    dpx_extr->Fill((vecDetBegin_extr.Px() - vecDetEnd_extr.Px()) * 0.001, weight);
    dpy_extr->Fill((vecDetBegin_extr.Py() - vecDetEnd_extr.Py()) * 0.001, weight);
    dpz_extr->Fill((vecDetBegin_extr.Pz() - vecDetEnd_extr.Pz()) * 0.001, weight);
    // 2D
    dp_vs_p->Fill(vecDetBegin.Mag() * 0.001, (vecDetBegin.Mag() - vecDetEnd.Mag()) * 0.001, weight);
    dp_vs_phi->Fill(vecDetBegin.Phi(), (vecDetBegin.Mag() - vecDetEnd.Mag()) * 0.001, weight);
    dp_vs_eta->Fill(vecDetBegin.Eta(), (vecDetBegin.Mag() - vecDetEnd.Mag()) * 0.001, weight);
    // not sure how to weight this
    dp_vs_eta_phi->Fill(vecDetBegin.Eta(), vecDetBegin.Phi(), (vecDetBegin.Mag() - vecDetEnd.Mag()) * 0.001);
    p_vs_p->Fill(vecDetBegin.Mag() * 0.001, vecDetEnd.Mag() * 0.001, weight);
    p_extr_vs_p_extr->Fill(vecDetBegin_extr.Mag() * 0.001, vecDetEnd_extr.Mag() * 0.001, weight);
    dp_extr_vs_eta->Fill(vecDetBegin_extr.Eta(), (vecDetBegin_extr.Mag() - vecDetEnd_extr.Mag()) * 0.001, weight);
 
    // multiple scattering
    dR->Fill(vecDetBegin.DeltaR(vecDetEnd),weight);
    // 2D
    dR_vs_p->Fill(vecDetBegin.Mag() * 0.001, vecDetBegin.DeltaR(vecDetEnd),weight);
    dR_vs_dp->Fill((vecDetBegin.Mag() - vecDetEnd.Mag()) * 0.001, vecDetBegin.DeltaR(vecDetEnd), weight);
    dR_vs_eta->Fill(vecDetBegin.Eta(), vecDetBegin.DeltaR(vecDetEnd), weight);
    dR_vs_phi->Fill(vecDetBegin.Phi(), vecDetBegin.DeltaR(vecDetEnd), weight);
    // not sure how to weight this
    dR_vs_eta_phi->Fill(vecDetBegin.Eta(), vecDetBegin.Phi(), vecDetBegin.DeltaR(vecDetEnd) );
 
    dphi->Fill(vecDetBegin.DeltaPhi(vecDetEnd), weight);
    dtheta->Fill(vecDetBegin.Theta() - vecDetEnd.Theta(), weight);
 
    dAngle->Fill(vecDetBegin.Angle(vecDetEnd), weight);
  }

finalize()

void PlotBase::finalize ( )

inherited

Definition at line 47 of file PlotBase.cxx.

                   {
  for (auto *subNode: m_vSubNodes) {
    subNode->finalize();
  }
  finalizePlots();
}

finalizePlots()

virtual void PlotBase::finalizePlots ( )

inlineprivatevirtualinherited

Reimplemented in IDTPM::DuplicateRatePlots, IDTPM::EfficiencyPlots, IDTPM::FakeRatePlots, IDTPM::HitsOnTracksPlots, IDTPM::NtracksPlots, IDTPM::OfflineElectronPlots, IDTPM::ResolutionPlots, IDTPM::SummaryPlots, IDTPM::TrackParametersPlots, IDTPM::VertexParametersPlots, InDetBasicPlot, InDetPerfPlot_Duplicate, InDetPerfPlot_Efficiency, InDetPerfPlot_FakeRate, InDetPerfPlot_Resolution, InDetPerfPlot_TrackParameters, InDetPerfPlot_TRTExtension, InDetPerfPlot_VertexTruthMatching, JetTagDQA::BTaggingValidationPlots, TCCPlots, Trk::EfficiencyPlots, ZeeValidation::ReconElectronsPlots, ZeeValidation::TrueElectronsPlots, and ZeeValidation::TrueFwdElectronsPlots.

Definition at line 92 of file PlotBase.h.

{;}

getDirectory()

const std::string & PlotBase::getDirectory ( )

inlineinherited

Definition at line 88 of file PlotBase.h.

{return m_sDirectory;}

initialize()

void PlotBase::initialize ( )

inherited

Definition at line 39 of file PlotBase.cxx.

                     {
  for (auto *subNode: m_vSubNodes) {
    subNode->initialize();
  }
  initializePlots();
}

initializePlots()

void Trk::ExtrRegionPlots::initializePlots ( )

privatevirtual

Reimplemented from PlotBase.

Definition at line 133 of file TruthTrkExtrapolationPlots.cxx.

                                   {
    // histogram definitions for scattering plots
    // variables within subdetector(s)
    // energy loss related quantities
    dp = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dp",
                m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth p: p_{begin} vs p_{end}; #Delta p=p_{begin}-p_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                20.);
    dpt = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dpt",
                 m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth pt;#Delta pt=pt_{begin}-pt_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                 20.);
    dpx = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dpx",
                 m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth px;#Delta px=px_{begin}-px_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                 20.);
    dpy = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dpy",
                 m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth py;#Delta py=py_{begin}-py_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                 20.);
    dpz = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dpz",
                 m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth pz;#Delta pz=pz_{begin}-pz_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                 20.);
 
    // extrapolated quantities: CaloEntryLayer_p_extr=extrapolated value from IP to CaloEntryLayer using e-loss
    // parametrization
    // if comparison to IP, truth_p_extr=truth_p
    dp_extr = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dp_extr",
                     m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth p, extrapolated from IP to current layer; #Delta p=p_extr_{begin}-p_extr_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                     20.);
    dpx_extr = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dpx_extr",
                      m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth px, extrapolated from IP to current layer; #Delta px=px_extr_{begin}-px_extr_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                      20.);
    dpy_extr = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dpy_extr",
                      m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth py, extrapolated from IP to current layer; #Delta py=py_extr_{begin}-py_extr_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                      20.);
    dpz_extr = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dpz_extr",
                      m_sRegionBegin + " vs " + m_sRegionEnd + " difference in truth pz, extrapolated from IP to current layer; #Delta pz=pz_extr_{begin}-pz_extr_{end} [GeV];Entries / 0.5 GeV", 80, -20.,
                      20.);
 
    // 2D plots
    dp_vs_p = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dp_vs_p",
                     m_sRegionBegin + " vs " + m_sRegionEnd + " difference in p (#Delta p=p_{begin}-p_{end}) vs p_{begin};p_{begin} [GeV];#Delta p=p_{begin}-p_{end} [GeV]", 250, 0., 250., 80, 0.,
                     20);
    dp_vs_phi = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dp_vs_phi",
                       m_sRegionBegin + " vs " + m_sRegionEnd + " difference in p(#Delta p=p_{begin}-p_{end}) vs phi_{begin};phi_{begin};#Delta p=p_{begin}-p_{end} [GeV]", 128, -3.2, 3.2, 80, 0.,
                       20);
    dp_vs_eta = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dp_vs_eta",
                       m_sRegionBegin + " vs " + m_sRegionEnd + " difference in p(#Delta p=p_{begin}-p_{end}) vs eta_{begin};#eta_{begin};#Delta p=p_{begin}-p_{end} [GeV]", 128, -3.2, 3.2, 80, 0.,
                       20);
    p_vs_p = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_p_vs_p",
                    m_sRegionBegin + " vs " + m_sRegionEnd + " p_{begin} vs p_{end};p_{begin} [GeV]; p_{end} [GeV]",
                    250, 0., 250, 250, 0., 250);
    p_extr_vs_p_extr = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_p_extr_vs_p_extr",
                              m_sRegionBegin + " vs " + m_sRegionEnd + " p_{begin extr} vs p_{end extr};p_{begin extr} [GeV]; p_{end extr} [GeV]", 250, 0., 250, 250, 0.,
                              250);
    dp_extr_vs_eta = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dp_extr_vs_eta",
                            m_sRegionBegin + " vs " + m_sRegionEnd + " difference in p_{extr}(#Delta p_{extr}=p_{begin extr}-p_{end extr}) vs eta_{begin};#eta__{begin};#Delta p_{extr}=p_{begin extr}-p_{end extr} [GeV]", 128, -3.2, 3.2, 80, 0.,
                            20);
 
    // 3D
    dp_vs_eta_phi = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dp_vs_eta_phi",
                           m_sRegionBegin + " vs " + m_sRegionEnd + " sum of bin entries: difference in p(#Delta p=p_{begin}-p_{end}) vs eta_{begin} and phi_{begin};#eta__{begin};#phi__{begin};#Delta p=p_{begin}-p_{end} [GeV]", 128, -3.2, 3.2, 128, -3.2,
                           -3.2);
 
    // measure for multiple scattering
    dR = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dR",
                m_sRegionBegin + " vs " + m_sRegionEnd + " dR between the two positions in the detector;pt_{begin}->dR(pt_{end});Entries / 0.05", 100, 0.,
                0.5);
    // 2D
    dR_vs_p = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dR_vs_p",
                     m_sRegionBegin + " vs " + m_sRegionEnd + " dR vs p_{begin};p_{begin} [GeV];_{begin}->dR(p_{end})",
                     250, 0., 250, 100, 0., 0.5);
    dR_vs_dp = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dR_vs_dp",
                      m_sRegionBegin + " vs " + m_sRegionEnd + " dR vs difference in p;#Delta p=p_{begin}-p_{end} [GeV];p_{begin}->dR(p_{end})", 80, 0., 20, 100, 0.,
                      0.5);
    dR_vs_eta = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dR_vs_eta",
                       m_sRegionBegin + " vs " + m_sRegionEnd + " dR vs eta;#eta_{begin};p_{begin}->dR(p_{end})", 128,
                       -3.2, 3.2, 100, 0., 0.5);
    dR_vs_phi = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dR_vs_phi",
                       m_sRegionBegin + " vs " + m_sRegionEnd + " dR vs phi;#phi_{begin};p_{begin}->dR(p_{end})", 128,
                       -3.2, 3.2, 100, 0., 0.5);
    dR_vs_eta_phi = Book2D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dR_vs_eta_phi",
                           m_sRegionBegin + " vs " + m_sRegionEnd + " sum of bin entries: dR vs eta and phi;#eta_{begin};#phi_{begin};#DeltaR", 128, -3.2, 3.2, 128, -3.2,
                           3.2);
 
    dphi = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dphi",
                  m_sRegionBegin + " vs " + m_sRegionEnd + " difference in phi;#Delta phi;Entries / 0.05", 100, -0.25,
                  0.25);
    dtheta = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dtheta",
                    m_sRegionBegin + " vs " + m_sRegionEnd + " difference in theta;#Delta theta;Entries / 0.002", 100,
                    -0.25, 0.25);
    dAngle = Book1D(m_sRegionBegin + "_vs_" + m_sRegionEnd + "_dAngle",
                    m_sRegionBegin + " vs " + m_sRegionEnd + " total angle difference;#Delta angle;Entries / 0.002",
                    100, -0.25, 0.25);
  }

RegisterSubPlot()

void PlotBase::RegisterSubPlot ( PlotBase * pPlotBase )

inlineinherited

Definition at line 41 of file PlotBase.h.

{m_vSubNodes.push_back(pPlotBase);}

retrieveBookedEfficiencies()

std::vector< EfficiencyData > PlotBase::retrieveBookedEfficiencies ( )

inherited

Retrieve all booked efficiency objects.

Definition at line 83 of file PlotBase.cxx.

                                     {
  std::vector<EfficiencyData> vBookedEfficiencies = m_vBookedEfficiencies;
  for (const auto &subNode: m_vSubNodes) {
    std::vector<EfficiencyData> subNodeHists = subNode->retrieveBookedEfficiencies();
    vBookedEfficiencies.insert(vBookedEfficiencies.end(), subNodeHists.begin(), subNodeHists.end());
  }
  return vBookedEfficiencies;
}

retrieveBookedHistograms()

std::vector< HistData > PlotBase::retrieveBookedHistograms ( )

inherited

Retrieve all booked histograms.

Definition at line 63 of file PlotBase.cxx.

                                   {
  std::vector<HistData> vBookedHistograms = m_vBookedHistograms;
  for (const auto &subNode: m_vSubNodes) {
    std::vector<HistData> subNodeHists = subNode->retrieveBookedHistograms();
    vBookedHistograms.insert(vBookedHistograms.end(), subNodeHists.begin(), subNodeHists.end());
  }
  return vBookedHistograms;
}

retrieveBookedTrees()

std::vector< TreeData > PlotBase::retrieveBookedTrees ( )

inherited

Retrieve all booked trees.

Definition at line 73 of file PlotBase.cxx.

                              {
  std::vector<TreeData> vBookedTrees = m_vBookedTrees;
  for (auto *subNode: m_vSubNodes) {
    std::vector<TreeData> subNodeTrees = subNode->retrieveBookedTrees();
    vBookedTrees.insert(vBookedTrees.end(), subNodeTrees.begin(), subNodeTrees.end());
  }
  return vBookedTrees;
}

setDetailLevel()

void PlotBase::setDetailLevel ( int iDetailLevel )

inherited

Definition at line 55 of file PlotBase.cxx.

                                         {
  for (auto *subNode: m_vSubNodes) {
    subNode->setDetailLevel(iDetailLevel);
  }
  m_iDetailLevel = iDetailLevel;
}

Member Data Documentation

dAngle

TH1* Trk::ExtrRegionPlots::dAngle

Definition at line 78 of file TruthTrkExtrapolationPlots.h.

dp

TH1* Trk::ExtrRegionPlots::dp

Definition at line 50 of file TruthTrkExtrapolationPlots.h.

dp_CaloEntry_IP_dp_CaloEntry_etxr

TH2* Trk::ExtrRegionPlots::dp_CaloEntry_IP_dp_CaloEntry_etxr

Definition at line 80 of file TruthTrkExtrapolationPlots.h.

dp_extr

TH1* Trk::ExtrRegionPlots::dp_extr

Definition at line 55 of file TruthTrkExtrapolationPlots.h.

dp_extr_vs_eta

TH2* Trk::ExtrRegionPlots::dp_extr_vs_eta

Definition at line 66 of file TruthTrkExtrapolationPlots.h.

dp_vs_eta

TH2* Trk::ExtrRegionPlots::dp_vs_eta

Definition at line 62 of file TruthTrkExtrapolationPlots.h.

dp_vs_eta_phi

TH2* Trk::ExtrRegionPlots::dp_vs_eta_phi

Definition at line 63 of file TruthTrkExtrapolationPlots.h.

dp_vs_p

TH2* Trk::ExtrRegionPlots::dp_vs_p

Definition at line 60 of file TruthTrkExtrapolationPlots.h.

dp_vs_phi

TH2* Trk::ExtrRegionPlots::dp_vs_phi

Definition at line 61 of file TruthTrkExtrapolationPlots.h.

dphi

TH1* Trk::ExtrRegionPlots::dphi

Definition at line 76 of file TruthTrkExtrapolationPlots.h.

dpt

TH1* Trk::ExtrRegionPlots::dpt

Definition at line 51 of file TruthTrkExtrapolationPlots.h.

dpx

TH1* Trk::ExtrRegionPlots::dpx

Definition at line 52 of file TruthTrkExtrapolationPlots.h.

dpx_extr

TH1* Trk::ExtrRegionPlots::dpx_extr

Definition at line 56 of file TruthTrkExtrapolationPlots.h.

dpy

TH1* Trk::ExtrRegionPlots::dpy

Definition at line 53 of file TruthTrkExtrapolationPlots.h.

dpy_extr

TH1* Trk::ExtrRegionPlots::dpy_extr

Definition at line 57 of file TruthTrkExtrapolationPlots.h.

dpz

TH1* Trk::ExtrRegionPlots::dpz

Definition at line 54 of file TruthTrkExtrapolationPlots.h.

dpz_extr

TH1* Trk::ExtrRegionPlots::dpz_extr

Definition at line 58 of file TruthTrkExtrapolationPlots.h.

dR

TH1* Trk::ExtrRegionPlots::dR

Definition at line 69 of file TruthTrkExtrapolationPlots.h.

dR_vs_dp

TH2* Trk::ExtrRegionPlots::dR_vs_dp

Definition at line 71 of file TruthTrkExtrapolationPlots.h.

dR_vs_eta

TH2* Trk::ExtrRegionPlots::dR_vs_eta

Definition at line 72 of file TruthTrkExtrapolationPlots.h.

dR_vs_eta_phi

TH2* Trk::ExtrRegionPlots::dR_vs_eta_phi

Definition at line 74 of file TruthTrkExtrapolationPlots.h.

dR_vs_p

TH2* Trk::ExtrRegionPlots::dR_vs_p

Definition at line 70 of file TruthTrkExtrapolationPlots.h.

dR_vs_phi

TH2* Trk::ExtrRegionPlots::dR_vs_phi

Definition at line 73 of file TruthTrkExtrapolationPlots.h.

dtheta

TH1* Trk::ExtrRegionPlots::dtheta

Definition at line 77 of file TruthTrkExtrapolationPlots.h.

m_iDetailLevel

int PlotBase::m_iDetailLevel

protectedinherited

Definition at line 101 of file PlotBase.h.

m_sDirectory

std::string PlotBase::m_sDirectory

protectedinherited

Definition at line 100 of file PlotBase.h.

m_sRegionBegin

std::string Trk::ExtrRegionPlots::m_sRegionBegin

Definition at line 47 of file TruthTrkExtrapolationPlots.h.

m_sRegionEnd

std::string Trk::ExtrRegionPlots::m_sRegionEnd

Definition at line 48 of file TruthTrkExtrapolationPlots.h.

m_vBookedEfficiencies

std::vector<EfficiencyData> PlotBase::m_vBookedEfficiencies

protectedinherited

Definition at line 99 of file PlotBase.h.

m_vBookedHistograms

std::vector<HistData> PlotBase::m_vBookedHistograms

protectedinherited

Definition at line 97 of file PlotBase.h.

m_vBookedTrees

std::vector<TreeData> PlotBase::m_vBookedTrees

protectedinherited

Definition at line 98 of file PlotBase.h.

m_vSubNodes

std::vector<PlotBase*> PlotBase::m_vSubNodes

protectedinherited

Definition at line 96 of file PlotBase.h.

p_extr_vs_p_extr

TH2* Trk::ExtrRegionPlots::p_extr_vs_p_extr

Definition at line 65 of file TruthTrkExtrapolationPlots.h.

p_vs_p

TH2* Trk::ExtrRegionPlots::p_vs_p

Definition at line 64 of file TruthTrkExtrapolationPlots.h.

---

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

• TruthTrkExtrapolationPlots.h
• TruthTrkExtrapolationPlots.cxx