CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
MuonTrackValidator Class Reference

#include <MuonTrackValidator.h>

Inheritance diagram for MuonTrackValidator:
DQMEDAnalyzer MuonTrackValidatorBase edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > > edm::stream::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &)
 Method called before the event loop. More...
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
 
void endRun (edm::Run const &, edm::EventSetup const &)
 Method called at the end of the event loop. More...
 
 MuonTrackValidator (const edm::ParameterSet &pset)
 Constructor. More...
 
virtual ~MuonTrackValidator ()
 Destructor. More...
 
- Public Member Functions inherited from DQMEDAnalyzer
virtual void beginRun (edm::Run const &, edm::EventSetup const &) final
 
virtual void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer (void)
 
virtual void endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
virtual void endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
uint32_t streamId () const
 
- Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
 EDAnalyzer ()=default
 
- Public Member Functions inherited from edm::stream::EDAnalyzerBase
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDAnalyzerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Member Functions

void getRecoMomentum (const reco::Track &track, double &pt, double &ptError, double &qoverp, double &qoverpError, double &lambda, double &lambdaError, double &phi, double &phiError) const
 retrieval of reconstructed momentum components from reco::Track (== mean values for GSF) More...
 
void getRecoMomentum (const reco::GsfTrack &gsfTrack, double &pt, double &ptError, double &qoverp, double &qoverpError, double &lambda, double &lambdaError, double &phi, double &phiError) const
 retrieval of reconstructed momentum components based on the mode of a reco::GsfTrack More...
 

Private Attributes

edm::EDGetTokenT
< SimHitTPAssociationProducer::SimHitTPAssociationList > 		_simHitTpMapTag
 
edm::InputTag associatormap
 
bool BiDirectional_RecoToSim_association
 
std::vector< MonitorElement * > chi2_vs_eta
 
std::vector< MonitorElement * > chi2_vs_nhits
 
std::vector< MonitorElement * > chi2_vs_phi
 
CosmicTrackingParticleSelector cosmictpSelector
 
std::vector< MonitorElement * > cotThetares_vs_eta
 
std::vector< MonitorElement * > cotThetares_vs_pt
 
std::string dirName_
 
std::vector< MonitorElement * > dxypull_vs_eta
 
std::vector< MonitorElement * > dxyres_vs_eta
 
std::vector< MonitorElement * > dxyres_vs_pt
 
std::vector< MonitorElement * > dzpull_vs_eta
 
std::vector< MonitorElement * > dzres_vs_eta
 
std::vector< MonitorElement * > dzres_vs_pt
 
std::vector< MonitorElement * > etares_vs_eta
 
std::vector< MonitorElement * > h_assochi2
 
std::vector< MonitorElement * > h_assochi2_prob
 
std::vector< MonitorElement * > h_chi2mean_vs_phi
 
std::vector< MonitorElement * > h_chi2meanh
 
std::vector< MonitorElement * > h_chi2meanhitsh
 
std::vector< MonitorElement * > h_dxypulleta
 
std::vector< MonitorElement * > h_dzpulleta
 
std::vector< MonitorElement * > h_hits_phi
 
std::vector< MonitorElement * > h_losthits
 
std::vector< MonitorElement * > h_losthits_eta
 
std::vector< MonitorElement * > h_nchi2
 
std::vector< MonitorElement * > h_nchi2_prob
 
std::vector< MonitorElement * > h_phipulleta
 
std::vector< MonitorElement * > h_phipullphi
 
std::vector< MonitorElement * > h_ptpulleta
 
std::vector< MonitorElement * > h_ptpullphi
 
std::vector< MonitorElement * > h_ptshifteta
 
std::vector< MonitorElement * > h_thetapulleta
 
std::vector< MonitorElement * > h_thetapullphi
 
bool MABH
 
double maxPhi
 
double minPhi
 
std::vector< MonitorElement * > nhits_vs_phi
 
int nintPhi
 
std::vector< MonitorElement * > nlosthits_vs_eta
 
std::vector< MonitorElement * > phimean_vs_eta_phi
 
std::vector< MonitorElement * > phipull_vs_eta
 
std::vector< MonitorElement * > phipull_vs_phi
 
std::vector< MonitorElement * > phires_vs_eta
 
std::vector< MonitorElement * > phires_vs_phi
 
std::vector< MonitorElement * > phires_vs_pt
 
std::vector< MonitorElement * > ptmean_vs_eta_phi
 
std::vector< MonitorElement * > ptpull_vs_eta
 
std::vector< MonitorElement * > ptpull_vs_phi
 
std::vector< MonitorElement * > ptres_vs_eta
 
std::vector< MonitorElement * > ptres_vs_phi
 
std::vector< MonitorElement * > ptres_vs_pt
 
edm::EDGetTokenT
< reco::RecoToSimCollection > 		recoToSimCollection_Token
 
edm::EDGetTokenT
< reco::SimToRecoCollection > 		simToRecoCollection_Token
 
std::vector< MonitorElement * > thetapull_vs_eta
 
std::vector< MonitorElement * > thetapull_vs_phi
 
TrackingParticleSelector tpSelector
 
bool UseAssociators
 
bool useGsf
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...> CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T...> HasAbility
 
typedef
CacheTypes::LuminosityBlockCache 		LuminosityBlockCache
 
typedef
LuminosityBlockContextT
< LuminosityBlockCache,
RunCache, GlobalCache > 		LuminosityBlockContext
 
typedef
CacheTypes::LuminosityBlockSummaryCache 		LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache,
GlobalCache > 		RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase ModuleType
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr
< dqmDetails::NoCache > 		globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
 
static std::shared_ptr
< dqmDetails::NoCache > 		globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
 
static void globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
 
static void globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
 
- Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
- Protected Member Functions inherited from MuonTrackValidatorBase
void BinLogX (TH1 *h)
 
virtual void doProfileX (TH2 *th2, MonitorElement *me)
 
virtual void doProfileX (MonitorElement *th2m, MonitorElement *me)
 
void fillPlotFromVector (MonitorElement *h, std::vector< int > &vec)
 
void fillPlotFromVectors (MonitorElement *h, std::vector< int > &numerator, std::vector< int > &denominator, std::string type)
 
virtual double getEta (double eta)
 
virtual double getPt (double pt)
 
 MuonTrackValidatorBase (const edm::ParameterSet &pset, edm::ConsumesCollector iC)
 Constructor. More...
 
 MuonTrackValidatorBase (const edm::ParameterSet &pset)
 
void setUpVectors ()
 
virtual ~MuonTrackValidatorBase ()
 Destructor. More...
 
- Protected Attributes inherited from MuonTrackValidatorBase
std::vector< std::string > associators
 
edm::InputTag bsSrc
 
edm::EDGetTokenT< reco::BeamSpot > bsSrc_Token
 
int cotThetaRes_nbin
 
double cotThetaRes_rangeMax
 
double cotThetaRes_rangeMin
 
DQMStore * dbe_
 
std::vector< std::vector
< double > > 		dxyintervals
 
int dxyRes_nbin
 
double dxyRes_rangeMax
 
double dxyRes_rangeMin
 
std::vector< std::vector
< double > > 		dzintervals
 
int dzRes_nbin
 
double dzRes_rangeMax
 
double dzRes_rangeMin
 
std::vector< std::vector
< double > > 		etaintervals
 
std::vector< MonitorElement * > h_assoc2dxy
 
std::vector< MonitorElement * > h_assoc2dz
 
std::vector< MonitorElement * > h_assoc2eta
 
std::vector< MonitorElement * > h_assoc2hit
 
std::vector< MonitorElement * > h_assoc2phi
 
std::vector< MonitorElement * > h_assoc2pT
 
std::vector< MonitorElement * > h_assocdxy
 
std::vector< MonitorElement * > h_assocdz
 
std::vector< MonitorElement * > h_assoceta
 
std::vector< MonitorElement * > h_assoceta_Quality05
 
std::vector< MonitorElement * > h_assoceta_Quality075
 
std::vector< MonitorElement * > h_assocFraction
 
std::vector< MonitorElement * > h_assochit
 
std::vector< MonitorElement * > h_assocphi
 
std::vector< MonitorElement * > h_assocphi_Quality05
 
std::vector< MonitorElement * > h_assocphi_Quality075
 
std::vector< MonitorElement * > h_assocpT
 
std::vector< MonitorElement * > h_assocpT_Quality05
 
std::vector< MonitorElement * > h_assocpT_Quality075
 
std::vector< MonitorElement * > h_assocSharedHit
 
std::vector< MonitorElement * > h_assocvertpos
 
std::vector< MonitorElement * > h_assoczpos
 
std::vector< MonitorElement * > h_charge
 
std::vector< MonitorElement * > h_CSChits_eta
 
std::vector< MonitorElement * > h_DThits_eta
 
std::vector< MonitorElement * > h_eta
 
std::vector< MonitorElement * > h_etaSIM
 
std::vector< MonitorElement * > h_fakes
 
std::vector< MonitorElement * > h_hits
 
std::vector< MonitorElement * > h_hits_eta
 
std::vector< MonitorElement * > h_pt
 
std::vector< MonitorElement * > h_ptSIM
 
std::vector< MonitorElement * > h_pullDxy
 
std::vector< MonitorElement * > h_pullDz
 
std::vector< MonitorElement * > h_pullPhi
 
std::vector< MonitorElement * > h_pullQoverp
 
std::vector< MonitorElement * > h_pullTheta
 
std::vector< MonitorElement * > h_PurityVsQuality
 
std::vector< MonitorElement * > h_recodxy
 
std::vector< MonitorElement * > h_recodz
 
std::vector< MonitorElement * > h_recoeta
 
std::vector< MonitorElement * > h_recohit
 
std::vector< MonitorElement * > h_recophi
 
std::vector< MonitorElement * > h_recopT
 
std::vector< MonitorElement * > h_RPChits_eta
 
std::vector< MonitorElement * > h_simuldxy
 
std::vector< MonitorElement * > h_simuldz
 
std::vector< MonitorElement * > h_simuleta
 
std::vector< MonitorElement * > h_simulhit
 
std::vector< MonitorElement * > h_simulphi
 
std::vector< MonitorElement * > h_simulpT
 
std::vector< MonitorElement * > h_simulvertpos
 
std::vector< MonitorElement * > h_simulzpos
 
std::vector< MonitorElement * > h_tracks
 
std::vector< MonitorElement * > h_tracksSIM
 
std::vector< MonitorElement * > h_vertposSIM
 
bool ignoremissingtkcollection_
 
std::vector< edm::InputTag > label
 
edm::InputTag label_tp_effic
 
edm::InputTag label_tp_fake
 
double max
 
double maxDxy
 
double maxDz
 
double maxHit
 
double maxPhi
 
double maxpT
 
double maxVertpos
 
double maxZpos
 
double min
 
double minDxy
 
double minDz
 
double minHit
 
double minPhi
 
double minpT
 
double minVertpos
 
double minZpos
 
std::vector< MonitorElement * > nCSChits_vs_eta
 
std::vector< MonitorElement * > nDThits_vs_eta
 
std::vector< MonitorElement * > nhits_vs_eta
 
int nint
 
int nintDxy
 
int nintDz
 
int nintHit
 
int nintPhi
 
int nintpT
 
int nintVertpos
 
int nintZpos
 
std::vector< MonitorElement * > nrec_vs_nsim
 
std::vector< MonitorElement * > nrecHit_vs_nsimHit_rec2sim
 
std::vector< MonitorElement * > nrecHit_vs_nsimHit_sim2rec
 
std::vector< MonitorElement * > nRPChits_vs_eta
 
std::string out
 
std::string parametersDefiner
 
std::vector< std::vector
< double > > 		phiintervals
 
int phiRes_nbin
 
double phiRes_rangeMax
 
double phiRes_rangeMin
 
std::vector< std::vector
< double > > 		pTintervals
 
int ptRes_nbin
 
double ptRes_rangeMax
 
double ptRes_rangeMin
 
edm::ESHandle< MagneticField > theMF
 
std::vector< std::vector< int > > totASS2_dxy
 
std::vector< std::vector< int > > totASS2_dz
 
std::vector< std::vector< int > > totASS2_hit
 
std::vector< std::vector< int > > totASS2_phi
 
std::vector< std::vector< int > > totASS2eta
 
std::vector< std::vector< int > > totASS2pT
 
std::vector< std::vector< int > > totASS_dxy
 
std::vector< std::vector< int > > totASS_dz
 
std::vector< std::vector< int > > totASS_hit
 
std::vector< std::vector< int > > totASS_phi
 
std::vector< std::vector< int > > totASS_phi_Quality05
 
std::vector< std::vector< int > > totASS_phi_Quality075
 
std::vector< std::vector< int > > totASS_vertpos
 
std::vector< std::vector< int > > totASS_zpos
 
std::vector< std::vector< int > > totASSeta
 
std::vector< std::vector< int > > totASSeta_Quality05
 
std::vector< std::vector< int > > totASSeta_Quality075
 
std::vector< std::vector< int > > totASSpT
 
std::vector< std::vector< int > > totASSpT_Quality05
 
std::vector< std::vector< int > > totASSpT_Quality075
 
std::vector< std::vector< int > > totREC_dxy
 
std::vector< std::vector< int > > totREC_dz
 
std::vector< std::vector< int > > totREC_hit
 
std::vector< std::vector< int > > totREC_phi
 
std::vector< std::vector< int > > totRECeta
 
std::vector< std::vector< int > > totRECpT
 
std::vector< std::vector< int > > totSIM_dxy
 
std::vector< std::vector< int > > totSIM_dz
 
std::vector< std::vector< int > > totSIM_hit
 
std::vector< std::vector< int > > totSIM_phi
 
std::vector< std::vector< int > > totSIM_vertpos
 
std::vector< std::vector< int > > totSIM_zpos
 
std::vector< std::vector< int > > totSIMeta
 
std::vector< std::vector< int > > totSIMpT
 
edm::EDGetTokenT
< TrackingParticleCollection > 		tp_effic_Token
 
edm::EDGetTokenT
< TrackingParticleCollection > 		tp_fake_Token
 
std::vector< edm::EDGetTokenT
< edm::View< reco::Track > > > 		track_Collection_Token
 
bool useFabs
 
bool useInvPt
 
bool useLogPt
 
bool usemuon
 
bool usetracker
 
std::vector< std::vector
< double > > 		vertposintervals
 
std::vector< std::vector
< double > > 		zposintervals
 

Detailed Description

Class that produces histograms to validate Muon Track Reconstruction performances

Definition at line 18 of file MuonTrackValidator.h.

Constructor & Destructor Documentation

MuonTrackValidator::MuonTrackValidator ( const edm::ParameterSet &  pset)
inline

Constructor.

Definition at line 21 of file MuonTrackValidator.h.

References _simHitTpMapTag, associatormap, MuonTrackValidatorBase::associators, BiDirectional_RecoToSim_association, MuonTrackValidatorBase::bsSrc, MuonTrackValidatorBase::bsSrc_Token, edm::EDConsumerBase::consumes(), cosmictpSelector, dirName_, edm::ParameterSet::dump(), edm::ParameterSet::getParameter(), HLT_25ns14e33_v1_cff::InputTag, MuonTrackValidatorBase::label, MuonTrackValidatorBase::label_tp_effic, MuonTrackValidatorBase::label_tp_fake, MABH, maxPhi, minPhi, nintPhi, recoToSimCollection_Token, simToRecoCollection_Token, AlCaHLTBitMon_QueryRunRegistry::string, MuonTrackValidatorBase::tp_effic_Token, MuonTrackValidatorBase::tp_fake_Token, tpSelector, MuonTrackValidatorBase::track_Collection_Token, UseAssociators, useGsf, MuonTrackValidatorBase::usemuon, and MuonTrackValidatorBase::usetracker.

