PrimaryVertexValidation Class Reference

#include <Alignment/OfflineValidation/plugins/PrimaryVertexValidation.cc>

Inheritance diagram for PrimaryVertexValidation:

Public Member Functions
	PrimaryVertexValidation (const edm::ParameterSet &)
	~PrimaryVertexValidation ()
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	add (std::map< std::string, TH1 * > &h, TH1 *hist)
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
std::vector< TH1F * >	bookResidualsHistogram (TFileDirectory dir, unsigned int theNOfBins, TString resType, TString varType)
std::map< std::string, TH1 * >	bookVertexHistograms (TFileDirectory dir)
virtual void	endJob ()
void	fill (std::map< std::string, TH1 * > &h, std::string s, double x)
void	fill (std::map< std::string, TH1 * > &h, std::string s, double x, double y)
void	fillByIndex (std::vector< TH1F * > &h, unsigned int index, double x)
void	fillMap (TH2F *trendMap, TH1F *residualsMapPlot[100][100], statmode::estimator fitPar_)
void	fillTrackHistos (std::map< std::string, TH1 * > &h, const std::string &ttype, const reco::TransientTrack *tt, const reco::Vertex &v, const reco::BeamSpot &beamSpot, double fBfield)
void	fillTrendPlot (TH1F *trendPlot, TH1F *residualsPlot[100], statmode::estimator fitPar_, TString var_)
void	fillTrendPlotByIndex (TH1F *trendPlot, std::vector< TH1F * > &h, statmode::estimator fitPar_)
std::pair< std::pair< Double_t, Double_t >, std::pair< Double_t, Double_t > >	fitResiduals (TH1 *hist)
std::pair< Double_t, Double_t >	getMAD (TH1F *histo)
std::pair< Double_t, Double_t >	getMedian (TH1F *histo)
bool	hasFirstLayerPixelHits (const reco::TransientTrack track)
bool	isBFieldConsistentWithMode (const edm::EventSetup &iSetup) const
bool	isHit2D (const TrackingRecHit &hit) const
bool	passesTrackCuts (const reco::Track &track, const reco::Vertex &vertex, std::string qualityString_, double dxyErrMax_, double dzErrMax_, double ptErrMax_)
std::pair< bool, bool >	pixelHitsCheck (const reco::TransientTrack track)
void	SetVarToZero ()
double	square (double x)

Static Private Member Functions
static bool	vtxSort (const reco::Vertex &a, const reco::Vertex &b)

Private Attributes
TH1F *	a_d3DEtaResiduals [nMaxBins_]
TH1F *	a_d3DPhiResiduals [nMaxBins_]
TH1F *	a_d3DResidualsMap [nMaxBins_][nMaxBins_]
TH1F *	a_dxEtaResiduals [nMaxBins_]
TH1F *	a_dxPhiResiduals [nMaxBins_]
TH1F *	a_dxyBiasResidualsMap [nMaxBins_][nMaxBins_]
TH1F *	a_dxyEtaBiasResiduals [nMaxBins_]
TH1F *	a_dxyEtaMADBiasTrend
TH1F *	a_dxyEtaMADTrend
TH1F *	a_dxyEtaMeanBiasTrend
TH1F *	a_dxyEtaMeanTrend
TH1F *	a_dxyEtaMedianBiasTrend
TH1F *	a_dxyEtaMedianTrend
TH1F *	a_dxyEtaResiduals [nMaxBins_]
TH1F *	a_dxyEtaWidthBiasTrend
TH1F *	a_dxyEtaWidthTrend
TH2F *	a_dxyMeanBiasMap
TH2F *	a_dxyMeanMap
TH1F *	a_dxyPhiBiasResiduals [nMaxBins_]
TH1F *	a_dxyPhiMADBiasTrend
TH1F *	a_dxyPhiMADTrend
TH1F *	a_dxyPhiMeanBiasTrend
TH1F *	a_dxyPhiMeanTrend
TH1F *	a_dxyPhiMedianBiasTrend
TH1F *	a_dxyPhiMedianTrend
TH1F *	a_dxyPhiResiduals [nMaxBins_]
TH1F *	a_dxyPhiWidthBiasTrend
TH1F *	a_dxyPhiWidthTrend
TH1F *	a_dxypTCentralMeanTrend
TH1F *	a_dxypTCentralWidthTrend
TH1F *	a_dxypTMeanTrend
TH1F *	a_dxypTWidthTrend
TH1F *	a_dxyResidualsMap [nMaxBins_][nMaxBins_]
TH2F *	a_dxyVsEta
TH2F *	a_dxyVsPhi
TH2F *	a_dxyWidthBiasMap
TH2F *	a_dxyWidthMap
TH1F *	a_dyEtaResiduals [nMaxBins_]
TH1F *	a_dyPhiResiduals [nMaxBins_]
TH1F *	a_dzBiasResidualsMap [nMaxBins_][nMaxBins_]
TH1F *	a_dzEtaBiasResiduals [nMaxBins_]
TH1F *	a_dzEtaMADBiasTrend
TH1F *	a_dzEtaMADTrend
TH1F *	a_dzEtaMeanBiasTrend
TH1F *	a_dzEtaMeanTrend
TH1F *	a_dzEtaMedianBiasTrend
TH1F *	a_dzEtaMedianTrend
TH1F *	a_dzEtaResiduals [nMaxBins_]
TH1F *	a_dzEtaWidthBiasTrend
TH1F *	a_dzEtaWidthTrend
TH2F *	a_dzMeanBiasMap
TH2F *	a_dzMeanMap
TH1F *	a_dzPhiBiasResiduals [nMaxBins_]
TH1F *	a_dzPhiMADBiasTrend
TH1F *	a_dzPhiMADTrend
TH1F *	a_dzPhiMeanBiasTrend
TH1F *	a_dzPhiMeanTrend
TH1F *	a_dzPhiMedianBiasTrend
TH1F *	a_dzPhiMedianTrend
TH1F *	a_dzPhiResiduals [nMaxBins_]
TH1F *	a_dzPhiWidthBiasTrend
TH1F *	a_dzPhiWidthTrend
TH1F *	a_dzpTCentralMeanTrend
TH1F *	a_dzpTCentralWidthTrend
TH1F *	a_dzpTMeanTrend
TH1F *	a_dzpTWidthTrend
TH1F *	a_dzResidualsMap [nMaxBins_][nMaxBins_]
TH2F *	a_dzVsEta
TH2F *	a_dzVsPhi
TH2F *	a_dzWidthBiasMap
TH2F *	a_dzWidthMap
TH1F *	a_IP2DEtaResiduals [nMaxBins_]
TH1F *	a_IP2DPhiResiduals [nMaxBins_]
TH1F *	a_IP3DEtaResiduals [nMaxBins_]
TH1F *	a_IP3DPhiResiduals [nMaxBins_]
TH1F *	a_reszEtaResiduals [nMaxBins_]
TH1F *	a_reszPhiResiduals [nMaxBins_]
bool	askFirstLayerHit_
double	Beamdxdz_
double	Beamsigmaz_
double	BeamWidthX_
double	BeamWidthY_
double	BSx0_
double	BSy0_
double	BSz0_
int	charge_ [nMaxtracks_]
double	chi2_ [nMaxtracks_]
double	chi2ndof_ [nMaxtracks_]
float	chi2normUnbiasedVertex_ [nMaxtracks_]
float	chi2ProbUnbiasedVertex_ [nMaxtracks_]
float	chi2UnbiasedVertex_ [nMaxtracks_]
double	d3DErrorFromMyVertex_ [nMaxtracks_]
double	d3DFromMyVertex_ [nMaxtracks_]
bool	debug_
bool	doBPix_
bool	doFPix_
float	DOFUnbiasedVertex_ [nMaxtracks_]
double	dxy_ [nMaxtracks_]
double	dxyBs_ [nMaxtracks_]
double	dxyErrorFromMyVertex_ [nMaxtracks_]
double	dxyFromMyVertex_ [nMaxtracks_]
double	dz_ [nMaxtracks_]
double	dzBs_ [nMaxtracks_]
double	dzErrorFromMyVertex_ [nMaxtracks_]
double	dzFromMyVertex_ [nMaxtracks_]
double	eta_ [nMaxtracks_]
double	etaOfProbe_
float	etaSect_
unsigned int	EventNumber_
TH2F *	h2_probeEtaPhi_
TH2F *	h2_probeEtaPt_
TH1F *	h_Beamsigmaz
TH1F *	h_BeamWidthX
TH1F *	h_BeamWidthY
TH1F *	h_BSx0
TH1F *	h_BSy0
TH1F *	h_BSz0
std::vector< TH1F * >	h_dxy_Central_pT_
std::vector< TH1F * >	h_dxy_pT_
std::vector< TH1F * >	h_dz_Central_pT_
std::vector< TH1F * >	h_dz_pT_
TH1F *	h_etaMax
TH1F *	h_fitVtxChi2_
TH1F *	h_fitVtxChi2ndf_
TH1F *	h_fitVtxChi2Prob_
TH1F *	h_fitVtxNdof_
TH1F *	h_fitVtxNtracks_
TH1F *	h_fitVtxTrackAverageWeight_
TH1F *	h_fitVtxTrackWeights_
TH1F *	h_lumiFromConfig
TH1F *	h_nbins
TH1F *	h_nClus
TH1F *	h_nOfflineVertices
std::vector< TH1F * >	h_norm_dxy_Central_pT_
std::vector< TH1F * >	h_norm_dxy_pT_
std::vector< TH1F * >	h_norm_dz_Central_pT_
std::vector< TH1F * >	h_norm_dz_pT_
TH1F *	h_nTracks
TH1F *	h_probeCharge_
TH1F *	h_probeChi2_
TH1F *	h_probed0RefitV_
TH1F *	h_probed3DRefitV_
TH1F *	h_probedxyRecoV_
TH1F *	h_probedxyRefitV_
TH1F *	h_probedzRecoV_
TH1F *	h_probedzRefitV_
TH1F *	h_probeEta_
TH1F *	h_probeHits1D_
TH1F *	h_probeHits2D_
TH1F *	h_probeHits_
TH1F *	h_probeHitsInBPIX_
TH1F *	h_probeHitsInFPIX_
TH1F *	h_probeHitsInTEC_
TH1F *	h_probeHitsInTIB_
TH1F *	h_probeHitsInTID_
TH1F *	h_probeHitsInTOB_
TH1F *	h_probeNormChi2_
TH1F *	h_probeP_
TH1F *	h_probePhi_
TH1F *	h_probePt_
TH1F *	h_probeQoverP_
TH1F *	h_probeRecoVSigXY_
TH1F *	h_probeRecoVSigZ_
TH1F *	h_probeRefitVLogSig3D_
TH1F *	h_probeRefitVSig3D_
TH1F *	h_probeRefitVSigResZ_
TH1F *	h_probeRefitVSigXY_
TH1F *	h_probeRefitVSigZ_
TH1F *	h_probereszRefitV_
TH1F *	h_probesignIP2DRefitV_
TH1F *	h_probez0RefitV_
TH1F *	h_recoVtxChi2ndf_
TH1F *	h_recoVtxChi2Prob_
TH1F *	h_recoVtxNtracks_
TH1F *	h_recoVtxSumPt_
TH1I *	h_runFromConfig
TH1I *	h_runFromEvent
TH1F *	h_runNumber
TH1F *	h_xErrOfflineVertex
TH1F *	h_xOfflineVertex
TH1F *	h_yErrOfflineVertex
TH1F *	h_yOfflineVertex
TH1F *	h_zErrOfflineVertex
TH1F *	h_zOfflineVertex
int	hasRecVertex_ [nMaxtracks_]
std::map< std::string, TH1 * >	hDA
double	intLumi_
double	IP3DsigFromMyVertex_ [nMaxtracks_]
double	IPLsigFromMyVertex_ [nMaxtracks_]
double	IPTsigFromMyVertex_ [nMaxtracks_]
int	isGoodTrack_ [nMaxtracks_]
int	isHighPurity_ [nMaxtracks_]
bool	isPhase1_
bool	lightNtupleSwitch_
unsigned int	LuminosityBlockNumber_
const float	mypT_bins_ [nPtBins_+1] = {0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3.0,3.1,3.2,3.3,3.4,3.5,3.6,3.7,3.8,3.9,4.0,4.1,4.25,4.5,4.75,5.0,5.5,6.0,7.0,8.0,9.0,11.0,14.0,20.}
TH1F *	n_d3DEtaResiduals [nMaxBins_]
TH1F *	n_d3DPhiResiduals [nMaxBins_]
TH1F *	n_d3DResidualsMap [nMaxBins_][nMaxBins_]
TH1F *	n_dxyBiasResidualsMap [nMaxBins_][nMaxBins_]
TH1F *	n_dxyEtaBiasResiduals [nMaxBins_]
TH1F *	n_dxyEtaMADBiasTrend
TH1F *	n_dxyEtaMADTrend
TH1F *	n_dxyEtaMeanBiasTrend
TH1F *	n_dxyEtaMeanTrend
TH1F *	n_dxyEtaMedianBiasTrend
TH1F *	n_dxyEtaMedianTrend
TH1F *	n_dxyEtaResiduals [nMaxBins_]
TH1F *	n_dxyEtaWidthBiasTrend
TH1F *	n_dxyEtaWidthTrend
TH2F *	n_dxyMeanBiasMap
TH2F *	n_dxyMeanMap
TH1F *	n_dxyPhiBiasResiduals [nMaxBins_]
TH1F *	n_dxyPhiMADBiasTrend
TH1F *	n_dxyPhiMADTrend
TH1F *	n_dxyPhiMeanBiasTrend
TH1F *	n_dxyPhiMeanTrend
TH1F *	n_dxyPhiMedianBiasTrend
TH1F *	n_dxyPhiMedianTrend
TH1F *	n_dxyPhiResiduals [nMaxBins_]
TH1F *	n_dxyPhiWidthBiasTrend
TH1F *	n_dxyPhiWidthTrend
TH1F *	n_dxypTCentralMeanTrend
TH1F *	n_dxypTCentralWidthTrend
TH1F *	n_dxypTMeanTrend
TH1F *	n_dxypTWidthTrend
TH1F *	n_dxyResidualsMap [nMaxBins_][nMaxBins_]
TH2F *	n_dxyVsEta
TH2F *	n_dxyVsPhi
TH2F *	n_dxyWidthBiasMap
TH2F *	n_dxyWidthMap
TH1F *	n_dzBiasResidualsMap [nMaxBins_][nMaxBins_]
TH1F *	n_dzEtaBiasResiduals [nMaxBins_]
TH1F *	n_dzEtaMADBiasTrend
TH1F *	n_dzEtaMADTrend
TH1F *	n_dzEtaMeanBiasTrend
TH1F *	n_dzEtaMeanTrend
TH1F *	n_dzEtaMedianBiasTrend
TH1F *	n_dzEtaMedianTrend
TH1F *	n_dzEtaResiduals [nMaxBins_]
TH1F *	n_dzEtaWidthBiasTrend
TH1F *	n_dzEtaWidthTrend
TH2F *	n_dzMeanBiasMap
TH2F *	n_dzMeanMap
TH1F *	n_dzPhiBiasResiduals [nMaxBins_]
TH1F *	n_dzPhiMADBiasTrend
TH1F *	n_dzPhiMADTrend
TH1F *	n_dzPhiMeanBiasTrend
TH1F *	n_dzPhiMeanTrend
TH1F *	n_dzPhiMedianBiasTrend
TH1F *	n_dzPhiMedianTrend
TH1F *	n_dzPhiResiduals [nMaxBins_]
TH1F *	n_dzPhiWidthBiasTrend
TH1F *	n_dzPhiWidthTrend
TH1F *	n_dzpTCentralMeanTrend
TH1F *	n_dzpTCentralWidthTrend
TH1F *	n_dzpTMeanTrend
TH1F *	n_dzpTWidthTrend
TH1F *	n_dzResidualsMap [nMaxBins_][nMaxBins_]
TH2F *	n_dzVsEta
TH2F *	n_dzVsPhi
TH2F *	n_dzWidthBiasMap
TH2F *	n_dzWidthMap
TH1F *	n_IP2DEtaResiduals [nMaxBins_]
TH1F *	n_IP2DPhiResiduals [nMaxBins_]
TH1F *	n_IP3DEtaResiduals [nMaxBins_]
TH1F *	n_IP3DPhiResiduals [nMaxBins_]
TH1F *	n_reszEtaResiduals [nMaxBins_]
TH1F *	n_reszPhiResiduals [nMaxBins_]
int	nBins_
int	nClus_
int	Nevt_
int	nhits1D_ [nMaxtracks_]
int	nhits2D_ [nMaxtracks_]
int	nhits_ [nMaxtracks_]
int	nhitsBPIX_ [nMaxtracks_]
int	nhitsFPIX_ [nMaxtracks_]
int	nhitsTEC_ [nMaxtracks_]
int	nhitsTIB_ [nMaxtracks_]
int	nhitsTID_ [nMaxtracks_]
int	nhitsTOB_ [nMaxtracks_]
int	nOfflineVertices_
int	nTracks_
int	nTracksPerClus_
double	p_ [nMaxtracks_]
double	phi_ [nMaxtracks_]
float	phiSect_
double	pt_ [nMaxtracks_]
double	ptOfProbe_
double	qoverp_ [nMaxtracks_]
TTree *	rootTree_
bool	runControl_
std::vector< unsigned int >	runControlNumbers_
unsigned int	RunNumber_
bool	storeNtuple_
float	sumOfWeightsUnbiasedVertex_ [nMaxtracks_]
edm::EDGetTokenT< reco::BeamSpot >	theBeamspotToken
edm::ParameterSet	theConfig
double	theta_ [nMaxtracks_]
TrackClusterizerInZ *	theTrackClusterizer_
edm::EDGetTokenT< reco::TrackCollection >	theTrackCollectionToken
TrackFilterForPVFindingBase *	theTrackFilter_
edm::EDGetTokenT< reco::VertexCollection >	theVertexCollectionToken
int	tracksUsedForVertexing_ [nMaxtracks_]
bool	useTracksFromRecoVtx_
double	vertexZMax_
double	wxy2_
double	xErrOfflineVertex_
double	xOfflineVertex_
double	xPCA_ [nMaxtracks_]
double	xUnbiasedVertex_ [nMaxtracks_]
double	yErrOfflineVertex_
double	yOfflineVertex_
double	yPCA_ [nMaxtracks_]
double	yUnbiasedVertex_ [nMaxtracks_]
double	zErrOfflineVertex_
double	zOfflineVertex_
double	zPCA_ [nMaxtracks_]
double	zUnbiasedVertex_ [nMaxtracks_]

Static Private Attributes
static const int	cmToum = 10000
static const int	nMaxBins_ = 100
static const int	nMaxtracks_ = 1000
static const int	nPtBins_ = 48

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: Validate alignment constants using unbiased vertex residuals

Implementation: <Notes on="" implementation>="">

Definition at line 75 of file PrimaryVertexValidation.h.

Constructor & Destructor Documentation

PrimaryVertexValidation::PrimaryVertexValidation ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 67 of file PrimaryVertexValidation.cc.

References etaOfProbe_, etaSect_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), TFileService::kSharedResource, nBins_, phiSect_, Pi, runControlNumbers_, AlCaHLTBitMon_QueryRunRegistry::string, theBeamspotToken, theTrackClusterizer_, theTrackCollectionToken, theTrackFilter_, and theVertexCollectionToken.

                                                                              :
  storeNtuple_(iConfig.getParameter<bool>("storeNtuple")),
  lightNtupleSwitch_(iConfig.getParameter<bool>("isLightNtuple")),
  useTracksFromRecoVtx_(iConfig.getParameter<bool>("useTracksFromRecoVtx")),
  vertexZMax_(iConfig.getUntrackedParameter<double>("vertexZMax",99.)),
  intLumi_(iConfig.getUntrackedParameter<double>("intLumi",0.)),
  askFirstLayerHit_(iConfig.getParameter<bool>("askFirstLayerHit")),
  doBPix_(iConfig.getUntrackedParameter<bool>("doBPix",true)),
  doFPix_(iConfig.getUntrackedParameter<bool>("doFPix",true)),
  ptOfProbe_(iConfig.getUntrackedParameter<double>("probePt",0.)),
  etaOfProbe_(iConfig.getUntrackedParameter<double>("probeEta",2.4)),
  nBins_(iConfig.getUntrackedParameter<int>("numberOfBins",24)),
  debug_(iConfig.getParameter<bool>("Debug")),
  runControl_(iConfig.getUntrackedParameter<bool>("runControl",false))
{
  
  // now do what ever initialization is needed
  // initialize phase space boundaries
  
  usesResource(TFileService::kSharedResource);

  std::vector<unsigned int> defaultRuns;
  defaultRuns.push_back(0);
  runControlNumbers_ = iConfig.getUntrackedParameter<std::vector<unsigned int> >("runControlNumber",defaultRuns);

  phiSect_ = (2*TMath::Pi())/nBins_;
  etaSect_ = (2*etaOfProbe_)/nBins_;

  edm::InputTag TrackCollectionTag_ = iConfig.getParameter<edm::InputTag>("TrackCollectionTag");
  theTrackCollectionToken = consumes<reco::TrackCollection>(TrackCollectionTag_);

  edm::InputTag VertexCollectionTag_ = iConfig.getParameter<edm::InputTag>("VertexCollectionTag");
  theVertexCollectionToken = consumes<reco::VertexCollection>(VertexCollectionTag_);

  edm::InputTag BeamspotTag_ = edm::InputTag("offlineBeamSpot");
  theBeamspotToken = consumes<reco::BeamSpot>(BeamspotTag_);

  // select and configure the track filter 
  theTrackFilter_= new TrackFilterForPVFinding(iConfig.getParameter<edm::ParameterSet>("TkFilterParameters") );
  // select and configure the track clusterizer  
  std::string clusteringAlgorithm=iConfig.getParameter<edm::ParameterSet>("TkClusParameters").getParameter<std::string>("algorithm");
  if (clusteringAlgorithm=="gap"){
    theTrackClusterizer_ = new GapClusterizerInZ(iConfig.getParameter<edm::ParameterSet>("TkClusParameters").getParameter<edm::ParameterSet>("TkGapClusParameters"));
  }else if(clusteringAlgorithm=="DA"){
    theTrackClusterizer_ = new DAClusterizerInZ(iConfig.getParameter<edm::ParameterSet>("TkClusParameters").getParameter<edm::ParameterSet>("TkDAClusParameters"));
    // provide the vectorized version of the clusterizer, if supported by the build
  } else if(clusteringAlgorithm=="DA_vect") {
    theTrackClusterizer_ = new DAClusterizerInZ_vect(iConfig.getParameter<edm::ParameterSet>("TkClusParameters").getParameter<edm::ParameterSet>("TkDAClusParameters"));
  }else{
    throw VertexException("PrimaryVertexProducerAlgorithm: unknown clustering algorithm: " + clusteringAlgorithm);  
  }
}

PrimaryVertexValidation::~PrimaryVertexValidation

(

)

Definition at line 121 of file PrimaryVertexValidation.cc.

{
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

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

private

Definition at line 2748 of file PrimaryVertexValidation.cc.

References create_public_lumi_plots::hist.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

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

void PrimaryVertexValidation::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

Definition at line 134 of file PrimaryVertexValidation.cc.

References a_d3DEtaResiduals, a_d3DPhiResiduals, a_d3DResidualsMap, a_dxEtaResiduals, a_dxPhiResiduals, a_dxyBiasResidualsMap, a_dxyEtaBiasResiduals, a_dxyEtaResiduals, a_dxyPhiBiasResiduals, a_dxyPhiResiduals, a_dxyResidualsMap, a_dxyVsEta, a_dxyVsPhi, a_dyEtaResiduals, a_dyPhiResiduals, a_dzBiasResidualsMap, a_dzEtaBiasResiduals, a_dzEtaResiduals, a_dzPhiBiasResiduals, a_dzPhiResiduals, a_dzResidualsMap, a_dzVsEta, a_dzVsPhi, a_IP2DEtaResiduals, a_IP2DPhiResiduals, a_IP3DEtaResiduals, a_IP3DPhiResiduals, a_reszEtaResiduals, a_reszPhiResiduals, funct::abs(), IPTools::absoluteImpactParameter3D(), askFirstLayerHit_, Beamdxdz_, Beamsigmaz_, ecalDrivenElectronSeedsParameters_cff::beamSpot, reco::BeamSpot::BeamWidthX(), BeamWidthX_, reco::BeamSpot::BeamWidthY(), BeamWidthY_, BSx0_, BSy0_, BSz0_, TransientTrackBuilder::build(), reco::TrackBase::charge(), charge_, reco::Vertex::chi2(), reco::TrackBase::chi2(), chi2_, chi2ndof_, chi2normUnbiasedVertex_, chi2ProbUnbiasedVertex_, chi2UnbiasedVertex_, TrackClusterizerInZ::clusterize(), fastPrimaryVertexProducer_cfi::clusters, cmToum, GlobalErrorBase< T, ErrorWeightType >::czz(), allConversions_cfi::d0, reco::TrackBase::d0Error(), d3DErrorFromMyVertex_, d3DFromMyVertex_, debug_, TransientVertex::degreesOfFreedom(), TrackerGeometry::dets(), doBPix_, doFPix_, DOFUnbiasedVertex_, reco::BeamSpot::dxdz(), reco::TrackBase::dxy(), dxy_, dxyBs_, dxyErrorFromMyVertex_, dxyFromMyVertex_, reco::TrackBase::dz(), dz_, dzBs_, reco::TrackBase::dzError(), dzErrorFromMyVertex_, dzFromMyVertex_, reco::TrackBase::eta(), eta_, etaOfProbe_, etaSect_, edm::EventID::event(), edm::Event::eventAuxiliary(), EventNumber_, fillByIndex(), fillTrackHistos(), edm::EventSetup::get(), edm::Event::getByToken(), h2_probeEtaPhi_, h2_probeEtaPt_, h_Beamsigmaz, h_BeamWidthX, h_BeamWidthY, h_BSx0, h_BSy0, h_BSz0, h_dxy_Central_pT_, h_dxy_pT_, h_dz_Central_pT_, h_dz_pT_, h_etaMax, h_fitVtxChi2_, h_fitVtxChi2ndf_, h_fitVtxChi2Prob_, h_fitVtxNdof_, h_fitVtxNtracks_, h_fitVtxTrackAverageWeight_, h_fitVtxTrackWeights_, h_nbins, h_nClus, h_nOfflineVertices, h_norm_dxy_Central_pT_, h_norm_dxy_pT_, h_norm_dz_Central_pT_, h_norm_dz_pT_, h_nTracks, h_probeCharge_, h_probeChi2_, h_probed0RefitV_, h_probed3DRefitV_, h_probedxyRecoV_, h_probedxyRefitV_, h_probedzRecoV_, h_probedzRefitV_, h_probeEta_, h_probeHits1D_, h_probeHits2D_, h_probeHits_, h_probeHitsInBPIX_, h_probeHitsInFPIX_, h_probeHitsInTEC_, h_probeHitsInTIB_, h_probeHitsInTID_, h_probeHitsInTOB_, h_probeNormChi2_, h_probeP_, h_probePhi_, h_probePt_, h_probeQoverP_, h_probeRecoVSigXY_, h_probeRecoVSigZ_, h_probeRefitVLogSig3D_, h_probeRefitVSig3D_, h_probeRefitVSigResZ_, h_probeRefitVSigXY_, h_probeRefitVSigZ_, h_probereszRefitV_, h_probesignIP2DRefitV_, h_probez0RefitV_, h_recoVtxChi2ndf_, h_recoVtxChi2Prob_, h_recoVtxNtracks_, h_recoVtxSumPt_, h_runFromEvent, h_runNumber, h_xErrOfflineVertex, h_xOfflineVertex, h_yErrOfflineVertex, h_yOfflineVertex, h_zErrOfflineVertex, h_zOfflineVertex, hasFirstLayerPixelHits(), hasRecVertex_, TransientVertex::hasTrackWeight(), hDA, reco::TrackBase::hitPattern(), hfClusterShapes_cfi::hits, mps_fire::i, edm::EventAuxiliary::id(), createfilelist::int, IP3DsigFromMyVertex_, IPLsigFromMyVertex_, IPTsigFromMyVertex_, isBFieldConsistentWithMode(), isGoodTrack_, isHighPurity_, isHit2D(), isPhase1_, TrackerGeometry::isThere(), reco::Vertex::isValid(), edm::HandleBase::isValid(), TransientVertex::isValid(), gen::k, LogTrace, PerigeeTrajectoryParameters::longitudinalImpactParameter(), PerigeeTrajectoryError::longitudinalImpactParameterError(), edm::EventAuxiliary::luminosityBlock(), LuminosityBlockNumber_, min(), mypT_bins_, n_d3DEtaResiduals, n_d3DPhiResiduals, n_d3DResidualsMap, n_dxyBiasResidualsMap, n_dxyEtaBiasResiduals, n_dxyEtaResiduals, n_dxyPhiBiasResiduals, n_dxyPhiResiduals, n_dxyResidualsMap, n_dxyVsEta, n_dxyVsPhi, n_dzBiasResidualsMap, n_dzEtaBiasResiduals, n_dzEtaResiduals, n_dzPhiBiasResiduals, n_dzPhiResiduals, n_dzResidualsMap, n_dzVsEta, n_dzVsPhi, n_IP2DEtaResiduals, n_IP2DPhiResiduals, n_IP3DEtaResiduals, n_IP3DPhiResiduals, n_reszEtaResiduals, n_reszPhiResiduals, nBins_, nClus_, reco::Vertex::ndof(), Nevt_, nhits1D_, nhits2D_, nhits_, nhitsBPIX_, nhitsFPIX_, nhitsTEC_, nhitsTIB_, nhitsTID_, nhitsTOB_, nMaxtracks_, nOfflineVertices_, TransientVertex::normalisedChiSquared(), reco::Vertex::normalizedChi2(), reco::TrackBase::normalizedChi2(), nPtBins_, nTracks_, nTracksPerClus_, reco::TrackBase::numberOfValidHits(), reco::HitPattern::numberOfValidPixelBarrelHits(), reco::HitPattern::numberOfValidPixelEndcapHits(), reco::HitPattern::numberOfValidStripTECHits(), reco::HitPattern::numberOfValidStripTIBHits(), reco::HitPattern::numberOfValidStripTIDHits(), reco::HitPattern::numberOfValidStripTOBHits(), reco::TrackBase::p(), GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, p_, TrajectoryStateClosestToPoint::perigeeError(), TrajectoryStateClosestToPoint::perigeeParameters(), reco::TrackBase::phi(), phi_, phiSect_, Pi, pixelHitsCheck(), point, FreeTrajectoryState::position(), TrajectoryStateClosestToPoint::position(), reco::Vertex::position(), TransientVertex::position(), position, TransientVertex::positionError(), funct::pow(), EnergyCorrector::pt, reco::TrackBase::pt(), pt_, ptOfProbe_, MetAnalyzer::pv(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), reco::TrackBase::qoverp(), qoverp_, reco::TrackBase::quality(), rootTree_, edm::EventAuxiliary::run(), runControl_, runControlNumbers_, RunNumber_, TrackFilterForPVFindingBase::select(), reco::TransientTrack::setBeamSpot(), SetVarToZero(), reco::BeamSpot::sigmaZ(), IPTools::signedImpactParameter3D(), IPTools::signedTransverseImpactParameter(), mathSSE::sqrt(), storeNtuple_, sumOfWeightsUnbiasedVertex_, funct::tan(), theBeamspotToken, TrajectoryStateClosestToPoint::theState(), reco::TrackBase::theta(), theta_, theTrackClusterizer_, theTrackCollectionToken, theTrackFilter_, theVertexCollectionToken, TransientVertex::totalChiSquared(), HiIsolationCommonParameters_cff::track, reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), reco::Vertex::tracksSize(), tracksUsedForVertexing_, TransientVertex::trackWeight(), PerigeeTrajectoryParameters::transverseImpactParameter(), groupFilesInBlocks::tt, useTracksFromRecoVtx_, VertexFitter< N >::vertex(), reco::TrackBase::vertex(), vertexZMax_, ecalDrivenElectronSeedsParameters_cff::vertices, vtxSort(), cms::Exception::what(), wxy2_, x, PV3DBase< T, PVType, FrameType >::x(), reco::BeamSpot::x0(), xErrOfflineVertex_, xOfflineVertex_, xPCA_, xUnbiasedVertex_, y, PV3DBase< T, PVType, FrameType >::y(), reco::BeamSpot::y0(), yErrOfflineVertex_, yOfflineVertex_, yPCA_, yUnbiasedVertex_, z, PV3DBase< T, PVType, FrameType >::z(), reco::BeamSpot::z0(), zErrOfflineVertex_, zOfflineVertex_, zPCA_, zUnbiasedVertex_, and CombinatorialSeedGeneratorForCosmicsRegionalReconstruction_cfi::zVertex.

