PVValHelper Namespace Reference

Classes
struct	histodetails

Typedefs
typedef std::tuple< std::string, std::string, std::string >	plotLabels

Enumerations
enum	detectorPhase { phase0 = 0, phase1 = 1, phase2 = 2, END_OF_PHASES = 3 }
enum	estimator {   MEAN = 1, WIDTH = 2, MEDIAN = 3, MAD = 4,   UNKWN = -1 }
enum	plotVariable {   phi = 1, eta = 2, pT = 3, pTCentral = 4,   ladder = 5, modZ = 6, END_OF_PLOTS = 7 }
enum	residualType {   dxy = 1, dx = 2, dy = 3, dz = 4,   IP2D = 5, resz = 6, IP3D = 7, d3D = 8,   norm_dxy = 9, norm_dx = 10, norm_dy = 11, norm_dz = 12,   norm_IP2D = 13, norm_resz = 14, norm_IP3D = 15, norm_d3D = 16,   END_OF_TYPES = 17 }

Functions
void	add (std::map< std::string, TH1 *> &h, TH1 *hist)
void	fill (std::map< std::string, TH1 *> &h, const std::string &s, double x)
void	fill (std::map< std::string, TH1 *> &h, const std::string &s, double x, double y)
void	fillByIndex (std::vector< TH1F *> &h, unsigned int index, double x, std::string tag="")
std::pair< Measurement1D, Measurement1D >	fitResiduals (TH1 *hist)
std::vector< float >	generateBins (int n, float start, float range)
Measurement1D	getMAD (TH1F *histo)
Measurement1D	getMedian (TH1F *histo)
plotLabels	getTypeString (residualType type)
plotLabels	getVarString (plotVariable var)
template<typename T , size_t N>
std::array< T, N+1 >	makeLogBins (const T &min, const T &max)
void	shrinkHistVectorToFit (std::vector< TH1F *> &h, unsigned int desired_size)

Variables
constexpr double	max_eta_phase0 = 2.5
constexpr double	max_eta_phase1 = 2.7
constexpr double	max_eta_phase2 = 4.0

Typedef Documentation

plotLabels

typedef std::tuple<std::string, std::string, std::string> PVValHelper::plotLabels

Definition at line 123 of file PVValidationHelpers.h.

Enumeration Type Documentation

detectorPhase

enum PVValHelper::detectorPhase

Enumerator
phase0
phase1
phase2
END_OF_PHASES

Definition at line 78 of file PVValidationHelpers.h.

                     {
    phase0 = 0,
    phase1 = 1,
    phase2 = 2,
    END_OF_PHASES = 3,
  };

estimator

enum PVValHelper::estimator

Enumerator
MEAN
WIDTH
MEDIAN
MAD
UNKWN

Definition at line 45 of file PVValidationHelpers.h.

{ MEAN = 1, WIDTH = 2, MEDIAN = 3, MAD = 4, UNKWN = -1 };

plotVariable

enum PVValHelper::plotVariable

Enumerator
phi
eta
pT
pTCentral
ladder
modZ
END_OF_PLOTS

Definition at line 68 of file PVValidationHelpers.h.

                    {
    phi = 1,
    eta = 2,
    pT = 3,
    pTCentral = 4,
    ladder = 5,
    modZ = 6,
    END_OF_PLOTS = 7,
  };

residualType

enum PVValHelper::residualType

Enumerator
dxy
dx
dy
dz
IP2D
resz
IP3D
d3D
norm_dxy
norm_dx
norm_dy
norm_dz
norm_IP2D
norm_resz
norm_IP3D
norm_d3D
END_OF_TYPES

Definition at line 47 of file PVValidationHelpers.h.

                    {
    dxy = 1,
    dx = 2,
    dy = 3,
    dz = 4,
    IP2D = 5,
    resz = 6,
    IP3D = 7,
    d3D = 8,

    norm_dxy = 9,
    norm_dx = 10,
    norm_dy = 11,
    norm_dz = 12,
    norm_IP2D = 13,
    norm_resz = 14,
    norm_IP3D = 15,
    norm_d3D = 16,
    END_OF_TYPES = 17,
  };

Function Documentation

add()

void PVValHelper::add	(	std::map< std::string, TH1 *> &	h,
		TH1 *	hist
	)

Definition at line 12 of file PVValidationHelpers.cc.

References h, and compareTotals::hist.

