CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
TkConvValidator Class Reference

#include <TkConvValidator.h>

Inheritance diagram for TkConvValidator:
DQMOneEDAnalyzer<> edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

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
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
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 &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

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

Private Attributes

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

Additional Inherited Members

- Public Types inherited from 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 DQMOneEDAnalyzer<>
virtual void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)=0
 
- Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType > produces ()
 declare what type of product will make and with which optional label More...
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces ()
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces ()
 
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)
 
template<Transition B>
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 (std::string instanceName)
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (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)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
- Protected Attributes inherited from DQMOneEDAnalyzer<>
edm::EDPutTokenT< DQMToken > runToken_
 

Detailed Description

$Id: TkConvValidator

Author
N.Marinelli - Univ. of Notre Dame

Definition at line 49 of file TkConvValidator.h.

Constructor & Destructor Documentation

◆ TkConvValidator()

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

Definition at line 99 of file TkConvValidator.cc.

99  {
100  fName_ = pset.getUntrackedParameter<std::string>("Name");
101  verbosity_ = pset.getUntrackedParameter<int>("Verbosity");
102  parameters_ = pset;
103 
104  photonCollectionProducer_ = pset.getParameter<std::string>("phoProducer");
105  photonCollection_ = pset.getParameter<std::string>("photonCollection");
106  photonCollectionPr_Token_ =
107  consumes<reco::PhotonCollection>(edm::InputTag(photonCollectionProducer_, photonCollection_));
108 
109  conversionCollectionProducer_ = pset.getParameter<std::string>("convProducer");
110  conversionCollection_ = pset.getParameter<std::string>("conversionCollection");
111  conversionCollectionPr_Token_ =
112  consumes<reco::ConversionCollection>(edm::InputTag(conversionCollectionProducer_, conversionCollection_));
113 
114  // conversionTrackProducer_ = pset.getParameter<std::string>("trackProducer");
115  dqmpath_ = pset.getParameter<std::string>("dqmpath");
116  minPhoEtCut_ = pset.getParameter<double>("minPhoEtCut");
117  generalTracksOnly_ = pset.getParameter<bool>("generalTracksOnly");
118  arbitratedMerged_ = pset.getParameter<bool>("arbitratedMerged");
119  arbitratedEcalSeeded_ = pset.getParameter<bool>("arbitratedEcalSeeded");
120  ecalalgotracks_ = pset.getParameter<bool>("ecalalgotracks");
121  highPurity_ = pset.getParameter<bool>("highPurity");
122  minProb_ = pset.getParameter<double>("minProb");
123  maxHitsBeforeVtx_ = pset.getParameter<uint>("maxHitsBeforeVtx");
124  minLxy_ = pset.getParameter<double>("minLxy");
125  isRunCentrally_ = pset.getParameter<bool>("isRunCentrally");
126 
127  offline_pvToken_ = consumes<reco::VertexCollection>(
128  pset.getUntrackedParameter<edm::InputTag>("offlinePV", edm::InputTag("offlinePrimaryVertices")));
129  beamspotToken_ =
130  consumes<reco::BeamSpot>(pset.getUntrackedParameter<edm::InputTag>("beamspot", edm::InputTag("offlineBeamSpot")));
131  g4_simTk_Token_ = consumes<edm::SimTrackContainer>(pset.getParameter<edm::InputTag>("simTracks"));
132  g4_simVtx_Token_ = consumes<edm::SimVertexContainer>(pset.getParameter<edm::InputTag>("simTracks"));
133 
134  tpSelForEff_Token_ = consumes<TrackingParticleRefVector>(edm::InputTag("tpSelecForEfficiency"));
135  tpSelForFake_Token_ = consumes<TrackingParticleRefVector>(edm::InputTag("tpSelecForFakeRate"));
136  //hepMC_Token_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));
137  //genjets_Token_ = consumes<reco::GenJetCollection>(
138  // edm::InputTag("ak4GenJets"));
139 
140  trackAssociator_Token_ =
141  consumes<reco::TrackToTrackingParticleAssociator>(edm::InputTag("trackAssociatorByHitsForConversionValidation"));
142 }

References muonDTDigis_cfi::pset, AlCaHLTBitMon_QueryRunRegistry::string, and parallelization::uint.

◆ ~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 1239 of file TkConvValidator.cc.

