#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 DQMOneEDAnalyzer<>
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

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::Accumulator, Args... >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () 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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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_

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

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_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magneticFieldToken_

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_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeometryToken_

const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord >	transientTrackBuilderToken_

double	trkIsolExtRadius_

double	trkIsolInnRadius_

double	trkPtLow_

double	trkPtSumCut_

int	verbosity_

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer<>
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
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Protected Attributes inherited from DQMOneEDAnalyzer<>
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

$Id: TkConvValidator

Author: N.Marinelli - Univ. of Notre Dame

Definition at line 51 of file TkConvValidator.h.

Constructor & Destructor Documentation

◆ TkConvValidator()

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

explicit

Definition at line 95 of file TkConvValidator.cc.

     : magneticFieldToken_{esConsumes<edm::Transition::BeginRun>()},
       caloGeometryToken_{esConsumes()},
       transientTrackBuilderToken_{esConsumes(edm::ESInputTag("", "TransientTrackBuilder"))},
       trackerGeometryToken_{esConsumes()} {
   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::~TkConvValidator ( )

override

Definition at line 144 of file TkConvValidator.cc.

144 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 1237 of file TkConvValidator.cc.

                                                                           {
   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:
   theCaloGeom_ = esup.getHandle(caloGeometryToken_);
 
   // Transform Track into TransientTrack (needed by the Vertex fitter)
   auto theTTB = esup.getHandle(transientTrackBuilderToken_);
 
   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);
   auto tracker = esup.getHandle(trackerGeometryToken_);
   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());
 
   // ################  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();
 
     float 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 (p1AtVtx.perp2() > 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.0;
     //    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 (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 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_));
 }

◆ bookHistograms()

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

overridevirtual

Histograms for efficiencies

Denominators

zooms

Implements DQMOneEDAnalyzer<>.

Definition at line 146 of file TkConvValidator.cc.

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

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 1226 of file TkConvValidator.cc.

References edm::EventSetup::getHandle(), magneticFieldToken_, theMF_, and thePhotonMCTruthFinder_.

                                                                                      {
   //get magnetic field
   edm::LogInfo("ConvertedPhotonProducer") << " get magnetic field"
                                           << "\n";
   theMF_ = theEventSetup.getHandle(magneticFieldToken_);
 
   thePhotonMCTruthFinder_ = new PhotonMCTruthFinder();
 }

◆ dqmEndRun()

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

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 1235 of file TkConvValidator.cc.

References thePhotonMCTruthFinder_.

1235 { delete thePhotonMCTruthFinder_; }

TkConvValidator::thePhotonMCTruthFinder_

PhotonMCTruthFinder * thePhotonMCTruthFinder_

Definition: TkConvValidator.h:121

◆ endJob()

void TkConvValidator::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 2229 of file TkConvValidator.cc.

References dbe_, edm::ParameterSet::getParameter(), isRunCentrally_, nEvt_, reco_skim_cfg_mod::outputFileName, parameters_, dqm::implementation::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;
 }

◆ etaTransformation()

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

private

Definition at line 2292 of file TkConvValidator.cc.

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

Referenced by analyze().

                                                                          {
   //---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
 }

◆ phiNormalization()

float TkConvValidator::phiNormalization ( float & a )

private

Definition at line 2276 of file TkConvValidator.cc.

References phi, PI, and TWOPI.

Referenced by analyze().

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

◆ recalculateMomentumAtFittedVertex()

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

private

Definition at line 2242 of file TkConvValidator.cc.

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

