TkConvValidator Class Reference

#include <TkConvValidator.h>

Inheritance diagram for TkConvValidator:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &r, edm::EventSetup const &theEventSetup) override
void	dqmEndRun (edm::Run const &r, edm::EventSetup const &es) override
void	endJob () override
	TkConvValidator (const edm::ParameterSet &)
	~TkConvValidator () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
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
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
float	etaTransformation (float a, float b)
float	phiNormalization (float &a)
math::XYZVector	recalculateMomentumAtFittedVertex (const MagneticField &mf, const TrackerGeometry &trackerGeom, const edm::RefToBase< reco::Track > &tk, const reco::Vertex &vtx)

Private Attributes
bool	arbitratedEcalSeeded_
bool	arbitratedMerged_
double	bcEtLow_
edm::EDGetTokenT< reco::BeamSpot >	beamspotToken_
std::string	conversionCollection_
edm::EDGetTokenT< reco::ConversionCollection >	conversionCollectionPr_Token_
std::string	conversionCollectionProducer_
std::string	conversionTrackProducer_
DQMStore *	dbe_
bool	dCotCutOn_
double	dCotCutValue_
double	dCotHardCutValue_
std::string	dqmpath_
bool	ecalalgotracks_
double	ecalEtSumCut_
double	ecalIsolRadius_
std::string	fName_
edm::EDGetTokenT< edm::SimTrackContainer >	g4_simTk_Token_
edm::EDGetTokenT< edm::SimVertexContainer >	g4_simVtx_Token_
bool	generalTracksOnly_
edm::EDGetTokenT< reco::GenJetCollection >	genjets_Token_
MonitorElement *	h2_Chi2VsEta_ [3]
MonitorElement *	h2_Chi2VsR_ [3]
MonitorElement *	h2_convVtxdRVsEta_
MonitorElement *	h2_convVtxdRVsR_
MonitorElement *	h2_convVtxRrecVsTrue_
MonitorElement *	h2_DCotTracksVsEta_
MonitorElement *	h2_DCotTracksVsR_
MonitorElement *	h2_DPhiTracksAtEcalVsEta_
MonitorElement *	h2_DPhiTracksAtEcalVsR_
MonitorElement *	h2_DPhiTracksAtVtxVsEta_
MonitorElement *	h2_DPhiTracksAtVtxVsR_
MonitorElement *	h2_dzPVVsR_
MonitorElement *	h2_photonPtRecVsPtSim_
MonitorElement *	h2_PtRecVsPtSim_ [3]
MonitorElement *	h2_TkPtPull_ [3]
MonitorElement *	h_AllSimConv_ [5]
	Denominator for efficiencies. More...
MonitorElement *	h_convEta2_ [3][3]
MonitorElement *	h_convEta_ [3][3]
MonitorElement *	h_convEtaMatchSC_ [3][3]
MonitorElement *	h_convPhi_ [3][3]
MonitorElement *	h_convPt_ [3][3]
MonitorElement *	h_convPtRes_ [3]
MonitorElement *	h_convR_ [3][3]
MonitorElement *	h_convRplot_
MonitorElement *	h_convSCdEta_ [3][3]
MonitorElement *	h_convSCdPhi_ [3][3]
MonitorElement *	h_convVtxdEta_
MonitorElement *	h_convVtxdPhi_
MonitorElement *	h_convVtxdR_
MonitorElement *	h_convVtxdR_barrel_
MonitorElement *	h_convVtxdR_endcap_
MonitorElement *	h_convVtxdX_
MonitorElement *	h_convVtxdX_barrel_
MonitorElement *	h_convVtxdX_endcap_
MonitorElement *	h_convVtxdY_
MonitorElement *	h_convVtxdY_barrel_
MonitorElement *	h_convVtxdY_endcap_
MonitorElement *	h_convVtxdZ_
MonitorElement *	h_convVtxdZ_barrel_
MonitorElement *	h_convVtxdZ_endcap_
MonitorElement *	h_convVtxRvsZ_ [3]
MonitorElement *	h_convVtxRvsZ_zoom_ [2]
MonitorElement *	h_convVtxYvsX_
MonitorElement *	h_convVtxYvsX_zoom_ [2]
MonitorElement *	h_convZ_ [3][3]
MonitorElement *	h_convZplot_
MonitorElement *	h_DCotTracks_ [3][3]
MonitorElement *	h_deltaExpectedHitsInner_ [3][3]
MonitorElement *	h_DEtaTracksAtEcal_ [3][3]
MonitorElement *	h_distMinAppTracks_ [3][3]
MonitorElement *	h_dlClosestHitToVtx_ [3]
MonitorElement *	h_dlClosestHitToVtxSig_ [3]
MonitorElement *	h_DPhiTracksAtEcal_ [3][3]
MonitorElement *	h_DPhiTracksAtVtx_ [3][3]
MonitorElement *	h_dzPVFromTracks_ [2]
MonitorElement *	h_EoverPTracks_ [3][3]
MonitorElement *	h_invMass_ [3][3]
MonitorElement *	h_leadExpectedHitsInner_ [3][3]
MonitorElement *	h_leadNHitsBeforeVtx_ [3][3]
MonitorElement *	h_lxybs_ [3][3]
MonitorElement *	h_match_
MonitorElement *	h_maxDlClosestHitToVtx_ [3][3]
MonitorElement *	h_maxDlClosestHitToVtxSig_ [3][3]
MonitorElement *	h_maxNHitsBeforeVtx_ [3][3]
MonitorElement *	h_nConv_ [3][3]
	info per conversion More...
MonitorElement *	h_nHitsBeforeVtx_ [3]
MonitorElement *	h_nSharedHits_ [3][3]
MonitorElement *	h_nSimConv_ [2]
MonitorElement *	h_RecoConvTwoMTracks_ [5]
MonitorElement *	h_RecoConvTwoTracks_ [5]
MonitorElement *	h_SimConvEtaPix_ [2]
MonitorElement *	h_SimConvOneMTracks_ [5]
MonitorElement *	h_SimConvOneTracks_ [5]
	Numerator for efficiencies. More...
MonitorElement *	h_SimConvTwoMTracks_ [5]
MonitorElement *	h_SimConvTwoMTracksAndVtxPGT0005_ [5]
MonitorElement *	h_SimConvTwoMTracksAndVtxPGT01_ [5]
MonitorElement *	h_SimConvTwoMTracksAndVtxPGT0_ [5]
MonitorElement *	h_SimConvTwoTracks_ [5]
MonitorElement *	h_simConvVtxRvsZ_ [4]
MonitorElement *	h_simConvVtxYvsX_
MonitorElement *	h_SimRecConvOneMTracks_ [5]
MonitorElement *	h_SimRecConvOneTracks_ [5]
MonitorElement *	h_SimRecConvTwoMTracks_ [5]
MonitorElement *	h_SimRecConvTwoTracks_ [5]
MonitorElement *	h_simTkEta_
MonitorElement *	h_simTkPt_
MonitorElement *	h_sumNHitsBeforeVtx_ [3][3]
MonitorElement *	h_tkChi2_ [3]
MonitorElement *	h_tkChi2Large_ [3]
MonitorElement *	h_TkD0_ [3]
MonitorElement *	h_TkPtPull_ [3]
MonitorElement *	h_trailNHitsBeforeVtx_ [3][3]
MonitorElement *	h_VisSimConv_ [6]
MonitorElement *	h_VisSimConvLarge_
MonitorElement *	h_vtxChi2Prob_ [3][3]
MonitorElement *	h_zPVFromTracks_ [2]
double	hcalEtSumCut_
double	hcalHitEtLow_
double	hcalIsolExtRadius_
double	hcalIsolInnRadius_
edm::EDGetTokenT< edm::HepMCProduct >	hepMC_Token_
bool	highPurity_
bool	isRunCentrally_
edm::InputTag	label_tp_
double	lip_
uint	maxHitsBeforeVtx_
double	maxPhoEtaForEffic
double	maxPhoEtaForPurity
double	maxPhoRForEffic
double	maxPhoRForPurity
double	maxPhoZForEffic
double	maxPhoZForPurity
double	mcConvEta_
double	mcConvPhi_
double	mcConvPt_
double	mcConvR_
double	mcConvX_
double	mcConvY_
double	mcConvZ_
double	mcEta_
double	mcJetEta_
double	mcJetPhi_
double	mcPhi_
	global variable for the MC photon More...
double	minLxy_
double	minPhoEtCut_
double	minPhoPtForEffic
double	minPhoPtForPurity
double	minProb_
int	nEntry_
int	nEvt_
MonitorElement *	nHits_ [3]
MonitorElement *	nHitsVsEta_ [3]
MonitorElement *	nHitsVsR_ [3]
int	nInvalidPCA_
int	nMatched_
int	nRecConv_
int	nRecConvAss_
int	nRecConvAssWithEcal_
int	nSimConv_ [2]
int	numOfTracksInCone_
edm::EDGetTokenT< reco::VertexCollection >	offline_pvToken_
MonitorElement *	p2_convVtxdRVsRZ_
MonitorElement *	p2_convVtxdZVsRZ_
MonitorElement *	p2_effRZ_
MonitorElement *	p_Chi2VsEta_ [3]
MonitorElement *	p_Chi2VsR_ [3]
MonitorElement *	p_convVtxdRVsEta_
MonitorElement *	p_convVtxdRVsR_
MonitorElement *	p_convVtxdXVsX_
MonitorElement *	p_convVtxdYVsY_
MonitorElement *	p_convVtxdZVsR_
MonitorElement *	p_convVtxdZVsZ_
MonitorElement *	p_DCotTracksVsEta_
MonitorElement *	p_DCotTracksVsR_
MonitorElement *	p_DPhiTracksAtEcalVsEta_
MonitorElement *	p_DPhiTracksAtEcalVsR_
MonitorElement *	p_DPhiTracksAtVtxVsEta_
MonitorElement *	p_DPhiTracksAtVtxVsR_
MonitorElement *	p_dzPVVsR_
MonitorElement *	p_nHitsVsEta_ [3]
MonitorElement *	p_nHitsVsR_ [3]
MonitorElement *	p_TkPtPull_ [3]
edm::ParameterSet	parameters_
std::string	photonCollection_
edm::EDGetTokenT< reco::PhotonCollection >	photonCollectionPr_Token_
std::string	photonCollectionProducer_
double	recMaxPt_
double	recMinPt_
	Global variables for reco Photon. More...
double	simMaxPt_
double	simMinPt_
edm::ESHandle< CaloGeometry >	theCaloGeom_
edm::ESHandle< CaloTopology >	theCaloTopo_
edm::RefVector< TrackingParticleCollection >	theConvTP_
edm::ESHandle< MagneticField >	theMF_
PhotonMCTruthFinder *	thePhotonMCTruthFinder_
edm::EDGetTokenT< TrackingParticleRefVector >	tpSelForEff_Token_
edm::EDGetTokenT< TrackingParticleRefVector >	tpSelForFake_Token_
edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator >	trackAssociator_Token_
double	trkIsolExtRadius_
double	trkIsolInnRadius_
double	trkPtLow_
double	trkPtSumCut_
int	verbosity_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

$Id: TkConvValidator

Author

N.Marinelli - Univ. of Notre Dame

Definition at line 49 of file TkConvValidator.h.

Constructor & Destructor Documentation

TkConvValidator::TkConvValidator ( const edm::ParameterSet &  pset )

explicit

Definition at line 101 of file TkConvValidator.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_2018_cff::InputTag, muonDTDigis_cfi::pset, AlCaHLTBitMon_QueryRunRegistry::string, and parallelization::uint.

                                                            {
  fName_ = pset.getUntrackedParameter<std::string>("Name");
  verbosity_ = pset.getUntrackedParameter<int>("Verbosity");
  parameters_ = pset;

  photonCollectionProducer_ = pset.getParameter<std::string>("phoProducer");
  photonCollection_ = pset.getParameter<std::string>("photonCollection");
  photonCollectionPr_Token_ =
      consumes<reco::PhotonCollection>(edm::InputTag(photonCollectionProducer_, photonCollection_));

  conversionCollectionProducer_ = pset.getParameter<std::string>("convProducer");
  conversionCollection_ = pset.getParameter<std::string>("conversionCollection");
  conversionCollectionPr_Token_ =
      consumes<reco::ConversionCollection>(edm::InputTag(conversionCollectionProducer_, conversionCollection_));

  // conversionTrackProducer_ = pset.getParameter<std::string>("trackProducer");
  dqmpath_ = pset.getParameter<std::string>("dqmpath");
  minPhoEtCut_ = pset.getParameter<double>("minPhoEtCut");
  generalTracksOnly_ = pset.getParameter<bool>("generalTracksOnly");
  arbitratedMerged_ = pset.getParameter<bool>("arbitratedMerged");
  arbitratedEcalSeeded_ = pset.getParameter<bool>("arbitratedEcalSeeded");
  ecalalgotracks_ = pset.getParameter<bool>("ecalalgotracks");
  highPurity_ = pset.getParameter<bool>("highPurity");
  minProb_ = pset.getParameter<double>("minProb");
  maxHitsBeforeVtx_ = pset.getParameter<uint>("maxHitsBeforeVtx");
  minLxy_ = pset.getParameter<double>("minLxy");
  isRunCentrally_ = pset.getParameter<bool>("isRunCentrally");

  offline_pvToken_ = consumes<reco::VertexCollection>(
      pset.getUntrackedParameter<edm::InputTag>("offlinePV", edm::InputTag("offlinePrimaryVertices")));
  beamspotToken_ =
      consumes<reco::BeamSpot>(pset.getUntrackedParameter<edm::InputTag>("beamspot", edm::InputTag("offlineBeamSpot")));
  g4_simTk_Token_ = consumes<edm::SimTrackContainer>(pset.getParameter<edm::InputTag>("simTracks"));
  g4_simVtx_Token_ = consumes<edm::SimVertexContainer>(pset.getParameter<edm::InputTag>("simTracks"));

  tpSelForEff_Token_ = consumes<TrackingParticleRefVector>(edm::InputTag("tpSelecForEfficiency"));
  tpSelForFake_Token_ = consumes<TrackingParticleRefVector>(edm::InputTag("tpSelecForFakeRate"));
  //hepMC_Token_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));
  //genjets_Token_ = consumes<reco::GenJetCollection>(
  //    edm::InputTag("ak4GenJets"));

  trackAssociator_Token_ =
      consumes<reco::TrackToTrackingParticleAssociator>(edm::InputTag("trackAssociatorByHitsForConversionValidation"));
}

TkConvValidator::~TkConvValidator ( )

override

Definition at line 146 of file TkConvValidator.cc.

{}

Member Function Documentation

void TkConvValidator::analyze ( const edm::Event &  e, const edm::EventSetup &  esup  )

override

Definition at line 1241 of file TkConvValidator.cc.