1239  {
1240  thePhotonMCTruthFinder_->clear();
1241  using namespace edm;
1242  // const float etaPhiDistance=0.01;
1243  // Fiducial region
1244  // const float TRK_BARL =0.9;
1245  const float BARL = 1.4442; // DAQ TDR p.290
1246  // const float END_LO = 1.566; // unused
1247  const float END_HI = 2.5;
1248  // Electron mass
1249  // const Float_t mElec= 0.000511; // unused
1250 
1251  Handle<reco::TrackToTrackingParticleAssociator> theTrackAssociator;
1252  e.getByToken(trackAssociator_Token_, theTrackAssociator);
1253 
1254  nEvt_++;
1255  LogInfo("TkConvValidator") << "TkConvValidator Analyzing event number: " << e.id() << " Global Counter " << nEvt_
1256  << "\n";
1257  // std::cout << "TkConvValidator Analyzing event number: " << e.id() << " Global Counter " << nEvt_ <<"\n";
1258 
1259  // get the geometry from the event setup:
1260  esup.get<CaloGeometryRecord>().get(theCaloGeom_);
1261 
1262  // Transform Track into TransientTrack (needed by the Vertex fitter)
1263  edm::ESHandle<TransientTrackBuilder> theTTB;
1264  esup.get<TransientTrackRecord>().get("TransientTrackBuilder", theTTB);
1265 
1267  Handle<reco::ConversionCollection> convHandle;
1268  e.getByToken(conversionCollectionPr_Token_, convHandle);
1269  const reco::ConversionCollection convCollection = *(convHandle.product());
1270  if (!convHandle.isValid()) {
1271  edm::LogError("ConversionsProducer") << "Error! Can't get the collection " << std::endl;
1272  return;
1273  }
1274 
1276  Handle<reco::PhotonCollection> photonHandle;
1277  e.getByToken(photonCollectionPr_Token_, photonHandle);
1278  const reco::PhotonCollection photonCollection = *(photonHandle.product());
1279  if (!photonHandle.isValid()) {
1280  edm::LogError("PhotonProducer") << "Error! Can't get the Photon collection " << std::endl;
1281  return;
1282  }
1283 
1284  // offline Primary vertex
1285  edm::Handle<reco::VertexCollection> vertexHandle;
1286  reco::VertexCollection vertexCollection;
1287  e.getByToken(offline_pvToken_, vertexHandle);
1288  if (!vertexHandle.isValid()) {
1289  edm::LogError("TrackerOnlyConversionProducer") << "Error! Can't get the product primary Vertex Collection "
1290  << "\n";
1291  } else {
1292  vertexCollection = *(vertexHandle.product());
1293  }
1294  reco::Vertex the_pvtx;
1295  bool valid_pvtx = false;
1296  if (!vertexCollection.empty()) {
1297  the_pvtx = *(vertexCollection.begin());
1298  //asking for one good vertex
1299  if (the_pvtx.isValid() && fabs(the_pvtx.position().z()) <= 15 && the_pvtx.position().Rho() <= 2) {
1300  valid_pvtx = true;
1301  }
1302  }
1303 
1304  edm::Handle<reco::BeamSpot> bsHandle;
1305  e.getByToken(beamspotToken_, bsHandle);
1306  if (!bsHandle.isValid()) {
1307  edm::LogError("TrackerOnlyConversionProducer") << "Error! Can't get the product primary Vertex Collection "
1308  << "\n";
1309  return;
1310  }
1311  const reco::BeamSpot& thebs = *bsHandle.product();
1312 
1313  //get tracker geometry for hits positions
1314  edm::ESHandle<TrackerGeometry> tracker;
1315  esup.get<TrackerDigiGeometryRecord>().get(tracker);
1316  const TrackerGeometry* trackerGeom = tracker.product();
1317 
1319  //get simtrack info
1320  //std::vector<SimTrack> theSimTracks;
1321  //std::vector<SimVertex> theSimVertices;
1322 
1323  edm::Handle<SimTrackContainer> SimTk;
1324  edm::Handle<SimVertexContainer> SimVtx;
1325  e.getByToken(g4_simTk_Token_, SimTk);
1326  e.getByToken(g4_simVtx_Token_, SimVtx);
1327 
1328  bool useTP = parameters_.getParameter<bool>("useTP");
1329  TrackingParticleRefVector dummy;
1330  edm::Handle<TrackingParticleRefVector> TPHandleForEff;
1331  edm::Handle<TrackingParticleRefVector> TPHandleForFakeRate;
1332  if (useTP) {
1333  e.getByToken(tpSelForEff_Token_, TPHandleForEff);
1334  e.getByToken(tpSelForFake_Token_, TPHandleForFakeRate);
1335  }
1336 
1337  const TrackingParticleRefVector& tpForEfficiency = useTP ? *(TPHandleForEff.product()) : dummy;
1338  const TrackingParticleRefVector& tpForFakeRate = useTP ? *(TPHandleForFakeRate.product()) : dummy;
1339 
1340  const std::vector<SimTrack>& theSimTracks = *SimTk;
1341  const std::vector<SimVertex>& theSimVertices = *SimVtx;
1342 
1343  //theSimTracks.insert(theSimTracks.end(),SimTk->begin(),SimTk->end());
1344  //theSimVertices.insert(theSimVertices.end(),SimVtx->begin(),SimVtx->end());
1345  std::vector<PhotonMCTruth> mcPhotons = thePhotonMCTruthFinder_->find(theSimTracks, theSimVertices);
1346 
1347  //edm::Handle<edm::HepMCProduct> hepMC;
1348  //e.getByToken(hepMC_Token_, hepMC);
1349  // const HepMC::GenEvent *myGenEvent = hepMC->GetEvent(); // unused
1350 
1351  // get generated jets
1352  //edm::Handle<reco::GenJetCollection> GenJetsHandle;
1353  //e.getByToken(genjets_Token_, GenJetsHandle);
1354  //const reco::GenJetCollection &genJetCollection = *(GenJetsHandle.product());
1355 
1356  // ################ SIM to RECO ######################### //
1357  std::map<const reco::Track*, TrackingParticleRef> myAss;
1358  std::map<const reco::Track*, TrackingParticleRef>::const_iterator itAss;
1359 
1360  for (std::vector<PhotonMCTruth>::const_iterator mcPho = mcPhotons.begin(); mcPho != mcPhotons.end(); mcPho++) {
1361  mcConvPt_ = (*mcPho).fourMomentum().et();
1362  float mcPhi = (*mcPho).fourMomentum().phi();
1363  mcPhi_ = phiNormalization(mcPhi);
1364  mcEta_ = (*mcPho).fourMomentum().pseudoRapidity();
1365  mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z());
1366  mcConvR_ = (*mcPho).vertex().perp();
1367  mcConvX_ = (*mcPho).vertex().x();
1368  mcConvY_ = (*mcPho).vertex().y();
1369  mcConvZ_ = (*mcPho).vertex().z();
1370  mcConvEta_ = (*mcPho).vertex().eta();
1371  mcConvPhi_ = (*mcPho).vertex().phi();
1372 
1373  if (fabs(mcEta_) > END_HI)
1374  continue;
1375 
1376  if (mcConvPt_ < minPhoPtForEffic)
1377  continue;
1378  if (fabs(mcEta_) > maxPhoEtaForEffic)
1379  continue;
1380  if (fabs(mcConvZ_) > maxPhoZForEffic)
1381  continue;
1382  if (mcConvR_ > maxPhoRForEffic)
1383  continue;
1385 
1386  bool goodSimConversion = false;
1387  bool visibleConversion = false;
1388  bool visibleConversionsWithTwoSimTracks = false;
1389  if ((*mcPho).isAConversion() == 1) {
1390  nSimConv_[0]++;
1391  h_AllSimConv_[0]->Fill(mcEta_);
1392  h_AllSimConv_[1]->Fill(mcPhi_);
1393  h_AllSimConv_[2]->Fill(mcConvR_);
1394  h_AllSimConv_[3]->Fill(mcConvZ_);
1395  h_AllSimConv_[4]->Fill((*mcPho).fourMomentum().et());
1396 
1397  if (mcConvR_ < 15)
1398  h_SimConvEtaPix_[0]->Fill(mcEta_);
1399 
1400  if ((fabs(mcEta_) <= BARL && mcConvR_ < 85) ||
1401  (fabs(mcEta_) > BARL && fabs(mcEta_) <= END_HI && fabs((*mcPho).vertex().z()) < 210))
1402  visibleConversion = true;
1403 
1404  theConvTP_.clear();
1405  // std::cout << " TkConvValidator TrackingParticles TrackingParticleCollection size "<< trackingParticles.size() << "\n";
1406  //duplicated TP collections for two associations
1407  for (const TrackingParticleRef tp : tpForEfficiency) {
1408  if (fabs(tp->vx() - (*mcPho).vertex().x()) < 0.0001 && fabs(tp->vy() - (*mcPho).vertex().y()) < 0.0001 &&
1409  fabs(tp->vz() - (*mcPho).vertex().z()) < 0.0001) {
1410  theConvTP_.push_back(tp);
1411  }
1412  }
1413  //std::cout << " TkConvValidator theConvTP_ size " << theConvTP_.size() << std::endl;
1414 
1415  if (theConvTP_.size() == 2)
1416  visibleConversionsWithTwoSimTracks = true;
1417  goodSimConversion = false;
1418 
1419  if (visibleConversion && visibleConversionsWithTwoSimTracks)
1420  goodSimConversion = true;
1421  if (goodSimConversion) {
1422  nSimConv_[1]++;
1423  h_VisSimConv_[0]->Fill(mcEta_);
1424  h_VisSimConv_[1]->Fill(mcPhi_);
1425  h_VisSimConv_[2]->Fill(mcConvR_);
1426  h_VisSimConv_[3]->Fill(mcConvZ_);
1427  h_VisSimConv_[4]->Fill((*mcPho).fourMomentum().et());
1428  }
1429 
1430  for (edm::RefVector<TrackingParticleCollection>::iterator iTrk = theConvTP_.begin(); iTrk != theConvTP_.end();
1431  ++iTrk) {
1432  h_simTkPt_->Fill((*iTrk)->pt());
1433  h_simTkEta_->Fill((*iTrk)->eta());
1434  }
1435 
1436  }
1437 
1438  if (!(visibleConversion && visibleConversionsWithTwoSimTracks))
1439  continue;
1440 
1441  h_simConvVtxRvsZ_[0]->Fill(fabs(mcConvZ_), mcConvR_);
1442  if (fabs(mcEta_) <= 1.) {
1443  h_simConvVtxRvsZ_[1]->Fill(fabs(mcConvZ_), mcConvR_);
1444  h_simConvVtxYvsX_->Fill(mcConvX_, mcConvY_);
1445  } else
1446  h_simConvVtxRvsZ_[2]->Fill(fabs(mcConvZ_), mcConvR_);
1447 
1448  //std::cout << " TkConvValidator theConvTP_ size " << theConvTP_.size() << std::endl;
1449  for (edm::RefVector<TrackingParticleCollection>::iterator iTP = theConvTP_.begin(); iTP != theConvTP_.end();
1450  iTP++) {
1451  // std::cout << " SIM to RECO TP vertex " << (*iTP)->vx() << " " << (*iTP)->vy() << " " << (*iTP)->vz() << " pt " << (*iTP)->pt() << std::endl;
1452  }
1453 
1454  bool recomatch = false;
1455  float chi2Prob = 0.;
1457  // cout << " size of conversions " << convHandle->size() << endl;
1458  for (reco::ConversionCollection::const_iterator conv = convHandle->begin(); conv != convHandle->end(); ++conv) {
1459  const reco::Conversion aConv = (*conv);
1460  if (arbitratedMerged_ && !aConv.quality(reco::Conversion::arbitratedMerged))
1461  continue;
1462  if (generalTracksOnly_ && !aConv.quality(reco::Conversion::generalTracksOnly))
1463  continue;
1464  if (arbitratedEcalSeeded_ && !aConv.quality(reco::Conversion::arbitratedEcalSeeded))
1465  continue;
1466 
1467  if (highPurity_ && !aConv.quality(reco::Conversion::highPurity))
1468  continue;
1469 
1470  //problematic?
1471  const std::vector<edm::RefToBase<reco::Track> >& tracks = aConv.tracks();
1472 
1473  const reco::Vertex& vtx = aConv.conversionVertex();
1474  //requires two tracks and a valid vertex
1475  if (tracks.size() != 2 || !(vtx.isValid()))
1476  continue;
1477 
1478  if (ChiSquaredProbability(aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof()) <= minProb_)
1479  continue;
1480  if (aConv.nHitsBeforeVtx().size() > 1 &&
1481  max(aConv.nHitsBeforeVtx().at(0), aConv.nHitsBeforeVtx().at(1)) > maxHitsBeforeVtx_)
1482  continue;
1483 
1484  //compute transverse decay length with respect to beamspot
1485  math::XYZVectorF themom = aConv.refittedPairMomentum();
1486  double dbsx = aConv.conversionVertex().x() - thebs.x0();
1487  double dbsy = aConv.conversionVertex().y() - thebs.y0();
1488  double lxy = (themom.x() * dbsx + themom.y() * dbsy) / themom.rho();
1489 
1490  if (lxy < minLxy_)
1491  continue;
1492 
1493  // bool phoIsInBarrel=false; // unused
1494  // bool phoIsInEndcap=false; // unused
1495  RefToBase<reco::Track> tfrb1 = aConv.tracks().front();
1496  RefToBase<reco::Track> tfrb2 = aConv.tracks().back();
1497 
1498  if (ecalalgotracks_ && (!(tfrb1->algo() == reco::TrackBase::outInEcalSeededConv ||
1499  tfrb1->algo() == reco::TrackBase::inOutEcalSeededConv) ||
1500  !(tfrb2->algo() == reco::TrackBase::outInEcalSeededConv ||
1501  tfrb2->algo() == reco::TrackBase::inOutEcalSeededConv)))
1502  continue;
1503 
1504  //reco::TrackRef tk1 = aConv.tracks().front();
1505  //reco::TrackRef tk2 = aConv.tracks().back();
1506  //std::cout << "SIM to RECO conversion track pt " << tk1->pt() << " " << tk2->pt() << endl;
1507  //
1508  //Use two RefToBaseVector and do two association actions to avoid that if two tracks from different collection
1509  RefToBaseVector<reco::Track> tc1, tc2;
1510  tc1.push_back(tfrb1);
1511  tc2.push_back(tfrb2);
1512  bool isAssociated = false;
1513  reco::SimToRecoCollection q1 = theTrackAssociator->associateSimToReco(tc1, theConvTP_);
1514  reco::SimToRecoCollection q2 = theTrackAssociator->associateSimToReco(tc2, theConvTP_);
1515  //try {
1516  std::vector<std::pair<RefToBase<reco::Track>, double> > trackV1, trackV2;
1517 
1518  int tp_1 = 0, tp_2 = 1; //the index of associated tp in theConvTP_ for two tracks
1519  if (q1.find(theConvTP_[0]) != q1.end()) {
1520  trackV1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >)q1[theConvTP_[0]];
1521  } else if (q1.find(theConvTP_[1]) != q1.end()) {
1522  trackV1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >)q1[theConvTP_[1]];
1523  tp_1 = 1;
1524  }
1525  if (q2.find(theConvTP_[1]) != q2.end()) {
1526  trackV2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >)q2[theConvTP_[1]];
1527  } else if (q2.find(theConvTP_[0]) != q2.end()) {
1528  trackV2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >)q2[theConvTP_[0]];
1529  tp_2 = 0;
1530  }
1531  if (!(!trackV1.empty() && !trackV2.empty()))
1532  continue;
1533  if (tp_1 == tp_2)
1534  continue;
1535 
1536  edm::RefToBase<reco::Track> tr1 = trackV1.front().first;
1537  edm::RefToBase<reco::Track> tr2 = trackV2.front().first;
1538  //std::cout << "associated tp1 " <<theConvTP_[0]->pt() << " to track with pT=" << tr1->pt() << " " << (tr1.get())->pt() << endl;
1539  //std::cout << "associated tp2 " <<theConvTP_[1]->pt() << " to track with pT=" << tr2->pt() << " " << (tr2.get())->pt() << endl;
1540  myAss.insert(std::make_pair(tr1.get(), theConvTP_[tp_1]));
1541  myAss.insert(std::make_pair(tr2.get(), theConvTP_[tp_2]));
1542 
1543  //} catch (Exception const& event) {
1544  //cout << "continue: " << event.what() << endl;
1545  // continue;
1546  //}
1547 
1548  isAssociated = true;
1549  recomatch = true;
1550  chi2Prob = ChiSquaredProbability(aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof());
1551 
1552  if (isAssociated) {
1553  h_SimRecConvTwoMTracks_[0]->Fill(mcEta_);
1554  h_SimRecConvTwoMTracks_[1]->Fill(mcPhi_);
1555  h_SimRecConvTwoMTracks_[2]->Fill(mcConvR_);
1556  h_SimRecConvTwoMTracks_[3]->Fill(mcConvZ_);
1557  h_SimRecConvTwoMTracks_[4]->Fill((*mcPho).fourMomentum().et());
1558  }
1559 
1560  // break;
1561  } // loop over reco conversions
1562  if (recomatch) {
1564  h_SimConvTwoMTracks_[0]->Fill(mcEta_);
1565  h_SimConvTwoMTracks_[1]->Fill(mcPhi_);
1566  h_SimConvTwoMTracks_[2]->Fill(mcConvR_);
1567  h_SimConvTwoMTracks_[3]->Fill(mcConvZ_);
1568  h_SimConvTwoMTracks_[4]->Fill((*mcPho).fourMomentum().et());
1569 
1570  if (chi2Prob > 0) {
1571  h_SimConvTwoMTracksAndVtxPGT0_[0]->Fill(mcEta_);
1572  h_SimConvTwoMTracksAndVtxPGT0_[1]->Fill(mcPhi_);
1573  h_SimConvTwoMTracksAndVtxPGT0_[2]->Fill(mcConvR_);
1574  h_SimConvTwoMTracksAndVtxPGT0_[3]->Fill(mcConvZ_);
1575  h_SimConvTwoMTracksAndVtxPGT0_[4]->Fill((*mcPho).fourMomentum().et());
1576  }
1577  if (chi2Prob > 0.0005) {
1578  h_SimConvTwoMTracksAndVtxPGT0005_[0]->Fill(mcEta_);
1579  h_SimConvTwoMTracksAndVtxPGT0005_[1]->Fill(mcPhi_);
1580  h_SimConvTwoMTracksAndVtxPGT0005_[2]->Fill(mcConvR_);
1581  h_SimConvTwoMTracksAndVtxPGT0005_[3]->Fill(mcConvZ_);
1582  h_SimConvTwoMTracksAndVtxPGT0005_[4]->Fill((*mcPho).fourMomentum().et());
1583  }
1584  }
1585 
1586  } //End loop over simulated conversions
1587 
1588  // ########################### RECO to SIM ############################## //
1589 
1590  for (reco::ConversionCollection::const_iterator conv = convHandle->begin(); conv != convHandle->end(); ++conv) {
1591  const reco::Conversion aConv = (*conv);
1592  if (arbitratedMerged_ && !aConv.quality(reco::Conversion::arbitratedMerged))
1593  continue;
1594  if (generalTracksOnly_ && !aConv.quality(reco::Conversion::generalTracksOnly))
1595  continue;
1596  if (arbitratedEcalSeeded_ && !aConv.quality(reco::Conversion::arbitratedEcalSeeded))
1597  continue;
1598 
1599  if (highPurity_ && !aConv.quality(reco::Conversion::highPurity))
1600  continue;
1601 
1602  //problematic?
1603  std::vector<edm::RefToBase<reco::Track> > tracks = aConv.tracks();
1604 
1605  const reco::Vertex& vtx = aConv.conversionVertex();
1606  //requires two tracks and a valid vertex
1607  if (tracks.size() != 2 || !(vtx.isValid()))
1608  continue;
1609  //if (tracks.size() !=2) continue;
1610 
1611  if (ChiSquaredProbability(aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof()) <= minProb_)
1612  continue;
1613  if (aConv.nHitsBeforeVtx().size() > 1 &&
1614  max(aConv.nHitsBeforeVtx().at(0), aConv.nHitsBeforeVtx().at(1)) > maxHitsBeforeVtx_)
1615  continue;
1616 
1617  //compute transverse decay length with respect to beamspot
1618  math::XYZVectorF themom = aConv.refittedPairMomentum();
1619  double dbsx = aConv.conversionVertex().x() - thebs.x0();
1620  double dbsy = aConv.conversionVertex().y() - thebs.y0();
1621  double lxy = (themom.x() * dbsx + themom.y() * dbsy) / themom.rho();
1622 
1623  if (lxy < minLxy_)
1624  continue;
1625 
1626  bool phoIsInBarrel = false;
1627  bool phoIsInEndcap = false;
1628  RefToBase<reco::Track> tk1 = aConv.tracks().front();
1629  RefToBase<reco::Track> tk2 = aConv.tracks().back();
1630  RefToBaseVector<reco::Track> tc1, tc2;
1631  tc1.push_back(tk1);
1632  tc2.push_back(tk2);
1633 
1634  if (ecalalgotracks_ &&
1635  (!(tk1->algo() == reco::TrackBase::outInEcalSeededConv ||
1636  tk1->algo() == reco::TrackBase::inOutEcalSeededConv) ||
1637  !(tk2->algo() == reco::TrackBase::outInEcalSeededConv || tk2->algo() == reco::TrackBase::inOutEcalSeededConv)))
1638  continue;
1639 
1640  //std::cout << " RECO to SIM conversion track pt " << tk1->pt() << " " << tk2->pt() << endl;
1641  const reco::Track refTk1 = aConv.conversionVertex().refittedTracks().front();
1642  const reco::Track refTk2 = aConv.conversionVertex().refittedTracks().back();
1643 
1644  float dPhiTracksAtVtx;
1645  // override with the phi calculated at the vertex
1646  math::XYZVector p1AtVtx = recalculateMomentumAtFittedVertex((*theMF_), *trackerGeom, tk1, aConv.conversionVertex());
1647  math::XYZVector p2AtVtx = recalculateMomentumAtFittedVertex((*theMF_), *trackerGeom, tk2, aConv.conversionVertex());
1648  if (p1AtVtx.perp2() > p2AtVtx.perp2())
1649  dPhiTracksAtVtx = p1AtVtx.phi() - p2AtVtx.phi();
1650  else
1651  dPhiTracksAtVtx = p2AtVtx.phi() - p1AtVtx.phi();
1652 
1653  math::XYZVectorF refittedMom = aConv.refittedPairMomentum();
1654 
1655  if (fabs(refittedMom.eta()) < 1.479) {
1656  phoIsInBarrel = true;
1657  } else {
1658  phoIsInEndcap = true;
1659  }
1660 
1661  nRecConv_++;
1662 
1663  // check matching with reco photon
1664  double Mindeltaeta = 999999;
1665  double Mindeltaphi = 999999;
1666  bool matchConvSC = false;
1667  reco::PhotonCollection::const_iterator iMatchingSC;
1668  for (reco::PhotonCollection::const_iterator iPho = photonCollection.begin(); iPho != photonCollection.end();
1669  iPho++) {
1670  reco::Photon aPho = reco::Photon(*iPho);
1671  const double deltaphi = reco::deltaPhi(aConv.refittedPairMomentum().phi(), aPho.superCluster()->position().phi());
1672  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
1673  double deltaeta = abs(aPho.superCluster()->position().eta() - ConvEta);
1674  if (abs(deltaeta) < abs(Mindeltaeta) && abs(deltaphi) < abs(Mindeltaphi)) {
1675  Mindeltaphi = abs(deltaphi);
1676  Mindeltaeta = abs(deltaeta);
1677  iMatchingSC = iPho;
1678  }
1679  }
1680  if (abs(Mindeltaeta) < 0.1 && abs(Mindeltaphi) < 0.1) {
1681  matchConvSC = true;
1682  }
1683 
1685  int match = 0;
1686  float invM = aConv.pairInvariantMass();
1687  float chi2Prob = ChiSquaredProbability(aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof());
1688  uint maxNHitsBeforeVtx =
1689  aConv.nHitsBeforeVtx().size() > 1 ? max(aConv.nHitsBeforeVtx().at(0), aConv.nHitsBeforeVtx().at(1)) : 0;
1690  uint sumNHitsBeforeVtx =
1691  aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(0) + aConv.nHitsBeforeVtx().at(1) : 0;
1692  float maxDlClosestHitToVtx = aConv.dlClosestHitToVtx().size() > 1 ? max(aConv.dlClosestHitToVtx().at(0).value(),
1693  aConv.dlClosestHitToVtx().at(1).value())
1694  : 0;
1695  float maxDlClosestHitToVtxSig =
1696  aConv.dlClosestHitToVtx().size() > 1
1697  ? max(aConv.dlClosestHitToVtx().at(0).value() / aConv.dlClosestHitToVtx().at(0).error(),
1698  aConv.dlClosestHitToVtx().at(1).value() / aConv.dlClosestHitToVtx().at(1).error())
1699  : 0;
1700 
1701  int ilead = 0, itrail = 1;
1702  if (tk2->pt() > tk1->pt()) {
1703  ilead = 1;
1704  itrail = 0;
1705  }
1706  RefToBase<reco::Track> tklead = aConv.tracks().at(ilead);
1707  RefToBase<reco::Track> tktrail = aConv.tracks().at(itrail);
1708 
1709  int deltaExpectedHitsInner = tklead->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS) -
1710  tktrail->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
1711  int leadExpectedHitsInner = tklead->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
1712  uint leadNHitsBeforeVtx = aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(ilead) : 0;
1713  uint trailNHitsBeforeVtx = aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(itrail) : 0;
1714 
1715  h_convEta_[match][0]->Fill(refittedMom.eta());
1716  h_convEta2_[match][0]->Fill(refittedMom.eta());
1717 
1718  h_convPhi_[match][0]->Fill(refittedMom.phi());
1719  h_convR_[match][0]->Fill(sqrt(aConv.conversionVertex().position().perp2()));
1720  h_convRplot_->Fill(sqrt(aConv.conversionVertex().position().perp2()));
1721  h_convZ_[match][0]->Fill(aConv.conversionVertex().position().z());
1722  h_convZplot_->Fill(aConv.conversionVertex().position().z());
1723  h_convPt_[match][0]->Fill(sqrt(refittedMom.perp2()));
1724  h_invMass_[match][0]->Fill(invM);
1725  h_vtxChi2Prob_[match][0]->Fill(chi2Prob);
1726  h_lxybs_[match][0]->Fill(lxy);
1727  h_maxNHitsBeforeVtx_[match][0]->Fill(maxNHitsBeforeVtx);
1728  h_leadNHitsBeforeVtx_[match][0]->Fill(leadNHitsBeforeVtx);
1729  h_trailNHitsBeforeVtx_[match][0]->Fill(trailNHitsBeforeVtx);
1730  h_sumNHitsBeforeVtx_[match][0]->Fill(sumNHitsBeforeVtx);
1731  h_deltaExpectedHitsInner_[match][0]->Fill(deltaExpectedHitsInner);
1732  h_leadExpectedHitsInner_[match][0]->Fill(leadExpectedHitsInner);
1733  h_maxDlClosestHitToVtx_[match][0]->Fill(maxDlClosestHitToVtx);
1734  h_maxDlClosestHitToVtxSig_[match][0]->Fill(maxDlClosestHitToVtxSig);
1735  h_nSharedHits_[match][0]->Fill(aConv.nSharedHits());
1736 
1737  if (matchConvSC) {
1738  h_convEtaMatchSC_[match][0]->Fill(refittedMom.eta());
1739  h_EoverPTracks_[match][0]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
1740  h_convSCdPhi_[match][0]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
1741  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
1742  h_convSCdEta_[match][0]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
1743  }
1744 
1745  h_distMinAppTracks_[match][0]->Fill(aConv.distOfMinimumApproach());
1746  h_DPhiTracksAtVtx_[match][0]->Fill(dPhiTracksAtVtx);
1747  h2_DPhiTracksAtVtxVsEta_->Fill(mcEta_, dPhiTracksAtVtx);
1748  h2_DPhiTracksAtVtxVsR_->Fill(mcConvR_, dPhiTracksAtVtx);
1749  p_DPhiTracksAtVtxVsEta_->Fill(mcEta_, dPhiTracksAtVtx);
1750  p_DPhiTracksAtVtxVsR_->Fill(mcConvR_, dPhiTracksAtVtx);
1751 
1752  h_DCotTracks_[match][0]->Fill(aConv.pairCotThetaSeparation());
1753  h2_DCotTracksVsEta_->Fill(mcEta_, aConv.pairCotThetaSeparation());
1754  h2_DCotTracksVsR_->Fill(mcConvR_, aConv.pairCotThetaSeparation());
1755  p_DCotTracksVsEta_->Fill(mcEta_, aConv.pairCotThetaSeparation());
1756  p_DCotTracksVsR_->Fill(mcConvR_, aConv.pairCotThetaSeparation());
1757 
1758  if (phoIsInBarrel) {
1759  h_invMass_[match][1]->Fill(invM);
1760  h_vtxChi2Prob_[match][1]->Fill(chi2Prob);
1761  h_distMinAppTracks_[match][1]->Fill(aConv.distOfMinimumApproach());
1762  h_DPhiTracksAtVtx_[match][1]->Fill(dPhiTracksAtVtx);
1763  h_DCotTracks_[match][1]->Fill(aConv.pairCotThetaSeparation());
1764  h_lxybs_[match][1]->Fill(lxy);
1765  h_maxNHitsBeforeVtx_[match][1]->Fill(maxNHitsBeforeVtx);
1766  h_leadNHitsBeforeVtx_[match][1]->Fill(leadNHitsBeforeVtx);
1767  h_trailNHitsBeforeVtx_[match][1]->Fill(trailNHitsBeforeVtx);
1768  h_sumNHitsBeforeVtx_[match][1]->Fill(sumNHitsBeforeVtx);
1769  h_deltaExpectedHitsInner_[match][1]->Fill(deltaExpectedHitsInner);
1770  h_leadExpectedHitsInner_[match][1]->Fill(leadExpectedHitsInner);
1771  h_maxDlClosestHitToVtx_[match][1]->Fill(maxDlClosestHitToVtx);
1772  h_maxDlClosestHitToVtxSig_[match][1]->Fill(maxDlClosestHitToVtxSig);
1773  h_nSharedHits_[match][1]->Fill(aConv.nSharedHits());
1774 
1775  /*
1776  if ( aConv.caloCluster().size() ) {
1777  h_convSCdPhi_[match][1]->Fill( aConv.caloCluster()[0]->phi() - refittedMom.phi() );
1778  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1779  h_convSCdEta_[match][1]->Fill( aConv.caloCluster()[0]->eta() - ConvEta );
1780  }
1781  */
1782 
1783  if (matchConvSC) {
1784  h_EoverPTracks_[match][1]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
1785  h_convSCdPhi_[match][1]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
1786  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
1787  h_convSCdEta_[match][1]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
1788  }
1789  }
1790 
1791  if (phoIsInEndcap) {
1792  h_invMass_[match][2]->Fill(invM);
1793  h_vtxChi2Prob_[match][2]->Fill(chi2Prob);
1794  h_distMinAppTracks_[match][2]->Fill(aConv.distOfMinimumApproach());
1795  h_DPhiTracksAtVtx_[match][2]->Fill(dPhiTracksAtVtx);
1796  h_DCotTracks_[match][2]->Fill(aConv.pairCotThetaSeparation());
1797  h_lxybs_[match][2]->Fill(lxy);
1798  h_maxNHitsBeforeVtx_[match][2]->Fill(maxNHitsBeforeVtx);
1799  h_leadNHitsBeforeVtx_[match][2]->Fill(leadNHitsBeforeVtx);
1800  h_trailNHitsBeforeVtx_[match][2]->Fill(trailNHitsBeforeVtx);
1801  h_sumNHitsBeforeVtx_[match][2]->Fill(sumNHitsBeforeVtx);
1802  h_deltaExpectedHitsInner_[match][2]->Fill(deltaExpectedHitsInner);
1803  h_leadExpectedHitsInner_[match][2]->Fill(leadExpectedHitsInner);
1804  h_maxDlClosestHitToVtx_[match][2]->Fill(maxDlClosestHitToVtx);
1805  h_maxDlClosestHitToVtxSig_[match][2]->Fill(maxDlClosestHitToVtxSig);
1806  h_nSharedHits_[match][2]->Fill(aConv.nSharedHits());
1807  if (matchConvSC) {
1808  h_EoverPTracks_[match][2]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
1809  h_convSCdPhi_[match][2]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
1810  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
1811  h_convSCdEta_[match][2]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
1812  }
1813  }
1814 
1815  h_convVtxRvsZ_[0]->Fill(fabs(aConv.conversionVertex().position().z()),
1816  sqrt(aConv.conversionVertex().position().perp2()));
1817  h_convVtxYvsX_->Fill(aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y());
1818  h_convVtxYvsX_zoom_[0]->Fill(aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y());
1819  h_convVtxYvsX_zoom_[1]->Fill(aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y());
1820 
1821  // quantities per track: all conversions
1822  for (unsigned int i = 0; i < tracks.size(); i++) {
1823  double d0;
1824  if (valid_pvtx) {
1825  d0 = -tracks[i]->dxy(the_pvtx.position());
1826  } else {
1827  d0 = tracks[i]->d0();
1828  }
1829  h_TkD0_[match]->Fill(d0 * tracks[i]->charge());
1830  h_nHitsBeforeVtx_[match]->Fill(aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(i) : 0);
1831  h_dlClosestHitToVtx_[match]->Fill(aConv.dlClosestHitToVtx().size() > 1 ? aConv.dlClosestHitToVtx().at(i).value()
1832  : 0);
1833  h_dlClosestHitToVtxSig_[match]->Fill(aConv.dlClosestHitToVtx().size() > 1
1834  ? aConv.dlClosestHitToVtx().at(i).value() /
1835  aConv.dlClosestHitToVtx().at(i).error()
1836  : 0);
1837 
1838  nHitsVsEta_[match]->Fill(mcEta_, float(tracks[i]->numberOfValidHits()));
1839  nHitsVsR_[match]->Fill(mcConvR_, float(tracks[i]->numberOfValidHits()));
1840  p_nHitsVsEta_[match]->Fill(mcEta_, float(tracks[i]->numberOfValidHits()) - 0.0001);
1841  p_nHitsVsR_[match]->Fill(mcConvR_, float(tracks[i]->numberOfValidHits()) - 0.0001);
1842  h_tkChi2_[match]->Fill(tracks[i]->normalizedChi2());
1843  h_tkChi2Large_[match]->Fill(tracks[i]->normalizedChi2());
1844  h2_Chi2VsEta_[match]->Fill(mcEta_, tracks[i]->normalizedChi2());
1845  h2_Chi2VsR_[match]->Fill(mcConvR_, tracks[i]->normalizedChi2());
1846  p_Chi2VsEta_[match]->Fill(mcEta_, tracks[i]->normalizedChi2());
1847  p_Chi2VsR_[match]->Fill(mcConvR_, tracks[i]->normalizedChi2());
1848  }
1849 
1850  bool associated = false;
1851  float mcConvPt_ = -99999999.0;
1852  // float mcPhi= 0; // unused
1853  float simPV_Z = 0;
1854  for (std::vector<PhotonMCTruth>::const_iterator mcPho = mcPhotons.begin(); mcPho != mcPhotons.end(); mcPho++) {
1855  mcConvPt_ = (*mcPho).fourMomentum().et();
1856  float mcPhi = (*mcPho).fourMomentum().phi();
1857  simPV_Z = (*mcPho).primaryVertex().z();
1858  mcPhi_ = phiNormalization(mcPhi);
1859  mcEta_ = (*mcPho).fourMomentum().pseudoRapidity();
1860  mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z());
1861  mcConvR_ = (*mcPho).vertex().perp();
1862  mcConvX_ = (*mcPho).vertex().x();
1863  mcConvY_ = (*mcPho).vertex().y();
1864  mcConvZ_ = (*mcPho).vertex().z();
1865  mcConvEta_ = (*mcPho).vertex().eta();
1866  mcConvPhi_ = (*mcPho).vertex().phi();
1867  if (fabs(mcEta_) > END_HI)
1868  continue;
1869  if (mcConvPt_ < minPhoPtForPurity)
1870  continue;
1871  if (fabs(mcEta_) > maxPhoEtaForPurity)
1872  continue;
1873  if (fabs(mcConvZ_) > maxPhoZForPurity)
1874  continue;
1875  if (mcConvR_ > maxPhoRForEffic)
1876  continue;
1877 
1878  if ((*mcPho).isAConversion() != 1)
1879  continue;
1880  if (!((fabs(mcEta_) <= BARL && mcConvR_ < 85) ||
1881  (fabs(mcEta_) > BARL && fabs(mcEta_) <= END_HI && fabs((*mcPho).vertex().z()) < 210)))
1882  continue;
1883 
1884  theConvTP_.clear();
1885  for (const TrackingParticleRef tp : tpForFakeRate) {
1886  if (fabs(tp->vx() - (*mcPho).vertex().x()) < 0.0001 && fabs(tp->vy() - (*mcPho).vertex().y()) < 0.0001 &&
1887  fabs(tp->vz() - (*mcPho).vertex().z()) < 0.0001) {
1888  theConvTP_.push_back(tp);
1889  }
1890  }
1891 
1892  if (theConvTP_.size() < 2)
1893  continue;
1894 
1895  //associated = false;
1896  reco::RecoToSimCollection const& p1 = theTrackAssociator->associateRecoToSim(tc1, theConvTP_);
1897  reco::RecoToSimCollection const& p2 = theTrackAssociator->associateRecoToSim(tc2, theConvTP_);
1898 
1899  auto itP1 = p1.find(tk1);
1900  auto itP2 = p2.find(tk2);
1901  bool good = (itP1 != p1.end()) and (not itP1->val.empty()) and (itP2 != p2.end()) and (not itP2->val.empty());
1902  if (not good) {
1903  itP1 = p1.find(tk2);
1904  itP2 = p2.find(tk1);
1905  good = (itP1 != p1.end()) and (not itP1->val.empty()) and (itP2 != p2.end()) and (not itP2->val.empty());
1906  }
1907 
1908  if (good) {
1909  std::vector<std::pair<TrackingParticleRef, double> > const& tp1 = itP1->val;
1910  std::vector<std::pair<TrackingParticleRef, double> > const& tp2 = itP2->val;
1911  TrackingParticleRef tpr1 = tp1.front().first;
1912  TrackingParticleRef tpr2 = tp2.front().first;
1913  if (abs(tpr1->pdgId()) == 11 && abs(tpr2->pdgId()) == 11 && tpr1->pdgId() * tpr2->pdgId() < 0) {
1914  if ((tpr1->parentVertex()->sourceTracks_end() - tpr1->parentVertex()->sourceTracks_begin() == 1) &&
1915  (tpr2->parentVertex()->sourceTracks_end() - tpr2->parentVertex()->sourceTracks_begin() == 1)) {
1916  if (tpr1->parentVertex().key() == tpr2->parentVertex().key() &&
1917  ((*tpr1->parentVertex()->sourceTracks_begin())->pdgId() == 22)) {
1918  mcConvR_ = sqrt(tpr1->parentVertex()->position().Perp2());
1919  mcConvZ_ = tpr1->parentVertex()->position().z();
1920  mcConvX_ = tpr1->parentVertex()->position().x();
1921  mcConvY_ = tpr1->parentVertex()->position().y();
1922  mcConvEta_ = tpr1->parentVertex()->position().eta();
1923  mcConvPhi_ = tpr1->parentVertex()->position().phi();
1924  mcConvPt_ = sqrt((*tpr1->parentVertex()->sourceTracks_begin())->momentum().Perp2());
1925  //std::cout << " Reco to Sim mcconvpt " << mcConvPt_ << std::endl;
1926  //cout << "associated track1 to " << tpr1->pdgId() << " with p=" << tpr1->p4() << " with pT=" << tpr1->pt() << endl;
1927  //cout << "associated track2 to " << tpr2->pdgId() << " with p=" << tpr2->p4() << " with pT=" << tpr2->pt() << endl;
1928  associated = true;
1929  break;
1930  }
1931  }
1932  }
1933  }
1934 
1935  } // end loop on sim photons
1936 
1937  if (false) {
1938  theConvTP_.clear();
1939  theConvTP_ = tpForFakeRate;
1940  reco::RecoToSimCollection p1incl = theTrackAssociator->associateRecoToSim(tc1, theConvTP_);
1941  reco::RecoToSimCollection p2incl = theTrackAssociator->associateRecoToSim(tc2, theConvTP_);
1942 
1943  for (std::vector<PhotonMCTruth>::const_iterator mcPho = mcPhotons.begin(); mcPho != mcPhotons.end(); mcPho++) {
1944  mcConvPt_ = (*mcPho).fourMomentum().et();
1945  float mcPhi = (*mcPho).fourMomentum().phi();
1946  simPV_Z = (*mcPho).primaryVertex().z();
1947  mcPhi_ = phiNormalization(mcPhi);
1948  mcEta_ = (*mcPho).fourMomentum().pseudoRapidity();
1949  mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z());
1950  mcConvR_ = (*mcPho).vertex().perp();
1951  mcConvX_ = (*mcPho).vertex().x();
1952  mcConvY_ = (*mcPho).vertex().y();
1953  mcConvZ_ = (*mcPho).vertex().z();
1954  mcConvEta_ = (*mcPho).vertex().eta();
1955  mcConvPhi_ = (*mcPho).vertex().phi();
1956  if (fabs(mcEta_) > END_HI)
1957  continue;
1958  if (mcConvPt_ < minPhoPtForPurity)
1959  continue;
1960  if (fabs(mcEta_) > maxPhoEtaForPurity)
1961  continue;
1962  if (fabs(mcConvZ_) > maxPhoZForPurity)
1963  continue;
1964  if (mcConvR_ > maxPhoRForEffic)
1965  continue;
1966 
1967  if ((*mcPho).isAConversion() != 1)
1968  continue;
1969  if (!((fabs(mcEta_) <= BARL && mcConvR_ < 85) ||
1970  (fabs(mcEta_) > BARL && fabs(mcEta_) <= END_HI && fabs((*mcPho).vertex().z()) < 210)))
1971  continue;
1972 
1973  theConvTP_.clear();
1974  for (TrackingParticleRef tp : tpForFakeRate) {
1975  if (fabs(tp->vx() - (*mcPho).vertex().x()) < 0.0001 && fabs(tp->vy() - (*mcPho).vertex().y()) < 0.0001 &&
1976  fabs(tp->vz() - (*mcPho).vertex().z()) < 0.0001) {
1977  theConvTP_.push_back(tp);
1978  }
1979  }
1980 
1981  if (theConvTP_.size() < 2)
1982  continue;
1983 
1984  //associated = false;
1985  reco::RecoToSimCollection p1 = theTrackAssociator->associateRecoToSim(tc1, theConvTP_);
1986  reco::RecoToSimCollection p2 = theTrackAssociator->associateRecoToSim(tc2, theConvTP_);
1987 
1988  if ((!p1incl.empty() && !p2incl.empty()) && (!p1.empty() || !p2.empty())) { // associated = true;
1989  try {
1990  std::vector<std::pair<TrackingParticleRef, double> > tp1 = p1incl[tk1];
1991  std::vector<std::pair<TrackingParticleRef, double> > tp2 = p2incl[tk2];
1992  if (!(!tp1.empty() && !tp2.empty())) {
1993  tp1 = p1[tk2];
1994  tp2 = p2[tk1];
1995  }
1996  if (!tp1.empty() && !tp2.empty()) {
1997  TrackingParticleRef tpr1 = tp1.front().first;
1998  TrackingParticleRef tpr2 = tp2.front().first;
1999  if (abs(tpr1->pdgId()) == 11 && abs(tpr2->pdgId()) == 11 && tpr1->pdgId() * tpr2->pdgId() < 0) {
2000  if (((tpr1->parentVertex()->sourceTracks_end() - tpr1->parentVertex()->sourceTracks_begin() >= 1) &&
2001  (*tpr1->parentVertex()->sourceTracks_begin())->pdgId() == 22) &&
2002  ((tpr2->parentVertex()->sourceTracks_end() - tpr2->parentVertex()->sourceTracks_begin() >= 1) &&
2003  (*tpr2->parentVertex()->sourceTracks_begin())->pdgId() == 22)) {
2004  // if ( fabs(tpr1->vx() - tpr2->vx()) < 0.1 && fabs(tpr1->vy() - tpr2->vy()) < 0.1 && fabs(tpr1->vz() - tpr2->vz()) < 0.1) {
2005  //if (((*tpr1->parentVertex()->sourceTracks_begin())->pdgId()==22) || ((*tpr2->parentVertex()->sourceTracks_begin())->pdgId()==22)) {
2006  // mcConvR_ = sqrt(tpr1->parentVertex()->position().Perp2());
2007  // mcConvZ_ = tpr1->parentVertex()->position().z();
2008  // mcConvX_ = tpr1->parentVertex()->position().x();
2009  // mcConvY_ = tpr1->parentVertex()->position().y();
2010  // mcConvEta_ = tpr1->parentVertex()->position().eta();
2011  // mcConvPhi_ = tpr1->parentVertex()->position().phi();
2012  // mcConvPt_ = sqrt((*tpr1->parentVertex()->sourceTracks_begin())->momentum().Perp2());
2013  //std::cout << " Reco to Sim mcconvpt " << mcConvPt_ << std::endl;
2014  //cout << "associated track1 to " << tpr1->pdgId() << " with p=" << tpr1->p4() << " with pT=" << tpr1->pt() << endl;
2015  //cout << "associated track2 to " << tpr2->pdgId() << " with p=" << tpr2->p4() << " with pT=" << tpr2->pt() << endl;
2016  associated = true;
2017  break;
2018  //}
2019  //}
2020  }
2021  }
2022  }
2023  } catch (Exception const& event) {
2024  //cout << "do not continue: " << event.what() << endl;
2025  //continue;
2026  }
2027  }
2028  }
2029  }
2030 
2031  if (associated)
2032  match = 1;
2033  else
2034  match = 2;
2035 
2036  h_match_->Fill(float(match));
2038  if (match == 1)
2039  nRecConvAss_++;
2040  h_convEta_[match][0]->Fill(refittedMom.eta());
2041  h_convEta_[match][1]->Fill(refittedMom.eta());
2042  if (matchConvSC)
2043  h_convEtaMatchSC_[match][0]->Fill(refittedMom.eta());
2044  h_convPhi_[match][0]->Fill(refittedMom.phi());
2045  h_convR_[match][0]->Fill(sqrt(aConv.conversionVertex().position().perp2()));
2046  h_convZ_[match][0]->Fill(aConv.conversionVertex().position().z());
2047  h_convPt_[match][0]->Fill(sqrt(refittedMom.perp2()));
2048  h_invMass_[match][0]->Fill(invM);
2049  h_vtxChi2Prob_[match][0]->Fill(chi2Prob);
2050  h_DPhiTracksAtVtx_[match][0]->Fill(dPhiTracksAtVtx);
2051  h_DCotTracks_[match][0]->Fill(aConv.pairCotThetaSeparation());
2052  h_distMinAppTracks_[match][0]->Fill(aConv.distOfMinimumApproach());
2053  h_lxybs_[match][0]->Fill(lxy);
2054  h_maxNHitsBeforeVtx_[match][0]->Fill(maxNHitsBeforeVtx);
2055  h_leadNHitsBeforeVtx_[match][0]->Fill(leadNHitsBeforeVtx);
2056  h_trailNHitsBeforeVtx_[match][0]->Fill(trailNHitsBeforeVtx);
2057  h_sumNHitsBeforeVtx_[match][0]->Fill(sumNHitsBeforeVtx);
2058  h_deltaExpectedHitsInner_[match][0]->Fill(deltaExpectedHitsInner);
2059  h_leadExpectedHitsInner_[match][0]->Fill(leadExpectedHitsInner);
2060  h_maxDlClosestHitToVtx_[match][0]->Fill(maxDlClosestHitToVtx);
2061  h_maxDlClosestHitToVtxSig_[match][0]->Fill(maxDlClosestHitToVtxSig);
2062  h_nSharedHits_[match][0]->Fill(aConv.nSharedHits());
2063  if (matchConvSC) {
2064  //h_EoverPTracks_[match][0] ->Fill (aConv.EoverPrefittedTracks());
2065  h_EoverPTracks_[match][0]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
2066  h_convSCdPhi_[match][0]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
2067  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
2068  h_convSCdEta_[match][0]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
2069  }
2070  if (match == 1) {
2071  h2_photonPtRecVsPtSim_->Fill(mcConvPt_, sqrt(refittedMom.perp2()));
2072  h_convPtRes_[0]->Fill(sqrt(refittedMom.perp2()) / mcConvPt_);
2073  }
2074 
2075  if (phoIsInBarrel) {
2076  h_invMass_[match][1]->Fill(invM);
2077  h_vtxChi2Prob_[match][1]->Fill(chi2Prob);
2078  h_DPhiTracksAtVtx_[match][1]->Fill(dPhiTracksAtVtx);
2079  h_DCotTracks_[match][1]->Fill(aConv.pairCotThetaSeparation());
2080  h_distMinAppTracks_[match][1]->Fill(aConv.distOfMinimumApproach());
2081  h_lxybs_[match][1]->Fill(lxy);
2082  h_maxNHitsBeforeVtx_[match][1]->Fill(maxNHitsBeforeVtx);
2083  h_leadNHitsBeforeVtx_[match][1]->Fill(leadNHitsBeforeVtx);
2084  h_trailNHitsBeforeVtx_[match][1]->Fill(trailNHitsBeforeVtx);
2085  h_sumNHitsBeforeVtx_[match][1]->Fill(sumNHitsBeforeVtx);
2086  h_deltaExpectedHitsInner_[match][1]->Fill(deltaExpectedHitsInner);
2087  h_leadExpectedHitsInner_[match][1]->Fill(leadExpectedHitsInner);
2088  h_maxDlClosestHitToVtx_[match][1]->Fill(maxDlClosestHitToVtx);
2089  h_maxDlClosestHitToVtxSig_[match][1]->Fill(maxDlClosestHitToVtxSig);
2090  h_nSharedHits_[match][1]->Fill(aConv.nSharedHits());
2091  if (matchConvSC) {
2092  // h_EoverPTracks_[match][1] ->Fill (aConv.EoverPrefittedTracks());
2093  h_EoverPTracks_[match][1]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
2094  h_convSCdPhi_[match][1]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
2095  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
2096  h_convSCdEta_[match][1]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
2097  }
2098  if (match == 1)
2099  h_convPtRes_[1]->Fill(sqrt(refittedMom.perp2()) / mcConvPt_);
2100  }
2101 
2102  if (phoIsInEndcap) {
2103  h_invMass_[match][2]->Fill(invM);
2104  h_vtxChi2Prob_[match][2]->Fill(chi2Prob);
2105  h_DPhiTracksAtVtx_[match][2]->Fill(dPhiTracksAtVtx);
2106  h_DCotTracks_[match][2]->Fill(aConv.pairCotThetaSeparation());
2107  h_distMinAppTracks_[match][2]->Fill(aConv.distOfMinimumApproach());
2108  h_lxybs_[match][2]->Fill(lxy);
2109  h_maxNHitsBeforeVtx_[match][2]->Fill(maxNHitsBeforeVtx);
2110  h_leadNHitsBeforeVtx_[match][2]->Fill(leadNHitsBeforeVtx);
2111  h_trailNHitsBeforeVtx_[match][2]->Fill(trailNHitsBeforeVtx);
2112  h_sumNHitsBeforeVtx_[match][2]->Fill(sumNHitsBeforeVtx);
2113  h_deltaExpectedHitsInner_[match][2]->Fill(deltaExpectedHitsInner);
2114  h_leadExpectedHitsInner_[match][2]->Fill(leadExpectedHitsInner);
2115  h_maxDlClosestHitToVtx_[match][2]->Fill(maxDlClosestHitToVtx);
2116  h_maxDlClosestHitToVtxSig_[match][2]->Fill(maxDlClosestHitToVtxSig);
2117  h_nSharedHits_[match][2]->Fill(aConv.nSharedHits());
2118  if (matchConvSC) {
2119  // h_EoverPTracks_[match][2] ->Fill (aConv.EoverPrefittedTracks());
2120  h_EoverPTracks_[match][2]->Fill(iMatchingSC->superCluster()->energy() / sqrt(refittedMom.mag2()));
2121  h_convSCdPhi_[match][2]->Fill(iMatchingSC->superCluster()->position().phi() - refittedMom.phi());
2122  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(), aConv.zOfPrimaryVertexFromTracks());
2123  h_convSCdEta_[match][2]->Fill(iMatchingSC->superCluster()->position().eta() - ConvEta);
2124  }
2125  if (match == 1)
2126  h_convPtRes_[2]->Fill(sqrt(refittedMom.perp2()) / mcConvPt_);
2127  }
2128 
2129  if (match == 1) {
2130  h_convVtxdX_->Fill(aConv.conversionVertex().position().x() - mcConvX_);
2131  h_convVtxdY_->Fill(aConv.conversionVertex().position().y() - mcConvY_);
2132  h_convVtxdZ_->Fill(aConv.conversionVertex().position().z() - mcConvZ_);
2133  h_convVtxdR_->Fill(sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
2134  h_convVtxdPhi_->Fill(aConv.conversionVertex().position().phi() - mcConvPhi_);
2135  h_convVtxdEta_->Fill(aConv.conversionVertex().position().eta() - mcConvEta_);
2136  h2_convVtxdRVsR_->Fill(mcConvR_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
2137  h2_convVtxdRVsEta_->Fill(mcEta_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
2138  p_convVtxdRVsR_->Fill(mcConvR_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
2139  p_convVtxdRVsEta_->Fill(mcEta_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
2140  p_convVtxdXVsX_->Fill(mcConvX_, aConv.conversionVertex().position().x() - mcConvX_);
2141  p_convVtxdYVsY_->Fill(mcConvY_, aConv.conversionVertex().position().y() - mcConvY_);
2142  p_convVtxdZVsZ_->Fill(mcConvZ_, aConv.conversionVertex().position().z() - mcConvZ_);
2143  p_convVtxdZVsR_->Fill(mcConvR_, aConv.conversionVertex().position().z() - mcConvZ_);
2144 
2145  float dR = sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_;
2146  float dZ = aConv.conversionVertex().position().z() - mcConvZ_;
2147  p2_convVtxdRVsRZ_->Fill(mcConvZ_, mcConvR_, dR);
2148  p2_convVtxdZVsRZ_->Fill(mcConvZ_, mcConvR_, dZ);
2149 
2150  h2_convVtxRrecVsTrue_->Fill(mcConvR_, sqrt(aConv.conversionVertex().position().perp2()));
2151 
2152  h_zPVFromTracks_[match]->Fill(aConv.zOfPrimaryVertexFromTracks());
2153  h_dzPVFromTracks_[match]->Fill(aConv.zOfPrimaryVertexFromTracks() - simPV_Z);
2154  h2_dzPVVsR_->Fill(mcConvR_, aConv.zOfPrimaryVertexFromTracks() - simPV_Z);
2155  p_dzPVVsR_->Fill(mcConvR_, aConv.zOfPrimaryVertexFromTracks() - simPV_Z);
2156 
2157  if (phoIsInBarrel) {
2158  h_convVtxdX_barrel_->Fill(aConv.conversionVertex().position().x() - mcConvX_);
2159  h_convVtxdY_barrel_->Fill(aConv.conversionVertex().position().y() - mcConvY_);
2160  h_convVtxdZ_barrel_->Fill(aConv.conversionVertex().position().z() - mcConvZ_);
2161  h_convVtxdR_barrel_->Fill(sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
2162  }
2163  if (phoIsInEndcap) {
2164  h_convVtxdX_endcap_->Fill(aConv.conversionVertex().position().x() - mcConvX_);
2165  h_convVtxdY_endcap_->Fill(aConv.conversionVertex().position().y() - mcConvY_);
2166  h_convVtxdZ_endcap_->Fill(aConv.conversionVertex().position().z() - mcConvZ_);
2167  h_convVtxdR_endcap_->Fill(sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
2168  }
2169  }
2170 
2172  for (unsigned int i = 0; i < tracks.size(); i++) {
2173  //std::cout << " Loop over tracks pt " << tracks[i]->pt() << std::endl;
2174  RefToBase<reco::Track> tfrb(aConv.tracks()[i]);
2175  itAss = myAss.find(tfrb.get());
2176 
2177  nHitsVsEta_[match]->Fill(mcEta_, float(tracks[i]->numberOfValidHits()));
2178  nHitsVsR_[match]->Fill(mcConvR_, float(tracks[i]->numberOfValidHits()));
2179  p_nHitsVsEta_[match]->Fill(mcEta_, float(tracks[i]->numberOfValidHits()) - 0.0001);
2180  p_nHitsVsR_[match]->Fill(mcConvR_, float(tracks[i]->numberOfValidHits()) - 0.0001);
2181  h_tkChi2_[match]->Fill(tracks[i]->normalizedChi2());
2182  h_tkChi2Large_[match]->Fill(tracks[i]->normalizedChi2());
2183  h2_Chi2VsEta_[match]->Fill(mcEta_, tracks[i]->normalizedChi2());
2184  h2_Chi2VsR_[match]->Fill(mcConvR_, tracks[i]->normalizedChi2());
2185  p_Chi2VsEta_[match]->Fill(mcEta_, tracks[i]->normalizedChi2());
2186  p_Chi2VsR_[match]->Fill(mcConvR_, tracks[i]->normalizedChi2());
2187  double d0;
2188  if (valid_pvtx) {
2189  d0 = -tracks[i]->dxy(the_pvtx.position());
2190  } else {
2191  d0 = tracks[i]->d0();
2192  }
2193  h_TkD0_[match]->Fill(d0 * tracks[i]->charge());
2194  h_nHitsBeforeVtx_[match]->Fill(aConv.nHitsBeforeVtx().size() > 1 ? aConv.nHitsBeforeVtx().at(i) : 0);
2195  h_dlClosestHitToVtx_[match]->Fill(aConv.dlClosestHitToVtx().size() > 1 ? aConv.dlClosestHitToVtx().at(i).value()
2196  : 0);
2197  h_dlClosestHitToVtxSig_[match]->Fill(aConv.dlClosestHitToVtx().size() > 1
2198  ? aConv.dlClosestHitToVtx().at(i).value() /
2199  aConv.dlClosestHitToVtx().at(i).error()
2200  : 0);
2201 
2202  if (itAss == myAss.end())
2203  continue;
2204  reco::Track refTrack = aConv.conversionVertex().refittedTracks()[i];
2205 
2206  float simPt = sqrt(((*itAss).second)->momentum().perp2());
2207  float recPt = refTrack.pt();
2208  float ptres = recPt - simPt;
2209  //float pterror = aConv.tracks()[i]->ptError();
2210  float pterror = aConv.conversionVertex().refittedTracks()[i].ptError();
2211  h2_PtRecVsPtSim_[0]->Fill(simPt, recPt);
2212  h_TkPtPull_[0]->Fill(ptres / pterror);
2213  h2_TkPtPull_[0]->Fill(mcEta_, ptres / pterror);
2214 
2215  if (phoIsInBarrel) {
2216  h_TkPtPull_[1]->Fill(ptres / pterror);
2217  h2_PtRecVsPtSim_[1]->Fill(simPt, recPt);
2218  }
2219  if (phoIsInEndcap) {
2220  h_TkPtPull_[2]->Fill(ptres / pterror);
2221  h2_PtRecVsPtSim_[2]->Fill(simPt, recPt);
2222  }
2223  } // end loop over track
2224 
2225  } // loop over reco conversions
2226 
2227  h_nConv_[0][0]->Fill(float(nRecConv_));
2228  h_nConv_[1][0]->Fill(float(nRecConvAss_));
2229 }

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

◆ bookHistograms()

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

Histograms for efficiencies

Denominators

zooms

Definition at line 146 of file TkConvValidator.cc.

146  {
147  nEvt_ = 0;
148  nEntry_ = 0;
149  nRecConv_ = 0;
150  nRecConvAss_ = 0;
151  nRecConvAssWithEcal_ = 0;
152 
153  nInvalidPCA_ = 0;
154 
155  dbe_ = nullptr;
156  dbe_ = edm::Service<DQMStore>().operator->();
157 
158  double etMin = parameters_.getParameter<double>("etMin");
159  double etMax = parameters_.getParameter<double>("etMax");
160  int etBin = parameters_.getParameter<int>("etBin");
161 
162  double resMin = parameters_.getParameter<double>("resMin");
163  double resMax = parameters_.getParameter<double>("resMax");
164  int resBin = parameters_.getParameter<int>("resBin");
165 
166  double etaMin = parameters_.getParameter<double>("etaMin");
167  double etaMax = parameters_.getParameter<double>("etaMax");
168  int etaBin = parameters_.getParameter<int>("etaBin");
169  int etaBin2 = parameters_.getParameter<int>("etaBin2");
170 
171  double phiMin = parameters_.getParameter<double>("phiMin");
172  double phiMax = parameters_.getParameter<double>("phiMax");
173  int phiBin = parameters_.getParameter<int>("phiBin");
174 
175  double rMin = parameters_.getParameter<double>("rMin");
176  double rMax = parameters_.getParameter<double>("rMax");
177  int rBin = parameters_.getParameter<int>("rBin");
178 
179  double zMin = parameters_.getParameter<double>("zMin");
180  double zMax = parameters_.getParameter<double>("zMax");
181  int zBin = parameters_.getParameter<int>("zBin");
182 
183  double dPhiTracksMin = parameters_.getParameter<double>("dPhiTracksMin");
184  double dPhiTracksMax = parameters_.getParameter<double>("dPhiTracksMax");
185  int dPhiTracksBin = parameters_.getParameter<int>("dPhiTracksBin");
186 
187  double eoverpMin = parameters_.getParameter<double>("eoverpMin");
188  double eoverpMax = parameters_.getParameter<double>("eoverpMax");
189  int eoverpBin = parameters_.getParameter<int>("eoverpBin");
190 
191  // double dEtaTracksMin = parameters_.getParameter<double>("dEtaTracksMin"); // unused
192  // double dEtaTracksMax = parameters_.getParameter<double>("dEtaTracksMax"); // unused
193  // int dEtaTracksBin = parameters_.getParameter<int>("dEtaTracksBin"); // unused
194 
195  double dCotTracksMin = parameters_.getParameter<double>("dCotTracksMin");
196  double dCotTracksMax = parameters_.getParameter<double>("dCotTracksMax");
197  int dCotTracksBin = parameters_.getParameter<int>("dCotTracksBin");
198 
199  double chi2Min = parameters_.getParameter<double>("chi2Min");
200  double chi2Max = parameters_.getParameter<double>("chi2Max");
201 
202  double rMinForXray = parameters_.getParameter<double>("rMinForXray");
203  double rMaxForXray = parameters_.getParameter<double>("rMaxForXray");
204  int rBinForXray = parameters_.getParameter<int>("rBinForXray");
205  double zMinForXray = parameters_.getParameter<double>("zMinForXray");
206  double zMaxForXray = parameters_.getParameter<double>("zMaxForXray");
207  int zBinForXray = parameters_.getParameter<int>("zBinForXray");
208  int zBin2ForXray = parameters_.getParameter<int>("zBin2ForXray");
209 
210  minPhoPtForEffic = parameters_.getParameter<double>("minPhoPtForEffic");
211  maxPhoEtaForEffic = parameters_.getParameter<double>("maxPhoEtaForEffic");
212  maxPhoZForEffic = parameters_.getParameter<double>("maxPhoZForEffic");
213  maxPhoRForEffic = parameters_.getParameter<double>("maxPhoRForEffic");
214  minPhoPtForPurity = parameters_.getParameter<double>("minPhoPtForPurity");
215  maxPhoEtaForPurity = parameters_.getParameter<double>("maxPhoEtaForPurity");
216  maxPhoZForPurity = parameters_.getParameter<double>("maxPhoZForPurity");
217  maxPhoRForPurity = parameters_.getParameter<double>("maxPhoRForPurity");
218 
219  if (dbe_) {
221  // SC from reco photons
222 
223  //TString simfolder = TString(
224  std::string simpath = dqmpath_ + "SimulationInfo";
225  iBooker.setCurrentFolder(simpath);
226  //
227  // simulation information about conversions
229  std::string histname = "nOfSimConversions";
230  h_nSimConv_[0] = iBooker.book1D(histname, "# of Sim conversions per event ", 20, -0.5, 19.5);
232  histname = "h_AllSimConvEta";
233  h_AllSimConv_[0] = iBooker.book1D(histname, " All conversions: simulated #eta", etaBin2, etaMin, etaMax);
234  histname = "h_AllSimConvPhi";
235  h_AllSimConv_[1] = iBooker.book1D(histname, " All conversions: simulated #phi", phiBin, phiMin, phiMax);
236  histname = "h_AllSimConvR";
237  h_AllSimConv_[2] = iBooker.book1D(histname, " All conversions: simulated R", rBin, rMin, rMax);
238  histname = "h_AllSimConvZ";
239  h_AllSimConv_[3] = iBooker.book1D(histname, " All conversions: simulated Z", zBin, zMin, zMax);
240  histname = "h_AllSimConvEt";
241  h_AllSimConv_[4] = iBooker.book1D(histname, " All conversions: simulated Et", etBin, etMin, etMax);
242  //
243  histname = "nOfVisSimConversions";
244  h_nSimConv_[1] = iBooker.book1D(histname, "# of Sim conversions per event ", 20, -0.5, 19.5);
245  histname = "h_VisSimConvEta";
246  h_VisSimConv_[0] = iBooker.book1D(histname, " All vis conversions: simulated #eta", etaBin2, etaMin, etaMax);
247  histname = "h_VisSimConvPhi";
248  h_VisSimConv_[1] = iBooker.book1D(histname, " All vis conversions: simulated #phi", phiBin, phiMin, phiMax);
249  histname = "h_VisSimConvR";
250  h_VisSimConv_[2] = iBooker.book1D(histname, " All vis conversions: simulated R", rBin, rMin, rMax);
251  histname = "h_VisSimConvZ";
252  h_VisSimConv_[3] = iBooker.book1D(histname, " All vis conversions: simulated Z", zBin, zMin, zMax);
253  histname = "h_VisSimConvEt";
254  h_VisSimConv_[4] = iBooker.book1D(histname, " All vis conversions: simulated Et", etBin, etMin, etMax);
255 
256  //
257  histname = "h_SimConvTwoMTracksEta";
258  h_SimConvTwoMTracks_[0] = iBooker.book1D(
259  histname, " All vis conversions with 2 reco-matching tracks: simulated #eta", etaBin2, etaMin, etaMax);
260  histname = "h_SimConvTwoMTracksPhi";
261  h_SimConvTwoMTracks_[1] = iBooker.book1D(
262  histname, " All vis conversions with 2 reco-matching tracks: simulated #phi", phiBin, phiMin, phiMax);
263  histname = "h_SimConvTwoMTracksR";
264  h_SimConvTwoMTracks_[2] =
265  iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated R", rBin, rMin, rMax);
266  histname = "h_SimConvTwoMTracksZ";
267  h_SimConvTwoMTracks_[3] =
268  iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated Z", zBin, zMin, zMax);
269  histname = "h_SimConvTwoMTracksEt";
270  h_SimConvTwoMTracks_[4] =
271  iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated Et", etBin, etMin, etMax);
272  //
273  histname = "h_SimConvTwoTracksEta";
274  h_SimConvTwoTracks_[0] =
275  iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated #eta", etaBin2, etaMin, etaMax);
276  histname = "h_SimConvTwoTracksPhi";
277  h_SimConvTwoTracks_[1] =
278  iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated #phi", phiBin, phiMin, phiMax);
279  histname = "h_SimConvTwoTracksR";
280  h_SimConvTwoTracks_[2] =
281  iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated R", rBin, rMin, rMax);
282  histname = "h_SimConvTwoTracksZ";
283  h_SimConvTwoTracks_[3] =
284  iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated Z", zBin, zMin, zMax);
285  histname = "h_SimConvTwoTracksEt";
286  h_SimConvTwoTracks_[4] =
287  iBooker.book1D(histname, " All vis conversions with 2 reco tracks: simulated Et", etBin, etMin, etMax);
288  //
289  histname = "h_SimConvTwoMTracksEtaAndVtxPGT0";
290  h_SimConvTwoMTracksAndVtxPGT0_[0] = iBooker.book1D(
291  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated #eta", etaBin2, etaMin, etaMax);
292  histname = "h_SimConvTwoMTracksPhiAndVtxPGT0";
293  h_SimConvTwoMTracksAndVtxPGT0_[1] = iBooker.book1D(
294  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated #phi", phiBin, phiMin, phiMax);
295  histname = "h_SimConvTwoMTracksRAndVtxPGT0";
296  h_SimConvTwoMTracksAndVtxPGT0_[2] = iBooker.book1D(
297  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated R", rBin, rMin, rMax);
298  histname = "h_SimConvTwoMTracksZAndVtxPGT0";
299  h_SimConvTwoMTracksAndVtxPGT0_[3] = iBooker.book1D(
300  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated Z", zBin, zMin, zMax);
301  histname = "h_SimConvTwoMTracksEtAndVtxPGT0";
302  h_SimConvTwoMTracksAndVtxPGT0_[4] = iBooker.book1D(
303  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated Et", etBin, etMin, etMax);
304 
305  //
306  histname = "h_SimConvTwoMTracksEtaAndVtxPGT0005";
307  h_SimConvTwoMTracksAndVtxPGT0005_[0] = iBooker.book1D(
308  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated #eta", etaBin2, etaMin, etaMax);
309  histname = "h_SimConvTwoMTracksPhiAndVtxPGT0005";
310  h_SimConvTwoMTracksAndVtxPGT0005_[1] = iBooker.book1D(
311  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated #phi", phiBin, phiMin, phiMax);
312  histname = "h_SimConvTwoMTracksRAndVtxPGT0005";
313  h_SimConvTwoMTracksAndVtxPGT0005_[2] = iBooker.book1D(
314  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated R", rBin, rMin, rMax);
315  histname = "h_SimConvTwoMTracksZAndVtxPGT0005";
316  h_SimConvTwoMTracksAndVtxPGT0005_[3] = iBooker.book1D(
317  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated Z", zBin, zMin, zMax);
318  histname = "h_SimConvTwoMTracksEtAndVtxPGT0005";
319  h_SimConvTwoMTracksAndVtxPGT0005_[4] = iBooker.book1D(
320  histname, " All vis conversions with 2 reco-matching tracks + vertex: simulated Et", etBin, etMin, etMax);
321 
322  histname = "h_SimRecConvTwoMTracksEta";
323  h_SimRecConvTwoMTracks_[0] = iBooker.book1D(
324  histname, " All vis conversions with 2 reco-matching tracks: simulated #eta", etaBin2, etaMin, etaMax);
325  histname = "h_SimRecConvTwoMTracksPhi";
326  h_SimRecConvTwoMTracks_[1] = iBooker.book1D(
327  histname, " All vis conversions with 2 reco-matching tracks: simulated #phi", phiBin, phiMin, phiMax);
328  histname = "h_SimRecConvTwoMTracksR";
329  h_SimRecConvTwoMTracks_[2] =
330  iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated R", rBin, rMin, rMax);
331  histname = "h_SimRecConvTwoMTracksZ";
332  h_SimRecConvTwoMTracks_[3] =
333  iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated Z", zBin, zMin, zMax);
334  histname = "h_SimRecConvTwoMTracksEt";
335  h_SimRecConvTwoMTracks_[4] =
336  iBooker.book1D(histname, " All vis conversions with 2 reco-matching tracks: simulated Et", etBin, etMin, etMax);
337  //
338 
339  h_SimConvEtaPix_[0] = iBooker.book1D("simConvEtaPix", " sim converted Photon Eta: Pix ", etaBin, etaMin, etaMax);
340  h_simTkPt_ = iBooker.book1D("simTkPt", "Sim conversion tracks pt ", etBin * 3, 0., etMax);
341  h_simTkEta_ = iBooker.book1D("simTkEta", "Sim conversion tracks eta ", etaBin, etaMin, etaMax);
342 
343  h_simConvVtxRvsZ_[0] = iBooker.book2D("simConvVtxRvsZAll",
344  " Photon Sim conversion vtx position",
345  zBinForXray,
346  zMinForXray,
347  zMaxForXray,
348  rBinForXray,
349  rMinForXray,
350  rMaxForXray);
351  h_simConvVtxRvsZ_[1] = iBooker.book2D("simConvVtxRvsZBarrel",
352  " Photon Sim conversion vtx position",
353  zBinForXray,
354  zMinForXray,
355  zMaxForXray,
356  rBinForXray,
357  rMinForXray,
358  rMaxForXray);
359  h_simConvVtxRvsZ_[2] = iBooker.book2D("simConvVtxRvsZEndcap",
360  " Photon Sim conversion vtx position",
361  zBin2ForXray,
362  zMinForXray,
363  zMaxForXray,
364  rBinForXray,
365  rMinForXray,
366  rMaxForXray);
367  h_simConvVtxRvsZ_[3] = iBooker.book2D("simConvVtxRvsZBarrel2",
368  " Photon Sim conversion vtx position when reco R<4cm",
369  zBinForXray,
370  zMinForXray,
371  zMaxForXray,
372  rBinForXray,
373  rMinForXray,
374  rMaxForXray);
375  h_simConvVtxYvsX_ = iBooker.book2D(
376  "simConvVtxYvsXTrkBarrel", " Photon Sim conversion vtx position, (x,y) eta<1 ", 100, -80., 80., 100, -80., 80.);
377 
378  std::string convpath = dqmpath_ + "ConversionInfo";
379  iBooker.setCurrentFolder(convpath);
380 
381  histname = "nConv";
382  h_nConv_[0][0] = iBooker.book1D(
383  histname + "All", "Number Of Conversions per isolated candidates per events: All Ecal ", 10, -0.5, 9.5);
384  h_nConv_[0][1] = iBooker.book1D(
385  histname + "Barrel", "Number Of Conversions per isolated candidates per events: Ecal Barrel ", 10, -0.5, 9.5);
386  h_nConv_[0][2] = iBooker.book1D(
387  histname + "Endcap", "Number Of Conversions per isolated candidates per events: Ecal Endcap ", 10, -0.5, 9.5);
388  h_nConv_[1][0] = iBooker.book1D(histname + "All_Ass",
389  "Number Of associated Conversions per isolated candidates per events: All Ecal ",
390  10,
391  -0.5,
392  9.5);
393 
394  h_convEta_[0][0] = iBooker.book1D("convEta", " converted Photon Eta ", etaBin, etaMin, etaMax);
395  h_convEtaMatchSC_[0][0] =
396  iBooker.book1D("convEtaMatchSC", " converted Photon Eta when SC is matched ", etaBin, etaMin, etaMax);
397  h_convEta2_[0][0] = iBooker.book1D("convEta2", " converted Photon Eta ", etaBin2, etaMin, etaMax);
398  h_convPhi_[0][0] = iBooker.book1D("convPhi", " converted Photon Phi ", phiBin, phiMin, phiMax);
399  h_convR_[0][0] = iBooker.book1D("convR", " converted photon R", rBin, rMin, rMax);
400  h_convZ_[0][0] = iBooker.book1D("convZ", " converted photon Z", zBin, zMin, zMax);
401  h_convPt_[0][0] = iBooker.book1D("convPt", " conversions Transverse Energy: all eta ", etBin, etMin, etMax);
402 
403  h_convEta_[1][0] = iBooker.book1D("convEtaAss2", " Matched converted Photon Eta ", etaBin2, etaMin, etaMax);
404  h_convEta_[1][1] = iBooker.book1D("convEtaAss", " Matched converted Photon Eta ", etaBin, etaMin, etaMax);
405  h_convEtaMatchSC_[1][0] =
406  iBooker.book1D("convEtaMatchSCAss", " converted Photon Eta when SC is matched ", etaBin, etaMin, etaMax);
407  h_convPhi_[1][0] = iBooker.book1D("convPhiAss", " Matched converted Photon Phi ", phiBin, phiMin, phiMax);
408  h_convR_[1][0] = iBooker.book1D("convRAss", " Matched converted photon R", rBin, rMin, rMax);
409  h_convZ_[1][0] = iBooker.book1D("convZAss", " Matched converted photon Z", zBin, zMin, zMax);
410  h_convPt_[1][0] =
411  iBooker.book1D("convPtAss", "Matched conversions Transverse Energy: all eta ", etBin, etMin, etMax);
412 
413  h_convEta_[2][0] = iBooker.book1D("convEtaFake2", " Fake converted Photon Eta ", etaBin2, etaMin, etaMax);
414  h_convEta_[2][1] = iBooker.book1D("convEtaFake", " Fake converted Photon Eta ", etaBin, etaMin, etaMax);
415  h_convEtaMatchSC_[2][0] =
416  iBooker.book1D("convEtaMatchSCFake", " converted Photon Eta when SC is matched ", etaBin, etaMin, etaMax);
417  h_convPhi_[2][0] = iBooker.book1D("convPhiFake", " Fake converted Photon Phi ", phiBin, phiMin, phiMax);
418  h_convR_[2][0] = iBooker.book1D("convRFake", " Fake converted photon R", rBin, rMin, rMax);
419  h_convZ_[2][0] = iBooker.book1D("convZFake", " Fake converted photon Z", zBin, zMin, zMax);
420  h_convPt_[2][0] = iBooker.book1D("convPtFake", "Fake conversions Transverse Energy: all eta ", etBin, etMin, etMax);
421 
422  h_convRplot_ = iBooker.book1D("convRplot", " converted photon R", 600, 0., 120.);
423  h_convZplot_ = iBooker.book1D("convZplot", " converted photon Z", 320, -160., 160.);
424 
425  histname = "convSCdPhi";
426  h_convSCdPhi_[0][0] = iBooker.book1D(histname + "All", "dPhi between SC and conversion", 100, -0.1, 0.1);
427  h_convSCdPhi_[0][1] =
428  iBooker.book1D(histname + "Barrel", " dPhi between SC and conversion: Barrel", 100, -0.1, 0.1);
429  h_convSCdPhi_[0][2] =
430  iBooker.book1D(histname + "Endcap", " dPhi between SC and conversion: Endcap", 100, -0.1, 0.1);
431  h_convSCdPhi_[1][0] = iBooker.book1D(histname + "All_Ass", "dPhi between SC and conversion", 100, -0.1, 0.1);
432  h_convSCdPhi_[1][1] =
433  iBooker.book1D(histname + "Barrel_Ass", " dPhi between SC and conversion: Barrel", 100, -0.1, 0.1);
434  h_convSCdPhi_[1][2] =
435  iBooker.book1D(histname + "Endcap_Ass", " dPhi between SC and conversion: Endcap", 100, -0.1, 0.1);
436  h_convSCdPhi_[2][0] = iBooker.book1D(histname + "All_Fakes", "dPhi between SC and conversion", 100, -0.1, 0.1);
437  h_convSCdPhi_[2][1] =
438  iBooker.book1D(histname + "Barrel_Fakes", " dPhi between SC and conversion: Barrel", 100, -0.1, 0.1);
439  h_convSCdPhi_[2][2] =
440  iBooker.book1D(histname + "Endcap_Fakes", " dPhi between SC and conversion: Endcap", 100, -0.1, 0.1);
441  histname = "convSCdEta";
442  h_convSCdEta_[0][0] = iBooker.book1D(histname + "All", " dEta between SC and conversion", 100, -0.1, 0.1);
443  h_convSCdEta_[0][1] =
444  iBooker.book1D(histname + "Barrel", " dEta between SC and conversion: Barrel", 100, -0.1, 0.1);
445  h_convSCdEta_[0][2] =
446  iBooker.book1D(histname + "Endcap", " dEta between SC and conversion: Endcap", 100, -0.1, 0.1);
447  h_convSCdEta_[1][0] = iBooker.book1D(histname + "All_Ass", " dEta between SC and conversion", 100, -0.1, 0.1);
448  h_convSCdEta_[1][1] =
449  iBooker.book1D(histname + "Barrel_Ass", " dEta between SC and conversion: Barrel", 100, -0.1, 0.1);
450  h_convSCdEta_[1][2] =
451  iBooker.book1D(histname + "Endcap_Ass", " dEta between SC and conversion: Endcap", 100, -0.1, 0.1);
452  h_convSCdEta_[2][0] = iBooker.book1D(histname + "All_Fakes", " dEta between SC and conversion", 100, -0.1, 0.1);
453  h_convSCdEta_[2][1] =
454  iBooker.book1D(histname + "Barrel_Fakes", " dEta between SC and conversion: Barrel", 100, -0.1, 0.1);
455  h_convSCdEta_[2][2] =
456  iBooker.book1D(histname + "Endcap_Fakes", " dEta between SC and conversion: Endcap", 100, -0.1, 0.1);
457 
458  histname = "convPtRes";
459  h_convPtRes_[0] = iBooker.book1D(histname + "All", " Conversion Pt rec/true : All ecal ", resBin, resMin, resMax);
460  h_convPtRes_[1] = iBooker.book1D(histname + "Barrel", " Conversion Pt rec/true : Barrel ", resBin, resMin, resMax);
461  h_convPtRes_[2] = iBooker.book1D(histname + "Endcap", " Conversion Pt rec/true : Endcap ", resBin, resMin, resMax);
462 
463  histname = "hInvMass";
464  h_invMass_[0][0] = iBooker.book1D(
465  histname + "All_AllTracks", " Photons:Tracks from conversion: Pair invariant mass: all Ecal ", 100, 0., 1.5);
466  h_invMass_[0][1] = iBooker.book1D(histname + "Barrel_AllTracks",
467  " Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",
468  100,
469  0.,
470  1.5);
471  h_invMass_[0][2] = iBooker.book1D(histname + "Endcap_AllTracks",
472  " Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",
473  100,
474  0.,
475  1.5);
476  //
477  h_invMass_[1][0] = iBooker.book1D(
478  histname + "All_AssTracks", " Photons:Tracks from conversion: Pair invariant mass: all Ecal ", 100, 0., 1.5);
479  h_invMass_[1][1] = iBooker.book1D(histname + "Barrel_AssTracks",
480  " Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",
481  100,
482  0.,
483  1.5);
484  h_invMass_[1][2] = iBooker.book1D(histname + "Endcap_AssTracks",
485  " Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",
486  100,
487  0.,
488  1.5);
489  //
490  h_invMass_[2][0] = iBooker.book1D(
491  histname + "All_FakeTracks", " Photons:Tracks from conversion: Pair invariant mass: all Ecal ", 100, 0., 1.5);
492  h_invMass_[2][1] = iBooker.book1D(histname + "Barrel_FakeTracks",
493  " Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",
494  100,
495  0.,
496  1.5);
497  h_invMass_[2][2] = iBooker.book1D(histname + "Endcap_FaleTracks",
498  " Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",
499  100,
500  0.,
501  1.5);
502 
503  histname = "hDPhiTracksAtVtx";
504  h_DPhiTracksAtVtx_[0][0] = iBooker.book1D(histname + "All",
505  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",
506  dPhiTracksBin,
507  dPhiTracksMin,
508  dPhiTracksMax);
509  h_DPhiTracksAtVtx_[0][1] =
510  iBooker.book1D(histname + "Barrel",
511  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",
512  dPhiTracksBin,
513  dPhiTracksMin,
514  dPhiTracksMax);
515  h_DPhiTracksAtVtx_[0][2] =
516  iBooker.book1D(histname + "Endcap",
517  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",
518  dPhiTracksBin,
519  dPhiTracksMin,
520  dPhiTracksMax);
521  h_DPhiTracksAtVtx_[1][0] = iBooker.book1D(histname + "All_Ass",
522  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",
523  dPhiTracksBin,
524  dPhiTracksMin,
525  dPhiTracksMax);
526  h_DPhiTracksAtVtx_[1][1] =
527  iBooker.book1D(histname + "Barrel_Ass",
528  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",
529  dPhiTracksBin,
530  dPhiTracksMin,
531  dPhiTracksMax);
532  h_DPhiTracksAtVtx_[1][2] =
533  iBooker.book1D(histname + "Endcap_Ass",
534  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",
535  dPhiTracksBin,
536  dPhiTracksMin,
537  dPhiTracksMax);
538  h_DPhiTracksAtVtx_[2][0] = iBooker.book1D(histname + "All_Fakes",
539  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",
540  dPhiTracksBin,
541  dPhiTracksMin,
542  dPhiTracksMax);
543  h_DPhiTracksAtVtx_[2][1] =
544  iBooker.book1D(histname + "Barrel_Fakes",
545  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",
546  dPhiTracksBin,
547  dPhiTracksMin,
548  dPhiTracksMax);
549  h_DPhiTracksAtVtx_[2][2] =
550  iBooker.book1D(histname + "Endcap_Fakes",
551  " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",
552  dPhiTracksBin,
553  dPhiTracksMin,
554  dPhiTracksMax);
555 
556  histname = "hDPhiTracksAtVtxVsEta";
557  h2_DPhiTracksAtVtxVsEta_ = iBooker.book2D(histname + "All",
558  " Photons:Tracks from conversions: #delta#phi Tracks at vertex vs #eta",
559  etaBin2,
560  etaMin,
561  etaMax,
562  100,
563  -0.5,
564  0.5);
565  histname = "pDPhiTracksAtVtxVsEta";
566  p_DPhiTracksAtVtxVsEta_ =
567  iBooker.bookProfile(histname + "All",
568  " Photons:Tracks from conversions: #delta#phi Tracks at vertex vs #eta ",
569  etaBin2,
570  etaMin,
571  etaMax,
572  100,
573  -0.5,
574  0.5,
575  "");
576 
577  histname = "hDPhiTracksAtVtxVsR";
578  h2_DPhiTracksAtVtxVsR_ = iBooker.book2D(histname + "All",
579  " Photons:Tracks from conversions: #delta#phi Tracks at vertex vs R",
580  rBin,
581  rMin,
582  rMax,
583  100,
584  -0.5,
585  0.5);
586  histname = "pDPhiTracksAtVtxVsR";
587  p_DPhiTracksAtVtxVsR_ = iBooker.bookProfile(histname + "All",
588  " Photons:Tracks from conversions: #delta#phi Tracks at vertex vs R ",
589  rBin,
590  rMin,
591  rMax,
592  100,
593  -0.5,
594  0.5,
595  "");
596 
597  histname = "hDCotTracks";
598  h_DCotTracks_[0][0] = iBooker.book1D(histname + "All",
599  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",
600  dCotTracksBin,
601  dCotTracksMin,
602  dCotTracksMax);
603  h_DCotTracks_[0][1] = iBooker.book1D(histname + "Barrel",
604  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",
605  dCotTracksBin,
606  dCotTracksMin,
607  dCotTracksMax);
608  h_DCotTracks_[0][2] = iBooker.book1D(histname + "Endcap",
609  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",
610  dCotTracksBin,
611  dCotTracksMin,
612  dCotTracksMax);
613  h_DCotTracks_[1][0] = iBooker.book1D(histname + "All_Ass",
614  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",
615  dCotTracksBin,
616  dCotTracksMin,
617  dCotTracksMax);
618  h_DCotTracks_[1][1] = iBooker.book1D(histname + "Barrel_Ass",
619  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",
620  dCotTracksBin,
621  dCotTracksMin,
622  dCotTracksMax);
623  h_DCotTracks_[1][2] = iBooker.book1D(histname + "Endcap_Ass",
624  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",
625  dCotTracksBin,
626  dCotTracksMin,
627  dCotTracksMax);
628  h_DCotTracks_[2][0] = iBooker.book1D(histname + "All_Fakes",
629  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",
630  dCotTracksBin,
631  dCotTracksMin,
632  dCotTracksMax);
633  h_DCotTracks_[2][1] = iBooker.book1D(histname + "Barrel_Fakes",
634  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",
635  dCotTracksBin,
636  dCotTracksMin,
637  dCotTracksMax);
638  h_DCotTracks_[2][2] = iBooker.book1D(histname + "Endcap_Fakes",
639  " Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",
640  dCotTracksBin,
641  dCotTracksMin,
642  dCotTracksMax);
643 
644  histname = "hDCotTracksVsEta";
645  h2_DCotTracksVsEta_ = iBooker.book2D(histname + "All",
646  " Photons:Tracks from conversions: #delta cotg(#Theta) Tracks vs #eta",
647  etaBin2,
648  etaMin,
649  etaMax,
650  100,
651  -0.2,
652  0.2);
653  histname = "pDCotTracksVsEta";
654  p_DCotTracksVsEta_ = iBooker.bookProfile(histname + "All",
655  " Photons:Tracks from conversions: #delta cotg(#Theta) Tracks vs #eta ",
656  etaBin2,
657  etaMin,
658  etaMax,
659  100,
660  -0.2,
661  0.2,
662  "");
663 
664  histname = "hDCotTracksVsR";
665  h2_DCotTracksVsR_ = iBooker.book2D(histname + "All",
666  " Photons:Tracks from conversions: #delta cotg(#Theta) Tracks at vertex vs R",
667  rBin,
668  rMin,
669  rMax,
670  100,
671  -0.2,
672  0.2);
673  histname = "pDCotTracksVsR";
674  p_DCotTracksVsR_ =
675  iBooker.bookProfile(histname + "All",
676  " Photons:Tracks from conversions: #delta cotg(#Theta) Tracks at vertex vs R ",
677  rBin,
678  rMin,
679  rMax,
680  100,
681  -0.2,
682  0.2,
683  "");
684 
685  histname = "hDistMinAppTracks";
686  h_distMinAppTracks_[0][0] = iBooker.book1D(
687  histname + "All", " Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ", 120, -0.5, 1.0);
688  h_distMinAppTracks_[0][1] = iBooker.book1D(
689  histname + "Barrel", " Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ", 120, -0.5, 1.0);
690  h_distMinAppTracks_[0][2] = iBooker.book1D(
691  histname + "Endcap", " Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ", 120, -0.5, 1.0);
692  h_distMinAppTracks_[1][0] = iBooker.book1D(
693  histname + "All_Ass", " Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ", 120, -0.5, 1.0);
694  h_distMinAppTracks_[1][1] =
695  iBooker.book1D(histname + "Barrel_Ass",
696  " Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ",
697  120,
698  -0.5,
699  1.0);
700  h_distMinAppTracks_[1][2] =
701  iBooker.book1D(histname + "Endcap_Ass",
702  " Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ",
703  120,
704  -0.5,
705  1.0);
706  h_distMinAppTracks_[2][0] = iBooker.book1D(
707  histname + "All_Fakes", " Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ", 120, -0.5, 1.0);
708  h_distMinAppTracks_[2][1] =
709  iBooker.book1D(histname + "Barrel_Fakes",
710  " Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ",
711  120,
712  -0.5,
713  1.0);
714  h_distMinAppTracks_[2][2] =
715  iBooker.book1D(histname + "Endcap_Fakes",
716  " Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ",
717  120,
718  -0.5,
719  1.0);
720 
721  h_convVtxRvsZ_[0] = iBooker.book2D("convVtxRvsZAll",
722  " Photon Reco conversion vtx position",
723  zBinForXray,
724  zMinForXray,
725  zMaxForXray,
726  rBinForXray,
727  rMinForXray,
728  rMaxForXray);
729  h_convVtxRvsZ_[1] = iBooker.book2D("convVtxRvsZBarrel",
730  " Photon Reco conversion vtx position",
731  zBinForXray,
732  zMinForXray,
733  zMaxForXray,
734  rBinForXray,
735  rMinForXray,
736  rMaxForXray);
737  h_convVtxRvsZ_[2] = iBooker.book2D("convVtxRvsZEndcap",
738  " Photon Reco conversion vtx position",
739  zBin2ForXray,
740  zMinForXray,
741  zMaxForXray,
742  rBinForXray,
743  rMinForXray,
744  rMaxForXray);
745  h_convVtxYvsX_ = iBooker.book2D(
746  "convVtxYvsXTrkBarrel", " Photon Reco conversion vtx position, (x,y) eta<1 ", 1000, -60., 60., 1000, -60., 60.);
748  h_convVtxRvsZ_zoom_[0] = iBooker.book2D("convVtxRvsZBarrelZoom1",
749  " Photon Reco conversion vtx position",
750  zBinForXray,
751  zMinForXray,
752  zMaxForXray,
753  rBinForXray,
754  -10.,
755  40.);
756  h_convVtxRvsZ_zoom_[1] = iBooker.book2D("convVtxRvsZBarrelZoom2",
757  " Photon Reco conversion vtx position",
758  zBinForXray,
759  zMinForXray,
760  zMaxForXray,
761  rBinForXray,
762  -10.,
763  20.);
764  h_convVtxYvsX_zoom_[0] = iBooker.book2D("convVtxYvsXTrkBarrelZoom1",
765  " Photon Reco conversion vtx position, (x,y) eta<1 ",
766  100,
767  -40.,
768  40.,
769  100,
770  -40.,
771  40.);
772  h_convVtxYvsX_zoom_[1] = iBooker.book2D("convVtxYvsXTrkBarrelZoom2",
773  " Photon Reco conversion vtx position, (x,y) eta<1 ",
774  100,
775  -20.,
776  20.,
777  100,
778  -20.,
779  20.);
780 
781  h_convVtxdR_ = iBooker.book1D("convVtxdR", " Photon Reco conversion vtx dR", 100, -10., 10.);
782  h_convVtxdX_ = iBooker.book1D("convVtxdX", " Photon Reco conversion vtx dX", 100, -10., 10.);
783  h_convVtxdY_ = iBooker.book1D("convVtxdY", " Photon Reco conversion vtx dY", 100, -10., 10.);
784  h_convVtxdZ_ = iBooker.book1D("convVtxdZ", " Photon Reco conversion vtx dZ", 100, -20., 20.);
785 
786  h_convVtxdPhi_ = iBooker.book1D("convVtxdPhi", " Photon Reco conversion vtx dPhi", 100, -0.01, 0.01);
787  h_convVtxdEta_ = iBooker.book1D("convVtxdEta", " Photon Reco conversion vtx dEta", 100, -0.5, 0.5);
788 
789  h_convVtxdR_barrel_ =
790  iBooker.book1D("convVtxdR_barrel", " Photon Reco conversion vtx dR, |eta|<=1.2", 100, -10., 10.);
791  h_convVtxdX_barrel_ =
792  iBooker.book1D("convVtxdX_barrel", " Photon Reco conversion vtx dX, |eta|<=1.2", 100, -10., 10.);
793  h_convVtxdY_barrel_ =
794  iBooker.book1D("convVtxdY_barrel", " Photon Reco conversion vtx dY, |eta|<=1.2 ", 100, -10., 10.);
795  h_convVtxdZ_barrel_ =
796  iBooker.book1D("convVtxdZ_barrel", " Photon Reco conversion vtx dZ, |eta|<=1.2,", 100, -20., 20.);
797 
798  h_convVtxdR_endcap_ =
799  iBooker.book1D("convVtxdR_endcap", " Photon Reco conversion vtx dR, |eta|>1.2 ", 100, -10., 10.);
800  h_convVtxdX_endcap_ =
801  iBooker.book1D("convVtxdX_endcap", " Photon Reco conversion vtx dX, |eta|>1.2", 100, -10., 10.);
802  h_convVtxdY_endcap_ =
803  iBooker.book1D("convVtxdY_endcap", " Photon Reco conversion vtx dY, |eta|>1.2", 100, -10., 10.);
804  h_convVtxdZ_endcap_ =
805  iBooker.book1D("convVtxdZ_endcap", " Photon Reco conversion vtx dZ, |eta|>1.2", 100, -20., 20.);
806 
807  h2_convVtxdRVsR_ = iBooker.book2D("h2ConvVtxdRVsR", " Conversion vtx dR vsR", rBin, rMin, rMax, 100, -20., 20.);
808  h2_convVtxdRVsEta_ =
809  iBooker.book2D("h2ConvVtxdRVsEta", "Conversion vtx dR vs Eta", etaBin2, etaMin, etaMax, 100, -20., 20.);
810 
811  p_convVtxdRVsR_ =
812  iBooker.bookProfile("pConvVtxdRVsR", " Conversion vtx dR vsR", rBin, rMin, rMax, 100, -20., 20., "");
813  p_convVtxdRVsEta_ =
814  iBooker.bookProfile("pConvVtxdRVsEta", "Conversion vtx dR vs Eta", etaBin2, etaMin, etaMax, 100, -20., 20., "");
815  p_convVtxdXVsX_ = iBooker.bookProfile("pConvVtxdXVsX", "Conversion vtx dX vs X", 120, -60, 60, 100, -20., 20., "");
816  p_convVtxdYVsY_ = iBooker.bookProfile("pConvVtxdYVsY", "Conversion vtx dY vs Y", 120, -60, 60, 100, -20., 20., "");
817  p_convVtxdZVsZ_ =
818  iBooker.bookProfile("pConvVtxdZVsZ", "Conversion vtx dZ vs Z", zBin, zMin, zMax, 100, -20., 20., "");
819 
820  p_convVtxdZVsR_ =
821  iBooker.bookProfile("pConvVtxdZVsR", "Conversion vtx dZ vs R", rBin, rMin, rMax, 100, -20., 20., "");
822  p2_convVtxdRVsRZ_ = iBooker.bookProfile2D(
823  "p2ConvVtxdRVsRZ", "Conversion vtx dR vs RZ", zBin, zMin, zMax, rBin, rMin, rMax, 100, 0., 20., "s");
824  p2_convVtxdZVsRZ_ = iBooker.bookProfile2D(
825  "p2ConvVtxdZVsRZ", "Conversion vtx dZ vs RZ", zBin, zMin, zMax, rBin, rMin, rMax, 100, 0., 20., "s");
826 
827  histname = "EoverPtracks";
828  h_EoverPTracks_[0][0] =
829  iBooker.book1D(histname + "All", " photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax);
830  h_EoverPTracks_[0][1] =
831  iBooker.book1D(histname + "Barrel", " photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
832  h_EoverPTracks_[0][2] =
833  iBooker.book1D(histname + "Endcap", " photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);
834  h_EoverPTracks_[1][0] =
835  iBooker.book1D(histname + "All_Ass", " photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax);
836  h_EoverPTracks_[1][1] = iBooker.book1D(
837  histname + "Barrel_Ass", " photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
838  h_EoverPTracks_[1][2] = iBooker.book1D(
839  histname + "Endcap_Ass", " photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);
840  h_EoverPTracks_[2][0] =
841  iBooker.book1D(histname + "All_Fakes", " photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax);
842  h_EoverPTracks_[2][1] = iBooker.book1D(
843  histname + "Barrel_Fakes", " photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
844  h_EoverPTracks_[2][2] = iBooker.book1D(
845  histname + "Endcap_Fakes", " photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);
846 
847  h2_convVtxRrecVsTrue_ = iBooker.book2D(
848  "h2ConvVtxRrecVsTrue", "Photon Reco conversion vtx R rec vs true", rBin, rMin, rMax, rBin, rMin, rMax);
849 
850  histname = "vtxChi2Prob";
851  h_vtxChi2Prob_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 100, 0., 1.);
852  h_vtxChi2Prob_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 100, 0., 1.);
853  h_vtxChi2Prob_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 100, 0., 1.);
854  h_vtxChi2Prob_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 100, 0., 1.);
855  h_vtxChi2Prob_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 100, 0., 1.);
856  h_vtxChi2Prob_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 100, 0., 1.);
857  h_vtxChi2Prob_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 100, 0., 1.);
858  h_vtxChi2Prob_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 100, 0., 1.);
859  h_vtxChi2Prob_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 100, 0., 1.);
860 
861  h_zPVFromTracks_[1] = iBooker.book1D("zPVFromTracks", " Photons: PV z from conversion tracks", 100, -25., 25.);
862  h_dzPVFromTracks_[1] =
863  iBooker.book1D("dzPVFromTracks", " Photons: PV Z_rec - Z_true from conversion tracks", 100, -5., 5.);
864  h2_dzPVVsR_ = iBooker.book2D("h2dzPVVsR", "Photon Reco conversions: dz(PV) vs R", rBin, rMin, rMax, 100, -3., 3.);
865  p_dzPVVsR_ =
866  iBooker.bookProfile("pdzPVVsR", "Photon Reco conversions: dz(PV) vs R", rBin, rMin, rMax, 100, -3., 3., "");
867 
868  histname = "lxybs";
869  h_lxybs_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 200, -100., 100.);
870  h_lxybs_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 200, -100., 100.);
871  h_lxybs_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 200, -100., 100.);
872  h_lxybs_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 200, -100., 100.);
873  h_lxybs_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 200, -100., 100.);
874  h_lxybs_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 200, -100., 100.);
875  h_lxybs_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 200, -100., 100.);
876  h_lxybs_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 200, -100., 100.);
877  h_lxybs_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 200, -100., 100.);
878 
879  histname = "maxNHitsBeforeVtx";
880  h_maxNHitsBeforeVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
881  h_maxNHitsBeforeVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
882  h_maxNHitsBeforeVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
883  h_maxNHitsBeforeVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
884  h_maxNHitsBeforeVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
885  h_maxNHitsBeforeVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
886  h_maxNHitsBeforeVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
887  h_maxNHitsBeforeVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
888  h_maxNHitsBeforeVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
889 
890  histname = "leadNHitsBeforeVtx";
891  h_leadNHitsBeforeVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
892  h_leadNHitsBeforeVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
893  h_leadNHitsBeforeVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
894  h_leadNHitsBeforeVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
895  h_leadNHitsBeforeVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
896  h_leadNHitsBeforeVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
897  h_leadNHitsBeforeVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
898  h_leadNHitsBeforeVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
899  h_leadNHitsBeforeVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
900 
901  histname = "trailNHitsBeforeVtx";
902  h_trailNHitsBeforeVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
903  h_trailNHitsBeforeVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
904  h_trailNHitsBeforeVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
905  h_trailNHitsBeforeVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
906  h_trailNHitsBeforeVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
907  h_trailNHitsBeforeVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
908  h_trailNHitsBeforeVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
909  h_trailNHitsBeforeVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
910  h_trailNHitsBeforeVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
911 
912  histname = "sumNHitsBeforeVtx";
913  h_sumNHitsBeforeVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
914  h_sumNHitsBeforeVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
915  h_sumNHitsBeforeVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
916  h_sumNHitsBeforeVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
917  h_sumNHitsBeforeVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
918  h_sumNHitsBeforeVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
919  h_sumNHitsBeforeVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
920  h_sumNHitsBeforeVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
921  h_sumNHitsBeforeVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
922 
923  histname = "maxDlClosestHitToVtx";
924  h_maxDlClosestHitToVtx_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 100, -10., 10.);
925  h_maxDlClosestHitToVtx_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 100, -10., 10.);
926  h_maxDlClosestHitToVtx_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 100, -10., 10.);
927  h_maxDlClosestHitToVtx_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 100, -10., 10.);
928  h_maxDlClosestHitToVtx_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 100, -10., 10.);
929  h_maxDlClosestHitToVtx_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 100, -10., 10.);
930  h_maxDlClosestHitToVtx_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 100, -10., 10.);
931  h_maxDlClosestHitToVtx_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 100, -10., 10.);
932  h_maxDlClosestHitToVtx_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 100, -10., 10.);
933 
934  histname = "maxDlClosestHitToVtxSig";
935  h_maxDlClosestHitToVtxSig_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 100, -8., 8.);
936  h_maxDlClosestHitToVtxSig_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 100, -8., 8.);
937  h_maxDlClosestHitToVtxSig_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 100, -8., 8.);
938  h_maxDlClosestHitToVtxSig_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 100, -8., 8.);
939  h_maxDlClosestHitToVtxSig_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 100, -8., 8.);
940  h_maxDlClosestHitToVtxSig_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 100, -8., 8.);
941  h_maxDlClosestHitToVtxSig_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 100, -8., 8.);
942  h_maxDlClosestHitToVtxSig_[2][1] =
943  iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 100, -8., 8.);
944  h_maxDlClosestHitToVtxSig_[2][2] =
945  iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 100, -8., 8.);
946 
947  histname = "deltaExpectedHitsInner";
948  h_deltaExpectedHitsInner_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 31, -15.5, 15.5);
949  h_deltaExpectedHitsInner_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 31, -15.5, 15.5);
950  h_deltaExpectedHitsInner_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 31, -15.5, 15.5);
951  h_deltaExpectedHitsInner_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 31, -15.5, 15.5);
952  h_deltaExpectedHitsInner_[1][1] =
953  iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 31, -15.5, 15.5);
954  h_deltaExpectedHitsInner_[1][2] =
955  iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 31, -15.5, 15.5);
956  h_deltaExpectedHitsInner_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 31, -15.5, 15.5);
957  h_deltaExpectedHitsInner_[2][1] =
958  iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 31, -15.5, 15.5);
959  h_deltaExpectedHitsInner_[2][2] =
960  iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 31, -15.5, 15.5);
961 
962  histname = "leadExpectedHitsInner";
963  h_leadExpectedHitsInner_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
964  h_leadExpectedHitsInner_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
965  h_leadExpectedHitsInner_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
966  h_leadExpectedHitsInner_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
967  h_leadExpectedHitsInner_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
968  h_leadExpectedHitsInner_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
969  h_leadExpectedHitsInner_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
970  h_leadExpectedHitsInner_[2][1] =
971  iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
972  h_leadExpectedHitsInner_[2][2] =
973  iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
974 
975  histname = "nSharedHits";
976  h_nSharedHits_[0][0] = iBooker.book1D(histname + "All", "vertex #chi^{2} all", 16, -0.5, 15.5);
977  h_nSharedHits_[0][1] = iBooker.book1D(histname + "Barrel", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
978  h_nSharedHits_[0][2] = iBooker.book1D(histname + "Endcap", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
979  h_nSharedHits_[1][0] = iBooker.book1D(histname + "All_Ass", "vertex #chi^{2} all", 16, -0.5, 15.5);
980  h_nSharedHits_[1][1] = iBooker.book1D(histname + "Barrel_Ass", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
981  h_nSharedHits_[1][2] = iBooker.book1D(histname + "Endcap_Ass", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
982  h_nSharedHits_[2][0] = iBooker.book1D(histname + "All_Fakes", "vertex #chi^{2} all", 16, -0.5, 15.5);
983  h_nSharedHits_[2][1] = iBooker.book1D(histname + "Barrel_Fakes", "vertex #chi^{2} barrel", 16, -0.5, 15.5);
984  h_nSharedHits_[2][2] = iBooker.book1D(histname + "Endcap_Fakes", "vertex #chi^{2} endcap", 16, -0.5, 15.5);
985 
987  histname = "nHits";
988  nHits_[0] = iBooker.book2D(histname + "AllTracks",
989  "Photons:Tracks from conversions: # of hits all tracks",
990  etaBin,
991  etaMin,
992  etaMax,
993  30,
994  0.,
995  30.);
996  nHits_[1] = iBooker.book2D(histname + "AllTracks_Ass",
997  "Photons:Tracks from conversions: # of hits all tracks ass",
998  etaBin,
999  etaMin,
1000  etaMax,
1001  30,
1002  0.,
1003  30.);
1004  nHits_[2] = iBooker.book2D(histname + "AllTracks_Fakes",
1005  "Photons:Tracks from conversions: # of hits all tracks fakes",
1006  etaBin,
1007  etaMin,
1008  etaMax,
1009  30,
1010  0.,
1011  30.);
1012 
1013  histname = "nHitsVsEta";
1014  nHitsVsEta_[0] = iBooker.book2D(histname + "AllTracks",
1015  "Photons:Tracks from conversions: # of hits vs #eta all tracks",
1016  etaBin,
1017  etaMin,
1018  etaMax,
1019  30,
1020  0.,
1021  30.);
1022  nHitsVsEta_[1] = iBooker.book2D(histname + "AllTracks_Ass",
1023  "Photons:Tracks from conversions: # of hits vs #eta all tracks",
1024  etaBin,
1025  etaMin,
1026  etaMax,
1027  30,
1028  0.,
1029  30.);
1030  nHitsVsEta_[2] = iBooker.book2D(histname + "AllTracks_Fakes",
1031  "Photons:Tracks from conversions: # of hits vs #eta all tracks",
1032  etaBin,
1033  etaMin,
1034  etaMax,
1035  30,
1036  0.,
1037  30.);
1038  histname = "h_nHitsVsEta";
1039  p_nHitsVsEta_[0] = iBooker.bookProfile(histname + "AllTracks",
1040  "Photons:Tracks from conversions: # of hits vs #eta all tracks",
1041  etaBin,
1042  etaMin,
1043  etaMax,
1044  30,
1045  -0.5,
1046  29.5,
1047  "");
1048  p_nHitsVsEta_[1] = iBooker.bookProfile(histname + "AllTracks_Ass",
1049  "Photons:Tracks from conversions: # of hits vs #eta all tracks",
1050  etaBin,
1051  etaMin,
1052  etaMax,
1053  30,
1054  -0.5,
1055  29.5,
1056  "");
1057  p_nHitsVsEta_[2] = iBooker.bookProfile(histname + "AllTracks_Fakes",
1058  "Photons:Tracks from conversions: # of hits vs #eta all tracks",
1059  etaBin,
1060  etaMin,
1061  etaMax,
1062  30,
1063  -0.5,
1064  29.5,
1065  "");
1066 
1067  histname = "nHitsVsR";
1068  nHitsVsR_[0] = iBooker.book2D(histname + "AllTracks",
1069  "Photons:Tracks from conversions: # of hits vs radius all tracks",
1070  rBin,
1071  rMin,
1072  rMax,
1073  30,
1074  0.,
1075  30.);
1076  nHitsVsR_[1] = iBooker.book2D(histname + "AllTracks_Ass",
1077  "Photons:Tracks from conversions: # of hits vs radius all tracks",
1078  rBin,
1079  rMin,
1080  rMax,
1081  30,
1082  0.,
1083  30.);
1084  nHitsVsR_[2] = iBooker.book2D(histname + "AllTracks_Fakes",
1085  "Photons:Tracks from conversions: # of hits vs radius all tracks",
1086  rBin,
1087  rMin,
1088  rMax,
1089  30,
1090  0.,
1091  30.);
1092 
1093  histname = "h_nHitsVsR";
1094  p_nHitsVsR_[0] = iBooker.bookProfile(histname + "AllTracks",
1095  "Photons:Tracks from conversions: # of hits vs radius all tracks",
1096  rBin,
1097  rMin,
1098  rMax,
1099  30,
1100  -0.5,
1101  29.5,
1102  "");
1103  p_nHitsVsR_[1] = iBooker.bookProfile(histname + "AllTracks_Ass",
1104  "Photons:Tracks from conversions: # of hits vs radius all tracks",
1105  rBin,
1106  rMin,
1107  rMax,
1108  30,
1109  -0.5,
1110  29.5,
1111  "");
1112  p_nHitsVsR_[2] = iBooker.bookProfile(histname + "AllTracks_Fakes",
1113  "Photons:Tracks from conversions: # of hits vs radius all tracks",
1114  rBin,
1115  rMin,
1116  rMax,
1117  30,
1118  -0.5,
1119  29.5,
1120  "");
1121 
1122  histname = "tkChi2";
1123  h_tkChi2_[0] = iBooker.book1D(
1124  histname + "AllTracks", "Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);
1125  h_tkChi2_[1] = iBooker.book1D(
1126  histname + "AllTracks_Ass", "Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);
1127  h_tkChi2_[2] = iBooker.book1D(
1128  histname + "AllTracks_Fakes", "Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);
1129 
1130  histname = "tkChi2Large";
1131  h_tkChi2Large_[0] = iBooker.book1D(
1132  histname + "AllTracks", "Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);
1133  h_tkChi2Large_[1] = iBooker.book1D(
1134  histname + "AllTracks_Ass", "Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);
1135  h_tkChi2Large_[2] = iBooker.book1D(
1136  histname + "AllTracks_Fakes", "Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);
1137 
1138  histname = "h2Chi2VsEta";
1139  h2_Chi2VsEta_[0] = iBooker.book2D(
1140  histname + "All", " Reco Track #chi^{2} vs #eta: All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max);
1141  h2_Chi2VsEta_[1] = iBooker.book2D(
1142  histname + "All_Ass", " Reco Track #chi^{2} vs #eta: All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max);
1143  h2_Chi2VsEta_[2] = iBooker.book2D(
1144  histname + "All_Fakes", " Reco Track #chi^{2} vs #eta: All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max);
1145  histname = "pChi2VsEta";
1146  p_Chi2VsEta_[0] = iBooker.bookProfile(
1147  histname + "All", " Reco Track #chi^{2} vs #eta : All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max, "");
1148  p_Chi2VsEta_[1] = iBooker.bookProfile(
1149  histname + "All_Ass", " Reco Track #chi^{2} vs #eta : All ", etaBin2, etaMin, etaMax, 100, chi2Min, chi2Max, "");
1150  p_Chi2VsEta_[2] = iBooker.bookProfile(histname + "All_Fakes",
1151  " Reco Track #chi^{2} vs #eta : All ",
1152  etaBin2,
1153  etaMin,
1154  etaMax,
1155  100,
1156  chi2Min,
1157  chi2Max,
1158  "");
1159 
1160  histname = "h2Chi2VsR";
1161  h2_Chi2VsR_[0] =
1162  iBooker.book2D(histname + "All", " Reco Track #chi^{2} vs R: All ", rBin, rMin, rMax, 100, chi2Min, chi2Max);
1163  h2_Chi2VsR_[1] = iBooker.book2D(
1164  histname + "All_Ass", " Reco Track #chi^{2} vs R: All ", rBin, rMin, rMax, 100, chi2Min, chi2Max);
1165  h2_Chi2VsR_[2] = iBooker.book2D(
1166  histname + "All_Fakes", " Reco Track #chi^{2} vs R: All ", rBin, rMin, rMax, 100, chi2Min, chi2Max);
1167  histname = "pChi2VsR";
1168  p_Chi2VsR_[0] = iBooker.bookProfile(
1169  histname + "All", " Reco Track #chi^{2} vas R : All ", rBin, rMin, rMax, 100, chi2Min, chi2Max, "");
1170  p_Chi2VsR_[1] = iBooker.bookProfile(
1171  histname + "All_Ass", " Reco Track #chi^{2} vas R : All ", rBin, rMin, rMax, 100, chi2Min, chi2Max, "");
1172  p_Chi2VsR_[2] = iBooker.bookProfile(
1173  histname + "All_Fakes", " Reco Track #chi^{2} vas R : All ", rBin, rMin, rMax, 100, chi2Min, chi2Max, "");
1174 
1175  histname = "hTkD0";
1176  h_TkD0_[0] = iBooker.book1D(histname + "All", " Reco Track D0*q: All ", 200, -0.1, 60);
1177  h_TkD0_[1] = iBooker.book1D(histname + "All_Ass", " Reco Track D0*q: Barrel ", 200, -0.1, 60);
1178  h_TkD0_[2] = iBooker.book1D(histname + "All_Fakes", " Reco Track D0*q: Endcap ", 200, -0.1, 60);
1179 
1180  histname = "hTkPtPull";
1181  h_TkPtPull_[0] = iBooker.book1D(histname + "All", " Reco Track Pt pull: All ", 100, -20., 10.);
1182  histname = "hTkPtPull";
1183  h_TkPtPull_[1] = iBooker.book1D(histname + "Barrel", " Reco Track Pt pull: Barrel ", 100, -20., 10.);
1184  histname = "hTkPtPull";
1185  h_TkPtPull_[2] = iBooker.book1D(histname + "Endcap", " Reco Track Pt pull: Endcap ", 100, -20., 10.);
1186 
1187  histname = "h2TkPtPullEta";
1188  h2_TkPtPull_[0] =
1189  iBooker.book2D(histname + "All", " Reco Track Pt pull: All ", etaBin2, etaMin, etaMax, 100, -20., 10.);
1190  histname = "pTkPtPullEta";
1191  p_TkPtPull_[0] = iBooker.bookProfile(
1192  histname + "All", " Reco Track Pt pull: All ", etaBin2, etaMin, etaMax, 100, -20., 10., " ");
1193 
1194  histname = "PtRecVsPtSim";
1195  h2_PtRecVsPtSim_[0] =
1196  iBooker.book2D(histname + "All", "Pt Rec vs Pt sim: All ", etBin, etMin, etMax, etBin, etMin, etMax);
1197  h2_PtRecVsPtSim_[1] =
1198  iBooker.book2D(histname + "Barrel", "Pt Rec vs Pt sim: Barrel ", etBin, etMin, etMax, etBin, etMin, etMax);
1199  h2_PtRecVsPtSim_[2] =
1200  iBooker.book2D(histname + "Endcap", "Pt Rec vs Pt sim: Endcap ", etBin, etMin, etMax, etBin, etMin, etMax);
1201 
1202  histname = "photonPtRecVsPtSim";
1203  h2_photonPtRecVsPtSim_ =
1204  iBooker.book2D(histname + "All", "Pt Rec vs Pt sim: All ", etBin, etMin, etMax, etBin, etMin, etMax);
1205 
1206  histname = "nHitsBeforeVtx";
1207  h_nHitsBeforeVtx_[0] = iBooker.book1D(histname + "All", "Pt Rec vs Pt sim: All ", 16, -0.5, 15.5);
1208  h_nHitsBeforeVtx_[1] = iBooker.book1D(histname + "Barrel", "Pt Rec vs Pt sim: Barrel ", 16, -0.5, 15.5);
1209  h_nHitsBeforeVtx_[2] = iBooker.book1D(histname + "Endcap", "Pt Rec vs Pt sim: Endcap ", 16, -0.5, 15.5);
1210 
1211  histname = "dlClosestHitToVtx";
1212  h_dlClosestHitToVtx_[0] = iBooker.book1D(histname + "All", "Pt Rec vs Pt sim: All ", 100, -10., 10.);
1213  h_dlClosestHitToVtx_[1] = iBooker.book1D(histname + "Barrel", "Pt Rec vs Pt sim: Barrel ", 100, -10., 10.);
1214  h_dlClosestHitToVtx_[2] = iBooker.book1D(histname + "Endcap", "Pt Rec vs Pt sim: Endcap ", 100, -10., 10.);
1215 
1216  histname = "dlClosestHitToVtxSig";
1217  h_dlClosestHitToVtxSig_[0] = iBooker.book1D(histname + "All", "Pt Rec vs Pt sim: All ", 100, -8., 8.);
1218  h_dlClosestHitToVtxSig_[1] = iBooker.book1D(histname + "Barrel", "Pt Rec vs Pt sim: Barrel ", 100, -8., 8.);
1219  h_dlClosestHitToVtxSig_[2] = iBooker.book1D(histname + "Endcap", "Pt Rec vs Pt sim: Endcap ", 100, -8., 8.);
1220 
1221  h_match_ = iBooker.book1D("h_match", " ", 3, -0.5, 2.5);
1222 
1223  } // if DQM
1224 }

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

◆ dqmBeginRun()

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

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 1226 of file TkConvValidator.cc.

1226  {
1227  //get magnetic field
1228  edm::LogInfo("ConvertedPhotonProducer") << " get magnetic field"
1229  << "\n";
1230  theEventSetup.get<IdealMagneticFieldRecord>().get(theMF_);
1231 
1232  thePhotonMCTruthFinder_ = new PhotonMCTruthFinder();
1233 }

References edm::EventSetup::get(), and get.

◆ dqmEndRun()

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

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 1235 of file TkConvValidator.cc.

1235  {
1236  delete thePhotonMCTruthFinder_;
1237 }

◆ endJob()

void TkConvValidator::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 2231 of file TkConvValidator.cc.

2231  {
2232  std::string outputFileName = parameters_.getParameter<std::string>("OutputFileName");
2233  if (!isRunCentrally_) {
2234  dbe_->save(outputFileName);
2235  }
2236 
2237  edm::LogInfo("TkConvValidator") << "Analyzed " << nEvt_ << "\n";
2238  // std::cout << "::endJob Analyzed " << nEvt_ << " events " << " with total " << nPho_ << " Photons " << "\n";
2239  // std::cout << "TkConvValidator::endJob Analyzed " << nEvt_ << " events " << "\n";
2240 
2241  return;
2242 }

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

◆ etaTransformation()

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

Definition at line 2294 of file TkConvValidator.cc.

2294  {
2295  //---Definitions
2296  const float PI = 3.1415927;
2297 
2298  //---Definitions for ECAL
2299  const float R_ECAL = 136.5;
2300  const float Z_Endcap = 328.0;
2301  const float etaBarrelEndcap = 1.479;
2302 
2303  //---ETA correction
2304 
2305  float Theta = 0.0;
2306  float ZEcal = R_ECAL * sinh(EtaParticle) + Zvertex;
2307 
2308  if (ZEcal != 0.0)
2309  Theta = atan(R_ECAL / ZEcal);
2310  if (Theta < 0.0)
2311  Theta = Theta + PI;
2312  float ETA = -log(tan(0.5 * Theta));
2313 
2314  if (fabs(ETA) > etaBarrelEndcap) {
2315  float Zend = Z_Endcap;
2316  if (EtaParticle < 0.0)
2317  Zend = -Zend;
2318  float Zlen = Zend - Zvertex;
2319  float RR = Zlen / sinh(EtaParticle);
2320  Theta = atan(RR / Zend);
2321  if (Theta < 0.0)
2322  Theta = Theta + PI;
2323  ETA = -log(tan(0.5 * Theta));
2324  }
2325  //---Return the result
2326  return ETA;
2327  //---end
2328 }

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

◆ phiNormalization()

float TkConvValidator::phiNormalization ( float &  a)
private

Definition at line 2278 of file TkConvValidator.cc.

2278  {
2279  //---Definitions
2280  const float PI = 3.1415927;
2281  const float TWOPI = 2.0 * PI;
2282 
2283  if (phi > PI) {
2284  phi = phi - TWOPI;
2285  }
2286  if (phi < -PI) {
2287  phi = phi + TWOPI;
2288  }
2289 
2290  // cout << " Float_t PHInormalization out " << PHI << endl;
2291  return phi;
2292 }

References PI, and TWOPI.

◆ 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 2244 of file TkConvValidator.cc.

2247  {
2248  math::XYZVector result;
2249  Surface::RotationType rot;
2250  auto scp = new SimpleCylinderBounds(sqrt(vtx.position().perp2()) - 0.001f,
2251  sqrt(vtx.position().perp2()) + 0.001f,
2252  -fabs(vtx.position().z()),
2253  fabs(vtx.position().z()));
2254  ReferenceCountingPointer<Cylinder> theBarrel_(
2255  new Cylinder(Cylinder::computeRadius(*scp), Surface::PositionType(0, 0, 0), rot, scp));
2256 
2257  ReferenceCountingPointer<Disk> theDisk_(
2258  new Disk(Surface::PositionType(0, 0, vtx.position().z()),
2259  rot,
2260  new SimpleDiskBounds(0, sqrt(vtx.position().perp2()), -0.001, 0.001)));
2261 
2262  const TrajectoryStateOnSurface myTSOS = trajectoryStateTransform::innerStateOnSurface(*tk, trackerGeom, &mf);
2263  PropagatorWithMaterial propag(anyDirection, 0.000511, &mf);
2264  TrajectoryStateOnSurface stateAtVtx;
2265  stateAtVtx = propag.propagate(myTSOS, *theBarrel_);
2266  if (!stateAtVtx.isValid()) {
2267  stateAtVtx = propag.propagate(myTSOS, *theDisk_);
2268  }
2269  if (stateAtVtx.isValid()) {
2270  return math::XYZVector(double(stateAtVtx.globalMomentum().x()),
2271  double(stateAtVtx.globalMomentum().y()),
2272  double(stateAtVtx.globalMomentum().z()));
2273  } else {
2274  return math::XYZVector(0., 0., 0.);
2275  }
2276 }

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

Member Data Documentation

◆ arbitratedEcalSeeded_

bool TkConvValidator::arbitratedEcalSeeded_
private

Definition at line 137 of file TkConvValidator.h.

◆ arbitratedMerged_

bool TkConvValidator::arbitratedMerged_
private

Definition at line 136 of file TkConvValidator.h.

◆ bcEtLow_

double TkConvValidator::bcEtLow_
private

Definition at line 124 of file TkConvValidator.h.

◆ beamspotToken_

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

Definition at line 101 of file TkConvValidator.h.

◆ conversionCollection_

std::string TkConvValidator::conversionCollection_
private

Definition at line 91 of file TkConvValidator.h.

◆ conversionCollectionPr_Token_

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

Definition at line 92 of file TkConvValidator.h.

◆ conversionCollectionProducer_

std::string TkConvValidator::conversionCollectionProducer_
private

Definition at line 90 of file TkConvValidator.h.

◆ conversionTrackProducer_

std::string TkConvValidator::conversionTrackProducer_
private

Definition at line 94 of file TkConvValidator.h.

◆ dbe_

DQMStore* TkConvValidator::dbe_
private

Definition at line 72 of file TkConvValidator.h.

◆ dCotCutOn_

bool TkConvValidator::dCotCutOn_
private

Definition at line 132 of file TkConvValidator.h.

◆ dCotCutValue_

double TkConvValidator::dCotCutValue_
private

Definition at line 133 of file TkConvValidator.h.

◆ dCotHardCutValue_

double TkConvValidator::dCotHardCutValue_
private

Definition at line 134 of file TkConvValidator.h.

◆ dqmpath_

std::string TkConvValidator::dqmpath_
private

Definition at line 110 of file TkConvValidator.h.

◆ ecalalgotracks_

bool TkConvValidator::ecalalgotracks_
private

Definition at line 138 of file TkConvValidator.h.

◆ ecalEtSumCut_

double TkConvValidator::ecalEtSumCut_
private

Definition at line 130 of file TkConvValidator.h.

◆ ecalIsolRadius_

double TkConvValidator::ecalIsolRadius_
private

Definition at line 123 of file TkConvValidator.h.

◆ fName_

std::string TkConvValidator::fName_
private

Definition at line 71 of file TkConvValidator.h.

◆ g4_simTk_Token_

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

Definition at line 102 of file TkConvValidator.h.

◆ g4_simVtx_Token_

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

Definition at line 103 of file TkConvValidator.h.

◆ generalTracksOnly_

bool TkConvValidator::generalTracksOnly_
private

Definition at line 135 of file TkConvValidator.h.

◆ genjets_Token_

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

Definition at line 107 of file TkConvValidator.h.

◆ h2_Chi2VsEta_

MonitorElement* TkConvValidator::h2_Chi2VsEta_[3]
private

Definition at line 317 of file TkConvValidator.h.

◆ h2_Chi2VsR_

MonitorElement* TkConvValidator::h2_Chi2VsR_[3]
private

Definition at line 319 of file TkConvValidator.h.

◆ h2_convVtxdRVsEta_

MonitorElement* TkConvValidator::h2_convVtxdRVsEta_
private

Definition at line 276 of file TkConvValidator.h.

◆ h2_convVtxdRVsR_

MonitorElement* TkConvValidator::h2_convVtxdRVsR_
private

Definition at line 275 of file TkConvValidator.h.

◆ h2_convVtxRrecVsTrue_

MonitorElement* TkConvValidator::h2_convVtxRrecVsTrue_
private

Definition at line 289 of file TkConvValidator.h.

◆ h2_DCotTracksVsEta_

MonitorElement* TkConvValidator::h2_DCotTracksVsEta_
private

Definition at line 238 of file TkConvValidator.h.

◆ h2_DCotTracksVsR_

MonitorElement* TkConvValidator::h2_DCotTracksVsR_
private

Definition at line 240 of file TkConvValidator.h.

◆ h2_DPhiTracksAtEcalVsEta_

MonitorElement* TkConvValidator::h2_DPhiTracksAtEcalVsEta_
private

Definition at line 248 of file TkConvValidator.h.

◆ h2_DPhiTracksAtEcalVsR_

MonitorElement* TkConvValidator::h2_DPhiTracksAtEcalVsR_
private

Definition at line 246 of file TkConvValidator.h.

◆ h2_DPhiTracksAtVtxVsEta_

MonitorElement* TkConvValidator::h2_DPhiTracksAtVtxVsEta_
private

Definition at line 232 of file TkConvValidator.h.

◆ h2_DPhiTracksAtVtxVsR_

MonitorElement* TkConvValidator::h2_DPhiTracksAtVtxVsR_
private

Definition at line 234 of file TkConvValidator.h.

◆ h2_dzPVVsR_

MonitorElement* TkConvValidator::h2_dzPVVsR_
private

Definition at line 295 of file TkConvValidator.h.

◆ h2_photonPtRecVsPtSim_

MonitorElement* TkConvValidator::h2_photonPtRecVsPtSim_
private

Definition at line 328 of file TkConvValidator.h.

◆ h2_PtRecVsPtSim_

MonitorElement* TkConvValidator::h2_PtRecVsPtSim_[3]
private

Definition at line 327 of file TkConvValidator.h.

◆ h2_TkPtPull_

MonitorElement* TkConvValidator::h2_TkPtPull_[3]
private

Definition at line 325 of file TkConvValidator.h.

◆ h_AllSimConv_

MonitorElement* TkConvValidator::h_AllSimConv_[5]
private

Denominator for efficiencies.

Definition at line 188 of file TkConvValidator.h.

◆ h_convEta2_

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

Definition at line 214 of file TkConvValidator.h.

◆ h_convEta_

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

Definition at line 213 of file TkConvValidator.h.

◆ h_convEtaMatchSC_

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

Definition at line 212 of file TkConvValidator.h.

◆ h_convPhi_

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

Definition at line 215 of file TkConvValidator.h.

◆ h_convPt_

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

Definition at line 218 of file TkConvValidator.h.

◆ h_convPtRes_

MonitorElement* TkConvValidator::h_convPtRes_[3]
private

Definition at line 227 of file TkConvValidator.h.

◆ h_convR_

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

Definition at line 216 of file TkConvValidator.h.

◆ h_convRplot_

MonitorElement* TkConvValidator::h_convRplot_
private

Definition at line 224 of file TkConvValidator.h.

◆ h_convSCdEta_

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

Definition at line 221 of file TkConvValidator.h.

◆ h_convSCdPhi_

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

Definition at line 222 of file TkConvValidator.h.

◆ h_convVtxdEta_

MonitorElement* TkConvValidator::h_convVtxdEta_
private

Definition at line 262 of file TkConvValidator.h.

◆ h_convVtxdPhi_

MonitorElement* TkConvValidator::h_convVtxdPhi_
private

Definition at line 263 of file TkConvValidator.h.

◆ h_convVtxdR_

MonitorElement* TkConvValidator::h_convVtxdR_
private

Definition at line 261 of file TkConvValidator.h.

◆ h_convVtxdR_barrel_

MonitorElement* TkConvValidator::h_convVtxdR_barrel_
private

Definition at line 268 of file TkConvValidator.h.

◆ h_convVtxdR_endcap_

MonitorElement* TkConvValidator::h_convVtxdR_endcap_
private

Definition at line 273 of file TkConvValidator.h.

◆ h_convVtxdX_

MonitorElement* TkConvValidator::h_convVtxdX_
private

Definition at line 258 of file TkConvValidator.h.

◆ h_convVtxdX_barrel_

MonitorElement* TkConvValidator::h_convVtxdX_barrel_
private

Definition at line 265 of file TkConvValidator.h.

◆ h_convVtxdX_endcap_

MonitorElement* TkConvValidator::h_convVtxdX_endcap_
private

Definition at line 270 of file TkConvValidator.h.

◆ h_convVtxdY_

MonitorElement* TkConvValidator::h_convVtxdY_
private

Definition at line 259 of file TkConvValidator.h.

◆ h_convVtxdY_barrel_

MonitorElement* TkConvValidator::h_convVtxdY_barrel_
private

Definition at line 266 of file TkConvValidator.h.

◆ h_convVtxdY_endcap_

MonitorElement* TkConvValidator::h_convVtxdY_endcap_
private

Definition at line 271 of file TkConvValidator.h.

◆ h_convVtxdZ_

MonitorElement* TkConvValidator::h_convVtxdZ_
private

Definition at line 260 of file TkConvValidator.h.

◆ h_convVtxdZ_barrel_

MonitorElement* TkConvValidator::h_convVtxdZ_barrel_
private

Definition at line 267 of file TkConvValidator.h.

◆ h_convVtxdZ_endcap_

MonitorElement* TkConvValidator::h_convVtxdZ_endcap_
private

Definition at line 272 of file TkConvValidator.h.

◆ h_convVtxRvsZ_

MonitorElement* TkConvValidator::h_convVtxRvsZ_[3]
private

Definition at line 253 of file TkConvValidator.h.

◆ h_convVtxRvsZ_zoom_

MonitorElement* TkConvValidator::h_convVtxRvsZ_zoom_[2]
private

Definition at line 255 of file TkConvValidator.h.

◆ h_convVtxYvsX_

MonitorElement* TkConvValidator::h_convVtxYvsX_
private

Definition at line 254 of file TkConvValidator.h.

◆ h_convVtxYvsX_zoom_

MonitorElement* TkConvValidator::h_convVtxYvsX_zoom_[2]
private

Definition at line 256 of file TkConvValidator.h.

◆ h_convZ_

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

Definition at line 217 of file TkConvValidator.h.

◆ h_convZplot_

MonitorElement* TkConvValidator::h_convZplot_
private

Definition at line 225 of file TkConvValidator.h.

◆ h_DCotTracks_

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

Definition at line 237 of file TkConvValidator.h.

◆ h_deltaExpectedHitsInner_

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

Definition at line 303 of file TkConvValidator.h.

◆ h_DEtaTracksAtEcal_

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

Definition at line 251 of file TkConvValidator.h.

◆ h_distMinAppTracks_

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

Definition at line 243 of file TkConvValidator.h.

◆ h_dlClosestHitToVtx_

MonitorElement* TkConvValidator::h_dlClosestHitToVtx_[3]
private

Definition at line 335 of file TkConvValidator.h.

◆ h_dlClosestHitToVtxSig_

MonitorElement* TkConvValidator::h_dlClosestHitToVtxSig_[3]
private

Definition at line 336 of file TkConvValidator.h.

◆ h_DPhiTracksAtEcal_

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

Definition at line 245 of file TkConvValidator.h.

◆ h_DPhiTracksAtVtx_

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

Definition at line 231 of file TkConvValidator.h.

◆ h_dzPVFromTracks_

MonitorElement* TkConvValidator::h_dzPVFromTracks_[2]
private

Definition at line 294 of file TkConvValidator.h.

◆ h_EoverPTracks_

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

Definition at line 219 of file TkConvValidator.h.

◆ h_invMass_

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

Definition at line 229 of file TkConvValidator.h.

◆ h_leadExpectedHitsInner_

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

Definition at line 304 of file TkConvValidator.h.

◆ h_leadNHitsBeforeVtx_

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

Definition at line 300 of file TkConvValidator.h.

◆ h_lxybs_

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

Definition at line 298 of file TkConvValidator.h.

◆ h_match_

MonitorElement* TkConvValidator::h_match_
private

Definition at line 330 of file TkConvValidator.h.

◆ h_maxDlClosestHitToVtx_

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

Definition at line 305 of file TkConvValidator.h.

◆ h_maxDlClosestHitToVtxSig_

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

Definition at line 306 of file TkConvValidator.h.

◆ h_maxNHitsBeforeVtx_

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

Definition at line 299 of file TkConvValidator.h.

◆ h_nConv_

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

info per conversion

Definition at line 211 of file TkConvValidator.h.

◆ h_nHitsBeforeVtx_

MonitorElement* TkConvValidator::h_nHitsBeforeVtx_[3]
private

Definition at line 334 of file TkConvValidator.h.

◆ h_nSharedHits_

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

Definition at line 307 of file TkConvValidator.h.

◆ h_nSimConv_

MonitorElement* TkConvValidator::h_nSimConv_[2]
private

Definition at line 178 of file TkConvValidator.h.

◆ h_RecoConvTwoMTracks_

MonitorElement* TkConvValidator::h_RecoConvTwoMTracks_[5]
private

Definition at line 208 of file TkConvValidator.h.

◆ h_RecoConvTwoTracks_

MonitorElement* TkConvValidator::h_RecoConvTwoTracks_[5]
private

Definition at line 206 of file TkConvValidator.h.

◆ h_SimConvEtaPix_

MonitorElement* TkConvValidator::h_SimConvEtaPix_[2]
private

Definition at line 179 of file TkConvValidator.h.

◆ h_SimConvOneMTracks_

MonitorElement* TkConvValidator::h_SimConvOneMTracks_[5]
private

Definition at line 193 of file TkConvValidator.h.

◆ h_SimConvOneTracks_

MonitorElement* TkConvValidator::h_SimConvOneTracks_[5]
private

Numerator for efficiencies.

Definition at line 192 of file TkConvValidator.h.

◆ h_SimConvTwoMTracks_

MonitorElement* TkConvValidator::h_SimConvTwoMTracks_[5]
private

Definition at line 195 of file TkConvValidator.h.

◆ h_SimConvTwoMTracksAndVtxPGT0005_

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT0005_[5]
private

Definition at line 197 of file TkConvValidator.h.

◆ h_SimConvTwoMTracksAndVtxPGT01_

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT01_[5]
private

Definition at line 198 of file TkConvValidator.h.

◆ h_SimConvTwoMTracksAndVtxPGT0_

MonitorElement* TkConvValidator::h_SimConvTwoMTracksAndVtxPGT0_[5]
private

Definition at line 196 of file TkConvValidator.h.

◆ h_SimConvTwoTracks_

MonitorElement* TkConvValidator::h_SimConvTwoTracks_[5]
private

Definition at line 194 of file TkConvValidator.h.

◆ h_simConvVtxRvsZ_

MonitorElement* TkConvValidator::h_simConvVtxRvsZ_[4]
private

Definition at line 184 of file TkConvValidator.h.

◆ h_simConvVtxYvsX_

MonitorElement* TkConvValidator::h_simConvVtxYvsX_
private

Definition at line 185 of file TkConvValidator.h.

◆ h_SimRecConvOneMTracks_

MonitorElement* TkConvValidator::h_SimRecConvOneMTracks_[5]
private

Definition at line 201 of file TkConvValidator.h.

◆ h_SimRecConvOneTracks_

MonitorElement* TkConvValidator::h_SimRecConvOneTracks_[5]
private

Definition at line 200 of file TkConvValidator.h.

◆ h_SimRecConvTwoMTracks_

MonitorElement* TkConvValidator::h_SimRecConvTwoMTracks_[5]
private

Definition at line 203 of file TkConvValidator.h.

◆ h_SimRecConvTwoTracks_

MonitorElement* TkConvValidator::h_SimRecConvTwoTracks_[5]
private

Definition at line 202 of file TkConvValidator.h.

◆ h_simTkEta_

MonitorElement* TkConvValidator::h_simTkEta_
private

Definition at line 182 of file TkConvValidator.h.

◆ h_simTkPt_

MonitorElement* TkConvValidator::h_simTkPt_
private

Definition at line 181 of file TkConvValidator.h.

◆ h_sumNHitsBeforeVtx_

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

Definition at line 302 of file TkConvValidator.h.

◆ h_tkChi2_

MonitorElement* TkConvValidator::h_tkChi2_[3]
private

Definition at line 315 of file TkConvValidator.h.

◆ h_tkChi2Large_

MonitorElement* TkConvValidator::h_tkChi2Large_[3]
private

Definition at line 316 of file TkConvValidator.h.

◆ h_TkD0_

MonitorElement* TkConvValidator::h_TkD0_[3]
private

Definition at line 322 of file TkConvValidator.h.

◆ h_TkPtPull_

MonitorElement* TkConvValidator::h_TkPtPull_[3]
private

Definition at line 324 of file TkConvValidator.h.

◆ h_trailNHitsBeforeVtx_

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

Definition at line 301 of file TkConvValidator.h.

◆ h_VisSimConv_

MonitorElement* TkConvValidator::h_VisSimConv_[6]
private

Definition at line 189 of file TkConvValidator.h.

◆ h_VisSimConvLarge_

MonitorElement* TkConvValidator::h_VisSimConvLarge_
private

Definition at line 190 of file TkConvValidator.h.

◆ h_vtxChi2Prob_

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

Definition at line 291 of file TkConvValidator.h.

◆ h_zPVFromTracks_

MonitorElement* TkConvValidator::h_zPVFromTracks_[2]
private

Definition at line 293 of file TkConvValidator.h.

◆ hcalEtSumCut_

double TkConvValidator::hcalEtSumCut_
private

Definition at line 131 of file TkConvValidator.h.

◆ hcalHitEtLow_

double TkConvValidator::hcalHitEtLow_
private

Definition at line 127 of file TkConvValidator.h.

◆ hcalIsolExtRadius_

double TkConvValidator::hcalIsolExtRadius_
private

Definition at line 125 of file TkConvValidator.h.

◆ hcalIsolInnRadius_

double TkConvValidator::hcalIsolInnRadius_
private

Definition at line 126 of file TkConvValidator.h.

◆ hepMC_Token_

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

Definition at line 106 of file TkConvValidator.h.

◆ highPurity_

bool TkConvValidator::highPurity_
private

Definition at line 139 of file TkConvValidator.h.

◆ isRunCentrally_

bool TkConvValidator::isRunCentrally_
private

Definition at line 116 of file TkConvValidator.h.

◆ label_tp_

edm::InputTag TkConvValidator::label_tp_
private

Definition at line 112 of file TkConvValidator.h.

◆ lip_

double TkConvValidator::lip_
private

Definition at line 122 of file TkConvValidator.h.

◆ maxHitsBeforeVtx_

uint TkConvValidator::maxHitsBeforeVtx_
private

Definition at line 141 of file TkConvValidator.h.

◆ maxPhoEtaForEffic

double TkConvValidator::maxPhoEtaForEffic
private

Definition at line 161 of file TkConvValidator.h.

◆ maxPhoEtaForPurity

double TkConvValidator::maxPhoEtaForPurity
private

Definition at line 165 of file TkConvValidator.h.

◆ maxPhoRForEffic

double TkConvValidator::maxPhoRForEffic
private

Definition at line 163 of file TkConvValidator.h.

◆ maxPhoRForPurity

double TkConvValidator::maxPhoRForPurity
private

Definition at line 167 of file TkConvValidator.h.

◆ maxPhoZForEffic

double TkConvValidator::maxPhoZForEffic
private

Definition at line 162 of file TkConvValidator.h.

◆ maxPhoZForPurity

double TkConvValidator::maxPhoZForPurity
private

Definition at line 166 of file TkConvValidator.h.

◆ mcConvEta_

double TkConvValidator::mcConvEta_
private

Definition at line 153 of file TkConvValidator.h.

◆ mcConvPhi_

double TkConvValidator::mcConvPhi_
private

Definition at line 152 of file TkConvValidator.h.

◆ mcConvPt_

double TkConvValidator::mcConvPt_
private

Definition at line 147 of file TkConvValidator.h.

◆ mcConvR_

double TkConvValidator::mcConvR_
private

Definition at line 148 of file TkConvValidator.h.

◆ mcConvX_

double TkConvValidator::mcConvX_
private

Definition at line 151 of file TkConvValidator.h.

◆ mcConvY_

double TkConvValidator::mcConvY_
private

Definition at line 150 of file TkConvValidator.h.

◆ mcConvZ_

double TkConvValidator::mcConvZ_
private

Definition at line 149 of file TkConvValidator.h.

◆ mcEta_

double TkConvValidator::mcEta_
private

Definition at line 146 of file TkConvValidator.h.

◆ mcJetEta_

double TkConvValidator::mcJetEta_
private

Definition at line 154 of file TkConvValidator.h.

◆ mcJetPhi_

double TkConvValidator::mcJetPhi_
private

Definition at line 155 of file TkConvValidator.h.

◆ mcPhi_

double TkConvValidator::mcPhi_
private

global variable for the MC photon

Definition at line 145 of file TkConvValidator.h.

◆ minLxy_

double TkConvValidator::minLxy_
private

Definition at line 142 of file TkConvValidator.h.

◆ minPhoEtCut_

double TkConvValidator::minPhoEtCut_
private

Definition at line 118 of file TkConvValidator.h.

◆ minPhoPtForEffic

double TkConvValidator::minPhoPtForEffic
private

Definition at line 160 of file TkConvValidator.h.

◆ minPhoPtForPurity

double TkConvValidator::minPhoPtForPurity
private

Definition at line 164 of file TkConvValidator.h.

◆ minProb_

double TkConvValidator::minProb_
private

Definition at line 140 of file TkConvValidator.h.

◆ nEntry_

int TkConvValidator::nEntry_
private

Definition at line 77 of file TkConvValidator.h.

◆ nEvt_

int TkConvValidator::nEvt_
private

Definition at line 76 of file TkConvValidator.h.

◆ nHits_

MonitorElement* TkConvValidator::nHits_[3]
private

Definition at line 310 of file TkConvValidator.h.

◆ nHitsVsEta_

MonitorElement* TkConvValidator::nHitsVsEta_[3]
private

Definition at line 312 of file TkConvValidator.h.

◆ nHitsVsR_

MonitorElement* TkConvValidator::nHitsVsR_[3]
private

Definition at line 314 of file TkConvValidator.h.

◆ nInvalidPCA_

int TkConvValidator::nInvalidPCA_
private

Definition at line 84 of file TkConvValidator.h.

◆ nMatched_

int TkConvValidator::nMatched_
private

Definition at line 79 of file TkConvValidator.h.

◆ nRecConv_

int TkConvValidator::nRecConv_
private

Definition at line 80 of file TkConvValidator.h.

◆ nRecConvAss_

int TkConvValidator::nRecConvAss_
private

Definition at line 81 of file TkConvValidator.h.

◆ nRecConvAssWithEcal_

int TkConvValidator::nRecConvAssWithEcal_
private

Definition at line 82 of file TkConvValidator.h.

◆ nSimConv_

int TkConvValidator::nSimConv_[2]
private

Definition at line 78 of file TkConvValidator.h.

◆ numOfTracksInCone_

int TkConvValidator::numOfTracksInCone_
private

Definition at line 128 of file TkConvValidator.h.

◆ offline_pvToken_

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

Definition at line 100 of file TkConvValidator.h.

◆ p2_convVtxdRVsRZ_

MonitorElement* TkConvValidator::p2_convVtxdRVsRZ_
private

Definition at line 286 of file TkConvValidator.h.

◆ p2_convVtxdZVsRZ_

MonitorElement* TkConvValidator::p2_convVtxdZVsRZ_
private

Definition at line 287 of file TkConvValidator.h.

◆ p2_effRZ_

MonitorElement* TkConvValidator::p2_effRZ_
private

Definition at line 332 of file TkConvValidator.h.

◆ p_Chi2VsEta_

MonitorElement* TkConvValidator::p_Chi2VsEta_[3]
private

Definition at line 318 of file TkConvValidator.h.

◆ p_Chi2VsR_

MonitorElement* TkConvValidator::p_Chi2VsR_[3]
private

Definition at line 320 of file TkConvValidator.h.

◆ p_convVtxdRVsEta_

MonitorElement* TkConvValidator::p_convVtxdRVsEta_
private

Definition at line 279 of file TkConvValidator.h.

◆ p_convVtxdRVsR_

MonitorElement* TkConvValidator::p_convVtxdRVsR_
private

Definition at line 278 of file TkConvValidator.h.

◆ p_convVtxdXVsX_

MonitorElement* TkConvValidator::p_convVtxdXVsX_
private

Definition at line 281 of file TkConvValidator.h.

◆ p_convVtxdYVsY_

MonitorElement* TkConvValidator::p_convVtxdYVsY_
private

Definition at line 282 of file TkConvValidator.h.

◆ p_convVtxdZVsR_

MonitorElement* TkConvValidator::p_convVtxdZVsR_
private

Definition at line 284 of file TkConvValidator.h.

◆ p_convVtxdZVsZ_

MonitorElement* TkConvValidator::p_convVtxdZVsZ_
private

Definition at line 283 of file TkConvValidator.h.

◆ p_DCotTracksVsEta_

MonitorElement* TkConvValidator::p_DCotTracksVsEta_
private

Definition at line 239 of file TkConvValidator.h.

◆ p_DCotTracksVsR_

MonitorElement* TkConvValidator::p_DCotTracksVsR_
private

Definition at line 241 of file TkConvValidator.h.

◆ p_DPhiTracksAtEcalVsEta_

MonitorElement* TkConvValidator::p_DPhiTracksAtEcalVsEta_
private

Definition at line 249 of file TkConvValidator.h.

◆ p_DPhiTracksAtEcalVsR_

MonitorElement* TkConvValidator::p_DPhiTracksAtEcalVsR_
private

Definition at line 247 of file TkConvValidator.h.

◆ p_DPhiTracksAtVtxVsEta_

MonitorElement* TkConvValidator::p_DPhiTracksAtVtxVsEta_
private

Definition at line 233 of file TkConvValidator.h.

◆ p_DPhiTracksAtVtxVsR_

MonitorElement* TkConvValidator::p_DPhiTracksAtVtxVsR_
private

Definition at line 235 of file TkConvValidator.h.

◆ p_dzPVVsR_

MonitorElement* TkConvValidator::p_dzPVVsR_
private

Definition at line 296 of file TkConvValidator.h.

◆ p_nHitsVsEta_

MonitorElement* TkConvValidator::p_nHitsVsEta_[3]
private

Definition at line 311 of file TkConvValidator.h.

◆ p_nHitsVsR_

MonitorElement* TkConvValidator::p_nHitsVsR_[3]
private

Definition at line 313 of file TkConvValidator.h.

◆ p_TkPtPull_

MonitorElement* TkConvValidator::p_TkPtPull_[3]
private

Definition at line 326 of file TkConvValidator.h.

◆ parameters_

edm::ParameterSet TkConvValidator::parameters_
private

Definition at line 86 of file TkConvValidator.h.

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

◆ photonCollection_

std::string TkConvValidator::photonCollection_
private

Definition at line 97 of file TkConvValidator.h.

◆ photonCollectionPr_Token_

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

Definition at line 98 of file TkConvValidator.h.

◆ photonCollectionProducer_

std::string TkConvValidator::photonCollectionProducer_
private

Definition at line 96 of file TkConvValidator.h.

◆ recMaxPt_

double TkConvValidator::recMaxPt_
private

Definition at line 174 of file TkConvValidator.h.

◆ recMinPt_

double TkConvValidator::recMinPt_
private

Global variables for reco Photon.

Definition at line 173 of file TkConvValidator.h.

◆ simMaxPt_

double TkConvValidator::simMaxPt_
private

Definition at line 170 of file TkConvValidator.h.

◆ simMinPt_

double TkConvValidator::simMinPt_
private

Definition at line 169 of file TkConvValidator.h.

◆ theCaloGeom_

edm::ESHandle<CaloGeometry> TkConvValidator::theCaloGeom_
private

Definition at line 87 of file TkConvValidator.h.

◆ theCaloTopo_

edm::ESHandle<CaloTopology> TkConvValidator::theCaloTopo_
private

Definition at line 88 of file TkConvValidator.h.

◆ theConvTP_

edm::RefVector<TrackingParticleCollection> TkConvValidator::theConvTP_
private

Definition at line 157 of file TkConvValidator.h.

◆ theMF_

edm::ESHandle<MagneticField> TkConvValidator::theMF_
private

Definition at line 73 of file TkConvValidator.h.

◆ thePhotonMCTruthFinder_

PhotonMCTruthFinder* TkConvValidator::thePhotonMCTruthFinder_
private

Definition at line 114 of file TkConvValidator.h.

◆ tpSelForEff_Token_

edm::EDGetTokenT<TrackingParticleRefVector> TkConvValidator::tpSelForEff_Token_
private

Definition at line 104 of file TkConvValidator.h.

◆ tpSelForFake_Token_

edm::EDGetTokenT<TrackingParticleRefVector> TkConvValidator::tpSelForFake_Token_
private

Definition at line 105 of file TkConvValidator.h.

◆ trackAssociator_Token_

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

Definition at line 108 of file TkConvValidator.h.

◆ trkIsolExtRadius_

double TkConvValidator::trkIsolExtRadius_
private

Definition at line 119 of file TkConvValidator.h.

◆ trkIsolInnRadius_

double TkConvValidator::trkIsolInnRadius_
private

Definition at line 120 of file TkConvValidator.h.

◆ trkPtLow_

double TkConvValidator::trkPtLow_
private

Definition at line 121 of file TkConvValidator.h.

◆ trkPtSumCut_

double TkConvValidator::trkPtSumCut_
private

Definition at line 129 of file TkConvValidator.h.

◆ verbosity_

int TkConvValidator::verbosity_
private

Definition at line 75 of file TkConvValidator.h.

TkConvValidator::h_convR_
MonitorElement * h_convR_[3][3]
Definition: TkConvValidator.h:216
TkConvValidator::p_convVtxdZVsZ_
MonitorElement * p_convVtxdZVsZ_
Definition: TkConvValidator.h:283
reco::Conversion::pairCotThetaSeparation
double pairCotThetaSeparation() const
Delta cot(Theta) where Theta is the angle in the (y,z) plane between the two tracks....
Definition: Conversion.cc:179
PI
Definition: PayloadInspector.h:19
TkConvValidator::h_convEta_
MonitorElement * h_convEta_[3][3]
Definition: TkConvValidator.h:213
TkRotation< float >
reco::Vertex::isValid
bool isValid() const
Tells whether the vertex is valid.
Definition: Vertex.h:72
TkConvValidator::h_SimRecConvTwoMTracks_
MonitorElement * h_SimRecConvTwoMTracks_[5]
Definition: TkConvValidator.h:203
reco::Conversion
Definition: Conversion.h:23
photonAnalyzer_cfi.rMax
rMax
Definition: photonAnalyzer_cfi.py:91
photonAnalyzer_cfi.etMin
etMin
Definition: photonAnalyzer_cfi.py:54
TkConvValidator::h_TkD0_
MonitorElement * h_TkD0_[3]
Definition: TkConvValidator.h:322
edm::RefVector::clear
void clear()
Clear the vector.
Definition: RefVector.h:142
TkConvValidator::maxPhoZForPurity
double maxPhoZForPurity
Definition: TkConvValidator.h:166
PDWG_EXOHSCP_cff.tracks
tracks
Definition: PDWG_EXOHSCP_cff.py:28
TkConvValidator::arbitratedEcalSeeded_
bool arbitratedEcalSeeded_
Definition: TkConvValidator.h:137
TkConvValidator::theMF_
edm::ESHandle< MagneticField > theMF_
Definition: TkConvValidator.h:73
ALCARECOTkAlBeamHalo_cff.etaMin
etaMin
GeV.
Definition: ALCARECOTkAlBeamHalo_cff.py:32
anyDirection
Definition: PropagationDirection.h:4
mps_fire.i
i
Definition: mps_fire.py:355
TkConvValidator::h_distMinAppTracks_
MonitorElement * h_distMinAppTracks_[3][3]
Definition: TkConvValidator.h:243
TkConvValidator::mcConvEta_
double mcConvEta_
Definition: TkConvValidator.h:153
TkConvValidator::p_DPhiTracksAtVtxVsEta_
MonitorElement * p_DPhiTracksAtVtxVsEta_
Definition: TkConvValidator.h:233
reco::Photon::superCluster
reco::SuperClusterRef superCluster() const override
Ref to SuperCluster.
reco::Conversion::generalTracksOnly
Definition: Conversion.h:28
photonAnalyzer_cfi.zMax
zMax
Definition: photonAnalyzer_cfi.py:95
TkConvValidator::h2_Chi2VsR_
MonitorElement * h2_Chi2VsR_[3]
Definition: TkConvValidator.h:319
TkConvValidator::h_convEta2_
MonitorElement * h_convEta2_[3][3]
Definition: TkConvValidator.h:214
TkConvValidator::h_convVtxdR_barrel_
MonitorElement * h_convVtxdR_barrel_
Definition: TkConvValidator.h:268
edm::Handle::product
T const * product() const
Definition: Handle.h:70
TkConvValidator::h_convVtxYvsX_
MonitorElement * h_convVtxYvsX_
Definition: TkConvValidator.h:254
photonAnalyzer_cfi.dPhiTracksMin
dPhiTracksMin
Definition: photonAnalyzer_cfi.py:102
etaBin
int etaBin(const l1t::HGCalMulticluster *cl)
Definition: L1EGammaEEProducer.cc:19
TkConvValidator::h_SimConvTwoTracks_
MonitorElement * h_SimConvTwoTracks_[5]
Definition: TkConvValidator.h:194
reco::TrackBase::inOutEcalSeededConv
Definition: TrackBase.h:106
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
photonAnalyzer_cfi.zBin
zBin
Definition: photonAnalyzer_cfi.py:93
TkConvValidator::h_convEtaMatchSC_
MonitorElement * h_convEtaMatchSC_[3][3]
Definition: TkConvValidator.h:212
reco::Conversion::tracks
std::vector< edm::RefToBase< reco::Track > > const & tracks() const
vector of track to base references
Definition: Conversion.cc:150
TkConvValidator::mcConvPt_
double mcConvPt_
Definition: TkConvValidator.h:147
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
photonAnalyzer_cfi.dPhiTracksBin
dPhiTracksBin
Definition: photonAnalyzer_cfi.py:101
PI
#define PI
Definition: QcdUeDQM.h:37
conv
static HepMC::IO_HEPEVT conv
Definition: BeamHaloProducer.cc:48
TkConvValidator::h_convSCdPhi_
MonitorElement * h_convSCdPhi_[3][3]
Definition: TkConvValidator.h:222
TkConvValidator::h_VisSimConv_
MonitorElement * h_VisSimConv_[6]
Definition: TkConvValidator.h:189
BeamMonitor_cff.phiBin
phiBin
Definition: BeamMonitor_cff.py:75
photonValidationSequence_cff.eoverpMax
eoverpMax
Definition: photonValidationSequence_cff.py:17
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
ChiSquaredProbability
float ChiSquaredProbability(double chiSquared, double nrDOF)
Definition: ChiSquaredProbability.cc:13
CaloGeometryRecord
Definition: CaloGeometryRecord.h:30
TkConvValidator::h_deltaExpectedHitsInner_
MonitorElement * h_deltaExpectedHitsInner_[3][3]
Definition: TkConvValidator.h:303
edm
HLT enums.
Definition: AlignableModifier.h:19
TkConvValidator::h_simConvVtxYvsX_
MonitorElement * h_simConvVtxYvsX_
Definition: TkConvValidator.h:185
reco::TrackToTrackingParticleAssociator::associateRecoToSim
reco::RecoToSimCollection associateRecoToSim(const edm::Handle< edm::View< reco::Track >> &tCH, const edm::Handle< TrackingParticleCollection > &tPCH) const
Definition: TrackToTrackingParticleAssociator.h:64
photonValidationSequence_cff.eoverpMin
eoverpMin
Definition: photonValidationSequence_cff.py:16
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
TkConvValidator::nRecConvAssWithEcal_
int nRecConvAssWithEcal_
Definition: TkConvValidator.h:82
edm::RefVector::begin
const_iterator begin() const
Initialize an iterator over the RefVector.
Definition: RefVector.h:223
TrackCollections2monitor_cff.numberOfValidHits
numberOfValidHits
Definition: TrackCollections2monitor_cff.py:248
TkConvValidator::h_convRplot_
MonitorElement * h_convRplot_
Definition: TkConvValidator.h:224
edm::LogInfo
Definition: MessageLogger.h:254
TkConvValidator::h2_DCotTracksVsEta_
MonitorElement * h2_DCotTracksVsEta_
Definition: TkConvValidator.h:238
TkConvValidator::h2_convVtxdRVsEta_
MonitorElement * h2_convVtxdRVsEta_
Definition: TkConvValidator.h:276
muonRecoAnalyzer_cfi.chi2Min
chi2Min
Definition: muonRecoAnalyzer_cfi.py:34
TkConvValidator::h_simTkEta_
MonitorElement * h_simTkEta_
Definition: TkConvValidator.h:182
TkConvValidator::nEntry_
int nEntry_
Definition: TkConvValidator.h:77
TkConvValidator::p_DCotTracksVsEta_
MonitorElement * p_DCotTracksVsEta_
Definition: TkConvValidator.h:239
TkConvValidator::h2_DPhiTracksAtVtxVsR_
MonitorElement * h2_DPhiTracksAtVtxVsR_
Definition: TkConvValidator.h:234
reco::Conversion::pairInvariantMass
double pairInvariantMass() const
if nTracks=2 returns the pair invariant mass. Original tracks are used here
Definition: Conversion.cc:159
reco::Conversion::zOfPrimaryVertexFromTracks
double zOfPrimaryVertexFromTracks(const math::XYZPoint &myBeamSpot=math::XYZPoint()) const
Definition: Conversion.h:136
TkConvValidator::tpSelForEff_Token_
edm::EDGetTokenT< TrackingParticleRefVector > tpSelForEff_Token_
Definition: TkConvValidator.h:104
reco::Conversion::dlClosestHitToVtx
const std::vector< Measurement1DFloat > & dlClosestHitToVtx() const
Vector of signed decay length with uncertainty from nearest hit on track to the conversion vtx positi...
Definition: Conversion.h:156
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:114
TkConvValidator::h_convVtxdPhi_
MonitorElement * h_convVtxdPhi_
Definition: TkConvValidator.h:263
TkConvValidator::h_EoverPTracks_
MonitorElement * h_EoverPTracks_[3][3]
Definition: TkConvValidator.h:219
photonValidator_cfi.dCotTracksMax
dCotTracksMax
Definition: photonValidator_cfi.py:100
TkConvValidator::p_TkPtPull_
MonitorElement * p_TkPtPull_[3]
Definition: TkConvValidator.h:326
TkConvValidator::h2_dzPVVsR_
MonitorElement * h2_dzPVVsR_
Definition: TkConvValidator.h:295
ReferenceCountingPointer< Cylinder >
TkConvValidator::h_convVtxdR_endcap_
MonitorElement * h_convVtxdR_endcap_
Definition: TkConvValidator.h:273
edm::RefVector< TrackingParticleCollection >
PhotonMCTruthFinder::find
std::vector< PhotonMCTruth > find(const std::vector< SimTrack > &simTracks, const std::vector< SimVertex > &simVertices)
Definition: PhotonMCTruthFinder.cc:17
photonValidator_cfi.zMinForXray
zMinForXray
Definition: photonValidator_cfi.py:122
TkConvValidator::h2_photonPtRecVsPtSim_
MonitorElement * h2_photonPtRecVsPtSim_
Definition: TkConvValidator.h:328
TkConvValidator::p_dzPVVsR_
MonitorElement * p_dzPVVsR_
Definition: TkConvValidator.h:296
reco::HitPattern::numberOfLostHits
int numberOfLostHits(HitCategory category) const
Definition: HitPattern.h:860
edm::Handle< reco::TrackToTrackingParticleAssociator >
parallelization.uint
uint
Definition: parallelization.py:124
TkConvValidator::parameters_
edm::ParameterSet parameters_
Definition: TkConvValidator.h:86
photonValidator_cfi.rMinForXray
rMinForXray
Definition: photonValidator_cfi.py:118
TkConvValidator::nHits_
MonitorElement * nHits_[3]
Definition: TkConvValidator.h:310
TkConvValidator::conversionCollectionPr_Token_
edm::EDGetTokenT< reco::ConversionCollection > conversionCollectionPr_Token_
Definition: TkConvValidator.h:92
TkConvValidator::h_TkPtPull_
MonitorElement * h_TkPtPull_[3]
Definition: TkConvValidator.h:324
TkConvValidator::beamspotToken_
edm::EDGetTokenT< reco::BeamSpot > beamspotToken_
Definition: TkConvValidator.h:101
TkConvValidator::h_SimConvTwoMTracksAndVtxPGT0_
MonitorElement * h_SimConvTwoMTracksAndVtxPGT0_[5]
Definition: TkConvValidator.h:196
TkConvValidator::h_nSimConv_
MonitorElement * h_nSimConv_[2]
Definition: TkConvValidator.h:178
edm::Ref< TrackingParticleCollection >
reco::ConversionCollection
std::vector< Conversion > ConversionCollection
collectin of Conversion objects
Definition: ConversionFwd.h:9
recPt
double recPt
Definition: PFJetBenchmarkAnalyzer.cc:75
TkConvValidator::isRunCentrally_
bool isRunCentrally_
Definition: TkConvValidator.h:116
edm::Exception
Definition: EDMException.h:77
photonAnalyzer_cfi.dPhiTracksMax
dPhiTracksMax
Definition: photonAnalyzer_cfi.py:103
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
SimpleCylinderBounds
Definition: SimpleCylinderBounds.h:19
PhotonMCTruthFinder::clear
void clear()
Definition: PhotonMCTruthFinder.h:27
TkConvValidator::fName_
std::string fName_
Definition: TkConvValidator.h:71
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:608
TkConvValidator::generalTracksOnly_
bool generalTracksOnly_
Definition: TkConvValidator.h:135
TkConvValidator::photonCollection_
std::string photonCollection_
Definition: TkConvValidator.h:97
IdealMagneticFieldRecord
Definition: IdealMagneticFieldRecord.h:11
SimpleDiskBounds
Definition: SimpleDiskBounds.h:11
conversionPostprocessing_cfi.etaBin2
etaBin2
Definition: conversionPostprocessing_cfi.py:19
TrajectoryStateOnSurface
Definition: TrajectoryStateOnSurface.h:16
reco::Conversion::refittedPairMomentum
math::XYZVectorF refittedPairMomentum() const
Conversion tracks momentum from the tracks refitted with vertex constraint.
Definition: Conversion.cc:206
TkConvValidator::tpSelForFake_Token_
edm::EDGetTokenT< TrackingParticleRefVector > tpSelForFake_Token_
Definition: TkConvValidator.h:105
reco::Conversion::arbitratedMerged
Definition: Conversion.h:30
TkConvValidator::h_dlClosestHitToVtxSig_
MonitorElement * h_dlClosestHitToVtxSig_[3]
Definition: TkConvValidator.h:336
edm::RefVector::end
const_iterator end() const
Termination of iteration.
Definition: RefVector.h:228
PropagatorWithMaterial
Definition: PropagatorWithMaterial.h:25
TkConvValidator::arbitratedMerged_
bool arbitratedMerged_
Definition: TkConvValidator.h:136
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
photonValidator_cfi.resMin
resMin
Definition: photonValidator_cfi.py:91
reco::Conversion::nHitsBeforeVtx
const std::vector< uint8_t > & nHitsBeforeVtx() const
Vector of the number of hits before the vertex along each track trajector.
Definition: Conversion.h:154
TkConvValidator::h2_convVtxRrecVsTrue_
MonitorElement * h2_convVtxRrecVsTrue_
Definition: TkConvValidator.h:289
photonValidationSequence_cff.useTP
useTP
Definition: photonValidationSequence_cff.py:14
edm::RefToBaseVector::push_back
void push_back(const RefToBase< T > &)
Definition: RefToBaseVector.h:217
reco::Conversion::arbitratedEcalSeeded
Definition: Conversion.h:29
TkConvValidator::h_convVtxdZ_endcap_
MonitorElement * h_convVtxdZ_endcap_
Definition: TkConvValidator.h:272
reco::TrackBase::outInEcalSeededConv
Definition: TrackBase.h:105
TkConvValidator::h_leadExpectedHitsInner_
MonitorElement * h_leadExpectedHitsInner_[3][3]
Definition: TkConvValidator.h:304
TkConvValidator::nRecConvAss_
int nRecConvAss_
Definition: TkConvValidator.h:81
TrackerDigiGeometryRecord
Definition: TrackerDigiGeometryRecord.h:15
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
TkConvValidator::p2_convVtxdZVsRZ_
MonitorElement * p2_convVtxdZVsRZ_
Definition: TkConvValidator.h:287
TkConvValidator::nHitsVsR_
MonitorElement * nHitsVsR_[3]
Definition: TkConvValidator.h:314
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
TkConvValidator::h_convVtxdX_
MonitorElement * h_convVtxdX_
Definition: TkConvValidator.h:258
AlignmentTrackSelector_cfi.phiMin
phiMin
Definition: AlignmentTrackSelector_cfi.py:18
photonValidator_cfi.rMaxForXray
rMaxForXray
Definition: photonValidator_cfi.py:119
reco::BeamSpot
Definition: BeamSpot.h:21
TkConvValidator::etaTransformation
float etaTransformation(float a, float b)
Definition: TkConvValidator.cc:2294
TransientTrackRecord
Definition: TransientTrackRecord.h:11
TkConvValidator::h_convSCdEta_
MonitorElement * h_convSCdEta_[3][3]
Definition: TkConvValidator.h:221
reco::Track
Definition: Track.h:27
edm::ESHandle< TransientTrackBuilder >
AlignmentTrackSelector_cfi.phiMax
phiMax
Definition: AlignmentTrackSelector_cfi.py:17
p2
double p2[4]
Definition: TauolaWrapper.h:90
TkConvValidator::h_convVtxdY_endcap_
MonitorElement * h_convVtxdY_endcap_
Definition: TkConvValidator.h:271
photonValidator_cfi.zMaxForXray
zMaxForXray
Definition: photonValidator_cfi.py:123
TkConvValidator::h2_Chi2VsEta_
MonitorElement * h2_Chi2VsEta_[3]
Definition: TkConvValidator.h:317
TkConvValidator::h_convPhi_
MonitorElement * h_convPhi_[3][3]
Definition: TkConvValidator.h:215
reco_skim_cfg_mod.outputFileName
outputFileName
Definition: reco_skim_cfg_mod.py:15
TkConvValidator::mcConvX_
double mcConvX_
Definition: TkConvValidator.h:151
TkConvValidator::h_sumNHitsBeforeVtx_
MonitorElement * h_sumNHitsBeforeVtx_[3][3]
Definition: TkConvValidator.h:302
Point3DBase< float, GlobalTag >
TkConvValidator::minPhoPtForPurity
double minPhoPtForPurity
Definition: TkConvValidator.h:164
dqm::implementation::DQMStore::save
DQM_DEPRECATED void save(std::string const &filename, std::string const &path="")
Definition: DQMStore.cc:760
PhotonMCTruthFinder
Definition: PhotonMCTruthFinder.h:20
TkConvValidator::p_nHitsVsR_
MonitorElement * p_nHitsVsR_[3]
Definition: TkConvValidator.h:313
edm::RefToBase::get
value_type const * get() const
Definition: RefToBase.h:209
q2
double q2[4]
Definition: TauolaWrapper.h:88
cmsswSequenceInfo.tp
tp
Definition: cmsswSequenceInfo.py:17
reco::Vertex::refittedTracks
const std::vector< Track > & refittedTracks() const
Returns the container of refitted tracks.
Definition: Vertex.h:183
TkConvValidator::recalculateMomentumAtFittedVertex
math::XYZVector recalculateMomentumAtFittedVertex(const MagneticField &mf, const TrackerGeometry &trackerGeom, const edm::RefToBase< reco::Track > &tk, const reco::Vertex &vtx)
Definition: TkConvValidator.cc:2244
TkConvValidator::h_lxybs_
MonitorElement * h_lxybs_[3][3]
Definition: TkConvValidator.h:298
photonValidator_cfi.eoverpBin
eoverpBin
Definition: photonValidator_cfi.py:106
badGlobalMuonTaggersAOD_cff.vtx
vtx
Definition: badGlobalMuonTaggersAOD_cff.py:5
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
TkConvValidator::verbosity_
int verbosity_
Definition: TkConvValidator.h:75
TkConvValidator::g4_simTk_Token_
edm::EDGetTokenT< edm::SimTrackContainer > g4_simTk_Token_
Definition: TkConvValidator.h:102
TkConvValidator::h_nHitsBeforeVtx_
MonitorElement * h_nHitsBeforeVtx_[3]
Definition: TkConvValidator.h:334
TkConvValidator::maxPhoEtaForPurity
double maxPhoEtaForPurity
Definition: TkConvValidator.h:165
TkConvValidator::nEvt_
int nEvt_
Definition: TkConvValidator.h:76
TkConvValidator::h_simConvVtxRvsZ_
MonitorElement * h_simConvVtxRvsZ_[4]
Definition: TkConvValidator.h:184
TkConvValidator::h_DPhiTracksAtVtx_
MonitorElement * h_DPhiTracksAtVtx_[3][3]
Definition: TkConvValidator.h:231
PbPb_ZMuSkimMuonDPG_cff.tracker
tracker
Definition: PbPb_ZMuSkimMuonDPG_cff.py:60
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
TkConvValidator::h_convVtxdX_barrel_
MonitorElement * h_convVtxdX_barrel_
Definition: TkConvValidator.h:265
TkConvValidator::h_simTkPt_
MonitorElement * h_simTkPt_
Definition: TkConvValidator.h:181
TkConvValidator::h_maxDlClosestHitToVtx_
MonitorElement * h_maxDlClosestHitToVtx_[3][3]
Definition: TkConvValidator.h:305
TkConvValidator::p_convVtxdRVsR_
MonitorElement * p_convVtxdRVsR_
Definition: TkConvValidator.h:278
TkConvValidator::mcConvZ_
double mcConvZ_
Definition: TkConvValidator.h:149
TkConvValidator::conversionCollectionProducer_
std::string conversionCollectionProducer_
Definition: TkConvValidator.h:90
q1
double q1[4]
Definition: TauolaWrapper.h:87
TkConvValidator::p_Chi2VsEta_
MonitorElement * p_Chi2VsEta_[3]
Definition: TkConvValidator.h:318
edm::LogError
Definition: MessageLogger.h:183
reco::Conversion::quality
bool quality(ConversionQuality q) const
Definition: Conversion.h:178
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
TWOPI
#define TWOPI
Definition: DQMSourcePi0.cc:36
TkConvValidator::photonCollectionProducer_
std::string photonCollectionProducer_
Definition: TkConvValidator.h:96
reco::Conversion::conversionVertex
const reco::Vertex & conversionVertex() const
returns the reco conversion vertex
Definition: Conversion.h:88
edm::AssociationMap< edm::OneToManyWithQualityGeneric< TrackingParticleCollection, edm::View< reco::Track >, double > >
TkConvValidator::dqmpath_
std::string dqmpath_
Definition: TkConvValidator.h:110
reco::Vertex::x
double x() const
x coordinate
Definition: Vertex.h:116
TkConvValidator::h_convPt_
MonitorElement * h_convPt_[3][3]
Definition: TkConvValidator.h:218
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:31
TkConvValidator::offline_pvToken_
edm::EDGetTokenT< reco::VertexCollection > offline_pvToken_
Definition: TkConvValidator.h:100
photonAnalyzer_cfi.zMin
zMin
Definition: photonAnalyzer_cfi.py:94
TkConvValidator::h_zPVFromTracks_
MonitorElement * h_zPVFromTracks_[2]
Definition: TkConvValidator.h:293
TkConvValidator::h_convVtxYvsX_zoom_
MonitorElement * h_convVtxYvsX_zoom_[2]
Definition: TkConvValidator.h:256
TkConvValidator::maxPhoRForPurity
double maxPhoRForPurity
Definition: TkConvValidator.h:167
reco::Conversion::highPurity
Definition: Conversion.h:33
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
photonValidator_cfi.dCotTracksMin
dCotTracksMin
Definition: photonValidator_cfi.py:99
edm::Service
Definition: Service.h:30
TkConvValidator::h_trailNHitsBeforeVtx_
MonitorElement * h_trailNHitsBeforeVtx_[3][3]
Definition: TkConvValidator.h:301
TkConvValidator::p_DPhiTracksAtVtxVsR_
MonitorElement * p_DPhiTracksAtVtxVsR_
Definition: TkConvValidator.h:235
photonAnalyzer_cfi.rBin
rBin
Definition: photonAnalyzer_cfi.py:89
muonRecoAnalyzer_cfi.chi2Max
chi2Max
Definition: muonRecoAnalyzer_cfi.py:35
reco::Conversion::distOfMinimumApproach
double distOfMinimumApproach() const
Definition: Conversion.h:116
TkConvValidator::h_convVtxRvsZ_zoom_
MonitorElement * h_convVtxRvsZ_zoom_[2]
Definition: TkConvValidator.h:255
TkConvValidator::maxPhoZForEffic
double maxPhoZForEffic
Definition: TkConvValidator.h:162
TkConvValidator::h_convVtxdEta_
MonitorElement * h_convVtxdEta_
Definition: TkConvValidator.h:262
photonValidator_cfi.resBin
resBin
Definition: photonValidator_cfi.py:90
p1
double p1[4]
Definition: TauolaWrapper.h:89
EgammaValidation_cff.pdgId
pdgId
Definition: EgammaValidation_cff.py:118
TkConvValidator::h_SimConvTwoMTracks_
MonitorElement * h_SimConvTwoMTracks_[5]
Definition: TkConvValidator.h:195
TkConvValidator::conversionCollection_
std::string conversionCollection_
Definition: TkConvValidator.h:91
TkConvValidator::p_convVtxdRVsEta_
MonitorElement * p_convVtxdRVsEta_
Definition: TkConvValidator.h:279
TkConvValidator::mcConvR_
double mcConvR_
Definition: TkConvValidator.h:148
photonValidator_cfi.zBin2ForXray
zBin2ForXray
Definition: photonValidator_cfi.py:121
TkConvValidator::h_DCotTracks_
MonitorElement * h_DCotTracks_[3][3]
Definition: TkConvValidator.h:237
TkConvValidator::h_dlClosestHitToVtx_
MonitorElement * h_dlClosestHitToVtx_[3]
Definition: TkConvValidator.h:335
get
#define get
edm::RefVector::push_back
void push_back(value_type const &ref)
Add a Ref<C, T> to the RefVector.
Definition: RefVector.h:67
edm::RefToBaseVector< reco::Track >
TrackCollections2monitor_cff.normalizedChi2
normalizedChi2
Definition: TrackCollections2monitor_cff.py:247
TkConvValidator::h_dzPVFromTracks_
MonitorElement * h_dzPVFromTracks_[2]
Definition: TkConvValidator.h:294
TkConvValidator::p_convVtxdYVsY_
MonitorElement * p_convVtxdYVsY_
Definition: TkConvValidator.h:282
nanoDQM_cff.Photon
Photon
Definition: nanoDQM_cff.py:63
TkConvValidator::h_match_
MonitorElement * h_match_
Definition: TkConvValidator.h:330
TkConvValidator::h_vtxChi2Prob_
MonitorElement * h_vtxChi2Prob_[3][3]
Definition: TkConvValidator.h:291
reco::Vertex::chi2
double chi2() const
chi-squares
Definition: Vertex.h:103
TkConvValidator::h_convVtxdY_
MonitorElement * h_convVtxdY_
Definition: TkConvValidator.h:259
TkConvValidator::dbe_
DQMStore * dbe_
Definition: TkConvValidator.h:72
TkConvValidator::mcConvY_
double mcConvY_
Definition: TkConvValidator.h:150
TkConvValidator::nRecConv_
int nRecConv_
Definition: TkConvValidator.h:80
TkConvValidator::theConvTP_
edm::RefVector< TrackingParticleCollection > theConvTP_
Definition: TkConvValidator.h:157
TkConvValidator::ecalalgotracks_
bool ecalalgotracks_
Definition: TkConvValidator.h:138
reco::TrackBase::algo
TrackAlgorithm algo() const
Definition: TrackBase.h:532
TkConvValidator::h_AllSimConv_
MonitorElement * h_AllSimConv_[5]
Denominator for efficiencies.
Definition: TkConvValidator.h:188
TkConvValidator::h_convVtxdR_
MonitorElement * h_convVtxdR_
Definition: TkConvValidator.h:261
DDAxes::phi
TkConvValidator::h2_TkPtPull_
MonitorElement * h2_TkPtPull_[3]
Definition: TkConvValidator.h:325
TkConvValidator::h_maxNHitsBeforeVtx_
MonitorElement * h_maxNHitsBeforeVtx_[3][3]
Definition: TkConvValidator.h:299
reco::Photon
Definition: Photon.h:21
TkConvValidator::p_DCotTracksVsR_
MonitorElement * p_DCotTracksVsR_
Definition: TkConvValidator.h:241
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
cosmicPhotonAnalyzer_cfi.etMax
etMax
Definition: cosmicPhotonAnalyzer_cfi.py:11
TkConvValidator::h_convZplot_
MonitorElement * h_convZplot_
Definition: TkConvValidator.h:225
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:489
reco::BeamSpot::x0
double x0() const
x coordinate
Definition: BeamSpot.h:61
edm::AssociationMap::empty
bool empty() const
return true if empty
Definition: AssociationMap.h:164
reco::Conversion::nSharedHits
uint8_t nSharedHits() const
number of shared hits btw the two track
Definition: Conversion.h:158
TkConvValidator::p2_convVtxdRVsRZ_
MonitorElement * p2_convVtxdRVsRZ_
Definition: TkConvValidator.h:286
TrajectoryStateOnSurface::globalMomentum
GlobalVector globalMomentum() const
Definition: TrajectoryStateOnSurface.h:66
TkConvValidator::p_nHitsVsEta_
MonitorElement * p_nHitsVsEta_[3]
Definition: TkConvValidator.h:311
TkConvValidator::h_convVtxdY_barrel_
MonitorElement * h_convVtxdY_barrel_
Definition: TkConvValidator.h:266
TkConvValidator::photonCollectionPr_Token_
edm::EDGetTokenT< reco::PhotonCollection > photonCollectionPr_Token_
Definition: TkConvValidator.h:98
Disk
Definition: BoundDisk.h:19
reco::Vertex::y
double y() const
y coordinate
Definition: Vertex.h:118
photonValidator_cfi.dCotTracksBin
dCotTracksBin
Definition: photonValidator_cfi.py:98
TkConvValidator::minPhoPtForEffic
double minPhoPtForEffic
Definition: TkConvValidator.h:160
TkConvValidator::h_convVtxdX_endcap_
MonitorElement * h_convVtxdX_endcap_
Definition: TkConvValidator.h:270
TkConvValidator::h_convVtxdZ_
MonitorElement * h_convVtxdZ_
Definition: TkConvValidator.h:260
spclusmultinvestigator_cfi.vertexCollection
vertexCollection
Definition: spclusmultinvestigator_cfi.py:4
TkConvValidator::mcEta_
double mcEta_
Definition: TkConvValidator.h:146
edm::RefVectorIterator
Definition: EDProductfwd.h:33
ETA
#define ETA
Definition: GenericBenchmark.cc:28
reco::HitPattern::MISSING_INNER_HITS
Definition: HitPattern.h:155
TkConvValidator::g4_simVtx_Token_
edm::EDGetTokenT< edm::SimVertexContainer > g4_simVtx_Token_
Definition: TkConvValidator.h:103
TkConvValidator::maxPhoEtaForEffic
double maxPhoEtaForEffic
Definition: TkConvValidator.h:161
makeMuonMisalignmentScenario.rot
rot
Definition: makeMuonMisalignmentScenario.py:322
math::XYZVectorF
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< float > > XYZVectorF
spatial vector with cartesian internal representation
Definition: Vector3D.h:16
TkConvValidator::h2_convVtxdRVsR_
MonitorElement * h2_convVtxdRVsR_
Definition: TkConvValidator.h:275
ALCARECOTkAlBeamHalo_cff.etaMax
etaMax
Definition: ALCARECOTkAlBeamHalo_cff.py:33
TkConvValidator::thePhotonMCTruthFinder_
PhotonMCTruthFinder * thePhotonMCTruthFinder_
Definition: TkConvValidator.h:114
TkConvValidator::h_convVtxdZ_barrel_
MonitorElement * h_convVtxdZ_barrel_
Definition: TkConvValidator.h:267
edm::RefToBase< reco::Track >
photonValidator_cfi.resMax
resMax
Definition: photonValidator_cfi.py:92
TkConvValidator::highPurity_
bool highPurity_
Definition: TkConvValidator.h:139
TkConvValidator::h_SimConvEtaPix_
MonitorElement * h_SimConvEtaPix_[2]
Definition: TkConvValidator.h:179
reco::TrackToTrackingParticleAssociator::associateSimToReco
reco::SimToRecoCollection associateSimToReco(const edm::Handle< edm::View< reco::Track >> &tCH, const edm::Handle< TrackingParticleCollection > &tPCH) const
Definition: TrackToTrackingParticleAssociator.h:71
TkConvValidator::maxHitsBeforeVtx_
uint maxHitsBeforeVtx_
Definition: TkConvValidator.h:141
TkConvValidator::p_Chi2VsR_
MonitorElement * p_Chi2VsR_[3]
Definition: TkConvValidator.h:320
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
photonValidator_cfi.rBinForXray
rBinForXray
Definition: photonValidator_cfi.py:117
TkConvValidator::h_nSharedHits_
MonitorElement * h_nSharedHits_[3][3]
Definition: TkConvValidator.h:307
Cylinder::computeRadius
static float computeRadius(Bounds const &bounds)
Definition: Cylinder.h:28
TkConvValidator::nInvalidPCA_
int nInvalidPCA_
Definition: TkConvValidator.h:84
edm::Ref::key
key_type key() const
Accessor for product key.
Definition: Ref.h:250
reco::PhotonCollection
std::vector< Photon > PhotonCollection
collectin of Photon objects
Definition: PhotonFwd.h:9
photonAnalyzer_cfi.rMin
rMin
Definition: photonAnalyzer_cfi.py:90
Cylinder
Definition: Cylinder.h:19
TkConvValidator::minLxy_
double minLxy_
Definition: TkConvValidator.h:142
mps_fire.result
result
Definition: mps_fire.py:303
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
TkConvValidator::h_SimConvTwoMTracksAndVtxPGT0005_
MonitorElement * h_SimConvTwoMTracksAndVtxPGT0005_[5]
Definition: TkConvValidator.h:197
TkConvValidator::h_nConv_
MonitorElement * h_nConv_[3][3]
info per conversion
Definition: TkConvValidator.h:211
HGC3DClusterGenMatchSelector_cfi.dR
dR
Definition: HGC3DClusterGenMatchSelector_cfi.py:7
TkConvValidator::theCaloGeom_
edm::ESHandle< CaloGeometry > theCaloGeom_
Definition: TkConvValidator.h:87
ExoticaDQM_cfi.photonCollection
photonCollection
Definition: ExoticaDQM_cfi.py:17
dummy
Definition: DummySelector.h:38
TkConvValidator::p_convVtxdXVsX_
MonitorElement * p_convVtxdXVsX_
Definition: TkConvValidator.h:281
TkConvValidator::p_convVtxdZVsR_
MonitorElement * p_convVtxdZVsR_
Definition: TkConvValidator.h:284
TkConvValidator::h_leadNHitsBeforeVtx_
MonitorElement * h_leadNHitsBeforeVtx_[3][3]
Definition: TkConvValidator.h:300
TkConvValidator::h2_PtRecVsPtSim_
MonitorElement * h2_PtRecVsPtSim_[3]
Definition: TkConvValidator.h:327
trajectoryStateTransform::innerStateOnSurface
TrajectoryStateOnSurface innerStateOnSurface(const reco::Track &tk, const TrackingGeometry &geom, const MagneticField *field, bool withErr=true)
Definition: TrajectoryStateTransform.cc:110
TkConvValidator::trackAssociator_Token_
edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator > trackAssociator_Token_
Definition: TkConvValidator.h:108
ZEcal
static constexpr float ZEcal
Definition: L1TkEmParticleProducer.cc:40
TkConvValidator::nHitsVsEta_
MonitorElement * nHitsVsEta_[3]
Definition: TkConvValidator.h:312
photonAnalyzer_cfi.etBin
etBin
Definition: photonAnalyzer_cfi.py:53
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
event
Definition: event.py:1
reco::BeamSpot::y0
double y0() const
y coordinate
Definition: BeamSpot.h:63
TkConvValidator::phiNormalization
float phiNormalization(float &a)
Definition: TkConvValidator.cc:2278
edm::RefVector::size
size_type size() const
Size of the RefVector.
Definition: RefVector.h:102
d0
static constexpr float d0
Definition: L1EGammaCrystalsEmulatorProducer.cc:84
TkConvValidator::mcPhi_
double mcPhi_
global variable for the MC photon
Definition: TkConvValidator.h:145
TkConvValidator::minProb_
double minProb_
Definition: TkConvValidator.h:140
reco::Vertex::ndof
double ndof() const
Definition: Vertex.h:110
TkConvValidator::h_maxDlClosestHitToVtxSig_
MonitorElement * h_maxDlClosestHitToVtxSig_[3][3]
Definition: TkConvValidator.h:306
edm::InputTag
Definition: InputTag.h:15
TkConvValidator::maxPhoRForEffic
double maxPhoRForEffic
Definition: TkConvValidator.h:163
Z_Endcap
static constexpr float Z_Endcap
Definition: ECALPositionCalculator.cc:11
TkConvValidator::nSimConv_
int nSimConv_[2]
Definition: TkConvValidator.h:78
TkConvValidator::h_tkChi2Large_
MonitorElement * h_tkChi2Large_[3]
Definition: TkConvValidator.h:316
reco::Vertex
Definition: Vertex.h:35
R_ECAL
static constexpr float R_ECAL
Definition: ECALPositionCalculator.cc:10
TrajectoryStateOnSurface::isValid
bool isValid() const
Definition: TrajectoryStateOnSurface.h:54
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
TkConvValidator::minPhoEtCut_
double minPhoEtCut_
Definition: TkConvValidator.h:118
TkConvValidator::h_invMass_
MonitorElement * h_invMass_[3][3]
Definition: TkConvValidator.h:229
TkConvValidator::h_tkChi2_
MonitorElement * h_tkChi2_[3]
Definition: TkConvValidator.h:315
TkConvValidator::h_convZ_
MonitorElement * h_convZ_[3][3]
Definition: TkConvValidator.h:217
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
TkConvValidator::mcConvPhi_
double mcConvPhi_
Definition: TkConvValidator.h:152
TkConvValidator::h_convVtxRvsZ_
MonitorElement * h_convVtxRvsZ_[3]
Definition: TkConvValidator.h:253
TkConvValidator::h_convPtRes_
MonitorElement * h_convPtRes_[3]
Definition: TkConvValidator.h:227
TrackerGeometry
Definition: TrackerGeometry.h:14
etaBarrelEndcap
static constexpr float etaBarrelEndcap
Definition: ECALPositionCalculator.cc:12
TkConvValidator::h2_DPhiTracksAtVtxVsEta_
MonitorElement * h2_DPhiTracksAtVtxVsEta_
Definition: TkConvValidator.h:232
photonValidator_cfi.zBinForXray
zBinForXray
Definition: photonValidator_cfi.py:120
TkConvValidator::h2_DCotTracksVsR_
MonitorElement * h2_DCotTracksVsR_
Definition: TkConvValidator.h:240
Generated for CMSSW Reference Manual by doxygen 1.8.16