{

  using namespace std;
  using namespace reco;
  using namespace IPTools;

  if (!isBFieldConsistentWithMode(iSetup)) {
    edm::LogWarning("PrimaryVertexValidation") << "*********************************************************************************\n" 
                           << "* The configuration (ptOfProbe > " << ptOfProbe_ << "GeV) is not correctly set for current value of magnetic field \n" 
                           << "* Switching it to 0. !!! \n"
                           << "*********************************************************************************"<< std::endl;
    ptOfProbe_=0.;
  }

  if(nBins_!=24 && debug_){ 
    edm::LogInfo("PrimaryVertexValidation")<<"Using: "<<nBins_<<" bins plots";
  }
  
  bool passesRunControl = false;

  if(runControl_){
    for(unsigned int j=0;j<runControlNumbers_.size();j++){
      if(iEvent.eventAuxiliary().run() == runControlNumbers_[j]){ 
    if (debug_){
      edm::LogInfo("PrimaryVertexValidation")<<" run number: "<<iEvent.eventAuxiliary().run()<<" keeping run:"<<runControlNumbers_[j];
    }
    passesRunControl = true;
    break;
      }
    }
    if (!passesRunControl) return;
  }

  Nevt_++;  

  //=======================================================
  // Initialize Root-tuple variables
  //=======================================================

  SetVarToZero();
 
  //=======================================================
  // Retrieve the Magnetic Field information
  //=======================================================

  edm::ESHandle<MagneticField> theMGField;
  iSetup.get<IdealMagneticFieldRecord>().get( theMGField );

  //=======================================================
  // Retrieve the Tracking Geometry information
  //=======================================================

  edm::ESHandle<GlobalTrackingGeometry> theTrackingGeometry;
  iSetup.get<GlobalTrackingGeometryRecord>().get( theTrackingGeometry );

  //=======================================================
  // Retrieve geometry information
  //=======================================================
 
  edm::LogInfo("read tracker geometry...");
  edm::ESHandle<TrackerGeometry> pDD;
  iSetup.get<TrackerDigiGeometryRecord>().get( pDD );
  edm::LogInfo("tracker geometry read")<<"There are: "<< pDD->dets().size() <<" detectors";

  // switch on the phase1 
  if( (pDD->isThere(GeomDetEnumerators::P1PXB)) || 
      (pDD->isThere(GeomDetEnumerators::P1PXEC)) ) {
    isPhase1_ = true;
    if (debug_){
      edm::LogInfo("PrimaryVertexValidation")<<" pixel phase1 setup ";
    }
  } else {
    isPhase1_ = false;
    if (debug_){
      edm::LogInfo("PrimaryVertexValidation")<<" pixel phase0 setup ";
    }
  }

  if(isPhase1_){
    etaOfProbe_ = std::min(etaOfProbe_,2.7);
  } else {
    etaOfProbe_ = std::min(etaOfProbe_,2.5);
  }

  if(h_etaMax->GetEntries()==0.){
    h_etaMax->SetBinContent(1,etaOfProbe_);
    h_nbins->SetBinContent(1,nBins_);
  }

  //=======================================================
  // Retrieve the Transient Track Builder information
  //=======================================================

  edm::ESHandle<TransientTrackBuilder> theB_;
  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",theB_);
  double fBfield_=((*theB_).field()->inTesla(GlobalPoint(0.,0.,0.))).z();

  //=======================================================
  // Retrieve the Track information
  //=======================================================
  
  edm::Handle<TrackCollection>  trackCollectionHandle;
  iEvent.getByToken(theTrackCollectionToken, trackCollectionHandle);
  
  //=======================================================
  // Retrieve offline vartex information (only for reco)
  //=======================================================
 
  //edm::Handle<VertexCollection> vertices;
  edm::Handle<std::vector<Vertex> > vertices;

  try {
    iEvent.getByToken(theVertexCollectionToken, vertices);
  } catch ( cms::Exception& er ) {
    LogTrace("PrimaryVertexValidation")<<"caught std::exception "<<er.what()<<std::endl;  
  }

  std::vector<Vertex> vsorted = *(vertices);
  // sort the vertices by number of tracks in descending order
  // use chi2 as tiebreaker
  std::sort( vsorted.begin(), vsorted.end(), PrimaryVertexValidation::vtxSort );
  
  // skip events with no PV, this should not happen
  if( vsorted.size() == 0) return;
  // skip events failing vertex cut
  if( fabs(vsorted[0].z()) > vertexZMax_ ) return; 
  
  if ( vsorted[0].isValid() ) {
    xOfflineVertex_ = (vsorted)[0].x();
    yOfflineVertex_ = (vsorted)[0].y();
    zOfflineVertex_ = (vsorted)[0].z();

    xErrOfflineVertex_ = (vsorted)[0].xError();
    yErrOfflineVertex_ = (vsorted)[0].yError();
    zErrOfflineVertex_ = (vsorted)[0].zError();
  }

  h_xOfflineVertex->Fill(xOfflineVertex_);    
  h_yOfflineVertex->Fill(yOfflineVertex_);     
  h_zOfflineVertex->Fill(zOfflineVertex_);     
  h_xErrOfflineVertex->Fill(xErrOfflineVertex_);  
  h_yErrOfflineVertex->Fill(yErrOfflineVertex_);  
  h_zErrOfflineVertex->Fill(zErrOfflineVertex_);  

  unsigned int vertexCollectionSize = vsorted.size();
  int nvvertex = 0;
  
  for (unsigned int i=0; i<vertexCollectionSize; i++) {
    const Vertex& vertex = vsorted.at(i);
    if (vertex.isValid()) nvvertex++;
  }

  nOfflineVertices_ = nvvertex;
  h_nOfflineVertices->Fill(nvvertex);

  if ( vsorted.size() && useTracksFromRecoVtx_ ) {
   
    double sumpt    = 0;
    size_t ntracks  = 0;
    double chi2ndf  = 0.; 
    double chi2prob = 0.;

    if (!vsorted.at(0).isFake()) {
      
      Vertex pv = vsorted.at(0);
      
      ntracks  = pv.tracksSize();
      chi2ndf  = pv.normalizedChi2();
      chi2prob = TMath::Prob(pv.chi2(),(int)pv.ndof());
      
      h_recoVtxNtracks_->Fill(ntracks);   
      h_recoVtxChi2ndf_->Fill(chi2ndf);    
      h_recoVtxChi2Prob_->Fill(chi2prob);   
      
      for (Vertex::trackRef_iterator itrk = pv.tracks_begin();itrk != pv.tracks_end(); ++itrk) {
    double pt = (**itrk).pt();
    sumpt += pt*pt;
    
    const math::XYZPoint myVertex(pv.position().x(),pv.position().y(),pv.position().z());
    
    double dxyRes = (**itrk).dxy(myVertex);
    double dzRes  = (**itrk).dz(myVertex);
    
    double dxy_err = (**itrk).dxyError();
    double dz_err  = (**itrk).dzError();
    
    float trackphi = ((**itrk).phi())*(180/TMath::Pi());
    float tracketa = (**itrk).eta();
    
    for(int i=0; i<nBins_; i++){
      
      float phiF = (-TMath::Pi()+i*phiSect_)*(180/TMath::Pi());
      float phiL = (-TMath::Pi()+(i+1)*phiSect_)*(180/TMath::Pi());
      
      float etaF=-etaOfProbe_+i*etaSect_;
      float etaL=-etaOfProbe_+(i+1)*etaSect_;
      
      if(tracketa >= etaF && tracketa < etaL ){

        a_dxyEtaBiasResiduals[i]->Fill(dxyRes*cmToum);
        a_dzEtaBiasResiduals[i]->Fill(dzRes*cmToum); 
        n_dxyEtaBiasResiduals[i]->Fill((dxyRes)/dxy_err);
        n_dzEtaBiasResiduals[i]->Fill((dzRes)/dz_err);

      }

      if(trackphi >= phiF && trackphi < phiL ){ 

        a_dxyPhiBiasResiduals[i]->Fill(dxyRes*cmToum);
        a_dzPhiBiasResiduals[i]->Fill(dzRes*cmToum); 
        n_dxyPhiBiasResiduals[i]->Fill((dxyRes)/dxy_err);
        n_dzPhiBiasResiduals[i]->Fill((dzRes)/dz_err); 
        
        for(int j=0; j<nBins_; j++){
          
          float etaJ=-etaOfProbe_+j*etaSect_;
          float etaK=-etaOfProbe_+(j+1)*etaSect_;
          
          if(tracketa >= etaJ && tracketa < etaK ){
        
        a_dxyBiasResidualsMap[i][j]->Fill(dxyRes*cmToum); 
        a_dzBiasResidualsMap[i][j]->Fill(dzRes*cmToum);   
        
        n_dxyBiasResidualsMap[i][j]->Fill((dxyRes)/dxy_err); 
        n_dzBiasResidualsMap[i][j]->Fill((dzRes)/dz_err);  
        
          }
        }
      }     
    }
      }
      
      h_recoVtxSumPt_->Fill(sumpt);   

    }
  }

  //=======================================================
  // Retrieve Beamspot information
  //=======================================================

  BeamSpot beamSpot;
  edm::Handle<BeamSpot> beamSpotHandle;
  iEvent.getByToken(theBeamspotToken, beamSpotHandle);
    
  if ( beamSpotHandle.isValid() ) {
    beamSpot = *beamSpotHandle;
    BSx0_    = beamSpot.x0();
    BSy0_    = beamSpot.y0();
    BSz0_    = beamSpot.z0();
    Beamsigmaz_ = beamSpot.sigmaZ();    
    Beamdxdz_   = beamSpot.dxdz();       
    BeamWidthX_ = beamSpot.BeamWidthX();
    BeamWidthY_ = beamSpot.BeamWidthY();

    wxy2_=TMath::Power(BeamWidthX_,2)+TMath::Power(BeamWidthY_,2);

  } else {
    edm::LogWarning("PrimaryVertexValidation")<<"No BeamSpot found!";
  }

  h_BSx0->Fill(BSx0_);      
  h_BSy0->Fill(BSy0_);     
  h_BSz0->Fill(BSz0_);      
  h_Beamsigmaz->Fill(Beamsigmaz_);
  h_BeamWidthX->Fill(BeamWidthX_);
  h_BeamWidthY->Fill(BeamWidthY_);

  if(debug_)
    edm::LogInfo("PrimaryVertexValidation")<<"Beamspot x:" <<BSx0_<<" y:"<<BSy0_<<" z:"<<BSz0_;
   
  //=======================================================
  // Starts here ananlysis
  //=======================================================
  
  RunNumber_=iEvent.eventAuxiliary().run();
  h_runNumber->Fill(RunNumber_);

  if(h_runFromEvent->GetEntries()==0){
    h_runFromEvent->SetBinContent(1,RunNumber_);
  }

  LuminosityBlockNumber_=iEvent.eventAuxiliary().luminosityBlock();
  EventNumber_=iEvent.eventAuxiliary().id().event();
    
  if(debug_)
    edm::LogInfo("PrimaryVertexValidation")<<" looping over "<<trackCollectionHandle->size()<< "tracks";

  h_nTracks->Fill(trackCollectionHandle->size()); 

  //======================================================
  // Interface RECO tracks to vertex reconstruction
  //======================================================
 
  std::vector<TransientTrack> t_tks;
  unsigned int k = 0;   
  for(TrackCollection::const_iterator track = trackCollectionHandle->begin(); track!= trackCollectionHandle->end(); ++track, ++k){
  
    TransientTrack tt = theB_->build(&(*track));  
    tt.setBeamSpot(beamSpot);
    t_tks.push_back(tt);
  
  }
  
  if(debug_) {
    edm::LogInfo("PrimaryVertexValidation") << "Found: " << t_tks.size() << " reconstructed tracks";
  }
  
  //======================================================
  // select the tracks
  //======================================================

  std::vector<TransientTrack> seltks = theTrackFilter_->select(t_tks);
    
  //======================================================
  // clusterize tracks in Z
  //======================================================

  vector< vector<TransientTrack> > clusters = theTrackClusterizer_->clusterize(seltks);
  
  if (debug_){
    edm::LogInfo("PrimaryVertexValidation")<<" looping over: "<< clusters.size() << " clusters  from " << t_tks.size() << " selected tracks";
  }
  
  nClus_=clusters.size();  
  h_nClus->Fill(nClus_);

  //======================================================
  // Starts loop on clusters 
  //======================================================

  for (vector< vector<TransientTrack> >::const_iterator iclus = clusters.begin(); iclus != clusters.end(); iclus++) {

    nTracksPerClus_=0;

    unsigned int i = 0;   
    for(vector<TransientTrack>::const_iterator theTTrack = iclus->begin(); theTTrack!= iclus->end(); ++theTTrack, ++i)
      {
    if ( nTracks_ >= nMaxtracks_ ) {
      edm::LogError("PrimaryVertexValidation")<<" Warning - Number of tracks: " << nTracks_ << " , greater than " << nMaxtracks_;
      continue;
    }
    
    const Track & theTrack = theTTrack->track();

    pt_[nTracks_]       = theTrack.pt();
    p_[nTracks_]        = theTrack.p();
    nhits_[nTracks_]    = theTrack.numberOfValidHits();
    eta_[nTracks_]      = theTrack.eta();
    theta_[nTracks_]    = theTrack.theta();
    phi_[nTracks_]      = theTrack.phi();
    chi2_[nTracks_]     = theTrack.chi2();
    chi2ndof_[nTracks_] = theTrack.normalizedChi2();
    charge_[nTracks_]   = theTrack.charge();
    qoverp_[nTracks_]   = theTrack.qoverp();
    dz_[nTracks_]       = theTrack.dz();
    dxy_[nTracks_]      = theTrack.dxy();
    
    TrackBase::TrackQuality _trackQuality = TrackBase::qualityByName("highPurity");
    isHighPurity_[nTracks_] = theTrack.quality(_trackQuality);
    
    math::XYZPoint point(BSx0_,BSy0_,BSz0_);
    dxyBs_[nTracks_]    = theTrack.dxy(point);
    dzBs_[nTracks_]     = theTrack.dz(point);

    xPCA_[nTracks_]     = theTrack.vertex().x();
    yPCA_[nTracks_]     = theTrack.vertex().y();
    zPCA_[nTracks_]     = theTrack.vertex().z(); 
    
    //=======================================================
    // Retrieve rechit information
    //=======================================================  
    
    const reco::HitPattern& hits = theTrack.hitPattern();
    
    int nRecHit1D=0;
    int nRecHit2D=0;
    int nhitinTIB  = hits.numberOfValidStripTIBHits(); 
    int nhitinTOB  = hits.numberOfValidStripTOBHits(); 
    int nhitinTID  = hits.numberOfValidStripTIDHits(); 
    int nhitinTEC  = hits.numberOfValidStripTECHits();
    int nhitinBPIX = hits.numberOfValidPixelBarrelHits();
    int nhitinFPIX = hits.numberOfValidPixelEndcapHits();
    
    for (trackingRecHit_iterator iHit = theTTrack->recHitsBegin(); iHit != theTTrack->recHitsEnd(); ++iHit) {
      if((*iHit)->isValid()) {  
        
        if (this->isHit2D(**iHit)) {++nRecHit2D;}
        else {++nRecHit1D; }
      }
    } 
    
    nhits1D_[nTracks_]     = nRecHit1D;
    nhits2D_[nTracks_]     = nRecHit2D;
    nhitsBPIX_[nTracks_]   = nhitinBPIX;
    nhitsFPIX_[nTracks_]   = nhitinFPIX;
    nhitsTIB_[nTracks_]    = nhitinTIB;
    nhitsTID_[nTracks_]    = nhitinTID;
    nhitsTOB_[nTracks_]    = nhitinTOB;
    nhitsTEC_[nTracks_]    = nhitinTEC;
    
    //=======================================================
    // Good tracks for vertexing selection
    //=======================================================  

    bool pass = true;
    if(askFirstLayerHit_) pass = this->hasFirstLayerPixelHits((*theTTrack));
    if (pass && (theTrack.pt() >=ptOfProbe_) && fabs(theTrack.eta()) <= etaOfProbe_){
      isGoodTrack_[nTracks_]=1;
    }
      
    //=======================================================
    // Fit unbiased vertex
    //=======================================================
    
    vector<TransientTrack> theFinalTracks;
    theFinalTracks.clear();

    for(vector<TransientTrack>::const_iterator tk = iclus->begin(); tk!= iclus->end(); ++tk){
      
      pass = this->hasFirstLayerPixelHits((*tk));
      if (pass){
        if( tk == theTTrack ) continue;
        else {
          theFinalTracks.push_back((*tk));
        }
      }
    }
    
    if(theFinalTracks.size() > 1){
        
      if(debug_)
        edm::LogInfo("PrimaryVertexValidation")<<"Transient Track Collection size: "<<theFinalTracks.size();
      try{
          
        VertexFitter<5>* theFitter = new AdaptiveVertexFitter;
        TransientVertex theFittedVertex = theFitter->vertex(theFinalTracks);

        //AdaptiveVertexFitter* theFitter = new AdaptiveVertexFitter;
        //TransientVertex theFittedVertex = theFitter->vertex(theFinalTracks,beamSpot);  // if you want the beam constraint

        double totalTrackWeights=0;
        if(theFittedVertex.isValid ()){

          
          if(theFittedVertex.hasTrackWeight()){
        for(size_t rtracks= 0; rtracks < theFinalTracks.size(); rtracks++){
          sumOfWeightsUnbiasedVertex_[nTracks_] += theFittedVertex.trackWeight(theFinalTracks[rtracks]);
          totalTrackWeights+= theFittedVertex.trackWeight(theFinalTracks[rtracks]);
          h_fitVtxTrackWeights_->Fill(theFittedVertex.trackWeight(theFinalTracks[rtracks]));
        }
          }
          
          h_fitVtxTrackAverageWeight_->Fill(totalTrackWeights/theFinalTracks.size());

          const math::XYZPoint theRecoVertex(xOfflineVertex_,yOfflineVertex_,zOfflineVertex_);      
          const math::XYZPoint myVertex(theFittedVertex.position().x(),theFittedVertex.position().y(),theFittedVertex.position().z());

          const Vertex vertex = theFittedVertex;
          fillTrackHistos(hDA,"all",&(*theTTrack),vertex,beamSpot,fBfield_);

          hasRecVertex_[nTracks_]    = 1;
          xUnbiasedVertex_[nTracks_] = theFittedVertex.position().x();
          yUnbiasedVertex_[nTracks_] = theFittedVertex.position().y();
          zUnbiasedVertex_[nTracks_] = theFittedVertex.position().z();
          
          chi2normUnbiasedVertex_[nTracks_] = theFittedVertex.normalisedChiSquared();
          chi2UnbiasedVertex_[nTracks_]     = theFittedVertex.totalChiSquared();
          DOFUnbiasedVertex_[nTracks_]      = theFittedVertex.degreesOfFreedom();   
          chi2ProbUnbiasedVertex_[nTracks_] = TMath::Prob(theFittedVertex.totalChiSquared(),(int)theFittedVertex.degreesOfFreedom());
          tracksUsedForVertexing_[nTracks_] = theFinalTracks.size();

          h_fitVtxNtracks_->Fill(theFinalTracks.size());        
          h_fitVtxChi2_->Fill(theFittedVertex.totalChiSquared());
          h_fitVtxNdof_->Fill(theFittedVertex.degreesOfFreedom());
          h_fitVtxChi2ndf_->Fill(theFittedVertex.normalisedChiSquared());        
          h_fitVtxChi2Prob_->Fill(TMath::Prob(theFittedVertex.totalChiSquared(),(int)theFittedVertex.degreesOfFreedom()));       
                
          // from my Vertex
          double dxyFromMyVertex = theTrack.dxy(myVertex);
          double dzFromMyVertex  = theTrack.dz(myVertex);

          GlobalPoint vert(theFittedVertex.position().x(),theFittedVertex.position().y(),theFittedVertex.position().z());

          //FreeTrajectoryState theTrackNearVertex = (*theTTrack).trajectoryStateClosestToPoint(vert).theState();
          //double dz_err = sqrt(theFittedVertex.positionError().czz() + theTrackNearVertex.cartesianError().position().czz());      
          //double dz_err = hypot(theTrack.dzError(),theFittedVertex.positionError().czz());
          
          double dz_err = sqrt(std::pow(theTrack.dzError(), 2) + theFittedVertex.positionError().czz());


          // PV2D 
          std::pair<bool,Measurement1D> s_ip2dpv = signedTransverseImpactParameter(*theTTrack,
                                               GlobalVector(theTrack.px(),
                                                    theTrack.py(),
                                                    theTrack.pz()),
                                               theFittedVertex);            
          
          double s_ip2dpv_corr = s_ip2dpv.second.value();
          double s_ip2dpv_err  = s_ip2dpv.second.error();
          
          // PV3D
          std::pair<bool, Measurement1D> s_ip3dpv = signedImpactParameter3D(*theTTrack,         
                                        GlobalVector(theTrack.px(),  
                                                 theTrack.py(),  
                                                 theTrack.pz()), 
                                        theFittedVertex);            
          
          double s_ip3dpv_corr = s_ip3dpv.second.value();
          double s_ip3dpv_err  = s_ip3dpv.second.error();

          // PV3D absolute
          std::pair<bool,Measurement1D> ip3dpv = absoluteImpactParameter3D(*theTTrack,theFittedVertex);
          double ip3d_corr = ip3dpv.second.value(); 
          double ip3d_err  = ip3dpv.second.error(); 
          
          // with respect to any specified vertex, such as primary vertex
          TrajectoryStateClosestToPoint traj = (*theTTrack).trajectoryStateClosestToPoint(vert);

          GlobalPoint refPoint = traj.position();
          GlobalPoint cPToVtx  = traj.theState().position();

          float my_dx = refPoint.x() - myVertex.x();
          float my_dy = refPoint.y() - myVertex.y();

          float my_dx2 = cPToVtx.x() - myVertex.x();
          float my_dy2 = cPToVtx.y() - myVertex.y();

          float my_dxy = std::sqrt(my_dx*my_dx + my_dy*my_dy);

          double d0 = traj.perigeeParameters().transverseImpactParameter();
          //double d0_error = traj.perigeeError().transverseImpactParameterError();
          double z0 = traj.perigeeParameters().longitudinalImpactParameter();
          double z0_error = traj.perigeeError().longitudinalImpactParameterError();

          if(debug_){
        edm::LogInfo("PrimaryVertexValidation")<< "my_dx:"   << my_dx  
                               << " my_dy:"  << my_dy  
                               << " my_dxy:" << my_dxy
                               << " my_dx2:" << my_dx2  
                               << " my_dy2:" << my_dy2
                               << " d0: "    << d0
                               << " dxyFromVtx:" << dxyFromMyVertex << "\n"
                               << " ============================== "<< "\n"
                               << "diff1:"    << std::abs(d0) - std::abs(my_dxy) << "\n"
                               << "diff2:"    << std::abs(d0) - std::abs(dxyFromMyVertex) << "\n"
                               << "diff3:"    << (my_dx - my_dx2) << " " << (my_dy - my_dy2) << "\n" 
                               << std::endl;    
          }

          // define IPs
         
          dxyFromMyVertex_[nTracks_]      = dxyFromMyVertex; 
          dxyErrorFromMyVertex_[nTracks_] = s_ip2dpv_err;
          IPTsigFromMyVertex_[nTracks_]   = dxyFromMyVertex/s_ip2dpv_err;
 
          dzFromMyVertex_[nTracks_]       = dzFromMyVertex;  
          dzErrorFromMyVertex_[nTracks_]  = dz_err;
          IPLsigFromMyVertex_[nTracks_]   = dzFromMyVertex/dz_err;

          d3DFromMyVertex_[nTracks_]      = ip3d_corr;
          d3DErrorFromMyVertex_[nTracks_] = ip3d_err;
          IP3DsigFromMyVertex_[nTracks_]  = (ip3d_corr/ip3d_err);
                  
          // fill directly the histograms of residuals
          
          float trackphi = (theTrack.phi())*(180/TMath::Pi());
          float tracketa = theTrack.eta();
          float trackpt  = theTrack.pt();

          
          // filling the pT-binned distributions

          for(int ipTBin=0; ipTBin<nPtBins_; ipTBin++){
        
        float pTF = mypT_bins_[ipTBin];
        float pTL = mypT_bins_[ipTBin+1];
        
        if(debug_)
          edm::LogInfo("PrimaryVertexValidation")<<"ipTBin:"<<ipTBin<< " "<<mypT_bins_[ipTBin]<< " < pT < "<<mypT_bins_[ipTBin+1]<<std::endl;
        
        if( fabs(tracketa)<1.5 && (trackpt >= pTF && trackpt < pTL) ){
          
          if(debug_)
            edm::LogInfo("PrimaryVertexValidation")<<"passes this cut: "<<mypT_bins_[ipTBin]<<std::endl;
          fillByIndex(h_dxy_pT_,ipTBin,dxyFromMyVertex*cmToum);
          fillByIndex(h_dz_pT_,ipTBin,dzFromMyVertex*cmToum);
          fillByIndex(h_norm_dxy_pT_,ipTBin,dxyFromMyVertex/s_ip2dpv_err);
          fillByIndex(h_norm_dz_pT_,ipTBin,dzFromMyVertex/dz_err);
          
          if(fabs(tracketa)<1.){
            
            if(debug_)
              edm::LogInfo("PrimaryVertexValidation")<<"passes tight eta cut: "<<mypT_bins_[ipTBin]<<std::endl;
            fillByIndex(h_dxy_Central_pT_,ipTBin,dxyFromMyVertex*cmToum);
            fillByIndex(h_dz_Central_pT_,ipTBin,dzFromMyVertex*cmToum);
            fillByIndex(h_norm_dxy_Central_pT_,ipTBin,dxyFromMyVertex/s_ip2dpv_err);
            fillByIndex(h_norm_dz_Central_pT_,ipTBin,dzFromMyVertex/dz_err);
          }
        }
          }
          
          // checks on the probe track quality
          if(trackpt >= ptOfProbe_ && fabs(tracketa)<= etaOfProbe_){

        std::pair<bool,bool> pixelOcc = pixelHitsCheck((*theTTrack));

        /*
          if(pixelOcc.first == true)
          std::cout<<"has BPIx hits"<<std::endl;
          if(pixelOcc.second == true)
          std::cout<<"has FPix hits"<<std::endl;
        */        

        if(!doBPix_ && (pixelOcc.first == true))  continue;
        if(!doFPix_ && (pixelOcc.second == true)) continue;
    
        //std::cout<<"track passed"<<std::endl;
    
        fillTrackHistos(hDA,"sel",&(*theTTrack),vertex,beamSpot,fBfield_);

        // probe checks
        h_probePt_->Fill(theTrack.pt());
        h_probeP_->Fill(theTrack.p());
        h_probeEta_->Fill(theTrack.eta());
        h_probePhi_->Fill(theTrack.phi());
        h2_probeEtaPhi_->Fill(theTrack.eta(),theTrack.phi());
        h2_probeEtaPt_->Fill(theTrack.eta(),theTrack.pt());

        h_probeChi2_->Fill(theTrack.chi2());
        h_probeNormChi2_->Fill(theTrack.normalizedChi2());
        h_probeCharge_->Fill(theTrack.charge());
        h_probeQoverP_->Fill(theTrack.qoverp());
        h_probeHits_->Fill(theTrack.numberOfValidHits());       
        h_probeHits1D_->Fill(nRecHit1D);
        h_probeHits2D_->Fill(nRecHit2D);
        h_probeHitsInTIB_->Fill(nhitinBPIX);
        h_probeHitsInTOB_->Fill(nhitinFPIX);
        h_probeHitsInTID_->Fill(nhitinTIB);
        h_probeHitsInTEC_->Fill(nhitinTID);
        h_probeHitsInBPIX_->Fill(nhitinTOB);
        h_probeHitsInFPIX_->Fill(nhitinTEC);
        
        float dxyRecoV = theTrack.dz(theRecoVertex);
        float dzRecoV  = theTrack.dxy(theRecoVertex);
        float dxysigmaRecoV = TMath::Sqrt(theTrack.d0Error()*theTrack.d0Error()+xErrOfflineVertex_*yErrOfflineVertex_);
        float dzsigmaRecoV  = TMath::Sqrt(theTrack.dzError()*theTrack.dzError()+zErrOfflineVertex_*zErrOfflineVertex_);

        double zTrack=(theTTrack->stateAtBeamLine().trackStateAtPCA()).position().z();
        double zVertex=theFittedVertex.position().z();
        double tantheta=tan((theTTrack->stateAtBeamLine().trackStateAtPCA()).momentum().theta());

        double dz2= pow(theTrack.dzError(),2)+wxy2_/pow(tantheta,2);
        double restrkz   = zTrack-zVertex;
        double pulltrkz  = (zTrack-zVertex)/TMath::Sqrt(dz2);

        h_probedxyRecoV_->Fill(dxyRecoV);
        h_probedzRecoV_->Fill(dzRecoV);
    
        h_probedzRefitV_->Fill(dxyFromMyVertex);
        h_probedxyRefitV_->Fill(dzFromMyVertex);

        h_probed0RefitV_->Fill(d0);
        h_probez0RefitV_->Fill(z0);

        h_probesignIP2DRefitV_->Fill(s_ip2dpv_corr);
        h_probed3DRefitV_->Fill(ip3d_corr);
        h_probereszRefitV_->Fill(restrkz);

        h_probeRecoVSigZ_->Fill(dzRecoV/dzsigmaRecoV);
        h_probeRecoVSigXY_->Fill(dxyRecoV/dxysigmaRecoV);   
        h_probeRefitVSigZ_->Fill(dzFromMyVertex/dz_err);
        h_probeRefitVSigXY_->Fill(dxyFromMyVertex/s_ip2dpv_err); 
        h_probeRefitVSig3D_->Fill(ip3d_corr/ip3d_err);
        h_probeRefitVLogSig3D_->Fill(log10(ip3d_corr/ip3d_err));
        h_probeRefitVSigResZ_->Fill(pulltrkz);
        
        a_dxyVsPhi->Fill(trackphi,dxyFromMyVertex*cmToum);
        a_dzVsPhi->Fill(trackphi,z0*cmToum);  
        n_dxyVsPhi->Fill(trackphi,dxyFromMyVertex/s_ip2dpv_err); 
        n_dzVsPhi->Fill(trackphi,z0/z0_error);
  
        a_dxyVsEta->Fill(tracketa,dxyFromMyVertex*cmToum); 
        a_dzVsEta->Fill(tracketa,z0*cmToum);  
        n_dxyVsEta->Fill(tracketa,dxyFromMyVertex/s_ip2dpv_err); 
        n_dzVsEta->Fill(tracketa,z0/z0_error); 
 
        // filling the binned distributions
        for(int i=0; i<nBins_; i++){
          
          float phiF = (-TMath::Pi()+i*phiSect_)*(180/TMath::Pi());
          float phiL = (-TMath::Pi()+(i+1)*phiSect_)*(180/TMath::Pi());
          
          float etaF=-etaOfProbe_+i*etaSect_;
          float etaL=-etaOfProbe_+(i+1)*etaSect_;

          if(tracketa >= etaF && tracketa < etaL ){

            a_dxyEtaResiduals[i]->Fill(dxyFromMyVertex*cmToum);
            a_dxEtaResiduals[i]->Fill(my_dx*cmToum);
            a_dyEtaResiduals[i]->Fill(my_dy*cmToum);
            a_dzEtaResiduals[i]->Fill(dzFromMyVertex*cmToum);   
            n_dxyEtaResiduals[i]->Fill(dxyFromMyVertex/s_ip2dpv_err);
            n_dzEtaResiduals[i]->Fill(dzFromMyVertex/dz_err);       
            a_IP2DEtaResiduals[i]->Fill(s_ip2dpv_corr*cmToum);
            n_IP2DEtaResiduals[i]->Fill(s_ip2dpv_corr/s_ip2dpv_err);
            a_reszEtaResiduals[i]->Fill(restrkz*cmToum);
            n_reszEtaResiduals[i]->Fill(pulltrkz);
            a_d3DEtaResiduals[i]->Fill(ip3d_corr*cmToum);   
            n_d3DEtaResiduals[i]->Fill(ip3d_corr/ip3d_err);
            a_IP3DEtaResiduals[i]->Fill(s_ip3dpv_corr*cmToum);
            n_IP3DEtaResiduals[i]->Fill(s_ip3dpv_corr/s_ip3dpv_err);

          }
              
          if(trackphi >= phiF && trackphi < phiL ){
            a_dxyPhiResiduals[i]->Fill(dxyFromMyVertex*cmToum);
            a_dxPhiResiduals[i]->Fill(my_dx*cmToum);
            a_dyPhiResiduals[i]->Fill(my_dy*cmToum);
            a_dzPhiResiduals[i]->Fill(dzFromMyVertex*cmToum); 
            n_dxyPhiResiduals[i]->Fill(dxyFromMyVertex/s_ip2dpv_err);
            n_dzPhiResiduals[i]->Fill(dzFromMyVertex/dz_err); 
            a_IP2DPhiResiduals[i]->Fill(s_ip2dpv_corr*cmToum);
            n_IP2DPhiResiduals[i]->Fill(s_ip2dpv_corr/s_ip2dpv_err); 
            a_reszPhiResiduals[i]->Fill(restrkz*cmToum);
            n_reszPhiResiduals[i]->Fill(pulltrkz);
            a_d3DPhiResiduals[i]->Fill(ip3d_corr*cmToum);   
            n_d3DPhiResiduals[i]->Fill(ip3d_corr/ip3d_err);
            a_IP3DPhiResiduals[i]->Fill(s_ip3dpv_corr*cmToum);
            n_IP3DPhiResiduals[i]->Fill(s_ip3dpv_corr/s_ip3dpv_err);

            for(int j=0; j<nBins_; j++){

              float etaJ=-etaOfProbe_+j*etaSect_;
              float etaK=-etaOfProbe_+(j+1)*etaSect_;

              if(tracketa >= etaJ && tracketa < etaK ){
            a_dxyResidualsMap[i][j]->Fill(dxyFromMyVertex*cmToum); 
            a_dzResidualsMap[i][j]->Fill(dzFromMyVertex*cmToum);        
            n_dxyResidualsMap[i][j]->Fill(dxyFromMyVertex/s_ip2dpv_err); 
            n_dzResidualsMap[i][j]->Fill(dzFromMyVertex/dz_err);  
            a_d3DResidualsMap[i][j]->Fill(ip3d_corr*cmToum);   
            n_d3DResidualsMap[i][j]->Fill(ip3d_corr/ip3d_err); 

              }
            }
          }     
        }
          }
              
          if(debug_){
        edm::LogInfo("PrimaryVertexValidation")<<" myVertex.x()= "<<myVertex.x()<<"\n"
                               <<" myVertex.y()= "<<myVertex.y()<<" \n"
                               <<" myVertex.z()= "<<myVertex.z()<<" \n"
                               <<" theTrack.dz(myVertex)= "<<theTrack.dz(myVertex)<<" \n"
                               <<" zPCA -myVertex.z() = "<<(theTrack.vertex().z() -myVertex.z());
        
          }// ends if debug_
        } // ends if the fitted vertex is Valid

        delete theFitter;

      }  catch ( cms::Exception& er ) {
        LogTrace("PrimaryVertexValidation")<<"caught std::exception "<<er.what()<<std::endl;
      }
        
    } //ends if theFinalTracks.size() > 2
    
    else {
      if(debug_)
        edm::LogInfo("PrimaryVertexValidation")<<"Not enough tracks to make a vertex.  Returns no vertex info";
    }
      
    ++nTracks_;  
    ++nTracksPerClus_;

    if(debug_)
      edm::LogInfo("PrimaryVertexValidation")<<"Track "<<i<<" : pT = "<<theTrack.pt();
      
      }// for loop on tracks
  } // for loop on track clusters
  
  // Fill the TTree if needed

  if(storeNtuple_){  
    rootTree_->Fill();
  }

  
} 

