PVValHelper Namespace Reference

Classes
struct	histodetails

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

Enumerations
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)

Typedef Documentation

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

Definition at line 112 of file PVValidationHelpers.h.

Enumeration Type Documentation

enum PVValHelper::estimator

Enumerator
MEAN
WIDTH
MEDIAN
MAD
UNKWN

Definition at line 34 of file PVValidationHelpers.h.

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

enum PVValHelper::plotVariable

Enumerator
phi
eta
pT
pTCentral
ladder
modZ
END_OF_PLOTS

Definition at line 63 of file PVValidationHelpers.h.

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

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 42 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

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

Definition at line 12 of file PVValidationHelpers.cc.

References create_public_lumi_plots::hist.

Referenced by BPHRecoBuilder::add(), SiPixelQuality::addDisabledModule(), edm::eventsetup::EventSetupRecord::cacheIdentifier(), Folder::cd(), FastLineRecognition::Cluster::Cluster(), popcon::PopConSourceHandler< CastorPedestals >::convertFromOld(), CosmicHitPairGenerator::CosmicHitPairGenerator(), CosmicHitTripletGenerator::CosmicHitTripletGenerator(), edmStreamStallGrapher::createPDFImage(), BPHWriteSpecificDecay::fill(), evf::FastMonitoringService::Encoding::fillReserved(), CaloHitResponse::finalizeHits(), TrackerMap::find_layer(), symbols::get_libraries(), python.rootplot.core::get_plot_inputs(), symbols::get_symbols(), L1GctHfBitCountsLut::getThresholdsGct(), PixelForwardLayer::groupedCompatibleDetsV(), PixelForwardLayerPhase1::groupedCompatibleDetsV(), TIBRing::groupedCompatibleDetsV(), MillePedeVariables::increaseHitsX(), MillePedeVariables::increaseHitsY(), EBHitResponse::keepBlank(), CTPPSIncludeAlignments::Merge(), CSCCLCTData::nlines(), MtvClusterizer1D< T >::operator()(), FsmwClusterizer1D< T >::operator()(), muonisolation::NominalEfficiencyThresholds::locless::operator()(), amc13::Packet::Packet(), python.rootplot.core::plot_hists_mpl(), L1GctLut< NAddressBits, NDataBits >::printLine(), objects.VertexAnalyzer.VertexAnalyzer::process(), core.PileUpAnalyzer.PileUpAnalyzer::process(), fwlite::RunLumiSelector::RunLumiSelector(), TIBLayer::searchNeighbors(), CompositeTECWedge::searchNeighbors(), PixelForwardLayer::searchNeighbors(), TECLayer::searchNeighbors(), PixelForwardLayerPhase1::searchNeighbors(), TIDRing::searchNeighbors(), Phase2EndcapRing::searchNeighbors(), TOBRod::searchNeighbors(), TBPLayer::searchNeighbors(), PixelBlade::searchNeighbors(), Phase1PixelBlade::searchNeighbors(), Phase2OTBarrelRod::searchNeighbors(), TIBRing::searchNeighbors(), CompositeTECPetal::searchNeighbors(), AlignmentParameterSelector::selectedParameters(), edm::EventSetup::setKnownRecordsSupplier(), CSCBaseElectronicsSim::setSignalTimeRange(), MatrixReader.MatrixReader::showRaw(), CSCAnodeData2007::sizeInWords(), CSCAnodeData2006::sizeInWords(), CSCALCTs2006::sizeInWords(), conddb_migrate::tags_in_gts(), evf::FastMonitoringService::Encoding::update(), AlignmentParameterBuilder::~AlignmentParameterBuilder(), DTVDriftCalibration::cellInfo::~cellInfo(), PhysicsTools::Calibration::MVAComputerContainer::~MVAComputerContainer(), and pat::PATUserDataMerger< pat::pat::MET, pat::helper::AddUserInt >::~PATUserDataMerger().

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

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

Definition at line 20 of file PVValidationHelpers.cc.

References alignCSCRings::s.

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);
}

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 alignCSCRings::s.

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

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

Definition at line 42 of file PVValidationHelpers.cc.

References diffTreeTool::index.

Referenced by PrimaryVertexValidation::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;
  }
}

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

(

TH1 *

hist

)

Definition at line 230 of file PVValidationHelpers.cc.

References RecoEcal_EventContent_cff::func, SiStripPI::mean, and mps_fire::result.

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

{
  //float fitResult(9999);
  //if (hist->GetEntries() < 20) return ;
  
  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;
}

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

Definition at line 164 of file PVValidationHelpers.cc.

References a, begin, end, and findQualityFiles::v.

Referenced by PrimaryVertexValidation::PrimaryVertexValidation().

{

  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;

}

Measurement1D PVValHelper::getMAD

(

TH1F *

histo

)

Definition at line 200 of file PVValidationHelpers.cc.

References funct::abs(), getMedian(), pileupCalc::nbins, mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, Measurement1D::value(), and create_public_lumi_plots::weights.

Referenced by PrimaryVertexValidation::fillMap(), PrimaryVertexValidation::fillTrendPlot(), and PrimaryVertexValidation::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;
}

Measurement1D PVValHelper::getMedian

(

TH1F *

histo

)

Definition at line 183 of file PVValidationHelpers.cc.

References lumiQueryAPI::q, and mps_fire::result.

Referenced by PrimaryVertexValidation::fillMap(), PrimaryVertexValidation::fillTrendPlot(), PrimaryVertexValidation::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;

}

PVValHelper::plotLabels PVValHelper::getTypeString

(

PVValHelper::residualType

type

)

Definition at line 62 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;
  
}

PVValHelper::plotLabels PVValHelper::getVarString

(

PVValHelper::plotVariable

var

)

Definition at line 131 of file PVValidationHelpers.cc.

References eta, ladder, modZ, phi, pT, and pTCentral.

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;
  
}

template<typename T , size_t N>

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

Definition at line 19 of file PVValidationHelpers.h.

References mps_fire::i, N, funct::pow(), 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;
  }

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

Definition at line 55 of file PVValidationHelpers.cc.

Referenced by PrimaryVertexValidation::analyze().

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