References funct::abs(), reco::TrackBase::algo(), reco::Conversion::arbitratedEcalSeeded, reco::Conversion::arbitratedMerged, reco::TrackToTrackingParticleAssociator::associateRecoToSim(), reco::TrackToTrackingParticleAssociator::associateSimToReco(), ALCARECOTkAlJpsiMuMu_cff::charge, reco::Vertex::chi2(), ChiSquaredProbability(), conv, reco::Conversion::conversionVertex(), HLTMuonOfflineAnalyzer_cfi::d0, reco::deltaPhi(), reco::Conversion::distOfMinimumApproach(), reco::Conversion::dlClosestHitToVtx(), reco::Conversion::dPhiTracksAtVtx(), HGC3DClusterGenMatchSelector_cfi::dR, edm::AssociationMap< Tag >::empty(), edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), reco::Conversion::generalTracksOnly, edm::EventSetup::get(), edm::RefToBase< T >::get(), edm::Event::getByToken(), reco::Conversion::highPurity, reco::TrackBase::hitPattern(), mps_fire::i, edm::EventBase::id(), reco::TrackBase::inOutEcalSeededConv, edm::HandleBase::isValid(), reco::Vertex::isValid(), edm::Ref< C, T, F >::key(), match(), SiStripPI::max, tkConvValidator_cfi::maxPhoEtaForEffic, tkConvValidator_cfi::maxPhoEtaForPurity, tkConvValidator_cfi::maxPhoRForEffic, tkConvValidator_cfi::maxPhoZForEffic, tkConvValidator_cfi::maxPhoZForPurity, tkConvValidator_cfi::minPhoPtForEffic, tkConvValidator_cfi::minPhoPtForPurity, reco::HitPattern::MISSING_INNER_HITS, reco::Vertex::ndof(), reco::Conversion::nHitsBeforeVtx(), TrackCollections2monitor_cff::normalizedChi2, reco::Conversion::nSharedHits(), reco::HitPattern::numberOfLostHits(), TrackCollections2monitor_cff::numberOfValidHits, reco::TrackBase::outInEcalSeededConv, p1, p2, reco::Conversion::pairCotThetaSeparation(), reco::Conversion::pairInvariantMass(), EgammaValidation_cff::pdgId, nanoDQM_cff::Photon, ExoticaDQM_cfi::photonCollection, reco::Vertex::position(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), reco::TrackBase::pt(), edm::RefToBaseVector< T >::push_back(), q1, q2, reco::Conversion::quality(), recPt, reco::Conversion::refittedPairMomentum(), reco::Vertex::refittedTracks(), mathSSE::sqrt(), reco::Photon::superCluster(), PbPb_ZMuSkimMuonDPG_cff::tracker, PDWG_EXOHSCP_cff::tracks, reco::Conversion::tracks(), parallelization::uint, photonValidationSequence_cff::useTP, spclusmultinvestigator_cfi::vertexCollection, badGlobalMuonTaggersAOD_cff::vtx, reco::Vertex::x(), reco::BeamSpot::x0(), reco::Vertex::y(), reco::BeamSpot::y0(), and reco::Conversion::zOfPrimaryVertexFromTracks().

                                                                          {
  thePhotonMCTruthFinder_->clear();
  using namespace edm;
  //  const float etaPhiDistance=0.01;
  // Fiducial region
  // const float TRK_BARL =0.9;
  const float BARL = 1.4442;  // DAQ TDR p.290
  //  const float END_LO = 1.566; // unused
  const float END_HI = 2.5;
  // Electron mass
  //  const Float_t mElec= 0.000511; // unused

  Handle<reco::TrackToTrackingParticleAssociator> theTrackAssociator;
  e.getByToken(trackAssociator_Token_, theTrackAssociator);

  nEvt_++;
  LogInfo("TkConvValidator") << "TkConvValidator Analyzing event number: " << e.id() << " Global Counter " << nEvt_
                             << "\n";
  //  std::cout << "TkConvValidator Analyzing event number: "  << e.id() << " Global Counter " << nEvt_ <<"\n";

  // get the geometry from the event setup:
  esup.get<CaloGeometryRecord>().get(theCaloGeom_);

  // Transform Track into TransientTrack (needed by the Vertex fitter)
  edm::ESHandle<TransientTrackBuilder> theTTB;
  esup.get<TransientTrackRecord>().get("TransientTrackBuilder", theTTB);

  Handle<reco::ConversionCollection> convHandle;
  e.getByToken(conversionCollectionPr_Token_, convHandle);
  const reco::ConversionCollection convCollection = *(convHandle.product());
  if (!convHandle.isValid()) {
    edm::LogError("ConversionsProducer") << "Error! Can't get the  collection " << std::endl;
    return;
  }

  Handle<reco::PhotonCollection> photonHandle;
  e.getByToken(photonCollectionPr_Token_, photonHandle);
  const reco::PhotonCollection photonCollection = *(photonHandle.product());
  if (!photonHandle.isValid()) {
    edm::LogError("PhotonProducer") << "Error! Can't get the Photon collection " << std::endl;
    return;
  }

  // offline  Primary vertex
  edm::Handle<reco::VertexCollection> vertexHandle;
  reco::VertexCollection vertexCollection;
  e.getByToken(offline_pvToken_, vertexHandle);
  if (!vertexHandle.isValid()) {
    edm::LogError("TrackerOnlyConversionProducer") << "Error! Can't get the product primary Vertex Collection "
                                                   << "\n";
  } else {
    vertexCollection = *(vertexHandle.product());
  }
  reco::Vertex the_pvtx;
  bool valid_pvtx = false;
  if (!vertexCollection.empty()) {
    the_pvtx = *(vertexCollection.begin());
    //asking for one good vertex
    if (the_pvtx.isValid() && fabs(the_pvtx.position().z()) <= 15 && the_pvtx.position().Rho() <= 2) {
      valid_pvtx = true;
    }
  }

  edm::Handle<reco::BeamSpot> bsHandle;
  e.getByToken(beamspotToken_, bsHandle);
  if (!bsHandle.isValid()) {
    edm::LogError("TrackerOnlyConversionProducer") << "Error! Can't get the product primary Vertex Collection "
                                                   << "\n";
    return;
  }
  const reco::BeamSpot& thebs = *bsHandle.product();

  //get tracker geometry for hits positions
  edm::ESHandle<TrackerGeometry> tracker;
  esup.get<TrackerDigiGeometryRecord>().get(tracker);
  const TrackerGeometry* trackerGeom = tracker.product();

  //get simtrack info
  //std::vector<SimTrack> theSimTracks;
  //std::vector<SimVertex> theSimVertices;

  edm::Handle<SimTrackContainer> SimTk;
  edm::Handle<SimVertexContainer> SimVtx;
  e.getByToken(g4_simTk_Token_, SimTk);
  e.getByToken(g4_simVtx_Token_, SimVtx);

  bool useTP = parameters_.getParameter<bool>("useTP");
  TrackingParticleRefVector dummy;
  edm::Handle<TrackingParticleRefVector> TPHandleForEff;
  edm::Handle<TrackingParticleRefVector> TPHandleForFakeRate;
  if (useTP) {
    e.getByToken(tpSelForEff_Token_, TPHandleForEff);
    e.getByToken(tpSelForFake_Token_, TPHandleForFakeRate);
  }

  const TrackingParticleRefVector& tpForEfficiency = useTP ? *(TPHandleForEff.product()) : dummy;
  const TrackingParticleRefVector& tpForFakeRate = useTP ? *(TPHandleForFakeRate.product()) : dummy;

  const std::vector<SimTrack>& theSimTracks = *SimTk;
  const std::vector<SimVertex>& theSimVertices = *SimVtx;

  //theSimTracks.insert(theSimTracks.end(),SimTk->begin(),SimTk->end());
  //theSimVertices.insert(theSimVertices.end(),SimVtx->begin(),SimVtx->end());
  std::vector<PhotonMCTruth> mcPhotons = thePhotonMCTruthFinder_->find(theSimTracks, theSimVertices);

  //edm::Handle<edm::HepMCProduct> hepMC;
  //e.getByToken(hepMC_Token_, hepMC);
  //  const HepMC::GenEvent *myGenEvent = hepMC->GetEvent(); // unused

  // get generated jets
  //edm::Handle<reco::GenJetCollection> GenJetsHandle;
  //e.getByToken(genjets_Token_, GenJetsHandle);
  //const reco::GenJetCollection &genJetCollection = *(GenJetsHandle.product());

  ConversionHitChecker hitChecker;

  // ################  SIM to RECO ######################### //
  std::map<const reco::Track*, TrackingParticleRef> myAss;
  std::map<const reco::Track*, TrackingParticleRef>::const_iterator itAss;

  for (std::vector<PhotonMCTruth>::const_iterator mcPho = mcPhotons.begin(); mcPho != mcPhotons.end(); mcPho++) {
    mcConvPt_ = (*mcPho).fourMomentum().et();
    float mcPhi = (*mcPho).fourMomentum().phi();
    mcPhi_ = phiNormalization(mcPhi);
    mcEta_ = (*mcPho).fourMomentum().pseudoRapidity();
    mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z());
    mcConvR_ = (*mcPho).vertex().perp();
    mcConvX_ = (*mcPho).vertex().x();
    mcConvY_ = (*mcPho).vertex().y();
    mcConvZ_ = (*mcPho).vertex().z();
    mcConvEta_ = (*mcPho).vertex().eta();
    mcConvPhi_ = (*mcPho).vertex().phi();

    if (fabs(mcEta_) > END_HI)
      continue;

    if (mcConvPt_ < minPhoPtForEffic)
      continue;
    if (fabs(mcEta_) > maxPhoEtaForEffic)
      continue;
    if (fabs(mcConvZ_) > maxPhoZForEffic)
      continue;
    if (mcConvR_ > maxPhoRForEffic)
      continue;

    bool goodSimConversion = false;
    bool visibleConversion = false;
    bool visibleConversionsWithTwoSimTracks = false;
    if ((*mcPho).isAConversion() == 1) {
      nSimConv_[0]++;
      h_AllSimConv_[0]->Fill(mcEta_);
      h_AllSimConv_[1]->Fill(mcPhi_);
      h_AllSimConv_[2]->Fill(mcConvR_);
      h_AllSimConv_[3]->Fill(mcConvZ_);
      h_AllSimConv_[4]->Fill((*mcPho).fourMomentum().et());

      if (mcConvR_ < 15)
        h_SimConvEtaPix_[0]->Fill(mcEta_);

      if ((fabs(mcEta_) <= BARL && mcConvR_ < 85) ||
          (fabs(mcEta_) > BARL && fabs(mcEta_) <= END_HI && fabs((*mcPho).vertex().z()) < 210))
        visibleConversion = true;

      theConvTP_.clear();
      //      std::cout << " TkConvValidator TrackingParticles   TrackingParticleCollection size "<<  trackingParticles.size() <<  "\n";
      //duplicated TP collections for two associations
      for (const TrackingParticleRef tp : tpForEfficiency) {
        if (fabs(tp->vx() - (*mcPho).vertex().x()) < 0.0001 && fabs(tp->vy() - (*mcPho).vertex().y()) < 0.0001 &&
            fabs(tp->vz() - (*mcPho).vertex().z()) < 0.0001) {
          theConvTP_.push_back(tp);
        }
      }
      //std::cout << " TkConvValidator  theConvTP_ size " <<   theConvTP_.size() << std::endl;

      if (theConvTP_.size() == 2)
        visibleConversionsWithTwoSimTracks = true;
      goodSimConversion = false;

      if (visibleConversion && visibleConversionsWithTwoSimTracks)
        goodSimConversion = true;
      if (goodSimConversion) {
        nSimConv_[1]++;
        h_VisSimConv_[0]->Fill(mcEta_);
        h_VisSimConv_[1]->Fill(mcPhi_);
        h_VisSimConv_[2]->Fill(mcConvR_);
        h_VisSimConv_[3]->Fill(mcConvZ_);
        h_VisSimConv_[4]->Fill((*mcPho).fourMomentum().et());
      }

      for (edm::RefVector<TrackingParticleCollection>::iterator iTrk = theConvTP_.begin(); iTrk != theConvTP_.end();
           ++iTrk) {
        h_simTkPt_->Fill((*iTrk)->pt());
        h_simTkEta_->Fill((*iTrk)->eta());
      }

    }  

    if (!(visibleConversion && visibleConversionsWithTwoSimTracks))
      continue;

    h_simConvVtxRvsZ_[0]->Fill(fabs(mcConvZ_), mcConvR_);
    if (fabs(mcEta_) <= 1.) {
      h_simConvVtxRvsZ_[1]->Fill(fabs(mcConvZ_), mcConvR_);
      h_simConvVtxYvsX_->Fill(mcConvX_, mcConvY_);
    } else
      h_simConvVtxRvsZ_[2]->Fill(fabs(mcConvZ_), mcConvR_);

    //std::cout << " TkConvValidator  theConvTP_ size " <<   theConvTP_.size() << std::endl;
    for (edm::RefVector<TrackingParticleCollection>::iterator iTP = theConvTP_.begin(); iTP != theConvTP_.end();
         iTP++) {
      //  std::cout << " SIM to RECO TP vertex " << (*iTP)->vx() << " " <<  (*iTP)->vy() << " " << (*iTP)->vz() << " pt " << (*iTP)->pt() << std::endl;
    }

    bool recomatch = false;
    float chi2Prob = 0.;
    // cout << " size of conversions " << convHandle->size() << endl;
    for (reco::ConversionCollection::const_iterator conv = convHandle->begin(); conv != convHandle->end(); ++conv) {
      const reco::Conversion aConv = (*conv);
      if (arbitratedMerged_ && !aConv.quality(reco::Conversion::arbitratedMerged))
        continue;
      if (generalTracksOnly_ && !aConv.quality(reco::Conversion::generalTracksOnly))
        continue;
      if (arbitratedEcalSeeded_ && !aConv.quality(reco::Conversion::arbitratedEcalSeeded))
        continue;

      if (highPurity_ && !aConv.quality(reco::Conversion::highPurity))
        continue;

      //problematic?
      const std::vector<edm::RefToBase<reco::Track> >& tracks = aConv.tracks();

      const reco::Vertex& vtx = aConv.conversionVertex();
      //requires two tracks and a valid vertex
      if (tracks.size() != 2 || !(vtx.isValid()))
        continue;

      if (ChiSquaredProbability(aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof()) <= minProb_)
        continue;
      if (aConv.nHitsBeforeVtx().size() > 1 &&
          max(aConv.nHitsBeforeVtx().at(0), aConv.nHitsBeforeVtx().at(1)) > maxHitsBeforeVtx_)
        continue;

      //compute transverse decay length with respect to beamspot
      math::XYZVectorF themom = aConv.refittedPairMomentum();
      double dbsx = aConv.conversionVertex().x() - thebs.x0();
      double dbsy = aConv.conversionVertex().y() - thebs.y0();
      double lxy = (themom.x() * dbsx + themom.y() * dbsy) / themom.rho();

      if (lxy < minLxy_)
        continue;

      //    bool  phoIsInBarrel=false; // unused
      //    bool  phoIsInEndcap=false; // unused
      RefToBase<reco::Track> tfrb1 = aConv.tracks().front();
      RefToBase<reco::Track> tfrb2 = aConv.tracks().back();

      if (ecalalgotracks_ && (!(tfrb1->algo() == reco::TrackBase::outInEcalSeededConv ||
                                tfrb1->algo() == reco::TrackBase::inOutEcalSeededConv) ||
                              !(tfrb2->algo() == reco::TrackBase::outInEcalSeededConv ||
                                tfrb2->algo() == reco::TrackBase::inOutEcalSeededConv)))
        continue;

      //reco::TrackRef tk1 = aConv.tracks().front();
      //reco::TrackRef tk2 = aConv.tracks().back();
      //std::cout << "SIM to RECO  conversion track pt " << tk1->pt() << " " << tk2->pt() << endl;
      //
      //Use two RefToBaseVector and do two association actions to avoid that if two tracks from different collection
      RefToBaseVector<reco::Track> tc1, tc2;
      tc1.push_back(tfrb1);
      tc2.push_back(tfrb2);
      bool isAssociated = false;
      reco::SimToRecoCollection q1 = theTrackAssociator->associateSimToReco(tc1, theConvTP_);
      reco::SimToRecoCollection q2 = theTrackAssociator->associateSimToReco(tc2, theConvTP_);
      //try {
      std::vector<std::pair<RefToBase<reco::Track>, double> > trackV1, trackV2;

      int tp_1 = 0, tp_2 = 1;  //the index of associated tp in theConvTP_ for two tracks
      if (q1.find(theConvTP_[0]) != q1.end()) {
        trackV1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >)q1[theConvTP_[0]];
      } else if (q1.find(theConvTP_[1]) != q1.end()) {
        trackV1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >)q1[theConvTP_[1]];
        tp_1 = 1;
      }
      if (q2.find(theConvTP_[1]) != q2.end()) {
        trackV2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >)q2[theConvTP_[1]];
      } else if (q2.find(theConvTP_[0]) != q2.end()) {
        trackV2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >)q2[theConvTP_[0]];
        tp_2 = 0;
      }
      if (!(!trackV1.empty() && !trackV2.empty()))
        continue;
      if (tp_1 == tp_2)
        continue;

      edm::RefToBase<reco::Track> tr1 = trackV1.front().first;
      edm::RefToBase<reco::Track> tr2 = trackV2.front().first;
      //std::cout << "associated tp1 " <<theConvTP_[0]->pt() << " to track with  pT=" << tr1->pt() << " " << (tr1.get())->pt() << endl;
      //std::cout << "associated tp2 " <<theConvTP_[1]->pt() << " to track with  pT=" << tr2->pt() << " " << (tr2.get())->pt() << endl;
      myAss.insert(std::make_pair(tr1.get(), theConvTP_[tp_1]));
      myAss.insert(std::make_pair(tr2.get(), theConvTP_[tp_2]));

      //} catch (Exception const& event) {
      //cout << "continue: " << event.what()  << endl;
      //  continue;
      //}

      isAssociated = true;
      recomatch = true;
      chi2Prob = ChiSquaredProbability(aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof());

      if (isAssociated) {
        h_SimRecConvTwoMTracks_[0]->Fill(mcEta_);
        h_SimRecConvTwoMTracks_[1]->Fill(mcPhi_);
        h_SimRecConvTwoMTracks_[2]->Fill(mcConvR_);
        h_SimRecConvTwoMTracks_[3]->Fill(mcConvZ_);
        h_SimRecConvTwoMTracks_[4]->Fill((*mcPho).fourMomentum().et());
      }

      // break;
    }  // loop over reco conversions
    if (recomatch) {
      h_SimConvTwoMTracks_[0]->Fill(mcEta_);
      h_SimConvTwoMTracks_[1]->Fill(mcPhi_);
      h_SimConvTwoMTracks_[2]->Fill(mcConvR_);
      h_SimConvTwoMTracks_[3]->Fill(mcConvZ_);
      h_SimConvTwoMTracks_[4]->Fill((*mcPho).fourMomentum().et());

      if (chi2Prob > 0) {
        h_SimConvTwoMTracksAndVtxPGT0_[0]->Fill(mcEta_);
        h_SimConvTwoMTracksAndVtxPGT0_[1]->Fill(mcPhi_);
        h_SimConvTwoMTracksAndVtxPGT0_[2]->Fill(mcConvR_);
        h_SimConvTwoMTracksAndVtxPGT0_[3]->Fill(mcConvZ_);
        h_SimConvTwoMTracksAndVtxPGT0_[4]->Fill((*mcPho).fourMomentum().et());
      }
      if (chi2Prob > 0.0005) {
        h_SimConvTwoMTracksAndVtxPGT0005_[0]->Fill(mcEta_);
        h_SimConvTwoMTracksAndVtxPGT0005_[1]->Fill(mcPhi_);
        h_SimConvTwoMTracksAndVtxPGT0005_[2]->Fill(mcConvR_);
        h_SimConvTwoMTracksAndVtxPGT0005_[3]->Fill(mcConvZ_);
        h_SimConvTwoMTracksAndVtxPGT0005_[4]->Fill((*mcPho).fourMomentum().et());
      }
    }

  }  //End loop over simulated conversions

  // ########################### RECO to SIM ############################## //

  for (reco::ConversionCollection::const_iterator conv = convHandle->begin(); conv != convHandle->end(); ++conv) {
    const reco::Conversion aConv = (*conv);
    if (arbitratedMerged_ && !aConv.quality(reco::Conversion::arbitratedMerged))
      continue;
    if (generalTracksOnly_ && !aConv.quality(reco::Conversion::generalTracksOnly))
      continue;
    if (arbitratedEcalSeeded_ && !aConv.quality(reco::Conversion::arbitratedEcalSeeded))
      continue;

    if (highPurity_ && !aConv.quality(reco::Conversion::highPurity))
      continue;

    //problematic?
    std::vector<edm::RefToBase<reco::Track> > tracks = aConv.tracks();

    const reco::Vertex& vtx = aConv.conversionVertex();
    //requires two tracks and a valid vertex
    if (tracks.size() != 2 || !(vtx.isValid()))
      continue;
    //if (tracks.size() !=2) continue;

    if (ChiSquaredProbability(aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof()) <= minProb_)
      continue;
    if (aConv.nHitsBeforeVtx().size() > 1 &&
        max(aConv.nHitsBeforeVtx().at(0), aConv.nHitsBeforeVtx().at(1)) > maxHitsBeforeVtx_)
      continue;

    //compute transverse decay length with respect to beamspot
    math::XYZVectorF themom = aConv.refittedPairMomentum();
    double dbsx = aConv.conversionVertex().x() - thebs.x0();
    double dbsy = aConv.conversionVertex().y() - thebs.y0();
    double lxy = (themom.x() * dbsx + themom.y() * dbsy) / themom.rho();

    if (lxy < minLxy_)
      continue;

    bool phoIsInBarrel = false;
    bool phoIsInEndcap = false;
    RefToBase<reco::Track> tk1 = aConv.tracks().front();
    RefToBase<reco::Track> tk2 = aConv.tracks().back();
    RefToBaseVector<reco::Track> tc1, tc2;
    tc1.push_back(tk1);
    tc2.push_back(tk2);

    if (ecalalgotracks_ &&
        (!(tk1->algo() == reco::TrackBase::outInEcalSeededConv ||
           tk1->algo() == reco::TrackBase::inOutEcalSeededConv) ||
         !(tk2->algo() == reco::TrackBase::outInEcalSeededConv || tk2->algo() == reco::TrackBase::inOutEcalSeededConv)))
      continue;

    //std::cout << " RECO to SIM conversion track pt " << tk1->pt() << " " << tk2->pt() << endl;
    const reco::Track refTk1 = aConv.conversionVertex().refittedTracks().front();
    const reco::Track refTk2 = aConv.conversionVertex().refittedTracks().back();

    //TODO replace it with phi at vertex
    float dPhiTracksAtVtx = aConv.dPhiTracksAtVtx();
    // override with the phi calculated at the vertex
    math::XYZVector p1AtVtx = recalculateMomentumAtFittedVertex((*theMF_), *trackerGeom, tk1, aConv.conversionVertex());
    math::XYZVector p2AtVtx = recalculateMomentumAtFittedVertex((*theMF_), *trackerGeom, tk2, aConv.conversionVertex());
    if (sqrt(p1AtVtx.perp2()) > sqrt(p2AtVtx.perp2()))
      dPhiTracksAtVtx = p1AtVtx.phi() - p2AtVtx.phi();
    else
      dPhiTracksAtVtx = p2AtVtx.phi() - p1AtVtx.phi();

    math::XYZVectorF refittedMom = aConv.refittedPairMomentum();

    if (fabs(refittedMom.eta()) < 1.479) {
      phoIsInBarrel = true;
    } else {
      phoIsInEndcap = true;
    }

    nRecConv_++;

    // check matching with reco photon
    double Mindeltaeta = 999999;
    double Mindeltaphi = 999999;
    bool matchConvSC = false;
    reco::PhotonCollection::const_iterator iMatchingSC;
    for (reco::PhotonCollection::const_iterator iPho = photonCollection.begin(); iPho != photonCollection.end();
         iPho++) {
      reco::Photon aPho = reco::Photon(*iPho);
      const double deltaphi = reco::deltaPhi(aConv.refittedPairMomentum().phi(), aPho.superCluster()->position().phi());
      double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
      double deltaeta = abs(aPho.superCluster()->position().eta() - ConvEta);
      if (abs(deltaeta) < abs(Mindeltaeta) && abs(deltaphi) < abs(Mindeltaphi)) {
        Mindeltaphi = abs(deltaphi);
        Mindeltaeta = abs(deltaeta);
        iMatchingSC = iPho;
      }
    }
    if (abs(Mindeltaeta) < 0.1 && abs(Mindeltaphi) < 0.1) {
      matchConvSC = true;
    }

    int match = 0;
    float invM = aConv.pairInvariantMass();
    float chi2Prob = ChiSquaredProbability(aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof());
    uint maxNHitsBeforeVtx =
        aConv.nHitsBeforeVtx().size() > 1 ? max(aConv.nHitsBeforeVtx().at(0), aConv.nHitsBeforeVtx().at(1)) : 0;
    uint sumNHitsBeforeVtx =
        aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(0) + aConv.nHitsBeforeVtx().at(1) : 0;
    float maxDlClosestHitToVtx = aConv.dlClosestHitToVtx().size() > 1 ? max(aConv.dlClosestHitToVtx().at(0).value(),
                                                                            aConv.dlClosestHitToVtx().at(1).value())
                                                                      : 0;
    float maxDlClosestHitToVtxSig =
        aConv.dlClosestHitToVtx().size() > 1
            ? max(aConv.dlClosestHitToVtx().at(0).value() / aConv.dlClosestHitToVtx().at(0).error(),
                  aConv.dlClosestHitToVtx().at(1).value() / aConv.dlClosestHitToVtx().at(1).error())
            : 0;

    int ilead = 0, itrail = 1;
    if (tk2->pt() > tk1->pt()) {
      ilead = 1;
      itrail = 0;
    }
    RefToBase<reco::Track> tklead = aConv.tracks().at(ilead);
    RefToBase<reco::Track> tktrail = aConv.tracks().at(itrail);

    int deltaExpectedHitsInner = tklead->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS) -
                                 tktrail->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
    int leadExpectedHitsInner = tklead->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
    uint leadNHitsBeforeVtx = aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(ilead) : 0;
    uint trailNHitsBeforeVtx = aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(itrail) : 0;

    h_convEta_[match][0]->Fill(refittedMom.eta());
    h_convEta2_[match][0]->Fill(refittedMom.eta());

    h_convPhi_[match][0]->Fill(refittedMom.phi());
    h_convR_[match][0]->Fill(sqrt(aConv.conversionVertex().position().perp2()));
    h_convRplot_->Fill(sqrt(aConv.conversionVertex().position().perp2()));
    h_convZ_[match][0]->Fill(aConv.conversionVertex().position().z());
    h_convZplot_->Fill(aConv.conversionVertex().position().z());
    h_convPt_[match][0]->Fill(sqrt(refittedMom.perp2()));
    h_invMass_[match][0]->Fill(invM);
    h_vtxChi2Prob_[match][0]->Fill(chi2Prob);
    h_lxybs_[match][0]->Fill(lxy);
    h_maxNHitsBeforeVtx_[match][0]->Fill(maxNHitsBeforeVtx);
    h_leadNHitsBeforeVtx_[match][0]->Fill(leadNHitsBeforeVtx);
    h_trailNHitsBeforeVtx_[match][0]->Fill(trailNHitsBeforeVtx);
    h_sumNHitsBeforeVtx_[match][0]->Fill(sumNHitsBeforeVtx);
    h_deltaExpectedHitsInner_[match][0]->Fill(deltaExpectedHitsInner);
    h_leadExpectedHitsInner_[match][0]->Fill(leadExpectedHitsInner);
    h_maxDlClosestHitToVtx_[match][0]->Fill(maxDlClosestHitToVtx);
    h_maxDlClosestHitToVtxSig_[match][0]->Fill(maxDlClosestHitToVtxSig);
    h_nSharedHits_[match][0]->Fill(aConv.nSharedHits());

    if (matchConvSC) {
      h_convEtaMatchSC_[match][0]->Fill(refittedMom.eta());
      h_EoverPTracks_[match][0]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
      h_convSCdPhi_[match][0]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
      double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
      h_convSCdEta_[match][0]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
    }

    h_distMinAppTracks_[match][0]->Fill(aConv.distOfMinimumApproach());
    h_DPhiTracksAtVtx_[match][0]->Fill(dPhiTracksAtVtx);
    h2_DPhiTracksAtVtxVsEta_->Fill(mcEta_, dPhiTracksAtVtx);
    h2_DPhiTracksAtVtxVsR_->Fill(mcConvR_, dPhiTracksAtVtx);
    p_DPhiTracksAtVtxVsEta_->Fill(mcEta_, dPhiTracksAtVtx);
    p_DPhiTracksAtVtxVsR_->Fill(mcConvR_, dPhiTracksAtVtx);

    h_DCotTracks_[match][0]->Fill(aConv.pairCotThetaSeparation());
    h2_DCotTracksVsEta_->Fill(mcEta_, aConv.pairCotThetaSeparation());
    h2_DCotTracksVsR_->Fill(mcConvR_, aConv.pairCotThetaSeparation());
    p_DCotTracksVsEta_->Fill(mcEta_, aConv.pairCotThetaSeparation());
    p_DCotTracksVsR_->Fill(mcConvR_, aConv.pairCotThetaSeparation());

    if (phoIsInBarrel) {
      h_invMass_[match][1]->Fill(invM);
      h_vtxChi2Prob_[match][1]->Fill(chi2Prob);
      h_distMinAppTracks_[match][1]->Fill(aConv.distOfMinimumApproach());
      h_DPhiTracksAtVtx_[match][1]->Fill(dPhiTracksAtVtx);
      h_DCotTracks_[match][1]->Fill(aConv.pairCotThetaSeparation());
      h_lxybs_[match][1]->Fill(lxy);
      h_maxNHitsBeforeVtx_[match][1]->Fill(maxNHitsBeforeVtx);
      h_leadNHitsBeforeVtx_[match][1]->Fill(leadNHitsBeforeVtx);
      h_trailNHitsBeforeVtx_[match][1]->Fill(trailNHitsBeforeVtx);
      h_sumNHitsBeforeVtx_[match][1]->Fill(sumNHitsBeforeVtx);
      h_deltaExpectedHitsInner_[match][1]->Fill(deltaExpectedHitsInner);
      h_leadExpectedHitsInner_[match][1]->Fill(leadExpectedHitsInner);
      h_maxDlClosestHitToVtx_[match][1]->Fill(maxDlClosestHitToVtx);
      h_maxDlClosestHitToVtxSig_[match][1]->Fill(maxDlClosestHitToVtxSig);
      h_nSharedHits_[match][1]->Fill(aConv.nSharedHits());

      /*
      if ( aConv.caloCluster().size() ) {
    h_convSCdPhi_[match][1]->Fill(  aConv.caloCluster()[0]->phi() - refittedMom.phi() );
    double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
    h_convSCdEta_[match][1]->Fill( aConv.caloCluster()[0]->eta() - ConvEta );
      }
      */

      if (matchConvSC) {
        h_EoverPTracks_[match][1]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
        h_convSCdPhi_[match][1]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
        double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
        h_convSCdEta_[match][1]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
      }
    }

    if (phoIsInEndcap) {
      h_invMass_[match][2]->Fill(invM);
      h_vtxChi2Prob_[match][2]->Fill(chi2Prob);
      h_distMinAppTracks_[match][2]->Fill(aConv.distOfMinimumApproach());
      h_DPhiTracksAtVtx_[match][2]->Fill(dPhiTracksAtVtx);
      h_DCotTracks_[match][2]->Fill(aConv.pairCotThetaSeparation());
      h_lxybs_[match][2]->Fill(lxy);
      h_maxNHitsBeforeVtx_[match][2]->Fill(maxNHitsBeforeVtx);
      h_leadNHitsBeforeVtx_[match][2]->Fill(leadNHitsBeforeVtx);
      h_trailNHitsBeforeVtx_[match][2]->Fill(trailNHitsBeforeVtx);
      h_sumNHitsBeforeVtx_[match][2]->Fill(sumNHitsBeforeVtx);
      h_deltaExpectedHitsInner_[match][2]->Fill(deltaExpectedHitsInner);
      h_leadExpectedHitsInner_[match][2]->Fill(leadExpectedHitsInner);
      h_maxDlClosestHitToVtx_[match][2]->Fill(maxDlClosestHitToVtx);
      h_maxDlClosestHitToVtxSig_[match][2]->Fill(maxDlClosestHitToVtxSig);
      h_nSharedHits_[match][2]->Fill(aConv.nSharedHits());
      if (matchConvSC) {
        h_EoverPTracks_[match][2]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
        h_convSCdPhi_[match][2]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
        double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
        h_convSCdEta_[match][2]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
      }
    }

    h_convVtxRvsZ_[0]->Fill(fabs(aConv.conversionVertex().position().z()),
                            sqrt(aConv.conversionVertex().position().perp2()));
    h_convVtxYvsX_->Fill(aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y());
    h_convVtxYvsX_zoom_[0]->Fill(aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y());
    h_convVtxYvsX_zoom_[1]->Fill(aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y());

    // quantities per track: all conversions
    for (unsigned int i = 0; i < tracks.size(); i++) {
      double d0;
      if (valid_pvtx) {
        d0 = -tracks[i]->dxy(the_pvtx.position());
      } else {
        d0 = tracks[i]->d0();
      }
      h_TkD0_[match]->Fill(d0 * tracks[i]->charge());
      h_nHitsBeforeVtx_[match]->Fill(aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(i) : 0);
      h_dlClosestHitToVtx_[match]->Fill(aConv.dlClosestHitToVtx().size() > 1 ? aConv.dlClosestHitToVtx().at(i).value()
                                                                             : 0);
      h_dlClosestHitToVtxSig_[match]->Fill(aConv.dlClosestHitToVtx().size() > 1
                                               ? aConv.dlClosestHitToVtx().at(i).value() /
                                                     aConv.dlClosestHitToVtx().at(i).error()
                                               : 0);

      nHitsVsEta_[match]->Fill(mcEta_, float(tracks[i]->numberOfValidHits()));
      nHitsVsR_[match]->Fill(mcConvR_, float(tracks[i]->numberOfValidHits()));
      p_nHitsVsEta_[match]->Fill(mcEta_, float(tracks[i]->numberOfValidHits()) - 0.0001);
      p_nHitsVsR_[match]->Fill(mcConvR_, float(tracks[i]->numberOfValidHits()) - 0.0001);
      h_tkChi2_[match]->Fill(tracks[i]->normalizedChi2());
      h_tkChi2Large_[match]->Fill(tracks[i]->normalizedChi2());
      h2_Chi2VsEta_[match]->Fill(mcEta_, tracks[i]->normalizedChi2());
      h2_Chi2VsR_[match]->Fill(mcConvR_, tracks[i]->normalizedChi2());
      p_Chi2VsEta_[match]->Fill(mcEta_, tracks[i]->normalizedChi2());
      p_Chi2VsR_[match]->Fill(mcConvR_, tracks[i]->normalizedChi2());
    }

    bool associated = false;
    float mcConvPt_ = -99999999;
    //    float mcPhi= 0; // unused
    float simPV_Z = 0;
    for (std::vector<PhotonMCTruth>::const_iterator mcPho = mcPhotons.begin(); mcPho != mcPhotons.end(); mcPho++) {
      mcConvPt_ = (*mcPho).fourMomentum().et();
      float mcPhi = (*mcPho).fourMomentum().phi();
      simPV_Z = (*mcPho).primaryVertex().z();
      mcPhi_ = phiNormalization(mcPhi);
      mcEta_ = (*mcPho).fourMomentum().pseudoRapidity();
      mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z());
      mcConvR_ = (*mcPho).vertex().perp();
      mcConvX_ = (*mcPho).vertex().x();
      mcConvY_ = (*mcPho).vertex().y();
      mcConvZ_ = (*mcPho).vertex().z();
      mcConvEta_ = (*mcPho).vertex().eta();
      mcConvPhi_ = (*mcPho).vertex().phi();
      if (fabs(mcEta_) > END_HI)
        continue;
      if (mcConvPt_ < minPhoPtForPurity)
        continue;
      if (fabs(mcEta_) > maxPhoEtaForPurity)
        continue;
      if (fabs(mcConvZ_) > maxPhoZForPurity)
        continue;
      if (mcConvR_ > maxPhoRForEffic)
        continue;

      if ((*mcPho).isAConversion() != 1)
        continue;
      if (!((fabs(mcEta_) <= BARL && mcConvR_ < 85) ||
            (fabs(mcEta_) > BARL && fabs(mcEta_) <= END_HI && fabs((*mcPho).vertex().z()) < 210)))
        continue;

      theConvTP_.clear();
      for (const TrackingParticleRef tp : tpForFakeRate) {
        if (fabs(tp->vx() - (*mcPho).vertex().x()) < 0.0001 && fabs(tp->vy() - (*mcPho).vertex().y()) < 0.0001 &&
            fabs(tp->vz() - (*mcPho).vertex().z()) < 0.0001) {
          theConvTP_.push_back(tp);
        }
      }

      if (theConvTP_.size() < 2)
        continue;

      //associated = false;
      reco::RecoToSimCollection const& p1 = theTrackAssociator->associateRecoToSim(tc1, theConvTP_);
      reco::RecoToSimCollection const& p2 = theTrackAssociator->associateRecoToSim(tc2, theConvTP_);

      auto itP1 = p1.find(tk1);
      auto itP2 = p2.find(tk2);
      bool good = (itP1 != p1.end()) and (not itP1->val.empty()) and (itP2 != p2.end()) and (not itP2->val.empty());
      if (not good) {
        itP1 = p1.find(tk2);
        itP2 = p2.find(tk1);
        good = (itP1 != p1.end()) and (not itP1->val.empty()) and (itP2 != p2.end()) and (not itP2->val.empty());
      }

      if (good) {
        std::vector<std::pair<TrackingParticleRef, double> > const& tp1 = itP1->val;
        std::vector<std::pair<TrackingParticleRef, double> > const& tp2 = itP2->val;
        TrackingParticleRef tpr1 = tp1.front().first;
        TrackingParticleRef tpr2 = tp2.front().first;
        if (abs(tpr1->pdgId()) == 11 && abs(tpr2->pdgId()) == 11 && tpr1->pdgId() * tpr2->pdgId() < 0) {
          if ((tpr1->parentVertex()->sourceTracks_end() - tpr1->parentVertex()->sourceTracks_begin() == 1) &&
              (tpr2->parentVertex()->sourceTracks_end() - tpr2->parentVertex()->sourceTracks_begin() == 1)) {
            if (tpr1->parentVertex().key() == tpr2->parentVertex().key() &&
                ((*tpr1->parentVertex()->sourceTracks_begin())->pdgId() == 22)) {
              mcConvR_ = sqrt(tpr1->parentVertex()->position().Perp2());
              mcConvZ_ = tpr1->parentVertex()->position().z();
              mcConvX_ = tpr1->parentVertex()->position().x();
              mcConvY_ = tpr1->parentVertex()->position().y();
              mcConvEta_ = tpr1->parentVertex()->position().eta();
              mcConvPhi_ = tpr1->parentVertex()->position().phi();
              mcConvPt_ = sqrt((*tpr1->parentVertex()->sourceTracks_begin())->momentum().Perp2());
              //std::cout << " Reco to Sim mcconvpt " << mcConvPt_ << std::endl;
              //cout << "associated track1 to " << tpr1->pdgId() << " with p=" << tpr1->p4() << " with pT=" << tpr1->pt() << endl;
              //cout << "associated track2 to " << tpr2->pdgId() << " with p=" << tpr2->p4() << " with pT=" << tpr2->pt() << endl;
              associated = true;
              break;
            }
          }
        }
      }

    }  // end loop on sim photons

    if (false) {
      theConvTP_.clear();
      theConvTP_ = tpForFakeRate;
      reco::RecoToSimCollection p1incl = theTrackAssociator->associateRecoToSim(tc1, theConvTP_);
      reco::RecoToSimCollection p2incl = theTrackAssociator->associateRecoToSim(tc2, theConvTP_);

      for (std::vector<PhotonMCTruth>::const_iterator mcPho = mcPhotons.begin(); mcPho != mcPhotons.end(); mcPho++) {
        mcConvPt_ = (*mcPho).fourMomentum().et();
        float mcPhi = (*mcPho).fourMomentum().phi();
        simPV_Z = (*mcPho).primaryVertex().z();
        mcPhi_ = phiNormalization(mcPhi);
        mcEta_ = (*mcPho).fourMomentum().pseudoRapidity();
        mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z());
        mcConvR_ = (*mcPho).vertex().perp();
        mcConvX_ = (*mcPho).vertex().x();
        mcConvY_ = (*mcPho).vertex().y();
        mcConvZ_ = (*mcPho).vertex().z();
        mcConvEta_ = (*mcPho).vertex().eta();
        mcConvPhi_ = (*mcPho).vertex().phi();
        if (fabs(mcEta_) > END_HI)
          continue;
        if (mcConvPt_ < minPhoPtForPurity)
          continue;
        if (fabs(mcEta_) > maxPhoEtaForPurity)
          continue;
        if (fabs(mcConvZ_) > maxPhoZForPurity)
          continue;
        if (mcConvR_ > maxPhoRForEffic)
          continue;

        if ((*mcPho).isAConversion() != 1)
          continue;
        if (!((fabs(mcEta_) <= BARL && mcConvR_ < 85) ||
              (fabs(mcEta_) > BARL && fabs(mcEta_) <= END_HI && fabs((*mcPho).vertex().z()) < 210)))
          continue;

        theConvTP_.clear();
        for (TrackingParticleRef tp : tpForFakeRate) {
          if (fabs(tp->vx() - (*mcPho).vertex().x()) < 0.0001 && fabs(tp->vy() - (*mcPho).vertex().y()) < 0.0001 &&
              fabs(tp->vz() - (*mcPho).vertex().z()) < 0.0001) {
            theConvTP_.push_back(tp);
          }
        }

        if (theConvTP_.size() < 2)
          continue;

        //associated = false;
        reco::RecoToSimCollection p1 = theTrackAssociator->associateRecoToSim(tc1, theConvTP_);
        reco::RecoToSimCollection p2 = theTrackAssociator->associateRecoToSim(tc2, theConvTP_);

        if ((!p1incl.empty() && !p2incl.empty()) && (!p1.empty() || !p2.empty())) {  // associated = true;
          try {
            std::vector<std::pair<TrackingParticleRef, double> > tp1 = p1incl[tk1];
            std::vector<std::pair<TrackingParticleRef, double> > tp2 = p2incl[tk2];
            if (!(!tp1.empty() && !tp2.empty())) {
              tp1 = p1[tk2];
              tp2 = p2[tk1];
            }
            if (!tp1.empty() && !tp2.empty()) {
              TrackingParticleRef tpr1 = tp1.front().first;
              TrackingParticleRef tpr2 = tp2.front().first;
              if (abs(tpr1->pdgId()) == 11 && abs(tpr2->pdgId()) == 11 && tpr1->pdgId() * tpr2->pdgId() < 0) {
                if (((tpr1->parentVertex()->sourceTracks_end() - tpr1->parentVertex()->sourceTracks_begin() >= 1) &&
                     (*tpr1->parentVertex()->sourceTracks_begin())->pdgId() == 22) &&
                    ((tpr2->parentVertex()->sourceTracks_end() - tpr2->parentVertex()->sourceTracks_begin() >= 1) &&
                     (*tpr2->parentVertex()->sourceTracks_begin())->pdgId() == 22)) {
                  // if ( fabs(tpr1->vx() - tpr2->vx()) < 0.1 && fabs(tpr1->vy() - tpr2->vy()) < 0.1 && fabs(tpr1->vz() - tpr2->vz()) < 0.1) {
                  //if (((*tpr1->parentVertex()->sourceTracks_begin())->pdgId()==22) || ((*tpr2->parentVertex()->sourceTracks_begin())->pdgId()==22)) {
                  //                       mcConvR_ = sqrt(tpr1->parentVertex()->position().Perp2());
                  //                       mcConvZ_ = tpr1->parentVertex()->position().z();
                  //                       mcConvX_ = tpr1->parentVertex()->position().x();
                  //                       mcConvY_ = tpr1->parentVertex()->position().y();
                  //                       mcConvEta_ = tpr1->parentVertex()->position().eta();
                  //                       mcConvPhi_ = tpr1->parentVertex()->position().phi();
                  //                       mcConvPt_ = sqrt((*tpr1->parentVertex()->sourceTracks_begin())->momentum().Perp2());
                  //std::cout << " Reco to Sim mcconvpt " << mcConvPt_ << std::endl;
                  //cout << "associated track1 to " << tpr1->pdgId() << " with p=" << tpr1->p4() << " with pT=" << tpr1->pt() << endl;
                  //cout << "associated track2 to " << tpr2->pdgId() << " with p=" << tpr2->p4() << " with pT=" << tpr2->pt() << endl;
                  associated = true;
                  break;
                  //}
                  //}
                }
              }
            }
          } catch (Exception const& event) {
            //cout << "do not continue: " << event.what()  << endl;
            //continue;
          }
        }
      }
    }

    if (associated)
      match = 1;
    else
      match = 2;

    h_match_->Fill(float(match));
    if (match == 1)
      nRecConvAss_++;
    h_convEta_[match][0]->Fill(refittedMom.eta());
    h_convEta_[match][1]->Fill(refittedMom.eta());
    if (matchConvSC)
      h_convEtaMatchSC_[match][0]->Fill(refittedMom.eta());
    h_convPhi_[match][0]->Fill(refittedMom.phi());
    h_convR_[match][0]->Fill(sqrt(aConv.conversionVertex().position().perp2()));
    h_convZ_[match][0]->Fill(aConv.conversionVertex().position().z());
    h_convPt_[match][0]->Fill(sqrt(refittedMom.perp2()));
    h_invMass_[match][0]->Fill(invM);
    h_vtxChi2Prob_[match][0]->Fill(chi2Prob);
    h_DPhiTracksAtVtx_[match][0]->Fill(dPhiTracksAtVtx);
    h_DCotTracks_[match][0]->Fill(aConv.pairCotThetaSeparation());
    h_distMinAppTracks_[match][0]->Fill(aConv.distOfMinimumApproach());
    h_lxybs_[match][0]->Fill(lxy);
    h_maxNHitsBeforeVtx_[match][0]->Fill(maxNHitsBeforeVtx);
    h_leadNHitsBeforeVtx_[match][0]->Fill(leadNHitsBeforeVtx);
    h_trailNHitsBeforeVtx_[match][0]->Fill(trailNHitsBeforeVtx);
    h_sumNHitsBeforeVtx_[match][0]->Fill(sumNHitsBeforeVtx);
    h_deltaExpectedHitsInner_[match][0]->Fill(deltaExpectedHitsInner);
    h_leadExpectedHitsInner_[match][0]->Fill(leadExpectedHitsInner);
    h_maxDlClosestHitToVtx_[match][0]->Fill(maxDlClosestHitToVtx);
    h_maxDlClosestHitToVtxSig_[match][0]->Fill(maxDlClosestHitToVtxSig);
    h_nSharedHits_[match][0]->Fill(aConv.nSharedHits());
    if (matchConvSC) {
      //h_EoverPTracks_[match][0] ->Fill (aConv.EoverPrefittedTracks());
      h_EoverPTracks_[match][0]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
      h_convSCdPhi_[match][0]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
      double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
      h_convSCdEta_[match][0]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
    }
    if (match == 1) {
      h2_photonPtRecVsPtSim_->Fill(mcConvPt_, sqrt(refittedMom.perp2()));
      h_convPtRes_[0]->Fill(sqrt(refittedMom.perp2()) / mcConvPt_);
    }

    if (phoIsInBarrel) {
      h_invMass_[match][1]->Fill(invM);
      h_vtxChi2Prob_[match][1]->Fill(chi2Prob);
      h_DPhiTracksAtVtx_[match][1]->Fill(dPhiTracksAtVtx);
      h_DCotTracks_[match][1]->Fill(aConv.pairCotThetaSeparation());
      h_distMinAppTracks_[match][1]->Fill(aConv.distOfMinimumApproach());
      h_lxybs_[match][1]->Fill(lxy);
      h_maxNHitsBeforeVtx_[match][1]->Fill(maxNHitsBeforeVtx);
      h_leadNHitsBeforeVtx_[match][1]->Fill(leadNHitsBeforeVtx);
      h_trailNHitsBeforeVtx_[match][1]->Fill(trailNHitsBeforeVtx);
      h_sumNHitsBeforeVtx_[match][1]->Fill(sumNHitsBeforeVtx);
      h_deltaExpectedHitsInner_[match][1]->Fill(deltaExpectedHitsInner);
      h_leadExpectedHitsInner_[match][1]->Fill(leadExpectedHitsInner);
      h_maxDlClosestHitToVtx_[match][1]->Fill(maxDlClosestHitToVtx);
      h_maxDlClosestHitToVtxSig_[match][1]->Fill(maxDlClosestHitToVtxSig);
      h_nSharedHits_[match][1]->Fill(aConv.nSharedHits());
      if (matchConvSC) {
        //  h_EoverPTracks_[match][1] ->Fill (aConv.EoverPrefittedTracks());
        h_EoverPTracks_[match][1]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
        h_convSCdPhi_[match][1]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
        double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
        h_convSCdEta_[match][1]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
      }
      if (match == 1)
        h_convPtRes_[1]->Fill(sqrt(refittedMom.perp2()) / mcConvPt_);
    }

    if (phoIsInEndcap) {
      h_invMass_[match][2]->Fill(invM);
      h_vtxChi2Prob_[match][2]->Fill(chi2Prob);
      h_DPhiTracksAtVtx_[match][2]->Fill(dPhiTracksAtVtx);
      h_DCotTracks_[match][2]->Fill(aConv.pairCotThetaSeparation());
      h_distMinAppTracks_[match][2]->Fill(aConv.distOfMinimumApproach());
      h_lxybs_[match][2]->Fill(lxy);
      h_maxNHitsBeforeVtx_[match][2]->Fill(maxNHitsBeforeVtx);
      h_leadNHitsBeforeVtx_[match][2]->Fill(leadNHitsBeforeVtx);
      h_trailNHitsBeforeVtx_[match][2]->Fill(trailNHitsBeforeVtx);
      h_sumNHitsBeforeVtx_[match][2]->Fill(sumNHitsBeforeVtx);
      h_deltaExpectedHitsInner_[match][2]->Fill(deltaExpectedHitsInner);
      h_leadExpectedHitsInner_[match][2]->Fill(leadExpectedHitsInner);
      h_maxDlClosestHitToVtx_[match][2]->Fill(maxDlClosestHitToVtx);
      h_maxDlClosestHitToVtxSig_[match][2]->Fill(maxDlClosestHitToVtxSig);
      h_nSharedHits_[match][2]->Fill(aConv.nSharedHits());
      if (matchConvSC) {
        //  h_EoverPTracks_[match][2] ->Fill (aConv.EoverPrefittedTracks());
        h_EoverPTracks_[match][2]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
        h_convSCdPhi_[match][2]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
        double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
        h_convSCdEta_[match][2]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
      }
      if (match == 1)
        h_convPtRes_[2]->Fill(sqrt(refittedMom.perp2()) / mcConvPt_);
    }

    if (match == 1) {
      h_convVtxdX_->Fill(aConv.conversionVertex().position().x() - mcConvX_);
      h_convVtxdY_->Fill(aConv.conversionVertex().position().y() - mcConvY_);
      h_convVtxdZ_->Fill(aConv.conversionVertex().position().z() - mcConvZ_);
      h_convVtxdR_->Fill(sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
      h_convVtxdPhi_->Fill(aConv.conversionVertex().position().phi() - mcConvPhi_);
      h_convVtxdEta_->Fill(aConv.conversionVertex().position().eta() - mcConvEta_);
      h2_convVtxdRVsR_->Fill(mcConvR_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
      h2_convVtxdRVsEta_->Fill(mcEta_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
      p_convVtxdRVsR_->Fill(mcConvR_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
      p_convVtxdRVsEta_->Fill(mcEta_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
      p_convVtxdXVsX_->Fill(mcConvX_, aConv.conversionVertex().position().x() - mcConvX_);
      p_convVtxdYVsY_->Fill(mcConvY_, aConv.conversionVertex().position().y() - mcConvY_);
      p_convVtxdZVsZ_->Fill(mcConvZ_, aConv.conversionVertex().position().z() - mcConvZ_);
      p_convVtxdZVsR_->Fill(mcConvR_, aConv.conversionVertex().position().z() - mcConvZ_);

      float dR = sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_;
      float dZ = aConv.conversionVertex().position().z() - mcConvZ_;
      p2_convVtxdRVsRZ_->Fill(mcConvZ_, mcConvR_, dR);
      p2_convVtxdZVsRZ_->Fill(mcConvZ_, mcConvR_, dZ);

      h2_convVtxRrecVsTrue_->Fill(mcConvR_, sqrt(aConv.conversionVertex().position().perp2()));

      h_zPVFromTracks_[match]->Fill(aConv.zOfPrimaryVertexFromTracks());
      h_dzPVFromTracks_[match]->Fill(aConv.zOfPrimaryVertexFromTracks() - simPV_Z);
      h2_dzPVVsR_->Fill(mcConvR_, aConv.zOfPrimaryVertexFromTracks() - simPV_Z);
      p_dzPVVsR_->Fill(mcConvR_, aConv.zOfPrimaryVertexFromTracks() - simPV_Z);

      if (phoIsInBarrel) {
        h_convVtxdX_barrel_->Fill(aConv.conversionVertex().position().x() - mcConvX_);
        h_convVtxdY_barrel_->Fill(aConv.conversionVertex().position().y() - mcConvY_);
        h_convVtxdZ_barrel_->Fill(aConv.conversionVertex().position().z() - mcConvZ_);
        h_convVtxdR_barrel_->Fill(sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
      }
      if (phoIsInEndcap) {
        h_convVtxdX_endcap_->Fill(aConv.conversionVertex().position().x() - mcConvX_);
        h_convVtxdY_endcap_->Fill(aConv.conversionVertex().position().y() - mcConvY_);
        h_convVtxdZ_endcap_->Fill(aConv.conversionVertex().position().z() - mcConvZ_);
        h_convVtxdR_endcap_->Fill(sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
      }
    }

    for (unsigned int i = 0; i < tracks.size(); i++) {
      //std::cout << " Loop over tracks  pt " << tracks[i]->pt() << std::endl;
      RefToBase<reco::Track> tfrb(aConv.tracks()[i]);
      itAss = myAss.find(tfrb.get());

      nHitsVsEta_[match]->Fill(mcEta_, float(tracks[i]->numberOfValidHits()));
      nHitsVsR_[match]->Fill(mcConvR_, float(tracks[i]->numberOfValidHits()));
      p_nHitsVsEta_[match]->Fill(mcEta_, float(tracks[i]->numberOfValidHits()) - 0.0001);
      p_nHitsVsR_[match]->Fill(mcConvR_, float(tracks[i]->numberOfValidHits()) - 0.0001);
      h_tkChi2_[match]->Fill(tracks[i]->normalizedChi2());
      h_tkChi2Large_[match]->Fill(tracks[i]->normalizedChi2());
      h2_Chi2VsEta_[match]->Fill(mcEta_, tracks[i]->normalizedChi2());
      h2_Chi2VsR_[match]->Fill(mcConvR_, tracks[i]->normalizedChi2());
      p_Chi2VsEta_[match]->Fill(mcEta_, tracks[i]->normalizedChi2());
      p_Chi2VsR_[match]->Fill(mcConvR_, tracks[i]->normalizedChi2());
      double d0;
      if (valid_pvtx) {
        d0 = -tracks[i]->dxy(the_pvtx.position());
      } else {
        d0 = tracks[i]->d0();
      }
      h_TkD0_[match]->Fill(d0 * tracks[i]->charge());
      h_nHitsBeforeVtx_[match]->Fill(aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(i) : 0);
      h_dlClosestHitToVtx_[match]->Fill(aConv.dlClosestHitToVtx().size() > 1 ? aConv.dlClosestHitToVtx().at(i).value()
                                                                             : 0);
      h_dlClosestHitToVtxSig_[match]->Fill(aConv.dlClosestHitToVtx().size() > 1
                                               ? aConv.dlClosestHitToVtx().at(i).value() /
                                                     aConv.dlClosestHitToVtx().at(i).error()
                                               : 0);

      if (itAss == myAss.end())
        continue;
      reco::Track refTrack = aConv.conversionVertex().refittedTracks()[i];

      float simPt = sqrt(((*itAss).second)->momentum().perp2());
      float recPt = refTrack.pt();
      float ptres = recPt - simPt;
      //float pterror = aConv.tracks()[i]->ptError();
      float pterror = aConv.conversionVertex().refittedTracks()[i].ptError();
      h2_PtRecVsPtSim_[0]->Fill(simPt, recPt);
      h_TkPtPull_[0]->Fill(ptres / pterror);
      h2_TkPtPull_[0]->Fill(mcEta_, ptres / pterror);

      if (phoIsInBarrel) {
        h_TkPtPull_[1]->Fill(ptres / pterror);
        h2_PtRecVsPtSim_[1]->Fill(simPt, recPt);
      }
      if (phoIsInEndcap) {
        h_TkPtPull_[2]->Fill(ptres / pterror);
        h2_PtRecVsPtSim_[2]->Fill(simPt, recPt);
      }
    }  // end loop over track

  }  // loop over reco conversions

  h_nConv_[0][0]->Fill(float(nRecConv_));
  h_nConv_[1][0]->Fill(float(nRecConvAss_));
}

void TkConvValidator::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  run, edm::EventSetup const &  es  )

overridevirtual

Histograms for efficiencies

Denominators

zooms

Implements DQMEDAnalyzer.

Definition at line 148 of file TkConvValidator.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile2D(), muonRecoAnalyzer_cfi::chi2Max, muonRecoAnalyzer_cfi::chi2Min, dbe_, photonValidator_cfi::dCotTracksBin, photonValidator_cfi::dCotTracksMax, photonValidator_cfi::dCotTracksMin, photonAnalyzer_cfi::dPhiTracksBin, photonAnalyzer_cfi::dPhiTracksMax, photonAnalyzer_cfi::dPhiTracksMin, photonValidator_cfi::eoverpBin, photonValidationSequence_cff::eoverpMax, photonValidationSequence_cff::eoverpMin, muonRecoAnalyzer_cfi::etaBin, conversionPostprocessing_cfi::etaBin2, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, photonAnalyzer_cfi::etBin, cosmicPhotonAnalyzer_cfi::etMax, photonAnalyzer_cfi::etMin, tkConvValidator_cfi::maxPhoEtaForEffic, tkConvValidator_cfi::maxPhoEtaForPurity, tkConvValidator_cfi::maxPhoRForEffic, tkConvValidator_cfi::maxPhoRForPurity, tkConvValidator_cfi::maxPhoZForEffic, tkConvValidator_cfi::maxPhoZForPurity, tkConvValidator_cfi::minPhoPtForEffic, tkConvValidator_cfi::minPhoPtForPurity, Utilities::operator, BeamMonitor_cff::phiBin, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, photonAnalyzer_cfi::rBin, photonValidator_cfi::rBinForXray, photonValidator_cfi::resBin, photonValidator_cfi::resMax, photonValidator_cfi::resMin, photonAnalyzer_cfi::rMax, photonValidator_cfi::rMaxForXray, photonAnalyzer_cfi::rMin, photonValidator_cfi::rMinForXray, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, photonAnalyzer_cfi::zBin, photonValidator_cfi::zBin2ForXray, photonValidator_cfi::zBinForXray, photonAnalyzer_cfi::zMax, photonValidator_cfi::zMaxForXray, photonAnalyzer_cfi::zMin, and photonValidator_cfi::zMinForXray.

                                                                                                         {
  nEvt_ = 0;
  nEntry_ = 0;
  nRecConv_ = 0;
  nRecConvAss_ = 0;
  nRecConvAssWithEcal_ = 0;

  nInvalidPCA_ = 0;

  dbe_ = nullptr;
  dbe_ = edm::Service<DQMStore>().operator->();

  double etMin = parameters_.getParameter<double>("etMin");
  double etMax = parameters_.getParameter<double>("etMax");
  int etBin = parameters_.getParameter<int>("etBin");

  double resMin = parameters_.getParameter<double>("resMin");
  double resMax = parameters_.getParameter<double>("resMax");
  int resBin = parameters_.getParameter<int>("resBin");

  double etaMin = parameters_.getParameter<double>("etaMin");
  double etaMax = parameters_.getParameter<double>("etaMax");
  int etaBin = parameters_.getParameter<int>("etaBin");
  int etaBin2 = parameters_.getParameter<int>("etaBin2");

  double phiMin = parameters_.getParameter<double>("phiMin");
  double phiMax = parameters_.getParameter<double>("phiMax");
  int phiBin = parameters_.getParameter<int>("phiBin");

  double rMin = parameters_.getParameter<double>("rMin");
  double rMax = parameters_.getParameter<double>("rMax");
  int rBin = parameters_.getParameter<int>("rBin");

  double zMin = parameters_.getParameter<double>("zMin");
  double zMax = parameters_.getParameter<double>("zMax");
  int zBin = parameters_.getParameter<int>("zBin");

  double dPhiTracksMin = parameters_.getParameter<double>("dPhiTracksMin");
  double dPhiTracksMax = parameters_.getParameter<double>("dPhiTracksMax");
  int dPhiTracksBin = parameters_.getParameter<int>("dPhiTracksBin");

  double eoverpMin = parameters_.getParameter<double>("eoverpMin");
  double eoverpMax = parameters_.getParameter<double>("eoverpMax");
  int eoverpBin = parameters_.getParameter<int>("eoverpBin");

  //  double dEtaTracksMin = parameters_.getParameter<double>("dEtaTracksMin");  // unused
  //  double dEtaTracksMax = parameters_.getParameter<double>("dEtaTracksMax"); // unused
  //  int    dEtaTracksBin = parameters_.getParameter<int>("dEtaTracksBin");  // unused

  double dCotTracksMin = parameters_.getParameter<double>("dCotTracksMin");
  double dCotTracksMax = parameters_.getParameter<double>("dCotTracksMax");
  int dCotTracksBin = parameters_.getParameter<int>("dCotTracksBin");

  double chi2Min = parameters_.getParameter<double>("chi2Min");
  double chi2Max = parameters_.getParameter<double>("chi2Max");

  double rMinForXray = parameters_.getParameter<double>("rMinForXray");
  double rMaxForXray = parameters_.getParameter<double>("rMaxForXray");
  int rBinForXray = parameters_.getParameter<int>("rBinForXray");
  double zMinForXray = parameters_.getParameter<double>("zMinForXray");
  double zMaxForXray = parameters_.getParameter<double>("zMaxForXray");
  int zBinForXray = parameters_.getParameter<int>("zBinForXray");
  int zBin2ForXray = parameters_.getParameter<int>("zBin2ForXray");

  minPhoPtForEffic = parameters_.getParameter<double>("minPhoPtForEffic");
  maxPhoEtaForEffic = parameters_.getParameter<double>("maxPhoEtaForEffic");
  maxPhoZForEffic = parameters_.getParameter<double>("maxPhoZForEffic");
  maxPhoRForEffic = parameters_.getParameter<double>("maxPhoRForEffic");
  minPhoPtForPurity = parameters_.getParameter<double>("minPhoPtForPurity");
  maxPhoEtaForPurity = parameters_.getParameter<double>("maxPhoEtaForPurity");
  maxPhoZForPurity = parameters_.getParameter<double>("maxPhoZForPurity");
  maxPhoRForPurity = parameters_.getParameter<double>("maxPhoRForPurity");

  if (dbe_) {
    // SC from reco photons

    //TString simfolder = TString(
    std::string simpath = dqmpath_ + "SimulationInfo";
    iBooker.setCurrentFolder(simpath);
    //
    // simulation information about conversions
    std::string histname = "nOfSimConversions";
    h_nSimConv_[0] = iBooker.book1D(histname, "# of Sim conversions per event ", 20, -0.5, 19.5);
    histname = "h_AllSimConvEta";
    h_AllSimConv_[0] = iBooker.book1D(histname, " All conversions: simulated #eta", etaBin2, etaMin, etaMax);
    histname = "h_AllSimConvPhi";
    h_AllSimConv_[1] = iBooker.book1D(histname, " All conversions: simulated #phi", phiBin, phiMin, phiMax);
    histname = "h_AllSimConvR";
    h_AllSimConv_[2] = iBooker.book1D(histname, " All conversions: simulated R", rBin, rMin, rMax);
    histname = "h_AllSimConvZ";
    h_AllSimConv_[3] = iBooker.book1D(histname, " All conversions: simulated Z", zBin, zMin, zMax);
    histname = "h_AllSimConvEt";
    h_AllSimConv_[4] = iBooker.book1D(histname, " All conversions: simulated Et", etBin, etMin, etMax);
    //
    histname = "nOfVisSimConversions";
    h_nSimConv_[1] = iBooker.book1D(histname, "# of Sim conversions per event ", 20, -0.5, 19.5);
    histname = "h_VisSimConvEta";
    h_VisSimConv_[0] = iBooker.book1D(histname, " All vis conversions: simulated #eta", etaBin2, etaMin, etaMax);
    histname = "h_VisSimConvPhi";
    h_VisSimConv_[1] = iBooker.book1D(histname, " All vis conversions: simulated #phi", phiBin, phiMin, phiMax);
    histname = "h_VisSimConvR";
    h_VisSimConv_[2] = iBooker.book1D(histname, " All vis conversions: simulated R", rBin, rMin, rMax);
    histname = "h_VisSimConvZ";
    h_VisSimConv_[3] = iBooker.book1D(histname, " All vis conversions: simulated Z", zBin, zMin, zMax);
    histname = "h_VisSimConvEt";
    h_VisSimConv_[4] = iBooker.book1D(histname, " All vis conversions: simulated Et", etBin, etMin, etMax);

    //
    histname = "h_SimConvTwoMTracksEta";
    h_SimConvTwoMTracks_[0] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks: simulated #eta", etaBin2, etaMin, etaMax);
    histname = "h_SimConvTwoMTracksPhi";
    h_SimConvTwoMTracks_[1] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks: simulated #phi", phiBin, phiMin, phiMax);
    histname = "h_SimConvTwoMTracksR";
    h_SimConvTwoMTracks_[2] =
        iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated R", rBin, rMin, rMax);
    histname = "h_SimConvTwoMTracksZ";
    h_SimConvTwoMTracks_[3] =
        iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated Z", zBin, zMin, zMax);
    histname = "h_SimConvTwoMTracksEt";
    h_SimConvTwoMTracks_[4] =
        iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated Et", etBin, etMin, etMax);
    //
    histname = "h_SimConvTwoTracksEta";
    h_SimConvTwoTracks_[0] =
        iBooker.book1D(histname, " All vis conversions with 2 reco  tracks: simulated #eta", etaBin2, etaMin, etaMax);
    histname = "h_SimConvTwoTracksPhi";
    h_SimConvTwoTracks_[1] =
        iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated #phi", phiBin, phiMin, phiMax);
    histname = "h_SimConvTwoTracksR";
    h_SimConvTwoTracks_[2] =
        iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated R", rBin, rMin, rMax);
    histname = "h_SimConvTwoTracksZ";
    h_SimConvTwoTracks_[3] =
        iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated Z", zBin, zMin, zMax);
    histname = "h_SimConvTwoTracksEt";
    h_SimConvTwoTracks_[4] =
        iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated Et", etBin, etMin, etMax);
    //
    histname = "h_SimConvTwoMTracksEtaAndVtxPGT0";
    h_SimConvTwoMTracksAndVtxPGT0_[0] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated #eta", etaBin2, etaMin, etaMax);
    histname = "h_SimConvTwoMTracksPhiAndVtxPGT0";
    h_SimConvTwoMTracksAndVtxPGT0_[1] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated #phi", phiBin, phiMin, phiMax);
    histname = "h_SimConvTwoMTracksRAndVtxPGT0";
    h_SimConvTwoMTracksAndVtxPGT0_[2] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated R", rBin, rMin, rMax);
    histname = "h_SimConvTwoMTracksZAndVtxPGT0";
    h_SimConvTwoMTracksAndVtxPGT0_[3] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated Z", zBin, zMin, zMax);
    histname = "h_SimConvTwoMTracksEtAndVtxPGT0";
    h_SimConvTwoMTracksAndVtxPGT0_[4] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated Et", etBin, etMin, etMax);

    //
    histname = "h_SimConvTwoMTracksEtaAndVtxPGT0005";
    h_SimConvTwoMTracksAndVtxPGT0005_[0] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated #eta", etaBin2, etaMin, etaMax);
    histname = "h_SimConvTwoMTracksPhiAndVtxPGT0005";
    h_SimConvTwoMTracksAndVtxPGT0005_[1] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated #phi", phiBin, phiMin, phiMax);
    histname = "h_SimConvTwoMTracksRAndVtxPGT0005";
    h_SimConvTwoMTracksAndVtxPGT0005_[2] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated R", rBin, rMin, rMax);
    histname = "h_SimConvTwoMTracksZAndVtxPGT0005";
    h_SimConvTwoMTracksAndVtxPGT0005_[3] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated Z", zBin, zMin, zMax);
    histname = "h_SimConvTwoMTracksEtAndVtxPGT0005";
    h_SimConvTwoMTracksAndVtxPGT0005_[4] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated Et", etBin, etMin, etMax);

    histname = "h_SimRecConvTwoMTracksEta";
    h_SimRecConvTwoMTracks_[0] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks: simulated #eta", etaBin2, etaMin, etaMax);
    histname = "h_SimRecConvTwoMTracksPhi";
    h_SimRecConvTwoMTracks_[1] = iBooker.book1D(
        histname, " All vis conversions with 2 reco-matching tracks: simulated #phi", phiBin, phiMin, phiMax);
    histname = "h_SimRecConvTwoMTracksR";
    h_SimRecConvTwoMTracks_[2] =
        iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated R", rBin, rMin, rMax);
    histname = "h_SimRecConvTwoMTracksZ";
    h_SimRecConvTwoMTracks_[3] =
        iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated Z", zBin, zMin, zMax);
    histname = "h_SimRecConvTwoMTracksEt";
    h_SimRecConvTwoMTracks_[4] =
        iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated Et", etBin, etMin, etMax);
    //

    h_SimConvEtaPix_[0] = iBooker.book1D("simConvEtaPix", " sim converted Photon Eta: Pix ", etaBin, etaMin, etaMax);
    h_simTkPt_ = iBooker.book1D("simTkPt", "Sim conversion tracks pt ", etBin * 3, 0., etMax);
    h_simTkEta_ = iBooker.book1D("simTkEta", "Sim conversion tracks eta ", etaBin, etaMin, etaMax);

    h_simConvVtxRvsZ_[0] = iBooker.book2D("simConvVtxRvsZAll",
                                          " Photon Sim conversion vtx position",
                                          zBinForXray,
                                          zMinForXray,
                                          zMaxForXray,
                                          rBinForXray,
                                          rMinForXray,
                                          rMaxForXray);
    h_simConvVtxRvsZ_[1] = iBooker.book2D("simConvVtxRvsZBarrel",
                                          " Photon Sim conversion vtx position",
                                          zBinForXray,
                                          zMinForXray,
                                          zMaxForXray,
                                          rBinForXray,
                                          rMinForXray,
                                          rMaxForXray);
    h_simConvVtxRvsZ_[2] = iBooker.book2D("simConvVtxRvsZEndcap",
                                          " Photon Sim conversion vtx position",
                                          zBin2ForXray,
                                          zMinForXray,
                                          zMaxForXray,
                                          rBinForXray,
                                          rMinForXray,
                                          rMaxForXray);
    h_simConvVtxRvsZ_[3] = iBooker.book2D("simConvVtxRvsZBarrel2",
                                          " Photon Sim conversion vtx position when reco R<4cm",
                                          zBinForXray,
                                          zMinForXray,
                                          zMaxForXray,
                                          rBinForXray,
                                          rMinForXray,
                                          rMaxForXray);
    h_simConvVtxYvsX_ = iBooker.book2D(
        "simConvVtxYvsXTrkBarrel", " Photon Sim conversion vtx position, (x,y) eta<1 ", 100, -80., 80., 100, -80., 80.);

    std::string convpath = dqmpath_ + "ConversionInfo";
    iBooker.setCurrentFolder(convpath);

    histname = "nConv";
    h_nConv_[0][0] = iBooker.book1D(
        histname + "All", "Number Of Conversions per isolated candidates per events: All Ecal  ", 10, -0.5, 9.5);
    h_nConv_[0][1] = iBooker.book1D(
        histname + "Barrel", "Number Of Conversions per isolated candidates per events: Ecal Barrel  ", 10, -0.5, 9.5);
    h_nConv_[0][2] = iBooker.book1D(
        histname + "Endcap", "Number Of Conversions per isolated candidates per events: Ecal Endcap ", 10, -0.5, 9.5);
    h_nConv_[1][0] = iBooker.book1D(histname + "All_Ass",
                                    "Number Of associated Conversions per isolated candidates per events: All Ecal  ",
                                    10,
                                    -0.5,
                                    9.5);

    h_convEta_[0][0] = iBooker.book1D("convEta", " converted Photon  Eta ", etaBin, etaMin, etaMax);
    h_convEtaMatchSC_[0][0] =
        iBooker.book1D("convEtaMatchSC", " converted Photon  Eta when SC is matched ", etaBin, etaMin, etaMax);
    h_convEta2_[0][0] = iBooker.book1D("convEta2", " converted Photon  Eta ", etaBin2, etaMin, etaMax);
    h_convPhi_[0][0] = iBooker.book1D("convPhi", " converted Photon  Phi ", phiBin, phiMin, phiMax);
    h_convR_[0][0] = iBooker.book1D("convR", " converted photon R", rBin, rMin, rMax);
    h_convZ_[0][0] = iBooker.book1D("convZ", " converted photon Z", zBin, zMin, zMax);
    h_convPt_[0][0] = iBooker.book1D("convPt", "  conversions Transverse Energy: all eta ", etBin, etMin, etMax);

    h_convEta_[1][0] = iBooker.book1D("convEtaAss2", " Matched converted Photon  Eta ", etaBin2, etaMin, etaMax);
    h_convEta_[1][1] = iBooker.book1D("convEtaAss", " Matched converted Photon  Eta ", etaBin, etaMin, etaMax);
    h_convEtaMatchSC_[1][0] =
        iBooker.book1D("convEtaMatchSCAss", " converted Photon  Eta when SC is matched ", etaBin, etaMin, etaMax);
    h_convPhi_[1][0] = iBooker.book1D("convPhiAss", " Matched converted Photon  Phi ", phiBin, phiMin, phiMax);
    h_convR_[1][0] = iBooker.book1D("convRAss", " Matched converted photon R", rBin, rMin, rMax);
    h_convZ_[1][0] = iBooker.book1D("convZAss", " Matched converted photon Z", zBin, zMin, zMax);
    h_convPt_[1][0] =
        iBooker.book1D("convPtAss", "Matched conversions Transverse Energy: all eta ", etBin, etMin, etMax);

    h_convEta_[2][0] = iBooker.book1D("convEtaFake2", " Fake converted Photon  Eta ", etaBin2, etaMin, etaMax);
    h_convEta_[2][1] = iBooker.book1D("convEtaFake", " Fake converted Photon  Eta ", etaBin, etaMin, etaMax);
    h_convEtaMatchSC_[2][0] =
        iBooker.book1D("convEtaMatchSCFake", " converted Photon  Eta when SC is matched ", etaBin, etaMin, etaMax);
    h_convPhi_[2][0] = iBooker.book1D("convPhiFake", " Fake converted Photon  Phi ", phiBin, phiMin, phiMax);
    h_convR_[2][0] = iBooker.book1D("convRFake", " Fake converted photon R", rBin, rMin, rMax);
    h_convZ_[2][0] = iBooker.book1D("convZFake", " Fake converted photon Z", zBin, zMin, zMax);
    h_convPt_[2][0] = iBooker.book1D("convPtFake", "Fake conversions Transverse Energy: all eta ", etBin, etMin, etMax);

    h_convRplot_ = iBooker.book1D("convRplot", " converted photon R", 600, 0., 120.);
    h_convZplot_ = iBooker.book1D("convZplot", " converted photon Z", 320, -160., 160.);

    histname = "convSCdPhi";
    h_convSCdPhi_[0][0] = iBooker.book1D(histname + "All", "dPhi between SC and conversion", 100, -0.1, 0.1);
    h_convSCdPhi_[0][1] =
        iBooker.book1D(histname + "Barrel", " dPhi between SC and conversion: Barrel", 100, -0.1, 0.1);
    h_convSCdPhi_[0][2] =
        iBooker.book1D(histname + "Endcap", " dPhi between SC and conversion: Endcap", 100, -0.1, 0.1);
    h_convSCdPhi_[1][0] = iBooker.book1D(histname + "All_Ass", "dPhi between SC and conversion", 100, -0.1, 0.1);
    h_convSCdPhi_[1][1] =
        iBooker.book1D(histname + "Barrel_Ass", " dPhi between SC and conversion: Barrel", 100, -0.1, 0.1);
    h_convSCdPhi_[1][2] =
        iBooker.book1D(histname + "Endcap_Ass", " dPhi between SC and conversion: Endcap", 100, -0.1, 0.1);
    h_convSCdPhi_[2][0] = iBooker.book1D(histname + "All_Fakes", "dPhi between SC and conversion", 100, -0.1, 0.1);
    h_convSCdPhi_[2][1] =
        iBooker.book1D(histname + "Barrel_Fakes", " dPhi between SC and conversion: Barrel", 100, -0.1, 0.1);
    h_convSCdPhi_[2][2] =
        iBooker.book1D(histname + "Endcap_Fakes", " dPhi between SC and conversion: Endcap", 100, -0.1, 0.1);
    histname = "convSCdEta";
    h_convSCdEta_[0][0] = iBooker.book1D(histname + "All", " dEta between SC and conversion", 100, -0.1, 0.1);
    h_convSCdEta_[0][1] =
        iBooker.book1D(histname + "Barrel", " dEta between SC and conversion: Barrel", 100, -0.1, 0.1);
    h_convSCdEta_[0][2] =
        iBooker.book1D(histname + "Endcap", " dEta between SC and conversion: Endcap", 100, -0.1, 0.1);
    h_convSCdEta_[1][0] = iBooker.book1D(histname + "All_Ass", " dEta between SC and conversion", 100, -0.1, 0.1);
    h_convSCdEta_[1][1] =
        iBooker.book1D(histname + "Barrel_Ass", " dEta between SC and conversion: Barrel", 100, -0.1, 0.1);
    h_convSCdEta_[1][2] =
        iBooker.book1D(histname + "Endcap_Ass", " dEta between SC and conversion: Endcap", 100, -0.1, 0.1);
    h_convSCdEta_[2][0] = iBooker.book1D(histname + "All_Fakes", " dEta between SC and conversion", 100, -0.1, 0.1);
    h_convSCdEta_[2][1] =
        iBooker.book1D(histname + "Barrel_Fakes", " dEta between SC and conversion: Barrel", 100, -0.1, 0.1);
    h_convSCdEta_[2][2] =
        iBooker.book1D(histname + "Endcap_Fakes", " dEta between SC and conversion: Endcap", 100, -0.1, 0.1);

    histname = "convPtRes";
    h_convPtRes_[0] = iBooker.book1D(histname + "All", " Conversion Pt rec/true : All ecal ", resBin, resMin, resMax);
    h_convPtRes_[1] = iBooker.book1D(histname + "Barrel", " Conversion Pt rec/true : Barrel ", resBin, resMin, resMax);
    h_convPtRes_[2] = iBooker.book1D(histname + "Endcap", " Conversion Pt rec/true : Endcap ", resBin, resMin, resMax);

    histname = "hInvMass";
    h_invMass_[0][0] = iBooker.book1D(
        histname + "All_AllTracks", " Photons:Tracks from conversion: Pair invariant mass: all Ecal ", 100, 0., 1.5);
    h_invMass_[0][1] = iBooker.book1D(histname + "Barrel_AllTracks",
                                      " Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",
                                      100,
                                      0.,
                                      1.5);
    h_invMass_[0][2] = iBooker.book1D(histname + "Endcap_AllTracks",
                                      " Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",
                                      100,
                                      0.,
                                      1.5);
    //
    h_invMass_[1][0] = iBooker.book1D(
        histname + "All_AssTracks", " Photons:Tracks from conversion: Pair invariant mass: all Ecal ", 100, 0., 1.5);
    h_invMass_[1][1] = iBooker.book1D(histname + "Barrel_AssTracks",
                                      " Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",
                                      100,
                                      0.,
                                      1.5);
    h_invMass_[1][2] = iBooker.book1D(histname + "Endcap_AssTracks",
                                      " Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",
                                      100,
                                      0.,
                                      1.5);
    //
    h_invMass_[2][0] = iBooker.book1D(
        histname + "All_FakeTracks", " Photons:Tracks from conversion: Pair invariant mass: all Ecal ", 100, 0., 1.5);
    h_invMass_[2][1] = iBooker.book1D(histname + "Barrel_FakeTracks",
                                      " Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",
                                      100,
                                      0.,
                                      1.5);
    h_invMass_[2][2] = iBooker.book1D(histname + "Endcap_FaleTracks",
                                      " Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",
                                      100,
                                      0.,
                                      1.5);

    histname = "hDPhiTracksAtVtx";
    h_DPhiTracksAtVtx_[0][0] = iBooker.book1D(histname + "All",
                                              " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",
                                              dPhiTracksBin,
                                              dPhiTracksMin,
                                              dPhiTracksMax);
    h_DPhiTracksAtVtx_[0][1] =
        iBooker.book1D(histname + "Barrel",
                       " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",
                       dPhiTracksBin,
                       dPhiTracksMin,
                       dPhiTracksMax);
    h_DPhiTracksAtVtx_[0][2] =
        iBooker.book1D(histname + "Endcap",
                       " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",
                       dPhiTracksBin,
                       dPhiTracksMin,
                       dPhiTracksMax);
    h_DPhiTracksAtVtx_[1][0] = iBooker.book1D(histname + "All_Ass",
                                              " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",
                                              dPhiTracksBin,
                                              dPhiTracksMin,
                                              dPhiTracksMax);
    h_DPhiTracksAtVtx_[1][1] =
        iBooker.book1D(histname + "Barrel_Ass",
                       " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",
                       dPhiTracksBin,
                       dPhiTracksMin,
                       dPhiTracksMax);
    h_DPhiTracksAtVtx_[1][2] =
        iBooker.book1D(histname + "Endcap_Ass",
                       " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",
                       dPhiTracksBin,
                       dPhiTracksMin,
                       dPhiTracksMax);
    h_DPhiTracksAtVtx_[2][0] = iBooker.book1D(histname + "All_Fakes",
                                              " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",
                                              dPhiTracksBin,
                                              dPhiTracksMin,
                                              dPhiTracksMax);
    h_DPhiTracksAtVtx_[2][1] =
        iBooker.book1D(histname + "Barrel_Fakes",
                       " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",
                       dPhiTracksBin,
                       dPhiTracksMin,
                       dPhiTracksMax);
    h_DPhiTracksAtVtx_[2][2] =
        iBooker.book1D(histname + "Endcap_Fakes",
                       " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",
                       dPhiTracksBin,
                       dPhiTracksMin,
                       dPhiTracksMax);

    histname = "hDPhiTracksAtVtxVsEta";
    h2_DPhiTracksAtVtxVsEta_ = iBooker.book2D(histname + "All",
                                              "  Photons:Tracks from conversions: #delta#phi Tracks at vertex vs #eta",
                                              etaBin2,
                                              etaMin,
                                              etaMax,
                                              100,
                                              -0.5,
                                              0.5);
    histname = "pDPhiTracksAtVtxVsEta";
    p_DPhiTracksAtVtxVsEta_ =
        iBooker.bookProfile(histname + "All",
                            " Photons:Tracks from conversions: #delta#phi Tracks at vertex vs #eta ",
                            etaBin2,
                            etaMin,
                            etaMax,
                            100,
                            -0.5,
                            0.5,
                            "");

    histname = "hDPhiTracksAtVtxVsR";
    h2_DPhiTracksAtVtxVsR_ = iBooker.book2D(histname + "All",
                                            "  Photons:Tracks from conversions: #delta#phi Tracks at vertex vs R",
                                            rBin,
                                            rMin,
                                            rMax,
                                            100,
                                            -0.5,
                                            0.5);
    histname = "pDPhiTracksAtVtxVsR";
    p_DPhiTracksAtVtxVsR_ = iBooker.bookProfile(histname + "All",
                                                " Photons:Tracks from conversions: #delta#phi Tracks at vertex vs R ",
                                                rBin,
                                                rMin,
                                                rMax,
                                                100,
                                                -0.5,
                                                0.5,
                                                "");

    histname = "hDCotTracks";
    h_DCotTracks_[0][0] = iBooker.book1D(histname + "All",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);
    h_DCotTracks_[0][1] = iBooker.book1D(histname + "Barrel",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);
    h_DCotTracks_[0][2] = iBooker.book1D(histname + "Endcap",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);
    h_DCotTracks_[1][0] = iBooker.book1D(histname + "All_Ass",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);
    h_DCotTracks_[1][1] = iBooker.book1D(histname + "Barrel_Ass",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);
    h_DCotTracks_[1][2] = iBooker.book1D(histname + "Endcap_Ass",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);
    h_DCotTracks_[2][0] = iBooker.book1D(histname + "All_Fakes",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);
    h_DCotTracks_[2][1] = iBooker.book1D(histname + "Barrel_Fakes",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);
    h_DCotTracks_[2][2] = iBooker.book1D(histname + "Endcap_Fakes",
                                         " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",
                                         dCotTracksBin,
                                         dCotTracksMin,
                                         dCotTracksMax);

    histname = "hDCotTracksVsEta";
    h2_DCotTracksVsEta_ = iBooker.book2D(histname + "All",
                                         "  Photons:Tracks from conversions:  #delta cotg(#Theta) Tracks vs #eta",
                                         etaBin2,
                                         etaMin,
                                         etaMax,
                                         100,
                                         -0.2,
                                         0.2);
    histname = "pDCotTracksVsEta";
    p_DCotTracksVsEta_ = iBooker.bookProfile(histname + "All",
                                             " Photons:Tracks from conversions:  #delta cotg(#Theta) Tracks vs #eta ",
                                             etaBin2,
                                             etaMin,
                                             etaMax,
                                             100,
                                             -0.2,
                                             0.2,
                                             "");

    histname = "hDCotTracksVsR";
    h2_DCotTracksVsR_ = iBooker.book2D(histname + "All",
                                       "  Photons:Tracks from conversions:  #delta cotg(#Theta)  Tracks at vertex vs R",
                                       rBin,
                                       rMin,
                                       rMax,
                                       100,
                                       -0.2,
                                       0.2);
    histname = "pDCotTracksVsR";
    p_DCotTracksVsR_ =
        iBooker.bookProfile(histname + "All",
                            " Photons:Tracks from conversions:  #delta cotg(#Theta) Tracks at vertex vs R ",
                            rBin,
                            rMin,
                            rMax,
                            100,
                            -0.2,
                            0.2,
                            "");

    histname = "hDistMinAppTracks";
    h_distMinAppTracks_[0][0] = iBooker.book1D(
        histname + "All", " Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ", 120, -0.5, 1.0);
    h_distMinAppTracks_[0][1] = iBooker.book1D(
        histname + "Barrel", " Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ", 120, -0.5, 1.0);
    h_distMinAppTracks_[0][2] = iBooker.book1D(
        histname + "Endcap", " Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ", 120, -0.5, 1.0);
    h_distMinAppTracks_[1][0] = iBooker.book1D(
        histname + "All_Ass", " Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ", 120, -0.5, 1.0);
    h_distMinAppTracks_[1][1] =
        iBooker.book1D(histname + "Barrel_Ass",
                       " Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ",
                       120,
                       -0.5,
                       1.0);
    h_distMinAppTracks_[1][2] =
        iBooker.book1D(histname + "Endcap_Ass",
                       " Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ",
                       120,
                       -0.5,
                       1.0);
    h_distMinAppTracks_[2][0] = iBooker.book1D(
        histname + "All_Fakes", " Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ", 120, -0.5, 1.0);
    h_distMinAppTracks_[2][1] =
        iBooker.book1D(histname + "Barrel_Fakes",
                       " Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ",
                       120,
                       -0.5,
                       1.0);
    h_distMinAppTracks_[2][2] =
        iBooker.book1D(histname + "Endcap_Fakes",
                       " Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ",
                       120,
                       -0.5,
                       1.0);

    h_convVtxRvsZ_[0] = iBooker.book2D("convVtxRvsZAll",
                                       " Photon Reco conversion vtx position",
                                       zBinForXray,
                                       zMinForXray,
                                       zMaxForXray,
                                       rBinForXray,
                                       rMinForXray,
                                       rMaxForXray);
    h_convVtxRvsZ_[1] = iBooker.book2D("convVtxRvsZBarrel",
                                       " Photon Reco conversion vtx position",
                                       zBinForXray,
                                       zMinForXray,
                                       zMaxForXray,
                                       rBinForXray,
                                       rMinForXray,
                                       rMaxForXray);
    h_convVtxRvsZ_[2] = iBooker.book2D("convVtxRvsZEndcap",
                                       " Photon Reco conversion vtx position",
                                       zBin2ForXray,
                                       zMinForXray,
                                       zMaxForXray,
                                       rBinForXray,
                                       rMinForXray,
                                       rMaxForXray);
    h_convVtxYvsX_ = iBooker.book2D(
        "convVtxYvsXTrkBarrel", " Photon Reco conversion vtx position, (x,y) eta<1 ", 1000, -60., 60., 1000, -60., 60.);
    h_convVtxRvsZ_zoom_[0] = iBooker.book2D("convVtxRvsZBarrelZoom1",
                                            " Photon Reco conversion vtx position",
                                            zBinForXray,
                                            zMinForXray,
                                            zMaxForXray,
                                            rBinForXray,
                                            -10.,
                                            40.);
    h_convVtxRvsZ_zoom_[1] = iBooker.book2D("convVtxRvsZBarrelZoom2",
                                            " Photon Reco conversion vtx position",
                                            zBinForXray,
                                            zMinForXray,
                                            zMaxForXray,
                                            rBinForXray,
                                            -10.,
                                            20.);
    h_convVtxYvsX_zoom_[0] = iBooker.book2D("convVtxYvsXTrkBarrelZoom1",
                                            " Photon Reco conversion vtx position, (x,y) eta<1 ",
                                            100,
                                            -40.,
                                            40.,
                                            100,
                                            -40.,
                                            40.);
    h_convVtxYvsX_zoom_[1] = iBooker.book2D("convVtxYvsXTrkBarrelZoom2",
                                            " Photon Reco conversion vtx position, (x,y) eta<1 ",
                                            100,
                                            -20.,
                                            20.,
                                            100,
                                            -20.,
                                            20.);

    h_convVtxdR_ = iBooker.book1D("convVtxdR", " Photon Reco conversion vtx dR", 100, -10., 10.);
    h_convVtxdX_ = iBooker.book1D("convVtxdX", " Photon Reco conversion vtx dX", 100, -10., 10.);
    h_convVtxdY_ = iBooker.book1D("convVtxdY", " Photon Reco conversion vtx dY", 100, -10., 10.);
    h_convVtxdZ_ = iBooker.book1D("convVtxdZ", " Photon Reco conversion vtx dZ", 100, -20., 20.);

    h_convVtxdPhi_ = iBooker.book1D("convVtxdPhi", " Photon Reco conversion vtx dPhi", 100, -0.01, 0.01);
    h_convVtxdEta_ = iBooker.book1D("convVtxdEta", " Photon Reco conversion vtx dEta", 100, -0.5, 0.5);

    h_convVtxdR_barrel_ =
        iBooker.book1D("convVtxdR_barrel", " Photon Reco conversion vtx dR, |eta|<=1.2", 100, -10., 10.);
    h_convVtxdX_barrel_ =
        iBooker.book1D("convVtxdX_barrel", " Photon Reco conversion vtx dX, |eta|<=1.2", 100, -10., 10.);
    h_convVtxdY_barrel_ =
        iBooker.book1D("convVtxdY_barrel", " Photon Reco conversion vtx dY, |eta|<=1.2 ", 100, -10., 10.);
    h_convVtxdZ_barrel_ =
        iBooker.book1D("convVtxdZ_barrel", " Photon Reco conversion vtx dZ, |eta|<=1.2,", 100, -20., 20.);

    h_convVtxdR_endcap_ =
        iBooker.book1D("convVtxdR_endcap", " Photon Reco conversion vtx dR,  |eta|>1.2 ", 100, -10., 10.);
    h_convVtxdX_endcap_ =
        iBooker.book1D("convVtxdX_endcap", " Photon Reco conversion vtx dX,  |eta|>1.2", 100, -10., 10.);
    h_convVtxdY_endcap_ =
        iBooker.book1D("convVtxdY_endcap", " Photon Reco conversion vtx dY,  |eta|>1.2", 100, -10., 10.);
    h_convVtxdZ_endcap_ =
        iBooker.book1D("convVtxdZ_endcap", " Photon Reco conversion vtx dZ,  |eta|>1.2", 100, -20., 20.);

    h2_convVtxdRVsR_ = iBooker.book2D("h2ConvVtxdRVsR", " Conversion vtx dR vsR", rBin, rMin, rMax, 100, -20., 20.);
    h2_convVtxdRVsEta_ =
        iBooker.book2D("h2ConvVtxdRVsEta", "Conversion vtx dR vs Eta", etaBin2, etaMin, etaMax, 100, -20., 20.);

    p_convVtxdRVsR_ =
        iBooker.bookProfile("pConvVtxdRVsR", " Conversion vtx dR vsR", rBin, rMin, rMax, 100, -20., 20., "");
    p_convVtxdRVsEta_ =
        iBooker.bookProfile("pConvVtxdRVsEta", "Conversion vtx dR vs Eta", etaBin2, etaMin, etaMax, 100, -20., 20., "");
    p_convVtxdXVsX_ = iBooker.bookProfile("pConvVtxdXVsX", "Conversion vtx dX vs X", 120, -60, 60, 100, -20., 20., "");
    p_convVtxdYVsY_ = iBooker.bookProfile("pConvVtxdYVsY", "Conversion vtx dY vs Y", 120, -60, 60, 100, -20., 20., "");
    p_convVtxdZVsZ_ =
        iBooker.bookProfile("pConvVtxdZVsZ", "Conversion vtx dZ vs Z", zBin, zMin, zMax, 100, -20., 20., "");

    p_convVtxdZVsR_ =
        iBooker.bookProfile("pConvVtxdZVsR", "Conversion vtx dZ vs R", rBin, rMin, rMax, 100, -20., 20., "");
    p2_convVtxdRVsRZ_ = iBooker.bookProfile2D(
        "p2ConvVtxdRVsRZ", "Conversion vtx dR vs RZ", zBin, zMin, zMax, rBin, rMin, rMax, 100, 0., 20., "s");
    p2_convVtxdZVsRZ_ = iBooker.bookProfile2D(
        "p2ConvVtxdZVsRZ", "Conversion vtx dZ vs RZ", zBin, zMin, zMax, rBin, rMin, rMax, 100, 0., 20., "s");

    histname = "EoverPtracks";
    h_EoverPTracks_[0][0] =
        iBooker.book1D(histname + "All", " photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax);
    h_EoverPTracks_[0][1] =
        iBooker.book1D(histname + "Barrel", " photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
    h_EoverPTracks_[0][2] =
        iBooker.book1D(histname + "Endcap", " photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);
    h_EoverPTracks_[1][0] =
        iBooker.book1D(histname + "All_Ass", " photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax);
    h_EoverPTracks_[1][1] = iBooker.book1D(
        histname + "Barrel_Ass", " photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
    h_EoverPTracks_[1][2] = iBooker.book1D(
        histname + "Endcap_Ass", " photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);
    h_EoverPTracks_[2][0] =
        iBooker.book1D(histname + "All_Fakes", " photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax);
    h_EoverPTracks_[2][1] = iBooker.book1D(
        histname + "Barrel_Fakes", " photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
    h_EoverPTracks_[2][2] = iBooker.book1D(
        histname + "Endcap_Fakes", " photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);

    h2_convVtxRrecVsTrue_ = iBooker.book2D(
        "h2ConvVtxRrecVsTrue", "Photon Reco conversion vtx R rec vs true", rBin, rMin, rMax, rBin, rMin, rMax);

    histname = "vtxChi2Prob";
    h_vtxChi2Prob_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 100, 0., 1.);
    h_vtxChi2Prob_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 100, 0., 1.);
    h_vtxChi2Prob_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 100, 0., 1.);
    h_vtxChi2Prob_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 100, 0., 1.);
    h_vtxChi2Prob_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 100, 0., 1.);
    h_vtxChi2Prob_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 100, 0., 1.);
    h_vtxChi2Prob_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 100, 0., 1.);
    h_vtxChi2Prob_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 100, 0., 1.);
    h_vtxChi2Prob_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 100, 0., 1.);

    h_zPVFromTracks_[1] = iBooker.book1D("zPVFromTracks", " Photons: PV z from conversion tracks", 100, -25., 25.);
    h_dzPVFromTracks_[1] =
        iBooker.book1D("dzPVFromTracks", " Photons: PV Z_rec - Z_true from conversion tracks", 100, -5., 5.);
    h2_dzPVVsR_ = iBooker.book2D("h2dzPVVsR", "Photon Reco conversions: dz(PV) vs R", rBin, rMin, rMax, 100, -3., 3.);
    p_dzPVVsR_ =
        iBooker.bookProfile("pdzPVVsR", "Photon Reco conversions: dz(PV) vs R", rBin, rMin, rMax, 100, -3., 3., "");

    histname = "lxybs";
    h_lxybs_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 200, -100., 100.);
    h_lxybs_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 200, -100., 100.);
    h_lxybs_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 200, -100., 100.);
    h_lxybs_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 200, -100., 100.);
    h_lxybs_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 200, -100., 100.);
    h_lxybs_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 200, -100., 100.);
    h_lxybs_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 200, -100., 100.);
    h_lxybs_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 200, -100., 100.);
    h_lxybs_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 200, -100., 100.);

    histname = "maxNHitsBeforeVtx";
    h_maxNHitsBeforeVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_maxNHitsBeforeVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_maxNHitsBeforeVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_maxNHitsBeforeVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_maxNHitsBeforeVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_maxNHitsBeforeVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_maxNHitsBeforeVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_maxNHitsBeforeVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_maxNHitsBeforeVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);

    histname = "leadNHitsBeforeVtx";
    h_leadNHitsBeforeVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_leadNHitsBeforeVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_leadNHitsBeforeVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_leadNHitsBeforeVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_leadNHitsBeforeVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_leadNHitsBeforeVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_leadNHitsBeforeVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_leadNHitsBeforeVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_leadNHitsBeforeVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);

    histname = "trailNHitsBeforeVtx";
    h_trailNHitsBeforeVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_trailNHitsBeforeVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_trailNHitsBeforeVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_trailNHitsBeforeVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_trailNHitsBeforeVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_trailNHitsBeforeVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_trailNHitsBeforeVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_trailNHitsBeforeVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_trailNHitsBeforeVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);

    histname = "sumNHitsBeforeVtx";
    h_sumNHitsBeforeVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_sumNHitsBeforeVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_sumNHitsBeforeVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_sumNHitsBeforeVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_sumNHitsBeforeVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_sumNHitsBeforeVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_sumNHitsBeforeVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_sumNHitsBeforeVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_sumNHitsBeforeVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);

    histname = "maxDlClosestHitToVtx";
    h_maxDlClosestHitToVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 100, -10., 10.);
    h_maxDlClosestHitToVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 100, -10., 10.);
    h_maxDlClosestHitToVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 100, -10., 10.);
    h_maxDlClosestHitToVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 100, -10., 10.);
    h_maxDlClosestHitToVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 100, -10., 10.);
    h_maxDlClosestHitToVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 100, -10., 10.);
    h_maxDlClosestHitToVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 100, -10., 10.);
    h_maxDlClosestHitToVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 100, -10., 10.);
    h_maxDlClosestHitToVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 100, -10., 10.);

    histname = "maxDlClosestHitToVtxSig";
    h_maxDlClosestHitToVtxSig_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 100, -8., 8.);
    h_maxDlClosestHitToVtxSig_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 100, -8., 8.);
    h_maxDlClosestHitToVtxSig_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 100, -8., 8.);
    h_maxDlClosestHitToVtxSig_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 100, -8., 8.);
    h_maxDlClosestHitToVtxSig_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 100, -8., 8.);
    h_maxDlClosestHitToVtxSig_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 100, -8., 8.);
    h_maxDlClosestHitToVtxSig_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 100, -8., 8.);
    h_maxDlClosestHitToVtxSig_[2][1] =
        iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 100, -8., 8.);
    h_maxDlClosestHitToVtxSig_[2][2] =
        iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 100, -8., 8.);

    histname = "deltaExpectedHitsInner";
    h_deltaExpectedHitsInner_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 31, -15.5, 15.5);
    h_deltaExpectedHitsInner_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 31, -15.5, 15.5);
    h_deltaExpectedHitsInner_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 31, -15.5, 15.5);
    h_deltaExpectedHitsInner_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 31, -15.5, 15.5);
    h_deltaExpectedHitsInner_[1][1] =
        iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 31, -15.5, 15.5);
    h_deltaExpectedHitsInner_[1][2] =
        iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 31, -15.5, 15.5);
    h_deltaExpectedHitsInner_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 31, -15.5, 15.5);
    h_deltaExpectedHitsInner_[2][1] =
        iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 31, -15.5, 15.5);
    h_deltaExpectedHitsInner_[2][2] =
        iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 31, -15.5, 15.5);

    histname = "leadExpectedHitsInner";
    h_leadExpectedHitsInner_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_leadExpectedHitsInner_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_leadExpectedHitsInner_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_leadExpectedHitsInner_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_leadExpectedHitsInner_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_leadExpectedHitsInner_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_leadExpectedHitsInner_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_leadExpectedHitsInner_[2][1] =
        iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_leadExpectedHitsInner_[2][2] =
        iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);

    histname = "nSharedHits";
    h_nSharedHits_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_nSharedHits_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_nSharedHits_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_nSharedHits_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_nSharedHits_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_nSharedHits_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
    h_nSharedHits_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
    h_nSharedHits_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
    h_nSharedHits_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);

    histname = "nHits";
    nHits_[0] = iBooker.book2D(histname + "AllTracks",
                               "Photons:Tracks from conversions: # of hits all tracks",
                               etaBin,
                               etaMin,
                               etaMax,
                               30,
                               0.,
                               30.);
    nHits_[1] = iBooker.book2D(histname + "AllTracks_Ass",
                               "Photons:Tracks from conversions: # of hits  all tracks ass",
                               etaBin,
                               etaMin,
                               etaMax,
                               30,
                               0.,
                               30.);
    nHits_[2] = iBooker.book2D(histname + "AllTracks_Fakes",
                               "Photons:Tracks from conversions: # of hits all tracks fakes",
                               etaBin,
                               etaMin,
                               etaMax,
                               30,
                               0.,
                               30.);

    histname = "nHitsVsEta";
    nHitsVsEta_[0] = iBooker.book2D(histname + "AllTracks",
                                    "Photons:Tracks from conversions: # of hits vs #eta all tracks",
                                    etaBin,
                                    etaMin,
                                    etaMax,
                                    30,
                                    0.,
                                    30.);
    nHitsVsEta_[1] = iBooker.book2D(histname + "AllTracks_Ass",
                                    "Photons:Tracks from conversions: # of hits vs #eta all tracks",
                                    etaBin,
                                    etaMin,
                                    etaMax,
                                    30,
                                    0.,
                                    30.);
    nHitsVsEta_[2] = iBooker.book2D(histname + "AllTracks_Fakes",
                                    "Photons:Tracks from conversions: # of hits vs #eta all tracks",
                                    etaBin,
                                    etaMin,
                                    etaMax,
                                    30,
                                    0.,
                                    30.);
    histname = "h_nHitsVsEta";
    p_nHitsVsEta_[0] = iBooker.bookProfile(histname + "AllTracks",
                                           "Photons:Tracks from conversions: # of hits vs #eta all tracks",
                                           etaBin,
                                           etaMin,
                                           etaMax,
                                           30,
                                           -0.5,
                                           29.5,
                                           "");
    p_nHitsVsEta_[1] = iBooker.bookProfile(histname + "AllTracks_Ass",
                                           "Photons:Tracks from conversions: # of hits vs #eta all tracks",
                                           etaBin,
                                           etaMin,
                                           etaMax,
                                           30,
                                           -0.5,
                                           29.5,
                                           "");
    p_nHitsVsEta_[2] = iBooker.bookProfile(histname + "AllTracks_Fakes",
                                           "Photons:Tracks from conversions: # of hits vs #eta all tracks",
                                           etaBin,
                                           etaMin,
                                           etaMax,
                                           30,
                                           -0.5,
                                           29.5,
                                           "");

    histname = "nHitsVsR";
    nHitsVsR_[0] = iBooker.book2D(histname + "AllTracks",
                                  "Photons:Tracks from conversions: # of hits vs radius all tracks",
                                  rBin,
                                  rMin,
                                  rMax,
                                  30,
                                  0.,
                                  30.);
    nHitsVsR_[1] = iBooker.book2D(histname + "AllTracks_Ass",
                                  "Photons:Tracks from conversions: # of hits vs radius all tracks",
                                  rBin,
                                  rMin,
                                  rMax,
                                  30,
                                  0.,
                                  30.);
    nHitsVsR_[2] = iBooker.book2D(histname + "AllTracks_Fakes",
                                  "Photons:Tracks from conversions: # of hits vs radius all tracks",
                                  rBin,
                                  rMin,
                                  rMax,
                                  30,
                                  0.,
                                  30.);

    histname = "h_nHitsVsR";
    p_nHitsVsR_[0] = iBooker.bookProfile(histname + "AllTracks",
                                         "Photons:Tracks from conversions: # of hits vs radius all tracks",
                                         rBin,
                                         rMin,
                                         rMax,
                                         30,
                                         -0.5,
                                         29.5,
                                         "");
    p_nHitsVsR_[1] = iBooker.bookProfile(histname + "AllTracks_Ass",
                                         "Photons:Tracks from conversions: # of hits vs radius all tracks",
                                         rBin,
                                         rMin,
                                         rMax,
                                         30,
                                         -0.5,
                                         29.5,
                                         "");
    p_nHitsVsR_[2] = iBooker.bookProfile(histname + "AllTracks_Fakes",
                                         "Photons:Tracks from conversions: # of hits vs radius all tracks",
                                         rBin,
                                         rMin,
                                         rMax,
                                         30,
                                         -0.5,
                                         29.5,
                                         "");

    histname = "tkChi2";
    h_tkChi2_[0] = iBooker.book1D(
        histname + "AllTracks", "Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);
    h_tkChi2_[1] = iBooker.book1D(
        histname + "AllTracks_Ass", "Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);
    h_tkChi2_[2] = iBooker.book1D(
        histname + "AllTracks_Fakes", "Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);

    histname = "tkChi2Large";
    h_tkChi2Large_[0] = iBooker.book1D(
        histname + "AllTracks", "Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);
    h_tkChi2Large_[1] = iBooker.book1D(
        histname + "AllTracks_Ass", "Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);
    h_tkChi2Large_[2] = iBooker.book1D(
        histname + "AllTracks_Fakes", "Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);

    histname = "h2Chi2VsEta";
    h2_Chi2VsEta_[0] = iBooker.book2D(
        histname + "All", " Reco Track  #chi^{2} vs #eta: All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max);
    h2_Chi2VsEta_[1] = iBooker.book2D(
        histname + "All_Ass", " Reco Track  #chi^{2} vs #eta: All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max);
    h2_Chi2VsEta_[2] = iBooker.book2D(
        histname + "All_Fakes", " Reco Track  #chi^{2} vs #eta: All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max);
    histname = "pChi2VsEta";
    p_Chi2VsEta_[0] = iBooker.bookProfile(
        histname + "All", " Reco Track #chi^{2} vs #eta : All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max, "");
    p_Chi2VsEta_[1] = iBooker.bookProfile(
        histname + "All_Ass", " Reco Track #chi^{2} vs #eta : All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max, "");
    p_Chi2VsEta_[2] = iBooker.bookProfile(histname + "All_Fakes",
                                          " Reco Track #chi^{2} vs #eta : All ",
                                          etaBin2,
                                          etaMin,
                                          etaMax,
                                          100,
                                          chi2Min,
                                          chi2Max,
                                          "");

    histname = "h2Chi2VsR";
    h2_Chi2VsR_[0] =
        iBooker.book2D(histname + "All", " Reco Track  #chi^{2} vs R: All ", rBin, rMin, rMax, 100, chi2Min, chi2Max);
    h2_Chi2VsR_[1] = iBooker.book2D(
        histname + "All_Ass", " Reco Track  #chi^{2} vs R: All ", rBin, rMin, rMax, 100, chi2Min, chi2Max);
    h2_Chi2VsR_[2] = iBooker.book2D(
        histname + "All_Fakes", " Reco Track  #chi^{2} vs R: All ", rBin, rMin, rMax, 100, chi2Min, chi2Max);
    histname = "pChi2VsR";
    p_Chi2VsR_[0] = iBooker.bookProfile(
        histname + "All", " Reco Track #chi^{2} vas R : All ", rBin, rMin, rMax, 100, chi2Min, chi2Max, "");
    p_Chi2VsR_[1] = iBooker.bookProfile(
        histname + "All_Ass", " Reco Track #chi^{2} vas R : All ", rBin, rMin, rMax, 100, chi2Min, chi2Max, "");
    p_Chi2VsR_[2] = iBooker.bookProfile(
        histname + "All_Fakes", " Reco Track #chi^{2} vas R : All ", rBin, rMin, rMax, 100, chi2Min, chi2Max, "");

    histname = "hTkD0";
    h_TkD0_[0] = iBooker.book1D(histname + "All", " Reco Track D0*q: All ", 200, -0.1, 60);
    h_TkD0_[1] = iBooker.book1D(histname + "All_Ass", " Reco Track D0*q: Barrel ", 200, -0.1, 60);
    h_TkD0_[2] = iBooker.book1D(histname + "All_Fakes", " Reco Track D0*q: Endcap ", 200, -0.1, 60);

    histname = "hTkPtPull";
    h_TkPtPull_[0] = iBooker.book1D(histname + "All", " Reco Track Pt pull: All ", 100, -20., 10.);
    histname = "hTkPtPull";
    h_TkPtPull_[1] = iBooker.book1D(histname + "Barrel", " Reco Track Pt pull: Barrel ", 100, -20., 10.);
    histname = "hTkPtPull";
    h_TkPtPull_[2] = iBooker.book1D(histname + "Endcap", " Reco Track Pt pull: Endcap ", 100, -20., 10.);

    histname = "h2TkPtPullEta";
    h2_TkPtPull_[0] =
        iBooker.book2D(histname + "All", " Reco Track Pt pull: All ", etaBin2, etaMin, etaMax, 100, -20., 10.);
    histname = "pTkPtPullEta";
    p_TkPtPull_[0] = iBooker.bookProfile(
        histname + "All", " Reco Track Pt pull: All ", etaBin2, etaMin, etaMax, 100, -20., 10., " ");

    histname = "PtRecVsPtSim";
    h2_PtRecVsPtSim_[0] =
        iBooker.book2D(histname + "All", "Pt Rec vs Pt sim: All ", etBin, etMin, etMax, etBin, etMin, etMax);
    h2_PtRecVsPtSim_[1] =
        iBooker.book2D(histname + "Barrel", "Pt Rec vs Pt sim: Barrel ", etBin, etMin, etMax, etBin, etMin, etMax);
    h2_PtRecVsPtSim_[2] =
        iBooker.book2D(histname + "Endcap", "Pt Rec vs Pt sim: Endcap ", etBin, etMin, etMax, etBin, etMin, etMax);

    histname = "photonPtRecVsPtSim";
    h2_photonPtRecVsPtSim_ =
        iBooker.book2D(histname + "All", "Pt Rec vs Pt sim: All ", etBin, etMin, etMax, etBin, etMin, etMax);

    histname = "nHitsBeforeVtx";
    h_nHitsBeforeVtx_[0] = iBooker.book1D(histname + "All", "Pt Rec vs Pt sim: All ", 16, -0.5, 15.5);
    h_nHitsBeforeVtx_[1] = iBooker.book1D(histname + "Barrel", "Pt Rec vs Pt sim: Barrel ", 16, -0.5, 15.5);
    h_nHitsBeforeVtx_[2] = iBooker.book1D(histname + "Endcap", "Pt Rec vs Pt sim: Endcap ", 16, -0.5, 15.5);

    histname = "dlClosestHitToVtx";
    h_dlClosestHitToVtx_[0] = iBooker.book1D(histname + "All", "Pt Rec vs Pt sim: All ", 100, -10., 10.);
    h_dlClosestHitToVtx_[1] = iBooker.book1D(histname + "Barrel", "Pt Rec vs Pt sim: Barrel ", 100, -10., 10.);
    h_dlClosestHitToVtx_[2] = iBooker.book1D(histname + "Endcap", "Pt Rec vs Pt sim: Endcap ", 100, -10., 10.);

    histname = "dlClosestHitToVtxSig";
    h_dlClosestHitToVtxSig_[0] = iBooker.book1D(histname + "All", "Pt Rec vs Pt sim: All ", 100, -8., 8.);
    h_dlClosestHitToVtxSig_[1] = iBooker.book1D(histname + "Barrel", "Pt Rec vs Pt sim: Barrel ", 100, -8., 8.);
    h_dlClosestHitToVtxSig_[2] = iBooker.book1D(histname + "Endcap", "Pt Rec vs Pt sim: Endcap ", 100, -8., 8.);

    h_match_ = iBooker.book1D("h_match", " ", 3, -0.5, 2.5);

  }  // if DQM
}

void TkConvValidator::dqmBeginRun ( edm::Run const &  r, edm::EventSetup const &  theEventSetup  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 1228 of file TkConvValidator.cc.

References edm::EventSetup::get().

                                                                                     {
  //get magnetic field
  edm::LogInfo("ConvertedPhotonProducer") << " get magnetic field"
                                          << "\n";
  theEventSetup.get<IdealMagneticFieldRecord>().get(theMF_);

  thePhotonMCTruthFinder_ = new PhotonMCTruthFinder();
}

void TkConvValidator::dqmEndRun ( edm::Run const &  r, edm::EventSetup const &  es  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 1237 of file TkConvValidator.cc.

                                                                                   {
  delete thePhotonMCTruthFinder_;
}

void TkConvValidator::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 2236 of file TkConvValidator.cc.

References dbe_, reco_skim_cfg_mod::outputFileName, dqm::dqmstoreimpl::DQMStore::save(), and AlCaHLTBitMon_QueryRunRegistry::string.

                             {
  std::string outputFileName = parameters_.getParameter<std::string>("OutputFileName");
  if (!isRunCentrally_) {
    dbe_->save(outputFileName);
  }

  edm::LogInfo("TkConvValidator") << "Analyzed " << nEvt_ << "\n";
  // std::cout  << "::endJob Analyzed " << nEvt_ << " events " << " with total " << nPho_ << " Photons " << "\n";
  //  std::cout  << "TkConvValidator::endJob Analyzed " << nEvt_ << " events " << "\n";

  return;
}

float TkConvValidator::etaTransformation ( float  a, float  b  )

private

Definition at line 2299 of file TkConvValidator.cc.

References ETA, etaBarrelEndcap, dqm-mbProfile::log, PI, R_ECAL, funct::tan(), and Z_Endcap.

                                                                         {
  //---Definitions
  const float PI = 3.1415927;

  //---Definitions for ECAL
  const float R_ECAL = 136.5;
  const float Z_Endcap = 328.0;
  const float etaBarrelEndcap = 1.479;

  //---ETA correction

  float Theta = 0.0;
  float ZEcal = R_ECAL * sinh(EtaParticle) + Zvertex;

  if (ZEcal != 0.0)
    Theta = atan(R_ECAL / ZEcal);
  if (Theta < 0.0)
    Theta = Theta + PI;
  float ETA = -log(tan(0.5 * Theta));

  if (fabs(ETA) > etaBarrelEndcap) {
    float Zend = Z_Endcap;
    if (EtaParticle < 0.0)
      Zend = -Zend;
    float Zlen = Zend - Zvertex;
    float RR = Zlen / sinh(EtaParticle);
    Theta = atan(RR / Zend);
    if (Theta < 0.0)
      Theta = Theta + PI;
    ETA = -log(tan(0.5 * Theta));
  }
  //---Return the result
  return ETA;
  //---end
}

float TkConvValidator::phiNormalization ( float &  a )

private

Definition at line 2283 of file TkConvValidator.cc.

References PI, and TWOPI.

                                                  {
  //---Definitions
  const float PI = 3.1415927;
  const float TWOPI = 2.0 * PI;

  if (phi > PI) {
    phi = phi - TWOPI;
  }
  if (phi < -PI) {
    phi = phi + TWOPI;
  }

  //  cout << " Float_t PHInormalization out " << PHI << endl;
  return phi;
}

math::XYZVector TkConvValidator::recalculateMomentumAtFittedVertex ( const MagneticField &  mf, const TrackerGeometry &  trackerGeom, const edm::RefToBase< reco::Track > &  tk, const reco::Vertex &  vtx  )

private

Definition at line 2249 of file TkConvValidator.cc.

References anyDirection, Cylinder::computeRadius(), f, TrajectoryStateOnSurface::globalMomentum(), trajectoryStateTransform::innerStateOnSurface(), TrajectoryStateOnSurface::isValid(), reco::Vertex::position(), Propagator::propagate(), mps_fire::result, makeMuonMisalignmentScenario::rot, mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                          {
  math::XYZVector result;
  Surface::RotationType rot;
  auto scp = new SimpleCylinderBounds(sqrt(vtx.position().perp2()) - 0.001f,
                                      sqrt(vtx.position().perp2()) + 0.001f,
                                      -fabs(vtx.position().z()),
                                      fabs(vtx.position().z()));
  ReferenceCountingPointer<Cylinder> theBarrel_(
      new Cylinder(Cylinder::computeRadius(*scp), Surface::PositionType(0, 0, 0), rot, scp));

  ReferenceCountingPointer<Disk> theDisk_(
      new Disk(Surface::PositionType(0, 0, vtx.position().z()),
               rot,
               new SimpleDiskBounds(0, sqrt(vtx.position().perp2()), -0.001, 0.001)));

  const TrajectoryStateOnSurface myTSOS = trajectoryStateTransform::innerStateOnSurface(*tk, trackerGeom, &mf);
  PropagatorWithMaterial propag(anyDirection, 0.000511, &mf);
  TrajectoryStateOnSurface stateAtVtx;
  stateAtVtx = propag.propagate(myTSOS, *theBarrel_);
  if (!stateAtVtx.isValid()) {
    stateAtVtx = propag.propagate(myTSOS, *theDisk_);
  }
  if (stateAtVtx.isValid()) {
    return math::XYZVector(double(stateAtVtx.globalMomentum().x()),
                           double(stateAtVtx.globalMomentum().y()),
                           double(stateAtVtx.globalMomentum().z()));
  } else {
    return math::XYZVector(0., 0., 0.);
  }
}

Member Data Documentation

bool TkConvValidator::arbitratedEcalSeeded_

private

Definition at line 137 of file TkConvValidator.h.

bool TkConvValidator::arbitratedMerged_

private

Definition at line 136 of file TkConvValidator.h.

double TkConvValidator::bcEtLow_

private

Definition at line 124 of file TkConvValidator.h.

edm::EDGetTokenT<reco::BeamSpot> TkConvValidator::beamspotToken_

private

Definition at line 101 of file TkConvValidator.h.

std::string TkConvValidator::conversionCollection_

private

Definition at line 91 of file TkConvValidator.h.

edm::EDGetTokenT<reco::ConversionCollection> TkConvValidator::conversionCollectionPr_Token_

private

Definition at line 92 of file TkConvValidator.h.

std::string TkConvValidator::conversionCollectionProducer_

private

Definition at line 90 of file TkConvValidator.h.

std::string TkConvValidator::conversionTrackProducer_

private

Definition at line 94 of file TkConvValidator.h.

DQMStore* TkConvValidator::dbe_

private

Definition at line 72 of file TkConvValidator.h.

bool TkConvValidator::dCotCutOn_

private

Definition at line 132 of file TkConvValidator.h.

double TkConvValidator::dCotCutValue_

private

Definition at line 133 of file TkConvValidator.h.

double TkConvValidator::dCotHardCutValue_

private

Definition at line 134 of file TkConvValidator.h.

std::string TkConvValidator::dqmpath_

private

Definition at line 110 of file TkConvValidator.h.

bool TkConvValidator::ecalalgotracks_

private

Definition at line 138 of file TkConvValidator.h.

double TkConvValidator::ecalEtSumCut_

private

Definition at line 130 of file TkConvValidator.h.

double TkConvValidator::ecalIsolRadius_

private

Definition at line 123 of file TkConvValidator.h.

std::string TkConvValidator::fName_

private

Definition at line 71 of file TkConvValidator.h.

edm::EDGetTokenT<edm::SimTrackContainer> TkConvValidator::g4_simTk_Token_

private

Definition at line 102 of file TkConvValidator.h.

edm::EDGetTokenT<edm::SimVertexContainer> TkConvValidator::g4_simVtx_Token_

private

Definition at line 103 of file TkConvValidator.h.

bool TkConvValidator::generalTracksOnly_

private

Definition at line 135 of file TkConvValidator.h.

edm::EDGetTokenT<reco::GenJetCollection> TkConvValidator::genjets_Token_

private

Definition at line 107 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_Chi2VsEta_[3]

private

Definition at line 317 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_Chi2VsR_[3]

private

Definition at line 319 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_convVtxdRVsEta_

private

Definition at line 276 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_convVtxdRVsR_

private

Definition at line 275 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_convVtxRrecVsTrue_

private

Definition at line 289 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_DCotTracksVsEta_

private

Definition at line 238 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_DCotTracksVsR_

private

Definition at line 240 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_DPhiTracksAtEcalVsEta_

private

Definition at line 248 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_DPhiTracksAtEcalVsR_

private

Definition at line 246 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_DPhiTracksAtVtxVsEta_

private

Definition at line 232 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_DPhiTracksAtVtxVsR_

private

Definition at line 234 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_dzPVVsR_

private

Definition at line 295 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_photonPtRecVsPtSim_

private

Definition at line 328 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_PtRecVsPtSim_[3]

private

Definition at line 327 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h2_TkPtPull_[3]

private

Definition at line 325 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_AllSimConv_[5]

private

Denominator for efficiencies.

Definition at line 188 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convEta2_[3][3]

private

Definition at line 214 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convEta_[3][3]

private

Definition at line 213 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convEtaMatchSC_[3][3]

private

Definition at line 212 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convPhi_[3][3]

private

Definition at line 215 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convPt_[3][3]

private

Definition at line 218 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convPtRes_[3]

private

Definition at line 227 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convR_[3][3]

private

Definition at line 216 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convRplot_

private

Definition at line 224 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convSCdEta_[3][3]

private

Definition at line 221 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convSCdPhi_[3][3]

private

Definition at line 222 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdEta_

private

Definition at line 262 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdPhi_

private

Definition at line 263 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdR_

private

Definition at line 261 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdR_barrel_

private

Definition at line 268 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdR_endcap_

private

Definition at line 273 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdX_

private

Definition at line 258 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdX_barrel_

private

Definition at line 265 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdX_endcap_

private

Definition at line 270 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdY_

private

Definition at line 259 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdY_barrel_

private

Definition at line 266 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdY_endcap_

private

Definition at line 271 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdZ_

private

Definition at line 260 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdZ_barrel_

private

Definition at line 267 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxdZ_endcap_

private

Definition at line 272 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxRvsZ_[3]

private

Definition at line 253 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxRvsZ_zoom_[2]

private

Definition at line 255 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxYvsX_

private

Definition at line 254 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convVtxYvsX_zoom_[2]

private

Definition at line 256 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convZ_[3][3]

private

Definition at line 217 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_convZplot_

private

Definition at line 225 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_DCotTracks_[3][3]

private

Definition at line 237 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_deltaExpectedHitsInner_[3][3]

private

Definition at line 303 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_DEtaTracksAtEcal_[3][3]

private

Definition at line 251 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_distMinAppTracks_[3][3]

private

Definition at line 243 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_dlClosestHitToVtx_[3]

private

Definition at line 335 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_dlClosestHitToVtxSig_[3]

private

Definition at line 336 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_DPhiTracksAtEcal_[3][3]

private

Definition at line 245 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_DPhiTracksAtVtx_[3][3]

private

Definition at line 231 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_dzPVFromTracks_[2]

private

Definition at line 294 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_EoverPTracks_[3][3]

private

Definition at line 219 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_invMass_[3][3]

private

Definition at line 229 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_leadExpectedHitsInner_[3][3]

private

Definition at line 304 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_leadNHitsBeforeVtx_[3][3]

private

Definition at line 300 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_lxybs_[3][3]

private

Definition at line 298 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_match_

private

Definition at line 330 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_maxDlClosestHitToVtx_[3][3]

private

Definition at line 305 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_maxDlClosestHitToVtxSig_[3][3]

private

Definition at line 306 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_maxNHitsBeforeVtx_[3][3]

private

Definition at line 299 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_nConv_[3][3]

private

info per conversion

Definition at line 211 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_nHitsBeforeVtx_[3]

private

Definition at line 334 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_nSharedHits_[3][3]

private

Definition at line 307 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_nSimConv_[2]

private

Definition at line 178 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_RecoConvTwoMTracks_[5]

private

Definition at line 208 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_RecoConvTwoTracks_[5]

private

Definition at line 206 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimConvEtaPix_[2]

private

Definition at line 179 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimConvOneMTracks_[5]

private

Definition at line 193 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimConvOneTracks_[5]

private

Numerator for efficiencies.

Definition at line 192 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimConvTwoMTracks_[5]

private

Definition at line 195 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT0005_[5]

private

Definition at line 197 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT01_[5]

private

Definition at line 198 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT0_[5]

private

Definition at line 196 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimConvTwoTracks_[5]

private

Definition at line 194 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_simConvVtxRvsZ_[4]

private

Definition at line 184 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_simConvVtxYvsX_

private

Definition at line 185 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimRecConvOneMTracks_[5]

private

Definition at line 201 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimRecConvOneTracks_[5]

private

Definition at line 200 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimRecConvTwoMTracks_[5]

private

Definition at line 203 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_SimRecConvTwoTracks_[5]

private

Definition at line 202 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_simTkEta_

private

Definition at line 182 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_simTkPt_

private

Definition at line 181 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_sumNHitsBeforeVtx_[3][3]

private

Definition at line 302 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_tkChi2_[3]

private

Definition at line 315 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_tkChi2Large_[3]

private

Definition at line 316 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_TkD0_[3]

private

Definition at line 322 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_TkPtPull_[3]

private

Definition at line 324 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_trailNHitsBeforeVtx_[3][3]

private

Definition at line 301 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_VisSimConv_[6]

private

Definition at line 189 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_VisSimConvLarge_

private

Definition at line 190 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_vtxChi2Prob_[3][3]

private

Definition at line 291 of file TkConvValidator.h.

MonitorElement* TkConvValidator::h_zPVFromTracks_[2]

private

Definition at line 293 of file TkConvValidator.h.

double TkConvValidator::hcalEtSumCut_

private

Definition at line 131 of file TkConvValidator.h.

double TkConvValidator::hcalHitEtLow_

private

Definition at line 127 of file TkConvValidator.h.

double TkConvValidator::hcalIsolExtRadius_

private

Definition at line 125 of file TkConvValidator.h.

double TkConvValidator::hcalIsolInnRadius_

private

Definition at line 126 of file TkConvValidator.h.

edm::EDGetTokenT<edm::HepMCProduct> TkConvValidator::hepMC_Token_

private

Definition at line 106 of file TkConvValidator.h.

bool TkConvValidator::highPurity_

private

Definition at line 139 of file TkConvValidator.h.

bool TkConvValidator::isRunCentrally_

private

Definition at line 116 of file TkConvValidator.h.

edm::InputTag TkConvValidator::label_tp_

private

Definition at line 112 of file TkConvValidator.h.

double TkConvValidator::lip_

private

Definition at line 122 of file TkConvValidator.h.

uint TkConvValidator::maxHitsBeforeVtx_

private

Definition at line 141 of file TkConvValidator.h.

double TkConvValidator::maxPhoEtaForEffic

private

Definition at line 161 of file TkConvValidator.h.

double TkConvValidator::maxPhoEtaForPurity

private

Definition at line 165 of file TkConvValidator.h.

double TkConvValidator::maxPhoRForEffic

private

Definition at line 163 of file TkConvValidator.h.

double TkConvValidator::maxPhoRForPurity

private

Definition at line 167 of file TkConvValidator.h.

double TkConvValidator::maxPhoZForEffic

private

Definition at line 162 of file TkConvValidator.h.

double TkConvValidator::maxPhoZForPurity

private

Definition at line 166 of file TkConvValidator.h.

double TkConvValidator::mcConvEta_

private

Definition at line 153 of file TkConvValidator.h.

double TkConvValidator::mcConvPhi_

private

Definition at line 152 of file TkConvValidator.h.

double TkConvValidator::mcConvPt_

private

Definition at line 147 of file TkConvValidator.h.

double TkConvValidator::mcConvR_

private

Definition at line 148 of file TkConvValidator.h.

double TkConvValidator::mcConvX_

private

Definition at line 151 of file TkConvValidator.h.

double TkConvValidator::mcConvY_

private

Definition at line 150 of file TkConvValidator.h.

double TkConvValidator::mcConvZ_

private

Definition at line 149 of file TkConvValidator.h.

double TkConvValidator::mcEta_

private

Definition at line 146 of file TkConvValidator.h.

double TkConvValidator::mcJetEta_

private

Definition at line 154 of file TkConvValidator.h.

double TkConvValidator::mcJetPhi_

private

Definition at line 155 of file TkConvValidator.h.

double TkConvValidator::mcPhi_

private

global variable for the MC photon

Definition at line 145 of file TkConvValidator.h.

double TkConvValidator::minLxy_

private

Definition at line 142 of file TkConvValidator.h.

double TkConvValidator::minPhoEtCut_

private

Definition at line 118 of file TkConvValidator.h.

double TkConvValidator::minPhoPtForEffic

private

Definition at line 160 of file TkConvValidator.h.

double TkConvValidator::minPhoPtForPurity

private

Definition at line 164 of file TkConvValidator.h.

double TkConvValidator::minProb_

private

Definition at line 140 of file TkConvValidator.h.

int TkConvValidator::nEntry_

private

Definition at line 77 of file TkConvValidator.h.

int TkConvValidator::nEvt_

private

Definition at line 76 of file TkConvValidator.h.

MonitorElement* TkConvValidator::nHits_[3]

private

Definition at line 310 of file TkConvValidator.h.

MonitorElement* TkConvValidator::nHitsVsEta_[3]

private

Definition at line 312 of file TkConvValidator.h.

MonitorElement* TkConvValidator::nHitsVsR_[3]

private

Definition at line 314 of file TkConvValidator.h.

int TkConvValidator::nInvalidPCA_

private

Definition at line 84 of file TkConvValidator.h.

int TkConvValidator::nMatched_

private

Definition at line 79 of file TkConvValidator.h.

int TkConvValidator::nRecConv_

private

Definition at line 80 of file TkConvValidator.h.

int TkConvValidator::nRecConvAss_

private

Definition at line 81 of file TkConvValidator.h.

int TkConvValidator::nRecConvAssWithEcal_

private

Definition at line 82 of file TkConvValidator.h.

int TkConvValidator::nSimConv_[2]

private

Definition at line 78 of file TkConvValidator.h.

int TkConvValidator::numOfTracksInCone_

private

Definition at line 128 of file TkConvValidator.h.

edm::EDGetTokenT<reco::VertexCollection> TkConvValidator::offline_pvToken_

private

Definition at line 100 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p2_convVtxdRVsRZ_

private

Definition at line 286 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p2_convVtxdZVsRZ_

private

Definition at line 287 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p2_effRZ_

private

Definition at line 332 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_Chi2VsEta_[3]

private

Definition at line 318 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_Chi2VsR_[3]

private

Definition at line 320 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_convVtxdRVsEta_

private

Definition at line 279 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_convVtxdRVsR_

private

Definition at line 278 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_convVtxdXVsX_

private

Definition at line 281 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_convVtxdYVsY_

private

Definition at line 282 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_convVtxdZVsR_

private

Definition at line 284 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_convVtxdZVsZ_

private

Definition at line 283 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_DCotTracksVsEta_

private

Definition at line 239 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_DCotTracksVsR_

private

Definition at line 241 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_DPhiTracksAtEcalVsEta_

private

Definition at line 249 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_DPhiTracksAtEcalVsR_

private

Definition at line 247 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_DPhiTracksAtVtxVsEta_

private

Definition at line 233 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_DPhiTracksAtVtxVsR_

private

Definition at line 235 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_dzPVVsR_

private

Definition at line 296 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_nHitsVsEta_[3]

private

Definition at line 311 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_nHitsVsR_[3]

private

Definition at line 313 of file TkConvValidator.h.

MonitorElement* TkConvValidator::p_TkPtPull_[3]

private

Definition at line 326 of file TkConvValidator.h.

edm::ParameterSet TkConvValidator::parameters_

private

Definition at line 86 of file TkConvValidator.h.

Referenced by Mixins._TypedParameterizable::clone(), Types.PSet::clone(), Types.EDAlias::clone(), and Mixins._TypedParameterizable::copy().

std::string TkConvValidator::photonCollection_

private

Definition at line 97 of file TkConvValidator.h.

edm::EDGetTokenT<reco::PhotonCollection> TkConvValidator::photonCollectionPr_Token_

private

Definition at line 98 of file TkConvValidator.h.

std::string TkConvValidator::photonCollectionProducer_

private

Definition at line 96 of file TkConvValidator.h.

double TkConvValidator::recMaxPt_

private

Definition at line 174 of file TkConvValidator.h.

double TkConvValidator::recMinPt_

private

Global variables for reco Photon.

Definition at line 173 of file TkConvValidator.h.

double TkConvValidator::simMaxPt_

private

Definition at line 170 of file TkConvValidator.h.

double TkConvValidator::simMinPt_

private

Definition at line 169 of file TkConvValidator.h.

edm::ESHandle<CaloGeometry> TkConvValidator::theCaloGeom_

private

Definition at line 87 of file TkConvValidator.h.

edm::ESHandle<CaloTopology> TkConvValidator::theCaloTopo_

private

Definition at line 88 of file TkConvValidator.h.

edm::RefVector<TrackingParticleCollection> TkConvValidator::theConvTP_

private

Definition at line 157 of file TkConvValidator.h.

edm::ESHandle<MagneticField> TkConvValidator::theMF_

private

Definition at line 73 of file TkConvValidator.h.

PhotonMCTruthFinder* TkConvValidator::thePhotonMCTruthFinder_

private

Definition at line 114 of file TkConvValidator.h.

edm::EDGetTokenT<TrackingParticleRefVector> TkConvValidator::tpSelForEff_Token_

private

Definition at line 104 of file TkConvValidator.h.

edm::EDGetTokenT<TrackingParticleRefVector> TkConvValidator::tpSelForFake_Token_

private

Definition at line 105 of file TkConvValidator.h.

edm::EDGetTokenT<reco::TrackToTrackingParticleAssociator> TkConvValidator::trackAssociator_Token_

private

Definition at line 108 of file TkConvValidator.h.

double TkConvValidator::trkIsolExtRadius_

private

Definition at line 119 of file TkConvValidator.h.

double TkConvValidator::trkIsolInnRadius_

private

Definition at line 120 of file TkConvValidator.h.

double TkConvValidator::trkPtLow_

private

Definition at line 121 of file TkConvValidator.h.

double TkConvValidator::trkPtSumCut_

private

Definition at line 129 of file TkConvValidator.h.

int TkConvValidator::verbosity_

private

Definition at line 75 of file TkConvValidator.h.