Referenced by helpers.CloneTaskVisitor::__appendToTopTask(), Matriplex::__attribute__(), Options.Options::__getitem__(), SequenceTypes._ModuleSequenceType::_replaceIfHeldDirectly(), SequenceTypes.Schedule::_replaceIfHeldDirectly(), LHCInfoPopConSourceHandler::addEmptyPayload(), LHCInfoPerLSPopConSourceHandler::addEmptyPayload(), LHCInfoPerFillPopConSourceHandler::addEmptyPayload(), custom_jme_cff::AddJetID(), custom_jme_cff::AddPileUpJetIDVars(), custom_jme_cff::AddQGLTaggerVars(), helpers::addTaskToProcess(), SimTrackstersProducer::addTrackster(), custom_muon_cff::AddVariablesForMuon(), Config.ProcessModifier::apply(), SequenceTypes._ModuleSequenceType::associate(), objects.VertexAnalyzer.VertexAnalyzer::beginLoop(), core.PileUpAnalyzer.PileUpAnalyzer::beginLoop(), CosmicHitPairGenerator::CosmicHitPairGenerator(), CosmicHitTripletGenerator::CosmicHitTripletGenerator(), o2olib.O2OTool::create(), edmStreamStallGrapher::createPDFImage(), evf::FastMonEncoding::encode(), evf::FastMonEncoding::encodeString(), SequenceVisitors.NodeVisitor::enter(), helpers.CloneSequenceVisitor::enter(), SequenceTypes.NodeNameVisitor::enter(), BPHWriteSpecificDecay::fill(), Phase2ITValidateCluster::fillITHistos(), createIOVlist::fillJson(), Phase2OTValidateCluster::fillOTHistos(), evf::FastMonEncoding::fillReserved(), TrackerMap::find_layer(), symbols::get_libraries(), python.rootplot.core::get_plot_inputs(), symbols::get_symbols(), L1TPhase2CorrelatorOffline::getQuantile(), L1GctHfBitCountsLut::getThresholdsGct(), PixelForwardLayer::groupedCompatibleDetsV(), PixelForwardLayerPhase1::groupedCompatibleDetsV(), TIBRing::groupedCompatibleDetsV(), MillePedeVariables::increaseHitsX(), MillePedeVariables::increaseHitsY(), CTPPSRPAlignmentCorrectionsDataESSourceXMLCommon::Merge(), edmTracerLogToSimpleConfig.PathParser::parse(), edmTracerLogToSimpleConfig.ConsumesParser::parse(), python.rootplot.core::plot_hists_mpl(), L1GctLut< 16, 12 >::printLine(), core.JSONAnalyzer.JSONAnalyzer::process(), core.PileUpAnalyzer.PileUpAnalyzer::process(), objects.VertexAnalyzer.VertexAnalyzer::process(), TIBLayer::searchNeighbors(), TECLayer::searchNeighbors(), CompositeTECWedge::searchNeighbors(), PixelForwardLayer::searchNeighbors(), PixelForwardLayerPhase1::searchNeighbors(), Phase2EndcapSingleRing::searchNeighbors(), TIDRing::searchNeighbors(), Phase2EndcapRing::searchNeighbors(), TOBRod::searchNeighbors(), TBPLayer::searchNeighbors(), PixelBlade::searchNeighbors(), Phase2OTBarrelRod::searchNeighbors(), Phase1PixelBlade::searchNeighbors(), TIBRing::searchNeighbors(), CompositeTECPetal::searchNeighbors(), MatrixReader.MatrixReader::showRaw(), conddb_migrate::tags_in_gts(), theLHCInfoPerLSImpl::transferPayloads(), theLHCInfoPerFillImpl::transferPayloads(), theLHCInfoImpl::transferPayloads(), evf::FastMonEncoding::update(), evf::FastMonEncoding::updatePreinit(), evf::FastMonEncoding::updateReserved(), and querying.connection::write().

{
  h[hist->GetName()] = hist;
  hist->StatOverflows(kTRUE);
}

fill() [1/2]

void PVValHelper::fill	(	std::map< std::string, TH1 *> &	h,
		const std::string &	s,
		double	x
	)

Definition at line 20 of file PVValidationHelpers.cc.

References h, alignCSCRings::s, and x.

Referenced by PrimaryVertexValidation::fillTrackHistos().