Referenced by analyze().

                                                                                           {
   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

◆ arbitratedEcalSeeded_

bool TkConvValidator::arbitratedEcalSeeded_

private

Definition at line 144 of file TkConvValidator.h.

Referenced by analyze().

◆ arbitratedMerged_

bool TkConvValidator::arbitratedMerged_

private

Definition at line 143 of file TkConvValidator.h.

Referenced by analyze().

◆ bcEtLow_

double TkConvValidator::bcEtLow_

private

Definition at line 131 of file TkConvValidator.h.

◆ beamspotToken_

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

private

Definition at line 103 of file TkConvValidator.h.

Referenced by analyze().

◆ caloGeometryToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> TkConvValidator::caloGeometryToken_

private

Definition at line 113 of file TkConvValidator.h.

Referenced by analyze().

◆ conversionCollection_

std::string TkConvValidator::conversionCollection_

private

Definition at line 93 of file TkConvValidator.h.

◆ conversionCollectionPr_Token_

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

private

Definition at line 94 of file TkConvValidator.h.

Referenced by analyze().

◆ conversionCollectionProducer_

std::string TkConvValidator::conversionCollectionProducer_

private

Definition at line 92 of file TkConvValidator.h.

◆ conversionTrackProducer_

std::string TkConvValidator::conversionTrackProducer_

private

Definition at line 96 of file TkConvValidator.h.

◆ dbe_

DQMStore* TkConvValidator::dbe_

private

Definition at line 74 of file TkConvValidator.h.

Referenced by bookHistograms(), and endJob().

◆ dCotCutOn_

bool TkConvValidator::dCotCutOn_

private

Definition at line 139 of file TkConvValidator.h.

◆ dCotCutValue_

double TkConvValidator::dCotCutValue_

private

Definition at line 140 of file TkConvValidator.h.

◆ dCotHardCutValue_

double TkConvValidator::dCotHardCutValue_

private

Definition at line 141 of file TkConvValidator.h.

◆ dqmpath_

std::string TkConvValidator::dqmpath_

private

Definition at line 117 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ ecalalgotracks_

bool TkConvValidator::ecalalgotracks_

private

Definition at line 145 of file TkConvValidator.h.

Referenced by analyze().

◆ ecalEtSumCut_

double TkConvValidator::ecalEtSumCut_

private

Definition at line 137 of file TkConvValidator.h.

◆ ecalIsolRadius_

double TkConvValidator::ecalIsolRadius_

private

Definition at line 130 of file TkConvValidator.h.

◆ fName_

std::string TkConvValidator::fName_

private

Definition at line 73 of file TkConvValidator.h.

◆ g4_simTk_Token_

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

private

Definition at line 104 of file TkConvValidator.h.

Referenced by analyze().

◆ g4_simVtx_Token_

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

private

Definition at line 105 of file TkConvValidator.h.

Referenced by analyze().

◆ generalTracksOnly_

bool TkConvValidator::generalTracksOnly_

private

Definition at line 142 of file TkConvValidator.h.

Referenced by analyze().

◆ genjets_Token_

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

private

Definition at line 109 of file TkConvValidator.h.

◆ h2_Chi2VsEta_

MonitorElement* TkConvValidator::h2_Chi2VsEta_[3]

private

Definition at line 324 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_Chi2VsR_

MonitorElement* TkConvValidator::h2_Chi2VsR_[3]

private

Definition at line 326 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_convVtxdRVsEta_

MonitorElement* TkConvValidator::h2_convVtxdRVsEta_

private

Definition at line 283 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_convVtxdRVsR_

MonitorElement* TkConvValidator::h2_convVtxdRVsR_

private

Definition at line 282 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_convVtxRrecVsTrue_

MonitorElement* TkConvValidator::h2_convVtxRrecVsTrue_

private

Definition at line 296 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_DCotTracksVsEta_

MonitorElement* TkConvValidator::h2_DCotTracksVsEta_

private

Definition at line 245 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_DCotTracksVsR_

MonitorElement* TkConvValidator::h2_DCotTracksVsR_

private

Definition at line 247 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_DPhiTracksAtEcalVsEta_

MonitorElement* TkConvValidator::h2_DPhiTracksAtEcalVsEta_

private

Definition at line 255 of file TkConvValidator.h.

◆ h2_DPhiTracksAtEcalVsR_

MonitorElement* TkConvValidator::h2_DPhiTracksAtEcalVsR_

private

Definition at line 253 of file TkConvValidator.h.

◆ h2_DPhiTracksAtVtxVsEta_

MonitorElement* TkConvValidator::h2_DPhiTracksAtVtxVsEta_

private

Definition at line 239 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_DPhiTracksAtVtxVsR_

MonitorElement* TkConvValidator::h2_DPhiTracksAtVtxVsR_

private

Definition at line 241 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_dzPVVsR_

MonitorElement* TkConvValidator::h2_dzPVVsR_

private

Definition at line 302 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_photonPtRecVsPtSim_

MonitorElement* TkConvValidator::h2_photonPtRecVsPtSim_

private

Definition at line 335 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_PtRecVsPtSim_

MonitorElement* TkConvValidator::h2_PtRecVsPtSim_[3]

private

Definition at line 334 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h2_TkPtPull_

MonitorElement* TkConvValidator::h2_TkPtPull_[3]

private

Definition at line 332 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_AllSimConv_

MonitorElement* TkConvValidator::h_AllSimConv_[5]

private

Denominator for efficiencies.

Definition at line 195 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convEta2_

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

private

Definition at line 221 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convEta_

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

private

Definition at line 220 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convEtaMatchSC_

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

private

Definition at line 219 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convPhi_

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

private

Definition at line 222 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convPt_

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

private

Definition at line 225 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convPtRes_

MonitorElement* TkConvValidator::h_convPtRes_[3]

private

Definition at line 234 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convR_

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

private

Definition at line 223 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convRplot_

MonitorElement* TkConvValidator::h_convRplot_

private

Definition at line 231 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convSCdEta_

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

private

Definition at line 228 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convSCdPhi_

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

private

Definition at line 229 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdEta_

MonitorElement* TkConvValidator::h_convVtxdEta_

private

Definition at line 269 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdPhi_

MonitorElement* TkConvValidator::h_convVtxdPhi_

private

Definition at line 270 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdR_

MonitorElement* TkConvValidator::h_convVtxdR_

private

Definition at line 268 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdR_barrel_

MonitorElement* TkConvValidator::h_convVtxdR_barrel_

private

Definition at line 275 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdR_endcap_

MonitorElement* TkConvValidator::h_convVtxdR_endcap_

private

Definition at line 280 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdX_

MonitorElement* TkConvValidator::h_convVtxdX_

private

Definition at line 265 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdX_barrel_

MonitorElement* TkConvValidator::h_convVtxdX_barrel_

private

Definition at line 272 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdX_endcap_

MonitorElement* TkConvValidator::h_convVtxdX_endcap_

private

Definition at line 277 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdY_

MonitorElement* TkConvValidator::h_convVtxdY_

private

Definition at line 266 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdY_barrel_

MonitorElement* TkConvValidator::h_convVtxdY_barrel_

private

Definition at line 273 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdY_endcap_

MonitorElement* TkConvValidator::h_convVtxdY_endcap_

private

Definition at line 278 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdZ_

MonitorElement* TkConvValidator::h_convVtxdZ_

private

Definition at line 267 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdZ_barrel_

MonitorElement* TkConvValidator::h_convVtxdZ_barrel_

private

Definition at line 274 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxdZ_endcap_

MonitorElement* TkConvValidator::h_convVtxdZ_endcap_

private

Definition at line 279 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxRvsZ_

MonitorElement* TkConvValidator::h_convVtxRvsZ_[3]

private

Definition at line 260 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxRvsZ_zoom_

MonitorElement* TkConvValidator::h_convVtxRvsZ_zoom_[2]

private

Definition at line 262 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ h_convVtxYvsX_

MonitorElement* TkConvValidator::h_convVtxYvsX_

private

Definition at line 261 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convVtxYvsX_zoom_

MonitorElement* TkConvValidator::h_convVtxYvsX_zoom_[2]

private

Definition at line 263 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convZ_

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

private

Definition at line 224 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_convZplot_

MonitorElement* TkConvValidator::h_convZplot_

private

Definition at line 232 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_DCotTracks_

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

private

Definition at line 244 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_deltaExpectedHitsInner_

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

private

Definition at line 310 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_DEtaTracksAtEcal_

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

private

Definition at line 258 of file TkConvValidator.h.

◆ h_distMinAppTracks_

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

private

Definition at line 250 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_dlClosestHitToVtx_

MonitorElement* TkConvValidator::h_dlClosestHitToVtx_[3]

private

Definition at line 342 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_dlClosestHitToVtxSig_

MonitorElement* TkConvValidator::h_dlClosestHitToVtxSig_[3]

private

Definition at line 343 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_DPhiTracksAtEcal_

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

private

Definition at line 252 of file TkConvValidator.h.

◆ h_DPhiTracksAtVtx_

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

private

Definition at line 238 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_dzPVFromTracks_

MonitorElement* TkConvValidator::h_dzPVFromTracks_[2]

private

Definition at line 301 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_EoverPTracks_

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

private

Definition at line 226 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_invMass_

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

private

Definition at line 236 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_leadExpectedHitsInner_

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

private

Definition at line 311 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_leadNHitsBeforeVtx_

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

private

Definition at line 307 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_lxybs_

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

private

Definition at line 305 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_match_

MonitorElement* TkConvValidator::h_match_

private

Definition at line 337 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_maxDlClosestHitToVtx_

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

private

Definition at line 312 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_maxDlClosestHitToVtxSig_

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

private

Definition at line 313 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_maxNHitsBeforeVtx_

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

private

Definition at line 306 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_nConv_

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

private

info per conversion

Definition at line 218 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_nHitsBeforeVtx_

MonitorElement* TkConvValidator::h_nHitsBeforeVtx_[3]

private

Definition at line 341 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_nSharedHits_

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

private

Definition at line 314 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_nSimConv_

MonitorElement* TkConvValidator::h_nSimConv_[2]

private

Definition at line 185 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ h_RecoConvTwoMTracks_

MonitorElement* TkConvValidator::h_RecoConvTwoMTracks_[5]

private

Definition at line 215 of file TkConvValidator.h.

◆ h_RecoConvTwoTracks_

MonitorElement* TkConvValidator::h_RecoConvTwoTracks_[5]

private

Definition at line 213 of file TkConvValidator.h.

◆ h_SimConvEtaPix_

MonitorElement* TkConvValidator::h_SimConvEtaPix_[2]

private

Definition at line 186 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_SimConvOneMTracks_

MonitorElement* TkConvValidator::h_SimConvOneMTracks_[5]

private

Definition at line 200 of file TkConvValidator.h.

◆ h_SimConvOneTracks_

MonitorElement* TkConvValidator::h_SimConvOneTracks_[5]

private

Numerator for efficiencies.

Definition at line 199 of file TkConvValidator.h.

◆ h_SimConvTwoMTracks_

MonitorElement* TkConvValidator::h_SimConvTwoMTracks_[5]

private

Definition at line 202 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_SimConvTwoMTracksAndVtxPGT0005_

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT0005_[5]

private

Definition at line 204 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_SimConvTwoMTracksAndVtxPGT01_

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT01_[5]

private

Definition at line 205 of file TkConvValidator.h.

◆ h_SimConvTwoMTracksAndVtxPGT0_

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT0_[5]

private

Definition at line 203 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_SimConvTwoTracks_

MonitorElement* TkConvValidator::h_SimConvTwoTracks_[5]

private

Definition at line 201 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ h_simConvVtxRvsZ_

MonitorElement* TkConvValidator::h_simConvVtxRvsZ_[4]

private

Definition at line 191 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_simConvVtxYvsX_

MonitorElement* TkConvValidator::h_simConvVtxYvsX_

private

Definition at line 192 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_SimRecConvOneMTracks_

MonitorElement* TkConvValidator::h_SimRecConvOneMTracks_[5]

private

Definition at line 208 of file TkConvValidator.h.

◆ h_SimRecConvOneTracks_

MonitorElement* TkConvValidator::h_SimRecConvOneTracks_[5]

private

Definition at line 207 of file TkConvValidator.h.

◆ h_SimRecConvTwoMTracks_

MonitorElement* TkConvValidator::h_SimRecConvTwoMTracks_[5]

private

Definition at line 210 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_SimRecConvTwoTracks_

MonitorElement* TkConvValidator::h_SimRecConvTwoTracks_[5]

private

Definition at line 209 of file TkConvValidator.h.

◆ h_simTkEta_

MonitorElement* TkConvValidator::h_simTkEta_

private

Definition at line 189 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_simTkPt_

MonitorElement* TkConvValidator::h_simTkPt_

private

Definition at line 188 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_sumNHitsBeforeVtx_

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

private

Definition at line 309 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_tkChi2_

MonitorElement* TkConvValidator::h_tkChi2_[3]

private

Definition at line 322 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_tkChi2Large_

MonitorElement* TkConvValidator::h_tkChi2Large_[3]

private

Definition at line 323 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_TkD0_

MonitorElement* TkConvValidator::h_TkD0_[3]

private

Definition at line 329 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_TkPtPull_

MonitorElement* TkConvValidator::h_TkPtPull_[3]

private

Definition at line 331 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_trailNHitsBeforeVtx_

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

private

Definition at line 308 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_VisSimConv_

MonitorElement* TkConvValidator::h_VisSimConv_[6]

private

Definition at line 196 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_VisSimConvLarge_

MonitorElement* TkConvValidator::h_VisSimConvLarge_

private

Definition at line 197 of file TkConvValidator.h.

◆ h_vtxChi2Prob_

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

private

Definition at line 298 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ h_zPVFromTracks_

MonitorElement* TkConvValidator::h_zPVFromTracks_[2]

private

Definition at line 300 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ hcalEtSumCut_

double TkConvValidator::hcalEtSumCut_

private

Definition at line 138 of file TkConvValidator.h.

◆ hcalHitEtLow_

double TkConvValidator::hcalHitEtLow_

private

Definition at line 134 of file TkConvValidator.h.

◆ hcalIsolExtRadius_

double TkConvValidator::hcalIsolExtRadius_

private

Definition at line 132 of file TkConvValidator.h.

◆ hcalIsolInnRadius_

double TkConvValidator::hcalIsolInnRadius_

private

Definition at line 133 of file TkConvValidator.h.

◆ hepMC_Token_

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

private

Definition at line 108 of file TkConvValidator.h.

◆ highPurity_

bool TkConvValidator::highPurity_

private

Definition at line 146 of file TkConvValidator.h.

Referenced by analyze().

◆ isRunCentrally_

bool TkConvValidator::isRunCentrally_

private

Definition at line 123 of file TkConvValidator.h.

Referenced by endJob().

◆ label_tp_

edm::InputTag TkConvValidator::label_tp_

private

Definition at line 119 of file TkConvValidator.h.

◆ lip_

double TkConvValidator::lip_

private

Definition at line 129 of file TkConvValidator.h.

◆ magneticFieldToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> TkConvValidator::magneticFieldToken_

private

Definition at line 112 of file TkConvValidator.h.

Referenced by dqmBeginRun().

◆ maxHitsBeforeVtx_

uint TkConvValidator::maxHitsBeforeVtx_

private

Definition at line 148 of file TkConvValidator.h.

Referenced by analyze().

◆ maxPhoEtaForEffic

double TkConvValidator::maxPhoEtaForEffic

private

Definition at line 168 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ maxPhoEtaForPurity

double TkConvValidator::maxPhoEtaForPurity

private

Definition at line 172 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ maxPhoRForEffic

double TkConvValidator::maxPhoRForEffic

private

Definition at line 170 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ maxPhoRForPurity

double TkConvValidator::maxPhoRForPurity

private

Definition at line 174 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ maxPhoZForEffic

double TkConvValidator::maxPhoZForEffic

private

Definition at line 169 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ maxPhoZForPurity

double TkConvValidator::maxPhoZForPurity

private

Definition at line 173 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ mcConvEta_

double TkConvValidator::mcConvEta_

private

Definition at line 160 of file TkConvValidator.h.

Referenced by analyze().

◆ mcConvPhi_

double TkConvValidator::mcConvPhi_

private

Definition at line 159 of file TkConvValidator.h.

Referenced by analyze().

◆ mcConvPt_

double TkConvValidator::mcConvPt_

private

Definition at line 154 of file TkConvValidator.h.

Referenced by analyze().

◆ mcConvR_

double TkConvValidator::mcConvR_

private

Definition at line 155 of file TkConvValidator.h.

Referenced by analyze().

◆ mcConvX_

double TkConvValidator::mcConvX_

private

Definition at line 158 of file TkConvValidator.h.

Referenced by analyze().

◆ mcConvY_

double TkConvValidator::mcConvY_

private

Definition at line 157 of file TkConvValidator.h.

Referenced by analyze().

◆ mcConvZ_

double TkConvValidator::mcConvZ_

private

Definition at line 156 of file TkConvValidator.h.

Referenced by analyze().

◆ mcEta_

double TkConvValidator::mcEta_

private

Definition at line 153 of file TkConvValidator.h.

Referenced by analyze().

◆ mcJetEta_

double TkConvValidator::mcJetEta_

private

Definition at line 161 of file TkConvValidator.h.

◆ mcJetPhi_

double TkConvValidator::mcJetPhi_

private

Definition at line 162 of file TkConvValidator.h.

◆ mcPhi_

double TkConvValidator::mcPhi_

private

global variable for the MC photon

Definition at line 152 of file TkConvValidator.h.

Referenced by analyze().

◆ minLxy_

double TkConvValidator::minLxy_

private

Definition at line 149 of file TkConvValidator.h.

Referenced by analyze().

◆ minPhoEtCut_

double TkConvValidator::minPhoEtCut_

private

Definition at line 125 of file TkConvValidator.h.

◆ minPhoPtForEffic

double TkConvValidator::minPhoPtForEffic

private

Definition at line 167 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ minPhoPtForPurity

double TkConvValidator::minPhoPtForPurity

private

Definition at line 171 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ minProb_

double TkConvValidator::minProb_

private

Definition at line 147 of file TkConvValidator.h.

Referenced by analyze().

◆ nEntry_

int TkConvValidator::nEntry_

private

Definition at line 79 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ nEvt_

int TkConvValidator::nEvt_

private

Definition at line 78 of file TkConvValidator.h.

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

◆ nHits_

MonitorElement* TkConvValidator::nHits_[3]

private

Definition at line 317 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ nHitsVsEta_

MonitorElement* TkConvValidator::nHitsVsEta_[3]

private

Definition at line 319 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ nHitsVsR_

MonitorElement* TkConvValidator::nHitsVsR_[3]

private

Definition at line 321 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ nInvalidPCA_

int TkConvValidator::nInvalidPCA_

private

Definition at line 86 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ nMatched_

int TkConvValidator::nMatched_

private

Definition at line 81 of file TkConvValidator.h.

◆ nRecConv_

int TkConvValidator::nRecConv_

private

Definition at line 82 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ nRecConvAss_

int TkConvValidator::nRecConvAss_

private

Definition at line 83 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ nRecConvAssWithEcal_

int TkConvValidator::nRecConvAssWithEcal_

private

Definition at line 84 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ nSimConv_

int TkConvValidator::nSimConv_[2]

private

Definition at line 80 of file TkConvValidator.h.

Referenced by analyze().

◆ numOfTracksInCone_

int TkConvValidator::numOfTracksInCone_

private

Definition at line 135 of file TkConvValidator.h.

◆ offline_pvToken_

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

private

Definition at line 102 of file TkConvValidator.h.

Referenced by analyze().

◆ p2_convVtxdRVsRZ_

MonitorElement* TkConvValidator::p2_convVtxdRVsRZ_

private

Definition at line 293 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p2_convVtxdZVsRZ_

MonitorElement* TkConvValidator::p2_convVtxdZVsRZ_

private

Definition at line 294 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p2_effRZ_

MonitorElement* TkConvValidator::p2_effRZ_

private

Definition at line 339 of file TkConvValidator.h.

◆ p_Chi2VsEta_

MonitorElement* TkConvValidator::p_Chi2VsEta_[3]

private

Definition at line 325 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_Chi2VsR_

MonitorElement* TkConvValidator::p_Chi2VsR_[3]

private

Definition at line 327 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_convVtxdRVsEta_

MonitorElement* TkConvValidator::p_convVtxdRVsEta_

private

Definition at line 286 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_convVtxdRVsR_

MonitorElement* TkConvValidator::p_convVtxdRVsR_

private

Definition at line 285 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_convVtxdXVsX_

MonitorElement* TkConvValidator::p_convVtxdXVsX_

private

Definition at line 288 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_convVtxdYVsY_

MonitorElement* TkConvValidator::p_convVtxdYVsY_

private

Definition at line 289 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_convVtxdZVsR_

MonitorElement* TkConvValidator::p_convVtxdZVsR_

private

Definition at line 291 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_convVtxdZVsZ_

MonitorElement* TkConvValidator::p_convVtxdZVsZ_

private

Definition at line 290 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_DCotTracksVsEta_

MonitorElement* TkConvValidator::p_DCotTracksVsEta_

private

Definition at line 246 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_DCotTracksVsR_

MonitorElement* TkConvValidator::p_DCotTracksVsR_

private

Definition at line 248 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_DPhiTracksAtEcalVsEta_

MonitorElement* TkConvValidator::p_DPhiTracksAtEcalVsEta_

private

Definition at line 256 of file TkConvValidator.h.

◆ p_DPhiTracksAtEcalVsR_

MonitorElement* TkConvValidator::p_DPhiTracksAtEcalVsR_

private

Definition at line 254 of file TkConvValidator.h.

◆ p_DPhiTracksAtVtxVsEta_

MonitorElement* TkConvValidator::p_DPhiTracksAtVtxVsEta_

private

Definition at line 240 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_DPhiTracksAtVtxVsR_

MonitorElement* TkConvValidator::p_DPhiTracksAtVtxVsR_

private

Definition at line 242 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_dzPVVsR_

MonitorElement* TkConvValidator::p_dzPVVsR_

private

Definition at line 303 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_nHitsVsEta_

MonitorElement* TkConvValidator::p_nHitsVsEta_[3]

private

Definition at line 318 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_nHitsVsR_

MonitorElement* TkConvValidator::p_nHitsVsR_[3]

private

Definition at line 320 of file TkConvValidator.h.

Referenced by analyze(), and bookHistograms().

◆ p_TkPtPull_

MonitorElement* TkConvValidator::p_TkPtPull_[3]

private

Definition at line 333 of file TkConvValidator.h.

Referenced by bookHistograms().

◆ parameters_

edm::ParameterSet TkConvValidator::parameters_

private

Definition at line 88 of file TkConvValidator.h.

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

◆ photonCollection_

std::string TkConvValidator::photonCollection_

private

Definition at line 99 of file TkConvValidator.h.

◆ photonCollectionPr_Token_

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

private

Definition at line 100 of file TkConvValidator.h.

Referenced by analyze().

◆ photonCollectionProducer_

std::string TkConvValidator::photonCollectionProducer_

private

Definition at line 98 of file TkConvValidator.h.

◆ recMaxPt_

double TkConvValidator::recMaxPt_

private

Definition at line 181 of file TkConvValidator.h.

◆ recMinPt_

double TkConvValidator::recMinPt_

private

Global variables for reco Photon.

Definition at line 180 of file TkConvValidator.h.

◆ simMaxPt_

double TkConvValidator::simMaxPt_

private

Definition at line 177 of file TkConvValidator.h.

◆ simMinPt_

double TkConvValidator::simMinPt_

private

Definition at line 176 of file TkConvValidator.h.

◆ theCaloGeom_

edm::ESHandle<CaloGeometry> TkConvValidator::theCaloGeom_

private

Definition at line 89 of file TkConvValidator.h.

Referenced by analyze().

◆ theCaloTopo_

edm::ESHandle<CaloTopology> TkConvValidator::theCaloTopo_

private

Definition at line 90 of file TkConvValidator.h.

◆ theConvTP_

edm::RefVector<TrackingParticleCollection> TkConvValidator::theConvTP_

private

Definition at line 164 of file TkConvValidator.h.

Referenced by analyze().

◆ theMF_

edm::ESHandle<MagneticField> TkConvValidator::theMF_

private

Definition at line 75 of file TkConvValidator.h.

Referenced by analyze(), and dqmBeginRun().

◆ thePhotonMCTruthFinder_

PhotonMCTruthFinder* TkConvValidator::thePhotonMCTruthFinder_

private

Definition at line 121 of file TkConvValidator.h.

Referenced by analyze(), dqmBeginRun(), and dqmEndRun().

◆ tpSelForEff_Token_

edm::EDGetTokenT<TrackingParticleRefVector> TkConvValidator::tpSelForEff_Token_

private

Definition at line 106 of file TkConvValidator.h.

Referenced by analyze().

◆ tpSelForFake_Token_

edm::EDGetTokenT<TrackingParticleRefVector> TkConvValidator::tpSelForFake_Token_

private

Definition at line 107 of file TkConvValidator.h.

Referenced by analyze().

◆ trackAssociator_Token_

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

private

Definition at line 110 of file TkConvValidator.h.

Referenced by analyze().

◆ trackerGeometryToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TkConvValidator::trackerGeometryToken_

private

Definition at line 115 of file TkConvValidator.h.

Referenced by analyze().

◆ transientTrackBuilderToken_

const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> TkConvValidator::transientTrackBuilderToken_

private

Definition at line 114 of file TkConvValidator.h.

Referenced by analyze().

◆ trkIsolExtRadius_

double TkConvValidator::trkIsolExtRadius_

private

Definition at line 126 of file TkConvValidator.h.

◆ trkIsolInnRadius_

double TkConvValidator::trkIsolInnRadius_

private

Definition at line 127 of file TkConvValidator.h.

◆ trkPtLow_

double TkConvValidator::trkPtLow_

private

Definition at line 128 of file TkConvValidator.h.

◆ trkPtSumCut_

double TkConvValidator::trkPtSumCut_

private

Definition at line 136 of file TkConvValidator.h.

◆ verbosity_

int TkConvValidator::verbosity_

private

Definition at line 77 of file TkConvValidator.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TkConvValidator()

◆ ~TkConvValidator()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

◆ dqmEndRun()

◆ endJob()

◆ etaTransformation()

◆ phiNormalization()

◆ recalculateMomentumAtFittedVertex()

Member Data Documentation

◆ arbitratedEcalSeeded_

◆ arbitratedMerged_

◆ bcEtLow_

◆ beamspotToken_

◆ caloGeometryToken_

◆ conversionCollection_

◆ conversionCollectionPr_Token_

◆ conversionCollectionProducer_

◆ conversionTrackProducer_

◆ dbe_

◆ dCotCutOn_

◆ dCotCutValue_

◆ dCotHardCutValue_

◆ dqmpath_

◆ ecalalgotracks_

◆ ecalEtSumCut_

◆ ecalIsolRadius_

◆ fName_

◆ g4_simTk_Token_

◆ g4_simVtx_Token_

◆ generalTracksOnly_

◆ genjets_Token_

◆ h2_Chi2VsEta_

◆ h2_Chi2VsR_

◆ h2_convVtxdRVsEta_

◆ h2_convVtxdRVsR_

◆ h2_convVtxRrecVsTrue_

◆ h2_DCotTracksVsEta_

◆ h2_DCotTracksVsR_

◆ h2_DPhiTracksAtEcalVsEta_

◆ h2_DPhiTracksAtEcalVsR_

◆ h2_DPhiTracksAtVtxVsEta_

◆ h2_DPhiTracksAtVtxVsR_

◆ h2_dzPVVsR_

◆ h2_photonPtRecVsPtSim_

◆ h2_PtRecVsPtSim_

◆ h2_TkPtPull_

◆ h_AllSimConv_

◆ h_convEta2_

◆ h_convEta_

◆ h_convEtaMatchSC_

◆ h_convPhi_

◆ h_convPt_

◆ h_convPtRes_

◆ h_convR_

◆ h_convRplot_

◆ h_convSCdEta_

◆ h_convSCdPhi_

◆ h_convVtxdEta_

◆ h_convVtxdPhi_

◆ h_convVtxdR_

◆ h_convVtxdR_barrel_

◆ h_convVtxdR_endcap_

◆ h_convVtxdX_

◆ h_convVtxdX_barrel_

◆ h_convVtxdX_endcap_

◆ h_convVtxdY_

◆ h_convVtxdY_barrel_

◆ h_convVtxdY_endcap_

◆ h_convVtxdZ_

◆ h_convVtxdZ_barrel_

◆ h_convVtxdZ_endcap_