void PrimaryVertexValidation::beginJob ( void  )

privatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 994 of file PrimaryVertexValidation.cc.

References a_d3DEtaResiduals, a_d3DPhiResiduals, a_d3DResidualsMap, a_dxEtaResiduals, a_dxPhiResiduals, a_dxyBiasResidualsMap, a_dxyEtaBiasResiduals, a_dxyEtaMADBiasTrend, a_dxyEtaMADTrend, a_dxyEtaMeanBiasTrend, a_dxyEtaMeanTrend, a_dxyEtaMedianBiasTrend, a_dxyEtaMedianTrend, a_dxyEtaResiduals, a_dxyEtaWidthBiasTrend, a_dxyEtaWidthTrend, a_dxyMeanBiasMap, a_dxyMeanMap, a_dxyPhiBiasResiduals, a_dxyPhiMADBiasTrend, a_dxyPhiMADTrend, a_dxyPhiMeanBiasTrend, a_dxyPhiMeanTrend, a_dxyPhiMedianBiasTrend, a_dxyPhiMedianTrend, a_dxyPhiResiduals, a_dxyPhiWidthBiasTrend, a_dxyPhiWidthTrend, a_dxypTCentralMeanTrend, a_dxypTCentralWidthTrend, a_dxypTMeanTrend, a_dxypTWidthTrend, a_dxyResidualsMap, a_dxyVsEta, a_dxyVsPhi, a_dxyWidthBiasMap, a_dxyWidthMap, a_dyEtaResiduals, a_dyPhiResiduals, a_dzBiasResidualsMap, a_dzEtaBiasResiduals, a_dzEtaMADBiasTrend, a_dzEtaMADTrend, a_dzEtaMeanBiasTrend, a_dzEtaMeanTrend, a_dzEtaMedianBiasTrend, a_dzEtaMedianTrend, a_dzEtaResiduals, a_dzEtaWidthBiasTrend, a_dzEtaWidthTrend, a_dzMeanBiasMap, a_dzMeanMap, a_dzPhiBiasResiduals, a_dzPhiMADBiasTrend, a_dzPhiMADTrend, a_dzPhiMeanBiasTrend, a_dzPhiMeanTrend, a_dzPhiMedianBiasTrend, a_dzPhiMedianTrend, a_dzPhiResiduals, a_dzPhiWidthBiasTrend, a_dzPhiWidthTrend, a_dzpTCentralMeanTrend, a_dzpTCentralWidthTrend, a_dzpTMeanTrend, a_dzpTWidthTrend, a_dzResidualsMap, a_dzVsEta, a_dzVsPhi, a_dzWidthBiasMap, a_dzWidthMap, a_IP2DEtaResiduals, a_IP2DPhiResiduals, a_IP3DEtaResiduals, a_IP3DPhiResiduals, a_reszEtaResiduals, a_reszPhiResiduals, Beamdxdz_, Beamsigmaz_, BeamWidthX_, BeamWidthY_, bookResidualsHistogram(), bookVertexHistograms(), BSx0_, BSy0_, BSz0_, charge_, chi2_, chi2ndof_, chi2normUnbiasedVertex_, chi2ProbUnbiasedVertex_, chi2UnbiasedVertex_, d3DFromMyVertex_, DOFUnbiasedVertex_, dxy_, dxyBs_, dxyErrorFromMyVertex_, dxyFromMyVertex_, dz_, dzBs_, dzErrorFromMyVertex_, dzFromMyVertex_, eta_, etaOfProbe_, etaSect_, EventNumber_, h2_probeEtaPhi_, h2_probeEtaPt_, h_Beamsigmaz, h_BeamWidthX, h_BeamWidthY, h_BSx0, h_BSy0, h_BSz0, h_dxy_Central_pT_, h_dxy_pT_, h_dz_Central_pT_, h_dz_pT_, h_etaMax, h_fitVtxChi2_, h_fitVtxChi2ndf_, h_fitVtxChi2Prob_, h_fitVtxNdof_, h_fitVtxNtracks_, h_fitVtxTrackAverageWeight_, h_fitVtxTrackWeights_, h_lumiFromConfig, h_nbins, h_nClus, h_nOfflineVertices, h_norm_dxy_Central_pT_, h_norm_dxy_pT_, h_norm_dz_Central_pT_, h_norm_dz_pT_, h_nTracks, h_probeCharge_, h_probeChi2_, h_probed0RefitV_, h_probed3DRefitV_, h_probedxyRecoV_, h_probedxyRefitV_, h_probedzRecoV_, h_probedzRefitV_, h_probeEta_, h_probeHits1D_, h_probeHits2D_, h_probeHits_, h_probeHitsInBPIX_, h_probeHitsInFPIX_, h_probeHitsInTEC_, h_probeHitsInTIB_, h_probeHitsInTID_, h_probeHitsInTOB_, h_probeNormChi2_, h_probeP_, h_probePhi_, h_probePt_, h_probeQoverP_, h_probeRecoVSigXY_, h_probeRecoVSigZ_, h_probeRefitVLogSig3D_, h_probeRefitVSig3D_, h_probeRefitVSigResZ_, h_probeRefitVSigXY_, h_probeRefitVSigZ_, h_probereszRefitV_, h_probesignIP2DRefitV_, h_probez0RefitV_, h_recoVtxChi2ndf_, h_recoVtxChi2Prob_, h_recoVtxNtracks_, h_recoVtxSumPt_, h_runFromConfig, h_runFromEvent, h_runNumber, h_xErrOfflineVertex, h_xOfflineVertex, h_yErrOfflineVertex, h_yOfflineVertex, h_zErrOfflineVertex, h_zOfflineVertex, hasRecVertex_, hDA, mps_fire::i, intLumi_, IP3DsigFromMyVertex_, IPLsigFromMyVertex_, IPTsigFromMyVertex_, isGoodTrack_, isHighPurity_, lightNtupleSwitch_, LuminosityBlockNumber_, TFileDirectory::make(), TFileService::make(), TFileService::mkdir(), mypT_bins_, n_d3DEtaResiduals, n_d3DPhiResiduals, n_d3DResidualsMap, n_dxyBiasResidualsMap, n_dxyEtaBiasResiduals, n_dxyEtaMADBiasTrend, n_dxyEtaMADTrend, n_dxyEtaMeanBiasTrend, n_dxyEtaMeanTrend, n_dxyEtaMedianBiasTrend, n_dxyEtaMedianTrend, n_dxyEtaResiduals, n_dxyEtaWidthBiasTrend, n_dxyEtaWidthTrend, n_dxyMeanBiasMap, n_dxyMeanMap, n_dxyPhiBiasResiduals, n_dxyPhiMADBiasTrend, n_dxyPhiMADTrend, n_dxyPhiMeanBiasTrend, n_dxyPhiMeanTrend, n_dxyPhiMedianBiasTrend, n_dxyPhiMedianTrend, n_dxyPhiResiduals, n_dxyPhiWidthBiasTrend, n_dxyPhiWidthTrend, n_dxypTCentralMeanTrend, n_dxypTCentralWidthTrend, n_dxypTMeanTrend, n_dxypTWidthTrend, n_dxyResidualsMap, n_dxyVsEta, n_dxyVsPhi, n_dxyWidthBiasMap, n_dxyWidthMap, n_dzBiasResidualsMap, n_dzEtaBiasResiduals, n_dzEtaMADBiasTrend, n_dzEtaMADTrend, n_dzEtaMeanBiasTrend, n_dzEtaMeanTrend, n_dzEtaMedianBiasTrend, n_dzEtaMedianTrend, n_dzEtaResiduals, n_dzEtaWidthBiasTrend, n_dzEtaWidthTrend, n_dzMeanBiasMap, n_dzMeanMap, n_dzPhiBiasResiduals, n_dzPhiMADBiasTrend, n_dzPhiMADTrend, n_dzPhiMeanBiasTrend, n_dzPhiMeanTrend, n_dzPhiMedianBiasTrend, n_dzPhiMedianTrend, n_dzPhiResiduals, n_dzPhiWidthBiasTrend, n_dzPhiWidthTrend, n_dzpTCentralMeanTrend, n_dzpTCentralWidthTrend, n_dzpTMeanTrend, n_dzpTWidthTrend, n_dzResidualsMap, n_dzVsEta, n_dzVsPhi, n_dzWidthBiasMap, n_dzWidthMap, n_IP2DEtaResiduals, n_IP2DPhiResiduals, n_IP3DEtaResiduals, n_IP3DPhiResiduals, n_reszEtaResiduals, n_reszPhiResiduals, nBins_, nClus_, Nevt_, nhits1D_, nhits2D_, nhits_, nhitsBPIX_, nhitsFPIX_, nhitsTEC_, nhitsTIB_, nhitsTID_, nhitsTOB_, nOfflineVertices_, nPtBins_, nTracks_, nTracksPerClus_, p_, phi_, phiSect_, Pi, pt_, qoverp_, alignCSCRings::r, rootTree_, runControlNumbers_, RunNumber_, sumOfWeightsUnbiasedVertex_, theta_, tracksUsedForVertexing_, useTracksFromRecoVtx_, xOfflineVertex_, xPCA_, xUnbiasedVertex_, yOfflineVertex_, yPCA_, yUnbiasedVertex_, zOfflineVertex_, zPCA_, and zUnbiasedVertex_.