{
  if (h.count(s) == 0) {
    edm::LogWarning("PVValidationHelpers") << "Trying to fill non-existing Histogram named " << s << std::endl;
    return;
  }
  h[s]->Fill(x);
}

fill() [2/2]

void PVValHelper::fill	(	std::map< std::string, TH1 *> &	h,
		const std::string &	s,
		double	x,
		double	y
	)

Definition at line 31 of file PVValidationHelpers.cc.

References h, alignCSCRings::s, x, and detailsBasic3DVector::y.

{
  if (h.count(s) == 0) {
    edm::LogWarning("PVValidationHelpers") << "Trying to fill non-existing Histogram named " << s << std::endl;
    return;
  }
  h[s]->Fill(x, y);
}

fillByIndex()

void PVValHelper::fillByIndex	(	std::vector< TH1F *> &	h,
		unsigned int	index,
		double	x,
		std::string	tag = ""
	)

Definition at line 42 of file PVValidationHelpers.cc.

References cms::cuda::assert(), h, makeGlobalPositionRcd_cfg::tag, and x.

Referenced by PrimaryVertexValidation::analyze(), SplitVertexResolution::analyze(), and DMRChecker::analyze().

{
  assert(!h.empty());
  if (index < h.size()) {
    h[index]->Fill(x);
  } else {
    edm::LogWarning("PVValidationHelpers") << "Trying to fill non-existing Histogram with index " << index
                                           << " for array with size: " << h.size() << " tag: " << tag << std::endl;
    return;
  }
}

fitResiduals()

std::pair< Measurement1D, Measurement1D > PVValHelper::fitResiduals

(

TH1 *

hist

)

Definition at line 225 of file PVValidationHelpers.cc.

References EcalMonitorTask_cff::func, compareTotals::hist, SiStripPI::mean, and mps_fire::result.

Referenced by PrimaryVertexValidation::fillMap(), PrimaryVertexValidation::fillTrendPlot(), and PrimaryVertexValidation::fillTrendPlotByIndex().

{
  //float fitResult(9999);
  if (hist->GetEntries() < 1) {
    return std::make_pair(Measurement1D(0., 0.), Measurement1D(0., 0.));
  };

  float mean = hist->GetMean();
  float sigma = hist->GetRMS();

  TF1 func("tmp", "gaus", mean - 1.5 * sigma, mean + 1.5 * sigma);
  if (0 == hist->Fit(&func, "QNR")) {  // N: do not blow up file by storing fit!
    mean = func.GetParameter(1);
    sigma = func.GetParameter(2);
    // second fit: three sigma of first fit around mean of first fit
    func.SetRange(mean - 2 * sigma, mean + 2 * sigma);
    // I: integral gives more correct results if binning is too wide
    // L: Likelihood can treat empty bins correctly (if hist not weighted...)
    if (0 == hist->Fit(&func, "Q0LR")) {
      if (hist->GetFunction(func.GetName())) {  // Take care that it is later on drawn:
        hist->GetFunction(func.GetName())->ResetBit(TF1::kNotDraw);
      }
    }
  }

  float res_mean = func.GetParameter(1);
  float res_width = func.GetParameter(2);

  float res_mean_err = func.GetParError(1);
  float res_width_err = func.GetParError(2);

  Measurement1D resultM(res_mean, res_mean_err);
  Measurement1D resultW(res_width, res_width_err);

  std::pair<Measurement1D, Measurement1D> result;

  result = std::make_pair(resultM, resultW);
  return result;
}

generateBins()

std::vector< float > PVValHelper::generateBins	(	int	n,
		float	start,
		float	range
	)

Definition at line 161 of file PVValidationHelpers.cc.

References a, mps_fire::end, writeEcalDQMStatus::interval, dqmiodumpmetadata::n, FastTimerService_cff::range, and findQualityFiles::v.

Referenced by PrimaryVertexValidation::PrimaryVertexValidation(), and SplitVertexResolution::SplitVertexResolution().

{
  std::vector<float> v(n);
  float interval = range / (n - 1);
  std::iota(v.begin(), v.end(), 1.);

  /*
    std::cout<<" interval:"<<interval<<std::endl;
    for(float &a : v) { std::cout<< a << " ";  }
    std::cout<< "\n";
  */

  std::for_each(begin(v), end(v), [&](float& a) { a = start + ((a - 1) * interval); });

  return v;
}

getMAD()