21  :MuonTrackValidatorBase(pset){
22  dirName_ = pset.getParameter<std::string>("dirName");
23  associatormap = pset.getParameter< edm::InputTag >("associatormap");
24  UseAssociators = pset.getParameter< bool >("UseAssociators");
25  tpSelector = TrackingParticleSelector(pset.getParameter<double>("ptMinTP"),
26  pset.getParameter<double>("minRapidityTP"),
27  pset.getParameter<double>("maxRapidityTP"),
28  pset.getParameter<double>("tipTP"),
29  pset.getParameter<double>("lipTP"),
30  pset.getParameter<int>("minHitTP"),
31  pset.getParameter<bool>("signalOnlyTP"),
32  pset.getParameter<bool>("chargedOnlyTP"),
33  pset.getParameter<bool>("stableOnlyTP"),
34  pset.getParameter<std::vector<int> >("pdgIdTP"));
35  cosmictpSelector = CosmicTrackingParticleSelector(pset.getParameter<double>("ptMinTP"),
36  pset.getParameter<double>("minRapidityTP"),
37  pset.getParameter<double>("maxRapidityTP"),
38  pset.getParameter<double>("tipTP"),
39  pset.getParameter<double>("lipTP"),
40  pset.getParameter<int>("minHitTP"),
41  pset.getParameter<bool>("chargedOnlyTP"),
42  pset.getParameter<std::vector<int> >("pdgIdTP"));
43 
44  minPhi = pset.getParameter<double>("minPhi");
45  maxPhi = pset.getParameter<double>("maxPhi");
46  nintPhi = pset.getParameter<int>("nintPhi");
47  useGsf = pset.getParameter<bool>("useGsf");
48  BiDirectional_RecoToSim_association = pset.getParameter<bool>("BiDirectional_RecoToSim_association");
49 
50  // dump cfg parameters
51  edm::LogVerbatim("MuonTrackValidator") << "constructing MuonTrackValidator: " << pset.dump();
52 
53  // Declare consumes (also for the base class)
54  bsSrc_Token = consumes<reco::BeamSpot>(bsSrc);
55  tp_effic_Token = consumes<TrackingParticleCollection>(label_tp_effic);
56  tp_fake_Token = consumes<TrackingParticleCollection>(label_tp_fake);
57  for (unsigned int www=0;www<label.size();www++){
58  track_Collection_Token.push_back(consumes<edm::View<reco::Track> >(label[www]));
59  }
60  simToRecoCollection_Token = consumes<reco::SimToRecoCollection>(associatormap);
61  recoToSimCollection_Token = consumes<reco::RecoToSimCollection>(associatormap);
62 
63  _simHitTpMapTag = mayConsume<SimHitTPAssociationProducer::SimHitTPAssociationList>(pset.getParameter<edm::InputTag>("simHitTpMapTag"));
64 
65  MABH = false;
66  if (!UseAssociators) {
67  // flag MuonAssociatorByHits
68  if (associators[0] == "MuonAssociationByHits") MABH = true;
69  // reset string associators to the map label
70  associators.clear();
71  associators.push_back(associatormap.label());
72  edm::LogVerbatim("MuonTrackValidator") << "--> associators reset to: " <<associators[0];
73  } else {
74  for (auto const& associator :associators) {
75  consumes<reco::TrackToTrackingParticleAssociator>(edm::InputTag(associator));
76  }
77  }
78 
79  // inform on which SimHits will be counted
80  if (usetracker) edm::LogVerbatim("MuonTrackValidator")
81  <<"\n usetracker = TRUE : Tracker SimHits WILL be counted";
82  else edm::LogVerbatim("MuonTrackValidator")
83  <<"\n usetracker = FALSE : Tracker SimHits WILL NOT be counted";
84  if (usemuon) edm::LogVerbatim("MuonTrackValidator")
85  <<" usemuon = TRUE : Muon SimHits WILL be counted";
86  else edm::LogVerbatim("MuonTrackValidator")
87  <<" usemuon = FALSE : Muon SimHits WILL NOT be counted"<<std::endl;
88 
89  // loop over the reco::Track collections to validate: check for inconsistent input settings
90  for (unsigned int www=0;www<label.size();www++) {
91  std::string recoTracksLabel = label[www].label();
92  std::string recoTracksInstance = label[www].instance();
93 
94  // tracks with hits only on tracker
95  if (recoTracksLabel=="generalTracks" ||
96  (recoTracksLabel.find("cutsRecoTracks") != std::string::npos) ||
97  recoTracksLabel=="ctfWithMaterialTracksP5LHCNavigation" ||
98  recoTracksLabel=="hltL3TkTracksFromL2" ||
99  (recoTracksLabel=="hltL3Muons" && recoTracksInstance=="L2Seeded"))
100  {
101  if (usemuon) {
102  edm::LogWarning("MuonTrackValidator")
103  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
104  <<"\n with usemuon == true"<<"\n ---> please change to usemuon == false ";
105  }
106  if (!usetracker) {
107  edm::LogWarning("MuonTrackValidator")
108  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
109  <<"\n with usetracker == false"<<"\n ---> please change to usetracker == true ";
110  }
111  }
112 
113  // tracks with hits only on muon detectors
114  else if (recoTracksLabel=="standAloneMuons" ||
115  recoTracksLabel=="standAloneSETMuons" ||
116  recoTracksLabel=="cosmicMuons" ||
117  recoTracksLabel=="hltL2Muons")
118  {
119  if (usetracker) {
120  edm::LogWarning("MuonTrackValidator")
121  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
122  <<"\n with usetracker == true"<<"\n ---> please change to usetracker == false ";
123  }
124  if (!usemuon) {
125  edm::LogWarning("MuonTrackValidator")
126  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
127  <<"\n with usemuon == false"<<"\n ---> please change to usemuon == true ";
128  }
129  }
130 
131  } // for (unsigned int www=0;www<label.size();www++)
132  }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
MuonTrackValidatorBase::track_Collection_Token
std::vector< edm::EDGetTokenT< edm::View< reco::Track > > > track_Collection_Token
Definition: MuonTrackValidatorBase.h:358
edm::ParameterSet::dump
std::string dump(unsigned int indent=0) const
Definition: ParameterSet.cc:922
HLT_25ns14e33_v1_cff.InputTag
tuple InputTag
Definition: HLT_25ns14e33_v1_cff.py:45333
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:98
MuonTrackValidator::_simHitTpMapTag
edm::EDGetTokenT< SimHitTPAssociationProducer::SimHitTPAssociationList > _simHitTpMapTag
Definition: MuonTrackValidator.h:160
reco::modules::TrackingParticleSelector
SingleObjectSelector< TrackingParticleCollection,::TrackingParticleSelector > TrackingParticleSelector
Definition: TrackingParticleSelector.cc:18
MuonTrackValidator::cosmictpSelector
CosmicTrackingParticleSelector cosmictpSelector
Definition: MuonTrackValidator.h:169
MuonTrackValidatorBase::tp_fake_Token
edm::EDGetTokenT< TrackingParticleCollection > tp_fake_Token
Definition: MuonTrackValidatorBase.h:361
MuonTrackValidatorBase::label
std::vector< edm::InputTag > label
Definition: MuonTrackValidatorBase.h:349
MuonTrackValidatorBase::MuonTrackValidatorBase
MuonTrackValidatorBase(const edm::ParameterSet &pset, edm::ConsumesCollector iC)
Constructor.
Definition: MuonTrackValidatorBase.h:41
reco::CosmicTrackingParticleSelector
ObjectSelector< CosmicTrackingParticleSelector > CosmicTrackingParticleSelector
Definition: CosmicTrackingParticleSelector.cc:16
MuonTrackValidator::recoToSimCollection_Token
edm::EDGetTokenT< reco::RecoToSimCollection > recoToSimCollection_Token
Definition: MuonTrackValidator.h:159
MuonTrackValidator::simToRecoCollection_Token
edm::EDGetTokenT< reco::SimToRecoCollection > simToRecoCollection_Token
Definition: MuonTrackValidator.h:158
MuonTrackValidator::tpSelector
TrackingParticleSelector tpSelector
Definition: MuonTrackValidator.h:168
MuonTrackValidator::associatormap
edm::InputTag associatormap
Definition: MuonTrackValidator.h:157
MuonTrackValidatorBase::associators
std::vector< std::string > associators
Definition: MuonTrackValidatorBase.h:355
MuonTrackValidatorBase::tp_effic_Token
edm::EDGetTokenT< TrackingParticleCollection > tp_effic_Token
Definition: MuonTrackValidatorBase.h:360
MuonTrackValidatorBase::bsSrc_Token
edm::EDGetTokenT< reco::BeamSpot > bsSrc_Token
Definition: MuonTrackValidatorBase.h:359
virtual MuonTrackValidator::~MuonTrackValidator ( )
inlinevirtual

Destructor.

Definition at line 135 of file MuonTrackValidator.h.

135 { }

Member Function Documentation

void MuonTrackValidator::analyze ( const edm::Event &  event,
const edm::EventSetup &  setup 
)
virtual

Method called before the event loop.

Method called once per event

Implements edm::stream::EDAnalyzerBase.

Definition at line 227 of file MuonTrackValidator.cc.

References funct::abs(), asciidump::at, reco::TrackBase::charge(), reco::TrackBase::chi2(), funct::cos(), reco::TrackBase::dxy(), reco::TrackBase::dxyError(), reco::TrackBase::dz(), reco::TrackBase::dzError(), edm::AssociationMap< Tag >::end(), reco::TrackBase::eta(), f, edm::AssociationMap< Tag >::find(), reco::Track::found(), HLT_25ns14e33_v1_cff::fraction, edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), reco::TrackBase::hitPattern(), i, edm::RefToBase< T >::key(), edm::Ref< C, T, F >::key(), diffTwoXMLs::label, LogTrace, M_PI, min(), TrackingParticle::momentum(), reco::TrackBase::momentum(), reco::TrackBase::ndof(), reco::TrackBase::normalizedChi2(), reco::TrackBase::numberOfLostHits(), reco::TrackBase::numberOfValidHits(), reco::HitPattern::numberOfValidMuonCSCHits(), reco::HitPattern::numberOfValidMuonDTHits(), reco::HitPattern::numberOfValidMuonRPCHits(), reco::TrackBase::parameters(), edm::AssociationMap< Tag >::post_insert(), edm::Handle< T >::product(), reco::TrackBase::pt(), HLT_25ns14e33_v1_cff::quality, dt_dqm_sourceclient_common_cff::reco, funct::sin(), mathSSE::sqrt(), reco::TrackBase::theta(), tmp, HLT_50ns_5e33_v1_cff::trackCollection, TrackingParticle::vertex(), and w.

