Collaboration diagram for HistFamily:

Public Member Functions
	HistFamily (const string &)
void	fill (const xAOD::TrackParticle &, const xAOD::Vertex *)
	HistFamily (const HistFamily &)=delete
HistFamily &	operator= (const HistFamily &)=delete

Private Member Functions
uint8_t	getSum (const xAOD::TrackParticle &, xAOD::SummaryType) const

Private Attributes
TH2 *	h_pt = nullptr
TH2 *	h_si_hits_phys = nullptr
TH2 *	h_si_hits = nullptr
TH2 *	h_si_holes = nullptr
TH2 *	h_si_shared = nullptr
TH2 *	h_ibl_hits = nullptr
TH2 *	h_bl_hits = nullptr
TH2 *	h_ibl_expected = nullptr
TH2 *	h_bl_expected = nullptr
TH2 *	h_pix_hits_phys = nullptr
TH2 *	h_pix_hits = nullptr
TH2 *	h_pix_holes = nullptr
TH2 *	h_pix_shared = nullptr
TH2 *	h_sct_hits_phys = nullptr
TH2 *	h_sct_hits = nullptr
TH2 *	h_sct_holes = nullptr
TH2 *	h_sct_shared = nullptr
TH2 *	h_trt_hits = nullptr
TH2 *	h_trt_outlier_fraction = nullptr
TH2 *	h_d0 = nullptr
TH2 *	h_z0sintheta = nullptr

Detailed Description

Definition at line 33 of file InDetTrackSelectionToolTester.cxx.

Constructor & Destructor Documentation

◆ HistFamily() [1/2]

HistFamily::HistFamily ( const string & cut_name )

Definition at line 187 of file InDetTrackSelectionToolTester.cxx.

{
  auto* currentFile = TFile::CurrentFile();
  currentFile->mkdir(cut_name.data())->cd(); // create a directory for this cut type and cd into it
#define HIST_INIT( NAME, AXIS_LABEL, AXIS_N, AXIS_XL, AXIS_XH )     \
  do{                                   \
    h_##NAME = new TH2F(#NAME, #NAME ";#eta;" AXIS_LABEL, 50,-2.5,2.5,  \
            AXIS_N, AXIS_XL, AXIS_XH);          \
  } while (false)
 
  HIST_INIT( pt, "p_{T}  [GeV]", 100, 0, 20 );
  HIST_INIT( si_hits_phys, "Si physical hits", 24, 0, 24 );
  HIST_INIT( si_hits, "Si hits", 24, 0, 24 );
  HIST_INIT( si_holes, "Si holes", 6, 0, 6 );
  HIST_INIT( si_shared, "Si shared hits", 6, 0, 6 );
  HIST_INIT( pix_hits_phys, "Pixel physical hits", 10, 0, 10 );
  HIST_INIT( pix_hits, "Pixel hits", 10, 0, 10 );
  HIST_INIT( pix_holes, "Pixel holes", 4, 0, 4 );
  HIST_INIT( pix_shared, "Pixel shared hits", 4, 0, 4 );
  HIST_INIT( sct_hits_phys, "SCT physical hits", 16, 0, 16 );
  HIST_INIT( sct_hits, "SCT hits", 16, 0, 16 );
  HIST_INIT( sct_holes, "SCT holes", 4, 0, 4 );
  HIST_INIT( sct_shared, "SCT shared hits", 4, 0, 4 );
  HIST_INIT( ibl_hits, "IBL hits", 4, 0, 4 );
  HIST_INIT( ibl_expected, "expect IBL hit", 2, 0, 2 );
  HIST_INIT( bl_hits, "BLayer hits", 4, 0, 4 );
  HIST_INIT( bl_expected, "expect BLayer hit", 2, 0, 2 );
  HIST_INIT( trt_hits, "TRT hits + outliers", 60, 0, 60 );
  HIST_INIT( trt_outlier_fraction, "TRT outlier fraction", 25, 0., 1. );
  
  h_d0 = new TH2F( "d0", "d_{0}^{BL};d_{0}^{BL}  [mm];#sigma_{d_{0}}  [mm]", 60, -3., 3., 60, 0., 3. );
  h_z0sintheta = new TH2F( "z0sintheta", "z_{0}^{PV} sin #theta;z_{0}^{PV} sin #theta  [mm];#sigma_{z_{0} sin #theta}  [mm]", 100, -5., 5., 50, 0., 5. );
 