Measurement1D PVValHelper::getMAD

(

TH1F *

histo

)

Definition at line 196 of file PVValidationHelpers.cc.

References funct::abs(), getMedian(), timingPdfMaker::histo, dqmiolumiharvest::j, median(), LaserClient_cfi::nbins, mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, Measurement1D::value(), and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

Referenced by PrimaryVertexValidation::fillMap(), PrimaryVertexValidation::fillTrendPlot(), PrimaryVertexValidation::fillTrendPlotByIndex(), and SplitVertexResolution::fillTrendPlotByIndex().

{
  int nbins = histo->GetNbinsX();
  double median = getMedian(histo).value();
  double x_lastBin = histo->GetBinLowEdge(nbins + 1);
  const char* HistoName = histo->GetName();
  std::string Finalname = "resMed_";
  Finalname.append(HistoName);
  TH1F* newHisto = new TH1F(Finalname.c_str(), Finalname.c_str(), nbins, 0., x_lastBin);
  double* residuals = new double[nbins];
  const float* weights = histo->GetArray();

  for (int j = 0; j < nbins; j++) {
    residuals[j] = std::abs(median - histo->GetBinCenter(j + 1));
    newHisto->Fill(residuals[j], weights[j]);
  }

  double theMAD = (PVValHelper::getMedian(newHisto).value()) * 1.4826;

  delete[] residuals;
  residuals = nullptr;
  newHisto->Delete("");

  Measurement1D result(theMAD, theMAD / histo->GetEntries());
  return result;
}

getMedian()

Measurement1D PVValHelper::getMedian

(

TH1F *

histo

)

Definition at line 180 of file PVValidationHelpers.cc.

References timingPdfMaker::histo, median(), submitPVResolutionJobs::q, and mps_fire::result.

Referenced by PrimaryVertexValidation::fillMap(), PrimaryVertexValidation::fillTrendPlot(), PrimaryVertexValidation::fillTrendPlotByIndex(), SplitVertexResolution::fillTrendPlotByIndex(), and getMAD().

{
  double median = 0.;
  double q = 0.5;  // 0.5 quantile for "median"
  // protect against empty histograms
  if (histo->Integral() != 0) {
    histo->GetQuantiles(1, &median, &q);
  }

  Measurement1D result(median, median / TMath::Sqrt(histo->GetEntries()));

  return result;
}

getTypeString()

PVValHelper::plotLabels PVValHelper::getTypeString

(

PVValHelper::residualType

type

)

Definition at line 63 of file PVValidationHelpers.cc.

References d3D, dx, dxy, dy, dz, IP2D, IP3D, norm_d3D, norm_dx, norm_dxy, norm_dy, norm_dz, norm_IP2D, norm_IP3D, norm_resz, resz, and reco::returnType().

Referenced by PrimaryVertexValidation::bookResidualsHistogram(), and PrimaryVertexValidation::PrimaryVertexValidation().

{
  PVValHelper::plotLabels returnType;
  switch (type) {
      // absoulte

    case PVValHelper::dxy:
      returnType = std::make_tuple("dxy", "d_{xy}", "[#mum]");
      break;
    case PVValHelper::dx:
      returnType = std::make_tuple("dx", "d_{x}", "[#mum]");
      break;
    case PVValHelper::dy:
      returnType = std::make_tuple("dy", "d_{y}", "[#mum]");
      break;
    case PVValHelper::dz:
      returnType = std::make_tuple("dz", "d_{z}", "[#mum]");
      break;
    case PVValHelper::IP2D:
      returnType = std::make_tuple("IP2D", "IP_{2D}", "[#mum]");
      break;
    case PVValHelper::resz:
      returnType = std::make_tuple("resz", "z_{trk}-z_{vtx}", "[#mum]");
      break;
    case PVValHelper::IP3D:
      returnType = std::make_tuple("IP3D", "IP_{3D}", "[#mum]");
      break;
    case PVValHelper::d3D:
      returnType = std::make_tuple("d3D", "d_{3D}", "[#mum]");
      break;

      // normalized

    case PVValHelper::norm_dxy:
      returnType = std::make_tuple("norm_dxy", "d_{xy}/#sigma_{d_{xy}}", "");
      break;
    case PVValHelper::norm_dx:
      returnType = std::make_tuple("norm_dx", "d_{x}/#sigma_{d_{x}}", "");
      break;
    case PVValHelper::norm_dy:
      returnType = std::make_tuple("norm_dy", "d_{y}/#sigma_{d_{y}}", "");
      break;
    case PVValHelper::norm_dz:
      returnType = std::make_tuple("norm_dz", "d_{z}/#sigma_{d_{z}}", "");
      break;
    case PVValHelper::norm_IP2D:
      returnType = std::make_tuple("norm_IP2D", "IP_{2D}/#sigma_{IP_{2D}}", "");
      break;
    case PVValHelper::norm_resz:
      returnType = std::make_tuple("norm_resz", "z_{trk}-z_{vtx}/#sigma_{res_{z}}", "");
      break;
    case PVValHelper::norm_IP3D:
      returnType = std::make_tuple("norm_IP3D", "IP_{3D}/#sigma_{IP_{3D}}", "");
      break;
    case PVValHelper::norm_d3D:
      returnType = std::make_tuple("norm_d3D", "d_{3D}/#sigma_{d_{3D}}", "");
      break;

    default:
      edm::LogWarning("PVValidationHelpers") << " getTypeString() unknown residual type: " << type << std::endl;
  }

  return returnType;
}