227  {
228  using namespace reco;
229 
230  edm::LogInfo("MuonTrackValidator") << "\n====================================================" << "\n"
231  << "Analyzing new event" << "\n"
232  << "====================================================\n" << "\n";
233 
234  edm::ESHandle<ParametersDefinerForTP> Lhc_parametersDefinerTP;
235  std::unique_ptr<ParametersDefinerForTP> Cosmic_parametersDefinerTP;
236 
237  if(parametersDefiner=="LhcParametersDefinerForTP") {
238  setup.get<TrackAssociatorRecord>().get(parametersDefiner, Lhc_parametersDefinerTP);
239  }
240  else if(parametersDefiner=="CosmicParametersDefinerForTP") {
241  edm::ESHandle<CosmicParametersDefinerForTP> _Cosmic_parametersDefinerTP;
242  setup.get<TrackAssociatorRecord>().get(parametersDefiner, _Cosmic_parametersDefinerTP);
243 
244  //Since we modify the object, we must clone it
245  Cosmic_parametersDefinerTP = _Cosmic_parametersDefinerTP->clone();
246 
247  edm::Handle<SimHitTPAssociationProducer::SimHitTPAssociationList> simHitsTPAssoc;
248  //warning: make sure the TP collection used in the map is the same used here
249  event.getByToken(_simHitTpMapTag,simHitsTPAssoc);
250  Cosmic_parametersDefinerTP->initEvent(simHitsTPAssoc);
251  cosmictpSelector.initEvent(simHitsTPAssoc);
252  }
253  else {
254  edm::LogError("MuonTrackValidator")
255  << "Unexpected label: parametersDefiner = "<< parametersDefiner.c_str() << "\n";
256  }
257 
258  edm::Handle<TrackingParticleCollection> TPCollectionHeff ;
259  event.getByToken(tp_effic_Token,TPCollectionHeff);
260  TrackingParticleCollection const & tPCeff = *(TPCollectionHeff.product());
261 
262  edm::Handle<TrackingParticleCollection> TPCollectionHfake ;
263  event.getByToken(tp_fake_Token,TPCollectionHfake);
264 
265  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
266  event.getByToken(bsSrc_Token,recoBeamSpotHandle);
267  reco::BeamSpot bs = *recoBeamSpotHandle;
268 
269  std::vector<const reco::TrackToTrackingParticleAssociator*> associator;
270  if (UseAssociators) {
271  edm::Handle<reco::TrackToTrackingParticleAssociator> theAssociator;
272  for (unsigned int w=0;w<associators.size();w++) {
273  event.getByLabel(associators[w],theAssociator);
274  associator.push_back( theAssociator.product() );
275  }
276  }
277 
278  int w=0;
279  for (unsigned int ww=0;ww<associators.size();ww++){
280  for (unsigned int www=0;www<label.size();www++){
281  //
282  //get collections from the event
283  //
284  edm::Handle<edm::View<Track> > trackCollection;
285 
286  reco::RecoToSimCollection recSimColl;
287  reco::SimToRecoCollection simRecColl;
288  unsigned int trackCollectionSize = 0;
289 
290  if(!event.getByToken(track_Collection_Token[www], trackCollection)&&ignoremissingtkcollection_) {
291 
292  recSimColl.post_insert();
293  simRecColl.post_insert();
294 
295  }
296 
297  else {
298 
299  trackCollectionSize = trackCollection->size();
300 
301  //associate tracks
302  if(UseAssociators){
303  edm::LogVerbatim("MuonTrackValidator") << "Analyzing "
304  << label[www].process()<<":"
305  << label[www].label()<<":"
306  << label[www].instance()<<" with "
307  << associators[ww].c_str() <<"\n";
308 
309  LogTrace("MuonTrackValidator") << "Calling associateRecoToSim method" << "\n";
310  recSimColl=associator[ww]->associateRecoToSim(trackCollection,
311  TPCollectionHfake);
312  LogTrace("MuonTrackValidator") << "Calling associateSimToReco method" << "\n";
313  simRecColl=associator[ww]->associateSimToReco(trackCollection,
314  TPCollectionHeff);
315  }
316  else{
317  edm::LogVerbatim("MuonTrackValidator") << "Analyzing "
318  << label[www].process()<<":"
319  << label[www].label()<<":"
320  << label[www].instance()<<" with "
321  << associatormap.process()<<":"
322  << associatormap.label()<<":"
323  << associatormap.instance()<<"\n";
324 
325  Handle<reco::SimToRecoCollection > simtorecoCollectionH;
326  event.getByToken(simToRecoCollection_Token,simtorecoCollectionH);
327  simRecColl= *(simtorecoCollectionH.product());
328 
329  Handle<reco::RecoToSimCollection > recotosimCollectionH;
330  event.getByToken(recoToSimCollection_Token,recotosimCollectionH);
331  recSimColl= *(recotosimCollectionH.product());
332  }
333 
334  }
335 
336  //
337  //fill simulation histograms
338  //compute number of tracks per eta interval
339  //
340  edm::LogVerbatim("MuonTrackValidator") << "\n# of TrackingParticles: " << tPCeff.size() << "\n";
341  int ats = 0;
342  int st = 0;
343  for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
344  bool TP_is_matched = false;
345  double quality = 0.;
346  bool Quality05 = false;
347  bool Quality075 = false;
348 
349  TrackingParticleRef tpr(TPCollectionHeff, i);
350  TrackingParticle* tp = const_cast<TrackingParticle*>(tpr.get());
351 
352  TrackingParticle::Vector momentumTP;
353  TrackingParticle::Point vertexTP;
354  double dxySim = 0;
355  double dzSim = 0;
356 
357  //If the TrackingParticle is collision-like, get the momentum and vertex at production state
358  //and the impact parameters w.r.t. PCA
359  if(parametersDefiner=="LhcParametersDefinerForTP")
360  {
361  LogTrace("MuonTrackValidator") <<"TrackingParticle "<< i;
362  if(! tpSelector(*tp)) continue;
363  momentumTP = tp->momentum();
364  vertexTP = tp->vertex();
365  TrackingParticle::Vector momentum = Lhc_parametersDefinerTP->momentum(event,setup,tpr);
366  TrackingParticle::Point vertex = Lhc_parametersDefinerTP->vertex(event,setup,tpr);
367  dxySim = (-vertex.x()*sin(momentum.phi())+vertex.y()*cos(momentum.phi()));
368  dzSim = vertex.z() - (vertex.x()*momentum.x()+vertex.y()*momentum.y()) /
369  sqrt(momentum.perp2()) * momentum.z()/sqrt(momentum.perp2());
370  }
371  //for cosmics get the momentum and vertex at PCA
372  else if(parametersDefiner=="CosmicParametersDefinerForTP")
373  {
374  edm::LogVerbatim("MuonTrackValidator") <<"TrackingParticle "<< i;
375  if(! cosmictpSelector(tpr,&bs,event,setup)) continue;
376  momentumTP = Cosmic_parametersDefinerTP->momentum(event,setup,tpr);
377  vertexTP = Cosmic_parametersDefinerTP->vertex(event,setup,tpr);
378  dxySim = (-vertexTP.x()*sin(momentumTP.phi())+vertexTP.y()*cos(momentumTP.phi()));
379  dzSim = vertexTP.z() - (vertexTP.x()*momentumTP.x()+vertexTP.y()*momentumTP.y()) /
380  sqrt(momentumTP.perp2()) * momentumTP.z()/sqrt(momentumTP.perp2());
381  }
382  edm::LogVerbatim("MuonTrackValidator") <<"--------------------Selected TrackingParticle #"<<tpr.key();
383  edm::LogVerbatim("MuonTrackValidator") <<"momentumTP: pt = "<<sqrt(momentumTP.perp2())<<", pz = "<<momentumTP.z()
384  <<", \t vertexTP: radius = "<<sqrt(vertexTP.perp2())<< ", z = "<<vertexTP.z() <<"\n";
385  st++;
386 
387  h_ptSIM[w]->Fill(sqrt(momentumTP.perp2()));
388  h_etaSIM[w]->Fill(momentumTP.eta());
389  h_vertposSIM[w]->Fill(sqrt(vertexTP.perp2()));
390 
391  std::vector<std::pair<RefToBase<Track>, double> > rt;
392  if(simRecColl.find(tpr) != simRecColl.end()){
393  rt = (std::vector<std::pair<RefToBase<Track>, double> >) simRecColl[tpr];
394  if (rt.size()!=0) {
395  RefToBase<Track> assoc_recoTrack = rt.begin()->first;
396  edm::LogVerbatim("MuonTrackValidator")<<"-----------------------------associated Track #"<<assoc_recoTrack.key();
397  TP_is_matched = true;
398  ats++;
399  quality = rt.begin()->second;
400  edm::LogVerbatim("MuonTrackValidator") << "TrackingParticle #" <<tpr.key()
401  << " with pt=" << sqrt(momentumTP.perp2())
402  << " associated with quality:" << quality <<"\n";
403  if (MABH) {
404  if (quality > 0.75) {
405  Quality075 = true;
406  Quality05 = true;
407  }
408  else if (quality > 0.5) {
409  Quality05 = true;
410  }
411  }
412  }
413  }else{
414  edm::LogVerbatim("MuonTrackValidator")
415  << "TrackingParticle #" << tpr.key()
416  << " with pt,eta,phi: "
417  << sqrt(momentumTP.perp2()) << " , "
418  << momentumTP.eta() << " , "
419  << momentumTP.phi() << " , "
420  << " NOT associated to any reco::Track" << "\n";
421  }
422 
423  for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
424  if (getEta(momentumTP.eta())>etaintervals[w][f]&&
425  getEta(momentumTP.eta())<etaintervals[w][f+1]) {
426  totSIMeta[w][f]++;
427  if (TP_is_matched) {
428  totASSeta[w][f]++;
429 
430  if (MABH) {
431  if (Quality075) {
432  totASSeta_Quality075[w][f]++;
433  totASSeta_Quality05[w][f]++;
434  }
435  else if (Quality05) {
436  totASSeta_Quality05[w][f]++;
437  }
438  }
439  }
440  }
441  } // END for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
442 
443  for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
444  if (momentumTP.phi() > phiintervals[w][f]&&
445  momentumTP.phi() <phiintervals[w][f+1]) {
446  totSIM_phi[w][f]++;
447  if (TP_is_matched) {
448  totASS_phi[w][f]++;
449 
450  if (MABH) {
451  if (Quality075) {
452  totASS_phi_Quality075[w][f]++;
453  totASS_phi_Quality05[w][f]++;
454  }
455  else if (Quality05) {
456  totASS_phi_Quality05[w][f]++;
457  }
458  }
459  }
460  }
461  } // END for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
462 
463  for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
464  if (getPt(sqrt(momentumTP.perp2()))>pTintervals[w][f]&&
465  getPt(sqrt(momentumTP.perp2()))<pTintervals[w][f+1]) {
466  totSIMpT[w][f]++;
467  if (TP_is_matched) {
468  totASSpT[w][f]++;
469 
470  if (MABH) {
471  if (Quality075) {
472  totASSpT_Quality075[w][f]++;
473  totASSpT_Quality05[w][f]++;
474  }
475  else if (Quality05) {
476  totASSpT_Quality05[w][f]++;
477  }
478  }
479  }
480  }
481  } // END for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
482 
483  for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
484  if (dxySim>dxyintervals[w][f]&&
485  dxySim<dxyintervals[w][f+1]) {
486  totSIM_dxy[w][f]++;
487  if (TP_is_matched) {
488  totASS_dxy[w][f]++;
489  }
490  }
491  } // END for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
492 
493  for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
494  if (dzSim>dzintervals[w][f]&&
495  dzSim<dzintervals[w][f+1]) {
496  totSIM_dz[w][f]++;
497  if (TP_is_matched) {
498  totASS_dz[w][f]++;
499  }
500  }
501  } // END for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
502 
503  for (unsigned int f=0; f<vertposintervals[w].size()-1; f++){
504  if (sqrt(vertexTP.perp2())>vertposintervals[w][f]&&
505  sqrt(vertexTP.perp2())<vertposintervals[w][f+1]) {
506  totSIM_vertpos[w][f]++;
507  if (TP_is_matched) {
508  totASS_vertpos[w][f]++;
509  }
510  }
511  } // END for (unsigned int f=0; f<vertposintervals[w].size()-1; f++){
512 
513  for (unsigned int f=0; f<zposintervals[w].size()-1; f++){
514  if (vertexTP.z()>zposintervals[w][f]&&
515  vertexTP.z()<zposintervals[w][f+1]) {
516  totSIM_zpos[w][f]++;
517  if (TP_is_matched) {
518  totASS_zpos[w][f]++;
519  }
520  }
521  } // END for (unsigned int f=0; f<zposintervals[w].size()-1; f++){
522 
523  int nSimHits = 0;
524  if (usetracker && usemuon) {
525  nSimHits= tpr.get()->numberOfHits();
526  }
527  else if (!usetracker && usemuon) {
528  nSimHits= tpr.get()->numberOfHits() - tpr.get()->numberOfTrackerHits();
529  }
530  else if (usetracker && !usemuon) {
531  nSimHits=tpr.get()->numberOfTrackerHits();
532  }
533 
534  int tmp = std::min(nSimHits,int(maxHit-1));
535  edm::LogVerbatim("MuonTrackValidator") << "\t N simhits = "<< nSimHits<<"\n";
536 
537  totSIM_hit[w][tmp]++;
538  if (TP_is_matched) totASS_hit[w][tmp]++;
539 
540  if (TP_is_matched)
541  {
542  RefToBase<Track> assoctrack = rt.begin()->first;
543  nrecHit_vs_nsimHit_sim2rec[w]->Fill( assoctrack->numberOfValidHits(),nSimHits);
544  }
545  } // End for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
546  if (st!=0) h_tracksSIM[w]->Fill(st);
547 
548  //
549  //fill reconstructed track histograms
550  //
551  edm::LogVerbatim("MuonTrackValidator") << "\n# of reco::Tracks with "
552  << label[www].process()<<":"
553  << label[www].label()<<":"
554  << label[www].instance()
555  << ": " << trackCollectionSize << "\n";
556  int at = 0;
557  int rT = 0;
558  for(edm::View<Track>::size_type i=0; i<trackCollectionSize; ++i){
559  bool Track_is_matched = false;
560  RefToBase<Track> track(trackCollection, i);
561  rT++;
562 
563  std::vector<std::pair<TrackingParticleRef, double> > tp;
564  TrackingParticleRef tpr;
565 
566  // new logic (bidirectional)
567  if (BiDirectional_RecoToSim_association) {
568  edm::LogVerbatim("MuonTrackValidator")<<"----------------------------------------Track #"<< track.key();
569 
570  if(recSimColl.find(track) != recSimColl.end()) {
571  tp = recSimColl[track];
572  if (tp.size() != 0) {
573  tpr = tp.begin()->first;
574  // RtS and StR must associate the same pair !
575  if(simRecColl.find(tpr) != simRecColl.end()) {
576  std::vector<std::pair<RefToBase<Track>, double> > track_checkback = simRecColl[tpr];
577  RefToBase<Track> assoc_track_checkback;
578  assoc_track_checkback = track_checkback.begin()->first;
579 
580  if ( assoc_track_checkback.key() == track.key() ) {
581  edm::LogVerbatim("MuonTrackValidator")<<"------------------associated TrackingParticle #"<<tpr.key();
582  Track_is_matched = true;
583  at++;
584  double Purity = tp.begin()->second;
585  double Quality = track_checkback.begin()->second;
586  edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
587  << " associated with quality:" << Purity <<"\n";
588  if (MABH) h_PurityVsQuality[w]->Fill(Quality,Purity);
589  }
590  }
591  }
592  }
593 
594  if (!Track_is_matched)
595  edm::LogVerbatim("MuonTrackValidator")
596  << "reco::Track #" << track.key() << " with pt=" << track->pt() << " NOT associated to any TrackingParticle" << "\n";
597  }
598  // old logic (bugged for collision scenario, still valid for cosmics 2 legs reco)
599  else {
600  if(recSimColl.find(track) != recSimColl.end()){
601  tp = recSimColl[track];
602  if (tp.size()!=0) {
603  Track_is_matched = true;
604  tpr = tp.begin()->first;
605  at++;
606  edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
607  << " associated with quality:" << tp.begin()->second <<"\n";
608  }
609  } else {
610  edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
611  << " NOT associated to any TrackingParticle" << "\n";
612  }
613  }
614 
615  //Compute fake rate vs eta
616  for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
617  if (getEta(track->momentum().eta())>etaintervals[w][f]&&
618  getEta(track->momentum().eta())<etaintervals[w][f+1]) {
619  totRECeta[w][f]++;
620  if (Track_is_matched) {
621  totASS2eta[w][f]++;
622  }
623  }
624  } // End for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
625 
626  for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
627  if (track->momentum().phi()>phiintervals[w][f]&&
628  track->momentum().phi()<phiintervals[w][f+1]) {
629  totREC_phi[w][f]++;
630  if (Track_is_matched) {
631  totASS2_phi[w][f]++;
632  }
633  }
634  } // End for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
635 
636  for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
637  if (getPt(sqrt(track->momentum().perp2()))>pTintervals[w][f]&&
638  getPt(sqrt(track->momentum().perp2()))<pTintervals[w][f+1]) {
639  totRECpT[w][f]++;
640  if (Track_is_matched) {
641  totASS2pT[w][f]++;
642  }
643  }
644  } // End for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
645 
646  for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
647  if (track->dxy(bs.position())>dxyintervals[w][f]&&
648  track->dxy(bs.position())<dxyintervals[w][f+1]) {
649  totREC_dxy[w][f]++;
650  if (Track_is_matched) {
651  totASS2_dxy[w][f]++;
652  }
653  }
654  } // End for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
655 
656  for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
657  if (track->dz(bs.position())>dzintervals[w][f]&&
658  track->dz(bs.position())<dzintervals[w][f+1]) {
659  totREC_dz[w][f]++;
660  if (Track_is_matched) {
661  totASS2_dz[w][f]++;
662  }
663  }
664  } // End for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
665 
666  int tmp = std::min((int)track->found(),int(maxHit-1));
667  totREC_hit[w][tmp]++;
668  if (Track_is_matched) totASS2_hit[w][tmp]++;
669 
670  edm::LogVerbatim("MuonTrackValidator") << "\t N valid rechits = "<< (int)track->found() <<"\n";
671 
672  // Fill other histos
673  TrackingParticle* tpp = const_cast<TrackingParticle*>(tpr.get());
674  // TrackingParticle parameters at point of closest approach to the beamline
675  TrackingParticle::Vector momentumTP;
676  TrackingParticle::Point vertexTP;
677 
678  if (parametersDefiner=="LhcParametersDefinerForTP") {
679  // following reco plots are made only from tracks associated to selected signal TPs
680  if (! (Track_is_matched && tpSelector(*tpp)) ) continue;
681  else {
682  momentumTP = Lhc_parametersDefinerTP->momentum(event,setup,tpr) ;
683  vertexTP = Lhc_parametersDefinerTP->vertex(event,setup,tpr);
684  }
685  }
686  else if (parametersDefiner=="CosmicParametersDefinerForTP") {
687  // following reco plots are made only from tracks associated to selected signal TPs
688  if (! (Track_is_matched && cosmictpSelector(tpr,&bs,event,setup)) ) continue;
689  else {
690  momentumTP = Cosmic_parametersDefinerTP->momentum(event,setup,tpr) ;
691  vertexTP = Cosmic_parametersDefinerTP->vertex(event,setup,tpr);
692  }
693  }
694 
695  if (associators[ww]=="trackAssociatorByChi2"){
696  //association chi2
697  double assocChi2 = -tp.begin()->second;//in association map is stored -chi2
698  h_assochi2[www]->Fill(assocChi2);
699  h_assochi2_prob[www]->Fill(TMath::Prob((assocChi2)*5,5));
700  }
701  else if (associators[ww]=="trackAssociatorByHits"){
702  double fraction = tp.begin()->second;
703  h_assocFraction[www]->Fill(fraction);
704  h_assocSharedHit[www]->Fill(fraction*track->numberOfValidHits());
705  }
706 
707  //nchi2 and hits global distributions
708  h_nchi2[w]->Fill(track->normalizedChi2());
709  h_nchi2_prob[w]->Fill(TMath::Prob(track->chi2(),(int)track->ndof()));
710  h_hits[w]->Fill(track->numberOfValidHits());
711  h_losthits[w]->Fill(track->numberOfLostHits());
712  chi2_vs_nhits[w]->Fill(track->numberOfValidHits(),track->normalizedChi2());
713  h_charge[w]->Fill( track->charge() );
714 
715  double ptSim = sqrt(momentumTP.perp2());
716  double qoverpSim = tpr->charge()/sqrt(momentumTP.x()*momentumTP.x()+momentumTP.y()*momentumTP.y()+momentumTP.z()*momentumTP.z());
717  double thetaSim = momentumTP.theta();
718  double lambdaSim = M_PI/2-momentumTP.theta();
719  double phiSim = momentumTP.phi();
720  double dxySim = (-vertexTP.x()*sin(momentumTP.phi())+vertexTP.y()*cos(momentumTP.phi()));
721  double dzSim = vertexTP.z() - (vertexTP.x()*momentumTP.x()+vertexTP.y()*momentumTP.y()) /
722  sqrt(momentumTP.perp2()) * momentumTP.z()/sqrt(momentumTP.perp2());
723 
724  // removed unused variable, left this in case it has side effects
725  track->parameters();
726 
727  double qoverpRec(0);
728  double qoverpErrorRec(0);
729  double ptRec(0);
730  double ptErrorRec(0);
731  double lambdaRec(0);
732  double lambdaErrorRec(0);
733  double phiRec(0);
734  double phiErrorRec(0);
735 
736  //loop to decide whether to take gsfTrack (utilisation of mode-function) or common track
737  const GsfTrack* gsfTrack(0);
738  if(useGsf){
739  gsfTrack = dynamic_cast<const GsfTrack*>(&(*track));
740  if (gsfTrack==0) edm::LogInfo("MuonTrackValidator") << "Trying to access mode for a non-GsfTrack";
741  }
742 
743  if (gsfTrack) {
744  // get values from mode
745  getRecoMomentum(*gsfTrack, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
746  lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
747  }
748 
749  else {
750  // get values from track (without mode)
751  getRecoMomentum(*track, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
752  lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
753  }
754 
755  double thetaRec = track->theta();
756  double ptError = ptErrorRec;
757  double ptres = ptRec - ptSim;
758  double etares = track->eta()-momentumTP.Eta();
759  double dxyRec = track->dxy(bs.position());
760  double dzRec = track->dz(bs.position());
761  // eta residue; pt, k, theta, phi, dxy, dz pulls
762  double qoverpPull=(qoverpRec-qoverpSim)/qoverpErrorRec;
763  double thetaPull=(lambdaRec-lambdaSim)/lambdaErrorRec;
764  double phiDiff = phiRec - phiSim;
765  if (abs(phiDiff) > M_PI) {
766  if (phiDiff >0.) phiDiff = phiDiff - 2.*M_PI;
767  else phiDiff = phiDiff + 2.*M_PI;
768  }
769  double phiPull=phiDiff/phiErrorRec;
770  double dxyPull=(dxyRec-dxySim)/track->dxyError();
771  double dzPull=(dzRec-dzSim)/track->dzError();
772 
773  double contrib_Qoverp = ((qoverpRec-qoverpSim)/qoverpErrorRec)*((qoverpRec-qoverpSim)/qoverpErrorRec)/5;
774  double contrib_dxy = ((dxyRec-dxySim)/track->dxyError())*((dxyRec-dxySim)/track->dxyError())/5;
775  double contrib_dz = ((dzRec-dzSim)/track->dzError())*((dzRec-dzSim)/track->dzError())/5;
776  double contrib_theta = ((lambdaRec-lambdaSim)/lambdaErrorRec)*((lambdaRec-lambdaSim)/lambdaErrorRec)/5;
777  double contrib_phi = (phiDiff/phiErrorRec)*(phiDiff/phiErrorRec)/5;
778 
779  edm::LogVerbatim("MuonTrackValidator") << "assocChi2=" << tp.begin()->second << "\n"
780  << "" << "\n"
781  << "ptREC=" << ptRec << "\n"
782  << "etaREC=" << track->eta() << "\n"
783  << "qoverpREC=" << qoverpRec << "\n"
784  << "dxyREC=" << dxyRec << "\n"
785  << "dzREC=" << dzRec << "\n"
786  << "thetaREC=" << track->theta() << "\n"
787  << "phiREC=" << phiRec << "\n"
788  << "" << "\n"
789  << "qoverpError()=" << qoverpErrorRec << "\n"
790  << "dxyError()=" << track->dxyError() << "\n"
791  << "dzError()=" << track->dzError() << "\n"
792  << "thetaError()=" << lambdaErrorRec << "\n"
793  << "phiError()=" << phiErrorRec << "\n"
794  << "" << "\n"
795  << "ptSIM=" << ptSim << "\n"
796  << "etaSIM=" << momentumTP.Eta() << "\n"
797  << "qoverpSIM=" << qoverpSim << "\n"
798  << "dxySIM=" << dxySim << "\n"
799  << "dzSIM=" << dzSim << "\n"
800  << "thetaSIM=" << M_PI/2-lambdaSim << "\n"
801  << "phiSIM=" << phiSim << "\n"
802  << "" << "\n"
803  << "contrib_Qoverp=" << contrib_Qoverp << "\n"
804  << "contrib_dxy=" << contrib_dxy << "\n"
805  << "contrib_dz=" << contrib_dz << "\n"
806  << "contrib_theta=" << contrib_theta << "\n"
807  << "contrib_phi=" << contrib_phi << "\n"
808  << "" << "\n"
809  <<"chi2PULL="<<contrib_Qoverp+contrib_dxy+contrib_dz+contrib_theta+contrib_phi<<"\n";
810 
811  h_pullQoverp[w]->Fill(qoverpPull);
812  h_pullTheta[w]->Fill(thetaPull);
813  h_pullPhi[w]->Fill(phiPull);
814  h_pullDxy[w]->Fill(dxyPull);
815  h_pullDz[w]->Fill(dzPull);
816 
817  h_pt[w]->Fill(ptres/ptError);
818  h_eta[w]->Fill(etares);
819  etares_vs_eta[w]->Fill(getEta(track->eta()),etares);
820 
821  //chi2 and #hit vs eta: fill 2D histos
822  chi2_vs_eta[w]->Fill(getEta(track->eta()),track->normalizedChi2());
823  nhits_vs_eta[w]->Fill(getEta(track->eta()),track->numberOfValidHits());
824  nDThits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonDTHits());
825  nCSChits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonCSCHits());
826  nRPChits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonRPCHits());
827 
828  nlosthits_vs_eta[w]->Fill(getEta(track->eta()),track->numberOfLostHits());
829 
830  //resolution of track params: fill 2D histos
831  dxyres_vs_eta[w]->Fill(getEta(track->eta()),dxyRec-dxySim);
832  ptres_vs_eta[w]->Fill(getEta(track->eta()),(ptRec-ptSim)/ptRec);
833  dzres_vs_eta[w]->Fill(getEta(track->eta()),dzRec-dzSim);
834  phires_vs_eta[w]->Fill(getEta(track->eta()),phiDiff);
835  cotThetares_vs_eta[w]->Fill(getEta(track->eta()), cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
836 
837  //same as before but vs pT
838  dxyres_vs_pt[w]->Fill(getPt(ptRec),dxyRec-dxySim);
839  ptres_vs_pt[w]->Fill(getPt(ptRec),(ptRec-ptSim)/ptRec);
840  dzres_vs_pt[w]->Fill(getPt(ptRec),dzRec-dzSim);
841  phires_vs_pt[w]->Fill(getPt(ptRec),phiDiff);
842  cotThetares_vs_pt[w]->Fill(getPt(ptRec), cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
843 
844  //pulls of track params vs eta: fill 2D histos
845  dxypull_vs_eta[w]->Fill(getEta(track->eta()),dxyPull);
846  ptpull_vs_eta[w]->Fill(getEta(track->eta()),ptres/ptError);
847  dzpull_vs_eta[w]->Fill(getEta(track->eta()),dzPull);
848  phipull_vs_eta[w]->Fill(getEta(track->eta()),phiPull);
849  thetapull_vs_eta[w]->Fill(getEta(track->eta()),thetaPull);
850 
851  //plots vs phi
852  nhits_vs_phi[w]->Fill(phiRec,track->numberOfValidHits());
853  chi2_vs_phi[w]->Fill(phiRec,track->normalizedChi2());
854  ptmean_vs_eta_phi[w]->Fill(phiRec,getEta(track->eta()),ptRec);
855  phimean_vs_eta_phi[w]->Fill(phiRec,getEta(track->eta()),phiRec);
856  ptres_vs_phi[w]->Fill(phiRec,(ptRec-ptSim)/ptRec);
857  phires_vs_phi[w]->Fill(phiRec,phiDiff);
858  ptpull_vs_phi[w]->Fill(phiRec,ptres/ptError);
859  phipull_vs_phi[w]->Fill(phiRec,phiPull);
860  thetapull_vs_phi[w]->Fill(phiRec,thetaPull);
861 
862  int nSimHits = 0;
863  if (usetracker && usemuon) {
864  nSimHits= tpr.get()->numberOfHits();
865  }
866  else if (!usetracker && usemuon) {
867  nSimHits= tpr.get()->numberOfHits() - tpr.get()->numberOfTrackerHits();
868  }
869  else if (usetracker && !usemuon) {
870  nSimHits=tpr.get()->numberOfTrackerHits();
871  }
872 
873  nrecHit_vs_nsimHit_rec2sim[w]->Fill(track->numberOfValidHits(), nSimHits);
874 
875  } // End of for(View<Track>::size_type i=0; i<trackCollectionSize; ++i){
876 
877  if (at!=0) h_tracks[w]->Fill(at);
878  h_fakes[w]->Fill(rT-at);
879  edm::LogVerbatim("MuonTrackValidator") << "Total Simulated: " << st << "\n"
880  << "Total Associated (simToReco): " << ats << "\n"
881  << "Total Reconstructed: " << rT << "\n"
882  << "Total Associated (recoToSim): " << at << "\n"
883  << "Total Fakes: " << rT-at << "\n";
884  nrec_vs_nsim[w]->Fill(rT,st);
885  w++;
886  } // End of for (unsigned int www=0;www<label.size();www++){
887  } //END of for (unsigned int ww=0;ww<associators.size();ww++){
888 }
MuonTrackValidatorBase::nrecHit_vs_nsimHit_sim2rec
std::vector< MonitorElement * > nrecHit_vs_nsimHit_sim2rec
Definition: MuonTrackValidatorBase.h:412
edm::View::size_type
unsigned int size_type
Definition: View.h:85
MuonTrackValidatorBase::totASS_zpos
std::vector< std::vector< int > > totASS_zpos
Definition: MuonTrackValidatorBase.h:439
i
int i
Definition: DBlmapReader.cc:9
MuonTrackValidatorBase::h_PurityVsQuality
std::vector< MonitorElement * > h_PurityVsQuality
Definition: MuonTrackValidatorBase.h:442
MuonTrackValidatorBase::totASSeta_Quality05
std::vector< std::vector< int > > totASSeta_Quality05
Definition: MuonTrackValidatorBase.h:443
HLT_25ns14e33_v1_cff.fraction
tuple fraction
Definition: HLT_25ns14e33_v1_cff.py:28703
MuonTrackValidatorBase::totREC_phi
std::vector< std::vector< int > > totREC_phi
Definition: MuonTrackValidatorBase.h:436
MuonTrackValidator::etares_vs_eta
std::vector< MonitorElement * > etares_vs_eta
Definition: MuonTrackValidator.h:180
MuonTrackValidatorBase::nrecHit_vs_nsimHit_rec2sim
std::vector< MonitorElement * > nrecHit_vs_nsimHit_rec2sim
Definition: MuonTrackValidatorBase.h:413
MuonTrackValidator::ptres_vs_phi
std::vector< MonitorElement * > ptres_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::totASS2_hit
std::vector< std::vector< int > > totASS2_hit
Definition: MuonTrackValidatorBase.h:435
w
const double w
Definition: UKUtility.cc:23
MuonTrackValidatorBase::track_Collection_Token
std::vector< edm::EDGetTokenT< edm::View< reco::Track > > > track_Collection_Token
Definition: MuonTrackValidatorBase.h:358
MuonTrackValidatorBase::totRECpT
std::vector< std::vector< int > > totRECpT
Definition: MuonTrackValidatorBase.h:434
TrackingParticleCollection
std::vector< TrackingParticle > TrackingParticleCollection
Definition: HiEvtPlaneFlatProducer.cc:102
MuonTrackValidatorBase::totREC_dz
std::vector< std::vector< int > > totREC_dz
Definition: MuonTrackValidatorBase.h:438
MuonTrackValidatorBase::nrec_vs_nsim
std::vector< MonitorElement * > nrec_vs_nsim
Definition: MuonTrackValidatorBase.h:411
edm::AssociationMap::end
const_iterator end() const
last iterator over the map (read only)
Definition: AssociationMap.h:89
MuonTrackValidatorBase::h_pullDz
std::vector< MonitorElement * > h_pullDz
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidatorBase::h_etaSIM
std::vector< MonitorElement * > h_etaSIM
Definition: MuonTrackValidatorBase.h:392
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
MuonTrackValidatorBase::h_fakes
std::vector< MonitorElement * > h_fakes
Definition: MuonTrackValidatorBase.h:395
MuonTrackValidator::ptres_vs_eta
std::vector< MonitorElement * > ptres_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidator::cotThetares_vs_eta
std::vector< MonitorElement * > cotThetares_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidatorBase::getEta
virtual double getEta(double eta)
Definition: MuonTrackValidatorBase.h:131
MuonTrackValidatorBase::h_eta
std::vector< MonitorElement * > h_eta
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidator::phimean_vs_eta_phi
std::vector< MonitorElement * > phimean_vs_eta_phi
Definition: MuonTrackValidator.h:185
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
MuonTrackValidatorBase::totASS_dz
std::vector< std::vector< int > > totASS_dz
Definition: MuonTrackValidatorBase.h:438
edm::AssociationMap::find
const_iterator find(const key_type &k) const
find element with specified reference key
Definition: AssociationMap.h:91
MuonTrackValidatorBase::etaintervals
std::vector< std::vector< double > > etaintervals
Definition: MuonTrackValidatorBase.h:426
MuonTrackValidatorBase::h_pt
std::vector< MonitorElement * > h_pt
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidator::phires_vs_pt
std::vector< MonitorElement * > phires_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidatorBase::h_pullQoverp
std::vector< MonitorElement * > h_pullQoverp
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidator::h_assochi2_prob
std::vector< MonitorElement * > h_assochi2_prob
Definition: MuonTrackValidator.h:188
MuonTrackValidatorBase::totSIM_zpos
std::vector< std::vector< int > > totSIM_zpos
Definition: MuonTrackValidatorBase.h:439
MuonTrackValidatorBase::totASS2_dz
std::vector< std::vector< int > > totASS2_dz
Definition: MuonTrackValidatorBase.h:438
MuonTrackValidator::cotThetares_vs_pt
std::vector< MonitorElement * > cotThetares_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidator::dxyres_vs_eta
std::vector< MonitorElement * > dxyres_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidatorBase::totASS_phi_Quality05
std::vector< std::vector< int > > totASS_phi_Quality05
Definition: MuonTrackValidatorBase.h:445
MuonTrackValidator::dzres_vs_pt
std::vector< MonitorElement * > dzres_vs_pt
Definition: MuonTrackValidator.h:198
edm::Ref::key
key_type key() const
Accessor for product key.
Definition: Ref.h:266
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
MuonTrackValidatorBase::h_tracks
std::vector< MonitorElement * > h_tracks
Definition: MuonTrackValidatorBase.h:395
MuonTrackValidator::dzres_vs_eta
std::vector< MonitorElement * > dzres_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidator::getRecoMomentum
void getRecoMomentum(const reco::Track &track, double &pt, double &ptError, double &qoverp, double &qoverpError, double &lambda, double &lambdaError, double &phi, double &phiError) const
retrieval of reconstructed momentum components from reco::Track (== mean values for GSF) ...
Definition: MuonTrackValidator.cc:961
MuonTrackValidatorBase::dxyintervals
std::vector< std::vector< double > > dxyintervals
Definition: MuonTrackValidatorBase.h:429
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:19
MuonTrackValidatorBase::nDThits_vs_eta
std::vector< MonitorElement * > nDThits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::totASS_phi
std::vector< std::vector< int > > totASS_phi
Definition: MuonTrackValidatorBase.h:436
MuonTrackValidator::thetapull_vs_phi
std::vector< MonitorElement * > thetapull_vs_phi
Definition: MuonTrackValidator.h:202
MuonTrackValidatorBase::totASS_dxy
std::vector< std::vector< int > > totASS_dxy
Definition: MuonTrackValidatorBase.h:437
MuonTrackValidator::_simHitTpMapTag
edm::EDGetTokenT< SimHitTPAssociationProducer::SimHitTPAssociationList > _simHitTpMapTag
Definition: MuonTrackValidator.h:160
MuonTrackValidator::h_assochi2
std::vector< MonitorElement * > h_assochi2
Definition: MuonTrackValidator.h:188
MuonTrackValidator::chi2_vs_phi
std::vector< MonitorElement * > chi2_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::totSIM_dxy
std::vector< std::vector< int > > totSIM_dxy
Definition: MuonTrackValidatorBase.h:437
MuonTrackValidator::chi2_vs_nhits
std::vector< MonitorElement * > chi2_vs_nhits
Definition: MuonTrackValidator.h:180
TrackingParticle::Point
math::XYZPointD Point
point in the space
Definition: TrackingParticle.h:36
MuonTrackValidatorBase::totASSeta_Quality075
std::vector< std::vector< int > > totASSeta_Quality075
Definition: MuonTrackValidatorBase.h:443
MuonTrackValidatorBase::nhits_vs_eta
std::vector< MonitorElement * > nhits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidator::phires_vs_phi
std::vector< MonitorElement * > phires_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::getPt
virtual double getPt(double pt)
Definition: MuonTrackValidatorBase.h:136
MuonTrackValidator::cosmictpSelector
CosmicTrackingParticleSelector cosmictpSelector
Definition: MuonTrackValidator.h:169
MuonTrackValidator::nhits_vs_phi
std::vector< MonitorElement * > nhits_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::totASS2pT
std::vector< std::vector< int > > totASS2pT
Definition: MuonTrackValidatorBase.h:434
MuonTrackValidatorBase::h_vertposSIM
std::vector< MonitorElement * > h_vertposSIM
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidatorBase::tp_fake_Token
edm::EDGetTokenT< TrackingParticleCollection > tp_fake_Token
Definition: MuonTrackValidatorBase.h:361
MuonTrackValidator::dxyres_vs_pt
std::vector< MonitorElement * > dxyres_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidatorBase::label
std::vector< edm::InputTag > label
Definition: MuonTrackValidatorBase.h:349
MuonTrackValidatorBase::h_assocSharedHit
std::vector< MonitorElement * > h_assocSharedHit
Definition: MuonTrackValidatorBase.h:416
MuonTrackValidator::h_nchi2
std::vector< MonitorElement * > h_nchi2
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::pTintervals
std::vector< std::vector< double > > pTintervals
Definition: MuonTrackValidatorBase.h:427
edm::RefToBase::key
size_t key() const
Definition: RefToBase.h:228
edm::AssociationMap::post_insert
void post_insert()
post insert action
Definition: AssociationMap.h:132
MuonTrackValidatorBase::totASS2eta
std::vector< std::vector< int > > totASS2eta
Definition: MuonTrackValidatorBase.h:433
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
MuonTrackValidatorBase::totSIMpT
std::vector< std::vector< int > > totSIMpT
Definition: MuonTrackValidatorBase.h:434
MuonTrackValidatorBase::totREC_dxy
std::vector< std::vector< int > > totREC_dxy
Definition: MuonTrackValidatorBase.h:437
MuonTrackValidatorBase::h_pullPhi
std::vector< MonitorElement * > h_pullPhi
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidatorBase::totSIM_hit
std::vector< std::vector< int > > totSIM_hit
Definition: MuonTrackValidatorBase.h:435
MuonTrackValidator::h_nchi2_prob
std::vector< MonitorElement * > h_nchi2_prob
Definition: MuonTrackValidator.h:177
MuonTrackValidator::ptpull_vs_eta
std::vector< MonitorElement * > ptpull_vs_eta
Definition: MuonTrackValidator.h:201
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
MuonTrackValidatorBase::totASSpT
std::vector< std::vector< int > > totASSpT
Definition: MuonTrackValidatorBase.h:434
MuonTrackValidator::dxypull_vs_eta
std::vector< MonitorElement * > dxypull_vs_eta
Definition: MuonTrackValidator.h:201
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
MuonTrackValidatorBase::totASS_phi_Quality075
std::vector< std::vector< int > > totASS_phi_Quality075
Definition: MuonTrackValidatorBase.h:445
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
CosmicTrackingParticleSelector::initEvent
void initEvent(edm::Handle< SimHitTPAssociationProducer::SimHitTPAssociationList > simHitsTPAssocToSet) const
Definition: CosmicTrackingParticleSelector.h:81
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:868
MuonTrackValidatorBase::h_pullTheta
std::vector< MonitorElement * > h_pullTheta
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidatorBase::totSIM_vertpos
std::vector< std::vector< int > > totSIM_vertpos
Definition: MuonTrackValidatorBase.h:439
MuonTrackValidatorBase::totASSpT_Quality05
std::vector< std::vector< int > > totASSpT_Quality05
Definition: MuonTrackValidatorBase.h:444
MuonTrackValidatorBase::totREC_hit
std::vector< std::vector< int > > totREC_hit
Definition: MuonTrackValidatorBase.h:435
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:242
min
T min(T a, T b)
Definition: MathUtil.h:58
MuonTrackValidator::phipull_vs_eta
std::vector< MonitorElement * > phipull_vs_eta
Definition: MuonTrackValidator.h:201
MuonTrackValidatorBase::totSIMeta
std::vector< std::vector< int > > totSIMeta
Definition: MuonTrackValidatorBase.h:433
MuonTrackValidatorBase::totASS2_dxy
std::vector< std::vector< int > > totASS2_dxy
Definition: MuonTrackValidatorBase.h:437
MuonTrackValidator::thetapull_vs_eta
std::vector< MonitorElement * > thetapull_vs_eta
Definition: MuonTrackValidator.h:201
MuonTrackValidatorBase::h_charge
std::vector< MonitorElement * > h_charge
Definition: MuonTrackValidatorBase.h:395
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
M_PI
#define M_PI
Definition: readTestVector.cc:25
MuonTrackValidator::ptmean_vs_eta_phi
std::vector< MonitorElement * > ptmean_vs_eta_phi
Definition: MuonTrackValidator.h:185
MuonTrackValidatorBase::totASS2_phi
std::vector< std::vector< int > > totASS2_phi
Definition: MuonTrackValidatorBase.h:436
MuonTrackValidatorBase::phiintervals
std::vector< std::vector< double > > phiintervals
Definition: MuonTrackValidatorBase.h:428
MuonTrackValidator::recoToSimCollection_Token
edm::EDGetTokenT< reco::RecoToSimCollection > recoToSimCollection_Token
Definition: MuonTrackValidator.h:159
MuonTrackValidator::simToRecoCollection_Token
edm::EDGetTokenT< reco::SimToRecoCollection > simToRecoCollection_Token
Definition: MuonTrackValidator.h:158
MuonTrackValidatorBase::totASSpT_Quality075
std::vector< std::vector< int > > totASSpT_Quality075
Definition: MuonTrackValidatorBase.h:444
MuonTrackValidatorBase::h_pullDxy
std::vector< MonitorElement * > h_pullDxy
Definition: MuonTrackValidatorBase.h:403
edm::Handle::product
T const * product() const
Definition: Handle.h:81
MuonTrackValidatorBase::h_tracksSIM
std::vector< MonitorElement * > h_tracksSIM
Definition: MuonTrackValidatorBase.h:392
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
MuonTrackValidator::chi2_vs_eta
std::vector< MonitorElement * > chi2_vs_eta
Definition: MuonTrackValidator.h:191
MuonTrackValidator::tpSelector
TrackingParticleSelector tpSelector
Definition: MuonTrackValidator.h:168
MuonTrackValidatorBase::totSIM_phi
std::vector< std::vector< int > > totSIM_phi
Definition: MuonTrackValidatorBase.h:436
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:42
TrackingParticle::vertex
Point vertex() const
Parent vertex position.
Definition: TrackingParticle.h:120
MuonTrackValidatorBase::h_hits
std::vector< MonitorElement * > h_hits
Definition: MuonTrackValidatorBase.h:395
MuonTrackValidator::h_losthits
std::vector< MonitorElement * > h_losthits
Definition: MuonTrackValidator.h:177
MuonTrackValidator::ptpull_vs_phi
std::vector< MonitorElement * > ptpull_vs_phi
Definition: MuonTrackValidator.h:202
edm::InputTag::process
std::string const & process() const
Definition: InputTag.h:46
tmp
std::vector< std::vector< double > > tmp
Definition: MVATrainer.cc:100
MuonTrackValidator::associatormap
edm::InputTag associatormap
Definition: MuonTrackValidator.h:157
MuonTrackValidatorBase::totASSeta
std::vector< std::vector< int > > totASSeta
Definition: MuonTrackValidatorBase.h:433
MuonTrackValidator::ptres_vs_pt
std::vector< MonitorElement * > ptres_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidatorBase::totRECeta
std::vector< std::vector< int > > totRECeta
Definition: MuonTrackValidatorBase.h:433
MuonTrackValidatorBase::vertposintervals
std::vector< std::vector< double > > vertposintervals
Definition: MuonTrackValidatorBase.h:431
MuonTrackValidatorBase::dzintervals
std::vector< std::vector< double > > dzintervals
Definition: MuonTrackValidatorBase.h:430
MuonTrackValidator::phires_vs_eta
std::vector< MonitorElement * > phires_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidatorBase::h_ptSIM
std::vector< MonitorElement * > h_ptSIM
Definition: MuonTrackValidatorBase.h:392
TrackingParticle::momentum
Vector momentum() const
spatial momentum vector
Definition: TrackingParticle.cc:122
HLT_25ns14e33_v1_cff.quality
tuple quality
Definition: HLT_25ns14e33_v1_cff.py:33032
TrackingParticle
Monte Carlo truth information used for tracking validation.
Definition: TrackingParticle.h:29
MuonTrackValidatorBase::zposintervals
std::vector< std::vector< double > > zposintervals
Definition: MuonTrackValidatorBase.h:432
MuonTrackValidatorBase::h_assocFraction
std::vector< MonitorElement * > h_assocFraction
Definition: MuonTrackValidatorBase.h:416
TrackingParticle::Vector
math::XYZVectorD Vector
point in the space
Definition: TrackingParticle.h:37
HLT_50ns_5e33_v1_cff.trackCollection
tuple trackCollection
Definition: HLT_50ns_5e33_v1_cff.py:45691
MuonTrackValidatorBase::totSIM_dz
std::vector< std::vector< int > > totSIM_dz
Definition: MuonTrackValidatorBase.h:438
MuonTrackValidatorBase::associators
std::vector< std::string > associators
Definition: MuonTrackValidatorBase.h:355
MuonTrackValidator::nlosthits_vs_eta
std::vector< MonitorElement * > nlosthits_vs_eta
Definition: MuonTrackValidator.h:191
MuonTrackValidator::dzpull_vs_eta
std::vector< MonitorElement * > dzpull_vs_eta
Definition: MuonTrackValidator.h:201
MuonTrackValidatorBase::tp_effic_Token
edm::EDGetTokenT< TrackingParticleCollection > tp_effic_Token
Definition: MuonTrackValidatorBase.h:360
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:43
MuonTrackValidatorBase::nRPChits_vs_eta
std::vector< MonitorElement * > nRPChits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::totASS_vertpos
std::vector< std::vector< int > > totASS_vertpos
Definition: MuonTrackValidatorBase.h:439
MuonTrackValidatorBase::totASS_hit
std::vector< std::vector< int > > totASS_hit
Definition: MuonTrackValidatorBase.h:435
MuonTrackValidatorBase::nCSChits_vs_eta
std::vector< MonitorElement * > nCSChits_vs_eta
Definition: MuonTrackValidatorBase.h:419
asciidump.at
list at
Definition: asciidump.py:428
MuonTrackValidator::phipull_vs_phi
std::vector< MonitorElement * > phipull_vs_phi
Definition: MuonTrackValidator.h:202
MuonTrackValidatorBase::bsSrc_Token
edm::EDGetTokenT< reco::BeamSpot > bsSrc_Token
Definition: MuonTrackValidatorBase.h:359
void MuonTrackValidator::bookHistograms ( DQMStore::IBooker &  ibooker,
edm::Run const &  ,
edm::EventSetup const &  setup 
)
virtual