{
  edm::LogInfo("PrimaryVertexValidation") 
    <<"######################################\n"
    <<"Begin Job \n" 
    <<"######################################";

  // Define TTree for output
  Nevt_    = 0;
  
  //  rootFile_ = new TFile(filename_.c_str(),"recreate");
  edm::Service<TFileService> fs;
  rootTree_ = fs->make<TTree>("tree","PV Validation tree");
 
  // Track Paramters 

  if(lightNtupleSwitch_){ 

    rootTree_->Branch("EventNumber",&EventNumber_,"EventNumber/i");
    rootTree_->Branch("RunNumber",&RunNumber_,"RunNumber/i");
    rootTree_->Branch("LuminosityBlockNumber",&LuminosityBlockNumber_,"LuminosityBlockNumber/i");
    rootTree_->Branch("nOfflineVertices",&nOfflineVertices_,"nOfflineVertices/I");
    rootTree_->Branch("nTracks",&nTracks_,"nTracks/I");
    rootTree_->Branch("phi",&phi_,"phi[nTracks]/D");
    rootTree_->Branch("eta",&eta_,"eta[nTracks]/D");
    rootTree_->Branch("pt",&pt_,"pt[nTracks]/D");
    rootTree_->Branch("dxyFromMyVertex",&dxyFromMyVertex_,"dxyFromMyVertex[nTracks]/D");
    rootTree_->Branch("dzFromMyVertex",&dzFromMyVertex_,"dzFromMyVertex[nTracks]/D"); 
    rootTree_->Branch("d3DFromMyVertex",&d3DFromMyVertex_,"d3DFromMyVertex[nTracks]/D"); 
    rootTree_->Branch("IPTsigFromMyVertex",&IPTsigFromMyVertex_,"IPTsigFromMyVertex_[nTracks]/D");    
    rootTree_->Branch("IPLsigFromMyVertex",&IPLsigFromMyVertex_,"IPLsigFromMyVertex_[nTracks]/D"); 
    rootTree_->Branch("IP3DsigFromMyVertex",&IP3DsigFromMyVertex_,"IP3DsigFromMyVertex_[nTracks]/D"); 
    rootTree_->Branch("hasRecVertex",&hasRecVertex_,"hasRecVertex[nTracks]/I");
    rootTree_->Branch("isGoodTrack",&isGoodTrack_,"isGoodTrack[nTracks]/I");
    rootTree_->Branch("isHighPurity",&isHighPurity_,"isHighPurity_[nTracks]/I");
    
  } else {
    
    rootTree_->Branch("nTracks",&nTracks_,"nTracks/I");
    rootTree_->Branch("nTracksPerClus",&nTracksPerClus_,"nTracksPerClus/I");
    rootTree_->Branch("nClus",&nClus_,"nClus/I");
    rootTree_->Branch("xOfflineVertex",&xOfflineVertex_,"xOfflineVertex/D");
    rootTree_->Branch("yOfflineVertex",&yOfflineVertex_,"yOfflineVertex/D");
    rootTree_->Branch("zOfflineVertex",&zOfflineVertex_,"zOfflineVertex/D");
    rootTree_->Branch("BSx0",&BSx0_,"BSx0/D");
    rootTree_->Branch("BSy0",&BSy0_,"BSy0/D");
    rootTree_->Branch("BSz0",&BSz0_,"BSz0/D");
    rootTree_->Branch("Beamsigmaz",&Beamsigmaz_,"Beamsigmaz/D");
    rootTree_->Branch("Beamdxdz",&Beamdxdz_,"Beamdxdz/D");
    rootTree_->Branch("BeamWidthX",&BeamWidthX_,"BeamWidthX/D");
    rootTree_->Branch("BeamWidthY",&BeamWidthY_,"BeamWidthY/D");
    rootTree_->Branch("pt",&pt_,"pt[nTracks]/D");
    rootTree_->Branch("p",&p_,"p[nTracks]/D");
    rootTree_->Branch("nhits",&nhits_,"nhits[nTracks]/I");
    rootTree_->Branch("nhits1D",&nhits1D_,"nhits1D[nTracks]/I");
    rootTree_->Branch("nhits2D",&nhits2D_,"nhits2D[nTracks]/I");
    rootTree_->Branch("nhitsBPIX",&nhitsBPIX_,"nhitsBPIX[nTracks]/I");
    rootTree_->Branch("nhitsFPIX",&nhitsFPIX_,"nhitsFPIX[nTracks]/I");
    rootTree_->Branch("nhitsTIB",&nhitsTIB_,"nhitsTIB[nTracks]/I");
    rootTree_->Branch("nhitsTID",&nhitsTID_,"nhitsTID[nTracks]/I");
    rootTree_->Branch("nhitsTOB",&nhitsTOB_,"nhitsTOB[nTracks]/I");
    rootTree_->Branch("nhitsTEC",&nhitsTEC_,"nhitsTEC[nTracks]/I");
    rootTree_->Branch("eta",&eta_,"eta[nTracks]/D");
    rootTree_->Branch("theta",&theta_,"theta[nTracks]/D");
    rootTree_->Branch("phi",&phi_,"phi[nTracks]/D");
    rootTree_->Branch("chi2",&chi2_,"chi2[nTracks]/D");
    rootTree_->Branch("chi2ndof",&chi2ndof_,"chi2ndof[nTracks]/D");
    rootTree_->Branch("charge",&charge_,"charge[nTracks]/I");
    rootTree_->Branch("qoverp",&qoverp_,"qoverp[nTracks]/D");
    rootTree_->Branch("dz",&dz_,"dz[nTracks]/D");
    rootTree_->Branch("dxy",&dxy_,"dxy[nTracks]/D");
    rootTree_->Branch("dzBs",&dzBs_,"dzBs[nTracks]/D");
    rootTree_->Branch("dxyBs",&dxyBs_,"dxyBs[nTracks]/D");
    rootTree_->Branch("xPCA",&xPCA_,"xPCA[nTracks]/D");
    rootTree_->Branch("yPCA",&yPCA_,"yPCA[nTracks]/D");
    rootTree_->Branch("zPCA",&zPCA_,"zPCA[nTracks]/D");
    rootTree_->Branch("xUnbiasedVertex",&xUnbiasedVertex_,"xUnbiasedVertex[nTracks]/D");
    rootTree_->Branch("yUnbiasedVertex",&yUnbiasedVertex_,"yUnbiasedVertex[nTracks]/D");
    rootTree_->Branch("zUnbiasedVertex",&zUnbiasedVertex_,"zUnbiasedVertex[nTracks]/D");
    rootTree_->Branch("chi2normUnbiasedVertex",&chi2normUnbiasedVertex_,"chi2normUnbiasedVertex[nTracks]/F");
    rootTree_->Branch("chi2UnbiasedVertex",&chi2UnbiasedVertex_,"chi2UnbiasedVertex[nTracks]/F");
    rootTree_->Branch("DOFUnbiasedVertex",&DOFUnbiasedVertex_," DOFUnbiasedVertex[nTracks]/F");
    rootTree_->Branch("chi2ProbUnbiasedVertex",&chi2ProbUnbiasedVertex_,"chi2ProbUnbiasedVertex[nTracks]/F");
    rootTree_->Branch("sumOfWeightsUnbiasedVertex",&sumOfWeightsUnbiasedVertex_,"sumOfWeightsUnbiasedVertex[nTracks]/F");
    rootTree_->Branch("tracksUsedForVertexing",&tracksUsedForVertexing_,"tracksUsedForVertexing[nTracks]/I");
    rootTree_->Branch("dxyFromMyVertex",&dxyFromMyVertex_,"dxyFromMyVertex[nTracks]/D");
    rootTree_->Branch("dzFromMyVertex",&dzFromMyVertex_,"dzFromMyVertex[nTracks]/D");
    rootTree_->Branch("dxyErrorFromMyVertex",&dxyErrorFromMyVertex_,"dxyErrorFromMyVertex_[nTracks]/D"); 
    rootTree_->Branch("dzErrorFromMyVertex",&dzErrorFromMyVertex_,"dzErrorFromMyVertex_[nTracks]/D");  
    rootTree_->Branch("IPTsigFromMyVertex",&IPTsigFromMyVertex_,"IPTsigFromMyVertex_[nTracks]/D");    
    rootTree_->Branch("IPLsigFromMyVertex",&IPLsigFromMyVertex_,"IPLsigFromMyVertex_[nTracks]/D"); 
    rootTree_->Branch("hasRecVertex",&hasRecVertex_,"hasRecVertex[nTracks]/I");
    rootTree_->Branch("isGoodTrack",&isGoodTrack_,"isGoodTrack[nTracks]/I");   
  }

  // event histograms
  TFileDirectory EventFeatures = fs->mkdir("EventFeatures");

  TH1F::SetDefaultSumw2(kTRUE);

  h_lumiFromConfig    = EventFeatures.make<TH1F>("h_lumiFromConfig","luminosity from config;;luminosity of present run",1,-0.5,0.5);
  h_lumiFromConfig->SetBinContent(1,intLumi_);

  h_runFromConfig     = EventFeatures.make<TH1I>("h_runFromConfig","run number from config;;run number (from configuration)",
                         runControlNumbers_.size(),0.,runControlNumbers_.size());
  for(unsigned int r=0;r<runControlNumbers_.size();r++){
    h_runFromConfig->SetBinContent(r+1,runControlNumbers_[r]);
  }
  
  h_runFromEvent      = EventFeatures.make<TH1I>("h_runFromEvent","run number from config;;run number (from event)",1,-0.5,0.5);

  h_nTracks           = EventFeatures.make<TH1F>("h_nTracks","number of tracks per event;n_{tracks}/event;n_{events}",300,-0.5,299.5);       
  h_nClus             = EventFeatures.make<TH1F>("h_nClus","number of track clusters;n_{clusters}/event;n_{events}",50,-0.5,49.5);       
  h_nOfflineVertices  = EventFeatures.make<TH1F>("h_nOfflineVertices","number of offline reconstructed vertices;n_{vertices}/event;n_{events}",50,-0.5,49.5);  
  h_runNumber         = EventFeatures.make<TH1F>("h_runNumber","run number;run number;n_{events}",100000,150000.,250000.);       
  h_xOfflineVertex    = EventFeatures.make<TH1F>("h_xOfflineVertex","x-coordinate of offline vertex;x_{vertex};n_{events}",100,-0.1,0.1);    
  h_yOfflineVertex    = EventFeatures.make<TH1F>("h_yOfflineVertex","y-coordinate of offline vertex;y_{vertex};n_{events}",100,-0.1,0.1);    
  h_zOfflineVertex    = EventFeatures.make<TH1F>("h_zOfflineVertex","z-coordinate of offline vertex;z_{vertex};n_{events}",100,-30.,30.);    
  h_xErrOfflineVertex = EventFeatures.make<TH1F>("h_xErrOfflineVertex","x-coordinate error of offline vertex;err_{x}^{vtx};n_{events}",100,0.,0.01); 
  h_yErrOfflineVertex = EventFeatures.make<TH1F>("h_yErrOfflineVertex","y-coordinate error of offline vertex;err_{y}^{vtx};n_{events}",100,0.,0.01); 
  h_zErrOfflineVertex = EventFeatures.make<TH1F>("h_zErrOfflineVertex","z-coordinate error of offline vertex;err_{z}^{vtx};n_{events}",100,0.,10.); 
  h_BSx0              = EventFeatures.make<TH1F>("h_BSx0","x-coordinate of reco beamspot;x^{BS}_{0};n_{events}",100,-0.1,0.1);           
  h_BSy0              = EventFeatures.make<TH1F>("h_BSy0","y-coordinate of reco beamspot;y^{BS}_{0};n_{events}",100,-0.1,0.1);           
  h_BSz0              = EventFeatures.make<TH1F>("h_BSz0","z-coordinate of reco beamspot;z^{BS}_{0};n_{events}",100,-1.,1.);             
  h_Beamsigmaz        = EventFeatures.make<TH1F>("h_Beamsigmaz","z-coordinate beam width;#sigma_{Z}^{beam};n_{events}",100,0.,1.);               
  h_BeamWidthX        = EventFeatures.make<TH1F>("h_BeamWidthX","x-coordinate beam width;#sigma_{X}^{beam};n_{events}",100,0.,0.01);         
  h_BeamWidthY        = EventFeatures.make<TH1F>("h_BeamWidthY","y-coordinate beam width;#sigma_{Y}^{beam};n_{events}",100,0.,0.01);        

  h_etaMax            = EventFeatures.make<TH1F>("etaMax","etaMax",1,-0.5,0.5);
  h_nbins             = EventFeatures.make<TH1F>("nbins","nbins",1,-0.5,0.5);

  // probe track histograms
  TFileDirectory ProbeFeatures = fs->mkdir("ProbeTrackFeatures");

  h_probePt_         = ProbeFeatures.make<TH1F>("h_probePt","p_{T} of probe track;track p_{T} (GeV); tracks",100,0.,50.);   
  h_probeP_          = ProbeFeatures.make<TH1F>("h_probeP","momentum of probe track;track p (GeV); tracks",100,0.,100.);   
  h_probeEta_        = ProbeFeatures.make<TH1F>("h_probeEta","#eta of the probe track;track #eta;tracks",54,-2.8,2.8);  
  h_probePhi_        = ProbeFeatures.make<TH1F>("h_probePhi","#phi of probe track;track #phi (rad);tracks",100,-3.15,3.15);  

  h2_probeEtaPhi_    = ProbeFeatures.make<TH2F>("h2_probeEtaPhi","probe track #phi vs #eta;#eta of probe track;track #phi of probe track (rad); tracks",54,-2.8,2.8,100,-3.15,3.15);  
  h2_probeEtaPt_     = ProbeFeatures.make<TH2F>("h2_probeEtaPt","probe track p_{T} vs #eta;#eta of probe track;track p_{T} (GeV); tracks",54,-2.8,2.8,100,0.,50.);    

  h_probeChi2_       = ProbeFeatures.make<TH1F>("h_probeChi2","#chi^{2} of probe track;track #chi^{2}; tracks",100,0.,100.); 
  h_probeNormChi2_   = ProbeFeatures.make<TH1F>("h_probeNormChi2"," normalized #chi^{2} of probe track;track #chi^{2}/ndof; tracks",100,0.,10.);
  h_probeCharge_     = ProbeFeatures.make<TH1F>("h_probeCharge","charge of profe track;track charge Q;tracks",3,-1.5,1.5);
  h_probeQoverP_     = ProbeFeatures.make<TH1F>("h_probeQoverP","q/p of probe track; track Q/p (GeV^{-1});tracks",200,-1.,1.);
  h_probedzRecoV_    = ProbeFeatures.make<TH1F>("h_probedzRecoV","d_{z}(V_{offline}) of probe track;track d_{z}(V_{off}) (cm);tracks",200,-1.,1.);  
  h_probedxyRecoV_   = ProbeFeatures.make<TH1F>("h_probedxyRecoV","d_{xy}(V_{offline}) of probe track;track d_{xy}(V_{off}) (cm);tracks",200,-1.,1.);      
  h_probedzRefitV_   = ProbeFeatures.make<TH1F>("h_probedzRefitV","d_{z}(V_{refit}) of probe track;track d_{z}(V_{fit}) (cm);tracks",200,-0.5,0.5);
  h_probesignIP2DRefitV_ = ProbeFeatures.make<TH1F>("h_probesignIPRefitV","ip_{2D}(V_{refit}) of probe track;track ip_{2D}(V_{fit}) (cm);tracks",200,-1.,1.);
  h_probedxyRefitV_  = ProbeFeatures.make<TH1F>("h_probedxyRefitV","d_{xy}(V_{refit}) of probe track;track d_{xy}(V_{fit}) (cm);tracks",200,-0.5,0.5); 

  h_probez0RefitV_   = ProbeFeatures.make<TH1F>("h_probez0RefitV","z_{0}(V_{refit}) of probe track;track z_{0}(V_{fit}) (cm);tracks",200,-1.,1.);
  h_probed0RefitV_   = ProbeFeatures.make<TH1F>("h_probed0RefitV","d_{0}(V_{refit}) of probe track;track d_{0}(V_{fit}) (cm);tracks",200,-1.,1.);

  h_probed3DRefitV_  = ProbeFeatures.make<TH1F>("h_probed3DRefitV","d_{3D}(V_{refit}) of probe track;track d_{3D}(V_{fit}) (cm);tracks",200,0.,1.); 
  h_probereszRefitV_ = ProbeFeatures.make<TH1F>("h_probeReszRefitV","z_{track} -z_{V_{refit}};track res_{z}(V_{refit}) (cm);tracks",200,-1.,1.); 

  h_probeRecoVSigZ_  = ProbeFeatures.make<TH1F>("h_probeRecoVSigZ"  ,"Longitudinal DCA Significance (reco);d_{z}(V_{off})/#sigma_{dz};tracks",100,-8,8);
  h_probeRecoVSigXY_ = ProbeFeatures.make<TH1F>("h_probeRecoVSigXY" ,"Transverse DCA Significance (reco);d_{xy}(V_{off})/#sigma_{dxy};tracks",100,-8,8);
  h_probeRefitVSigZ_ = ProbeFeatures.make<TH1F>("h_probeRefitVSigZ" ,"Longitudinal DCA Significance (refit);d_{z}(V_{fit})/#sigma_{dz};tracks",100,-8,8);
  h_probeRefitVSigXY_= ProbeFeatures.make<TH1F>("h_probeRefitVSigXY","Transverse DCA Significance (refit);d_{xy}(V_{fit})/#sigma_{dxy};tracks",100,-8,8);
  h_probeRefitVSig3D_= ProbeFeatures.make<TH1F>("h_probeRefitVSig3D","3D DCA Significance (refit);d_{3D}/#sigma_{3D};tracks",100,0.,20.); 
  h_probeRefitVLogSig3D_ = ProbeFeatures.make<TH1F>("h_probeRefitVLogSig3D","log_{10}(3D DCA-Significance) (refit);log_{10}(d_{3D}/#sigma_{3D});tracks",100,-5.,4.); 
  h_probeRefitVSigResZ_ = ProbeFeatures.make<TH1F>("h_probeRefitVSigResZ" ,"Longitudinal residual significance (refit);(z_{track} -z_{V_{fit}})/#sigma_{res_{z}};tracks",100,-8,8);

  h_probeHits_       = ProbeFeatures.make<TH1F>("h_probeNRechits"    ,"N_{hits}     ;N_{hits}    ;tracks",40,-0.5,39.5);
  h_probeHits1D_     = ProbeFeatures.make<TH1F>("h_probeNRechits1D"  ,"N_{hits} 1D  ;N_{hits} 1D ;tracks",40,-0.5,39.5);
  h_probeHits2D_     = ProbeFeatures.make<TH1F>("h_probeNRechits2D"  ,"N_{hits} 2D  ;N_{hits} 2D ;tracks",40,-0.5,39.5);
  h_probeHitsInTIB_  = ProbeFeatures.make<TH1F>("h_probeNRechitsTIB" ,"N_{hits} TIB ;N_{hits} TIB;tracks",40,-0.5,39.5);
  h_probeHitsInTOB_  = ProbeFeatures.make<TH1F>("h_probeNRechitsTOB" ,"N_{hits} TOB ;N_{hits} TOB;tracks",40,-0.5,39.5);
  h_probeHitsInTID_  = ProbeFeatures.make<TH1F>("h_probeNRechitsTID" ,"N_{hits} TID ;N_{hits} TID;tracks",40,-0.5,39.5);
  h_probeHitsInTEC_  = ProbeFeatures.make<TH1F>("h_probeNRechitsTEC" ,"N_{hits} TEC ;N_{hits} TEC;tracks",40,-0.5,39.5);
  h_probeHitsInBPIX_ = ProbeFeatures.make<TH1F>("h_probeNRechitsBPIX","N_{hits} BPIX;N_{hits} BPIX;tracks",40,-0.5,39.5);
  h_probeHitsInFPIX_ = ProbeFeatures.make<TH1F>("h_probeNRechitsFPIX","N_{hits} FPIX;N_{hits} FPIX;tracks",40,-0.5,39.5);

  // refit vertex features
  TFileDirectory RefitVertexFeatures = fs->mkdir("RefitVertexFeatures");
  h_fitVtxNtracks_          = RefitVertexFeatures.make<TH1F>("h_fitVtxNtracks"  ,"N_{trks} used in vertex fit;N^{fit}_{tracks};vertices"        ,100,-0.5,99.5);
  h_fitVtxNdof_             = RefitVertexFeatures.make<TH1F>("h_fitVtxNdof"     ,"N_{DOF} of vertex fit;N_{DOF} of refit vertex;vertices"          ,100,-0.5,99.5);
  h_fitVtxChi2_             = RefitVertexFeatures.make<TH1F>("h_fitVtxChi2"     ,"#chi^{2} of vertex fit;vertex #chi^{2};vertices"            ,100,-0.5,99.5);
  h_fitVtxChi2ndf_          = RefitVertexFeatures.make<TH1F>("h_fitVtxChi2ndf"  ,"#chi^{2}/ndf of vertex fit;vertex #chi^{2}/ndf;vertices"    ,100,-0.5,9.5);
  h_fitVtxChi2Prob_         = RefitVertexFeatures.make<TH1F>("h_fitVtxChi2Prob" ,"Prob(#chi^{2},ndf) of vertex fit;Prob(#chi^{2},ndf);vertices",40,0.,1.);
  h_fitVtxTrackWeights_     = RefitVertexFeatures.make<TH1F>("h_fitVtxTrackWeights","track weights associated to track;track weights;tracks",40,0.,1.);
  h_fitVtxTrackAverageWeight_ = RefitVertexFeatures.make<TH1F>("h_fitVtxTrackAverageWeight_","average track weight per vertex;#LT track weight #GT;vertices",40,0.,1.);

  if(useTracksFromRecoVtx_) {

    TFileDirectory RecoVertexFeatures = fs->mkdir("RecoVertexFeatures");
    h_recoVtxNtracks_          = RecoVertexFeatures.make<TH1F>("h_recoVtxNtracks"  ,"N^{vtx}_{trks};N^{vtx}_{trks};vertices"        ,100,-0.5,99.5);
    h_recoVtxChi2ndf_          = RecoVertexFeatures.make<TH1F>("h_recoVtxChi2ndf"  ,"#chi^{2}/ndf vtx;#chi^{2}/ndf vtx;vertices"    ,10,-0.5,9.5);
    h_recoVtxChi2Prob_         = RecoVertexFeatures.make<TH1F>("h_recoVtxChi2Prob" ,"Prob(#chi^{2},ndf);Prob(#chi^{2},ndf);vertices",40,0.,1.);
    h_recoVtxSumPt_            = RecoVertexFeatures.make<TH1F>("h_recoVtxSumPt"    ,"Sum(p^{trks}_{T});Sum(p^{trks}_{T});vertices"  ,100,0.,200.);
    
  }


  TFileDirectory DA = fs->mkdir("DA");
  //DA.cd();
  hDA=bookVertexHistograms(DA);
  //for(std::map<std::string,TH1*>::const_iterator hist=hDA.begin(); hist!=hDA.end(); hist++){
  //hist->second->SetDirectory(DA);
  // DA.make<TH1F>(hist->second);
  // }

  // initialize the residuals histograms 

  float dxymax_phi = 2000; 
  float dzmax_phi  = 2000; 
  float dxymax_eta = 3000; 
  float dzmax_eta  = 3000;

  float d3Dmax_phi = hypot(dxymax_phi,dzmax_phi);
  float d3Dmax_eta = hypot(dxymax_eta,dzmax_eta);

  const int mybins_ = 500;

  //  
  // Unbiased track-to-vertex residuals
  // The vertex is refit without the probe track
  //

  TFileDirectory AbsTransPhiRes  = fs->mkdir("Abs_Transv_Phi_Residuals");
  TFileDirectory AbsTransEtaRes  = fs->mkdir("Abs_Transv_Eta_Residuals");
                              
  TFileDirectory AbsLongPhiRes   = fs->mkdir("Abs_Long_Phi_Residuals");
  TFileDirectory AbsLongEtaRes   = fs->mkdir("Abs_Long_Eta_Residuals");
                  
  TFileDirectory Abs3DPhiRes     = fs->mkdir("Abs_3D_Phi_Residuals");
  TFileDirectory Abs3DEtaRes     = fs->mkdir("Abs_3D_Eta_Residuals");

  TFileDirectory NormTransPhiRes = fs->mkdir("Norm_Transv_Phi_Residuals");
  TFileDirectory NormTransEtaRes = fs->mkdir("Norm_Transv_Eta_Residuals");
                              
  TFileDirectory NormLongPhiRes  = fs->mkdir("Norm_Long_Phi_Residuals");
  TFileDirectory NormLongEtaRes  = fs->mkdir("Norm_Long_Eta_Residuals");

  TFileDirectory Norm3DPhiRes    = fs->mkdir("Norm_3D_Phi_Residuals");
  TFileDirectory Norm3DEtaRes    = fs->mkdir("Norm_3D_Eta_Residuals");

  TFileDirectory AbsDoubleDiffRes   = fs->mkdir("Abs_DoubleDiffResiduals");
  TFileDirectory NormDoubleDiffRes  = fs->mkdir("Norm_DoubleDiffResiduals");
  
  // book residuals vs pT histograms
  
  TFileDirectory AbsTranspTRes  = fs->mkdir("Abs_Transv_pT_Residuals"); 
  h_dxy_pT_      = bookResidualsHistogram(AbsTranspTRes,nPtBins_,"dxy","pT");          

  TFileDirectory AbsLongpTRes   = fs->mkdir("Abs_Long_pT_Residuals"); 
  h_dz_pT_       = bookResidualsHistogram(AbsLongpTRes,nPtBins_,"dz","pT");            

  TFileDirectory NormTranspTRes = fs->mkdir("Norm_Transv_pT_Residuals"); 
  h_norm_dxy_pT_ = bookResidualsHistogram(NormTranspTRes,nPtBins_,"norm_dxy","pT");        

  TFileDirectory NormLongpTRes  = fs->mkdir("Norm_Long_pT_Residuals"); 
  h_norm_dz_pT_  = bookResidualsHistogram(NormLongpTRes,nPtBins_,"norm_dz","pT");          

  // book residuals vs pT histograms in central region (|eta|<1.0)
               
  TFileDirectory AbsTranspTCentralRes  = fs->mkdir("Abs_Transv_pTCentral_Residuals"); 
  h_dxy_Central_pT_ = bookResidualsHistogram(AbsTranspTCentralRes,nPtBins_,"dxy","pTCentral");       

  TFileDirectory AbsLongpTCentralRes   = fs->mkdir("Abs_Long_pTCentral_Residuals"); 
  h_dz_Central_pT_  = bookResidualsHistogram(AbsLongpTCentralRes,nPtBins_,"dz","pTCentral");           

  TFileDirectory NormTranspTCentralRes = fs->mkdir("Norm_Transv_pTCentral_Residuals"); 
  h_norm_dxy_Central_pT_ = bookResidualsHistogram(NormTranspTCentralRes,nPtBins_,"norm_dxy","pTCentral");  

  TFileDirectory NormLongpTCentralRes  = fs->mkdir("Norm_Long_pTCentral_Residuals"); 
  h_norm_dz_Central_pT_  = bookResidualsHistogram(NormLongpTCentralRes,nPtBins_,"norm_dz","pTCentral");   

  // book residuals as function of phi and eta

  for ( int i=0; i<nBins_; ++i ) {

    float phiF = (-TMath::Pi()+i*phiSect_)*(180/TMath::Pi());
    float phiL = (-TMath::Pi()+(i+1)*phiSect_)*(180/TMath::Pi());
    
    float etaF=-etaOfProbe_+i*etaSect_;
    float etaL=-etaOfProbe_+(i+1)*etaSect_;
    
    // dxy vs phi and eta
     
    a_dxyPhiResiduals[i] = AbsTransPhiRes.make<TH1F>(Form("histo_dxy_phi_plot%i",i),
                             Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{xy} [#mum];tracks",phiF,phiL),
                             mybins_,-dxymax_phi,dxymax_phi);
    
    a_dxyEtaResiduals[i] = AbsTransEtaRes.make<TH1F>(Form("histo_dxy_eta_plot%i",i),
                             Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{xy} [#mum];tracks",etaF,etaL),
                             mybins_,-dxymax_eta,dxymax_eta);


    // dx vs phi and eta
     
    a_dxPhiResiduals[i] = AbsTransPhiRes.make<TH1F>(Form("histo_dx_phi_plot%i",i),
                            Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{x} [#mum];tracks",phiF,phiL),
                            mybins_,-dxymax_phi,dxymax_phi);
    
    a_dxEtaResiduals[i] = AbsTransEtaRes.make<TH1F>(Form("histo_dx_eta_plot%i",i),
                            Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{x} [#mum];tracks",etaF,etaL),
                            mybins_,-dxymax_eta,dxymax_eta);


    // dy vs phi and eta
     
    a_dyPhiResiduals[i] = AbsTransPhiRes.make<TH1F>(Form("histo_dy_phi_plot%i",i),
                            Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{y} [#mum];tracks",phiF,phiL),
                            mybins_,-dxymax_phi,dxymax_phi);
    
    a_dyEtaResiduals[i] = AbsTransEtaRes.make<TH1F>(Form("histo_dy_eta_plot%i",i),
                            Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{y} [#mum];tracks",etaF,etaL),
                            mybins_,-dxymax_eta,dxymax_eta);
    
    // IP2D vs phi and eta

    a_IP2DPhiResiduals[i] = AbsTransPhiRes.make<TH1F>(Form("histo_IP2D_phi_plot%i",i),
                             Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;IP_{2D} [#mum];tracks",phiF,phiL),
                             mybins_,-dxymax_phi,dxymax_phi);
    
    a_IP2DEtaResiduals[i] = AbsTransEtaRes.make<TH1F>(Form("histo_IP2D_eta_plot%i",i),
                             Form("%.2f<#eta^{probe}_{tk}<%.2f;IP_{2D} [#mum];tracks",etaF,etaL),
                             mybins_,-dxymax_eta,dxymax_eta);

    // IP3D vs phi and eta

    a_IP3DPhiResiduals[i] = Abs3DPhiRes.make<TH1F>(Form("histo_IP3D_phi_plot%i",i),
                           Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;IP_{3D} [#mum];tracks",phiF,phiL),
                           mybins_,-dxymax_phi,dxymax_phi);
    
    a_IP3DEtaResiduals[i] = Abs3DEtaRes.make<TH1F>(Form("histo_IP3D_eta_plot%i",i),
                           Form("%.2f<#eta^{probe}_{tk}<%.2f;IP_{3D} [#mum];tracks",etaF,etaL),
                           mybins_,-dxymax_eta,dxymax_eta);

    // dz vs phi and eta

    a_dzPhiResiduals[i]  = AbsLongPhiRes.make<TH1F>(Form("histo_dz_phi_plot%i",i),
                            Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{z} [#mum];tracks",phiF,phiL),
                            mybins_,-dzmax_phi,dzmax_phi);
    
    a_dzEtaResiduals[i]  = AbsLongEtaRes.make<TH1F>(Form("histo_dz_eta_plot%i",i),
                            Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{z} [#mum];tracks",etaF,etaL),
                            mybins_,-dzmax_eta,dzmax_eta);


    // resz vs phi and eta

    a_reszPhiResiduals[i]  = AbsLongPhiRes.make<TH1F>(Form("histo_resz_phi_plot%i",i),
                            Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;z_{trk} - z_{vtx} [#mum];tracks",phiF,phiL),
                            mybins_,-dzmax_phi,dzmax_phi);
    
    a_reszEtaResiduals[i]  = AbsLongEtaRes.make<TH1F>(Form("histo_resz_eta_plot%i",i),
                            Form("%.2f<#eta^{probe}_{tk}<%.2f;z_{trk} - z_{vtx} [#mum];tracks",etaF,etaL),
                            mybins_,-dzmax_eta,dzmax_eta);

    // d3D vs phi and eta

    a_d3DPhiResiduals[i] = Abs3DPhiRes.make<TH1F>(Form("histo_d3D_phi_plot%i",i),
                          Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{3D} [#mum];tracks",phiF,phiL),
                          mybins_,0.,d3Dmax_phi);
    
    a_d3DEtaResiduals[i] = Abs3DEtaRes.make<TH1F>(Form("histo_d3D_eta_plot%i",i),
                          Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{3D} [#mum];tracks",etaF,etaL),
                          mybins_,0.,d3Dmax_eta);
       
    //  _  _                    _ _           _   ___        _    _           _    
    // | \| |___ _ _ _ __  __ _| (_)______ __| | | _ \___ __(_)__| |_  _ __ _| |___
    // | .` / _ \ '_| '  \/ _` | | |_ / -_) _` | |   / -_|_-< / _` | || / _` | (_-<
    // |_|\_\___/_| |_|_|_\__,_|_|_/__\___\__,_| |_|_\___/__/_\__,_|\_,_\__,_|_/__/
    //
                                                                             
    // normalized dxy vs eta and phi
                
    n_dxyPhiResiduals[i] = NormTransPhiRes.make<TH1F>(Form("histo_norm_dxy_phi_plot%i",i),
                              Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{xy}/#sigma_{d_{xy}};tracks",phiF,phiL),
                              mybins_,-dxymax_phi/100.,dxymax_phi/100.);
    
    n_dxyEtaResiduals[i] = NormTransEtaRes.make<TH1F>(Form("histo_norm_dxy_eta_plot%i",i),
                              Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{xy}/#sigma_{d_{xy}};tracks",etaF,etaL),
                              mybins_,-dxymax_eta/100.,dxymax_eta/100.);
    
    // normalized IP2d vs eta and phi
    
    n_IP2DPhiResiduals[i] = NormTransPhiRes.make<TH1F>(Form("histo_norm_IP2D_phi_plot%i",i),
                               Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;IP_{2D}/#sigma_{IP_{2D}};tracks",phiF,phiL),
                               mybins_,-dxymax_phi/100.,dxymax_phi/100.);
    
    n_IP2DEtaResiduals[i] = NormTransEtaRes.make<TH1F>(Form("histo_norm_IP2D_eta_plot%i",i),
                               Form("%.2f<#eta^{probe}_{tk}<%.2f;IP_{2D}/#sigma_{IP_{2D}};tracks",etaF,etaL),
                               mybins_,-dxymax_eta/100.,dxymax_eta/100.);
    
    // normalized IP3d vs eta and phi
    
    n_IP3DPhiResiduals[i] = Norm3DPhiRes.make<TH1F>(Form("histo_norm_IP3D_phi_plot%i",i),
                            Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;IP_{3D}/#sigma_{IP_{3D}};tracks",phiF,phiL),
                            mybins_,-dxymax_phi/100.,dxymax_phi/100.);
    
    n_IP3DEtaResiduals[i] = Norm3DEtaRes.make<TH1F>(Form("histo_norm_IP3D_eta_plot%i",i),
                            Form("%.2f<#eta^{probe}_{tk}<%.2f;IP_{3D}/#sigma_{IP_{3D}};tracks",etaF,etaL),
                            mybins_,-dxymax_eta/100.,dxymax_eta/100.);

    // normalized dz vs phi and eta

    n_dzPhiResiduals[i]  = NormLongPhiRes.make<TH1F>(Form("histo_norm_dz_phi_plot%i",i),
                             Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{z}/#sigma_{d_{z}};tracks",phiF,phiL),
                             mybins_,-dzmax_phi/100.,dzmax_phi/100.);
    
    n_dzEtaResiduals[i]  = NormLongEtaRes.make<TH1F>(Form("histo_norm_dz_eta_plot%i",i),
                             Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{z}/#sigma_{d_{z}};tracks",etaF,etaL),
                             mybins_,-dzmax_eta/100.,dzmax_eta/100.);

    // pull of resz

    n_reszPhiResiduals[i]  = NormLongPhiRes.make<TH1F>(Form("histo_norm_resz_phi_plot%i",i),
                             Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;(z_{trk}-z_{vtx})/#sigma_{res_{z}};tracks",phiF,phiL),
                             mybins_,-dzmax_phi/100.,dzmax_phi/100.);
    
    n_reszEtaResiduals[i]  = NormLongEtaRes.make<TH1F>(Form("histo_norm_resz_eta_plot%i",i),
                             Form("%.2f<#eta^{probe}_{tk}<%.2f;(z_{trk}-z_{vtx})/#sigma_{res_{z}};tracks",etaF,etaL),
                             mybins_,-dzmax_eta/100.,dzmax_eta/100.);

    // normalized d3D vs phi and eta

    n_d3DPhiResiduals[i] = Norm3DPhiRes.make<TH1F>(Form("histo_norm_d3D_phi_plot%i",i),
                           Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{3D}/#sigma_{d_{3D}};tracks",phiF,phiL),
                           mybins_,0.,d3Dmax_phi/100.);
    
    n_d3DEtaResiduals[i] = Norm3DEtaRes.make<TH1F>(Form("histo_norm_d3D_eta_plot%i",i),
                           Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{3D}/#sigma_{d_{3D}};tracks",etaF,etaL),
                           mybins_,0.,d3Dmax_eta/100.);

    //  ___           _    _     ___  _  __  __   ___        _    _           _    
    // |   \ ___ _  _| |__| |___|   \(_)/ _|/ _| | _ \___ __(_)__| |_  _ __ _| |___
    // | |) / _ \ || | '_ \ / -_) |) | |  _|  _| |   / -_|_-< / _` | || / _` | (_-<
    // |___/\___/\_,_|_.__/_\___|___/|_|_| |_|   |_|_\___/__/_\__,_|\_,_\__,_|_/__/
    
    for ( int j=0; j<nBins_; ++j ) {
 
      a_dxyResidualsMap[i][j] = AbsDoubleDiffRes.make<TH1F>(Form("histo_dxy_eta_plot%i_phi_plot%i",i,j),
                                Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{xy};tracks",etaF,etaL,phiF,phiL),
                                mybins_,-dzmax_eta,dzmax_eta);
      
      a_dzResidualsMap[i][j]  = AbsDoubleDiffRes.make<TH1F>(Form("histo_dz_eta_plot%i_phi_plot%i",i,j),
                                Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{z};tracks",etaF,etaL,phiF,phiL),
                                mybins_,-dzmax_eta,dzmax_eta);
      
      a_d3DResidualsMap[i][j] = AbsDoubleDiffRes.make<TH1F>(Form("histo_d3D_eta_plot%i_phi_plot%i",i,j),
                                Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{3D};tracks",etaF,etaL,phiF,phiL),
                                mybins_,0.,d3Dmax_eta);
      
      n_dxyResidualsMap[i][j] = NormDoubleDiffRes.make<TH1F>(Form("histo_norm_dxy_eta_plot%i_phi_plot%i",i,j),
                                 Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{xy}/#sigma_{d_{xy}};tracks",etaF,etaL,phiF,phiL),
                                 mybins_,-dzmax_eta/100,dzmax_eta/100);

      n_dzResidualsMap[i][j]  = NormDoubleDiffRes.make<TH1F>(Form("histo_norm_dz_eta_plot%i_phi_plot%i",i,j),
                                 Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{z}/#sigma_{d_{z}};tracks",etaF,etaL,phiF,phiL),
                                 mybins_,-dzmax_eta/100,dzmax_eta/100);

      n_d3DResidualsMap[i][j] = NormDoubleDiffRes.make<TH1F>(Form("histo_norm_d3D_eta_plot%i_phi_plot%i",i,j),
                                 Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{3D}/#sigma_{d_{3D}};tracks",etaF,etaL,phiF,phiL),
                                 mybins_,0.,d3Dmax_eta);

    }
  }

  // declaration of the directories
  
  TFileDirectory BiasVsParameter = fs->mkdir("BiasVsParameter");

  a_dxyVsPhi = BiasVsParameter.make<TH2F>("h2_dxy_vs_phi","d_{xy} vs track #phi;track #phi [rad];track d_{xy}(PV) [#mum]",
                      48,-TMath::Pi(),TMath::Pi(),mybins_,-dxymax_phi,dxymax_phi); 
 
  a_dzVsPhi  = BiasVsParameter.make<TH2F>("h2_dz_vs_phi","d_{z} vs track #phi;track #phi [rad];track d_{z}(PV) [#mum]",
                      48,-TMath::Pi(),TMath::Pi(),mybins_,-dzmax_phi,dzmax_phi);   
               
  n_dxyVsPhi = BiasVsParameter.make<TH2F>("h2_n_dxy_vs_phi","d_{xy}/#sigma_{d_{xy}} vs track #phi;track #phi [rad];track d_{xy}(PV)/#sigma_{d_{xy}}",
                      48,-TMath::Pi(),TMath::Pi(),mybins_,-dxymax_phi/100.,dxymax_phi/100.); 
  
  n_dzVsPhi  = BiasVsParameter.make<TH2F>("h2_n_dz_vs_phi","d_{z}/#sigma_{d_{z}} vs track #phi;track #phi [rad];track d_{z}(PV)/#sigma_{d_{z}}",
                      48,-TMath::Pi(),TMath::Pi(),mybins_,-dzmax_phi/100.,dzmax_phi/100.);   
               
  a_dxyVsEta = BiasVsParameter.make<TH2F>("h2_dxy_vs_eta","d_{xy} vs track #eta;track #eta;track d_{xy}(PV) [#mum]",
                      48,-etaOfProbe_,etaOfProbe_,mybins_,-dxymax_eta,dzmax_eta);
  
  a_dzVsEta  = BiasVsParameter.make<TH2F>("h2_dz_vs_eta","d_{z} vs track #eta;track #eta;track d_{z}(PV) [#mum]",
                      48,-etaOfProbe_,etaOfProbe_,mybins_,-dzmax_eta,dzmax_eta);   
               
  n_dxyVsEta = BiasVsParameter.make<TH2F>("h2_n_dxy_vs_eta","d_{xy}/#sigma_{d_{xy}} vs track #eta;track #eta;track d_{xy}(PV)/#sigma_{d_{xy}}",
                      48,-etaOfProbe_,etaOfProbe_,mybins_,-dxymax_eta/100.,dxymax_eta/100.);  

  n_dzVsEta  = BiasVsParameter.make<TH2F>("h2_n_dz_vs_eta","d_{z}/#sigma_{d_{z}} vs track #eta;track #eta;track d_{z}(PV)/#sigma_{d_{z}}",
                      48,-etaOfProbe_,etaOfProbe_,mybins_,-dzmax_eta/100.,dzmax_eta/100.);   

  TFileDirectory MeanTrendsDir   = fs->mkdir("MeanTrends");
  TFileDirectory WidthTrendsDir  = fs->mkdir("WidthTrends");
  TFileDirectory MedianTrendsDir = fs->mkdir("MedianTrends");
  TFileDirectory MADTrendsDir    = fs->mkdir("MADTrends");

  TFileDirectory Mean2DMapsDir   = fs->mkdir("MeanMaps");
  TFileDirectory Width2DMapsDir  = fs->mkdir("WidthMaps");

  Double_t highedge=nBins_-0.5;
  Double_t lowedge=-0.5;

  // means and widths from the fit

  a_dxyPhiMeanTrend  = MeanTrendsDir.make<TH1F> ("means_dxy_phi",
                         "#LT d_{xy} #GT vs #phi sector;#varphi (sector) [degrees];#LT d_{xy} #GT [#mum]",
                         nBins_,lowedge,highedge); 
  
  a_dxyPhiWidthTrend = WidthTrendsDir.make<TH1F>("widths_dxy_phi",
                         "#sigma_{d_{xy}} vs #phi sector;#varphi (sector) [degrees];#sigma_{d_{xy}} [#mum]",
                         nBins_,lowedge,highedge);
  
  a_dzPhiMeanTrend   = MeanTrendsDir.make<TH1F> ("means_dz_phi",
                         "#LT d_{z} #GT vs #phi sector;#varphi (sector) [degrees];#LT d_{z} #GT [#mum]",
                         nBins_,lowedge,highedge); 

  a_dzPhiWidthTrend  = WidthTrendsDir.make<TH1F>("widths_dz_phi","#sigma_{d_{z}} vs #phi sector;#varphi (sector) [degrees];#sigma_{d_{z}} [#mum]",
                         nBins_,lowedge,highedge);
              
  a_dxyEtaMeanTrend  = MeanTrendsDir.make<TH1F> ("means_dxy_eta",
                         "#LT d_{xy} #GT vs #eta sector;#eta (sector);#LT d_{xy} #GT [#mum]",
                         nBins_,lowedge,highedge);

  a_dxyEtaWidthTrend = WidthTrendsDir.make<TH1F>("widths_dxy_eta",
                         "#sigma_{d_{xy}} vs #eta sector;#eta (sector);#sigma_{d_{xy}} [#mum]",
                         nBins_,lowedge,highedge);
  
  a_dzEtaMeanTrend   = MeanTrendsDir.make<TH1F> ("means_dz_eta",
                         "#LT d_{z} #GT vs #eta sector;#eta (sector);#LT d_{z} #GT [#mum]"
                         ,nBins_,lowedge,highedge); 
  
  a_dzEtaWidthTrend  = WidthTrendsDir.make<TH1F>("widths_dz_eta",
                         "#sigma_{d_{xy}} vs #eta sector;#eta (sector);#sigma_{d_{z}} [#mum]",
                         nBins_,lowedge,highedge);
              
  n_dxyPhiMeanTrend  = MeanTrendsDir.make<TH1F> ("norm_means_dxy_phi",
                         "#LT d_{xy}/#sigma_{d_{xy}} #GT vs #phi sector;#varphi (sector) [degrees];#LT d_{xy}/#sigma_{d_{xy}} #GT",
                         nBins_,lowedge,highedge);
  
  n_dxyPhiWidthTrend = WidthTrendsDir.make<TH1F>("norm_widths_dxy_phi",
                         "width(d_{xy}/#sigma_{d_{xy}}) vs #phi sector;#varphi (sector) [degrees]; width(d_{xy}/#sigma_{d_{xy}})",
                         nBins_,lowedge,highedge);
  
  n_dzPhiMeanTrend   = MeanTrendsDir.make<TH1F> ("norm_means_dz_phi",
                         "#LT d_{z}/#sigma_{d_{z}} #GT vs #phi sector;#varphi (sector) [degrees];#LT d_{z}/#sigma_{d_{z}} #GT",
                         nBins_,lowedge,highedge); 
  
  n_dzPhiWidthTrend  = WidthTrendsDir.make<TH1F>("norm_widths_dz_phi",
                         "width(d_{z}/#sigma_{d_{z}}) vs #phi sector;#varphi (sector) [degrees];width(d_{z}/#sigma_{d_{z}})",
                         nBins_,lowedge,highedge);
                                            
  n_dxyEtaMeanTrend  = MeanTrendsDir.make<TH1F> ("norm_means_dxy_eta",
                         "#LT d_{xy}/#sigma_{d_{xy}} #GT vs #eta sector;#eta (sector);#LT d_{xy}/#sigma_{d_{z}} #GT",
                         nBins_,lowedge,highedge);
  
  n_dxyEtaWidthTrend = WidthTrendsDir.make<TH1F>("norm_widths_dxy_eta",
                         "width(d_{xy}/#sigma_{d_{xy}}) vs #eta sector;#eta (sector);width(d_{xy}/#sigma_{d_{z}})",
                         nBins_,lowedge,highedge);

  n_dzEtaMeanTrend   = MeanTrendsDir.make<TH1F> ("norm_means_dz_eta",
                         "#LT d_{z}/#sigma_{d_{z}} #GT vs #eta sector;#eta (sector);#LT d_{z}/#sigma_{d_{z}} #GT",
                         nBins_,lowedge,highedge);  
  
  n_dzEtaWidthTrend  = WidthTrendsDir.make<TH1F>("norm_widths_dz_eta",
                         "width(d_{z}/#sigma_{d_{z}}) vs #eta sector;#eta (sector);width(d_{z}/#sigma_{d_{z}})",
                         nBins_,lowedge,highedge);                        
  

  // means and widhts vs pT and pTCentral
  
  a_dxypTMeanTrend  = MeanTrendsDir.make<TH1F> ("means_dxy_pT",
                        "#LT d_{xy} #GT vs pT;p_{T} [GeV];#LT d_{xy} #GT [#mum]",
                        48,mypT_bins_); 
  
  a_dxypTWidthTrend = WidthTrendsDir.make<TH1F>("widths_dxy_pT",
                        "#sigma_{d_{xy}} vs pT;p_{T} [GeV];#sigma_{d_{xy}} [#mum]",
                        48,mypT_bins_);
  
  a_dzpTMeanTrend   = MeanTrendsDir.make<TH1F> ("means_dz_pT",
                        "#LT d_{z} #GT vs pT;p_{T} [GeV];#LT d_{z} #GT [#mum]",
                        48,mypT_bins_); 
  
  a_dzpTWidthTrend  = WidthTrendsDir.make<TH1F>("widths_dz_pT","#sigma_{d_{z}} vs pT;p_{T} [GeV];#sigma_{d_{z}} [#mum]",
                        48,mypT_bins_);
  
  a_dxypTCentralMeanTrend  = MeanTrendsDir.make<TH1F> ("means_dxy_pTCentral",
                               "#LT d_{xy} #GT vs p_{T};p_{T}(|#eta|<1.) [GeV];#LT d_{xy} #GT [#mum]",
                               48,mypT_bins_);
  
  a_dxypTCentralWidthTrend = WidthTrendsDir.make<TH1F>("widths_dxy_pTCentral",
                               "#sigma_{d_{xy}} vs p_{T};p_{T}(|#eta|<1.) [GeV];#sigma_{d_{xy}} [#mum]",
                               48,mypT_bins_);
  
  a_dzpTCentralMeanTrend   = MeanTrendsDir.make<TH1F> ("means_dz_pTCentral",
                               "#LT d_{z} #GT vs p_{T};p_{T}(|#eta|<1.) [GeV];#LT d_{z} #GT [#mum]"
                               ,48,mypT_bins_); 
  
  a_dzpTCentralWidthTrend  = WidthTrendsDir.make<TH1F>("widths_dz_pTCentral",
                               "#sigma_{d_{z}} vs p_{T};p_{T}(|#eta|<1.) [GeV];#sigma_{d_{z}} [#mum]",
                               48,mypT_bins_);
  
  n_dxypTMeanTrend  = MeanTrendsDir.make<TH1F> ("norm_means_dxy_pT",
                        "#LT d_{xy}/#sigma_{d_{xy}} #GT vs pT;p_{T} [GeV];#LT d_{xy}/#sigma_{d_{xy}} #GT",
                        48,mypT_bins_);
  
  n_dxypTWidthTrend = WidthTrendsDir.make<TH1F>("norm_widths_dxy_pT",
                        "width(d_{xy}/#sigma_{d_{xy}}) vs pT;p_{T} [GeV]; width(d_{xy}/#sigma_{d_{xy}})",
                        48,mypT_bins_);
  
  n_dzpTMeanTrend   = MeanTrendsDir.make<TH1F> ("norm_means_dz_pT",
                        "#LT d_{z}/#sigma_{d_{z}} #GT vs pT;p_{T} [GeV];#LT d_{z}/#sigma_{d_{z}} #GT",
                        48,mypT_bins_); 
  
  n_dzpTWidthTrend  = WidthTrendsDir.make<TH1F>("norm_widths_dz_pT",
                        "width(d_{z}/#sigma_{d_{z}}) vs pT;p_{T} [GeV];width(d_{z}/#sigma_{d_{z}})",
                        48,mypT_bins_);
  
  n_dxypTCentralMeanTrend  = MeanTrendsDir.make<TH1F> ("norm_means_dxy_pTCentral",
                               "#LT d_{xy}/#sigma_{d_{xy}} #GT vs p_{T};p_{T}(|#eta|<1.) [GeV];#LT d_{xy}/#sigma_{d_{z}} #GT",
                               48,mypT_bins_);
  
  n_dxypTCentralWidthTrend = WidthTrendsDir.make<TH1F>("norm_widths_dxy_pTCentral",
                               "width(d_{xy}/#sigma_{d_{xy}}) vs p_{T};p_{T}(|#eta|<1.) [GeV];width(d_{xy}/#sigma_{d_{z}})",
                               48,mypT_bins_);
  
  n_dzpTCentralMeanTrend   = MeanTrendsDir.make<TH1F> ("norm_means_dz_pTCentral",
                               "#LT d_{z}/#sigma_{d_{z}} #GT vs p_{T};p_{T}(|#eta|<1.) [GeV];#LT d_{z}/#sigma_{d_{z}} #GT",
                               48,mypT_bins_);  
  
  n_dzpTCentralWidthTrend  = WidthTrendsDir.make<TH1F>("norm_widths_dz_pTCentral",
                               "width(d_{z}/#sigma_{d_{z}}) vs p_{T};p_{T}(|#eta|<1.) [GeV];width(d_{z}/#sigma_{d_{z}})",
                               48,mypT_bins_); 

  // 2D maps

  a_dxyMeanMap       =  Mean2DMapsDir.make<TH2F>  ("means_dxy_map",
                           "#LT d_{xy} #GT map;#eta (sector);#varphi (sector) [degrees]",
                           nBins_,lowedge,highedge,nBins_,lowedge,highedge);
  
  a_dzMeanMap        =  Mean2DMapsDir.make<TH2F>  ("means_dz_map",
                           "#LT d_{z} #GT map;#eta (sector);#varphi (sector) [degrees]",
                           nBins_,lowedge,highedge,nBins_,lowedge,highedge);
             
  n_dxyMeanMap       =  Mean2DMapsDir.make<TH2F>  ("norm_means_dxy_map",
                           "#LT d_{xy}/#sigma_{d_{xy}} #GT map;#eta (sector);#varphi (sector) [degrees]",
                           nBins_,lowedge,highedge,nBins_,lowedge,highedge);

  n_dzMeanMap        =  Mean2DMapsDir.make<TH2F>  ("norm_means_dz_map",
                           "#LT d_{z}/#sigma_{d_{z}} #GT map;#eta (sector);#varphi (sector) [degrees]",
                           nBins_,lowedge,highedge,nBins_,lowedge,highedge);
             
  a_dxyWidthMap      =  Width2DMapsDir.make<TH2F> ("widths_dxy_map",
                           "#sigma_{d_{xy}} map;#eta (sector);#varphi (sector) [degrees]",
                           nBins_,lowedge,highedge,nBins_,lowedge,highedge);

  a_dzWidthMap       =  Width2DMapsDir.make<TH2F> ("widths_dz_map",
                           "#sigma_{d_{z}} map;#eta (sector);#varphi (sector) [degrees]",
                           nBins_,lowedge,highedge,nBins_,lowedge,highedge);
             
  n_dxyWidthMap      =  Width2DMapsDir.make<TH2F> ("norm_widths_dxy_map",
                           "width(d_{xy}/#sigma_{d_{xy}}) map;#eta (sector);#varphi (sector) [degrees]",
                           nBins_,lowedge,highedge,nBins_,lowedge,highedge);
  
  n_dzWidthMap       =  Width2DMapsDir.make<TH2F> ("norm_widths_dz_map",
                           "width(d_{z}/#sigma_{d_{z}}) map;#eta (sector);#varphi (sector) [degrees]",
                           nBins_,lowedge,highedge,nBins_,lowedge,highedge);

  // medians and MADs

  a_dxyPhiMedianTrend = MedianTrendsDir.make<TH1F>("medians_dxy_phi",
                           "Median of d_{xy} vs #phi sector;#varphi (sector) [degrees];#mu_{1/2}(d_{xy}) [#mum]",
                           nBins_,lowedge,highedge);        
         
  a_dxyPhiMADTrend    = MADTrendsDir.make<TH1F>   ("MADs_dxy_phi",
                           "Median absolute deviation of d_{xy} vs #phi sector;#varphi (sector) [degrees];MAD(d_{xy}) [#mum]",
                           nBins_,lowedge,highedge);
                    
  a_dzPhiMedianTrend  = MedianTrendsDir.make<TH1F>("medians_dz_phi",
                           "Median of d_{z} vs #phi sector;#varphi (sector) [degrees];#mu_{1/2}(d_{z}) [#mum]",
                           nBins_,lowedge,highedge); 
                    
  a_dzPhiMADTrend     = MADTrendsDir.make<TH1F>   ("MADs_dz_phi",
                           "Median absolute deviation of d_{z} vs #phi sector;#varphi (sector) [degrees];MAD(d_{z}) [#mum]",
                           nBins_,lowedge,highedge);                    
  
  a_dxyEtaMedianTrend = MedianTrendsDir.make<TH1F>("medians_dxy_eta",
                           "Median of d_{xy} vs #eta sector;#eta (sector);#mu_{1/2}(d_{xy}) [#mum]",
                           nBins_,lowedge,highedge);        
                  
  a_dxyEtaMADTrend    = MADTrendsDir.make<TH1F>   ("MADs_dxy_eta",
                           "Median absolute deviation of d_{xy} vs #eta sector;#eta (sector);MAD(d_{xy}) [#mum]",
                           nBins_,lowedge,highedge);    
                    
  a_dzEtaMedianTrend  = MedianTrendsDir.make<TH1F>("medians_dz_eta",
                           "Median of d_{z} vs #eta sector;#eta (sector);#mu_{1/2}(d_{z}) [#mum]",
                           nBins_,lowedge,highedge);    
                      
  a_dzEtaMADTrend     = MADTrendsDir.make<TH1F>   ("MADs_dz_eta",
                           "Median absolute deviation of d_{z} vs #eta sector;#eta (sector);MAD(d_{z}) [#mum]",
                           nBins_,lowedge,highedge);                          
  
  n_dxyPhiMedianTrend = MedianTrendsDir.make<TH1F>("norm_medians_dxy_phi",
                           "Median of d_{xy}/#sigma_{d_{xy}} vs #phi sector;#varphi (sector) [degrees];#mu_{1/2}(d_{xy}/#sigma_{d_{xy}})",
                           nBins_,lowedge,highedge); 

  n_dxyPhiMADTrend    = MADTrendsDir.make<TH1F>   ("norm_MADs_dxy_phi",
                           "Median absolute deviation of d_{xy}/#sigma_{d_{xy}} vs #phi sector;#varphi (sector) [degrees]; MAD(d_{xy}/#sigma_{d_{xy}})",
                           nBins_,lowedge,highedge);   

  n_dzPhiMedianTrend  = MedianTrendsDir.make<TH1F>("norm_medians_dz_phi",
                           "Median of d_{z}/#sigma_{d_{z}} vs #phi sector;#varphi (sector) [degrees];#mu_{1/2}(d_{z}/#sigma_{d_{z}})",
                           nBins_,lowedge,highedge);  
    
  n_dzPhiMADTrend     = MADTrendsDir.make<TH1F>   ("norm_MADs_dz_phi",
                           "Median absolute deviation of d_{z}/#sigma_{d_{z}} vs #phi sector;#varphi (sector) [degrees];MAD(d_{z}/#sigma_{d_{z}})",
                           nBins_,lowedge,highedge);        
  
  n_dxyEtaMedianTrend = MedianTrendsDir.make<TH1F>("norm_medians_dxy_eta",
                           "Median of d_{xy}/#sigma_{d_{xy}} vs #eta sector;#eta (sector);#mu_{1/2}(d_{xy}/#sigma_{d_{z}})",
                           nBins_,lowedge,highedge);
            
  n_dxyEtaMADTrend    = MADTrendsDir.make<TH1F>   ("norm_MADs_dxy_eta",
                           "Median absolute deviation of d_{xy}/#sigma_{d_{xy}} vs #eta sector;#eta (sector);MAD(d_{xy}/#sigma_{d_{z}})",
                           nBins_,lowedge,highedge);
            
  n_dzEtaMedianTrend  = MedianTrendsDir.make<TH1F>("norm_medians_dz_eta",
                           "Median of d_{z}/#sigma_{d_{z}} vs #eta sector;#eta (sector);#mu_{1/2}(d_{z}/#sigma_{d_{z}})",
                           nBins_,lowedge,highedge);    
        
  n_dzEtaMADTrend     = MADTrendsDir.make<TH1F>   ("norm_MADs_dz_eta",
                           "Median absolute deviation of d_{z}/#sigma_{d_{z}} vs #eta sector;#eta (sector);MAD(d_{z}/#sigma_{d_{z}})",
                           nBins_,lowedge,highedge);                      


  //  
  // plots of biased residuals
  // The vertex includes the probe track
  //

  if (useTracksFromRecoVtx_){

    TFileDirectory AbsTransPhiBiasRes  = fs->mkdir("Abs_Transv_Phi_BiasResiduals");
    TFileDirectory AbsTransEtaBiasRes  = fs->mkdir("Abs_Transv_Eta_BiasResiduals");
    
    TFileDirectory AbsLongPhiBiasRes   = fs->mkdir("Abs_Long_Phi_BiasResiduals");
    TFileDirectory AbsLongEtaBiasRes   = fs->mkdir("Abs_Long_Eta_BiasResiduals");
    
    TFileDirectory NormTransPhiBiasRes = fs->mkdir("Norm_Transv_Phi_BiasResiduals");
    TFileDirectory NormTransEtaBiasRes = fs->mkdir("Norm_Transv_Eta_BiasResiduals");
    
    TFileDirectory NormLongPhiBiasRes  = fs->mkdir("Norm_Long_Phi_BiasResiduals");
    TFileDirectory NormLongEtaBiasRes  = fs->mkdir("Norm_Long_Eta_BiasResiduals");
    
    TFileDirectory AbsDoubleDiffBiasRes   = fs->mkdir("Abs_DoubleDiffBiasResiduals");
    TFileDirectory NormDoubleDiffBiasRes  = fs->mkdir("Norm_DoubleDiffBiasResiduals");
    
    for ( int i=0; i<nBins_; ++i ) {
      
      float phiF = (-TMath::Pi()+i*phiSect_)*(180/TMath::Pi());
      float phiL = (-TMath::Pi()+(i+1)*phiSect_)*(180/TMath::Pi());
      
      float etaF=-etaOfProbe_+i*etaSect_;
      float etaL=-etaOfProbe_+(i+1)*etaSect_;
      
      a_dxyPhiBiasResiduals[i] = AbsTransPhiBiasRes.make<TH1F>(Form("histo_dxy_phi_plot%i",i),
                                   Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{xy} [#mum];tracks",phiF,phiL),
                                   mybins_,-dxymax_phi,dxymax_phi);

      a_dxyEtaBiasResiduals[i] = AbsTransEtaBiasRes.make<TH1F>(Form("histo_dxy_eta_plot%i",i),
                                   Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{xy} [#mum];tracks",etaF,etaL),
                                   mybins_,-dxymax_eta,dxymax_eta);
      
      a_dzPhiBiasResiduals[i]  = AbsLongPhiBiasRes.make<TH1F>(Form("histo_dz_phi_plot%i",i),
                                  Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f #circ;d_{z} [#mum];tracks",phiF,phiL),
                                  mybins_,-dzmax_phi,dzmax_phi);

      a_dzEtaBiasResiduals[i]  = AbsLongEtaBiasRes.make<TH1F>(Form("histo_dz_eta_plot%i",i),
                                  Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{z} [#mum];tracks",etaF,etaL),
                                  mybins_,-dzmax_eta,dzmax_eta);
      
      n_dxyPhiBiasResiduals[i] = NormTransPhiBiasRes.make<TH1F>(Form("histo_norm_dxy_phi_plot%i",i),
                                Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{xy}/#sigma_{d_{xy}};tracks",phiF,phiL),
                                mybins_,-dxymax_phi/100.,dxymax_phi/100.);

      n_dxyEtaBiasResiduals[i] = NormTransEtaBiasRes.make<TH1F>(Form("histo_norm_dxy_eta_plot%i",i),
                                Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{xy}/#sigma_{d_{xy}};tracks",etaF,etaL),
                                mybins_,-dxymax_eta/100.,dxymax_eta/100.);
      
      n_dzPhiBiasResiduals[i]  = NormLongPhiBiasRes.make<TH1F>(Form("histo_norm_dz_phi_plot%i",i),
                                   Form("%.2f#circ<#varphi^{probe}_{tk}<%.2f#circ;d_{z}/#sigma_{d_{z}};tracks",phiF,phiL),
                                   mybins_,-dzmax_phi/100.,dzmax_phi/100.);

      n_dzEtaBiasResiduals[i]  = NormLongEtaBiasRes.make<TH1F>(Form("histo_norm_dz_eta_plot%i",i),
                                   Form("%.2f<#eta^{probe}_{tk}<%.2f;d_{z}/#sigma_{d_{z}};tracks",etaF,etaL),
                                   mybins_,-dzmax_eta/100.,dzmax_eta/100.);
      
      for ( int j=0; j<nBins_; ++j ) {
    
    a_dxyBiasResidualsMap[i][j] = AbsDoubleDiffBiasRes.make<TH1F>(Form("histo_dxy_eta_plot%i_phi_plot%i",i,j),
                                      Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{xy} [#mum];tracks",etaF,etaL,phiF,phiL),
                                      mybins_,-dzmax_eta,dzmax_eta);
    
    a_dzBiasResidualsMap[i][j]  = AbsDoubleDiffBiasRes.make<TH1F>(Form("histo_dxy_eta_plot%i_phi_plot%i",i,j),
                                      Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{z} [#mum];tracks",etaF,etaL,phiF,phiL),
                                      mybins_,-dzmax_eta,dzmax_eta);
    
    n_dxyBiasResidualsMap[i][j] = NormDoubleDiffBiasRes.make<TH1F>(Form("histo_norm_dxy_eta_plot%i_phi_plot%i",i,j),
                                       Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{xy}/#sigma_{d_{xy}};tracks",etaF,etaL,phiF,phiL),
                                       mybins_,-dzmax_eta/100,dzmax_eta/100);

    n_dzBiasResidualsMap[i][j]  = NormDoubleDiffBiasRes.make<TH1F>(Form("histo_norm_dxy_eta_plot%i_phi_plot%i",i,j),
                                       Form("%.2f<#eta_{tk}<%.2f %.2f#circ<#varphi_{tk}<%.2f#circ;d_{z}/#sigma_{d_{z}};tracks",etaF,etaL,phiF,phiL),
                                       mybins_,-dzmax_eta/100,dzmax_eta/100);
      }
    }
    
    // declaration of the directories

    TFileDirectory MeanBiasTrendsDir   = fs->mkdir("MeanBiasTrends");
    TFileDirectory WidthBiasTrendsDir  = fs->mkdir("WidthBiasTrends");
    TFileDirectory MedianBiasTrendsDir = fs->mkdir("MedianBiasTrends");
    TFileDirectory MADBiasTrendsDir    = fs->mkdir("MADBiasTrends");
    
    TFileDirectory Mean2DBiasMapsDir   = fs->mkdir("MeanBiasMaps");
    TFileDirectory Width2DBiasMapsDir  = fs->mkdir("WidthBiasMaps");
    
    // means and widths from the fit
    
    a_dxyPhiMeanBiasTrend  = MeanBiasTrendsDir.make<TH1F> ("means_dxy_phi",
                               "#LT d_{xy} #GT vs #phi sector;#varphi (sector) [degrees];#LT d_{xy} #GT [#mum]",
                               nBins_,lowedge,highedge); 

    a_dxyPhiWidthBiasTrend = WidthBiasTrendsDir.make<TH1F>("widths_dxy_phi",
                               "#sigma_{d_{xy}} vs #phi sector;#varphi (sector) [degrees];#sigma_{d_{xy}} [#mum]",
                               nBins_,lowedge,highedge);

    a_dzPhiMeanBiasTrend   = MeanBiasTrendsDir.make<TH1F> ("means_dz_phi",
                               "#LT d_{z} #GT vs #phi sector;#varphi (sector) [degrees];#LT d_{z} #GT [#mum]",
                               nBins_,lowedge,highedge); 

    a_dzPhiWidthBiasTrend  = WidthBiasTrendsDir.make<TH1F>("widths_dz_phi",
                               "#sigma_{d_{z}} vs #phi sector;#varphi (sector) [degrees];#sigma_{d_{z}} [#mum]",
                               nBins_,lowedge,highedge);
    
    a_dxyEtaMeanBiasTrend  = MeanBiasTrendsDir.make<TH1F> ("means_dxy_eta",
                               "#LT d_{xy} #GT vs #eta sector;#eta (sector);#LT d_{xy} #GT [#mum]",
                               nBins_,lowedge,highedge);

    a_dxyEtaWidthBiasTrend = WidthBiasTrendsDir.make<TH1F>("widths_dxy_eta",
                               "#sigma_{d_{xy}} vs #eta sector;#eta (sector);#sigma_{d_{xy}} [#mum]",
                               nBins_,lowedge,highedge);

    a_dzEtaMeanBiasTrend   = MeanBiasTrendsDir.make<TH1F> ("means_dz_eta",
                               "#LT d_{z} #GT vs #eta sector;#eta (sector);#LT d_{z} #GT [#mum]",
                               nBins_,lowedge,highedge); 

    a_dzEtaWidthBiasTrend  = WidthBiasTrendsDir.make<TH1F>("widths_dz_eta",
                               "#sigma_{d_{xy}} vs #eta sector;#eta (sector);#sigma_{d_{z}} [#mum]",
                               nBins_,lowedge,highedge);
    
    n_dxyPhiMeanBiasTrend  = MeanBiasTrendsDir.make<TH1F> ("norm_means_dxy_phi",
                               "#LT d_{xy}/#sigma_{d_{xy}} #GT vs #phi sector;#varphi (sector) [degrees];#LT d_{xy}/#sigma_{d_{xy}} #GT",
                               nBins_,lowedge,highedge);

    n_dxyPhiWidthBiasTrend = WidthBiasTrendsDir.make<TH1F>("norm_widths_dxy_phi",
                               "width(d_{xy}/#sigma_{d_{xy}}) vs #phi sector;#varphi (sector) [degrees]; width(d_{xy}/#sigma_{d_{xy}})",
                               nBins_,lowedge,highedge);

    n_dzPhiMeanBiasTrend   = MeanBiasTrendsDir.make<TH1F> ("norm_means_dz_phi",
                               "#LT d_{z}/#sigma_{d_{z}} #GT vs #phi sector;#varphi (sector) [degrees];#LT d_{z}/#sigma_{d_{z}} #GT",
                               nBins_,lowedge,highedge); 

    n_dzPhiWidthBiasTrend  = WidthBiasTrendsDir.make<TH1F>("norm_widths_dz_phi",
                               "width(d_{z}/#sigma_{d_{z}}) vs #phi sector;#varphi (sector) [degrees];width(d_{z}/#sigma_{d_{z}})",
                               nBins_,lowedge,highedge);
    
    n_dxyEtaMeanBiasTrend  = MeanBiasTrendsDir.make<TH1F> ("norm_means_dxy_eta",
                               "#LT d_{xy}/#sigma_{d_{xy}} #GT vs #eta sector;#eta (sector);#LT d_{xy}/#sigma_{d_{z}} #GT",
                               nBins_,lowedge,highedge);

    n_dxyEtaWidthBiasTrend = WidthBiasTrendsDir.make<TH1F>("norm_widths_dxy_eta",
                               "width(d_{xy}/#sigma_{d_{xy}}) vs #eta sector;#eta (sector);width(d_{xy}/#sigma_{d_{z}})",
                               nBins_,lowedge,highedge);

    n_dzEtaMeanBiasTrend   = MeanBiasTrendsDir.make<TH1F> ("norm_means_dz_eta",
                               "#LT d_{z}/#sigma_{d_{z}} #GT vs #eta sector;#eta (sector);#LT d_{z}/#sigma_{d_{z}} #GT",
                               nBins_,lowedge,highedge);  

    n_dzEtaWidthBiasTrend  = WidthBiasTrendsDir.make<TH1F>("norm_widths_dz_eta",
                               "width(d_{z}/#sigma_{d_{z}}) vs #eta sector;#eta (sector);width(d_{z}/#sigma_{d_{z}})",
                               nBins_,lowedge,highedge);                        
    
    // 2D maps
    
    a_dxyMeanBiasMap       =  Mean2DBiasMapsDir.make<TH2F>  ("means_dxy_map",
                                 "#LT d_{xy} #GT map;#eta (sector);#varphi (sector) [degrees]",
                                 nBins_,lowedge,highedge,nBins_,lowedge,highedge);

    a_dzMeanBiasMap        =  Mean2DBiasMapsDir.make<TH2F>  ("means_dz_map",
                                 "#LT d_{z} #GT map;#eta (sector);#varphi (sector) [degrees]",
                                 nBins_,lowedge,highedge,nBins_,lowedge,highedge);
    
    n_dxyMeanBiasMap       =  Mean2DBiasMapsDir.make<TH2F>  ("norm_means_dxy_map",
                                 "#LT d_{xy}/#sigma_{d_{xy}} #GT map;#eta (sector);#varphi (sector) [degrees]",
                                 nBins_,lowedge,highedge,nBins_,lowedge,highedge);

    n_dzMeanBiasMap        =  Mean2DBiasMapsDir.make<TH2F>  ("norm_means_dz_map",
                                 "#LT d_{z}/#sigma_{d_{z}} #GT map;#eta (sector);#varphi (sector) [degrees]",
                                 nBins_,lowedge,highedge,nBins_,lowedge,highedge);
    
    a_dxyWidthBiasMap      =  Width2DBiasMapsDir.make<TH2F> ("widths_dxy_map",
                                 "#sigma_{d_{xy}} map;#eta (sector);#varphi (sector) [degrees]",
                                 nBins_,lowedge,highedge,nBins_,lowedge,highedge);

    a_dzWidthBiasMap       =  Width2DBiasMapsDir.make<TH2F> ("widths_dz_map",
                                 "#sigma_{d_{z}} map;#eta (sector);#varphi (sector) [degrees]",
                                 nBins_,lowedge,highedge,nBins_,lowedge,highedge);
    
    n_dxyWidthBiasMap      =  Width2DBiasMapsDir.make<TH2F> ("norm_widths_dxy_map",
                                 "width(d_{xy}/#sigma_{d_{xy}}) map;#eta (sector);#varphi (sector) [degrees]",
                                 nBins_,lowedge,highedge,nBins_,lowedge,highedge);

    n_dzWidthBiasMap       =  Width2DBiasMapsDir.make<TH2F> ("norm_widths_dz_map",
                                 "width(d_{z}/#sigma_{d_{z}}) map;#eta (sector);#varphi (sector) [degrees]",
                                 nBins_,lowedge,highedge,nBins_,lowedge,highedge);
    
    // medians and MADs
    
    a_dxyPhiMedianBiasTrend = MedianBiasTrendsDir.make<TH1F>("medians_dxy_phi",
                                 "Median of d_{xy} vs #phi sector;#varphi (sector) [degrees];#mu_{1/2}(d_{xy}) [#mum]",
                                 nBins_,lowedge,highedge);  
             
    a_dxyPhiMADBiasTrend    = MADBiasTrendsDir.make<TH1F>   ("MADs_dxy_phi",
                                 "Median absolute deviation of d_{xy} vs #phi sector;#varphi (sector) [degrees];MAD(d_{xy}) [#mum]",
                                 nBins_,lowedge,highedge);  
                
    a_dzPhiMedianBiasTrend  = MedianBiasTrendsDir.make<TH1F>("medians_dz_phi",
                                 "Median of d_{z} vs #phi sector;#varphi (sector) [degrees];#mu_{1/2}(d_{z}) [#mum]",
                                 nBins_,lowedge,highedge);      
            
    a_dzPhiMADBiasTrend     = MADBiasTrendsDir.make<TH1F>   ("MADs_dz_phi",
                                 "Median absolute deviation of d_{z} vs #phi sector;#varphi (sector) [degrees];MAD(d_{z}) [#mum]",
                                 nBins_,lowedge,highedge);                  
    
    a_dxyEtaMedianBiasTrend = MedianBiasTrendsDir.make<TH1F>("medians_dxy_eta",
                                 "Median of d_{xy} vs #eta sector;#eta (sector);#mu_{1/2}(d_{xy}) [#mum]",
                                 nBins_,lowedge,highedge);  
                      
    a_dxyEtaMADBiasTrend    = MADBiasTrendsDir.make<TH1F>   ("MADs_dxy_eta",
                                 "Median absolute deviation of d_{xy} vs #eta sector;#eta (sector);MAD(d_{xy}) [#mum]",
                                 nBins_,lowedge,highedge);      
                
    a_dzEtaMedianBiasTrend  = MedianBiasTrendsDir.make<TH1F>("medians_dz_eta",
                                 "Median of d_{z} vs #eta sector;#eta (sector);#mu_{1/2}(d_{z}) [#mum]",
                                 nBins_,lowedge,highedge);  
                      
    a_dzEtaMADBiasTrend     = MADBiasTrendsDir.make<TH1F>   ("MADs_dz_eta",
                                 "Median absolute deviation of d_{z} vs #eta sector;#eta (sector);MAD(d_{z}) [#mum]",
                                 nBins_,lowedge,highedge);                        
    
    n_dxyPhiMedianBiasTrend = MedianBiasTrendsDir.make<TH1F>("norm_medians_dxy_phi",
                                 "Median of d_{xy}/#sigma_{d_{xy}} vs #phi sector;#varphi (sector) [degrees];#mu_{1/2}(d_{xy}/#sigma_{d_{xy}})",
                                 nBins_,lowedge,highedge); 

    n_dxyPhiMADBiasTrend    = MADBiasTrendsDir.make<TH1F>   ("norm_MADs_dxy_phi",
                                 "Median absolute deviation of d_{xy}/#sigma_{d_{xy}} vs #phi sector;#varphi (sector) [degrees]; MAD(d_{xy}/#sigma_{d_{xy}})",
                                 nBins_,lowedge,highedge); 
  
    n_dzPhiMedianBiasTrend  = MedianBiasTrendsDir.make<TH1F>("norm_medians_dz_phi",
                                 "Median of d_{z}/#sigma_{d_{z}} vs #phi sector;#varphi (sector) [degrees];#mu_{1/2}(d_{z}/#sigma_{d_{z}})",
                                 nBins_,lowedge,highedge); 
     
    n_dzPhiMADBiasTrend     = MADBiasTrendsDir.make<TH1F>   ("norm_MADs_dz_phi",
                                 "Median absolute deviation of d_{z}/#sigma_{d_{z}} vs #phi sector;#varphi (sector) [degrees];MAD(d_{z}/#sigma_{d_{z}})",
                                 nBins_,lowedge,highedge);      
    
    n_dxyEtaMedianBiasTrend = MedianBiasTrendsDir.make<TH1F>("norm_medians_dxy_eta",
                                 "Median of d_{xy}/#sigma_{d_{xy}} vs #eta sector;#eta (sector);#mu_{1/2}(d_{xy}/#sigma_{d_{z}})",
                                 nBins_,lowedge,highedge);  
        
    n_dxyEtaMADBiasTrend    = MADBiasTrendsDir.make<TH1F>   ("norm_MADs_dxy_eta",
                                 "Median absolute deviation of d_{xy}/#sigma_{d_{xy}} vs #eta sector;#eta (sector);MAD(d_{xy}/#sigma_{d_{z}})",
                                 nBins_,lowedge,highedge);  
        
    n_dzEtaMedianBiasTrend  = MedianBiasTrendsDir.make<TH1F>("norm_medians_dz_eta",
                                 "Median of d_{z}/#sigma_{d_{z}} vs #eta sector;#eta (sector);#mu_{1/2}(d_{z}/#sigma_{d_{z}})",
                                 nBins_,lowedge,highedge);  
            
    n_dzEtaMADBiasTrend     = MADBiasTrendsDir.make<TH1F>   ("norm_MADs_dz_eta",
                                 "Median absolute deviation of d_{z}/#sigma_{d_{z}} vs #eta sector;#eta (sector);MAD(d_{z}/#sigma_{d_{z}})",
                                 nBins_,lowedge,highedge);                      
    
  }
}

std::vector< TH1F * > PrimaryVertexValidation::bookResidualsHistogram ( TFileDirectory  dir, unsigned int  theNOfBins, TString  resType, TString  varType  )

private

Definition at line 2603 of file PrimaryVertexValidation.cc.

References down, h, mps_fire::i, TFileDirectory::make(), and up.

Referenced by beginJob().

                                                   {
  TH1F::SetDefaultSumw2(kTRUE);
  
  Double_t up   = 1000;
  Double_t down = -up;
  
  if(resType.Contains("norm")){
    up = up*(1/100);
    down = down*(1/100);
  }
  
  std::vector<TH1F*> h;
  h.reserve(theNOfBins);
  
  const char* auxResType = (resType.ReplaceAll("_","")).Data();
  
  for(unsigned int i=0; i<theNOfBins;i++){
    TH1F* htemp = dir.make<TH1F>(Form("histo_%s_%s_plot%i",resType.Data(),varType.Data(),i),
                 Form("%s vs %s - bin %i;%s;tracks",auxResType,varType.Data(),i,auxResType),
                 500,down,up); 
    h.push_back(htemp);
  }
  
  return h;
  
}

std::map< std::string, TH1 * > PrimaryVertexValidation::bookVertexHistograms ( TFileDirectory  dir )

private

Definition at line 2549 of file PrimaryVertexValidation.cc.

References h, TFileDirectory::make(), Pi, AlCaHLTBitMon_QueryRunRegistry::string, lumiQTWidget::t, and cutBasedMuonId_MuonPOG_V0_cff::types.

Referenced by beginJob().

{

  TH1F::SetDefaultSumw2(kTRUE);

  std::map<std::string, TH1*> h;
  
  // histograms of track quality (Data and MC)
  std::string types[] = {"all","sel"};
  for(int t=0; t<2; t++){
    h["pseudorapidity_"+types[t]] =dir.make <TH1F>(("rapidity_"+types[t]).c_str(),"track pseudorapidity; track #eta; tracks",100,-3., 3.);
    h["z0_"+types[t]] = dir.make<TH1F>(("z0_"+types[t]).c_str(),"track z_{0};track z_{0} (cm);tracks",80,-40., 40.);
    h["phi_"+types[t]] = dir.make<TH1F>(("phi_"+types[t]).c_str(),"track #phi; track #phi;tracks",80,-TMath::Pi(), TMath::Pi());
    h["eta_"+types[t]] = dir.make<TH1F>(("eta_"+types[t]).c_str(),"track #eta; track #eta;tracks",80,-4., 4.);
    h["pt_"+types[t]] = dir.make<TH1F>(("pt_"+types[t]).c_str(),"track p_{T}; track p_{T} [GeV];tracks",100,0., 20.);
    h["p_"+types[t]] = dir.make<TH1F>(("p_"+types[t]).c_str(),"track p; track p [GeV];tracks",100,0., 20.);
    h["found_"+types[t]] = dir.make<TH1F>(("found_"+types[t]).c_str(),"n. found hits;n^{found}_{hits};tracks",30, 0., 30.);
    h["lost_"+types[t]] = dir.make<TH1F>(("lost_"+types[t]).c_str(),"n. lost hits;n^{lost}_{hits};tracks",20, 0., 20.);
    h["nchi2_"+types[t]] = dir.make<TH1F>(("nchi2_"+types[t]).c_str(),"normalized track #chi^{2};track #chi^{2}/ndf;tracks",100, 0., 20.);
    h["rstart_"+types[t]] = dir.make<TH1F>(("rstart_"+types[t]).c_str(),"track start radius; track innermost radius r (cm);tracks",100, 0., 20.);
    h["expectedInner_"+types[t]] = dir.make<TH1F>(("expectedInner_"+types[t]).c_str(),"n. expected inner hits;n^{expected}_{inner};tracks",10, 0., 10.);
    h["expectedOuter_"+types[t]] = dir.make<TH1F>(("expectedOuter_"+types[t]).c_str(),"n. expected outer hits;n^{expected}_{outer};tracks ",10, 0., 10.);
    h["logtresxy_"+types[t]] = dir.make<TH1F>(("logtresxy_"+types[t]).c_str(),"log10(track r-#phi resolution/#mum);log10(track r-#phi resolution/#mum);tracks",100, 0., 5.);
    h["logtresz_"+types[t]] = dir.make<TH1F>(("logtresz_"+types[t]).c_str(),"log10(track z resolution/#mum);log10(track z resolution/#mum);tracks",100, 0., 5.);
    h["tpullxy_"+types[t]] = dir.make<TH1F>(("tpullxy_"+types[t]).c_str(),"track r-#phi pull;pull_{r-#phi};tracks",100, -10., 10.);
    h["tpullz_"+types[t]] = dir.make<TH1F>(("tpullz_"+types[t]).c_str(),"track r-z pull;pull_{r-z};tracks",100, -50., 50.);
    h["tlogDCAxy_"+types[t]] = dir.make<TH1F>(("tlogDCAxy_"+types[t]).c_str(),"track log_{10}(DCA_{r-#phi});track log_{10}(DCA_{r-#phi});tracks",200, -5., 3.);
    h["tlogDCAz_"+types[t]] = dir.make<TH1F>(("tlogDCAz_"+types[t]).c_str(),"track log_{10}(DCA_{r-z});track log_{10}(DCA_{r-z});tracks",200, -5., 5.);
    h["lvseta_"+types[t]] = dir.make<TH2F>(("lvseta_"+types[t]).c_str(),"cluster length vs #eta;track #eta;cluster length",60,-3., 3., 20, 0., 20);
    h["lvstanlambda_"+types[t]] = dir.make<TH2F>(("lvstanlambda_"+types[t]).c_str(),"cluster length vs tan #lambda; tan#lambda;cluster length",60,-6., 6., 20, 0., 20);
    h["restrkz_"+types[t]] = dir.make<TH1F>(("restrkz_"+types[t]).c_str(),"z-residuals (track vs vertex);res_{z} (cm);tracks", 200, -5., 5.);
    h["restrkzvsphi_"+types[t]] = dir.make<TH2F>(("restrkzvsphi_"+types[t]).c_str(),"z-residuals (track - vertex) vs track #phi;track #phi;res_{z} (cm)", 12,-TMath::Pi(),TMath::Pi(),100, -0.5,0.5);
    h["restrkzvseta_"+types[t]] = dir.make<TH2F>(("restrkzvseta_"+types[t]).c_str(),"z-residuals (track - vertex) vs track #eta;track #eta;res_{z} (cm)", 12,-3.,3.,200, -0.5,0.5);
    h["pulltrkzvsphi_"+types[t]] = dir.make<TH2F>(("pulltrkzvsphi_"+types[t]).c_str(),"normalized z-residuals (track - vertex) vs track #phi;track #phi;res_{z}/#sigma_{res_{z}}", 12,-TMath::Pi(),TMath::Pi(),100, -5., 5.);
    h["pulltrkzvseta_"+types[t]] = dir.make<TH2F>(("pulltrkzvseta_"+types[t]).c_str(),"normalized z-residuals (track - vertex) vs track #eta;track #eta;res_{z}/#sigma_{res_{z}}", 12,-3.,3.,100, -5., 5.);
    h["pulltrkz_"+types[t]] = dir.make<TH1F>(("pulltrkz_"+types[t]).c_str(),"normalized z-residuals (track vs vertex);res_{z}/#sigma_{res_{z}};tracks", 100, -5., 5.);
    h["sigmatrkz0_"+types[t]] = dir.make<TH1F>(("sigmatrkz0_"+types[t]).c_str(),"z-resolution (excluding beam);#sigma^{trk}_{z_{0}} (cm);tracks", 100, 0., 5.);
    h["sigmatrkz_"+types[t]] = dir.make<TH1F>(("sigmatrkz_"+types[t]).c_str(),"z-resolution (including beam);#sigma^{trk}_{z} (cm);tracks", 100,0., 5.);
    h["nbarrelhits_"+types[t]] = dir.make<TH1F>(("nbarrelhits_"+types[t]).c_str(),"number of pixel barrel hits;n. hits Barrel Pixel;tracks", 10, 0., 10.);
    h["nbarrelLayers_"+types[t]] = dir.make<TH1F>(("nbarrelLayers_"+types[t]).c_str(),"number of pixel barrel layers;n. layers Barrel Pixel;tracks", 10, 0., 10.);
    h["nPxLayers_"+types[t]] = dir.make<TH1F>(("nPxLayers_"+types[t]).c_str(),"number of pixel layers (barrel+endcap);n. Pixel layers;tracks", 10, 0., 10.);
    h["nSiLayers_"+types[t]] = dir.make<TH1F>(("nSiLayers_"+types[t]).c_str(),"number of Tracker layers;n. Tracker layers;tracks", 20, 0., 20.);
    h["trackAlgo_"+types[t]] = dir.make<TH1F>(("trackAlgo_"+types[t]).c_str(),"track algorithm;track algo;tracks", 30, 0., 30.);
    h["trackQuality_"+types[t]] = dir.make<TH1F>(("trackQuality_"+types[t]).c_str(),"track quality;track quality;tracks", 7, -1., 6.);
  }

  return h;
  
}

void PrimaryVertexValidation::endJob ( void  )

privatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 1994 of file PrimaryVertexValidation.cc.

References a_dxyBiasResidualsMap, a_dxyEtaBiasResiduals, a_dxyEtaMADBiasTrend, a_dxyEtaMADTrend, a_dxyEtaMeanBiasTrend, a_dxyEtaMeanTrend, a_dxyEtaMedianBiasTrend, a_dxyEtaMedianTrend, a_dxyEtaResiduals, a_dxyEtaWidthBiasTrend, a_dxyEtaWidthTrend, a_dxyMeanBiasMap, a_dxyMeanMap, a_dxyPhiBiasResiduals, a_dxyPhiMADBiasTrend, a_dxyPhiMADTrend, a_dxyPhiMeanBiasTrend, a_dxyPhiMeanTrend, a_dxyPhiMedianBiasTrend, a_dxyPhiMedianTrend, a_dxyPhiResiduals, a_dxyPhiWidthBiasTrend, a_dxyPhiWidthTrend, a_dxypTCentralMeanTrend, a_dxypTCentralWidthTrend, a_dxypTMeanTrend, a_dxypTWidthTrend, a_dxyResidualsMap, a_dxyWidthBiasMap, a_dxyWidthMap, a_dzBiasResidualsMap, a_dzEtaBiasResiduals, a_dzEtaMADBiasTrend, a_dzEtaMADTrend, a_dzEtaMeanBiasTrend, a_dzEtaMeanTrend, a_dzEtaMedianBiasTrend, a_dzEtaMedianTrend, a_dzEtaResiduals, a_dzEtaWidthBiasTrend, a_dzEtaWidthTrend, a_dzMeanBiasMap, a_dzMeanMap, a_dzPhiBiasResiduals, a_dzPhiMADBiasTrend, a_dzPhiMADTrend, a_dzPhiMeanBiasTrend, a_dzPhiMeanTrend, a_dzPhiMedianBiasTrend, a_dzPhiMedianTrend, a_dzPhiResiduals, a_dzPhiWidthBiasTrend, a_dzPhiWidthTrend, a_dzpTCentralMeanTrend, a_dzpTCentralWidthTrend, a_dzpTMeanTrend, a_dzpTWidthTrend, a_dzResidualsMap, a_dzWidthBiasMap, a_dzWidthMap, fillMap(), fillTrendPlot(), fillTrendPlotByIndex(), h_dxy_Central_pT_, h_dxy_pT_, h_dz_Central_pT_, h_dz_pT_, h_norm_dxy_Central_pT_, h_norm_dxy_pT_, h_norm_dz_Central_pT_, h_norm_dz_pT_, statmode::MAD, statmode::MEAN, statmode::MEDIAN, n_dxyBiasResidualsMap, n_dxyEtaBiasResiduals, n_dxyEtaMADBiasTrend, n_dxyEtaMADTrend, n_dxyEtaMeanBiasTrend, n_dxyEtaMeanTrend, n_dxyEtaMedianBiasTrend, n_dxyEtaMedianTrend, n_dxyEtaResiduals, n_dxyEtaWidthBiasTrend, n_dxyEtaWidthTrend, n_dxyMeanBiasMap, n_dxyMeanMap, n_dxyPhiBiasResiduals, n_dxyPhiMADBiasTrend, n_dxyPhiMADTrend, n_dxyPhiMeanBiasTrend, n_dxyPhiMeanTrend, n_dxyPhiMedianBiasTrend, n_dxyPhiMedianTrend, n_dxyPhiResiduals, n_dxyPhiWidthBiasTrend, n_dxyPhiWidthTrend, n_dxypTCentralMeanTrend, n_dxypTCentralWidthTrend, n_dxypTMeanTrend, n_dxypTWidthTrend, n_dxyResidualsMap, n_dxyWidthBiasMap, n_dxyWidthMap, n_dzBiasResidualsMap, n_dzEtaBiasResiduals, n_dzEtaMADBiasTrend, n_dzEtaMADTrend, n_dzEtaMeanBiasTrend, n_dzEtaMeanTrend, n_dzEtaMedianBiasTrend, n_dzEtaMedianTrend, n_dzEtaResiduals, n_dzEtaWidthBiasTrend, n_dzEtaWidthTrend, n_dzMeanBiasMap, n_dzMeanMap, n_dzPhiBiasResiduals, n_dzPhiMADBiasTrend, n_dzPhiMADTrend, n_dzPhiMeanBiasTrend, n_dzPhiMeanTrend, n_dzPhiMedianBiasTrend, n_dzPhiMedianTrend, n_dzPhiResiduals, n_dzPhiWidthBiasTrend, n_dzPhiWidthTrend, n_dzpTCentralMeanTrend, n_dzpTCentralWidthTrend, n_dzpTMeanTrend, n_dzpTWidthTrend, n_dzResidualsMap, n_dzWidthBiasMap, n_dzWidthMap, Nevt_, useTracksFromRecoVtx_, and statmode::WIDTH.

Referenced by o2olib.O2ORunMgr::executeJob().

{

  edm::LogInfo("PrimaryVertexValidation")
    <<"######################################\n"
    <<"# PrimaryVertexValidation::endJob()\n" 
    <<"# Number of analyzed events: "<<Nevt_<<"\n"
    <<"######################################";

  if(useTracksFromRecoVtx_){

    fillTrendPlot(a_dxyPhiMeanBiasTrend ,a_dxyPhiBiasResiduals,statmode::MEAN,"phi");  
    fillTrendPlot(a_dxyPhiWidthBiasTrend,a_dxyPhiBiasResiduals,statmode::WIDTH,"phi");
    fillTrendPlot(a_dzPhiMeanBiasTrend  ,a_dzPhiBiasResiduals ,statmode::MEAN,"phi");   
    fillTrendPlot(a_dzPhiWidthBiasTrend ,a_dzPhiBiasResiduals ,statmode::WIDTH,"phi");  
    
    fillTrendPlot(a_dxyEtaMeanBiasTrend ,a_dxyEtaBiasResiduals,statmode::MEAN,"eta"); 
    fillTrendPlot(a_dxyEtaWidthBiasTrend,a_dxyEtaBiasResiduals,statmode::WIDTH,"eta");
    fillTrendPlot(a_dzEtaMeanBiasTrend  ,a_dzEtaBiasResiduals ,statmode::MEAN,"eta"); 
    fillTrendPlot(a_dzEtaWidthBiasTrend ,a_dzEtaBiasResiduals ,statmode::WIDTH,"eta");
    
    fillTrendPlot(n_dxyPhiMeanBiasTrend ,n_dxyPhiBiasResiduals,statmode::MEAN,"phi"); 
    fillTrendPlot(n_dxyPhiWidthBiasTrend,n_dxyPhiBiasResiduals,statmode::WIDTH,"phi");
    fillTrendPlot(n_dzPhiMeanBiasTrend  ,n_dzPhiBiasResiduals ,statmode::MEAN,"phi"); 
    fillTrendPlot(n_dzPhiWidthBiasTrend ,n_dzPhiBiasResiduals ,statmode::WIDTH,"phi");
    
    fillTrendPlot(n_dxyEtaMeanBiasTrend ,n_dxyEtaBiasResiduals,statmode::MEAN,"eta"); 
    fillTrendPlot(n_dxyEtaWidthBiasTrend,n_dxyEtaBiasResiduals,statmode::WIDTH,"eta");
    fillTrendPlot(n_dzEtaMeanBiasTrend  ,n_dzEtaBiasResiduals ,statmode::MEAN,"eta"); 
    fillTrendPlot(n_dzEtaWidthBiasTrend ,n_dzEtaBiasResiduals ,statmode::WIDTH,"eta");
    
    // medians and MADs   
    
    fillTrendPlot(a_dxyPhiMedianBiasTrend,a_dxyPhiBiasResiduals,statmode::MEDIAN,"phi");  
    fillTrendPlot(a_dxyPhiMADBiasTrend   ,a_dxyPhiBiasResiduals,statmode::MAD,"phi"); 
    fillTrendPlot(a_dzPhiMedianBiasTrend ,a_dzPhiBiasResiduals ,statmode::MEDIAN,"phi");  
    fillTrendPlot(a_dzPhiMADBiasTrend    ,a_dzPhiBiasResiduals ,statmode::MAD,"phi"); 
    
    fillTrendPlot(a_dxyEtaMedianBiasTrend,a_dxyEtaBiasResiduals,statmode::MEDIAN,"eta");  
    fillTrendPlot(a_dxyEtaMADBiasTrend   ,a_dxyEtaBiasResiduals,statmode::MAD,"eta"); 
    fillTrendPlot(a_dzEtaMedianBiasTrend ,a_dzEtaBiasResiduals ,statmode::MEDIAN,"eta");  
    fillTrendPlot(a_dzEtaMADBiasTrend    ,a_dzEtaBiasResiduals ,statmode::MAD,"eta"); 
    
    fillTrendPlot(n_dxyPhiMedianBiasTrend,n_dxyPhiBiasResiduals,statmode::MEDIAN,"phi");  
    fillTrendPlot(n_dxyPhiMADBiasTrend   ,n_dxyPhiBiasResiduals,statmode::MAD,"phi"); 
    fillTrendPlot(n_dzPhiMedianBiasTrend ,n_dzPhiBiasResiduals ,statmode::MEDIAN,"phi");  
    fillTrendPlot(n_dzPhiMADBiasTrend    ,n_dzPhiBiasResiduals ,statmode::MAD,"phi"); 
    
    fillTrendPlot(n_dxyEtaMedianBiasTrend,n_dxyEtaBiasResiduals,statmode::MEDIAN,"eta");  
    fillTrendPlot(n_dxyEtaMADBiasTrend   ,n_dxyEtaBiasResiduals,statmode::MAD,"eta"); 
    fillTrendPlot(n_dzEtaMedianBiasTrend ,n_dzEtaBiasResiduals ,statmode::MEDIAN,"eta");  
    fillTrendPlot(n_dzEtaMADBiasTrend    ,n_dzEtaBiasResiduals ,statmode::MAD,"eta"); 
   
    // 2d Maps

    fillMap(a_dxyMeanBiasMap ,a_dxyBiasResidualsMap,statmode::MEAN); 
    fillMap(a_dxyWidthBiasMap,a_dxyBiasResidualsMap,statmode::WIDTH);
    fillMap(a_dzMeanBiasMap  ,a_dzBiasResidualsMap,statmode::MEAN); 
    fillMap(a_dzWidthBiasMap ,a_dzBiasResidualsMap,statmode::WIDTH);

    fillMap(n_dxyMeanBiasMap ,n_dxyBiasResidualsMap,statmode::MEAN); 
    fillMap(n_dxyWidthBiasMap,n_dxyBiasResidualsMap,statmode::WIDTH);
    fillMap(n_dzMeanBiasMap  ,n_dzBiasResidualsMap,statmode::MEAN); 
    fillMap(n_dzWidthBiasMap ,n_dzBiasResidualsMap,statmode::WIDTH);
   
  }

  // do profiles

  fillTrendPlot(a_dxyPhiMeanTrend ,a_dxyPhiResiduals,statmode::MEAN,"phi");  
  fillTrendPlot(a_dxyPhiWidthTrend,a_dxyPhiResiduals,statmode::WIDTH,"phi");
  fillTrendPlot(a_dzPhiMeanTrend  ,a_dzPhiResiduals ,statmode::MEAN,"phi");   
  fillTrendPlot(a_dzPhiWidthTrend ,a_dzPhiResiduals ,statmode::WIDTH,"phi");  
  
  fillTrendPlot(a_dxyEtaMeanTrend ,a_dxyEtaResiduals,statmode::MEAN,"eta"); 
  fillTrendPlot(a_dxyEtaWidthTrend,a_dxyEtaResiduals,statmode::WIDTH,"eta");
  fillTrendPlot(a_dzEtaMeanTrend  ,a_dzEtaResiduals ,statmode::MEAN,"eta"); 
  fillTrendPlot(a_dzEtaWidthTrend ,a_dzEtaResiduals ,statmode::WIDTH,"eta");
  
  fillTrendPlot(n_dxyPhiMeanTrend ,n_dxyPhiResiduals,statmode::MEAN,"phi"); 
  fillTrendPlot(n_dxyPhiWidthTrend,n_dxyPhiResiduals,statmode::WIDTH,"phi");
  fillTrendPlot(n_dzPhiMeanTrend  ,n_dzPhiResiduals ,statmode::MEAN,"phi"); 
  fillTrendPlot(n_dzPhiWidthTrend ,n_dzPhiResiduals ,statmode::WIDTH,"phi");
  
  fillTrendPlot(n_dxyEtaMeanTrend ,n_dxyEtaResiduals,statmode::MEAN,"eta"); 
  fillTrendPlot(n_dxyEtaWidthTrend,n_dxyEtaResiduals,statmode::WIDTH,"eta");
  fillTrendPlot(n_dzEtaMeanTrend  ,n_dzEtaResiduals ,statmode::MEAN,"eta"); 
  fillTrendPlot(n_dzEtaWidthTrend ,n_dzEtaResiduals ,statmode::WIDTH,"eta");
    
  // vs transverse momentum

  fillTrendPlotByIndex(a_dxypTMeanTrend ,h_dxy_pT_,statmode::MEAN );  
  fillTrendPlotByIndex(a_dxypTWidthTrend,h_dxy_pT_,statmode::WIDTH);
  fillTrendPlotByIndex(a_dzpTMeanTrend  ,h_dz_pT_ ,statmode::MEAN );   
  fillTrendPlotByIndex(a_dzpTWidthTrend ,h_dz_pT_ ,statmode::WIDTH);  
  
  fillTrendPlotByIndex(a_dxypTCentralMeanTrend ,h_dxy_Central_pT_,statmode::MEAN ); 
  fillTrendPlotByIndex(a_dxypTCentralWidthTrend,h_dxy_Central_pT_,statmode::WIDTH);
  fillTrendPlotByIndex(a_dzpTCentralMeanTrend  ,h_dz_Central_pT_ ,statmode::MEAN ); 
  fillTrendPlotByIndex(a_dzpTCentralWidthTrend ,h_dz_Central_pT_ ,statmode::WIDTH);
  
  fillTrendPlotByIndex(n_dxypTMeanTrend ,h_norm_dxy_pT_,statmode::MEAN ); 
  fillTrendPlotByIndex(n_dxypTWidthTrend,h_norm_dxy_pT_,statmode::WIDTH);
  fillTrendPlotByIndex(n_dzpTMeanTrend  ,h_norm_dz_pT_ ,statmode::MEAN ); 
  fillTrendPlotByIndex(n_dzpTWidthTrend ,h_norm_dz_pT_ ,statmode::WIDTH);
  
  fillTrendPlotByIndex(n_dxypTCentralMeanTrend ,h_norm_dxy_Central_pT_,statmode::MEAN ); 
  fillTrendPlotByIndex(n_dxypTCentralWidthTrend,h_norm_dxy_Central_pT_,statmode::WIDTH);
  fillTrendPlotByIndex(n_dzpTCentralMeanTrend  ,h_norm_dz_Central_pT_ ,statmode::MEAN ); 
  fillTrendPlotByIndex(n_dzpTCentralWidthTrend ,h_norm_dz_Central_pT_ ,statmode::WIDTH);

  // medians and MADs     
  
  fillTrendPlot(a_dxyPhiMedianTrend,a_dxyPhiResiduals,statmode::MEDIAN,"phi");  
  fillTrendPlot(a_dxyPhiMADTrend   ,a_dxyPhiResiduals,statmode::MAD,"phi"); 
  fillTrendPlot(a_dzPhiMedianTrend ,a_dzPhiResiduals ,statmode::MEDIAN,"phi");  
  fillTrendPlot(a_dzPhiMADTrend    ,a_dzPhiResiduals ,statmode::MAD,"phi"); 
  
  fillTrendPlot(a_dxyEtaMedianTrend,a_dxyEtaResiduals,statmode::MEDIAN,"eta");  
  fillTrendPlot(a_dxyEtaMADTrend   ,a_dxyEtaResiduals,statmode::MAD,"eta"); 
  fillTrendPlot(a_dzEtaMedianTrend ,a_dzEtaResiduals ,statmode::MEDIAN,"eta");  
  fillTrendPlot(a_dzEtaMADTrend    ,a_dzEtaResiduals ,statmode::MAD,"eta"); 
  
  fillTrendPlot(n_dxyPhiMedianTrend,n_dxyPhiResiduals,statmode::MEDIAN,"phi");  
  fillTrendPlot(n_dxyPhiMADTrend   ,n_dxyPhiResiduals,statmode::MAD,"phi"); 
  fillTrendPlot(n_dzPhiMedianTrend ,n_dzPhiResiduals ,statmode::MEDIAN,"phi");  
  fillTrendPlot(n_dzPhiMADTrend    ,n_dzPhiResiduals ,statmode::MAD,"phi"); 
  
  fillTrendPlot(n_dxyEtaMedianTrend,n_dxyEtaResiduals,statmode::MEDIAN,"eta");  
  fillTrendPlot(n_dxyEtaMADTrend   ,n_dxyEtaResiduals,statmode::MAD,"eta"); 
  fillTrendPlot(n_dzEtaMedianTrend ,n_dzEtaResiduals ,statmode::MEDIAN,"eta");  
  fillTrendPlot(n_dzEtaMADTrend    ,n_dzEtaResiduals ,statmode::MAD,"eta"); 
    
  // 2D Maps

  fillMap(a_dxyMeanMap ,a_dxyResidualsMap,statmode::MEAN); 
  fillMap(a_dxyWidthMap,a_dxyResidualsMap,statmode::WIDTH);
  fillMap(a_dzMeanMap  ,a_dzResidualsMap,statmode::MEAN); 
  fillMap(a_dzWidthMap ,a_dzResidualsMap,statmode::WIDTH);
  
  fillMap(n_dxyMeanMap ,n_dxyResidualsMap,statmode::MEAN); 
  fillMap(n_dxyWidthMap,n_dxyResidualsMap,statmode::WIDTH);
  fillMap(n_dzMeanMap  ,n_dzResidualsMap,statmode::MEAN); 
  fillMap(n_dzWidthMap ,n_dzResidualsMap,statmode::WIDTH);

}

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

private

Definition at line 2756 of file PrimaryVertexValidation.cc.

References alignCSCRings::s.

Referenced by fillTrackHistos().

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

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

private

Definition at line 2767 of file PrimaryVertexValidation.cc.

References alignCSCRings::s.

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

void PrimaryVertexValidation::fillByIndex ( std::vector< TH1F * > &  h, unsigned int  index, double  x  )

private

Definition at line 2778 of file PrimaryVertexValidation.cc.

References DEFINE_FWK_MODULE, and diffTreeTool::index.

Referenced by analyze().

{
  if(index <= h.size()){
    h[index]->Fill(x);
  } else {
    edm::LogWarning("PrimaryVertexValidation") << "Trying to fill non-existing Histogram with index " << index << std::endl;
    return;
  }
}

void PrimaryVertexValidation::fillMap ( TH2F *  trendMap, TH1F *  residualsMapPlot[100][100], statmode::estimator  fitPar_  )

private

Definition at line 2449 of file PrimaryVertexValidation.cc.

References etaOfProbe_, etaSect_, fitResiduals(), getMAD(), getMedian(), mps_fire::i, statmode::MAD, statmode::MEAN, statmode::MEDIAN, nBins_, phiSect_, Pi, and statmode::WIDTH.

Referenced by endJob().

{
 
  for ( int i=0; i<nBins_; ++i ) {
    
    char phipositionString[129];
    float phiInterval = phiSect_*(180/TMath::Pi());
    float phiposition = (-180+i*phiInterval)+(phiInterval/2);
    sprintf(phipositionString,"%.f",phiposition);
    
    trendMap->GetYaxis()->SetBinLabel(i+1,phipositionString); 

    for ( int j=0; j<nBins_; ++j ) {

      char etapositionString[129];
      float etaposition = (-etaOfProbe_+j*etaSect_)+(etaSect_/2);
      sprintf(etapositionString,"%.1f",etaposition);

      if(i==0) { trendMap->GetXaxis()->SetBinLabel(j+1,etapositionString); }

      switch (fitPar_)
    { 
    case statmode::MEAN:
      {
        float mean_      = fitResiduals(residualsMapPlot[i][j]).first.first;
        float meanErr_   = fitResiduals(residualsMapPlot[i][j]).first.second;
        trendMap->SetBinContent(j+1,i+1,mean_);
        trendMap->SetBinError(j+1,i+1,meanErr_);
        break;
      }
    case statmode::WIDTH:
      {
        float width_     = fitResiduals(residualsMapPlot[i][j]).second.first;
        float widthErr_  = fitResiduals(residualsMapPlot[i][j]).second.second;
        trendMap->SetBinContent(j+1,i+1,width_);
        trendMap->SetBinError(j+1,i+1,widthErr_);
        break;
      }     
    case statmode::MEDIAN:
      {
        float median_    = getMedian(residualsMapPlot[i][j]).first;
        float medianErr_ = getMedian(residualsMapPlot[i][j]).second;
        trendMap->SetBinContent(j+1,i+1,median_);
        trendMap->SetBinError(j+1,i+1,medianErr_);
        break;
      }     
    case statmode::MAD:
      {
        float mad_       = getMAD(residualsMapPlot[i][j]).first; 
        float madErr_    = getMAD(residualsMapPlot[i][j]).second;
        trendMap->SetBinContent(j+1,i+1,mad_);
        trendMap->SetBinError(j+1,i+1,madErr_);
        break;
      }     
    default:
      edm::LogWarning("PrimaryVertexValidation:") <<" fillMap() "<<fitPar_<<" unknown estimator!"<<std::endl;
    }   
    } // closes loop on eta bins
  } // cloeses loop on phi bins
}

void PrimaryVertexValidation::fillTrackHistos ( std::map< std::string, TH1 * > &  h, const std::string &  ttype, const reco::TransientTrack *  tt, const reco::Vertex &  v, const reco::BeamSpot &  beamSpot, double  fBfield  )

private

Definition at line 2634 of file PrimaryVertexValidation.cc.

References reco::TrackBase::algo(), Reference_intrackfit_cff::barrel, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), funct::cos(), allConversions_cfi::d0, reco::TrackBase::d0Error(), reco::TrackBase::dxy(), allConversions_cfi::dz, reco::TrackBase::dz(), reco::TrackBase::dzError(), reco::TrackBase::eta(), fill(), objects.autophobj::float, reco::Track::found(), reco::TrackBase::hitPattern(), reco::Track::innerPosition(), edm::Ref< C, T, F >::isNonnull(), reco::Vertex::isValid(), kappa, reco::TrackBase::lambda(), cmsBatch::log, reco::Track::lost(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfHits(), reco::TrackBase::p(), reco::TrackBase::phi(), PixelSubdetector::PixelBarrel, reco::HitPattern::pixelBarrelLayersWithMeasurement(), reco::HitPattern::pixelLayersWithMeasurement(), reco::BeamSpot::position(), reco::Vertex::position(), position, funct::pow(), reco::TrackBase::pt(), lumiQueryAPI::q, reco::TrackBase::qoverp(), reco::TrackBase::qualityMask(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), funct::sin(), mathSSE::sqrt(), reco::TransientTrack::stateAtBeamLine(), funct::tan(), reco::TrackBase::theta(), reco::TransientTrack::track(), DetId::Tracker, reco::HitPattern::trackerLayersWithMeasurement(), TrajectoryStateClosestToBeamLine::trackStateAtPCA(), reco::TrackBase::vx(), reco::TrackBase::vy(), reco::TrackBase::vz(), wxy2_, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and z.

Referenced by analyze().

{

  using namespace reco;

  fill(h,"pseudorapidity_"+ttype,tt->track().eta());
  fill(h,"z0_"+ttype,tt->track().vz());
  fill(h,"phi_"+ttype,tt->track().phi());
  fill(h,"eta_"+ttype,tt->track().eta());
  fill(h,"pt_"+ttype,tt->track().pt());
  fill(h,"p_"+ttype,tt->track().p());
  fill(h,"found_"+ttype,tt->track().found());
  fill(h,"lost_"+ttype,tt->track().lost());
  fill(h,"nchi2_"+ttype,tt->track().normalizedChi2());
  fill(h,"rstart_"+ttype,(tt->track().innerPosition()).Rho());
  
  double d0Error=tt->track().d0Error();
  double d0=tt->track().dxy(beamSpot.position());
  double dz=tt->track().dz(beamSpot.position());
  if (d0Error>0){
    fill(h,"logtresxy_"+ttype,log(d0Error/0.0001)/log(10.));
    fill(h,"tpullxy_"+ttype,d0/d0Error);
    fill(h,"tlogDCAxy_"+ttype,log(fabs(d0/d0Error)));
    
  }
  //double z0=tt->track().vz();
  double dzError=tt->track().dzError();
  if(dzError>0){
    fill(h,"logtresz_"+ttype,log(dzError/0.0001)/log(10.));
    fill(h,"tpullz_"+ttype,dz/dzError);
    fill(h,"tlogDCAz_"+ttype,log(fabs(dz/dzError)));
  }
  
  //
  double wxy2_=pow(beamSpot.BeamWidthX(),2)+pow(beamSpot.BeamWidthY(),2);

  fill(h,"sigmatrkz_"+ttype,sqrt(pow(tt->track().dzError(),2)+wxy2_/pow(tan(tt->track().theta()),2)));
  fill(h,"sigmatrkz0_"+ttype,tt->track().dzError());
  
  // track vs vertex
  if( v.isValid()){ // && (v.ndof()<10.)) {
    // emulate clusterizer input
    //const TransientTrack & tt = theB_->build(&t); wrong !!!!
    //reco::TransientTrack tt = theB_->build(&t); 
    //ttt->track().setBeamSpot(beamSpot); // need the setBeamSpot !
    double z=(tt->stateAtBeamLine().trackStateAtPCA()).position().z();
    double tantheta=tan((tt->stateAtBeamLine().trackStateAtPCA()).momentum().theta());
    double dz2= pow(tt->track().dzError(),2)+wxy2_/pow(tantheta,2);
    
    fill(h,"restrkz_"+ttype,z-v.position().z());
    fill(h,"restrkzvsphi_"+ttype,tt->track().phi(), z-v.position().z());
    fill(h,"restrkzvseta_"+ttype,tt->track().eta(), z-v.position().z());
    fill(h,"pulltrkzvsphi_"+ttype,tt->track().phi(), (z-v.position().z())/sqrt(dz2));
    fill(h,"pulltrkzvseta_"+ttype,tt->track().eta(), (z-v.position().z())/sqrt(dz2));
    
    fill(h,"pulltrkz_"+ttype,(z-v.position().z())/sqrt(dz2));
    
    double x1=tt->track().vx()-beamSpot.x0(); double y1=tt->track().vy()-beamSpot.y0();
    
    double kappa=-0.002998*fBfield_*tt->track().qoverp()/cos(tt->track().theta());
    double D0=x1*sin(tt->track().phi())-y1*cos(tt->track().phi())-0.5*kappa*(x1*x1+y1*y1);
    double q=sqrt(1.-2.*kappa*D0);
    double s0=(x1*cos(tt->track().phi())+y1*sin(tt->track().phi()))/q;
    // double s1;
    if (fabs(kappa*s0)>0.001){
      //s1=asin(kappa*s0)/kappa;
    }else{
      //double ks02=(kappa*s0)*(kappa*s0);
      //s1=s0*(1.+ks02/6.+3./40.*ks02*ks02+5./112.*pow(ks02,3));
    }
    // sp.ddcap=-2.*D0/(1.+q);
    //double zdcap=tt->track().vz()-s1/tan(tt->track().theta());
    
  }
  //
  
  // collect some info on hits and clusters
  fill(h,"nbarrelLayers_"+ttype,tt->track().hitPattern().pixelBarrelLayersWithMeasurement());
  fill(h,"nPxLayers_"+ttype,tt->track().hitPattern().pixelLayersWithMeasurement());
  fill(h,"nSiLayers_"+ttype,tt->track().hitPattern().trackerLayersWithMeasurement());
  fill(h,"expectedInner_"+ttype,tt->track().hitPattern().numberOfHits(HitPattern::MISSING_INNER_HITS));
  fill(h,"expectedOuter_"+ttype,tt->track().hitPattern().numberOfHits(HitPattern::MISSING_OUTER_HITS));
  fill(h,"trackAlgo_"+ttype,tt->track().algo());
  fill(h,"trackQuality_"+ttype,tt->track().qualityMask());
  
  //
  int longesthit=0, nbarrel=0;
  for(trackingRecHit_iterator hit=tt->track().recHitsBegin(); hit!=tt->track().recHitsEnd(); hit++){
    if ((**hit).isValid() && (**hit).geographicalId().det() == DetId::Tracker ){
      bool barrel = DetId((**hit).geographicalId()).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
      //bool endcap = DetId::DetId((**hit).geographicalId()).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
      if (barrel){
    const SiPixelRecHit *pixhit = dynamic_cast<const SiPixelRecHit*>( &(**hit));
    edm::Ref<edmNew::DetSetVector<SiPixelCluster>, SiPixelCluster> const& clust = (*pixhit).cluster();
    if (clust.isNonnull()) {
      nbarrel++;
      if (clust->sizeY()-longesthit>0) longesthit=clust->sizeY();
      if (clust->sizeY()>20.){
        fill(h,"lvseta_"+ttype,tt->track().eta(), 19.9);
        fill(h,"lvstanlambda_"+ttype,tan(tt->track().lambda()), 19.9);
      }else{
        fill(h,"lvseta_"+ttype,tt->track().eta(), float(clust->sizeY()));
        fill(h,"lvstanlambda_"+ttype,tan(tt->track().lambda()), float(clust->sizeY()));
      }
    }
      }
    }
  }
  fill(h,"nbarrelhits_"+ttype,float(nbarrel));
  //------------------------------------------------------------------- 
}

void PrimaryVertexValidation::fillTrendPlot ( TH1F *  trendPlot, TH1F *  residualsPlot[100], statmode::estimator  fitPar_, TString  var_  )

private

Definition at line 2333 of file PrimaryVertexValidation.cc.

References gather_cfg::cout, etaOfProbe_, etaSect_, fitResiduals(), getMAD(), getMedian(), mps_fire::i, statmode::MAD, statmode::MEAN, statmode::MEDIAN, nBins_, phiSect_, Pi, and statmode::WIDTH.

Referenced by endJob().

{
   
  for ( int i=0; i<nBins_; ++i ) {
    
    char phipositionString[129];
    float phiInterval = phiSect_*(180/TMath::Pi());
    float phiposition = (-180+i*phiInterval)+(phiInterval/2);
    sprintf(phipositionString,"%.f",phiposition);
    
    char etapositionString[129];
    float etaposition = (-etaOfProbe_+i*etaSect_)+(etaSect_/2);
    sprintf(etapositionString,"%.1f",etaposition);
    
    switch(fitPar_)
      {
      case statmode::MEAN:
    {
      float mean_      = fitResiduals(residualsPlot[i]).first.first;
      float meanErr_   = fitResiduals(residualsPlot[i]).first.second;
      trendPlot->SetBinContent(i+1,mean_);
      trendPlot->SetBinError(i+1,meanErr_);
      break;
    } 
      case statmode::WIDTH:
    {
      float width_     = fitResiduals(residualsPlot[i]).second.first;
      float widthErr_  = fitResiduals(residualsPlot[i]).second.second;
      trendPlot->SetBinContent(i+1,width_);
      trendPlot->SetBinError(i+1,widthErr_);
      break;
    }
      case statmode::MEDIAN:
    {
      float median_    = getMedian(residualsPlot[i]).first;
      float medianErr_ = getMedian(residualsPlot[i]).second;
      trendPlot->SetBinContent(i+1,median_);
      trendPlot->SetBinError(i+1,medianErr_);
      break;
    } 
      case statmode::MAD:
    {
      float mad_       = getMAD(residualsPlot[i]).first; 
      float madErr_    = getMAD(residualsPlot[i]).second;
      trendPlot->SetBinContent(i+1,mad_);
      trendPlot->SetBinError(i+1,madErr_);
      break;
    } 
      default:
    edm::LogWarning("PrimaryVertexValidation")<<"fillTrendPlot() "<<fitPar_<<" unknown estimator!"<<std::endl;
    break;
      }

    if(var_=="eta"){
      trendPlot->GetXaxis()->SetBinLabel(i+1,etapositionString); 
    } else if(var_=="phi"){
      trendPlot->GetXaxis()->SetBinLabel(i+1,phipositionString); 
    } else {
      std::cout<<"PrimaryVertexValidation::fillTrendPlot() "<<var_<<" unknown track parameter!"<<std::endl;
    }
  }
}

void PrimaryVertexValidation::fillTrendPlotByIndex ( TH1F *  trendPlot, std::vector< TH1F * > &  h, statmode::estimator  fitPar_  )

private

Definition at line 2398 of file PrimaryVertexValidation.cc.

References stringResolutionProvider_cfi::bin, SoftLeptonByDistance_cfi::distance, fitResiduals(), getMAD(), getMedian(), statmode::MAD, statmode::MEAN, statmode::MEDIAN, and statmode::WIDTH.

Referenced by endJob().

{  

  for(auto iterator = h.begin(); iterator != h.end(); iterator++) {
    
    unsigned int bin = std::distance(h.begin(),iterator)+1;
    std::pair<std::pair<Double_t,Double_t>, std::pair<Double_t,Double_t>  > myFit = fitResiduals((*iterator));

    switch(fitPar_)
      {
      case statmode::MEAN: 
    {   
      float mean_      = myFit.first.first;
      float meanErr_   = myFit.first.second;
      trendPlot->SetBinContent(bin,mean_);
      trendPlot->SetBinError(bin,meanErr_);
      break;
    }
      case statmode::WIDTH:
    {
      float width_     = myFit.second.first;
      float widthErr_  = myFit.second.second;
      trendPlot->SetBinContent(bin,width_);
      trendPlot->SetBinError(bin,widthErr_);
      break;
    }
      case statmode::MEDIAN:
    {
      float median_    = getMedian(*iterator).first;
      float medianErr_ = getMedian(*iterator).second;
      trendPlot->SetBinContent(bin,median_);
      trendPlot->SetBinError(bin,medianErr_);
      break;
    }
      case statmode::MAD:
    {
      float mad_       = getMAD(*iterator).first; 
      float madErr_    = getMAD(*iterator).second;
      trendPlot->SetBinContent(bin,mad_);
      trendPlot->SetBinError(bin,madErr_);
      break;
    }
      default:
    edm::LogWarning("PrimaryVertexValidation")<<"fillTrendPlotByIndex() "<<fitPar_<<" unknown estimator!"<<std::endl;
    break;
      }
  }
}

std::pair< std::pair< Double_t, Double_t >, std::pair< Double_t, Double_t > > PrimaryVertexValidation::fitResiduals ( TH1 *  hist )

private

Definition at line 2290 of file PrimaryVertexValidation.cc.

References RecoEcal_EventContent_cff::func, RecoTauDiscriminantConfiguration::mean, mps_fire::result, and ctppsDiamondLocalTracks_cfi::sigma.

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

  std::pair<Double_t,Double_t> resultM;
  std::pair<Double_t,Double_t> resultW;

  resultM = std::make_pair(res_mean,res_mean_err);
  resultW = std::make_pair(res_width,res_width_err);

  std::pair<std::pair<Double_t,Double_t>, std::pair<Double_t,Double_t>  > result;
  
  result = std::make_pair(resultM,resultW);
  return result;
}

std::pair< Double_t, Double_t > PrimaryVertexValidation::getMAD ( TH1F *  histo )

private

Definition at line 2257 of file PrimaryVertexValidation.cc.

References Abs(), getMedian(), pileupCalc::nbins, mps_fire::result, and create_public_lumi_plots::weights.

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

{

  int nbins = histo->GetNbinsX();
  Double_t median = getMedian(histo).first;
  Double_t x_lastBin = histo->GetBinLowEdge(nbins+1);
  const char *HistoName =histo->GetName();
  TString Finalname = Form("resMed%s",HistoName);
  TH1F *newHisto = new TH1F(Finalname,Finalname,nbins,0.,x_lastBin);
  Double_t *residuals = new Double_t[nbins];
  Double_t *weights = new Double_t[nbins];

  for (int j = 0; j < nbins; j++) {
    residuals[j] = TMath::Abs(median - histo->GetBinCenter(j+1));
    weights[j]=histo->GetBinContent(j+1);
    newHisto->Fill(residuals[j],weights[j]);
  }
  
  Double_t theMAD = (getMedian(newHisto).first)*1.4826;
  
  delete[] residuals; residuals=0;
  delete[] weights; weights=0;
  newHisto->Delete("");
  
  std::pair<Double_t,Double_t> result;
  result = std::make_pair(theMAD,theMAD/histo->GetEntries());

  return result;

}

std::pair< Double_t, Double_t > PrimaryVertexValidation::getMedian ( TH1F *  histo )

private

Definition at line 2231 of file PrimaryVertexValidation.cc.

References pileupCalc::nbins, mps_fire::result, x, and y.

Referenced by fillMap(), fillTrendPlot(), fillTrendPlotByIndex(), and getMAD().

{
  Double_t median = 999;
  int nbins = histo->GetNbinsX();

  //extract median from histogram
  double *x = new double[nbins];
  double *y = new double[nbins];
  for (int j = 0; j < nbins; j++) {
    x[j] = histo->GetBinCenter(j+1);
    y[j] = histo->GetBinContent(j+1);
  }
  median = TMath::Median(nbins, x, y);
  
  delete[] x; x = 0;
  delete[] y; y = 0;  

  std::pair<Double_t,Double_t> result;
  result = std::make_pair(median,median/TMath::Sqrt(histo->GetEntries()));

  return result;

}

bool PrimaryVertexValidation::hasFirstLayerPixelHits ( const reco::TransientTrack  track )

private

Definition at line 976 of file PrimaryVertexValidation.cc.

References reco::HitPattern::getHitPattern(), reco::HitPattern::getLayer(), reco::TransientTrack::hitPattern(), mps_fire::i, reco::HitPattern::numberOfHits(), AlCaHLTBitMon_ParallelJobs::p, listBenchmarks::pattern, reco::HitPattern::pixelBarrelHitFilter(), reco::HitPattern::pixelEndcapHitFilter(), and reco::HitPattern::validHitFilter().

Referenced by analyze().

{
  using namespace reco;
  const HitPattern& p = track.hitPattern();      
  for (int i=0; i<p.numberOfHits(HitPattern::TRACK_HITS); i++) {
    uint32_t pattern = p.getHitPattern(HitPattern::TRACK_HITS, i);   
    if (p.pixelBarrelHitFilter(pattern) || p.pixelEndcapHitFilter(pattern) ) {
      if (p.getLayer(pattern) == 1) {
    if (p.validHitFilter(pattern)) {
      return true;
    }
      }
    }
  }
  return false;
} 

bool PrimaryVertexValidation::isBFieldConsistentWithMode ( const edm::EventSetup &  iSetup ) const

private

Definition at line 2142 of file PrimaryVertexValidation.cc.

References RunInfo::m_avg_current, edm::ESHandle< T >::product(), and ptOfProbe_.

Referenced by analyze().

{
  edm::ESHandle<RunInfo> runInfo;
  iSetup.get<RunInfoRcd>().get(runInfo);
  
  double average_current = runInfo.product()->m_avg_current;
  bool isOn = (average_current > 2000.);
  bool is0T = (ptOfProbe_==0.);
  
  return ( (isOn && !is0T) || (!isOn && is0T) );
}

bool PrimaryVertexValidation::isHit2D ( const TrackingRecHit &  hit ) const

private

Definition at line 927 of file PrimaryVertexValidation.cc.

References TrackingRecHit::dimension(), TrackingRecHit::geographicalId(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, and DetId::Tracker.

Referenced by analyze().

{
  if (hit.dimension() < 2) {
    return false; // some (muon...) stuff really has RecHit1D
  } else {
    const DetId detId(hit.geographicalId());
    if (detId.det() == DetId::Tracker) {
      if (detId.subdetId() == PixelSubdetector::PixelBarrel || detId.subdetId() == PixelSubdetector::PixelEndcap) {
        return true; // pixel is always 2D
      } else { // should be SiStrip now
        if (dynamic_cast<const SiStripRecHit2D*>(&hit)) return false; // normal hit
        else if (dynamic_cast<const SiStripMatchedRecHit2D*>(&hit)) return true; // matched is 2D
        else if (dynamic_cast<const ProjectedSiStripRecHit2D*>(&hit)) return false; // crazy hit...
        else {
          edm::LogError("UnkownType") << "@SUB=AlignmentTrackSelector::isHit2D"
                                      << "Tracker hit not in pixel and neither SiStripRecHit2D nor "
                                      << "SiStripMatchedRecHit2D nor ProjectedSiStripRecHit2D.";
          return false;
        }
      }
    } else { // not tracker??
      edm::LogWarning("DetectorMismatch") << "@SUB=AlignmentTrackSelector::isHit2D"
                                          << "Hit not in tracker with 'official' dimension >=2.";
      return true; // dimension() >= 2 so accept that...
    }
  }
  // never reached...
}

bool PrimaryVertexValidation::passesTrackCuts ( const reco::Track &  track, const reco::Vertex &  vertex, std::string  qualityString_, double  dxyErrMax_, double  dzErrMax_, double  ptErrMax_  )

private

Definition at line 2522 of file PrimaryVertexValidation.cc.

References reco::TrackBase::d0Error(), particleFlowDisplacedVertexCandidate_cfi::dxy, reco::TrackBase::dxy(), allConversions_cfi::dz, reco::TrackBase::dz(), reco::TrackBase::dzError(), reco::Vertex::position(), reco::TrackBase::pt(), reco::TrackBase::ptError(), reco::TrackBase::quality(), reco::TrackBase::qualityByName(), mathSSE::sqrt(), reco::Vertex::xError(), reco::Vertex::yError(), and reco::Vertex::zError().

{
 
   math::XYZPoint vtxPoint(0.0,0.0,0.0);
   double vzErr =0.0, vxErr=0.0, vyErr=0.0;
   vtxPoint=vertex.position();
   vzErr=vertex.zError();
   vxErr=vertex.xError();
   vyErr=vertex.yError();

   double dxy=0.0, dz=0.0, dxysigma=0.0, dzsigma=0.0;
   dxy = track.dxy(vtxPoint);
   dz = track.dz(vtxPoint);
   dxysigma = sqrt(track.d0Error()*track.d0Error()+vxErr*vyErr);
   dzsigma = sqrt(track.dzError()*track.dzError()+vzErr*vzErr);
 
   if(track.quality(reco::TrackBase::qualityByName(qualityString_)) != 1)return false;
   if(fabs(dxy/dxysigma) > dxyErrMax_) return false;
   if(fabs(dz/dzsigma) > dzErrMax_) return false;
   if(track.ptError() / track.pt() > ptErrMax_) return false;

   return true;
}

std::pair< bool, bool > PrimaryVertexValidation::pixelHitsCheck ( const reco::TransientTrack  track )

private

Definition at line 958 of file PrimaryVertexValidation.cc.

References reco::TransientTrack::hitPattern(), reco::HitPattern::numberOfValidPixelBarrelHits(), reco::HitPattern::numberOfValidPixelEndcapHits(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by analyze().

                                                                                        {
  
  bool hasBPixHits = false;
  bool hasFPixHits = false;
  
  const reco::HitPattern& p = track.hitPattern();
  if(p.numberOfValidPixelEndcapHits()!=0){
    hasFPixHits = true;
  } 
  if(p.numberOfValidPixelBarrelHits()!=0){
    hasBPixHits = true;
  }
  
  return std::make_pair(hasBPixHits,hasFPixHits);
}

void PrimaryVertexValidation::SetVarToZero ( )

private

Definition at line 2156 of file PrimaryVertexValidation.cc.

References Beamdxdz_, Beamsigmaz_, BeamWidthX_, BeamWidthY_, BSx0_, BSy0_, BSz0_, charge_, chi2_, chi2ndof_, chi2normUnbiasedVertex_, chi2ProbUnbiasedVertex_, chi2UnbiasedVertex_, d3DErrorFromMyVertex_, d3DFromMyVertex_, DOFUnbiasedVertex_, dxy_, dxyBs_, dxyErrorFromMyVertex_, dxyFromMyVertex_, dz_, dzBs_, dzErrorFromMyVertex_, dzFromMyVertex_, eta_, hasRecVertex_, mps_fire::i, IP3DsigFromMyVertex_, IPLsigFromMyVertex_, IPTsigFromMyVertex_, isGoodTrack_, isHighPurity_, LuminosityBlockNumber_, nClus_, nhits1D_, nhits2D_, nhits_, nhitsBPIX_, nhitsFPIX_, nhitsTEC_, nhitsTIB_, nhitsTID_, nhitsTOB_, nMaxtracks_, nOfflineVertices_, nTracks_, p_, phi_, pt_, qoverp_, RunNumber_, sumOfWeightsUnbiasedVertex_, theta_, tracksUsedForVertexing_, wxy2_, xErrOfflineVertex_, xOfflineVertex_, xPCA_, xUnbiasedVertex_, yErrOfflineVertex_, yOfflineVertex_, yPCA_, yUnbiasedVertex_, zErrOfflineVertex_, zOfflineVertex_, zPCA_, and zUnbiasedVertex_.

Referenced by analyze().

{
  nTracks_ = 0;
  nClus_ = 0;
  nOfflineVertices_=0;
  RunNumber_ =0;
  LuminosityBlockNumber_=0;
  xOfflineVertex_ =-999.;
  yOfflineVertex_ =-999.;
  zOfflineVertex_ =-999.;
  xErrOfflineVertex_=0.;
  yErrOfflineVertex_=0.;
  zErrOfflineVertex_=0.;
  BSx0_ = -999.;
  BSy0_ = -999.;
  BSz0_ = -999.;
  Beamsigmaz_=-999.;
  Beamdxdz_=-999.;   
  BeamWidthX_=-999.;
  BeamWidthY_=-999.;
  wxy2_=-999.;

  for ( int i=0; i<nMaxtracks_; ++i ) {
    
    pt_[i]        = 0;
    p_[i]         = 0;
    nhits_[i]     = 0;
    nhits1D_[i]   = 0;
    nhits2D_[i]   = 0;
    nhitsBPIX_[i] = 0;
    nhitsFPIX_[i] = 0;
    nhitsTIB_[i]  = 0;
    nhitsTID_[i]  = 0;
    nhitsTOB_[i]  = 0;
    nhitsTEC_[i]  = 0;
    isHighPurity_[i] = 0;
    eta_[i]       = 0;
    theta_[i]     = 0;
    phi_[i]       = 0;
    chi2_[i]      = 0;
    chi2ndof_[i]  = 0;
    charge_[i]    = 0;
    qoverp_[i]    = 0;
    dz_[i]        = 0;
    dxy_[i]       = 0;
    dzBs_[i]      = 0;
    dxyBs_[i]     = 0;
    xPCA_[i]      = 0;
    yPCA_[i]      = 0;
    zPCA_[i]      = 0;
    xUnbiasedVertex_[i] =0;    
    yUnbiasedVertex_[i] =0;
    zUnbiasedVertex_[i] =0;
    chi2normUnbiasedVertex_[i]=0;
    chi2UnbiasedVertex_[i]=0;
    chi2ProbUnbiasedVertex_[i]=0;
    DOFUnbiasedVertex_[i]=0;
    sumOfWeightsUnbiasedVertex_[i]=0;
    tracksUsedForVertexing_[i]=0;
    dxyFromMyVertex_[i]=0;
    dzFromMyVertex_[i]=0;
    d3DFromMyVertex_[i]=0;
    dxyErrorFromMyVertex_[i]=0; 
    dzErrorFromMyVertex_[i]=0;
    d3DErrorFromMyVertex_[i]=0;
    IPTsigFromMyVertex_[i]=0;   
    IPLsigFromMyVertex_[i]=0;  
    IP3DsigFromMyVertex_[i]=0;
    hasRecVertex_[i] = 0;
    isGoodTrack_[i]  = 0;
  } 
}

double PrimaryVertexValidation::square ( double  x )

inlineprivate

Definition at line 107 of file PrimaryVertexValidation.h.

                                {
    return x*x;
  }

bool PrimaryVertexValidation::vtxSort ( const reco::Vertex &  a, const reco::Vertex &  b  )

staticprivate

Definition at line 2512 of file PrimaryVertexValidation.cc.

References reco::Vertex::chi2(), funct::false, and reco::Vertex::tracksSize().

Referenced by analyze().

{
  if( a.tracksSize() != b.tracksSize() )
    return a.tracksSize() > b.tracksSize() ? true : false ;
  else
    return a.chi2() < b.chi2() ? true : false ;
}

Member Data Documentation

TH1F* PrimaryVertexValidation::a_d3DEtaResiduals[nMaxBins_]

private

Definition at line 274 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_d3DPhiResiduals[nMaxBins_]

private

Definition at line 273 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_d3DResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 300 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_dxEtaResiduals[nMaxBins_]

private

Definition at line 256 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_dxPhiResiduals[nMaxBins_]

private

Definition at line 255 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_dxyBiasResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 420 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaBiasResiduals[nMaxBins_]

private

Definition at line 405 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaMADBiasTrend

private

Definition at line 456 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaMADTrend

private

Definition at line 358 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaMeanBiasTrend

private

Definition at line 433 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaMeanTrend

private

Definition at line 313 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaMedianBiasTrend

private

Definition at line 455 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaMedianTrend

private

Definition at line 357 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaResiduals[nMaxBins_]

private

Definition at line 253 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaWidthBiasTrend

private

Definition at line 434 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyEtaWidthTrend

private

Definition at line 314 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::a_dxyMeanBiasMap

private

Definition at line 472 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::a_dxyMeanMap

private

Definition at line 388 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiBiasResiduals[nMaxBins_]

private

Definition at line 404 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiMADBiasTrend

private

Definition at line 451 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiMADTrend

private

Definition at line 353 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiMeanBiasTrend

private

Definition at line 428 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiMeanTrend

private

Definition at line 308 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiMedianBiasTrend

private

Definition at line 450 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiMedianTrend

private

Definition at line 352 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiResiduals[nMaxBins_]

private

Definition at line 252 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiWidthBiasTrend

private

Definition at line 429 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyPhiWidthTrend

private

Definition at line 309 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxypTCentralMeanTrend

private

Definition at line 335 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxypTCentralWidthTrend

private

Definition at line 336 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxypTMeanTrend

private

Definition at line 330 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxypTWidthTrend

private

Definition at line 331 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dxyResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 298 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* PrimaryVertexValidation::a_dxyVsEta

private

Definition at line 380 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::a_dxyVsPhi

private

Definition at line 374 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::a_dxyWidthBiasMap

private

Definition at line 478 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::a_dxyWidthMap

private

Definition at line 394 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dyEtaResiduals[nMaxBins_]

private

Definition at line 259 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_dyPhiResiduals[nMaxBins_]

private

Definition at line 258 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_dzBiasResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 421 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaBiasResiduals[nMaxBins_]

private

Definition at line 408 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaMADBiasTrend

private

Definition at line 458 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaMADTrend

private

Definition at line 360 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaMeanBiasTrend

private

Definition at line 435 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaMeanTrend

private

Definition at line 315 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaMedianBiasTrend

private

Definition at line 457 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaMedianTrend

private

Definition at line 359 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaResiduals[nMaxBins_]

private

Definition at line 262 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaWidthBiasTrend

private

Definition at line 436 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzEtaWidthTrend

private

Definition at line 316 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::a_dzMeanBiasMap

private

Definition at line 473 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::a_dzMeanMap

private

Definition at line 389 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiBiasResiduals[nMaxBins_]

private

Definition at line 407 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiMADBiasTrend

private

Definition at line 453 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiMADTrend

private

Definition at line 355 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiMeanBiasTrend

private

Definition at line 430 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiMeanTrend

private

Definition at line 310 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiMedianBiasTrend

private

Definition at line 452 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiMedianTrend

private

Definition at line 354 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiResiduals[nMaxBins_]

private

Definition at line 261 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiWidthBiasTrend

private

Definition at line 431 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzPhiWidthTrend

private

Definition at line 311 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzpTCentralMeanTrend

private

Definition at line 337 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzpTCentralWidthTrend

private

Definition at line 338 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzpTMeanTrend

private

Definition at line 332 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzpTWidthTrend

private

Definition at line 333 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_dzResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 299 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* PrimaryVertexValidation::a_dzVsEta

private

Definition at line 381 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::a_dzVsPhi

private

Definition at line 375 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::a_dzWidthBiasMap

private

Definition at line 479 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::a_dzWidthMap

private

Definition at line 395 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::a_IP2DEtaResiduals[nMaxBins_]

private

Definition at line 265 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_IP2DPhiResiduals[nMaxBins_]

private

Definition at line 264 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_IP3DEtaResiduals[nMaxBins_]

private

Definition at line 268 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_IP3DPhiResiduals[nMaxBins_]

private

Definition at line 267 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_reszEtaResiduals[nMaxBins_]

private

Definition at line 271 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::a_reszPhiResiduals[nMaxBins_]

private

Definition at line 270 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

bool PrimaryVertexValidation::askFirstLayerHit_

private

Definition at line 133 of file PrimaryVertexValidation.h.

Referenced by analyze().

double PrimaryVertexValidation::Beamdxdz_

private

Definition at line 187 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::Beamsigmaz_

private

Definition at line 186 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::BeamWidthX_

private

Definition at line 188 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::BeamWidthY_

private

Definition at line 189 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::BSx0_

private

Definition at line 183 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::BSy0_

private

Definition at line 184 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::BSz0_

private

Definition at line 185 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::charge_[nMaxtracks_]

private

Definition at line 210 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::chi2_[nMaxtracks_]

private

Definition at line 208 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::chi2ndof_[nMaxtracks_]

private

Definition at line 209 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

float PrimaryVertexValidation::chi2normUnbiasedVertex_[nMaxtracks_]

private

Definition at line 222 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

float PrimaryVertexValidation::chi2ProbUnbiasedVertex_[nMaxtracks_]

private

Definition at line 224 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

float PrimaryVertexValidation::chi2UnbiasedVertex_[nMaxtracks_]

private

Definition at line 223 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

const int PrimaryVertexValidation::cmToum = 10000

staticprivate

Definition at line 156 of file PrimaryVertexValidation.h.

Referenced by analyze().

double PrimaryVertexValidation::d3DErrorFromMyVertex_[nMaxtracks_]

private

Definition at line 235 of file PrimaryVertexValidation.h.

Referenced by analyze(), and SetVarToZero().

double PrimaryVertexValidation::d3DFromMyVertex_[nMaxtracks_]

private

Definition at line 231 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

bool PrimaryVertexValidation::debug_

private

Definition at line 142 of file PrimaryVertexValidation.h.

Referenced by analyze().

bool PrimaryVertexValidation::doBPix_

private

Definition at line 134 of file PrimaryVertexValidation.h.

Referenced by analyze().

bool PrimaryVertexValidation::doFPix_

private

Definition at line 135 of file PrimaryVertexValidation.h.

Referenced by analyze().

float PrimaryVertexValidation::DOFUnbiasedVertex_[nMaxtracks_]

private

Definition at line 225 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::dxy_[nMaxtracks_]

private

Definition at line 213 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::dxyBs_[nMaxtracks_]

private

Definition at line 214 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::dxyErrorFromMyVertex_[nMaxtracks_]

private

Definition at line 233 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::dxyFromMyVertex_[nMaxtracks_]

private

Definition at line 229 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::dz_[nMaxtracks_]

private

Definition at line 212 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::dzBs_[nMaxtracks_]

private

Definition at line 215 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::dzErrorFromMyVertex_[nMaxtracks_]

private

Definition at line 234 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::dzFromMyVertex_[nMaxtracks_]

private

Definition at line 230 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::eta_[nMaxtracks_]

private

Definition at line 205 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::etaOfProbe_

private

Definition at line 137 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), fillMap(), fillTrendPlot(), and PrimaryVertexValidation().

float PrimaryVertexValidation::etaSect_

private

Definition at line 160 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), fillMap(), fillTrendPlot(), and PrimaryVertexValidation().

unsigned int PrimaryVertexValidation::EventNumber_

private

Definition at line 173 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::h2_probeEtaPhi_

private

Definition at line 507 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::h2_probeEtaPt_

private

Definition at line 508 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_Beamsigmaz

private

Definition at line 501 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_BeamWidthX

private

Definition at line 502 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_BeamWidthY

private

Definition at line 503 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_BSx0

private

Definition at line 498 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_BSy0

private

Definition at line 499 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_BSz0

private

Definition at line 500 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

std::vector<TH1F*> PrimaryVertexValidation::h_dxy_Central_pT_

private

Definition at line 574 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> PrimaryVertexValidation::h_dxy_pT_

private

Definition at line 569 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> PrimaryVertexValidation::h_dz_Central_pT_

private

Definition at line 575 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> PrimaryVertexValidation::h_dz_pT_

private

Definition at line 570 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::h_etaMax

private

Definition at line 245 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_fitVtxChi2_

private

Definition at line 553 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_fitVtxChi2ndf_

private

Definition at line 555 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_fitVtxChi2Prob_

private

Definition at line 556 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_fitVtxNdof_

private

Definition at line 552 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_fitVtxNtracks_

private

Definition at line 554 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_fitVtxTrackAverageWeight_

private

Definition at line 558 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_fitVtxTrackWeights_

private

Definition at line 557 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_lumiFromConfig

private

Definition at line 489 of file PrimaryVertexValidation.h.

Referenced by beginJob().

TH1F* PrimaryVertexValidation::h_nbins

private

Definition at line 246 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_nClus

private

Definition at line 486 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_nOfflineVertices

private

Definition at line 487 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

std::vector<TH1F*> PrimaryVertexValidation::h_norm_dxy_Central_pT_

private

Definition at line 576 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> PrimaryVertexValidation::h_norm_dxy_pT_

private

Definition at line 571 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> PrimaryVertexValidation::h_norm_dz_Central_pT_

private

Definition at line 577 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> PrimaryVertexValidation::h_norm_dz_pT_

private

Definition at line 572 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::h_nTracks

private

Definition at line 485 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeCharge_

private

Definition at line 516 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeChi2_

private

Definition at line 514 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probed0RefitV_

private

Definition at line 525 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probed3DRefitV_

private

Definition at line 529 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probedxyRecoV_

private

Definition at line 520 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probedxyRefitV_

private

Definition at line 523 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probedzRecoV_

private

Definition at line 519 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probedzRefitV_

private

Definition at line 522 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeEta_

private

Definition at line 512 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHits1D_

private

Definition at line 541 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHits2D_

private

Definition at line 542 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHits_

private

Definition at line 540 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHitsInBPIX_

private

Definition at line 547 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHitsInFPIX_

private

Definition at line 548 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHitsInTEC_

private

Definition at line 546 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHitsInTIB_

private

Definition at line 543 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHitsInTID_

private

Definition at line 545 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeHitsInTOB_

private

Definition at line 544 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeNormChi2_

private

Definition at line 515 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeP_

private

Definition at line 510 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probePhi_

private

Definition at line 513 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probePt_

private

Definition at line 511 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeQoverP_

private

Definition at line 517 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeRecoVSigXY_

private

Definition at line 533 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeRecoVSigZ_

private

Definition at line 532 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeRefitVLogSig3D_

private

Definition at line 537 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeRefitVSig3D_

private

Definition at line 536 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeRefitVSigResZ_

private

Definition at line 538 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeRefitVSigXY_

private

Definition at line 535 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probeRefitVSigZ_

private

Definition at line 534 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probereszRefitV_

private

Definition at line 530 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probesignIP2DRefitV_

private

Definition at line 528 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_probez0RefitV_

private

Definition at line 526 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_recoVtxChi2ndf_

private

Definition at line 561 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_recoVtxChi2Prob_

private

Definition at line 562 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_recoVtxNtracks_

private

Definition at line 560 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_recoVtxSumPt_

private

Definition at line 563 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1I* PrimaryVertexValidation::h_runFromConfig

private

Definition at line 490 of file PrimaryVertexValidation.h.

Referenced by beginJob().

TH1I* PrimaryVertexValidation::h_runFromEvent

private

Definition at line 491 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_runNumber

private

Definition at line 488 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_xErrOfflineVertex

private

Definition at line 495 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_xOfflineVertex

private

Definition at line 492 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_yErrOfflineVertex

private

Definition at line 496 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_yOfflineVertex

private

Definition at line 493 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_zErrOfflineVertex

private

Definition at line 497 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::h_zOfflineVertex

private

Definition at line 494 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

int PrimaryVertexValidation::hasRecVertex_[nMaxtracks_]

private

Definition at line 241 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

std::map<std::string, TH1*> PrimaryVertexValidation::hDA

private

Definition at line 565 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

double PrimaryVertexValidation::intLumi_

private

Definition at line 130 of file PrimaryVertexValidation.h.

Referenced by beginJob().

double PrimaryVertexValidation::IP3DsigFromMyVertex_[nMaxtracks_]

private

Definition at line 239 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::IPLsigFromMyVertex_[nMaxtracks_]

private

Definition at line 238 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::IPTsigFromMyVertex_[nMaxtracks_]

private

Definition at line 237 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::isGoodTrack_[nMaxtracks_]

private

Definition at line 242 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::isHighPurity_[nMaxtracks_]

private

Definition at line 204 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

bool PrimaryVertexValidation::isPhase1_

private

Definition at line 138 of file PrimaryVertexValidation.h.

Referenced by analyze().

bool PrimaryVertexValidation::lightNtupleSwitch_

private

Definition at line 123 of file PrimaryVertexValidation.h.

Referenced by beginJob().

unsigned int PrimaryVertexValidation::LuminosityBlockNumber_

private

Definition at line 174 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

const float PrimaryVertexValidation::mypT_bins_[nPtBins_+1] = {0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3.0,3.1,3.2,3.3,3.4,3.5,3.6,3.7,3.8,3.9,4.0,4.1,4.25,4.5,4.75,5.0,5.5,6.0,7.0,8.0,9.0,11.0,14.0,20.}

private

Definition at line 165 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_d3DEtaResiduals[nMaxBins_]

private

Definition at line 294 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_d3DPhiResiduals[nMaxBins_]

private

Definition at line 293 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_d3DResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 304 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_dxyBiasResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 423 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaBiasResiduals[nMaxBins_]

private

Definition at line 413 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaMADBiasTrend

private

Definition at line 466 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaMADTrend

private

Definition at line 368 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaMeanBiasTrend

private

Definition at line 443 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaMeanTrend

private

Definition at line 323 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaMedianBiasTrend

private

Definition at line 465 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaMedianTrend

private

Definition at line 367 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaResiduals[nMaxBins_]

private

Definition at line 279 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaWidthBiasTrend

private

Definition at line 444 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyEtaWidthTrend

private

Definition at line 324 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::n_dxyMeanBiasMap

private

Definition at line 475 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::n_dxyMeanMap

private

Definition at line 391 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiBiasResiduals[nMaxBins_]

private

Definition at line 412 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiMADBiasTrend

private

Definition at line 461 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiMADTrend

private

Definition at line 363 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiMeanBiasTrend

private

Definition at line 438 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiMeanTrend

private

Definition at line 318 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiMedianBiasTrend

private

Definition at line 460 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiMedianTrend

private

Definition at line 362 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiResiduals[nMaxBins_]

private

Definition at line 278 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiWidthBiasTrend

private

Definition at line 439 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyPhiWidthTrend

private

Definition at line 319 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxypTCentralMeanTrend

private

Definition at line 345 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxypTCentralWidthTrend

private

Definition at line 346 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxypTMeanTrend

private

Definition at line 340 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxypTWidthTrend

private

Definition at line 341 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dxyResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 302 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* PrimaryVertexValidation::n_dxyVsEta

private

Definition at line 383 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::n_dxyVsPhi

private

Definition at line 377 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::n_dxyWidthBiasMap

private

Definition at line 481 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::n_dxyWidthMap

private

Definition at line 397 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzBiasResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 424 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaBiasResiduals[nMaxBins_]

private

Definition at line 416 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaMADBiasTrend

private

Definition at line 468 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaMADTrend

private

Definition at line 370 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaMeanBiasTrend

private

Definition at line 445 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaMeanTrend

private

Definition at line 325 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaMedianBiasTrend

private

Definition at line 467 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaMedianTrend

private

Definition at line 369 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaResiduals[nMaxBins_]

private

Definition at line 282 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaWidthBiasTrend

private

Definition at line 446 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzEtaWidthTrend

private

Definition at line 326 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::n_dzMeanBiasMap

private

Definition at line 476 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::n_dzMeanMap

private

Definition at line 392 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiBiasResiduals[nMaxBins_]

private

Definition at line 415 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiMADBiasTrend

private

Definition at line 463 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiMADTrend

private

Definition at line 365 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiMeanBiasTrend

private

Definition at line 440 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiMeanTrend

private

Definition at line 320 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiMedianBiasTrend

private

Definition at line 462 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiMedianTrend

private

Definition at line 364 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiResiduals[nMaxBins_]

private

Definition at line 281 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiWidthBiasTrend

private

Definition at line 441 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzPhiWidthTrend

private

Definition at line 321 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzpTCentralMeanTrend

private

Definition at line 347 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzpTCentralWidthTrend

private

Definition at line 348 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzpTMeanTrend

private

Definition at line 342 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzpTWidthTrend

private

Definition at line 343 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_dzResidualsMap[nMaxBins_][nMaxBins_]

private

Definition at line 303 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

TH2F* PrimaryVertexValidation::n_dzVsEta

private

Definition at line 384 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::n_dzVsPhi

private

Definition at line 378 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH2F* PrimaryVertexValidation::n_dzWidthBiasMap

private

Definition at line 482 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH2F* PrimaryVertexValidation::n_dzWidthMap

private

Definition at line 398 of file PrimaryVertexValidation.h.

Referenced by beginJob(), and endJob().

TH1F* PrimaryVertexValidation::n_IP2DEtaResiduals[nMaxBins_]

private

Definition at line 285 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_IP2DPhiResiduals[nMaxBins_]

private

Definition at line 284 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_IP3DEtaResiduals[nMaxBins_]

private

Definition at line 288 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_IP3DPhiResiduals[nMaxBins_]

private

Definition at line 287 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_reszEtaResiduals[nMaxBins_]

private

Definition at line 291 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

TH1F* PrimaryVertexValidation::n_reszPhiResiduals[nMaxBins_]

private

Definition at line 290 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

int PrimaryVertexValidation::nBins_

private

Definition at line 139 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), fillMap(), fillTrendPlot(), and PrimaryVertexValidation().

int PrimaryVertexValidation::nClus_

private

Definition at line 170 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::Nevt_

private

Definition at line 113 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

int PrimaryVertexValidation::nhits1D_[nMaxtracks_]

private

Definition at line 196 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nhits2D_[nMaxtracks_]

private

Definition at line 197 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nhits_[nMaxtracks_]

private

Definition at line 195 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nhitsBPIX_[nMaxtracks_]

private

Definition at line 198 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nhitsFPIX_[nMaxtracks_]

private

Definition at line 199 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nhitsTEC_[nMaxtracks_]

private

Definition at line 203 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nhitsTIB_[nMaxtracks_]

private

Definition at line 200 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nhitsTID_[nMaxtracks_]

private

Definition at line 201 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nhitsTOB_[nMaxtracks_]

private

Definition at line 202 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

const int PrimaryVertexValidation::nMaxBins_ = 100

staticprivate

Definition at line 119 of file PrimaryVertexValidation.h.

const int PrimaryVertexValidation::nMaxtracks_ = 1000

staticprivate

Definition at line 155 of file PrimaryVertexValidation.h.

Referenced by analyze(), and SetVarToZero().

int PrimaryVertexValidation::nOfflineVertices_

private

Definition at line 171 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

const int PrimaryVertexValidation::nPtBins_ = 48

staticprivate

Definition at line 157 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

int PrimaryVertexValidation::nTracks_

private

Definition at line 168 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

int PrimaryVertexValidation::nTracksPerClus_

private

Definition at line 169 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

double PrimaryVertexValidation::p_[nMaxtracks_]

private

Definition at line 194 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::phi_[nMaxtracks_]

private

Definition at line 207 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

float PrimaryVertexValidation::phiSect_

private

Definition at line 159 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), fillMap(), fillTrendPlot(), and PrimaryVertexValidation().

double PrimaryVertexValidation::pt_[nMaxtracks_]

private

Definition at line 193 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::ptOfProbe_

private

Definition at line 136 of file PrimaryVertexValidation.h.

Referenced by analyze(), and isBFieldConsistentWithMode().

double PrimaryVertexValidation::qoverp_[nMaxtracks_]

private

Definition at line 211 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

TTree* PrimaryVertexValidation::rootTree_

private

Definition at line 149 of file PrimaryVertexValidation.h.

Referenced by analyze(), and beginJob().

bool PrimaryVertexValidation::runControl_

private

Definition at line 143 of file PrimaryVertexValidation.h.

Referenced by analyze().

std::vector<unsigned int> PrimaryVertexValidation::runControlNumbers_

private

Definition at line 140 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and PrimaryVertexValidation().

unsigned int PrimaryVertexValidation::RunNumber_

private

Definition at line 172 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

bool PrimaryVertexValidation::storeNtuple_

private

Definition at line 122 of file PrimaryVertexValidation.h.

Referenced by analyze().

float PrimaryVertexValidation::sumOfWeightsUnbiasedVertex_[nMaxtracks_]

private

Definition at line 226 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

edm::EDGetTokenT<reco::BeamSpot> PrimaryVertexValidation::theBeamspotToken

private

Definition at line 147 of file PrimaryVertexValidation.h.

Referenced by analyze(), and PrimaryVertexValidation().

edm::ParameterSet PrimaryVertexValidation::theConfig

private

Definition at line 112 of file PrimaryVertexValidation.h.

double PrimaryVertexValidation::theta_[nMaxtracks_]

private

Definition at line 206 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

TrackClusterizerInZ* PrimaryVertexValidation::theTrackClusterizer_

private

Definition at line 116 of file PrimaryVertexValidation.h.

Referenced by analyze(), and PrimaryVertexValidation().

edm::EDGetTokenT<reco::TrackCollection> PrimaryVertexValidation::theTrackCollectionToken

private

Definition at line 145 of file PrimaryVertexValidation.h.

Referenced by analyze(), and PrimaryVertexValidation().

TrackFilterForPVFindingBase* PrimaryVertexValidation::theTrackFilter_

private

Definition at line 115 of file PrimaryVertexValidation.h.

Referenced by analyze(), and PrimaryVertexValidation().

edm::EDGetTokenT<reco::VertexCollection> PrimaryVertexValidation::theVertexCollectionToken

private

Definition at line 146 of file PrimaryVertexValidation.h.

Referenced by analyze(), and PrimaryVertexValidation().

int PrimaryVertexValidation::tracksUsedForVertexing_[nMaxtracks_]

private

Definition at line 227 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

bool PrimaryVertexValidation::useTracksFromRecoVtx_

private

Definition at line 124 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and endJob().

double PrimaryVertexValidation::vertexZMax_

private

Definition at line 127 of file PrimaryVertexValidation.h.

Referenced by analyze().

double PrimaryVertexValidation::wxy2_

private

Definition at line 190 of file PrimaryVertexValidation.h.

Referenced by analyze(), fillTrackHistos(), and SetVarToZero().

double PrimaryVertexValidation::xErrOfflineVertex_

private

Definition at line 179 of file PrimaryVertexValidation.h.

Referenced by analyze(), and SetVarToZero().

double PrimaryVertexValidation::xOfflineVertex_

private

Definition at line 175 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::xPCA_[nMaxtracks_]

private

Definition at line 216 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::xUnbiasedVertex_[nMaxtracks_]

private

Definition at line 219 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::yErrOfflineVertex_

private

Definition at line 180 of file PrimaryVertexValidation.h.

Referenced by analyze(), and SetVarToZero().

double PrimaryVertexValidation::yOfflineVertex_

private

Definition at line 176 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::yPCA_[nMaxtracks_]

private

Definition at line 217 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::yUnbiasedVertex_[nMaxtracks_]

private

Definition at line 220 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::zErrOfflineVertex_

private

Definition at line 181 of file PrimaryVertexValidation.h.

Referenced by analyze(), and SetVarToZero().

double PrimaryVertexValidation::zOfflineVertex_

private

Definition at line 177 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::zPCA_[nMaxtracks_]

private

Definition at line 218 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().

double PrimaryVertexValidation::zUnbiasedVertex_[nMaxtracks_]

private

Definition at line 221 of file PrimaryVertexValidation.h.

Referenced by analyze(), beginJob(), and SetVarToZero().