#undef HIST_INIT
 
  currentFile->cd(); // move back to the parent directory
}

◆ HistFamily() [2/2]

HistFamily::HistFamily ( const HistFamily & )

delete

Member Function Documentation

◆ fill()

void HistFamily::fill	(	const xAOD::TrackParticle &	trk,
		const xAOD::Vertex *	vtx )

Definition at line 225 of file InDetTrackSelectionToolTester.cxx.

{
  uint8_t nPixHits = getSum(trk, xAOD::numberOfPixelHits);
  uint8_t nPixDead = getSum(trk, xAOD::numberOfPixelDeadSensors);
  uint8_t nPixShared = getSum(trk, xAOD::numberOfPixelSharedHits);
  uint8_t nPixHoles = getSum(trk, xAOD::numberOfPixelHoles);
  uint8_t nSctHits = getSum(trk, xAOD::numberOfSCTHits);
  uint8_t nSctDead = getSum(trk, xAOD::numberOfSCTDeadSensors);
  uint8_t nSctShared = getSum(trk, xAOD::numberOfSCTSharedHits);
  uint8_t nSctHoles = getSum(trk, xAOD::numberOfSCTHoles);
  uint8_t nIblHits = getSum(trk, xAOD::numberOfInnermostPixelLayerHits);
  uint8_t expectIbl = getSum(trk, xAOD::expectInnermostPixelLayerHit);
  uint8_t nBlHits = getSum(trk, xAOD::numberOfNextToInnermostPixelLayerHits);
  uint8_t expectBl = getSum(trk, xAOD::expectNextToInnermostPixelLayerHit);
  uint8_t nTrtHits = getSum(trk, xAOD::numberOfTRTHits);
  uint8_t nTrtOutliers = getSum(trk, xAOD::numberOfTRTOutliers);
 
  double eta = trk.eta();
  h_pt->Fill(eta, trk.pt()*1e-3);
  h_ibl_hits->Fill(eta, nIblHits);
  h_ibl_expected->Fill(eta, expectIbl);
  h_bl_hits->Fill(eta, nBlHits);
  h_bl_expected->Fill(eta, expectBl);
  h_pix_hits_phys->Fill(eta, nPixHits);
  h_pix_hits->Fill(eta, nPixHits + nPixDead);
  h_pix_holes->Fill(eta, nPixHoles);
  h_pix_shared->Fill(eta, nPixShared);
  h_sct_hits_phys->Fill(eta, nSctHits);
  h_sct_hits->Fill(eta, nSctHits + nSctDead);
  h_sct_holes->Fill(eta, nSctHoles);
  h_sct_shared->Fill(eta, nSctShared);
  h_si_hits_phys->Fill(eta, nPixHits + nSctHits);
  h_si_hits->Fill(eta, nPixHits + nPixDead + nSctHits + nSctDead);
  h_si_holes->Fill(eta, nPixHoles + nSctHoles);
  h_si_shared->Fill(eta, nPixShared + nSctShared);
  h_trt_hits->Fill(eta, nTrtHits + nTrtOutliers);
  h_trt_outlier_fraction->Fill(eta, nTrtOutliers / static_cast<double>(nTrtHits + nTrtOutliers));
 
  const auto& covMatrix = trk.definingParametersCovMatrix();
  h_d0->Fill(trk.d0(), std::sqrt(covMatrix(0,0)) );
  double z0 = trk.z0();
  if (vtx != nullptr) {
    z0 += trk.vz() - vtx->z();
  }
  double sinTheta = std::sin(trk.theta());
  double cosTheta = std::cos(trk.theta());
  double varZ0SinTheta = 0.;
  varZ0SinTheta += covMatrix(1,1)*sinTheta*sinTheta;
  varZ0SinTheta += 2*covMatrix(1,3)*z0*sinTheta*cosTheta;
  varZ0SinTheta += covMatrix(3,3)*z0*z0*cosTheta*cosTheta;
  h_z0sintheta->Fill(z0*sinTheta, std::sqrt(varZ0SinTheta));
 
  return;
}

◆ getSum()

uint8_t HistFamily::getSum	(	const xAOD::TrackParticle &	trk,
		xAOD::SummaryType	sum_type ) const

private

Definition at line 280 of file InDetTrackSelectionToolTester.cxx.

{
  uint8_t sum_val = 0;
    if (!trk.summaryValue(sum_val, sum_type)) {
    Error( "HistFamily::getSum()", "Could not get summary type %i", sum_type );
  }
  return sum_val;
}