this are needed to calculate efficiency during tha harvesting for the automated validation

Implements DQMEDAnalyzer.

Definition at line 26 of file MuonTrackValidator.cc.

References HLT_25ns14e33_v1_cff::algo, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), DQMStore::IBooker::bookProfile2D(), DQMStore::IBooker::cd(), TrackerOfflineValidation_Dqm_cff::dirName, getTH1F(), getTH2F(), DQMStore::IBooker::goUp(), edm::InputTag::instance(), j, diffTwoXMLs::label, edm::InputTag::label(), bookConverter::max, gather_cfg::maxDxy, min(), HLT_25ns14e33_v1_cff::minHit, nint(), edm::InputTag::process(), python.rootplot.root2matplotlib::replace(), and DQMStore::IBooker::setCurrentFolder().

26  {
27 
28  int j=0;
29  for (unsigned int ww=0;ww<associators.size();ww++){
30  for (unsigned int www=0;www<label.size();www++){
31 
32  ibooker.cd();
33  InputTag algo = label[www];
34  string dirName=dirName_;
35  if (algo.process()!="")
36  dirName+=algo.process()+"_";
37  if(algo.label()!="")
38  dirName+=algo.label();
39  if(algo.instance()!="")
40  dirName+=("_"+algo.instance());
41  if (dirName.find("Tracks")<dirName.length()){
42  dirName.replace(dirName.find("Tracks"),6,"Trks");
43  }
44  if (dirName.find("UpdatedAtVtx")<dirName.length()){
45  dirName.replace(dirName.find("UpdatedAtVtx"),12,"UpdAtVtx");
46  }
47  string assoc= associators[ww];
48  if (assoc.find("tpToTkmuTrackAssociation")<assoc.length()){
49  dirName+="_TkAsso";
50  }
51  std::replace(dirName.begin(), dirName.end(), ':', '_');
52  ibooker.setCurrentFolder(dirName.c_str());
53 
54  setUpVectors();
55 
56  ibooker.goUp();
57  string subDirName = dirName + "/simulation";
58  ibooker.setCurrentFolder(subDirName.c_str());
59  h_ptSIM.push_back( ibooker.book1D("ptSIM", "generated p_{t}", 5500, 0, 110 ) );
60  h_etaSIM.push_back( ibooker.book1D("etaSIM", "generated pseudorapidity", 500, -2.5, 2.5 ) );
61  h_tracksSIM.push_back( ibooker.book1D("tracksSIM","number of simulated tracks",200,-0.5,99.5) );
62  h_vertposSIM.push_back( ibooker.book1D("vertposSIM","Transverse position of sim vertices",100,0.,120.) );
63 
64  ibooker.cd();
65  ibooker.setCurrentFolder(dirName.c_str());
66  h_tracks.push_back( ibooker.book1D("tracks","number of reconstructed tracks",200,-0.5,19.5) );
67  h_fakes.push_back( ibooker.book1D("fakes","number of fake reco tracks",20,-0.5,19.5) );
68  h_charge.push_back( ibooker.book1D("charge","charge",3,-1.5,1.5) );
69  h_hits.push_back( ibooker.book1D("hits", "number of hits per track", nintHit,minHit,maxHit ) );
70  h_losthits.push_back( ibooker.book1D("losthits", "number of lost hits per track", nintHit,minHit,maxHit) );
71  h_nchi2.push_back( ibooker.book1D("chi2", "normalized #chi^{2}", 200, 0, 20 ) );
72  h_nchi2_prob.push_back( ibooker.book1D("chi2_prob", "normalized #chi^{2} probability",100,0,1));
73 
75  h_recoeta.push_back( ibooker.book1D("num_reco_eta","N of reco track vs eta",nint,min,max) );
76  h_assoceta.push_back( ibooker.book1D("num_assoc(simToReco)_eta","N of associated tracks (simToReco) vs eta",nint,min,max) );
77  h_assoc2eta.push_back( ibooker.book1D("num_assoc(recoToSim)_eta","N of associated (recoToSim) tracks vs eta",nint,min,max) );
78  h_simuleta.push_back( ibooker.book1D("num_simul_eta","N of simulated tracks vs eta",nint,min,max) );
79  h_recopT.push_back( ibooker.book1D("num_reco_pT","N of reco track vs pT",nintpT,minpT,maxpT) );
80  h_assocpT.push_back( ibooker.book1D("num_assoc(simToReco)_pT","N of associated tracks (simToReco) vs pT",nintpT,minpT,maxpT) );
81  h_assoc2pT.push_back( ibooker.book1D("num_assoc(recoToSim)_pT","N of associated (recoToSim) tracks vs pT",nintpT,minpT,maxpT) );
82  h_simulpT.push_back( ibooker.book1D("num_simul_pT","N of simulated tracks vs pT",nintpT,minpT,maxpT) );
83  //
84  h_recohit.push_back( ibooker.book1D("num_reco_hit","N of reco track vs hit",nintHit,minHit,maxHit) );
85  h_assochit.push_back( ibooker.book1D("num_assoc(simToReco)_hit","N of associated tracks (simToReco) vs hit",nintHit,minHit,maxHit) );
86  h_assoc2hit.push_back( ibooker.book1D("num_assoc(recoToSim)_hit","N of associated (recoToSim) tracks vs hit",nintHit,minHit,maxHit) );
87  h_simulhit.push_back( ibooker.book1D("num_simul_hit","N of simulated tracks vs hit",nintHit,minHit,maxHit) );
88  //
89  h_recophi.push_back( ibooker.book1D("num_reco_phi","N of reco track vs phi",nintPhi,minPhi,maxPhi) );
90  h_assocphi.push_back( ibooker.book1D("num_assoc(simToReco)_phi","N of associated tracks (simToReco) vs phi",nintPhi,minPhi,maxPhi) );
91  h_assoc2phi.push_back( ibooker.book1D("num_assoc(recoToSim)_phi","N of associated (recoToSim) tracks vs phi",nintPhi,minPhi,maxPhi) );
92  h_simulphi.push_back( ibooker.book1D("num_simul_phi","N of simulated tracks vs phi",nintPhi,minPhi,maxPhi) );
93 
94  h_recodxy.push_back( ibooker.book1D("num_reco_dxy","N of reco track vs dxy",nintDxy,minDxy,maxDxy) );
95  h_assocdxy.push_back( ibooker.book1D("num_assoc(simToReco)_dxy","N of associated tracks (simToReco) vs dxy",nintDxy,minDxy,maxDxy) );
96  h_assoc2dxy.push_back( ibooker.book1D("num_assoc(recoToSim)_dxy","N of associated (recoToSim) tracks vs dxy",nintDxy,minDxy,maxDxy) );
97  h_simuldxy.push_back( ibooker.book1D("num_simul_dxy","N of simulated tracks vs dxy",nintDxy,minDxy,maxDxy) );
98 
99  h_recodz.push_back( ibooker.book1D("num_reco_dz","N of reco track vs dz",nintDz,minDz,maxDz) );
100  h_assocdz.push_back( ibooker.book1D("num_assoc(simToReco)_dz","N of associated tracks (simToReco) vs dz",nintDz,minDz,maxDz) );
101  h_assoc2dz.push_back( ibooker.book1D("num_assoc(recoToSim)_dz","N of associated (recoToSim) tracks vs dz",nintDz,minDz,maxDz) );
102  h_simuldz.push_back( ibooker.book1D("num_simul_dz","N of simulated tracks vs dz",nintDz,minDz,maxDz) );
103 
104  h_assocvertpos.push_back( ibooker.book1D("num_assoc(simToReco)_vertpos","N of associated tracks (simToReco) vs transverse vert position",nintVertpos,minVertpos,maxVertpos) );
105  h_simulvertpos.push_back( ibooker.book1D("num_simul_vertpos","N of simulated tracks vs transverse vert position",nintVertpos,minVertpos,maxVertpos) );
106 
107  h_assoczpos.push_back( ibooker.book1D("num_assoc(simToReco)_zpos","N of associated tracks (simToReco) vs z vert position",nintZpos,minZpos,maxZpos) );
108  h_simulzpos.push_back( ibooker.book1D("num_simul_zpos","N of simulated tracks vs z vert position",nintZpos,minZpos,maxZpos) );
109 
111 
112  h_eta.push_back( ibooker.book1D("eta", "pseudorapidity residue", 1000, -0.1, 0.1 ) );
113  h_pt.push_back( ibooker.book1D("pullPt", "pull of p_{t}", 100, -10, 10 ) );
114  h_pullTheta.push_back( ibooker.book1D("pullTheta","pull of #theta parameter",250,-25,25) );
115  h_pullPhi.push_back( ibooker.book1D("pullPhi","pull of #phi parameter",250,-25,25) );
116  h_pullDxy.push_back( ibooker.book1D("pullDxy","pull of dxy parameter",250,-25,25) );
117  h_pullDz.push_back( ibooker.book1D("pullDz","pull of dz parameter",250,-25,25) );
118  h_pullQoverp.push_back( ibooker.book1D("pullQoverp","pull of qoverp parameter",250,-25,25) );
119 
120  if (associators[ww]=="trackAssociatorByChi2"){
121  h_assochi2.push_back( ibooker.book1D("assocChi2","track association #chi^{2}",1000000,0,100000) );
122  h_assochi2_prob.push_back(ibooker.book1D("assocChi2_prob","probability of association #chi^{2}",100,0,1));
123  } else if (associators[ww]=="trackAssociatorByHits"){
124  h_assocFraction.push_back( ibooker.book1D("assocFraction","fraction of shared hits",200,0,2) );
125  h_assocSharedHit.push_back(ibooker.book1D("assocSharedHit","number of shared hits",20,0,20));
126  }
127 
128  chi2_vs_nhits.push_back( ibooker.book2D("chi2_vs_nhits","#chi^{2} vs nhits",25,0,25,100,0,10) );
129  h_chi2meanhitsh.push_back( ibooker.bookProfile("chi2mean_vs_nhits","mean #chi^{2} vs nhits",25,0,25,100,0,10) );
130 
131  etares_vs_eta.push_back( ibooker.book2D("etares_vs_eta","etaresidue vs eta",nint,min,max,200,-0.1,0.1) );
132  nrec_vs_nsim.push_back( ibooker.book2D("nrec_vs_nsim","nrec vs nsim",20,-0.5,19.5,20,-0.5,19.5) );
133 
134  chi2_vs_eta.push_back( ibooker.book2D("chi2_vs_eta","chi2_vs_eta",nint,min,max, 200, 0, 20 ));
135  h_chi2meanh.push_back( ibooker.bookProfile("chi2mean","mean #chi^{2} vs #eta",nint,min,max, 200, 0, 20) );
136  chi2_vs_phi.push_back( ibooker.book2D("chi2_vs_phi","#chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 200, 0, 20 ) );
137  h_chi2mean_vs_phi.push_back( ibooker.bookProfile("chi2mean_vs_phi","mean of #chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 200, 0, 20) );
138 
139  nhits_vs_eta.push_back( ibooker.book2D("nhits_vs_eta","nhits vs eta",nint,min,max,nintHit,minHit,maxHit) );
140  nDThits_vs_eta.push_back( ibooker.book2D("nDThits_vs_eta","# DT hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
141  nCSChits_vs_eta.push_back( ibooker.book2D("nCSChits_vs_eta","# CSC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
142  nRPChits_vs_eta.push_back( ibooker.book2D("nRPChits_vs_eta","# RPC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
143 
144  h_DThits_eta.push_back( ibooker.bookProfile("DThits_eta","mean # DT hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
145  h_CSChits_eta.push_back( ibooker.bookProfile("CSChits_eta","mean # CSC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
146  h_RPChits_eta.push_back( ibooker.bookProfile("RPChits_eta","mean # RPC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
147  h_hits_eta.push_back( ibooker.bookProfile("hits_eta","mean #hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
148  nhits_vs_phi.push_back( ibooker.book2D("nhits_vs_phi","#hits vs #phi",nintPhi,minPhi,maxPhi,nintHit,minHit,maxHit) );
149  h_hits_phi.push_back( ibooker.bookProfile("hits_phi","mean #hits vs #phi",nintPhi,minPhi,maxPhi, nintHit,minHit,maxHit) );
150 
151  nlosthits_vs_eta.push_back( ibooker.book2D("nlosthits_vs_eta","nlosthits vs eta",nint,min,max,nintHit,minHit,maxHit) );
152  h_losthits_eta.push_back( ibooker.bookProfile("losthits_eta","losthits_eta",nint,min,max,nintHit,minHit,maxHit) );
153 
154  ptres_vs_eta.push_back(ibooker.book2D("ptres_vs_eta","ptres_vs_eta",nint,min,max, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
155  ptres_vs_phi.push_back( ibooker.book2D("ptres_vs_phi","p_{t} res vs #phi",nintPhi,minPhi,maxPhi, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
156  ptres_vs_pt.push_back(ibooker.book2D("ptres_vs_pt","ptres_vs_pt",nintpT,minpT,maxpT, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
157 
158  cotThetares_vs_eta.push_back(ibooker.book2D("cotThetares_vs_eta","cotThetares_vs_eta",nint,min,max,cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
159  cotThetares_vs_pt.push_back(ibooker.book2D("cotThetares_vs_pt","cotThetares_vs_pt",nintpT,minpT,maxpT, cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
160 
161  phires_vs_eta.push_back(ibooker.book2D("phires_vs_eta","phires_vs_eta",nint,min,max, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
162  phires_vs_pt.push_back(ibooker.book2D("phires_vs_pt","phires_vs_pt",nintpT,minpT,maxpT, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
163  phires_vs_phi.push_back(ibooker.book2D("phires_vs_phi","#phi res vs #phi",nintPhi,minPhi,maxPhi,phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
164 
165  dxyres_vs_eta.push_back(ibooker.book2D("dxyres_vs_eta","dxyres_vs_eta",nint,min,max,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
166  dxyres_vs_pt.push_back( ibooker.book2D("dxyres_vs_pt","dxyres_vs_pt",nintpT,minpT,maxpT,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
167 
168  dzres_vs_eta.push_back(ibooker.book2D("dzres_vs_eta","dzres_vs_eta",nint,min,max,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
169  dzres_vs_pt.push_back(ibooker.book2D("dzres_vs_pt","dzres_vs_pt",nintpT,minpT,maxpT,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
170 
171  ptmean_vs_eta_phi.push_back(ibooker.bookProfile2D("ptmean_vs_eta_phi","mean p_{t} vs #eta and #phi",nintPhi,minPhi,maxPhi,nint,min,max,1000,0,1000));
172  phimean_vs_eta_phi.push_back(ibooker.bookProfile2D("phimean_vs_eta_phi","mean #phi vs #eta and #phi",nintPhi,minPhi,maxPhi,nint,min,max,nintPhi,minPhi,maxPhi));
173 
174  //pulls of track params vs eta: to be used with fitslicesytool
175  dxypull_vs_eta.push_back(ibooker.book2D("dxypull_vs_eta","dxypull_vs_eta",nint,min,max,100,-10,10));
176  ptpull_vs_eta.push_back(ibooker.book2D("ptpull_vs_eta","ptpull_vs_eta",nint,min,max,100,-10,10));
177  dzpull_vs_eta.push_back(ibooker.book2D("dzpull_vs_eta","dzpull_vs_eta",nint,min,max,100,-10,10));
178  phipull_vs_eta.push_back(ibooker.book2D("phipull_vs_eta","phipull_vs_eta",nint,min,max,100,-10,10));
179  thetapull_vs_eta.push_back(ibooker.book2D("thetapull_vs_eta","thetapull_vs_eta",nint,min,max,100,-10,10));
180 
181  //pulls of track params vs phi
182  ptpull_vs_phi.push_back(ibooker.book2D("ptpull_vs_phi","p_{t} pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
183  phipull_vs_phi.push_back(ibooker.book2D("phipull_vs_phi","#phi pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
184  thetapull_vs_phi.push_back(ibooker.book2D("thetapull_vs_phi","#theta pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
185 
186  nrecHit_vs_nsimHit_sim2rec.push_back( ibooker.book2D("nrecHit_vs_nsimHit_sim2rec","nrecHit vs nsimHit (Sim2RecAssoc)",nintHit,minHit,maxHit, nintHit,minHit,maxHit ));
187  nrecHit_vs_nsimHit_rec2sim.push_back( ibooker.book2D("nrecHit_vs_nsimHit_rec2sim","nrecHit vs nsimHit (Rec2simAssoc)",nintHit,minHit,maxHit, nintHit,minHit,maxHit ));
188 
189  if (MABH) {
190  h_PurityVsQuality.push_back
191  (ibooker.book2D("PurityVsQuality","Purity vs Quality (MABH)",20,0.01,1.01,20,0.01,1.01) );
192  h_assoceta_Quality05.push_back
193  (ibooker.book1D("num_assoc(simToReco)_eta_Q05","N of associated tracks (simToReco) vs eta (Quality>0.5)",nint,min,max) );
194  h_assoceta_Quality075.push_back
195  (ibooker.book1D("num_assoc(simToReco)_eta_Q075","N of associated tracks (simToReco) vs eta (Quality>0.75)",nint,min,max) );
196  h_assocpT_Quality05.push_back
197  (ibooker.book1D("num_assoc(simToReco)_pT_Q05","N of associated tracks (simToReco) vs pT (Quality>0.5)",nintpT,minpT,maxpT) );
198  h_assocpT_Quality075.push_back
199  (ibooker.book1D("num_assoc(simToReco)_pT_Q075","N of associated tracks (simToReco) vs pT (Quality>0.75)",nintpT,minpT,maxpT) );
200  h_assocphi_Quality05.push_back
201  (ibooker.book1D("num_assoc(simToReco)_phi_Q05","N of associated tracks (simToReco) vs phi (Quality>0.5)",nintPhi,minPhi,maxPhi) );
202  h_assocphi_Quality075.push_back
203  (ibooker.book1D("num_assoc(simToReco)_phi_Q075","N of associated tracks (simToReco) vs phi (Quality>0.75)",nintPhi,minPhi,maxPhi) );
204  }
205 
206  if(useLogPt){
207  BinLogX(dzres_vs_pt[j]->getTH2F());
208  BinLogX(dxyres_vs_pt[j]->getTH2F());
209  BinLogX(phires_vs_pt[j]->getTH2F());
210  BinLogX(cotThetares_vs_pt[j]->getTH2F());
211  BinLogX(ptres_vs_pt[j]->getTH2F());
212  BinLogX(h_recopT[j]->getTH1F());
213  BinLogX(h_assocpT[j]->getTH1F());
214  BinLogX(h_assoc2pT[j]->getTH1F());
215  BinLogX(h_simulpT[j]->getTH1F());
216  if (MABH) {
217  BinLogX(h_assocpT_Quality05[j]->getTH1F());
218  BinLogX(h_assocpT_Quality075[j]->getTH1F());
219  }
220  j++;
221  }
222 
223  }
224  }
225 }
MuonTrackValidatorBase::h_assoc2phi
std::vector< MonitorElement * > h_assoc2phi
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidatorBase::h_recoeta
std::vector< MonitorElement * > h_recoeta
Definition: MuonTrackValidatorBase.h:396
MuonTrackValidatorBase::nrecHit_vs_nsimHit_sim2rec
std::vector< MonitorElement * > nrecHit_vs_nsimHit_sim2rec
Definition: MuonTrackValidatorBase.h:412
MuonTrackValidatorBase::h_DThits_eta
std::vector< MonitorElement * > h_DThits_eta
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::h_assoc2hit
std::vector< MonitorElement * > h_assoc2hit
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidatorBase::h_PurityVsQuality
std::vector< MonitorElement * > h_PurityVsQuality
Definition: MuonTrackValidatorBase.h:442
MuonTrackValidatorBase::h_recopT
std::vector< MonitorElement * > h_recopT
Definition: MuonTrackValidatorBase.h:397
MuonTrackValidator::etares_vs_eta
std::vector< MonitorElement * > etares_vs_eta
Definition: MuonTrackValidator.h:180
MuonTrackValidatorBase::h_assoc2dxy
std::vector< MonitorElement * > h_assoc2dxy
Definition: MuonTrackValidatorBase.h:400
MuonTrackValidatorBase::nrecHit_vs_nsimHit_rec2sim
std::vector< MonitorElement * > nrecHit_vs_nsimHit_rec2sim
Definition: MuonTrackValidatorBase.h:413
MuonTrackValidator::ptres_vs_phi
std::vector< MonitorElement * > ptres_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidator::h_losthits_eta
std::vector< MonitorElement * > h_losthits_eta
Definition: MuonTrackValidator.h:192
MuonTrackValidatorBase::h_assocpT_Quality075
std::vector< MonitorElement * > h_assocpT_Quality075
Definition: MuonTrackValidatorBase.h:406
MuonTrackValidatorBase::nrec_vs_nsim
std::vector< MonitorElement * > nrec_vs_nsim
Definition: MuonTrackValidatorBase.h:411
MuonTrackValidatorBase::h_pullDz
std::vector< MonitorElement * > h_pullDz
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidatorBase::h_recohit
std::vector< MonitorElement * > h_recohit
Definition: MuonTrackValidatorBase.h:398
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
MuonTrackValidatorBase::h_etaSIM
std::vector< MonitorElement * > h_etaSIM
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidatorBase::h_assocvertpos
std::vector< MonitorElement * > h_assocvertpos
Definition: MuonTrackValidatorBase.h:402
DQMStore::IBooker::cd
void cd(void)
Definition: DQMStore.cc:266
MuonTrackValidatorBase::h_fakes
std::vector< MonitorElement * > h_fakes
Definition: MuonTrackValidatorBase.h:395
MuonTrackValidator::ptres_vs_eta
std::vector< MonitorElement * > ptres_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidator::cotThetares_vs_eta
std::vector< MonitorElement * > cotThetares_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidatorBase::h_eta
std::vector< MonitorElement * > h_eta
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidator::phimean_vs_eta_phi
std::vector< MonitorElement * > phimean_vs_eta_phi
Definition: MuonTrackValidator.h:185
MuonTrackValidatorBase::h_pt
std::vector< MonitorElement * > h_pt
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidator::phires_vs_pt
std::vector< MonitorElement * > phires_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidatorBase::h_pullQoverp
std::vector< MonitorElement * > h_pullQoverp
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidator::h_assochi2_prob
std::vector< MonitorElement * > h_assochi2_prob
Definition: MuonTrackValidator.h:188
MuonTrackValidatorBase::h_CSChits_eta
std::vector< MonitorElement * > h_CSChits_eta
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidator::cotThetares_vs_pt
std::vector< MonitorElement * > cotThetares_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidator::dxyres_vs_eta
std::vector< MonitorElement * > dxyres_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidator::dzres_vs_pt
std::vector< MonitorElement * > dzres_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidatorBase::h_assoceta
std::vector< MonitorElement * > h_assoceta
Definition: MuonTrackValidatorBase.h:396
MuonTrackValidatorBase::h_tracks
std::vector< MonitorElement * > h_tracks
Definition: MuonTrackValidatorBase.h:395
MuonTrackValidator::dzres_vs_eta
std::vector< MonitorElement * > dzres_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidatorBase::h_assoceta_Quality075
std::vector< MonitorElement * > h_assoceta_Quality075
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidatorBase::nDThits_vs_eta
std::vector< MonitorElement * > nDThits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidator::thetapull_vs_phi
std::vector< MonitorElement * > thetapull_vs_phi
Definition: MuonTrackValidator.h:202
MuonTrackValidator::h_assochi2
std::vector< MonitorElement * > h_assochi2
Definition: MuonTrackValidator.h:188
MuonTrackValidator::chi2_vs_phi
std::vector< MonitorElement * > chi2_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::h_simuldxy
std::vector< MonitorElement * > h_simuldxy
Definition: MuonTrackValidatorBase.h:400
MuonTrackValidatorBase::h_simuldz
std::vector< MonitorElement * > h_simuldz
Definition: MuonTrackValidatorBase.h:401
MuonTrackValidator::chi2_vs_nhits
std::vector< MonitorElement * > chi2_vs_nhits
Definition: MuonTrackValidator.h:180
MuonTrackValidatorBase::nhits_vs_eta
std::vector< MonitorElement * > nhits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidator::phires_vs_phi
std::vector< MonitorElement * > phires_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidator::h_chi2mean_vs_phi
std::vector< MonitorElement * > h_chi2mean_vs_phi
Definition: MuonTrackValidator.h:194
MuonTrackValidatorBase::h_assoczpos
std::vector< MonitorElement * > h_assoczpos
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidator::nhits_vs_phi
std::vector< MonitorElement * > nhits_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::h_simulphi
std::vector< MonitorElement * > h_simulphi
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidatorBase::h_vertposSIM
std::vector< MonitorElement * > h_vertposSIM
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidator::dxyres_vs_pt
std::vector< MonitorElement * > dxyres_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidatorBase::label
std::vector< edm::InputTag > label
Definition: MuonTrackValidatorBase.h:349
MuonTrackValidatorBase::h_assocSharedHit
std::vector< MonitorElement * > h_assocSharedHit
Definition: MuonTrackValidatorBase.h:416
MuonTrackValidator::h_nchi2
std::vector< MonitorElement * > h_nchi2
Definition: MuonTrackValidator.h:177
DQMStore::IBooker::bookProfile2D
MonitorElement * bookProfile2D(Args &&...args)
Definition: DQMStore.h:163
TrackerOfflineValidation_Dqm_cff.dirName
tuple dirName
Definition: TrackerOfflineValidation_Dqm_cff.py:56
MuonTrackValidatorBase::h_assoc2dz
std::vector< MonitorElement * > h_assoc2dz
Definition: MuonTrackValidatorBase.h:401
MuonTrackValidatorBase::h_pullPhi
std::vector< MonitorElement * > h_pullPhi
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidator::h_nchi2_prob
std::vector< MonitorElement * > h_nchi2_prob
Definition: MuonTrackValidator.h:177
MuonTrackValidator::ptpull_vs_eta
std::vector< MonitorElement * > ptpull_vs_eta
Definition: MuonTrackValidator.h:201
MuonTrackValidatorBase::h_recophi
std::vector< MonitorElement * > h_recophi
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidator::dxypull_vs_eta
std::vector< MonitorElement * > dxypull_vs_eta
Definition: MuonTrackValidator.h:201
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
j
int j
Definition: DBlmapReader.cc:9
getTH1F
TH1F * getTH1F(std::string name, std::string process, std::string rootfolder, DQMStore *dbe_, bool verb, bool clone)
Definition: HcalHistoUtils.h:57
MuonTrackValidatorBase::h_pullTheta
std::vector< MonitorElement * > h_pullTheta
Definition: MuonTrackValidatorBase.h:403
MuonTrackValidatorBase::h_simulhit
std::vector< MonitorElement * > h_simulhit
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidator::phipull_vs_eta
std::vector< MonitorElement * > phipull_vs_eta
Definition: MuonTrackValidator.h:201
MuonTrackValidator::thetapull_vs_eta
std::vector< MonitorElement * > thetapull_vs_eta
Definition: MuonTrackValidator.h:201
MuonTrackValidatorBase::h_charge
std::vector< MonitorElement * > h_charge
Definition: MuonTrackValidatorBase.h:395
MuonTrackValidatorBase::h_assocphi
std::vector< MonitorElement * > h_assocphi
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidator::ptmean_vs_eta_phi
std::vector< MonitorElement * > ptmean_vs_eta_phi
Definition: MuonTrackValidator.h:185
MuonTrackValidatorBase::h_assocpT
std::vector< MonitorElement * > h_assocpT
Definition: MuonTrackValidatorBase.h:397
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
DQMStore::IBooker::goUp
void goUp(void)
Definition: DQMStore.cc:278
MuonTrackValidatorBase::h_pullDxy
std::vector< MonitorElement * > h_pullDxy
Definition: MuonTrackValidatorBase.h:403
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
MuonTrackValidatorBase::h_simulpT
std::vector< MonitorElement * > h_simulpT
Definition: MuonTrackValidatorBase.h:397
MuonTrackValidatorBase::h_tracksSIM
std::vector< MonitorElement * > h_tracksSIM
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidatorBase::h_assocphi_Quality05
std::vector< MonitorElement * > h_assocphi_Quality05
Definition: MuonTrackValidatorBase.h:407
MuonTrackValidatorBase::h_simulzpos
std::vector< MonitorElement * > h_simulzpos
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidatorBase::h_recodz
std::vector< MonitorElement * > h_recodz
Definition: MuonTrackValidatorBase.h:401
MuonTrackValidator::chi2_vs_eta
std::vector< MonitorElement * > chi2_vs_eta
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_assocphi_Quality075
std::vector< MonitorElement * > h_assocphi_Quality075
Definition: MuonTrackValidatorBase.h:407
HLT_25ns14e33_v1_cff.algo
tuple algo
Definition: HLT_25ns14e33_v1_cff.py:1972
MuonTrackValidatorBase::h_hits_eta
std::vector< MonitorElement * > h_hits_eta
Definition: MuonTrackValidatorBase.h:422
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:42
MuonTrackValidatorBase::h_hits
std::vector< MonitorElement * > h_hits
Definition: MuonTrackValidatorBase.h:395
MuonTrackValidator::h_losthits
std::vector< MonitorElement * > h_losthits
Definition: MuonTrackValidator.h:177
MuonTrackValidator::ptpull_vs_phi
std::vector< MonitorElement * > ptpull_vs_phi
Definition: MuonTrackValidator.h:202
edm::InputTag::process
std::string const & process() const
Definition: InputTag.h:46
MuonTrackValidatorBase::h_simulvertpos
std::vector< MonitorElement * > h_simulvertpos
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidator::ptres_vs_pt
std::vector< MonitorElement * > ptres_vs_pt
Definition: MuonTrackValidator.h:198
MuonTrackValidator::h_chi2meanhitsh
std::vector< MonitorElement * > h_chi2meanhitsh
Definition: MuonTrackValidator.h:194
MuonTrackValidatorBase::h_recodxy
std::vector< MonitorElement * > h_recodxy
Definition: MuonTrackValidatorBase.h:400
MuonTrackValidator::phires_vs_eta
std::vector< MonitorElement * > phires_vs_eta
Definition: MuonTrackValidator.h:197
MuonTrackValidatorBase::h_ptSIM
std::vector< MonitorElement * > h_ptSIM
Definition: MuonTrackValidatorBase.h:392
getTH2F
TH2F * getTH2F(std::string name, std::string process, std::string rootfolder, DQMStore *dbe_, bool verb, bool clone)
Definition: HcalHistoUtils.h:85
MuonTrackValidatorBase::h_assocdz
std::vector< MonitorElement * > h_assocdz
Definition: MuonTrackValidatorBase.h:401
MuonTrackValidator::h_hits_phi
std::vector< MonitorElement * > h_hits_phi
Definition: MuonTrackValidator.h:193
MuonTrackValidatorBase::h_assocpT_Quality05
std::vector< MonitorElement * > h_assocpT_Quality05
Definition: MuonTrackValidatorBase.h:406
MuonTrackValidatorBase::h_RPChits_eta
std::vector< MonitorElement * > h_RPChits_eta
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::h_assoc2pT
std::vector< MonitorElement * > h_assoc2pT
Definition: MuonTrackValidatorBase.h:397
MuonTrackValidatorBase::h_assocFraction
std::vector< MonitorElement * > h_assocFraction
Definition: MuonTrackValidatorBase.h:416
MuonTrackValidatorBase::associators
std::vector< std::string > associators
Definition: MuonTrackValidatorBase.h:355
MuonTrackValidator::nlosthits_vs_eta
std::vector< MonitorElement * > nlosthits_vs_eta
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_assoceta_Quality05
std::vector< MonitorElement * > h_assoceta_Quality05
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidator::dzpull_vs_eta
std::vector< MonitorElement * > dzpull_vs_eta
Definition: MuonTrackValidator.h:201
MuonTrackValidatorBase::h_assocdxy
std::vector< MonitorElement * > h_assocdxy
Definition: MuonTrackValidatorBase.h:400
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:43
MuonTrackValidatorBase::nRPChits_vs_eta
std::vector< MonitorElement * > nRPChits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::h_assochit
std::vector< MonitorElement * > h_assochit
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidator::h_chi2meanh
std::vector< MonitorElement * > h_chi2meanh
Definition: MuonTrackValidator.h:192
MuonTrackValidatorBase::nCSChits_vs_eta
std::vector< MonitorElement * > nCSChits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::h_simuleta
std::vector< MonitorElement * > h_simuleta
Definition: MuonTrackValidatorBase.h:396
MuonTrackValidatorBase::h_assoc2eta
std::vector< MonitorElement * > h_assoc2eta
Definition: MuonTrackValidatorBase.h:396
MuonTrackValidator::phipull_vs_phi
std::vector< MonitorElement * > phipull_vs_phi
Definition: MuonTrackValidator.h:202
void MuonTrackValidator::endRun ( edm::Run const &  ,
edm::EventSetup const &   
)
virtual

Method called at the end of the event loop.

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 890 of file MuonTrackValidator.cc.

References dbe_, diffTwoXMLs::label, dbtoconf::out, and w.

890  {
891 
892  int w=0;
893  for (unsigned int ww=0;ww<associators.size();ww++){
894  for (unsigned int www=0;www<label.size();www++){
895 
896  //chi2 and #hit vs eta: get mean from 2D histos
897  doProfileX(chi2_vs_eta[w],h_chi2meanh[w]);
898  doProfileX(nhits_vs_eta[w],h_hits_eta[w]);
899  doProfileX(nDThits_vs_eta[w],h_DThits_eta[w]);
900  doProfileX(nCSChits_vs_eta[w],h_CSChits_eta[w]);
901  doProfileX(nRPChits_vs_eta[w],h_RPChits_eta[w]);
902 
903  doProfileX(nlosthits_vs_eta[w],h_losthits_eta[w]);
904  //vs phi
905  doProfileX(chi2_vs_nhits[w],h_chi2meanhitsh[w]);
906  doProfileX(chi2_vs_phi[w],h_chi2mean_vs_phi[w]);
907  doProfileX(nhits_vs_phi[w],h_hits_phi[w]);
908 
909  fillPlotFromVector(h_recoeta[w],totRECeta[w]);
910  fillPlotFromVector(h_simuleta[w],totSIMeta[w]);
911  fillPlotFromVector(h_assoceta[w],totASSeta[w]);
912  fillPlotFromVector(h_assoc2eta[w],totASS2eta[w]);
913 
914  fillPlotFromVector(h_recopT[w],totRECpT[w]);
915  fillPlotFromVector(h_simulpT[w],totSIMpT[w]);
916  fillPlotFromVector(h_assocpT[w],totASSpT[w]);
917  fillPlotFromVector(h_assoc2pT[w],totASS2pT[w]);
918 
919  fillPlotFromVector(h_recohit[w],totREC_hit[w]);
920  fillPlotFromVector(h_simulhit[w],totSIM_hit[w]);
921  fillPlotFromVector(h_assochit[w],totASS_hit[w]);
922  fillPlotFromVector(h_assoc2hit[w],totASS2_hit[w]);
923 
924  fillPlotFromVector(h_recophi[w],totREC_phi[w]);
925  fillPlotFromVector(h_simulphi[w],totSIM_phi[w]);
926  fillPlotFromVector(h_assocphi[w],totASS_phi[w]);
927  fillPlotFromVector(h_assoc2phi[w],totASS2_phi[w]);
928 
929  fillPlotFromVector(h_recodxy[w],totREC_dxy[w]);
930  fillPlotFromVector(h_simuldxy[w],totSIM_dxy[w]);
931  fillPlotFromVector(h_assocdxy[w],totASS_dxy[w]);
932  fillPlotFromVector(h_assoc2dxy[w],totASS2_dxy[w]);
933 
934  fillPlotFromVector(h_recodz[w],totREC_dz[w]);
935  fillPlotFromVector(h_simuldz[w],totSIM_dz[w]);
936  fillPlotFromVector(h_assocdz[w],totASS_dz[w]);
937  fillPlotFromVector(h_assoc2dz[w],totASS2_dz[w]);
938 
939  fillPlotFromVector(h_simulvertpos[w],totSIM_vertpos[w]);
940  fillPlotFromVector(h_assocvertpos[w],totASS_vertpos[w]);
941 
942  fillPlotFromVector(h_simulzpos[w],totSIM_zpos[w]);
943  fillPlotFromVector(h_assoczpos[w],totASS_zpos[w]);
944 
945  if (MABH) {
946  fillPlotFromVector(h_assoceta_Quality05[w] ,totASSeta_Quality05[w]);
947  fillPlotFromVector(h_assoceta_Quality075[w],totASSeta_Quality075[w]);
948  fillPlotFromVector(h_assocpT_Quality05[w] ,totASSpT_Quality05[w]);
949  fillPlotFromVector(h_assocpT_Quality075[w],totASSpT_Quality075[w]);
950  fillPlotFromVector(h_assocphi_Quality05[w] ,totASS_phi_Quality05[w]);
951  fillPlotFromVector(h_assocphi_Quality075[w],totASS_phi_Quality075[w]);
952  }
953 
954  w++;
955  }
956  }
957 
958  if ( out.size() != 0 && dbe_ ) dbe_->save(out);
959 }
MuonTrackValidatorBase::h_assoc2phi
std::vector< MonitorElement * > h_assoc2phi
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidatorBase::h_recoeta
std::vector< MonitorElement * > h_recoeta
Definition: MuonTrackValidatorBase.h:396
MuonTrackValidatorBase::h_DThits_eta
std::vector< MonitorElement * > h_DThits_eta
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::h_assoc2hit
std::vector< MonitorElement * > h_assoc2hit
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidatorBase::totASS_zpos
std::vector< std::vector< int > > totASS_zpos
Definition: MuonTrackValidatorBase.h:439
MuonTrackValidatorBase::totASSeta_Quality05
std::vector< std::vector< int > > totASSeta_Quality05
Definition: MuonTrackValidatorBase.h:443
MuonTrackValidatorBase::h_recopT
std::vector< MonitorElement * > h_recopT
Definition: MuonTrackValidatorBase.h:397
MuonTrackValidatorBase::totREC_phi
std::vector< std::vector< int > > totREC_phi
Definition: MuonTrackValidatorBase.h:436
MuonTrackValidatorBase::h_assoc2dxy
std::vector< MonitorElement * > h_assoc2dxy
Definition: MuonTrackValidatorBase.h:400
MuonTrackValidatorBase::totASS2_hit
std::vector< std::vector< int > > totASS2_hit
Definition: MuonTrackValidatorBase.h:435
w
const double w
Definition: UKUtility.cc:23
MuonTrackValidator::h_losthits_eta
std::vector< MonitorElement * > h_losthits_eta
Definition: MuonTrackValidator.h:192
MuonTrackValidatorBase::h_assocpT_Quality075
std::vector< MonitorElement * > h_assocpT_Quality075
Definition: MuonTrackValidatorBase.h:406
MuonTrackValidatorBase::totRECpT
std::vector< std::vector< int > > totRECpT
Definition: MuonTrackValidatorBase.h:434
MuonTrackValidatorBase::totREC_dz
std::vector< std::vector< int > > totREC_dz
Definition: MuonTrackValidatorBase.h:438
MuonTrackValidatorBase::h_recohit
std::vector< MonitorElement * > h_recohit
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidatorBase::h_assocvertpos
std::vector< MonitorElement * > h_assocvertpos
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidatorBase::totASS_dz
std::vector< std::vector< int > > totASS_dz
Definition: MuonTrackValidatorBase.h:438
MuonTrackValidatorBase::totSIM_zpos
std::vector< std::vector< int > > totSIM_zpos
Definition: MuonTrackValidatorBase.h:439
MuonTrackValidatorBase::h_CSChits_eta
std::vector< MonitorElement * > h_CSChits_eta
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::totASS2_dz
std::vector< std::vector< int > > totASS2_dz
Definition: MuonTrackValidatorBase.h:438
MuonTrackValidatorBase::totASS_phi_Quality05
std::vector< std::vector< int > > totASS_phi_Quality05
Definition: MuonTrackValidatorBase.h:445
MuonTrackValidatorBase::h_assoceta
std::vector< MonitorElement * > h_assoceta
Definition: MuonTrackValidatorBase.h:396
MuonTrackValidatorBase::h_assoceta_Quality075
std::vector< MonitorElement * > h_assoceta_Quality075
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidatorBase::nDThits_vs_eta
std::vector< MonitorElement * > nDThits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::totASS_phi
std::vector< std::vector< int > > totASS_phi
Definition: MuonTrackValidatorBase.h:436
MuonTrackValidatorBase::totASS_dxy
std::vector< std::vector< int > > totASS_dxy
Definition: MuonTrackValidatorBase.h:437
MuonTrackValidator::chi2_vs_phi
std::vector< MonitorElement * > chi2_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::h_simuldxy
std::vector< MonitorElement * > h_simuldxy
Definition: MuonTrackValidatorBase.h:400
MuonTrackValidatorBase::totSIM_dxy
std::vector< std::vector< int > > totSIM_dxy
Definition: MuonTrackValidatorBase.h:437
MuonTrackValidatorBase::h_simuldz
std::vector< MonitorElement * > h_simuldz
Definition: MuonTrackValidatorBase.h:401
MuonTrackValidator::chi2_vs_nhits
std::vector< MonitorElement * > chi2_vs_nhits
Definition: MuonTrackValidator.h:180
MuonTrackValidatorBase::totASSeta_Quality075
std::vector< std::vector< int > > totASSeta_Quality075
Definition: MuonTrackValidatorBase.h:443
MuonTrackValidatorBase::nhits_vs_eta
std::vector< MonitorElement * > nhits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidator::h_chi2mean_vs_phi
std::vector< MonitorElement * > h_chi2mean_vs_phi
Definition: MuonTrackValidator.h:194
MuonTrackValidatorBase::h_assoczpos
std::vector< MonitorElement * > h_assoczpos
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidator::nhits_vs_phi
std::vector< MonitorElement * > nhits_vs_phi
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::h_simulphi
std::vector< MonitorElement * > h_simulphi
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidatorBase::totASS2pT
std::vector< std::vector< int > > totASS2pT
Definition: MuonTrackValidatorBase.h:434
MuonTrackValidatorBase::doProfileX
virtual void doProfileX(TH2 *th2, MonitorElement *me)
Definition: MuonTrackValidatorBase.h:117
MuonTrackValidatorBase::label
std::vector< edm::InputTag > label
Definition: MuonTrackValidatorBase.h:349
MuonTrackValidatorBase::totASS2eta
std::vector< std::vector< int > > totASS2eta
Definition: MuonTrackValidatorBase.h:433
MuonTrackValidatorBase::totSIMpT
std::vector< std::vector< int > > totSIMpT
Definition: MuonTrackValidatorBase.h:434
MuonTrackValidatorBase::totREC_dxy
std::vector< std::vector< int > > totREC_dxy
Definition: MuonTrackValidatorBase.h:437
MuonTrackValidatorBase::h_assoc2dz
std::vector< MonitorElement * > h_assoc2dz
Definition: MuonTrackValidatorBase.h:401
MuonTrackValidatorBase::totSIM_hit
std::vector< std::vector< int > > totSIM_hit
Definition: MuonTrackValidatorBase.h:435
MuonTrackValidatorBase::h_recophi
std::vector< MonitorElement * > h_recophi
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidatorBase::totASSpT
std::vector< std::vector< int > > totASSpT
Definition: MuonTrackValidatorBase.h:434
MuonTrackValidatorBase::totASS_phi_Quality075
std::vector< std::vector< int > > totASS_phi_Quality075
Definition: MuonTrackValidatorBase.h:445
MuonTrackValidatorBase::totSIM_vertpos
std::vector< std::vector< int > > totSIM_vertpos
Definition: MuonTrackValidatorBase.h:439
MuonTrackValidatorBase::totASSpT_Quality05
std::vector< std::vector< int > > totASSpT_Quality05
Definition: MuonTrackValidatorBase.h:444
MuonTrackValidatorBase::totREC_hit
std::vector< std::vector< int > > totREC_hit
Definition: MuonTrackValidatorBase.h:435
MuonTrackValidatorBase::h_simulhit
std::vector< MonitorElement * > h_simulhit
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidatorBase::totSIMeta
std::vector< std::vector< int > > totSIMeta
Definition: MuonTrackValidatorBase.h:433
MuonTrackValidatorBase::totASS2_dxy
std::vector< std::vector< int > > totASS2_dxy
Definition: MuonTrackValidatorBase.h:437
MuonTrackValidatorBase::h_assocphi
std::vector< MonitorElement * > h_assocphi
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidatorBase::totASS2_phi
std::vector< std::vector< int > > totASS2_phi
Definition: MuonTrackValidatorBase.h:436
MuonTrackValidatorBase::totASSpT_Quality075
std::vector< std::vector< int > > totASSpT_Quality075
Definition: MuonTrackValidatorBase.h:444
MuonTrackValidatorBase::h_assocpT
std::vector< MonitorElement * > h_assocpT
Definition: MuonTrackValidatorBase.h:397
MuonTrackValidatorBase::h_simulpT
std::vector< MonitorElement * > h_simulpT
Definition: MuonTrackValidatorBase.h:397
MuonTrackValidatorBase::h_assocphi_Quality05
std::vector< MonitorElement * > h_assocphi_Quality05
Definition: MuonTrackValidatorBase.h:407
MuonTrackValidatorBase::h_simulzpos
std::vector< MonitorElement * > h_simulzpos
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidatorBase::h_recodz
std::vector< MonitorElement * > h_recodz
Definition: MuonTrackValidatorBase.h:401
MuonTrackValidator::chi2_vs_eta
std::vector< MonitorElement * > chi2_vs_eta
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_assocphi_Quality075
std::vector< MonitorElement * > h_assocphi_Quality075
Definition: MuonTrackValidatorBase.h:407
MuonTrackValidatorBase::totSIM_phi
std::vector< std::vector< int > > totSIM_phi
Definition: MuonTrackValidatorBase.h:436
MuonTrackValidatorBase::fillPlotFromVector
void fillPlotFromVector(MonitorElement *h, std::vector< int > &vec)
Definition: MuonTrackValidatorBase.h:141
MuonTrackValidatorBase::h_hits_eta
std::vector< MonitorElement * > h_hits_eta
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::totASSeta
std::vector< std::vector< int > > totASSeta
Definition: MuonTrackValidatorBase.h:433
MuonTrackValidatorBase::h_simulvertpos
std::vector< MonitorElement * > h_simulvertpos
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidatorBase::totRECeta
std::vector< std::vector< int > > totRECeta
Definition: MuonTrackValidatorBase.h:433
MuonTrackValidator::h_chi2meanhitsh
std::vector< MonitorElement * > h_chi2meanhitsh
Definition: MuonTrackValidator.h:194
MuonTrackValidatorBase::h_recodxy
std::vector< MonitorElement * > h_recodxy
Definition: MuonTrackValidatorBase.h:400
MuonTrackValidatorBase::h_assocdz
std::vector< MonitorElement * > h_assocdz
Definition: MuonTrackValidatorBase.h:401
MuonTrackValidator::h_hits_phi
std::vector< MonitorElement * > h_hits_phi
Definition: MuonTrackValidator.h:193
MuonTrackValidatorBase::h_assocpT_Quality05
std::vector< MonitorElement * > h_assocpT_Quality05
Definition: MuonTrackValidatorBase.h:406
MuonTrackValidatorBase::h_RPChits_eta
std::vector< MonitorElement * > h_RPChits_eta
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::h_assoc2pT
std::vector< MonitorElement * > h_assoc2pT
Definition: MuonTrackValidatorBase.h:397
MuonTrackValidatorBase::totSIM_dz
std::vector< std::vector< int > > totSIM_dz
Definition: MuonTrackValidatorBase.h:438
MuonTrackValidatorBase::associators
std::vector< std::string > associators
Definition: MuonTrackValidatorBase.h:355
MuonTrackValidator::nlosthits_vs_eta
std::vector< MonitorElement * > nlosthits_vs_eta
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_assoceta_Quality05
std::vector< MonitorElement * > h_assoceta_Quality05
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidatorBase::h_assocdxy
std::vector< MonitorElement * > h_assocdxy
Definition: MuonTrackValidatorBase.h:400
MuonTrackValidatorBase::nRPChits_vs_eta
std::vector< MonitorElement * > nRPChits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::totASS_vertpos
std::vector< std::vector< int > > totASS_vertpos
Definition: MuonTrackValidatorBase.h:439
MuonTrackValidatorBase::h_assochit
std::vector< MonitorElement * > h_assochit
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidatorBase::totASS_hit
std::vector< std::vector< int > > totASS_hit
Definition: MuonTrackValidatorBase.h:435
MuonTrackValidator::h_chi2meanh
std::vector< MonitorElement * > h_chi2meanh
Definition: MuonTrackValidator.h:192
MuonTrackValidatorBase::nCSChits_vs_eta
std::vector< MonitorElement * > nCSChits_vs_eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::h_simuleta
std::vector< MonitorElement * > h_simuleta
Definition: MuonTrackValidatorBase.h:396
MuonTrackValidatorBase::h_assoc2eta
std::vector< MonitorElement * > h_assoc2eta
Definition: MuonTrackValidatorBase.h:396
void MuonTrackValidator::getRecoMomentum ( const reco::Track &  track,
double &  pt,
double &  ptError,
double &  qoverp,
double &  qoverpError,
double &  lambda,
double &  lambdaError,
double &  phi,
double &  phiError 
) const
private

retrieval of reconstructed momentum components from reco::Track (== mean values for GSF)

Definition at line 961 of file MuonTrackValidator.cc.

References reco::TrackBase::lambda(), reco::TrackBase::lambdaError(), reco::TrackBase::phi(), reco::TrackBase::phiError(), reco::TrackBase::pt(), reco::TrackBase::ptError(), reco::TrackBase::qoverp(), and reco::TrackBase::qoverpError().

963  {
964  pt = track.pt();
965  ptError = track.ptError();
966  qoverp = track.qoverp();
967  qoverpError = track.qoverpError();
968  lambda = track.lambda();
969  lambdaError = track.lambdaError();
970  phi = track.phi();
971  phiError = track.phiError();
972 }
reco::TrackBase::qoverp
double qoverp() const
q / p
Definition: TrackBase.h:621
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:693
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:669
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:811
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:838
reco::TrackBase::lambda
double lambda() const
Lambda angle.
Definition: TrackBase.h:633
reco::TrackBase::qoverpError
double qoverpError() const
error on signed transverse curvature
Definition: TrackBase.h:805
reco::TrackBase::lambdaError
double lambdaError() const
error on lambda
Definition: TrackBase.h:826
phi
Definition: DDAxes.h:10
void MuonTrackValidator::getRecoMomentum ( const reco::GsfTrack &  gsfTrack,
double &  pt,
double &  ptError,
double &  qoverp,
double &  qoverpError,
double &  lambda,
double &  lambdaError,
double &  phi,
double &  phiError 
) const
private

retrieval of reconstructed momentum components based on the mode of a reco::GsfTrack

Definition at line 974 of file MuonTrackValidator.cc.

References reco::GsfTrack::lambdaMode(), reco::GsfTrack::lambdaModeError(), reco::GsfTrack::phiMode(), reco::GsfTrack::phiModeError(), reco::GsfTrack::ptMode(), reco::GsfTrack::ptModeError(), reco::GsfTrack::qoverpMode(), and reco::GsfTrack::qoverpModeError().

976  {
977  pt = gsfTrack.ptMode();
978  ptError = gsfTrack.ptModeError();
979  qoverp = gsfTrack.qoverpMode();
980  qoverpError = gsfTrack.qoverpModeError();
981  lambda = gsfTrack.lambdaMode();
982  lambdaError = gsfTrack.lambdaModeError();
983  phi = gsfTrack.phiMode();
984  phiError = gsfTrack.phiModeError();
985 }
reco::GsfTrack::phiModeError
double phiModeError() const
error on phi from mode
Definition: GsfTrack.h:94
reco::GsfTrack::lambdaMode
double lambdaMode() const
Lambda angle from mode.
Definition: GsfTrack.h:45
reco::GsfTrack::qoverpMode
double qoverpMode() const
q/p from mode
Definition: GsfTrack.h:41
reco::GsfTrack::ptModeError
double ptModeError() const
error on Pt (set to 1000 TeV if charge==0 for safety) from mode
Definition: GsfTrack.h:81
reco::GsfTrack::qoverpModeError
double qoverpModeError() const
error on signed transverse curvature from mode
Definition: GsfTrack.h:79
reco::GsfTrack::phiMode
double phiMode() const
azimuthal angle of momentum vector from mode
Definition: GsfTrack.h:57
reco::GsfTrack::lambdaModeError
double lambdaModeError() const
error on lambda from mode
Definition: GsfTrack.h:90
reco::GsfTrack::ptMode
double ptMode() const
track transverse momentum from mode
Definition: GsfTrack.h:49
phi
Definition: DDAxes.h:10

Member Data Documentation

edm::EDGetTokenT<SimHitTPAssociationProducer::SimHitTPAssociationList> MuonTrackValidator::_simHitTpMapTag
private

Definition at line 160 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

edm::InputTag MuonTrackValidator::associatormap
private

Definition at line 157 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::BiDirectional_RecoToSim_association
private

Definition at line 172 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

std::vector<MonitorElement*> MuonTrackValidator::chi2_vs_eta
private

Definition at line 191 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::chi2_vs_nhits
private

Definition at line 180 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::chi2_vs_phi
private

Definition at line 182 of file MuonTrackValidator.h.

CosmicTrackingParticleSelector MuonTrackValidator::cosmictpSelector
private

Definition at line 169 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

std::vector<MonitorElement*> MuonTrackValidator::cotThetares_vs_eta
private

Definition at line 197 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::cotThetares_vs_pt
private

Definition at line 198 of file MuonTrackValidator.h.

std::string MuonTrackValidator::dirName_
private

Definition at line 156 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

std::vector<MonitorElement*> MuonTrackValidator::dxypull_vs_eta
private

Definition at line 201 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::dxyres_vs_eta
private

Definition at line 197 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::dxyres_vs_pt
private

Definition at line 198 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::dzpull_vs_eta
private

Definition at line 201 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::dzres_vs_eta
private

Definition at line 197 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::dzres_vs_pt
private

Definition at line 198 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::etares_vs_eta
private

Definition at line 180 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_assochi2
private

Definition at line 188 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_assochi2_prob
private

Definition at line 188 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_chi2mean_vs_phi
private

Definition at line 194 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_chi2meanh
private

Definition at line 192 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_chi2meanhitsh
private

Definition at line 194 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_dxypulleta
private

Definition at line 203 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_dzpulleta
private

Definition at line 203 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_hits_phi
private

Definition at line 193 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_losthits
private

Definition at line 177 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_losthits_eta
private

Definition at line 192 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_nchi2
private

Definition at line 177 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_nchi2_prob
private

Definition at line 177 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_phipulleta
private

Definition at line 203 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_phipullphi
private

Definition at line 204 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_ptpulleta
private

Definition at line 203 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_ptpullphi
private

Definition at line 204 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_ptshifteta
private

Definition at line 181 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_thetapulleta
private

Definition at line 203 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::h_thetapullphi
private

Definition at line 204 of file MuonTrackValidator.h.

bool MuonTrackValidator::MABH
private

Definition at line 174 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

double MuonTrackValidator::maxPhi
private

Definition at line 163 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

double MuonTrackValidator::minPhi
private

Definition at line 163 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

std::vector<MonitorElement*> MuonTrackValidator::nhits_vs_phi
private

Definition at line 182 of file MuonTrackValidator.h.

int MuonTrackValidator::nintPhi
private

Definition at line 164 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

std::vector<MonitorElement*> MuonTrackValidator::nlosthits_vs_eta
private

Definition at line 191 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::phimean_vs_eta_phi
private

Definition at line 185 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::phipull_vs_eta
private

Definition at line 201 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::phipull_vs_phi
private

Definition at line 202 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::phires_vs_eta
private

Definition at line 197 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::phires_vs_phi
private

Definition at line 182 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::phires_vs_pt
private

Definition at line 198 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::ptmean_vs_eta_phi
private

Definition at line 185 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::ptpull_vs_eta
private

Definition at line 201 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::ptpull_vs_phi
private

Definition at line 202 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::ptres_vs_eta
private

Definition at line 197 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::ptres_vs_phi
private

Definition at line 182 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::ptres_vs_pt
private

Definition at line 198 of file MuonTrackValidator.h.

edm::EDGetTokenT<reco::RecoToSimCollection> MuonTrackValidator::recoToSimCollection_Token
private

Definition at line 159 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

edm::EDGetTokenT<reco::SimToRecoCollection> MuonTrackValidator::simToRecoCollection_Token
private

Definition at line 158 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

std::vector<MonitorElement*> MuonTrackValidator::thetapull_vs_eta
private

Definition at line 201 of file MuonTrackValidator.h.

std::vector<MonitorElement*> MuonTrackValidator::thetapull_vs_phi
private

Definition at line 202 of file MuonTrackValidator.h.

TrackingParticleSelector MuonTrackValidator::tpSelector
private

Definition at line 168 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::UseAssociators
private

Definition at line 162 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::useGsf
private

Definition at line 165 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

Generated for CMSSW Reference Manual by doxygen 1.8.5