getVarString()

PVValHelper::plotLabels PVValHelper::getVarString

(

PVValHelper::plotVariable

var

)

Definition at line 130 of file PVValidationHelpers.cc.

References eta, ladder, modZ, phi, pT, pTCentral, and trigObjTnPSource_cfi::var.

Referenced by PrimaryVertexValidation::bookResidualsHistogram(), and PrimaryVertexValidation::PrimaryVertexValidation().

{
  PVValHelper::plotLabels returnVar;
  switch (var) {
    case PVValHelper::phi:
      returnVar = std::make_tuple("phi", "#phi", "[rad]");
      break;
    case PVValHelper::eta:
      returnVar = std::make_tuple("eta", "#eta", "");
      break;
    case PVValHelper::pT:
      returnVar = std::make_tuple("pT", "p_{T}", "[GeV]");
      break;
    case PVValHelper::pTCentral:
      returnVar = std::make_tuple("pTCentral", "p_{T} |#eta|<1.", "[GeV]");
      break;
    case PVValHelper::ladder:
      returnVar = std::make_tuple("ladder", "ladder number", "");
      break;
    case PVValHelper::modZ:
      returnVar = std::make_tuple("modZ", "module number", "");
      break;
    default:
      edm::LogWarning("PVValidationHelpers") << " getVarString() unknown plot variable: " << var << std::endl;
  }

  return returnVar;
}

makeLogBins()

template<typename T , size_t N>

std::array<T, N + 1> PVValHelper::makeLogBins	(	const T &	min,
		const T &	max
	)

Definition at line 30 of file PVValidationHelpers.h.

References mps_fire::i, SiStripPI::max, SiStripPI::min, VarParsing::mult, N, funct::pow(), runTheMatrix::ret, and ApeEstimator_cff::width.

                                                             {
    const T minLog10 = std::log10(min);
    const T maxLog10 = std::log10(max);
    const T width = (maxLog10 - minLog10) / N;
    std::array<T, N + 1> ret;
    ret[0] = std::pow(10, minLog10);
    const T mult = std::pow(10, width);
    for (size_t i = 1; i <= N; ++i) {
      ret[i] = ret[i - 1] * mult;
    }
    return ret;
  }

shrinkHistVectorToFit()

void PVValHelper::shrinkHistVectorToFit	(	std::vector< TH1F *> &	h,
		unsigned int	desired_size
	)

Definition at line 56 of file PVValidationHelpers.cc.

References h.

Referenced by PrimaryVertexValidation::beginRun().

{
  h.erase(h.begin() + desired_size, h.end());
}

Variable Documentation

max_eta_phase0

constexpr double PVValHelper::max_eta_phase0 = 2.5

Definition at line 23 of file PVValidationHelpers.h.

Referenced by PrimaryVertexValidation::beginRun().

max_eta_phase1

constexpr double PVValHelper::max_eta_phase1 = 2.7

Definition at line 24 of file PVValidationHelpers.h.

Referenced by PrimaryVertexValidation::beginRun().

max_eta_phase2

constexpr double PVValHelper::max_eta_phase2 = 4.0

Definition at line 25 of file PVValidationHelpers.h.

Referenced by PrimaryVertexValidation::beginRun().