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:
edm::EDAnalyzer MuonTrackValidatorBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &)
 Method called once per event. More...
 
void beginRun (edm::Run const &, edm::EventSetup const &)
 Method called before the event loop. More...
 
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 edm::EDAnalyzer
 EDAnalyzer ()
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 

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::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_cotThetameanh
 
std::vector< MonitorElement * > h_cotThetameanhPt
 
std::vector< MonitorElement * > h_cotThetarmsh
 
std::vector< MonitorElement * > h_cotThetarmshPt
 
std::vector< MonitorElement * > h_dxymeanh
 
std::vector< MonitorElement * > h_dxymeanhPt
 
std::vector< MonitorElement * > h_dxypulleta
 
std::vector< MonitorElement * > h_dxypulletamean
 
std::vector< MonitorElement * > h_dxyrmsh
 
std::vector< MonitorElement * > h_dxyrmshPt
 
std::vector< MonitorElement * > h_dzmeanh
 
std::vector< MonitorElement * > h_dzmeanhPt
 
std::vector< MonitorElement * > h_dzpulleta
 
std::vector< MonitorElement * > h_dzpulletamean
 
std::vector< MonitorElement * > h_dzrmsh
 
std::vector< MonitorElement * > h_dzrmshPt
 
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_phimeanh
 
std::vector< MonitorElement * > h_phimeanhPhi
 
std::vector< MonitorElement * > h_phimeanhPt
 
std::vector< MonitorElement * > h_phipulleta
 
std::vector< MonitorElement * > h_phipulletamean
 
std::vector< MonitorElement * > h_phipullphi
 
std::vector< MonitorElement * > h_phipullphimean
 
std::vector< MonitorElement * > h_phirmsh
 
std::vector< MonitorElement * > h_phirmshPhi
 
std::vector< MonitorElement * > h_phirmshPt
 
std::vector< MonitorElement * > h_ptmeanh
 
std::vector< MonitorElement * > h_ptmeanhPhi
 
std::vector< MonitorElement * > h_ptmeanhPt
 
std::vector< MonitorElement * > h_ptpulleta
 
std::vector< MonitorElement * > h_ptpulletamean
 
std::vector< MonitorElement * > h_ptpullphi
 
std::vector< MonitorElement * > h_ptpullphimean
 
std::vector< MonitorElement * > h_ptrmsh
 
std::vector< MonitorElement * > h_ptrmshPhi
 
std::vector< MonitorElement * > h_ptrmshPt
 
std::vector< MonitorElement * > h_ptshifteta
 
std::vector< MonitorElement * > h_thetapulleta
 
std::vector< MonitorElement * > h_thetapulletamean
 
std::vector< MonitorElement * > h_thetapullphi
 
std::vector< MonitorElement * > h_thetapullphimean
 
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
 
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::EDAnalyzer
typedef EDAnalyzer ModuleType
 
typedef WorkerT< EDAnalyzer > WorkerType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const * currentContext () const
 
- 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)
 Constructor. More...
 
void setUpVectors ()
 
virtual ~MuonTrackValidatorBase ()
 Destructor. More...
 
- Protected Attributes inherited from MuonTrackValidatorBase
std::vector< const
TrackAssociatorBase * > 		associator
 
std::vector< std::string > associators
 
edm::InputTag bsSrc
 
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_effic
 
std::vector< MonitorElement * > h_effic_Quality05
 
std::vector< MonitorElement * > h_effic_Quality075
 
std::vector< MonitorElement * > h_effic_vs_dxy
 
std::vector< MonitorElement * > h_effic_vs_dz
 
std::vector< MonitorElement * > h_effic_vs_hit
 
std::vector< MonitorElement * > h_effic_vs_phi
 
std::vector< MonitorElement * > h_effic_vs_phi_Quality05
 
std::vector< MonitorElement * > h_effic_vs_phi_Quality075
 
std::vector< MonitorElement * > h_effic_vs_vertpos
 
std::vector< MonitorElement * > h_effic_vs_zpos
 
std::vector< MonitorElement * > h_efficPt
 
std::vector< MonitorElement * > h_efficPt_Quality05
 
std::vector< MonitorElement * > h_efficPt_Quality075
 
std::vector< MonitorElement * > h_eta
 
std::vector< MonitorElement * > h_etaSIM
 
std::vector< MonitorElement * > h_fake_vs_dxy
 
std::vector< MonitorElement * > h_fake_vs_dz
 
std::vector< MonitorElement * > h_fake_vs_hit
 
std::vector< MonitorElement * > h_fake_vs_phi
 
std::vector< MonitorElement * > h_fakerate
 
std::vector< MonitorElement * > h_fakeratePt
 
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
 
bool skipHistoFit
 
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
 
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

Date:
2010/03/25 10:27:54
Revision:
1.2

Definition at line 17 of file MuonTrackValidator.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 20 of file MuonTrackValidator.h.

References associatormap, MuonTrackValidatorBase::associators, BiDirectional_RecoToSim_association, cosmictpSelector, dirName_, edm::ParameterSet::dump(), edm::ParameterSet::getParameter(), MuonTrackValidatorBase::label, MABH, maxPhi, minPhi, nintPhi, ExpressReco_HICollisions_FallBack::recoTracksLabel, tpSelector, UseAssociators, useGsf, MuonTrackValidatorBase::usemuon, and MuonTrackValidatorBase::usetracker.

20  :MuonTrackValidatorBase(pset){
21  dirName_ = pset.getParameter<std::string>("dirName");
22  associatormap = pset.getParameter< edm::InputTag >("associatormap");
23  UseAssociators = pset.getParameter< bool >("UseAssociators");
24  tpSelector = TrackingParticleSelector(pset.getParameter<double>("ptMinTP"),
25  pset.getParameter<double>("minRapidityTP"),
26  pset.getParameter<double>("maxRapidityTP"),
27  pset.getParameter<double>("tipTP"),
28  pset.getParameter<double>("lipTP"),
29  pset.getParameter<int>("minHitTP"),
30  pset.getParameter<bool>("signalOnlyTP"),
31  pset.getParameter<bool>("chargedOnlyTP"),
32  pset.getParameter<std::vector<int> >("pdgIdTP"));
33  cosmictpSelector = CosmicTrackingParticleSelector(pset.getParameter<double>("ptMinTP"),
34  pset.getParameter<double>("minRapidityTP"),
35  pset.getParameter<double>("maxRapidityTP"),
36  pset.getParameter<double>("tipTP"),
37  pset.getParameter<double>("lipTP"),
38  pset.getParameter<int>("minHitTP"),
39  pset.getParameter<bool>("chargedOnlyTP"),
40  pset.getParameter<std::vector<int> >("pdgIdTP"));
41 
42  minPhi = pset.getParameter<double>("minPhi");
43  maxPhi = pset.getParameter<double>("maxPhi");
44  nintPhi = pset.getParameter<int>("nintPhi");
45  useGsf = pset.getParameter<bool>("useGsf");
46  BiDirectional_RecoToSim_association = pset.getParameter<bool>("BiDirectional_RecoToSim_association");
47 
48  // dump cfg parameters
49  edm::LogVerbatim("MuonTrackValidator") << "constructing MuonTrackValidator: " << pset.dump();
50 
51  MABH = false;
52  if (!UseAssociators) {
53  // flag MuonAssociatorByHits
54  if (associators[0] == "MuonAssociationByHits") MABH = true;
55  // reset string associators to the map label
56  associators.clear();
57  associators.push_back(associatormap.label());
58  edm::LogVerbatim("MuonTrackValidator") << "--> associators reset to: " <<associators[0];
59  }
60 
61  // inform on which SimHits will be counted
62  if (usetracker) edm::LogVerbatim("MuonTrackValidator")
63  <<"\n usetracker = TRUE : Tracker SimHits WILL be counted";
64  else edm::LogVerbatim("MuonTrackValidator")
65  <<"\n usetracker = FALSE : Tracker SimHits WILL NOT be counted";
66  if (usemuon) edm::LogVerbatim("MuonTrackValidator")
67  <<" usemuon = TRUE : Muon SimHits WILL be counted";
68  else edm::LogVerbatim("MuonTrackValidator")
69  <<" usemuon = FALSE : Muon SimHits WILL NOT be counted"<<std::endl;
70 
71  // loop over the reco::Track collections to validate: check for inconsistent input settings
72  for (unsigned int www=0;www<label.size();www++) {
73  std::string recoTracksLabel = label[www].label();
74  std::string recoTracksInstance = label[www].instance();
75 
76  // tracks with hits only on tracker
77  if (recoTracksLabel=="generalTracks" ||
78  (recoTracksLabel.find("cutsRecoTracks") != std::string::npos) ||
79  recoTracksLabel=="ctfWithMaterialTracksP5LHCNavigation" ||
80  recoTracksLabel=="hltL3TkTracksFromL2" ||
81  (recoTracksLabel=="hltL3Muons" && recoTracksInstance=="L2Seeded"))
82  {
83  if (usemuon) {
84  edm::LogWarning("MuonTrackValidator")
85  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
86  <<"\n with usemuon == true"<<"\n ---> please change to usemuon == false ";
87  }
88  if (!usetracker) {
89  edm::LogWarning("MuonTrackValidator")
90  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
91  <<"\n with usetracker == false"<<"\n ---> please change to usetracker == true ";
92  }
93  }
94 
95  // tracks with hits only on muon detectors
96  else if (recoTracksLabel=="standAloneMuons" ||
97  recoTracksLabel=="standAloneSETMuons" ||
98  recoTracksLabel=="cosmicMuons" ||
99  recoTracksLabel=="hltL2Muons")
100  {
101  if (usetracker) {
102  edm::LogWarning("MuonTrackValidator")
103  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
104  <<"\n with usetracker == true"<<"\n ---> please change to usetracker == false ";
105  }
106  if (!usemuon) {
107  edm::LogWarning("MuonTrackValidator")
108  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
109  <<"\n with usemuon == false"<<"\n ---> please change to usemuon == true ";
110  }
111  }
112 
113  } // for (unsigned int www=0;www<label.size();www++)
114  }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
ExpressReco_HICollisions_FallBack.recoTracksLabel
tuple recoTracksLabel
Definition: ExpressReco_HICollisions_FallBack.py:11702
reco::modules::TrackingParticleSelector
SingleObjectSelector< TrackingParticleCollection,::TrackingParticleSelector > TrackingParticleSelector
Definition: TrackingParticleSelector.cc:18
MuonTrackValidator::cosmictpSelector
CosmicTrackingParticleSelector cosmictpSelector
Definition: MuonTrackValidator.h:146
MuonTrackValidatorBase::label
std::vector< edm::InputTag > label
Definition: MuonTrackValidatorBase.h:338
reco::CosmicTrackingParticleSelector
ObjectSelector< CosmicTrackingParticleSelector > CosmicTrackingParticleSelector
Definition: CosmicTrackingParticleSelector.cc:16
MuonTrackValidator::tpSelector
TrackingParticleSelector tpSelector
Definition: MuonTrackValidator.h:145
MuonTrackValidator::associatormap
edm::InputTag associatormap
Definition: MuonTrackValidator.h:138
MuonTrackValidatorBase::MuonTrackValidatorBase
MuonTrackValidatorBase(const edm::ParameterSet &pset)
Constructor.
Definition: MuonTrackValidatorBase.h:39
MuonTrackValidatorBase::associators
std::vector< std::string > associators
Definition: MuonTrackValidatorBase.h:344
edm::ParameterSet::dump
std::string dump() const
Definition: ParameterSet.cc:817
virtual MuonTrackValidator::~MuonTrackValidator ( )
inlinevirtual

Destructor.

Definition at line 117 of file MuonTrackValidator.h.

117 { }

Member Function Documentation

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

Method called once per event.

Implements edm::EDAnalyzer.

Definition at line 332 of file MuonTrackValidator.cc.

References abs, asciidump::at, reco::TrackBase::charge(), reco::TrackBase::chi2(), funct::cos(), reco::TrackBase::dxy(), reco::TrackBase::dxyError(), reco::TrackBase::dz(), reco::TrackBase::dzError(), ExpressReco_HICollisions_FallBack::e, edm::AssociationMap< Tag >::end(), reco::TrackBase::eta(), f, edm::AssociationMap< Tag >::find(), reco::Track::found(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), reco::TrackBase::hitPattern(), i, edm::RefToBase< T >::key(), edm::Ref< C, T, F >::key(), label, LogTrace, M_PI, min, ParticleBase::momentum(), reco::TrackBase::momentum(), DetId::Muon, 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(), ExpressReco_HICollisions_FallBack::Quality, dt_offlineAnalysis_common_cff::reco, funct::sin(), mathSSE::sqrt(), reco::TrackBase::theta(), tmp, ExpressReco_HICollisions_FallBack::track, ExpressReco_HICollisions_FallBack::trackCollection, DetId::Tracker, TrackingParticle::trackPSimHit(), ParticleBase::vertex(), and cms::Exception::what().

332  {
333  using namespace reco;
334 
335  edm::LogInfo("MuonTrackValidator") << "\n====================================================" << "\n"
336  << "Analyzing new event" << "\n"
337  << "====================================================\n" << "\n";
338  edm::ESHandle<ParametersDefinerForTP> parametersDefinerTP;
339  setup.get<TrackAssociatorRecord>().get(parametersDefiner,parametersDefinerTP);
340 
341  edm::Handle<TrackingParticleCollection> TPCollectionHeff ;
342  event.getByLabel(label_tp_effic,TPCollectionHeff);
343  const TrackingParticleCollection tPCeff = *(TPCollectionHeff.product());
344 
345  edm::Handle<TrackingParticleCollection> TPCollectionHfake ;
346  event.getByLabel(label_tp_fake,TPCollectionHfake);
347  const TrackingParticleCollection tPCfake = *(TPCollectionHfake.product());
348 
349  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
350  event.getByLabel(bsSrc,recoBeamSpotHandle);
351  reco::BeamSpot bs = *recoBeamSpotHandle;
352 
353  int w=0;
354  for (unsigned int ww=0;ww<associators.size();ww++){
355  for (unsigned int www=0;www<label.size();www++){
356  //
357  //get collections from the event
358  //
359  edm::Handle<View<Track> > trackCollection;
360 
361  reco::RecoToSimCollection recSimColl;
362  reco::SimToRecoCollection simRecColl;
363  unsigned int trackCollectionSize = 0;
364 
365  // if(!event.getByLabel(label[www], trackCollection)&&ignoremissingtkcollection_) continue;
366  if(!event.getByLabel(label[www], trackCollection)&&ignoremissingtkcollection_) {
367 
368  recSimColl.post_insert();
369  simRecColl.post_insert();
370 
371  }
372 
373  else {
374 
375  trackCollectionSize = trackCollection->size();
376  //associate tracks
377  if(UseAssociators){
378  edm::LogVerbatim("MuonTrackValidator") << "Analyzing "
379  << label[www].process()<<":"
380  << label[www].label()<<":"
381  << label[www].instance()<<" with "
382  << associators[ww].c_str() <<"\n";
383 
384  LogTrace("MuonTrackValidator") << "Calling associateRecoToSim method" << "\n";
385  recSimColl=associator[ww]->associateRecoToSim(trackCollection,
386  TPCollectionHfake,
387  &event);
388  LogTrace("MuonTrackValidator") << "Calling associateSimToReco method" << "\n";
389  simRecColl=associator[ww]->associateSimToReco(trackCollection,
390  TPCollectionHeff,
391  &event);
392  }
393  else{
394  edm::LogVerbatim("MuonTrackValidator") << "Analyzing "
395  << label[www].process()<<":"
396  << label[www].label()<<":"
397  << label[www].instance()<<" with "
398  << associatormap.process()<<":"
399  << associatormap.label()<<":"
400  << associatormap.instance()<<"\n";
401 
402  Handle<reco::SimToRecoCollection > simtorecoCollectionH;
403  event.getByLabel(associatormap,simtorecoCollectionH);
404  simRecColl= *(simtorecoCollectionH.product());
405 
406  Handle<reco::RecoToSimCollection > recotosimCollectionH;
407  event.getByLabel(associatormap,recotosimCollectionH);
408  recSimColl= *(recotosimCollectionH.product());
409  }
410 
411  }
412 
413 
414  //
415  //fill simulation histograms
416  //compute number of tracks per eta interval
417  //
418  edm::LogVerbatim("MuonTrackValidator") << "\n# of TrackingParticles: " << tPCeff.size() << "\n";
419  int ats = 0;
420  int st=0;
421  for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
422  bool TP_is_matched = false;
423  double quality = 0.;
424  bool Quality05 = false;
425  bool Quality075 = false;
426 
427  TrackingParticleRef tpr(TPCollectionHeff, i);
428  TrackingParticle* tp=const_cast<TrackingParticle*>(tpr.get());
429  ParticleBase::Vector momentumTP;
430  ParticleBase::Point vertexTP;
431  double dxySim = 0;
432  double dzSim = 0;
433 
434  //If the TrackingParticle is collison like, get the momentum and vertex at production state
435  if(parametersDefiner=="LhcParametersDefinerForTP")
436  {
437  if(! tpSelector(*tp)) continue;
438  momentumTP = tp->momentum();
439  vertexTP = tp->vertex();
440  //Calcualte the impact parameters w.r.t. PCA
441  ParticleBase::Vector momentum = parametersDefinerTP->momentum(event,setup,*tp);
442  ParticleBase::Point vertex = parametersDefinerTP->vertex(event,setup,*tp);
443  dxySim = (-vertex.x()*sin(momentum.phi())+vertex.y()*cos(momentum.phi()));
444  dzSim = vertex.z() - (vertex.x()*momentum.x()+vertex.y()*momentum.y())/sqrt(momentum.perp2()) * momentum.z()/sqrt(momentum.perp2());
445  }
446  //If the TrackingParticle is comics, get the momentum and vertex at PCA
447  if(parametersDefiner=="CosmicParametersDefinerForTP")
448  {
449  if(! cosmictpSelector(*tp,&bs,event,setup)) continue;
450  momentumTP = parametersDefinerTP->momentum(event,setup,*tp);
451  vertexTP = parametersDefinerTP->vertex(event,setup,*tp);
452  dxySim = (-vertexTP.x()*sin(momentumTP.phi())+vertexTP.y()*cos(momentumTP.phi()));
453  dzSim = vertexTP.z() - (vertexTP.x()*momentumTP.x()+vertexTP.y()*momentumTP.y())/sqrt(momentumTP.perp2()) * momentumTP.z()/sqrt(momentumTP.perp2());
454  }
455  edm::LogVerbatim("MuonTrackValidator") <<"--------------------Selected TrackingParticle #"<<tpr.key();
456  st++;
457 
458  h_ptSIM[w]->Fill(sqrt(momentumTP.perp2()));
459  h_etaSIM[w]->Fill(momentumTP.eta());
460  h_vertposSIM[w]->Fill(sqrt(vertexTP.perp2()));
461 
462  std::vector<std::pair<RefToBase<Track>, double> > rt;
463  if(simRecColl.find(tpr) != simRecColl.end()){
464  rt = (std::vector<std::pair<RefToBase<Track>, double> >) simRecColl[tpr];
465  if (rt.size()!=0) {
466  RefToBase<Track> assoc_recoTrack = rt.begin()->first;
467  edm::LogVerbatim("MuonTrackValidator")<<"-----------------------------associated Track #"<<assoc_recoTrack.key();
468  TP_is_matched = true;
469  ats++;
470  quality = rt.begin()->second;
471  edm::LogVerbatim("MuonTrackValidator") << "TrackingParticle #" <<tpr.key()
472  << " with pt=" << sqrt(momentumTP.perp2())
473  << " associated with quality:" << quality <<"\n";
474  if (MABH) {
475  if (quality > 0.75) {
476  Quality075 = true;
477  Quality05 = true;
478  }
479  else if (quality > 0.5) {
480  Quality05 = true;
481  }
482  }
483  }
484  }else{
485  edm::LogVerbatim("MuonTrackValidator")
486  << "TrackingParticle #" << tpr.key()
487  << " with pt,eta,phi: "
488  << sqrt(momentumTP.perp2()) << " , "
489  << momentumTP.eta() << " , "
490  << momentumTP.phi() << " , "
491  << " NOT associated to any reco::Track" << "\n";
492  }
493 
494  for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
495  if (getEta(momentumTP.eta())>etaintervals[w][f]&&
496  getEta(momentumTP.eta())<etaintervals[w][f+1]) {
497  totSIMeta[w][f]++;
498  if (TP_is_matched) {
499  totASSeta[w][f]++;
500 
501  if (MABH) {
502  if (Quality075) {
503  totASSeta_Quality075[w][f]++;
504  totASSeta_Quality05[w][f]++;
505  }
506  else if (Quality05) {
507  totASSeta_Quality05[w][f]++;
508  }
509  }
510  }
511  }
512  } // END for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
513 
514  for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
515  if (momentumTP.phi() > phiintervals[w][f]&&
516  momentumTP.phi() <phiintervals[w][f+1]) {
517  totSIM_phi[w][f]++;
518  if (TP_is_matched) {
519  totASS_phi[w][f]++;
520 
521  if (MABH) {
522  if (Quality075) {
523  totASS_phi_Quality075[w][f]++;
524  totASS_phi_Quality05[w][f]++;
525  }
526  else if (Quality05) {
527  totASS_phi_Quality05[w][f]++;
528  }
529  }
530  }
531  }
532  } // END for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
533 
534 
535  for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
536  if (getPt(sqrt(momentumTP.perp2()))>pTintervals[w][f]&&
537  getPt(sqrt(momentumTP.perp2()))<pTintervals[w][f+1]) {
538  totSIMpT[w][f]++;
539  if (TP_is_matched) {
540  totASSpT[w][f]++;
541 
542  if (MABH) {
543  if (Quality075) {
544  totASSpT_Quality075[w][f]++;
545  totASSpT_Quality05[w][f]++;
546  }
547  else if (Quality05) {
548  totASSpT_Quality05[w][f]++;
549  }
550  }
551  }
552  }
553  } // END for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
554 
555  for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
556  if (dxySim>dxyintervals[w][f]&&
557  dxySim<dxyintervals[w][f+1]) {
558  totSIM_dxy[w][f]++;
559  if (TP_is_matched) {
560  totASS_dxy[w][f]++;
561  }
562  }
563  } // END for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
564 
565  for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
566  if (dzSim>dzintervals[w][f]&&
567  dzSim<dzintervals[w][f+1]) {
568  totSIM_dz[w][f]++;
569  if (TP_is_matched) {
570  totASS_dz[w][f]++;
571  }
572  }
573  } // END for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
574 
575  for (unsigned int f=0; f<vertposintervals[w].size()-1; f++){
576  if (sqrt(vertexTP.perp2())>vertposintervals[w][f]&&
577  sqrt(vertexTP.perp2())<vertposintervals[w][f+1]) {
578  totSIM_vertpos[w][f]++;
579  if (TP_is_matched) {
580  totASS_vertpos[w][f]++;
581  }
582  }
583  } // END for (unsigned int f=0; f<vertposintervals[w].size()-1; f++){
584 
585  for (unsigned int f=0; f<zposintervals[w].size()-1; f++){
586  if (vertexTP.z()>zposintervals[w][f]&&
587  vertexTP.z()<zposintervals[w][f+1]) {
588  totSIM_zpos[w][f]++;
589  if (TP_is_matched) {
590  totASS_zpos[w][f]++;
591  }
592  }
593  } // END for (unsigned int f=0; f<zposintervals[w].size()-1; f++){
594 
595  std::vector<PSimHit> simhits;
596 
597  if (usetracker && usemuon) {
598  simhits=tp->trackPSimHit();
599  }
600  else if (!usetracker && usemuon) {
601  simhits=tp->trackPSimHit(DetId::Muon);
602  }
603  else if (usetracker && !usemuon) {
604  simhits=tp->trackPSimHit(DetId::Tracker);
605  }
606 
607  int tmp = std::min((int)(simhits.end()-simhits.begin()),int(maxHit-1));
608  edm::LogVerbatim("MuonTrackValidator") << "\t N simhits = "<< (int)(simhits.end()-simhits.begin())<<"\n";
609 
610  totSIM_hit[w][tmp]++;
611  if (TP_is_matched) totASS_hit[w][tmp]++;
612 
613  if (TP_is_matched)
614  {
615  RefToBase<Track> assoctrack = rt.begin()->first;
616  nrecHit_vs_nsimHit_sim2rec[w]->Fill( assoctrack->numberOfValidHits(),(int)(simhits.end()-simhits.begin() ));
617  }
618  } // End for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
619  if (st!=0) h_tracksSIM[w]->Fill(st);
620 
621 
622  //
623  //fill reconstructed track histograms
624  //
625  edm::LogVerbatim("MuonTrackValidator") << "\n# of reco::Tracks with "
626  << label[www].process()<<":"
627  << label[www].label()<<":"
628  << label[www].instance()
629  << ": " << trackCollectionSize << "\n";
630  int at=0;
631  int rT=0;
632  for(View<Track>::size_type i=0; i<trackCollectionSize; ++i){
633  bool Track_is_matched = false;
634  RefToBase<Track> track(trackCollection, i);
635  rT++;
636 
637  std::vector<std::pair<TrackingParticleRef, double> > tp;
638  TrackingParticleRef tpr;
639 
640  // new logic (bidirectional)
641  if (BiDirectional_RecoToSim_association) {
642  edm::LogVerbatim("MuonTrackValidator")<<"----------------------------------------Track #"<< track.key();
643 
644  if(recSimColl.find(track) != recSimColl.end()) {
645  tp = recSimColl[track];
646  if (tp.size() != 0) {
647  tpr = tp.begin()->first;
648  // RtS and StR must associate the same pair !
649  if(simRecColl.find(tpr) != simRecColl.end()) {
650  std::vector<std::pair<RefToBase<Track>, double> > track_checkback = simRecColl[tpr];
651  RefToBase<Track> assoc_track_checkback;
652  assoc_track_checkback = track_checkback.begin()->first;
653 
654  if ( assoc_track_checkback.key() == track.key() ) {
655  edm::LogVerbatim("MuonTrackValidator")<<"------------------associated TrackingParticle #"<<tpr.key();
656  Track_is_matched = true;
657  at++;
658  double Purity = tp.begin()->second;
659  double Quality = track_checkback.begin()->second;
660  edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
661  << " associated with quality:" << Purity <<"\n";
662  if (MABH) h_PurityVsQuality[w]->Fill(Quality,Purity);
663  }
664  }
665  }
666  }
667 
668  if (!Track_is_matched) edm::LogVerbatim("MuonTrackValidator")
669  << "reco::Track #" << track.key() << " with pt=" << track->pt() << " NOT associated to any TrackingParticle" << "\n";
670  }
671  // old logic (bugged)
672  else {
673  if(recSimColl.find(track) != recSimColl.end()){
674  tp = recSimColl[track];
675  if (tp.size()!=0) {
676  Track_is_matched = true;
677  tpr = tp.begin()->first;
678  at++;
679  edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
680  << " associated with quality:" << tp.begin()->second <<"\n";
681  }
682  } else {
683  edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
684  << " NOT associated to any TrackingParticle" << "\n";
685  }
686  }
687 
688  //Compute fake rate vs eta
689  for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
690  if (getEta(track->momentum().eta())>etaintervals[w][f]&&
691  getEta(track->momentum().eta())<etaintervals[w][f+1]) {
692  totRECeta[w][f]++;
693  if (Track_is_matched) {
694  totASS2eta[w][f]++;
695  }
696  }
697  } // End for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
698 
699  for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
700  if (track->momentum().phi()>phiintervals[w][f]&&
701  track->momentum().phi()<phiintervals[w][f+1]) {
702  totREC_phi[w][f]++;
703  if (Track_is_matched) {
704  totASS2_phi[w][f]++;
705  }
706  }
707  } // End for (unsigned int f=0; f<phiintervals[w].size()-1; f++){
708 
709 
710  for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
711  if (getPt(sqrt(track->momentum().perp2()))>pTintervals[w][f]&&
712  getPt(sqrt(track->momentum().perp2()))<pTintervals[w][f+1]) {
713  totRECpT[w][f]++;
714  if (Track_is_matched) {
715  totASS2pT[w][f]++;
716  }
717  }
718  } // End for (unsigned int f=0; f<pTintervals[w].size()-1; f++){
719 
720  for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
721  if (track->dxy(bs.position())>dxyintervals[w][f]&&
722  track->dxy(bs.position())<dxyintervals[w][f+1]) {
723  totREC_dxy[w][f]++;
724  if (Track_is_matched) {
725  totASS2_dxy[w][f]++;
726  }
727  }
728  } // End for (unsigned int f=0; f<dxyintervals[w].size()-1; f++){
729 
730  for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
731  if (track->dz(bs.position())>dzintervals[w][f]&&
732  track->dz(bs.position())<dzintervals[w][f+1]) {
733  totREC_dz[w][f]++;
734  if (Track_is_matched) {
735  totASS2_dz[w][f]++;
736  }
737  }
738  } // End for (unsigned int f=0; f<dzintervals[w].size()-1; f++){
739 
740  int tmp = std::min((int)track->found(),int(maxHit-1));
741  totREC_hit[w][tmp]++;
742  if (Track_is_matched) totASS2_hit[w][tmp]++;
743 
744  edm::LogVerbatim("MuonTrackValidator") << "\t N valid rechits = "<< (int)track->found() <<"\n";
745 
746  //Fill other histos
747  try{
748  if (!Track_is_matched) continue;
749 
750  if (associators[ww]=="TrackAssociatorByChi2"){
751  //association chi2
752  double assocChi2 = -tp.begin()->second;//in association map is stored -chi2
753  h_assochi2[www]->Fill(assocChi2);
754  h_assochi2_prob[www]->Fill(TMath::Prob((assocChi2)*5,5));
755  }
756  else if (associators[ww]=="TrackAssociatorByHits"){
757  double fraction = tp.begin()->second;
758  h_assocFraction[www]->Fill(fraction);
759  h_assocSharedHit[www]->Fill(fraction*track->numberOfValidHits());
760  }
761 
762  //nchi2 and hits global distributions
763  h_nchi2[w]->Fill(track->normalizedChi2());
764  h_nchi2_prob[w]->Fill(TMath::Prob(track->chi2(),(int)track->ndof()));
765  h_hits[w]->Fill(track->numberOfValidHits());
766  h_losthits[w]->Fill(track->numberOfLostHits());
767  chi2_vs_nhits[w]->Fill(track->numberOfValidHits(),track->normalizedChi2());
768  h_charge[w]->Fill( track->charge() );
769 
770  //Get tracking particle parameters at point of closest approach to the beamline
771  ParticleBase::Vector momentumTP = parametersDefinerTP->momentum(event,setup,*(tpr.get()));
772  ParticleBase::Point vertexTP = parametersDefinerTP->vertex(event,setup,*(tpr.get()));
773  double ptSim = sqrt(momentumTP.perp2());
774  double qoverpSim = tpr->charge()/sqrt(momentumTP.x()*momentumTP.x()+momentumTP.y()*momentumTP.y()+momentumTP.z()*momentumTP.z());
775  double thetaSim = momentumTP.theta();
776  double lambdaSim = M_PI/2-momentumTP.theta();
777  double phiSim = momentumTP.phi();
778  double dxySim = (-vertexTP.x()*sin(momentumTP.phi())+vertexTP.y()*cos(momentumTP.phi()));
779  double dzSim = vertexTP.z() - (vertexTP.x()*momentumTP.x()+vertexTP.y()*momentumTP.y())/sqrt(momentumTP.perp2()) * momentumTP.z()/sqrt(momentumTP.perp2());
780 
781  TrackBase::ParameterVector rParameters = track->parameters();
782 
783  double qoverpRec(0);
784  double qoverpErrorRec(0);
785  double ptRec(0);
786  double ptErrorRec(0);
787  double lambdaRec(0);
788  double lambdaErrorRec(0);
789  double phiRec(0);
790  double phiErrorRec(0);
791 
792 
793  //loop to decide whether to take gsfTrack (utilisation of mode-function) or common track
794  const GsfTrack* gsfTrack(0);
795  if(useGsf){
796  gsfTrack = dynamic_cast<const GsfTrack*>(&(*track));
797  if (gsfTrack==0) edm::LogInfo("MuonTrackValidator") << "Trying to access mode for a non-GsfTrack";
798  }
799 
800  if (gsfTrack) {
801  // get values from mode
802  getRecoMomentum(*gsfTrack, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
803  lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
804  }
805 
806  else {
807  // get values from track (without mode)
808  getRecoMomentum(*track, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
809  lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
810  }
811 
812  double thetaRec = track->theta();
813  double ptError = ptErrorRec;
814  double ptres = ptRec - ptSim;
815  double etares = track->eta()-momentumTP.Eta();
816  double dxyRec = track->dxy(bs.position());
817  double dzRec = track->dz(bs.position());
818  // eta residue; pt, k, theta, phi, dxy, dz pulls
819  double qoverpPull=(qoverpRec-qoverpSim)/qoverpErrorRec;
820  double thetaPull=(lambdaRec-lambdaSim)/lambdaErrorRec;
821  double phiDiff = phiRec - phiSim;
822  if (abs(phiDiff) > M_PI) {
823  if (phiDiff >0.) phiDiff = phiDiff - 2.*M_PI;
824  else phiDiff = phiDiff + 2.*M_PI;
825  }
826  double phiPull=phiDiff/phiErrorRec;
827  double dxyPull=(dxyRec-dxySim)/track->dxyError();
828  double dzPull=(dzRec-dzSim)/track->dzError();
829 
830  double contrib_Qoverp = ((qoverpRec-qoverpSim)/qoverpErrorRec)*
831  ((qoverpRec-qoverpSim)/qoverpErrorRec)/5;
832  double contrib_dxy = ((dxyRec-dxySim)/track->dxyError())*((dxyRec-dxySim)/track->dxyError())/5;
833  double contrib_dz = ((dzRec-dzSim)/track->dzError())*((dzRec-dzSim)/track->dzError())/5;
834  double contrib_theta = ((lambdaRec-lambdaSim)/lambdaErrorRec)*
835  ((lambdaRec-lambdaSim)/lambdaErrorRec)/5;
836  double contrib_phi = (phiDiff/phiErrorRec)*(phiDiff/phiErrorRec)/5;
837 
838  LogTrace("MuonTrackValidator") << "assocChi2=" << tp.begin()->second << "\n"
839  << "" << "\n"
840  << "ptREC=" << ptRec << "\n"
841  << "etaREC=" << track->eta() << "\n"
842  << "qoverpREC=" << qoverpRec << "\n"
843  << "dxyREC=" << dxyRec << "\n"
844  << "dzREC=" << dzRec << "\n"
845  << "thetaREC=" << track->theta() << "\n"
846  << "phiREC=" << phiRec << "\n"
847  << "" << "\n"
848  << "qoverpError()=" << qoverpErrorRec << "\n"
849  << "dxyError()=" << track->dxyError() << "\n"
850  << "dzError()=" << track->dzError() << "\n"
851  << "thetaError()=" << lambdaErrorRec << "\n"
852  << "phiError()=" << phiErrorRec << "\n"
853  << "" << "\n"
854  << "ptSIM=" << ptSim << "\n"
855  << "etaSIM=" << momentumTP.Eta() << "\n"
856  << "qoverpSIM=" << qoverpSim << "\n"
857  << "dxySIM=" << dxySim << "\n"
858  << "dzSIM=" << dzSim << "\n"
859  << "thetaSIM=" << M_PI/2-lambdaSim << "\n"
860  << "phiSIM=" << phiSim << "\n"
861  << "" << "\n"
862  << "contrib_Qoverp=" << contrib_Qoverp << "\n"
863  << "contrib_dxy=" << contrib_dxy << "\n"
864  << "contrib_dz=" << contrib_dz << "\n"
865  << "contrib_theta=" << contrib_theta << "\n"
866  << "contrib_phi=" << contrib_phi << "\n"
867  << "" << "\n"
868  <<"chi2PULL="<<contrib_Qoverp+contrib_dxy+contrib_dz+contrib_theta+contrib_phi<<"\n";
869 
870  h_pullQoverp[w]->Fill(qoverpPull);
871  h_pullTheta[w]->Fill(thetaPull);
872  h_pullPhi[w]->Fill(phiPull);
873  h_pullDxy[w]->Fill(dxyPull);
874  h_pullDz[w]->Fill(dzPull);
875 
876 
877  h_pt[w]->Fill(ptres/ptError);
878  h_eta[w]->Fill(etares);
879  etares_vs_eta[w]->Fill(getEta(track->eta()),etares);
880 
881 
882  //chi2 and #hit vs eta: fill 2D histos
883  chi2_vs_eta[w]->Fill(getEta(track->eta()),track->normalizedChi2());
884  nhits_vs_eta[w]->Fill(getEta(track->eta()),track->numberOfValidHits());
885  nDThits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonDTHits());
886  nCSChits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonCSCHits());
887  nRPChits_vs_eta[w]->Fill(getEta(track->eta()),track->hitPattern().numberOfValidMuonRPCHits());
888 
889  nlosthits_vs_eta[w]->Fill(getEta(track->eta()),track->numberOfLostHits());
890 
891  //resolution of track params: fill 2D histos
892  dxyres_vs_eta[w]->Fill(getEta(track->eta()),dxyRec-dxySim);
893  ptres_vs_eta[w]->Fill(getEta(track->eta()),(ptRec-ptSim)/ptRec);
894  dzres_vs_eta[w]->Fill(getEta(track->eta()),dzRec-dzSim);
895  phires_vs_eta[w]->Fill(getEta(track->eta()),phiDiff);
896  cotThetares_vs_eta[w]->Fill(getEta(track->eta()), cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
897 
898  //same as before but vs pT
899  dxyres_vs_pt[w]->Fill(getPt(ptRec),dxyRec-dxySim);
900  ptres_vs_pt[w]->Fill(getPt(ptRec),(ptRec-ptSim)/ptRec);
901  dzres_vs_pt[w]->Fill(getPt(ptRec),dzRec-dzSim);
902  phires_vs_pt[w]->Fill(getPt(ptRec),phiDiff);
903  cotThetares_vs_pt[w]->Fill(getPt(ptRec), cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
904 
905  //pulls of track params vs eta: fill 2D histos
906  dxypull_vs_eta[w]->Fill(getEta(track->eta()),dxyPull);
907  ptpull_vs_eta[w]->Fill(getEta(track->eta()),ptres/ptError);
908  dzpull_vs_eta[w]->Fill(getEta(track->eta()),dzPull);
909  phipull_vs_eta[w]->Fill(getEta(track->eta()),phiPull);
910  thetapull_vs_eta[w]->Fill(getEta(track->eta()),thetaPull);
911 
912  //plots vs phi
913  nhits_vs_phi[w]->Fill(phiRec,track->numberOfValidHits());
914  chi2_vs_phi[w]->Fill(phiRec,track->normalizedChi2());
915  ptmean_vs_eta_phi[w]->Fill(phiRec,getEta(track->eta()),ptRec);
916  phimean_vs_eta_phi[w]->Fill(phiRec,getEta(track->eta()),phiRec);
917  ptres_vs_phi[w]->Fill(phiRec,(ptRec-ptSim)/ptRec);
918  phires_vs_phi[w]->Fill(phiRec,phiDiff);
919  ptpull_vs_phi[w]->Fill(phiRec,ptres/ptError);
920  phipull_vs_phi[w]->Fill(phiRec,phiPull);
921  thetapull_vs_phi[w]->Fill(phiRec,thetaPull);
922 
923  std::vector<PSimHit> simhits;
924 
925  if (usetracker && usemuon) {
926  simhits=tpr.get()->trackPSimHit();
927  }
928  else if (!usetracker && usemuon) {
929  simhits=tpr.get()->trackPSimHit(DetId::Muon);
930  }
931  else if (usetracker && !usemuon) {
932  simhits=tpr.get()->trackPSimHit(DetId::Tracker);
933  }
934 
935  nrecHit_vs_nsimHit_rec2sim[w]->Fill(track->numberOfValidHits(), (int)(simhits.end()-simhits.begin() ));
936 
937  } // End of try{
938  catch (cms::Exception e){
939  LogTrace("MuonTrackValidator") << "exception found: " << e.what() << "\n";
940  }
941  } // End of for(View<Track>::size_type i=0; i<trackCollectionSize; ++i){
942  if (at!=0) h_tracks[w]->Fill(at);
943  h_fakes[w]->Fill(rT-at);
944  edm::LogVerbatim("MuonTrackValidator") << "Total Simulated: " << st << "\n"
945  << "Total Associated (simToReco): " << ats << "\n"
946  << "Total Reconstructed: " << rT << "\n"
947  << "Total Associated (recoToSim): " << at << "\n"
948  << "Total Fakes: " << rT-at << "\n";
949  nrec_vs_nsim[w]->Fill(rT,st);
950  w++;
951  } // End of for (unsigned int www=0;www<label.size();www++){
952  } //END of for (unsigned int ww=0;ww<associators.size();ww++){
953 }
ExpressReco_HICollisions_FallBack.e
tuple e
Definition: ExpressReco_HICollisions_FallBack.py:1068
cms::Exception::what
virtual char const * what() const
Definition: Exception.cc:97
MuonTrackValidatorBase::nrecHit_vs_nsimHit_sim2rec
std::vector< MonitorElement * > nrecHit_vs_nsimHit_sim2rec
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidatorBase::totASS_zpos
std::vector< std::vector< int > > totASS_zpos
Definition: MuonTrackValidatorBase.h:425
i
int i
Definition: DBlmapReader.cc:9
MuonTrackValidatorBase::h_PurityVsQuality
std::vector< MonitorElement * > h_PurityVsQuality
Definition: MuonTrackValidatorBase.h:428
MuonTrackValidatorBase::totASSeta_Quality05
std::vector< std::vector< int > > totASSeta_Quality05
Definition: MuonTrackValidatorBase.h:429
ExpressReco_HICollisions_FallBack.track
tuple track
Definition: ExpressReco_HICollisions_FallBack.py:4641
MuonTrackValidatorBase::totREC_phi
std::vector< std::vector< int > > totREC_phi
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidator::etares_vs_eta
std::vector< MonitorElement * > etares_vs_eta
Definition: MuonTrackValidator.h:157
MuonTrackValidatorBase::nrecHit_vs_nsimHit_rec2sim
std::vector< MonitorElement * > nrecHit_vs_nsimHit_rec2sim
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidator::ptres_vs_phi
std::vector< MonitorElement * > ptres_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::totASS2_hit
std::vector< std::vector< int > > totASS2_hit
Definition: MuonTrackValidatorBase.h:421
ExpressReco_HICollisions_FallBack.trackCollection
tuple trackCollection
Definition: ExpressReco_HICollisions_FallBack.py:7877
MuonTrackValidatorBase::totRECpT
std::vector< std::vector< int > > totRECpT
Definition: MuonTrackValidatorBase.h:420
MuonTrackValidatorBase::totREC_dz
std::vector< std::vector< int > > totREC_dz
Definition: MuonTrackValidatorBase.h:424
MuonTrackValidatorBase::nrec_vs_nsim
std::vector< MonitorElement * > nrec_vs_nsim
Definition: MuonTrackValidatorBase.h:397
edm::AssociationMap::end
const_iterator end() const
last iterator over the map (read only)
Definition: AssociationMap.h:92
MuonTrackValidatorBase::h_pullDz
std::vector< MonitorElement * > h_pullDz
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidatorBase::h_etaSIM
std::vector< MonitorElement * > h_etaSIM
Definition: MuonTrackValidatorBase.h:379
TrackingParticle::trackPSimHit
const std::vector< PSimHit > & trackPSimHit() const
Definition: TrackingParticle.h:102
MuonTrackValidatorBase::h_fakes
std::vector< MonitorElement * > h_fakes
Definition: MuonTrackValidatorBase.h:382
MuonTrackValidator::ptres_vs_eta
std::vector< MonitorElement * > ptres_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidator::cotThetares_vs_eta
std::vector< MonitorElement * > cotThetares_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidatorBase::getEta
virtual double getEta(double eta)
Definition: MuonTrackValidatorBase.h:120
MuonTrackValidatorBase::h_eta
std::vector< MonitorElement * > h_eta
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidator::phimean_vs_eta_phi
std::vector< MonitorElement * > phimean_vs_eta_phi
Definition: MuonTrackValidator.h:162
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:424
edm::AssociationMap::find
const_iterator find(const key_type &k) const
find element with specified reference key
Definition: AssociationMap.h:94
MuonTrackValidatorBase::etaintervals
std::vector< std::vector< double > > etaintervals
Definition: MuonTrackValidatorBase.h:412
MuonTrackValidatorBase::h_pt
std::vector< MonitorElement * > h_pt
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidator::phires_vs_pt
std::vector< MonitorElement * > phires_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::h_pullQoverp
std::vector< MonitorElement * > h_pullQoverp
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidator::h_assochi2_prob
std::vector< MonitorElement * > h_assochi2_prob
Definition: MuonTrackValidator.h:165
abs
#define abs(x)
Definition: mlp_lapack.h:159
MuonTrackValidatorBase::totSIM_zpos
std::vector< std::vector< int > > totSIM_zpos
Definition: MuonTrackValidatorBase.h:425
MuonTrackValidatorBase::totASS2_dz
std::vector< std::vector< int > > totASS2_dz
Definition: MuonTrackValidatorBase.h:424
MuonTrackValidator::cotThetares_vs_pt
std::vector< MonitorElement * > cotThetares_vs_pt
Definition: MuonTrackValidator.h:175
min
#define min(a, b)
Definition: mlp_lapack.h:161
ParticleBase::Vector
math::XYZVectorD Vector
point in the space
Definition: ParticleBase.h:31
MuonTrackValidator::dxyres_vs_eta
std::vector< MonitorElement * > dxyres_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidatorBase::totASS_phi_Quality05
std::vector< std::vector< int > > totASS_phi_Quality05
Definition: MuonTrackValidatorBase.h:431
reco::TrackBase::ParameterVector
math::Vector< dimension >::type ParameterVector
parameter vector
Definition: TrackBase.h:69
MuonTrackValidator::dzres_vs_pt
std::vector< MonitorElement * > dzres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::h_tracks
std::vector< MonitorElement * > h_tracks
Definition: MuonTrackValidatorBase.h:382
MuonTrackValidator::dzres_vs_eta
std::vector< MonitorElement * > dzres_vs_eta
Definition: MuonTrackValidator.h:174
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:1124
MuonTrackValidatorBase::dxyintervals
std::vector< std::vector< double > > dxyintervals
Definition: MuonTrackValidatorBase.h:415
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:21
MuonTrackValidatorBase::nDThits_vs_eta
std::vector< MonitorElement * > nDThits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidatorBase::totASS_phi
std::vector< std::vector< int > > totASS_phi
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidator::thetapull_vs_phi
std::vector< MonitorElement * > thetapull_vs_phi
Definition: MuonTrackValidator.h:188
MuonTrackValidatorBase::totASS_dxy
std::vector< std::vector< int > > totASS_dxy
Definition: MuonTrackValidatorBase.h:423
MuonTrackValidator::h_assochi2
std::vector< MonitorElement * > h_assochi2
Definition: MuonTrackValidator.h:165
MuonTrackValidator::chi2_vs_phi
std::vector< MonitorElement * > chi2_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::totSIM_dxy
std::vector< std::vector< int > > totSIM_dxy
Definition: MuonTrackValidatorBase.h:423
MuonTrackValidator::chi2_vs_nhits
std::vector< MonitorElement * > chi2_vs_nhits
Definition: MuonTrackValidator.h:157
MuonTrackValidatorBase::totASSeta_Quality075
std::vector< std::vector< int > > totASSeta_Quality075
Definition: MuonTrackValidatorBase.h:429
MuonTrackValidatorBase::nhits_vs_eta
std::vector< MonitorElement * > nhits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidator::phires_vs_phi
std::vector< MonitorElement * > phires_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::getPt
virtual double getPt(double pt)
Definition: MuonTrackValidatorBase.h:125
MuonTrackValidator::cosmictpSelector
CosmicTrackingParticleSelector cosmictpSelector
Definition: MuonTrackValidator.h:146
ParticleBase::Point
math::XYZPointD Point
point in the space
Definition: ParticleBase.h:29
MuonTrackValidator::nhits_vs_phi
std::vector< MonitorElement * > nhits_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::totASS2pT
std::vector< std::vector< int > > totASS2pT
Definition: MuonTrackValidatorBase.h:420
MuonTrackValidatorBase::h_vertposSIM
std::vector< MonitorElement * > h_vertposSIM
Definition: MuonTrackValidatorBase.h:379
ExpressReco_HICollisions_FallBack.Quality
tuple Quality
Definition: ExpressReco_HICollisions_FallBack.py:11936
MuonTrackValidator::dxyres_vs_pt
std::vector< MonitorElement * > dxyres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::label
std::vector< edm::InputTag > label
Definition: MuonTrackValidatorBase.h:338
MuonTrackValidatorBase::h_assocSharedHit
std::vector< MonitorElement * > h_assocSharedHit
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidator::h_nchi2
std::vector< MonitorElement * > h_nchi2
Definition: MuonTrackValidator.h:154
MuonTrackValidatorBase::pTintervals
std::vector< std::vector< double > > pTintervals
Definition: MuonTrackValidatorBase.h:413
dt_offlineAnalysis_common_cff.reco
tuple reco
Definition: dt_offlineAnalysis_common_cff.py:54
edm::AssociationMap::post_insert
void post_insert()
post insert action
Definition: AssociationMap.h:135
MuonTrackValidatorBase::totASS2eta
std::vector< std::vector< int > > totASS2eta
Definition: MuonTrackValidatorBase.h:419
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
MuonTrackValidatorBase::totSIMpT
std::vector< std::vector< int > > totSIMpT
Definition: MuonTrackValidatorBase.h:420
MuonTrackValidatorBase::totREC_dxy
std::vector< std::vector< int > > totREC_dxy
Definition: MuonTrackValidatorBase.h:423
MuonTrackValidatorBase::h_pullPhi
std::vector< MonitorElement * > h_pullPhi
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidatorBase::totSIM_hit
std::vector< std::vector< int > > totSIM_hit
Definition: MuonTrackValidatorBase.h:421
MuonTrackValidator::h_nchi2_prob
std::vector< MonitorElement * > h_nchi2_prob
Definition: MuonTrackValidator.h:154
MuonTrackValidator::ptpull_vs_eta
std::vector< MonitorElement * > ptpull_vs_eta
Definition: MuonTrackValidator.h:187
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
MuonTrackValidatorBase::totASSpT
std::vector< std::vector< int > > totASSpT
Definition: MuonTrackValidatorBase.h:420
MuonTrackValidator::dxypull_vs_eta
std::vector< MonitorElement * > dxypull_vs_eta
Definition: MuonTrackValidator.h:187
edm::View::size_type
unsigned int size_type
Definition: View.h:90
MuonTrackValidatorBase::totASS_phi_Quality075
std::vector< std::vector< int > > totASS_phi_Quality075
Definition: MuonTrackValidatorBase.h:431
f
double f[11][100]
Definition: MuScleFitUtils.cc:79
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:225
MuonTrackValidatorBase::h_pullTheta
std::vector< MonitorElement * > h_pullTheta
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidatorBase::totSIM_vertpos
std::vector< std::vector< int > > totSIM_vertpos
Definition: MuonTrackValidatorBase.h:425
MuonTrackValidatorBase::totASSpT_Quality05
std::vector< std::vector< int > > totASSpT_Quality05
Definition: MuonTrackValidatorBase.h:430
MuonTrackValidatorBase::totREC_hit
std::vector< std::vector< int > > totREC_hit
Definition: MuonTrackValidatorBase.h:421
MuonTrackValidator::phipull_vs_eta
std::vector< MonitorElement * > phipull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidatorBase::totSIMeta
std::vector< std::vector< int > > totSIMeta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::totASS2_dxy
std::vector< std::vector< int > > totASS2_dxy
Definition: MuonTrackValidatorBase.h:423
MuonTrackValidator::thetapull_vs_eta
std::vector< MonitorElement * > thetapull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidatorBase::h_charge
std::vector< MonitorElement * > h_charge
Definition: MuonTrackValidatorBase.h:382
edm::RefToBase::key
size_t key() const
Definition: RefToBase.h:223
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:359
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:511
MuonTrackValidator::ptmean_vs_eta_phi
std::vector< MonitorElement * > ptmean_vs_eta_phi
Definition: MuonTrackValidator.h:162
MuonTrackValidatorBase::totASS2_phi
std::vector< std::vector< int > > totASS2_phi
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::phiintervals
std::vector< std::vector< double > > phiintervals
Definition: MuonTrackValidatorBase.h:414
MuonTrackValidatorBase::totASSpT_Quality075
std::vector< std::vector< int > > totASSpT_Quality075
Definition: MuonTrackValidatorBase.h:430
M_PI
#define M_PI
Definition: BFit3D.cc:3
MuonTrackValidatorBase::h_pullDxy
std::vector< MonitorElement * > h_pullDxy
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidatorBase::h_tracksSIM
std::vector< MonitorElement * > h_tracksSIM
Definition: MuonTrackValidatorBase.h:379
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::Ref::key
key_type key() const
Accessor for product key.
Definition: Ref.h:264
MuonTrackValidator::chi2_vs_eta
std::vector< MonitorElement * > chi2_vs_eta
Definition: MuonTrackValidator.h:168
MuonTrackValidator::tpSelector
TrackingParticleSelector tpSelector
Definition: MuonTrackValidator.h:145
ParticleBase::momentum
Vector momentum() const
spatial momentum vector
Definition: ParticleBase.h:87
edm::Handle::product
T const * product() const
Definition: Handle.h:74
MuonTrackValidatorBase::totSIM_phi
std::vector< std::vector< int > > totSIM_phi
Definition: MuonTrackValidatorBase.h:422
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:25
MuonTrackValidatorBase::h_hits
std::vector< MonitorElement * > h_hits
Definition: MuonTrackValidatorBase.h:382
MuonTrackValidator::h_losthits
std::vector< MonitorElement * > h_losthits
Definition: MuonTrackValidator.h:154
MuonTrackValidator::ptpull_vs_phi
std::vector< MonitorElement * > ptpull_vs_phi
Definition: MuonTrackValidator.h:188
edm::InputTag::process
std::string const & process() const
Definition: InputTag.h:29
tmp
std::vector< std::vector< double > > tmp
Definition: MVATrainer.cc:100
MuonTrackValidator::associatormap
edm::InputTag associatormap
Definition: MuonTrackValidator.h:138
MuonTrackValidatorBase::totASSeta
std::vector< std::vector< int > > totASSeta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidator::ptres_vs_pt
std::vector< MonitorElement * > ptres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::associator
std::vector< const TrackAssociatorBase * > associator
Definition: MuonTrackValidatorBase.h:376
MuonTrackValidatorBase::totRECeta
std::vector< std::vector< int > > totRECeta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::vertposintervals
std::vector< std::vector< double > > vertposintervals
Definition: MuonTrackValidatorBase.h:417
MuonTrackValidatorBase::dzintervals
std::vector< std::vector< double > > dzintervals
Definition: MuonTrackValidatorBase.h:416
MuonTrackValidator::phires_vs_eta
std::vector< MonitorElement * > phires_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidatorBase::h_ptSIM
std::vector< MonitorElement * > h_ptSIM
Definition: MuonTrackValidatorBase.h:379
MuonTrackValidatorBase::zposintervals
std::vector< std::vector< double > > zposintervals
Definition: MuonTrackValidatorBase.h:418
MuonTrackValidatorBase::h_assocFraction
std::vector< MonitorElement * > h_assocFraction
Definition: MuonTrackValidatorBase.h:402
TrackingParticleCollection
std::vector< TrackingParticle > TrackingParticleCollection
Definition: TrackingParticleFwd.h:8
MuonTrackValidatorBase::totSIM_dz
std::vector< std::vector< int > > totSIM_dz
Definition: MuonTrackValidatorBase.h:424
MuonTrackValidatorBase::associators
std::vector< std::string > associators
Definition: MuonTrackValidatorBase.h:344
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:239
ParticleBase::vertex
const Point & vertex() const
vertex position
Definition: ParticleBase.h:229
MuonTrackValidator::nlosthits_vs_eta
std::vector< MonitorElement * > nlosthits_vs_eta
Definition: MuonTrackValidator.h:168
MuonTrackValidator::dzpull_vs_eta
std::vector< MonitorElement * > dzpull_vs_eta
Definition: MuonTrackValidator.h:187
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:26
MuonTrackValidatorBase::nRPChits_vs_eta
std::vector< MonitorElement * > nRPChits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidatorBase::totASS_vertpos
std::vector< std::vector< int > > totASS_vertpos
Definition: MuonTrackValidatorBase.h:425
MuonTrackValidatorBase::totASS_hit
std::vector< std::vector< int > > totASS_hit
Definition: MuonTrackValidatorBase.h:421
MuonTrackValidatorBase::nCSChits_vs_eta
std::vector< MonitorElement * > nCSChits_vs_eta
Definition: MuonTrackValidatorBase.h:405
asciidump.at
list at
Definition: asciidump.py:428
MuonTrackValidator::phipull_vs_phi
std::vector< MonitorElement * > phipull_vs_phi
Definition: MuonTrackValidator.h:188
void MuonTrackValidator::beginRun ( edm::Run const &  ,
edm::EventSetup const &  setup 
)
virtual

Method called before the event loop.

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

Reimplemented from edm::EDAnalyzer.

Definition at line 29 of file MuonTrackValidator.cc.

References ExpressReco_HICollisions_FallBack::algo, DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), DQMStore::bookProfile2D(), DQMStore::cd(), dbe_, TrackerOfflineValidation_Dqm_cff::dirName, edm::EventSetup::get(), getTH1F(), getTH2F(), DQMStore::goUp(), edm::InputTag::instance(), j, edm::InputTag::label(), label, max(), min, nint(), edm::InputTag::process(), edm::ESHandle< class >::product(), linker::replace(), and DQMStore::setCurrentFolder().

29  {
30 
31  // dbe_->showDirStructure();
32 
33  int j=0;
34  for (unsigned int ww=0;ww<associators.size();ww++){
35  for (unsigned int www=0;www<label.size();www++){
36 
37  dbe_->cd();
38  InputTag algo = label[www];
39  string dirName=dirName_;
40  if (algo.process()!="")
41  dirName+=algo.process()+"_";
42  if(algo.label()!="")
43  dirName+=algo.label()+"_";
44  if(algo.instance()!="")
45  dirName+=algo.instance()+"_";
46  if (dirName.find("Tracks")<dirName.length()){
47  dirName.replace(dirName.find("Tracks"),6,"");
48  }
49  string assoc= associators[ww];
50  if (assoc.find("Track")<assoc.length()){
51  assoc.replace(assoc.find("Track"),5,"");
52  }
53  dirName+=assoc;
54  std::replace(dirName.begin(), dirName.end(), ':', '_');
55  dbe_->setCurrentFolder(dirName.c_str());
56 
57  setUpVectors();
58 
59  dbe_->goUp();
60  string subDirName = dirName + "/simulation";
61  dbe_->setCurrentFolder(subDirName.c_str());
62  h_ptSIM.push_back( dbe_->book1D("ptSIM", "generated p_{t}", 5500, 0, 110 ) );
63  h_etaSIM.push_back( dbe_->book1D("etaSIM", "generated pseudorapidity", 500, -2.5, 2.5 ) );
64  h_tracksSIM.push_back( dbe_->book1D("tracksSIM","number of simulated tracks",200,-0.5,99.5) );
65  h_vertposSIM.push_back( dbe_->book1D("vertposSIM","Transverse position of sim vertices",100,0.,120.) );
66 
67  dbe_->cd();
68  dbe_->setCurrentFolder(dirName.c_str());
69  h_tracks.push_back( dbe_->book1D("tracks","number of reconstructed tracks",200,-0.5,19.5) );
70  h_fakes.push_back( dbe_->book1D("fakes","number of fake reco tracks",20,-0.5,19.5) );
71  h_charge.push_back( dbe_->book1D("charge","charge",3,-1.5,1.5) );
72  h_hits.push_back( dbe_->book1D("hits", "number of hits per track", nintHit,minHit,maxHit ) );
73  h_losthits.push_back( dbe_->book1D("losthits", "number of lost hits per track", nintHit,minHit,maxHit) );
74  h_nchi2.push_back( dbe_->book1D("chi2", "normalized #chi^{2}", 200, 0, 20 ) );
75  h_nchi2_prob.push_back( dbe_->book1D("chi2_prob", "normalized #chi^{2} probability",100,0,1));
76 
77  h_effic.push_back( dbe_->book1D("effic","efficiency vs #eta",nint,min,max) );
78  h_efficPt.push_back( dbe_->book1D("efficPt","efficiency vs pT",nintpT,minpT,maxpT) );
79  h_effic_vs_hit.push_back( dbe_->book1D("effic_vs_hit","effic vs hit",nintHit,minHit,maxHit) );
80  h_effic_vs_phi.push_back( dbe_->book1D("effic_vs_phi","effic vs phi",nintPhi,minPhi,maxPhi) );
81  h_effic_vs_dxy.push_back( dbe_->book1D("effic_vs_dxy","effic vs dxy",nintDxy,minDxy,maxDxy) );
82  h_effic_vs_dz.push_back( dbe_->book1D("effic_vs_dz","effic vs dz",nintDz,minDz,maxDz) );
83  h_effic_vs_vertpos.push_back( dbe_->book1D("effic_vs_vertpos","effic vs vertpos",nintVertpos,minVertpos,maxVertpos) );
84  h_effic_vs_zpos.push_back( dbe_->book1D("effic_vs_zpos","effic vs zpos",nintZpos,minZpos,maxZpos) );
85 
86  h_fakerate.push_back( dbe_->book1D("fakerate","fake rate vs #eta",nint,min,max) );
87  h_fakeratePt.push_back( dbe_->book1D("fakeratePt","fake rate vs pT",nintpT,minpT,maxpT) );
88  h_fake_vs_hit.push_back( dbe_->book1D("fakerate_vs_hit","fake rate vs hit",nintHit,minHit,maxHit) );
89  h_fake_vs_phi.push_back( dbe_->book1D("fakerate_vs_phi","fake vs phi",nintPhi,minPhi,maxPhi) );
90  h_fake_vs_dxy.push_back( dbe_->book1D("fakerate_vs_dxy","fake rate vs dxy",nintDxy,minDxy,maxDxy) );
91  h_fake_vs_dz.push_back( dbe_->book1D("fakerate_vs_dz","fake vs dz",nintDz,minDz,maxDz) );
92 
93 
95  h_recoeta.push_back( dbe_->book1D("num_reco_eta","N of reco track vs eta",nint,min,max) );
96  h_assoceta.push_back( dbe_->book1D("num_assoc(simToReco)_eta","N of associated tracks (simToReco) vs eta",nint,min,max) );
97  h_assoc2eta.push_back( dbe_->book1D("num_assoc(recoToSim)_eta","N of associated (recoToSim) tracks vs eta",nint,min,max) );
98  h_simuleta.push_back( dbe_->book1D("num_simul_eta","N of simulated tracks vs eta",nint,min,max) );
99  h_recopT.push_back( dbe_->book1D("num_reco_pT","N of reco track vs pT",nintpT,minpT,maxpT) );
100  h_assocpT.push_back( dbe_->book1D("num_assoc(simToReco)_pT","N of associated tracks (simToReco) vs pT",nintpT,minpT,maxpT) );
101  h_assoc2pT.push_back( dbe_->book1D("num_assoc(recoToSim)_pT","N of associated (recoToSim) tracks vs pT",nintpT,minpT,maxpT) );
102  h_simulpT.push_back( dbe_->book1D("num_simul_pT","N of simulated tracks vs pT",nintpT,minpT,maxpT) );
103  //
104  h_recohit.push_back( dbe_->book1D("num_reco_hit","N of reco track vs hit",nintHit,minHit,maxHit) );
105  h_assochit.push_back( dbe_->book1D("num_assoc(simToReco)_hit","N of associated tracks (simToReco) vs hit",nintHit,minHit,maxHit) );
106  h_assoc2hit.push_back( dbe_->book1D("num_assoc(recoToSim)_hit","N of associated (recoToSim) tracks vs hit",nintHit,minHit,maxHit) );
107  h_simulhit.push_back( dbe_->book1D("num_simul_hit","N of simulated tracks vs hit",nintHit,minHit,maxHit) );
108  //
109  h_recophi.push_back( dbe_->book1D("num_reco_phi","N of reco track vs phi",nintPhi,minPhi,maxPhi) );
110  h_assocphi.push_back( dbe_->book1D("num_assoc(simToReco)_phi","N of associated tracks (simToReco) vs phi",nintPhi,minPhi,maxPhi) );
111  h_assoc2phi.push_back( dbe_->book1D("num_assoc(recoToSim)_phi","N of associated (recoToSim) tracks vs phi",nintPhi,minPhi,maxPhi) );
112  h_simulphi.push_back( dbe_->book1D("num_simul_phi","N of simulated tracks vs phi",nintPhi,minPhi,maxPhi) );
113 
114  h_recodxy.push_back( dbe_->book1D("num_reco_dxy","N of reco track vs dxy",nintDxy,minDxy,maxDxy) );
115  h_assocdxy.push_back( dbe_->book1D("num_assoc(simToReco)_dxy","N of associated tracks (simToReco) vs dxy",nintDxy,minDxy,maxDxy) );
116  h_assoc2dxy.push_back( dbe_->book1D("num_assoc(recoToSim)_dxy","N of associated (recoToSim) tracks vs dxy",nintDxy,minDxy,maxDxy) );
117  h_simuldxy.push_back( dbe_->book1D("num_simul_dxy","N of simulated tracks vs dxy",nintDxy,minDxy,maxDxy) );
118 
119  h_recodz.push_back( dbe_->book1D("num_reco_dz","N of reco track vs dz",nintDz,minDz,maxDz) );
120  h_assocdz.push_back( dbe_->book1D("num_assoc(simToReco)_dz","N of associated tracks (simToReco) vs dz",nintDz,minDz,maxDz) );
121  h_assoc2dz.push_back( dbe_->book1D("num_assoc(recoToSim)_dz","N of associated (recoToSim) tracks vs dz",nintDz,minDz,maxDz) );
122  h_simuldz.push_back( dbe_->book1D("num_simul_dz","N of simulated tracks vs dz",nintDz,minDz,maxDz) );
123 
124  h_assocvertpos.push_back( dbe_->book1D("num_assoc(simToReco)_vertpos","N of associated tracks (simToReco) vs transverse vert position",nintVertpos,minVertpos,maxVertpos) );
125  h_simulvertpos.push_back( dbe_->book1D("num_simul_vertpos","N of simulated tracks vs transverse vert position",nintVertpos,minVertpos,maxVertpos) );
126 
127  h_assoczpos.push_back( dbe_->book1D("num_assoc(simToReco)_zpos","N of associated tracks (simToReco) vs z vert position",nintZpos,minZpos,maxZpos) );
128  h_simulzpos.push_back( dbe_->book1D("num_simul_zpos","N of simulated tracks vs z vert position",nintZpos,minZpos,maxZpos) );
129 
130 
132 
133  h_eta.push_back( dbe_->book1D("eta", "pseudorapidity residue", 1000, -0.1, 0.1 ) );
134  h_pt.push_back( dbe_->book1D("pullPt", "pull of p_{t}", 100, -10, 10 ) );
135  h_pullTheta.push_back( dbe_->book1D("pullTheta","pull of #theta parameter",250,-25,25) );
136  h_pullPhi.push_back( dbe_->book1D("pullPhi","pull of #phi parameter",250,-25,25) );
137  h_pullDxy.push_back( dbe_->book1D("pullDxy","pull of dxy parameter",250,-25,25) );
138  h_pullDz.push_back( dbe_->book1D("pullDz","pull of dz parameter",250,-25,25) );
139  h_pullQoverp.push_back( dbe_->book1D("pullQoverp","pull of qoverp parameter",250,-25,25) );
140 
141  if (associators[ww]=="TrackAssociatorByChi2"){
142  h_assochi2.push_back( dbe_->book1D("assocChi2","track association #chi^{2}",1000000,0,100000) );
143  h_assochi2_prob.push_back(dbe_->book1D("assocChi2_prob","probability of association #chi^{2}",100,0,1));
144  } else if (associators[ww]=="TrackAssociatorByHits"){
145  h_assocFraction.push_back( dbe_->book1D("assocFraction","fraction of shared hits",200,0,2) );
146  h_assocSharedHit.push_back(dbe_->book1D("assocSharedHit","number of shared hits",20,0,20));
147  }
148 
149  chi2_vs_nhits.push_back( dbe_->book2D("chi2_vs_nhits","#chi^{2} vs nhits",25,0,25,100,0,10) );
150  h_chi2meanhitsh.push_back( dbe_->bookProfile("chi2mean_vs_nhits","mean #chi^{2} vs nhits",25,0,25,100,0,10) );
151 
152  etares_vs_eta.push_back( dbe_->book2D("etares_vs_eta","etaresidue vs eta",nint,min,max,200,-0.1,0.1) );
153  nrec_vs_nsim.push_back( dbe_->book2D("nrec_vs_nsim","nrec vs nsim",20,-0.5,19.5,20,-0.5,19.5) );
154 
155  chi2_vs_eta.push_back( dbe_->book2D("chi2_vs_eta","chi2_vs_eta",nint,min,max, 200, 0, 20 ));
156  h_chi2meanh.push_back( dbe_->bookProfile("chi2mean","mean #chi^{2} vs #eta",nint,min,max, 200, 0, 20) );
157  chi2_vs_phi.push_back( dbe_->book2D("chi2_vs_phi","#chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 200, 0, 20 ) );
158  h_chi2mean_vs_phi.push_back( dbe_->bookProfile("chi2mean_vs_phi","mean of #chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 200, 0, 20) );
159 
160  nhits_vs_eta.push_back( dbe_->book2D("nhits_vs_eta","nhits vs eta",nint,min,max,nintHit,minHit,maxHit) );
161  nDThits_vs_eta.push_back( dbe_->book2D("nDThits_vs_eta","# DT hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
162  nCSChits_vs_eta.push_back( dbe_->book2D("nCSChits_vs_eta","# CSC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
163  nRPChits_vs_eta.push_back( dbe_->book2D("nRPChits_vs_eta","# RPC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
164 
165  h_DThits_eta.push_back( dbe_->bookProfile("DThits_eta","mean # DT hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
166  h_CSChits_eta.push_back( dbe_->bookProfile("CSChits_eta","mean # CSC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
167  h_RPChits_eta.push_back( dbe_->bookProfile("RPChits_eta","mean # RPC hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
168  h_hits_eta.push_back( dbe_->bookProfile("hits_eta","mean #hits vs eta",nint,min,max,nintHit,minHit,maxHit) );
169  nhits_vs_phi.push_back( dbe_->book2D("nhits_vs_phi","#hits vs #phi",nintPhi,minPhi,maxPhi,nintHit,minHit,maxHit) );
170  h_hits_phi.push_back( dbe_->bookProfile("hits_phi","mean #hits vs #phi",nintPhi,minPhi,maxPhi, nintHit,minHit,maxHit) );
171 
172  nlosthits_vs_eta.push_back( dbe_->book2D("nlosthits_vs_eta","nlosthits vs eta",nint,min,max,nintHit,minHit,maxHit) );
173  h_losthits_eta.push_back( dbe_->bookProfile("losthits_eta","losthits_eta",nint,min,max,nintHit,minHit,maxHit) );
174 
175  ptres_vs_eta.push_back(dbe_->book2D("ptres_vs_eta","ptres_vs_eta",nint,min,max, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
176  h_ptrmsh.push_back( dbe_->book1D("ptres_vs_eta_Sigma","#sigma(#deltap_{t}/p_{t}) vs #eta",nint,min,max) );
177 
178  ptres_vs_phi.push_back( dbe_->book2D("ptres_vs_phi","p_{t} res vs #phi",nintPhi,minPhi,maxPhi, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
179  h_ptmeanhPhi.push_back( dbe_->book1D("ptres_vs_phi_Mean","mean of p_{t} resolution vs #phi",nintPhi,minPhi,maxPhi));
180  h_ptrmshPhi.push_back( dbe_->book1D("ptres_vs_phi_Sigma","#sigma(#deltap_{t}/p_{t}) vs #phi",nintPhi,minPhi,maxPhi) );
181 
182  ptres_vs_pt.push_back(dbe_->book2D("ptres_vs_pt","ptres_vs_pt",nintpT,minpT,maxpT, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
183  h_ptmeanhPt.push_back( dbe_->book1D("ptres_vs_pt_Mean","mean of p_{t} resolution vs p_{t}",nintpT,minpT,maxpT));
184  h_ptrmshPt.push_back( dbe_->book1D("ptres_vs_pt_Sigma","#sigma(#deltap_{t}/p_{t}) vs pT",nintpT,minpT,maxpT) );
185 
186  cotThetares_vs_eta.push_back(dbe_->book2D("cotThetares_vs_eta","cotThetares_vs_eta",nint,min,max,cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
187  h_cotThetameanh.push_back( dbe_->book1D("cotThetares_vs_eta_Mean","#sigma(cot(#theta)) vs #eta Mean",nint,min,max) );
188  h_cotThetarmsh.push_back( dbe_->book1D("cotThetares_vs_eta_Sigma","#sigma(cot(#theta)) vs #eta Sigma",nint,min,max) );
189 
190 
191  cotThetares_vs_pt.push_back(dbe_->book2D("cotThetares_vs_pt","cotThetares_vs_pt",nintpT,minpT,maxpT, cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
192  h_cotThetameanhPt.push_back( dbe_->book1D("cotThetares_vs_pt_Mean","#sigma(cot(#theta)) vs pT Mean",nintpT,minpT,maxpT) );
193  h_cotThetarmshPt.push_back( dbe_->book1D("cotThetares_vs_pt_Sigma","#sigma(cot(#theta)) vs pT Sigma",nintpT,minpT,maxpT) );
194 
195 
196  phires_vs_eta.push_back(dbe_->book2D("phires_vs_eta","phires_vs_eta",nint,min,max, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
197  h_phimeanh.push_back(dbe_->book1D("phires_vs_eta_Mean","mean of #phi res vs #eta",nint,min,max));
198  h_phirmsh.push_back( dbe_->book1D("phires_vs_eta_Sigma","#sigma(#delta#phi) vs #eta",nint,min,max) );
199 
200  phires_vs_pt.push_back(dbe_->book2D("phires_vs_pt","phires_vs_pt",nintpT,minpT,maxpT, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
201  h_phimeanhPt.push_back(dbe_->book1D("phires_vs_pt_Mean","mean of #phi res vs pT",nintpT,minpT,maxpT));
202  h_phirmshPt.push_back( dbe_->book1D("phires_vs_pt_Sigma","#sigma(#delta#phi) vs pT",nintpT,minpT,maxpT) );
203 
204  phires_vs_phi.push_back(dbe_->book2D("phires_vs_phi","#phi res vs #phi",nintPhi,minPhi,maxPhi,phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
205  h_phimeanhPhi.push_back(dbe_->book1D("phires_vs_phi_Mean","mean of #phi res vs #phi",nintPhi,minPhi,maxPhi));
206  h_phirmshPhi.push_back( dbe_->book1D("phires_vs_phi_Sigma","#sigma(#delta#phi) vs #phi",nintPhi,minPhi,maxPhi) );
207 
208  dxyres_vs_eta.push_back(dbe_->book2D("dxyres_vs_eta","dxyres_vs_eta",nint,min,max,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
209  h_dxymeanh.push_back( dbe_->book1D("dxyres_vs_eta_Mean","mean of dxyres vs #eta",nint,min,max) );
210  h_dxyrmsh.push_back( dbe_->book1D("dxyres_vs_eta_Sigma","#sigma(#deltadxy) vs #eta",nint,min,max) );
211 
212  dxyres_vs_pt.push_back( dbe_->book2D("dxyres_vs_pt","dxyres_vs_pt",nintpT,minpT,maxpT,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
213  h_dxymeanhPt.push_back( dbe_->book1D("dxyres_vs_pt_Mean","mean of dxyres vs pT",nintpT,minpT,maxpT) );
214  h_dxyrmshPt.push_back( dbe_->book1D("dxyres_vs_pt_Sigma","#sigmadxy vs pT",nintpT,minpT,maxpT) );
215 
216  dzres_vs_eta.push_back(dbe_->book2D("dzres_vs_eta","dzres_vs_eta",nint,min,max,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
217  h_dzmeanh.push_back( dbe_->book1D("dzres_vs_eta_Mean","mean of dzres vs #eta",nint,min,max) );
218  h_dzrmsh.push_back( dbe_->book1D("dzres_vs_eta_Sigma","#sigma(#deltadz) vs #eta",nint,min,max) );
219 
220  dzres_vs_pt.push_back(dbe_->book2D("dzres_vs_pt","dzres_vs_pt",nintpT,minpT,maxpT,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
221  h_dzmeanhPt.push_back( dbe_->book1D("dzres_vs_pt_Mean","mean of dzres vs pT",nintpT,minpT,maxpT) );
222  h_dzrmshPt.push_back( dbe_->book1D("dzres_vs_pt_Sigma","#sigma(#deltadz vs pT",nintpT,minpT,maxpT) );
223 
224  ptmean_vs_eta_phi.push_back(dbe_->bookProfile2D("ptmean_vs_eta_phi","mean p_{t} vs #eta and #phi",nintPhi,minPhi,maxPhi,nint,min,max,1000,0,1000));
225  phimean_vs_eta_phi.push_back(dbe_->bookProfile2D("phimean_vs_eta_phi","mean #phi vs #eta and #phi",nintPhi,minPhi,maxPhi,nint,min,max,nintPhi,minPhi,maxPhi));
226 
227  //pulls of track params vs eta: to be used with fitslicesytool
228  dxypull_vs_eta.push_back(dbe_->book2D("dxypull_vs_eta","dxypull_vs_eta",nint,min,max,100,-10,10));
229  ptpull_vs_eta.push_back(dbe_->book2D("ptpull_vs_eta","ptpull_vs_eta",nint,min,max,100,-10,10));
230  dzpull_vs_eta.push_back(dbe_->book2D("dzpull_vs_eta","dzpull_vs_eta",nint,min,max,100,-10,10));
231  phipull_vs_eta.push_back(dbe_->book2D("phipull_vs_eta","phipull_vs_eta",nint,min,max,100,-10,10));
232  thetapull_vs_eta.push_back(dbe_->book2D("thetapull_vs_eta","thetapull_vs_eta",nint,min,max,100,-10,10));
233 
234  h_dxypulletamean.push_back( dbe_->book1D("h_dxypulleta_Mean","mean of dxy pull vs #eta",nint,min,max) );
235  h_ptpulletamean.push_back( dbe_->book1D("h_ptpulleta_Mean","mean of p_{t} pull vs #eta",nint,min,max) );
236  h_dzpulletamean.push_back( dbe_->book1D("h_dzpulleta_Mean","mean of dz pull vs #eta",nint,min,max) );
237  h_phipulletamean.push_back( dbe_->book1D("h_phipulleta_Mean","mean of #phi pull vs #eta",nint,min,max) );
238  h_thetapulletamean.push_back( dbe_->book1D("h_thetapulleta_Mean","mean of #theta pull vs #eta",nint,min,max) );
239  // h_ptshiftetamean.push_back( dbe_->book1D("h_ptshifteta_Mean","<#deltapT/pT>[%] vs #eta",nint,min,max) );
240 
241  h_dxypulleta.push_back( dbe_->book1D("h_dxypulleta_Sigma","#sigma of dxy pull vs #eta",nint,min,max) );
242  h_ptpulleta.push_back( dbe_->book1D("h_ptpulleta_Sigma","#sigma of p_{t} pull vs #eta",nint,min,max) );
243  h_dzpulleta.push_back( dbe_->book1D("h_dzpulleta_Sigma","#sigma of dz pull vs #eta",nint,min,max) );
244  h_phipulleta.push_back( dbe_->book1D("h_phipulleta_Sigma","#sigma of #phi pull vs #eta",nint,min,max) );
245  h_thetapulleta.push_back( dbe_->book1D("h_thetapulleta_Sigma","#sigma of #theta pull vs #eta",nint,min,max) );
246  h_ptshifteta.push_back( dbe_->book1D("ptres_vs_eta_Mean","<#deltapT/pT>[%] vs #eta",nint,min,max) );
247 
248  //pulls of track params vs phi
249  ptpull_vs_phi.push_back(dbe_->book2D("ptpull_vs_phi","p_{t} pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
250  phipull_vs_phi.push_back(dbe_->book2D("phipull_vs_phi","#phi pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
251  thetapull_vs_phi.push_back(dbe_->book2D("thetapull_vs_phi","#theta pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
252 
253  h_ptpullphimean.push_back( dbe_->book1D("h_ptpullphi_Mean","mean of p_{t} pull vs #phi",nintPhi,minPhi,maxPhi) );
254  h_phipullphimean.push_back( dbe_->book1D("h_phipullphi_Mean","mean of #phi pull vs #phi",nintPhi,minPhi,maxPhi) );
255  h_thetapullphimean.push_back( dbe_->book1D("h_thetapullphi_Mean","mean of #theta pull vs #phi",nintPhi,minPhi,maxPhi) );
256 
257  h_ptpullphi.push_back( dbe_->book1D("h_ptpullphi_Sigma","#sigma of p_{t} pull vs #phi",nintPhi,minPhi,maxPhi) );
258  h_phipullphi.push_back( dbe_->book1D("h_phipullphi_Sigma","#sigma of #phi pull vs #phi",nintPhi,minPhi,maxPhi) );
259  h_thetapullphi.push_back( dbe_->book1D("h_thetapullphi_Sigma","#sigma of #theta pull vs #phi",nintPhi,minPhi,maxPhi) );
260 
261 
262  nrecHit_vs_nsimHit_sim2rec.push_back( dbe_->book2D("nrecHit_vs_nsimHit_sim2rec","nrecHit vs nsimHit (Sim2RecAssoc)",nintHit,minHit,maxHit, nintHit,minHit,maxHit ));
263  nrecHit_vs_nsimHit_rec2sim.push_back( dbe_->book2D("nrecHit_vs_nsimHit_rec2sim","nrecHit vs nsimHit (Rec2simAssoc)",nintHit,minHit,maxHit, nintHit,minHit,maxHit ));
264 
265  if (MABH) {
266  h_PurityVsQuality.push_back( dbe_->book2D("PurityVsQuality","Purity vs Quality (MABH)",20,0.01,1.01,20,0.01,1.01) );
267  h_effic_Quality05.push_back( dbe_->book1D("effic_Q05","efficiency vs #eta (Quality>0.5)",nint,min,max) );
268  h_effic_Quality075.push_back( dbe_->book1D("effic_Q075","efficiency vs #eta (Quality>0.75)",nint,min,max) );
269  h_efficPt_Quality05.push_back( dbe_->book1D("efficPt_Q05","efficiency vs pT (Quality>0.5)",nintpT,minpT,maxpT) );
270  h_efficPt_Quality075.push_back( dbe_->book1D("efficPt_Q075","efficiency vs pT (Quality>0.75)",nintpT,minpT,maxpT) );
271  h_effic_vs_phi_Quality05.push_back( dbe_->book1D("effic_vs_phi_Q05","effic vs phi (Quality>0.5)",nintPhi,minPhi,maxPhi) );
272  h_effic_vs_phi_Quality075.push_back( dbe_->book1D("effic_vs_phi_Q075","effic vs phi (Quality>0.75)",nintPhi,minPhi,maxPhi) );
273  //
274  h_assoceta_Quality05.push_back( dbe_->book1D("num_assoc(simToReco)_eta_Q05","N of associated tracks (simToReco) vs eta (Quality>0.5)",nint,min,max) );
275  h_assoceta_Quality075.push_back( dbe_->book1D("num_assoc(simToReco)_eta_Q075","N of associated tracks (simToReco) vs eta (Quality>0.75)",nint,min,max) );
276  h_assocpT_Quality05.push_back( dbe_->book1D("num_assoc(simToReco)_pT_Q05","N of associated tracks (simToReco) vs pT (Quality>0.5)",nintpT,minpT,maxpT) );
277  h_assocpT_Quality075.push_back( dbe_->book1D("num_assoc(simToReco)_pT_Q075","N of associated tracks (simToReco) vs pT (Quality>0.75)",nintpT,minpT,maxpT) );
278  h_assocphi_Quality05.push_back( dbe_->book1D("num_assoc(simToReco)_phi_Q05","N of associated tracks (simToReco) vs phi (Quality>0.5)",nintPhi,minPhi,maxPhi) );
279  h_assocphi_Quality075.push_back( dbe_->book1D("num_assoc(simToReco)_phi_Q075","N of associated tracks (simToReco) vs phi (Quality>0.75)",nintPhi,minPhi,maxPhi) );
280  }
281 
282  if(useLogPt){
283  BinLogX(dzres_vs_pt[j]->getTH2F());
284  BinLogX(h_dzmeanhPt[j]->getTH1F());
285  BinLogX(h_dzrmshPt[j]->getTH1F());
286 
287  BinLogX(dxyres_vs_pt[j]->getTH2F());
288  BinLogX(h_dxymeanhPt[j]->getTH1F());
289  BinLogX(h_dxyrmshPt[j]->getTH1F());
290 
291  BinLogX(phires_vs_pt[j]->getTH2F());
292  BinLogX(h_phimeanhPt[j]->getTH1F());
293  BinLogX(h_phirmshPt[j]->getTH1F());
294 
295  BinLogX(cotThetares_vs_pt[j]->getTH2F());
296  BinLogX(h_cotThetameanhPt[j]->getTH1F());
297  BinLogX(h_cotThetarmshPt[j]->getTH1F());
298 
299  BinLogX(ptres_vs_pt[j]->getTH2F());
300  BinLogX(h_ptmeanhPt[j]->getTH1F());
301  BinLogX(h_ptrmshPt[j]->getTH1F());
302 
303  BinLogX(h_efficPt[j]->getTH1F());
304  BinLogX(h_fakeratePt[j]->getTH1F());
305  BinLogX(h_recopT[j]->getTH1F());
306  BinLogX(h_assocpT[j]->getTH1F());
307  BinLogX(h_assoc2pT[j]->getTH1F());
308  BinLogX(h_simulpT[j]->getTH1F());
309 
310  if (MABH)
311  {
312  BinLogX(h_efficPt_Quality05[j]->getTH1F());
313  BinLogX(h_efficPt_Quality075[j]->getTH1F());
314  //
315  BinLogX(h_assocpT_Quality05[j]->getTH1F());
316  BinLogX(h_assocpT_Quality075[j]->getTH1F());
317  }
318 
319  j++;
320  }
321  }
322  }
323  if (UseAssociators) {
324  edm::ESHandle<TrackAssociatorBase> theAssociator;
325  for (unsigned int w=0;w<associators.size();w++) {
326  setup.get<TrackAssociatorRecord>().get(associators[w],theAssociator);
327  associator.push_back( theAssociator.product() );
328  }
329  }
330 }
MuonTrackValidator::h_dzmeanh
std::vector< MonitorElement * > h_dzmeanh
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::h_assoc2phi
std::vector< MonitorElement * > h_assoc2phi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidatorBase::h_recoeta
std::vector< MonitorElement * > h_recoeta
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidatorBase::nrecHit_vs_nsimHit_sim2rec
std::vector< MonitorElement * > nrecHit_vs_nsimHit_sim2rec
Definition: MuonTrackValidatorBase.h:398
MuonTrackValidatorBase::h_DThits_eta
std::vector< MonitorElement * > h_DThits_eta
Definition: MuonTrackValidatorBase.h:408
MuonTrackValidator::h_dxypulleta
std::vector< MonitorElement * > h_dxypulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidator::h_ptpulleta
std::vector< MonitorElement * > h_ptpulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidatorBase::h_assoc2hit
std::vector< MonitorElement * > h_assoc2hit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::h_PurityVsQuality
std::vector< MonitorElement * > h_PurityVsQuality
Definition: MuonTrackValidatorBase.h:428
MuonTrackValidatorBase::h_recopT
std::vector< MonitorElement * > h_recopT
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidator::etares_vs_eta
std::vector< MonitorElement * > etares_vs_eta
Definition: MuonTrackValidator.h:157
MuonTrackValidator::h_dxypulletamean
std::vector< MonitorElement * > h_dxypulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_assoc2dxy
std::vector< MonitorElement * > h_assoc2dxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidatorBase::nrecHit_vs_nsimHit_rec2sim
std::vector< MonitorElement * > nrecHit_vs_nsimHit_rec2sim
Definition: MuonTrackValidatorBase.h:399
MuonTrackValidator::ptres_vs_phi
std::vector< MonitorElement * > ptres_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::h_fake_vs_phi
std::vector< MonitorElement * > h_fake_vs_phi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidator::h_losthits_eta
std::vector< MonitorElement * > h_losthits_eta
Definition: MuonTrackValidator.h:169
MuonTrackValidatorBase::h_assocpT_Quality075
std::vector< MonitorElement * > h_assocpT_Quality075
Definition: MuonTrackValidatorBase.h:393
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:519
MuonTrackValidatorBase::nrec_vs_nsim
std::vector< MonitorElement * > nrec_vs_nsim
Definition: MuonTrackValidatorBase.h:397
MuonTrackValidatorBase::h_pullDz
std::vector< MonitorElement * > h_pullDz
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidatorBase::h_recohit
std::vector< MonitorElement * > h_recohit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::h_etaSIM
std::vector< MonitorElement * > h_etaSIM
Definition: MuonTrackValidatorBase.h:379
MuonTrackValidatorBase::h_assocvertpos
std::vector< MonitorElement * > h_assocvertpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::h_thetapullphi
std::vector< MonitorElement * > h_thetapullphi
Definition: MuonTrackValidator.h:190
MuonTrackValidator::h_ptpullphi
std::vector< MonitorElement * > h_ptpullphi
Definition: MuonTrackValidator.h:190
MuonTrackValidatorBase::h_fakes
std::vector< MonitorElement * > h_fakes
Definition: MuonTrackValidatorBase.h:382
MuonTrackValidator::ptres_vs_eta
std::vector< MonitorElement * > ptres_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidator::cotThetares_vs_eta
std::vector< MonitorElement * > cotThetares_vs_eta
Definition: MuonTrackValidator.h:174
DQMStore::cd
void cd(void)
go to top directory (ie. root)
Definition: DQMStore.cc:214
MuonTrackValidator::h_dzpulleta
std::vector< MonitorElement * > h_dzpulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidatorBase::h_eta
std::vector< MonitorElement * > h_eta
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidator::h_thetapulleta
std::vector< MonitorElement * > h_thetapulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidator::phimean_vs_eta_phi
std::vector< MonitorElement * > phimean_vs_eta_phi
Definition: MuonTrackValidator.h:162
MuonTrackValidatorBase::h_pt
std::vector< MonitorElement * > h_pt
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidator::phires_vs_pt
std::vector< MonitorElement * > phires_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::h_pullQoverp
std::vector< MonitorElement * > h_pullQoverp
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidatorBase::h_effic_vs_dz
std::vector< MonitorElement * > h_effic_vs_dz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidator::h_dxymeanhPt
std::vector< MonitorElement * > h_dxymeanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidator::h_assochi2_prob
std::vector< MonitorElement * > h_assochi2_prob
Definition: MuonTrackValidator.h:165
linker.replace
def replace
Definition: linker.py:10
MuonTrackValidatorBase::h_CSChits_eta
std::vector< MonitorElement * > h_CSChits_eta
Definition: MuonTrackValidatorBase.h:408
MuonTrackValidator::h_dxyrmsh
std::vector< MonitorElement * > h_dxyrmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidator::cotThetares_vs_pt
std::vector< MonitorElement * > cotThetares_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidator::h_ptshifteta
std::vector< MonitorElement * > h_ptshifteta
Definition: MuonTrackValidator.h:158
MuonTrackValidator::dxyres_vs_eta
std::vector< MonitorElement * > dxyres_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidator::dzres_vs_pt
std::vector< MonitorElement * > dzres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::h_assoceta
std::vector< MonitorElement * > h_assoceta
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidatorBase::h_tracks
std::vector< MonitorElement * > h_tracks
Definition: MuonTrackValidatorBase.h:382
MuonTrackValidator::dzres_vs_eta
std::vector< MonitorElement * > dzres_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidatorBase::h_assoceta_Quality075
std::vector< MonitorElement * > h_assoceta_Quality075
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidator::h_cotThetameanhPt
std::vector< MonitorElement * > h_cotThetameanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidator::h_cotThetarmsh
std::vector< MonitorElement * > h_cotThetarmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::nDThits_vs_eta
std::vector< MonitorElement * > nDThits_vs_eta
Definition: MuonTrackValidatorBase.h:405
ExpressReco_HICollisions_FallBack.algo
tuple algo
Definition: ExpressReco_HICollisions_FallBack.py:2684
MuonTrackValidatorBase::h_effic_vs_zpos
std::vector< MonitorElement * > h_effic_vs_zpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::thetapull_vs_phi
std::vector< MonitorElement * > thetapull_vs_phi
Definition: MuonTrackValidator.h:188
MuonTrackValidatorBase::h_efficPt_Quality05
std::vector< MonitorElement * > h_efficPt_Quality05
Definition: MuonTrackValidatorBase.h:393
MuonTrackValidatorBase::h_efficPt
std::vector< MonitorElement * > h_efficPt
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidator::h_assochi2
std::vector< MonitorElement * > h_assochi2
Definition: MuonTrackValidator.h:165
MuonTrackValidatorBase::h_fakerate
std::vector< MonitorElement * > h_fakerate
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidator::chi2_vs_phi
std::vector< MonitorElement * > chi2_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::h_simuldxy
std::vector< MonitorElement * > h_simuldxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidator::h_phipulleta
std::vector< MonitorElement * > h_phipulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidatorBase::h_simuldz
std::vector< MonitorElement * > h_simuldz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidator::chi2_vs_nhits
std::vector< MonitorElement * > chi2_vs_nhits
Definition: MuonTrackValidator.h:157
MuonTrackValidator::h_cotThetameanh
std::vector< MonitorElement * > h_cotThetameanh
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::nhits_vs_eta
std::vector< MonitorElement * > nhits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidator::phires_vs_phi
std::vector< MonitorElement * > phires_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidator::h_chi2mean_vs_phi
std::vector< MonitorElement * > h_chi2mean_vs_phi
Definition: MuonTrackValidator.h:171
MuonTrackValidator::h_phipullphi
std::vector< MonitorElement * > h_phipullphi
Definition: MuonTrackValidator.h:190
MuonTrackValidatorBase::h_assoczpos
std::vector< MonitorElement * > h_assoczpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::h_dzpulletamean
std::vector< MonitorElement * > h_dzpulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidator::nhits_vs_phi
std::vector< MonitorElement * > nhits_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::h_simulphi
std::vector< MonitorElement * > h_simulphi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidatorBase::h_vertposSIM
std::vector< MonitorElement * > h_vertposSIM
Definition: MuonTrackValidatorBase.h:379
MuonTrackValidator::dxyres_vs_pt
std::vector< MonitorElement * > dxyres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::label
std::vector< edm::InputTag > label
Definition: MuonTrackValidatorBase.h:338
MuonTrackValidatorBase::h_assocSharedHit
std::vector< MonitorElement * > h_assocSharedHit
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidator::h_ptrmsh
std::vector< MonitorElement * > h_ptrmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidator::h_nchi2
std::vector< MonitorElement * > h_nchi2
Definition: MuonTrackValidator.h:154
MuonTrackValidatorBase::h_efficPt_Quality075
std::vector< MonitorElement * > h_efficPt_Quality075
Definition: MuonTrackValidatorBase.h:393
TrackerOfflineValidation_Dqm_cff.dirName
tuple dirName
Definition: TrackerOfflineValidation_Dqm_cff.py:56
MuonTrackValidator::h_phirmshPhi
std::vector< MonitorElement * > h_phirmshPhi
Definition: MuonTrackValidator.h:179
MuonTrackValidatorBase::h_assoc2dz
std::vector< MonitorElement * > h_assoc2dz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidatorBase::h_pullPhi
std::vector< MonitorElement * > h_pullPhi
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidator::h_nchi2_prob
std::vector< MonitorElement * > h_nchi2_prob
Definition: MuonTrackValidator.h:154
MuonTrackValidator::ptpull_vs_eta
std::vector< MonitorElement * > ptpull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidator::h_dxyrmshPt
std::vector< MonitorElement * > h_dxyrmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidatorBase::h_effic_vs_phi_Quality05
std::vector< MonitorElement * > h_effic_vs_phi_Quality05
Definition: MuonTrackValidatorBase.h:394
MuonTrackValidatorBase::h_recophi
std::vector< MonitorElement * > h_recophi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidator::dxypull_vs_eta
std::vector< MonitorElement * > dxypull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidator::h_ptrmshPhi
std::vector< MonitorElement * > h_ptrmshPhi
Definition: MuonTrackValidator.h:179
j
int j
Definition: DBlmapReader.cc:9
MuonTrackValidatorBase::h_effic_Quality05
std::vector< MonitorElement * > h_effic_Quality05
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidatorBase::h_effic
std::vector< MonitorElement * > h_effic
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidatorBase::h_effic_vs_phi_Quality075
std::vector< MonitorElement * > h_effic_vs_phi_Quality075
Definition: MuonTrackValidatorBase.h:394
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:390
MuonTrackValidatorBase::h_simulhit
std::vector< MonitorElement * > h_simulhit
Definition: MuonTrackValidatorBase.h:385
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:833
MuonTrackValidator::h_cotThetarmshPt
std::vector< MonitorElement * > h_cotThetarmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidator::phipull_vs_eta
std::vector< MonitorElement * > phipull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidatorBase::h_effic_vs_dxy
std::vector< MonitorElement * > h_effic_vs_dxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidator::thetapull_vs_eta
std::vector< MonitorElement * > thetapull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidatorBase::h_charge
std::vector< MonitorElement * > h_charge
Definition: MuonTrackValidatorBase.h:382
MuonTrackValidatorBase::h_effic_vs_phi
std::vector< MonitorElement * > h_effic_vs_phi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidator::h_phirmshPt
std::vector< MonitorElement * > h_phirmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidatorBase::h_assocphi
std::vector< MonitorElement * > h_assocphi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidator::ptmean_vs_eta_phi
std::vector< MonitorElement * > ptmean_vs_eta_phi
Definition: MuonTrackValidator.h:162
MuonTrackValidator::h_dzrmshPt
std::vector< MonitorElement * > h_dzrmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidator::h_thetapullphimean
std::vector< MonitorElement * > h_thetapullphimean
Definition: MuonTrackValidator.h:192
MuonTrackValidatorBase::h_assocpT
std::vector< MonitorElement * > h_assocpT
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::h_pullDxy
std::vector< MonitorElement * > h_pullDxy
Definition: MuonTrackValidatorBase.h:390
MuonTrackValidatorBase::h_simulpT
std::vector< MonitorElement * > h_simulpT
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::h_tracksSIM
std::vector< MonitorElement * > h_tracksSIM
Definition: MuonTrackValidatorBase.h:379
MuonTrackValidatorBase::h_assocphi_Quality05
std::vector< MonitorElement * > h_assocphi_Quality05
Definition: MuonTrackValidatorBase.h:394
MuonTrackValidator::h_dzmeanhPt
std::vector< MonitorElement * > h_dzmeanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidatorBase::h_simulzpos
std::vector< MonitorElement * > h_simulzpos
Definition: MuonTrackValidatorBase.h:389
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
MuonTrackValidatorBase::h_recodz
std::vector< MonitorElement * > h_recodz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidatorBase::h_effic_Quality075
std::vector< MonitorElement * > h_effic_Quality075
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidator::chi2_vs_eta
std::vector< MonitorElement * > chi2_vs_eta
Definition: MuonTrackValidator.h:168
MuonTrackValidatorBase::h_assocphi_Quality075
std::vector< MonitorElement * > h_assocphi_Quality075
Definition: MuonTrackValidatorBase.h:394
MuonTrackValidator::h_phirmsh
std::vector< MonitorElement * > h_phirmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::h_hits_eta
std::vector< MonitorElement * > h_hits_eta
Definition: MuonTrackValidatorBase.h:408
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:25
MuonTrackValidatorBase::h_hits
std::vector< MonitorElement * > h_hits
Definition: MuonTrackValidatorBase.h:382
MuonTrackValidator::h_losthits
std::vector< MonitorElement * > h_losthits
Definition: MuonTrackValidator.h:154
MuonTrackValidator::ptpull_vs_phi
std::vector< MonitorElement * > ptpull_vs_phi
Definition: MuonTrackValidator.h:188
edm::InputTag::process
std::string const & process() const
Definition: InputTag.h:29
MuonTrackValidatorBase::h_fake_vs_dz
std::vector< MonitorElement * > h_fake_vs_dz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidator::h_ptrmshPt
std::vector< MonitorElement * > h_ptrmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidatorBase::h_effic_vs_vertpos
std::vector< MonitorElement * > h_effic_vs_vertpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::h_thetapulletamean
std::vector< MonitorElement * > h_thetapulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_simulvertpos
std::vector< MonitorElement * > h_simulvertpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::ptres_vs_pt
std::vector< MonitorElement * > ptres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::associator
std::vector< const TrackAssociatorBase * > associator
Definition: MuonTrackValidatorBase.h:376
MuonTrackValidator::h_chi2meanhitsh
std::vector< MonitorElement * > h_chi2meanhitsh
Definition: MuonTrackValidator.h:171
MuonTrackValidatorBase::h_recodxy
std::vector< MonitorElement * > h_recodxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidator::phires_vs_eta
std::vector< MonitorElement * > phires_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidatorBase::h_ptSIM
std::vector< MonitorElement * > h_ptSIM
Definition: MuonTrackValidatorBase.h:379
MuonTrackValidator::h_ptpulletamean
std::vector< MonitorElement * > h_ptpulletamean
Definition: MuonTrackValidator.h:191
getTH2F
TH2F * getTH2F(std::string name, std::string process, std::string rootfolder, DQMStore *dbe_, bool verb, bool clone)
Definition: HcalHistoUtils.h:85
MuonTrackValidator::h_dxymeanh
std::vector< MonitorElement * > h_dxymeanh
Definition: MuonTrackValidator.h:182
DQMStore::goUp
void goUp(void)
equivalent to &quot;cd ..&quot;
Definition: DQMStore.cc:248
MuonTrackValidatorBase::h_assocdz
std::vector< MonitorElement * > h_assocdz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidator::h_hits_phi
std::vector< MonitorElement * > h_hits_phi
Definition: MuonTrackValidator.h:170
MuonTrackValidatorBase::h_assocpT_Quality05
std::vector< MonitorElement * > h_assocpT_Quality05
Definition: MuonTrackValidatorBase.h:393
MuonTrackValidatorBase::h_RPChits_eta
std::vector< MonitorElement * > h_RPChits_eta
Definition: MuonTrackValidatorBase.h:408
MuonTrackValidator::h_ptmeanhPt
std::vector< MonitorElement * > h_ptmeanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidatorBase::h_assoc2pT
std::vector< MonitorElement * > h_assoc2pT
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::h_fake_vs_dxy
std::vector< MonitorElement * > h_fake_vs_dxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidatorBase::h_assocFraction
std::vector< MonitorElement * > h_assocFraction
Definition: MuonTrackValidatorBase.h:402
MuonTrackValidator::h_phimeanhPt
std::vector< MonitorElement * > h_phimeanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidator::h_ptpullphimean
std::vector< MonitorElement * > h_ptpullphimean
Definition: MuonTrackValidator.h:192
MuonTrackValidator::h_phimeanh
std::vector< MonitorElement * > h_phimeanh
Definition: MuonTrackValidator.h:182
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:647
MuonTrackValidatorBase::associators
std::vector< std::string > associators
Definition: MuonTrackValidatorBase.h:344
MuonTrackValidator::nlosthits_vs_eta
std::vector< MonitorElement * > nlosthits_vs_eta
Definition: MuonTrackValidator.h:168
MuonTrackValidatorBase::h_assoceta_Quality05
std::vector< MonitorElement * > h_assoceta_Quality05
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidator::dzpull_vs_eta
std::vector< MonitorElement * > dzpull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidatorBase::h_assocdxy
std::vector< MonitorElement * > h_assocdxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidator::h_phimeanhPhi
std::vector< MonitorElement * > h_phimeanhPhi
Definition: MuonTrackValidator.h:184
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:26
MuonTrackValidatorBase::nRPChits_vs_eta
std::vector< MonitorElement * > nRPChits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidator::h_phipullphimean
std::vector< MonitorElement * > h_phipullphimean
Definition: MuonTrackValidator.h:192
MuonTrackValidatorBase::h_assochit
std::vector< MonitorElement * > h_assochit
Definition: MuonTrackValidatorBase.h:385
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:237
MuonTrackValidator::h_chi2meanh
std::vector< MonitorElement * > h_chi2meanh
Definition: MuonTrackValidator.h:169
MuonTrackValidatorBase::h_fakeratePt
std::vector< MonitorElement * > h_fakeratePt
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::nCSChits_vs_eta
std::vector< MonitorElement * > nCSChits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidator::h_phipulletamean
std::vector< MonitorElement * > h_phipulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_fake_vs_hit
std::vector< MonitorElement * > h_fake_vs_hit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::h_effic_vs_hit
std::vector< MonitorElement * > h_effic_vs_hit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::h_simuleta
std::vector< MonitorElement * > h_simuleta
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidator::h_ptmeanhPhi
std::vector< MonitorElement * > h_ptmeanhPhi
Definition: MuonTrackValidator.h:184
MuonTrackValidator::h_dzrmsh
std::vector< MonitorElement * > h_dzrmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::h_assoc2eta
std::vector< MonitorElement * > h_assoc2eta
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidator::phipull_vs_phi
std::vector< MonitorElement * > phipull_vs_phi
Definition: MuonTrackValidator.h:188
DQMStore::bookProfile2D
MonitorElement * bookProfile2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, int nchZ, double lowZ, double highZ, const char *option="s")
Definition: DQMStore.cc:977
void MuonTrackValidator::endRun ( edm::Run const &  ,
edm::EventSetup const &   
)
virtual

Method called at the end of the event loop.

Reimplemented from edm::EDAnalyzer.

Definition at line 955 of file MuonTrackValidator.cc.

References dbe_, FitSlicesYTool::getFittedMeanWithError(), FitSlicesYTool::getFittedSigmaWithError(), label, dbtoconf::out, and DQMStore::save().

955  {
956 
957  int w=0;
958  for (unsigned int ww=0;ww<associators.size();ww++){
959  for (unsigned int www=0;www<label.size();www++){
960 
961  //resolution of track params: get sigma from 2D histos
962  if(!skipHistoFit){
963  FitSlicesYTool fsyt_dxy(dxyres_vs_eta[w]);
964  fsyt_dxy.getFittedSigmaWithError(h_dxyrmsh[w]);
965  fsyt_dxy.getFittedMeanWithError(h_dxymeanh[w]);
966  FitSlicesYTool fsyt_dxyPt(dxyres_vs_pt[w]);
967  fsyt_dxyPt.getFittedSigmaWithError(h_dxyrmshPt[w]);
968  fsyt_dxyPt.getFittedMeanWithError(h_dxymeanhPt[w]);
969  FitSlicesYTool fsyt_pt(ptres_vs_eta[w]);
970  fsyt_pt.getFittedSigmaWithError(h_ptrmsh[w]);
971  fsyt_pt.getFittedMeanWithError(h_ptshifteta[w]);
972  FitSlicesYTool fsyt_ptPt(ptres_vs_pt[w]);
973  fsyt_ptPt.getFittedSigmaWithError(h_ptrmshPt[w]);
974  fsyt_ptPt.getFittedMeanWithError(h_ptmeanhPt[w]);
975  FitSlicesYTool fsyt_ptPhi(ptres_vs_phi[w]);
976  fsyt_ptPhi.getFittedSigmaWithError(h_ptrmshPhi[w]);
977  fsyt_ptPhi.getFittedMeanWithError(h_ptmeanhPhi[w]);
978  FitSlicesYTool fsyt_dz(dzres_vs_eta[w]);
979  fsyt_dz.getFittedSigmaWithError(h_dzrmsh[w]);
980  fsyt_dz.getFittedMeanWithError(h_dzmeanh[w]);
981  FitSlicesYTool fsyt_dzPt(dzres_vs_pt[w]);
982  fsyt_dzPt.getFittedSigmaWithError(h_dzrmshPt[w]);
983  fsyt_dzPt.getFittedMeanWithError(h_dzmeanhPt[w]);
984  FitSlicesYTool fsyt_phi(phires_vs_eta[w]);
985  fsyt_phi.getFittedSigmaWithError(h_phirmsh[w]);
986  fsyt_phi.getFittedMeanWithError(h_phimeanh[w]);
987  FitSlicesYTool fsyt_phiPt(phires_vs_pt[w]);
988  fsyt_phiPt.getFittedSigmaWithError(h_phirmshPt[w]);
989  fsyt_phiPt.getFittedMeanWithError(h_phimeanhPt[w]);
990  FitSlicesYTool fsyt_phiPhi(phires_vs_phi[w]);
991  fsyt_phiPhi.getFittedSigmaWithError(h_phirmshPhi[w]);
992  fsyt_phiPhi.getFittedMeanWithError(h_phimeanhPhi[w]);
993  FitSlicesYTool fsyt_cotTheta(cotThetares_vs_eta[w]);
994  fsyt_cotTheta.getFittedSigmaWithError(h_cotThetarmsh[w]);
995  fsyt_cotTheta.getFittedMeanWithError(h_cotThetameanh[w]);
996  FitSlicesYTool fsyt_cotThetaPt(cotThetares_vs_pt[w]);
997  fsyt_cotThetaPt.getFittedSigmaWithError(h_cotThetarmshPt[w]);
998  fsyt_cotThetaPt.getFittedMeanWithError(h_cotThetameanhPt[w]);
999  }
1000  //chi2 and #hit vs eta: get mean from 2D histos
1001  doProfileX(chi2_vs_eta[w],h_chi2meanh[w]);
1002  doProfileX(nhits_vs_eta[w],h_hits_eta[w]);
1003  doProfileX(nDThits_vs_eta[w],h_DThits_eta[w]);
1004  doProfileX(nCSChits_vs_eta[w],h_CSChits_eta[w]);
1005  doProfileX(nRPChits_vs_eta[w],h_RPChits_eta[w]);
1006 
1007  doProfileX(nlosthits_vs_eta[w],h_losthits_eta[w]);
1008  //vs phi
1009  doProfileX(chi2_vs_nhits[w],h_chi2meanhitsh[w]);
1010  // doProfileX(ptres_vs_eta[w],h_ptresmean_vs_eta[w]);
1011  // doProfileX(phires_vs_eta[w],h_phiresmean_vs_eta[w]);
1012  doProfileX(chi2_vs_phi[w],h_chi2mean_vs_phi[w]);
1013  doProfileX(nhits_vs_phi[w],h_hits_phi[w]);
1014  // doProfileX(ptres_vs_phi[w],h_ptresmean_vs_phi[w]);
1015  // doProfileX(phires_vs_phi[w],h_phiresmean_vs_phi[w]);
1016  if(!skipHistoFit){
1017  //pulls of track params vs eta: get sigma from 2D histos
1018  FitSlicesYTool fsyt_dxyp(dxypull_vs_eta[w]);
1019  fsyt_dxyp.getFittedSigmaWithError(h_dxypulleta[w]);
1020  fsyt_dxyp.getFittedMeanWithError(h_dxypulletamean[w]);
1021  FitSlicesYTool fsyt_ptp(ptpull_vs_eta[w]);
1022  fsyt_ptp.getFittedSigmaWithError(h_ptpulleta[w]);
1023  fsyt_ptp.getFittedMeanWithError(h_ptpulletamean[w]);
1024  FitSlicesYTool fsyt_dzp(dzpull_vs_eta[w]);
1025  fsyt_dzp.getFittedSigmaWithError(h_dzpulleta[w]);
1026  fsyt_dzp.getFittedMeanWithError(h_dzpulletamean[w]);
1027  FitSlicesYTool fsyt_phip(phipull_vs_eta[w]);
1028  fsyt_phip.getFittedSigmaWithError(h_phipulleta[w]);
1029  fsyt_phip.getFittedMeanWithError(h_phipulletamean[w]);
1030  FitSlicesYTool fsyt_thetap(thetapull_vs_eta[w]);
1031  fsyt_thetap.getFittedSigmaWithError(h_thetapulleta[w]);
1032  fsyt_thetap.getFittedMeanWithError(h_thetapulletamean[w]);
1033  //vs phi
1034  FitSlicesYTool fsyt_ptpPhi(ptpull_vs_phi[w]);
1035  fsyt_ptpPhi.getFittedSigmaWithError(h_ptpullphi[w]);
1036  fsyt_ptpPhi.getFittedMeanWithError(h_ptpullphimean[w]);
1037  FitSlicesYTool fsyt_phipPhi(phipull_vs_phi[w]);
1038  fsyt_phipPhi.getFittedSigmaWithError(h_phipullphi[w]);
1039  fsyt_phipPhi.getFittedMeanWithError(h_phipullphimean[w]);
1040  FitSlicesYTool fsyt_thetapPhi(thetapull_vs_phi[w]);
1041  fsyt_thetapPhi.getFittedSigmaWithError(h_thetapullphi[w]);
1042  fsyt_thetapPhi.getFittedMeanWithError(h_thetapullphimean[w]);
1043  //effic&fake
1044  fillPlotFromVectors(h_effic[w],totASSeta[w],totSIMeta[w],"effic");
1045  fillPlotFromVectors(h_fakerate[w],totASS2eta[w],totRECeta[w],"fakerate");
1046  fillPlotFromVectors(h_efficPt[w],totASSpT[w],totSIMpT[w],"effic");
1047  fillPlotFromVectors(h_fakeratePt[w],totASS2pT[w],totRECpT[w],"fakerate");
1048  fillPlotFromVectors(h_effic_vs_hit[w],totASS_hit[w],totSIM_hit[w],"effic");
1049  fillPlotFromVectors(h_fake_vs_hit[w],totASS2_hit[w],totREC_hit[w],"fakerate");
1050  fillPlotFromVectors(h_effic_vs_phi[w],totASS_phi[w],totSIM_phi[w],"effic");
1051  fillPlotFromVectors(h_fake_vs_phi[w],totASS2_phi[w],totREC_phi[w],"fakerate");
1052  fillPlotFromVectors(h_effic_vs_dxy[w],totASS_dxy[w],totSIM_dxy[w],"effic");
1053  fillPlotFromVectors(h_fake_vs_dxy[w],totASS2_dxy[w],totREC_dxy[w],"fakerate");
1054  fillPlotFromVectors(h_effic_vs_dz[w],totASS_dz[w],totSIM_dz[w],"effic");
1055  fillPlotFromVectors(h_fake_vs_dz[w],totASS2_dz[w],totREC_dz[w],"fakerate");
1056  fillPlotFromVectors(h_effic_vs_vertpos[w],totASS_vertpos[w],totSIM_vertpos[w],"effic");
1057  fillPlotFromVectors(h_effic_vs_zpos[w],totASS_zpos[w],totSIM_zpos[w],"effic");
1058 
1059  if (MABH)
1060  {
1061  fillPlotFromVectors(h_effic_Quality05[w] ,totASSeta_Quality05[w] ,totSIMeta[w],"effic");
1062  fillPlotFromVectors(h_effic_Quality075[w],totASSeta_Quality075[w],totSIMeta[w],"effic");
1063  fillPlotFromVectors(h_efficPt_Quality05[w] ,totASSpT_Quality05[w] ,totSIMpT[w],"effic");
1064  fillPlotFromVectors(h_efficPt_Quality075[w],totASSpT_Quality075[w],totSIMpT[w],"effic");
1065  fillPlotFromVectors(h_effic_vs_phi_Quality05[w] ,totASS_phi_Quality05[w],totSIM_phi[w],"effic");
1066  fillPlotFromVectors(h_effic_vs_phi_Quality075[w],totASS_phi_Quality075[w],totSIM_phi[w],"effic");
1067  }
1068  }
1069 
1070  fillPlotFromVector(h_recoeta[w],totRECeta[w]);
1071  fillPlotFromVector(h_simuleta[w],totSIMeta[w]);
1072  fillPlotFromVector(h_assoceta[w],totASSeta[w]);
1073  fillPlotFromVector(h_assoc2eta[w],totASS2eta[w]);
1074 
1075  fillPlotFromVector(h_recopT[w],totRECpT[w]);
1076  fillPlotFromVector(h_simulpT[w],totSIMpT[w]);
1077  fillPlotFromVector(h_assocpT[w],totASSpT[w]);
1078  fillPlotFromVector(h_assoc2pT[w],totASS2pT[w]);
1079 
1080  fillPlotFromVector(h_recohit[w],totREC_hit[w]);
1081  fillPlotFromVector(h_simulhit[w],totSIM_hit[w]);
1082  fillPlotFromVector(h_assochit[w],totASS_hit[w]);
1083  fillPlotFromVector(h_assoc2hit[w],totASS2_hit[w]);
1084 
1085  fillPlotFromVector(h_recophi[w],totREC_phi[w]);
1086  fillPlotFromVector(h_simulphi[w],totSIM_phi[w]);
1087  fillPlotFromVector(h_assocphi[w],totASS_phi[w]);
1088  fillPlotFromVector(h_assoc2phi[w],totASS2_phi[w]);
1089 
1090  fillPlotFromVector(h_recodxy[w],totREC_dxy[w]);
1091  fillPlotFromVector(h_simuldxy[w],totSIM_dxy[w]);
1092  fillPlotFromVector(h_assocdxy[w],totASS_dxy[w]);
1093  fillPlotFromVector(h_assoc2dxy[w],totASS2_dxy[w]);
1094 
1095  fillPlotFromVector(h_recodz[w],totREC_dz[w]);
1096  fillPlotFromVector(h_simuldz[w],totSIM_dz[w]);
1097  fillPlotFromVector(h_assocdz[w],totASS_dz[w]);
1098  fillPlotFromVector(h_assoc2dz[w],totASS2_dz[w]);
1099 
1100  fillPlotFromVector(h_simulvertpos[w],totSIM_vertpos[w]);
1101  fillPlotFromVector(h_assocvertpos[w],totASS_vertpos[w]);
1102 
1103  fillPlotFromVector(h_simulzpos[w],totSIM_zpos[w]);
1104  fillPlotFromVector(h_assoczpos[w],totASS_zpos[w]);
1105 
1106  if (MABH) {
1107  fillPlotFromVector(h_assoceta_Quality05[w] ,totASSeta_Quality05[w]);
1108  fillPlotFromVector(h_assoceta_Quality075[w],totASSeta_Quality075[w]);
1109  fillPlotFromVector(h_assocpT_Quality05[w] ,totASSpT_Quality05[w]);
1110  fillPlotFromVector(h_assocpT_Quality075[w],totASSpT_Quality075[w]);
1111  fillPlotFromVector(h_assocphi_Quality05[w] ,totASS_phi_Quality05[w]);
1112  fillPlotFromVector(h_assocphi_Quality075[w],totASS_phi_Quality075[w]);
1113  }
1114 
1115  w++;
1116  }
1117  }
1118 
1119  if ( out.size() != 0 && dbe_ ) dbe_->save(out);
1120 }
MuonTrackValidator::h_dzmeanh
std::vector< MonitorElement * > h_dzmeanh
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::h_assoc2phi
std::vector< MonitorElement * > h_assoc2phi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidatorBase::h_recoeta
std::vector< MonitorElement * > h_recoeta
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidatorBase::h_DThits_eta
std::vector< MonitorElement * > h_DThits_eta
Definition: MuonTrackValidatorBase.h:408
MuonTrackValidator::h_dxypulleta
std::vector< MonitorElement * > h_dxypulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidator::h_ptpulleta
std::vector< MonitorElement * > h_ptpulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidatorBase::h_assoc2hit
std::vector< MonitorElement * > h_assoc2hit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::totASS_zpos
std::vector< std::vector< int > > totASS_zpos
Definition: MuonTrackValidatorBase.h:425
MuonTrackValidatorBase::totASSeta_Quality05
std::vector< std::vector< int > > totASSeta_Quality05
Definition: MuonTrackValidatorBase.h:429
MuonTrackValidatorBase::h_recopT
std::vector< MonitorElement * > h_recopT
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::totREC_phi
std::vector< std::vector< int > > totREC_phi
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidator::h_dxypulletamean
std::vector< MonitorElement * > h_dxypulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_assoc2dxy
std::vector< MonitorElement * > h_assoc2dxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidator::ptres_vs_phi
std::vector< MonitorElement * > ptres_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::totASS2_hit
std::vector< std::vector< int > > totASS2_hit
Definition: MuonTrackValidatorBase.h:421
MuonTrackValidatorBase::h_fake_vs_phi
std::vector< MonitorElement * > h_fake_vs_phi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidator::h_losthits_eta
std::vector< MonitorElement * > h_losthits_eta
Definition: MuonTrackValidator.h:169
MuonTrackValidatorBase::fillPlotFromVectors
void fillPlotFromVectors(MonitorElement *h, std::vector< int > &numerator, std::vector< int > &denominator, std::string type)
Definition: MuonTrackValidatorBase.h:136
MuonTrackValidatorBase::h_assocpT_Quality075
std::vector< MonitorElement * > h_assocpT_Quality075
Definition: MuonTrackValidatorBase.h:393
MuonTrackValidatorBase::totRECpT
std::vector< std::vector< int > > totRECpT
Definition: MuonTrackValidatorBase.h:420
MuonTrackValidatorBase::totREC_dz
std::vector< std::vector< int > > totREC_dz
Definition: MuonTrackValidatorBase.h:424
MuonTrackValidatorBase::h_recohit
std::vector< MonitorElement * > h_recohit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::h_assocvertpos
std::vector< MonitorElement * > h_assocvertpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::h_thetapullphi
std::vector< MonitorElement * > h_thetapullphi
Definition: MuonTrackValidator.h:190
MuonTrackValidator::h_ptpullphi
std::vector< MonitorElement * > h_ptpullphi
Definition: MuonTrackValidator.h:190
MuonTrackValidator::ptres_vs_eta
std::vector< MonitorElement * > ptres_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidator::cotThetares_vs_eta
std::vector< MonitorElement * > cotThetares_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidator::h_dzpulleta
std::vector< MonitorElement * > h_dzpulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidator::h_thetapulleta
std::vector< MonitorElement * > h_thetapulleta
Definition: MuonTrackValidator.h:189
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:1883
MuonTrackValidatorBase::totASS_dz
std::vector< std::vector< int > > totASS_dz
Definition: MuonTrackValidatorBase.h:424
MuonTrackValidator::phires_vs_pt
std::vector< MonitorElement * > phires_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::h_effic_vs_dz
std::vector< MonitorElement * > h_effic_vs_dz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidator::h_dxymeanhPt
std::vector< MonitorElement * > h_dxymeanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidatorBase::totSIM_zpos
std::vector< std::vector< int > > totSIM_zpos
Definition: MuonTrackValidatorBase.h:425
MuonTrackValidatorBase::h_CSChits_eta
std::vector< MonitorElement * > h_CSChits_eta
Definition: MuonTrackValidatorBase.h:408
MuonTrackValidator::h_dxyrmsh
std::vector< MonitorElement * > h_dxyrmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::totASS2_dz
std::vector< std::vector< int > > totASS2_dz
Definition: MuonTrackValidatorBase.h:424
MuonTrackValidator::cotThetares_vs_pt
std::vector< MonitorElement * > cotThetares_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidator::h_ptshifteta
std::vector< MonitorElement * > h_ptshifteta
Definition: MuonTrackValidator.h:158
MuonTrackValidator::dxyres_vs_eta
std::vector< MonitorElement * > dxyres_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidatorBase::totASS_phi_Quality05
std::vector< std::vector< int > > totASS_phi_Quality05
Definition: MuonTrackValidatorBase.h:431
MuonTrackValidator::dzres_vs_pt
std::vector< MonitorElement * > dzres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::h_assoceta
std::vector< MonitorElement * > h_assoceta
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidator::dzres_vs_eta
std::vector< MonitorElement * > dzres_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidatorBase::h_assoceta_Quality075
std::vector< MonitorElement * > h_assoceta_Quality075
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidator::h_cotThetameanhPt
std::vector< MonitorElement * > h_cotThetameanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidator::h_cotThetarmsh
std::vector< MonitorElement * > h_cotThetarmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::nDThits_vs_eta
std::vector< MonitorElement * > nDThits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidatorBase::totASS_phi
std::vector< std::vector< int > > totASS_phi
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::h_effic_vs_zpos
std::vector< MonitorElement * > h_effic_vs_zpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::thetapull_vs_phi
std::vector< MonitorElement * > thetapull_vs_phi
Definition: MuonTrackValidator.h:188
MuonTrackValidatorBase::h_efficPt_Quality05
std::vector< MonitorElement * > h_efficPt_Quality05
Definition: MuonTrackValidatorBase.h:393
MuonTrackValidatorBase::totASS_dxy
std::vector< std::vector< int > > totASS_dxy
Definition: MuonTrackValidatorBase.h:423
MuonTrackValidatorBase::h_efficPt
std::vector< MonitorElement * > h_efficPt
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::h_fakerate
std::vector< MonitorElement * > h_fakerate
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidator::chi2_vs_phi
std::vector< MonitorElement * > chi2_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::h_simuldxy
std::vector< MonitorElement * > h_simuldxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidatorBase::totSIM_dxy
std::vector< std::vector< int > > totSIM_dxy
Definition: MuonTrackValidatorBase.h:423
MuonTrackValidator::h_phipulleta
std::vector< MonitorElement * > h_phipulleta
Definition: MuonTrackValidator.h:189
MuonTrackValidatorBase::h_simuldz
std::vector< MonitorElement * > h_simuldz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidator::chi2_vs_nhits
std::vector< MonitorElement * > chi2_vs_nhits
Definition: MuonTrackValidator.h:157
MuonTrackValidatorBase::totASSeta_Quality075
std::vector< std::vector< int > > totASSeta_Quality075
Definition: MuonTrackValidatorBase.h:429
MuonTrackValidator::h_cotThetameanh
std::vector< MonitorElement * > h_cotThetameanh
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::nhits_vs_eta
std::vector< MonitorElement * > nhits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidator::phires_vs_phi
std::vector< MonitorElement * > phires_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidator::h_chi2mean_vs_phi
std::vector< MonitorElement * > h_chi2mean_vs_phi
Definition: MuonTrackValidator.h:171
MuonTrackValidator::h_phipullphi
std::vector< MonitorElement * > h_phipullphi
Definition: MuonTrackValidator.h:190
MuonTrackValidatorBase::h_assoczpos
std::vector< MonitorElement * > h_assoczpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::h_dzpulletamean
std::vector< MonitorElement * > h_dzpulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidator::nhits_vs_phi
std::vector< MonitorElement * > nhits_vs_phi
Definition: MuonTrackValidator.h:159
MuonTrackValidatorBase::h_simulphi
std::vector< MonitorElement * > h_simulphi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidatorBase::totASS2pT
std::vector< std::vector< int > > totASS2pT
Definition: MuonTrackValidatorBase.h:420
MuonTrackValidator::dxyres_vs_pt
std::vector< MonitorElement * > dxyres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::doProfileX
virtual void doProfileX(TH2 *th2, MonitorElement *me)
Definition: MuonTrackValidatorBase.h:106
MuonTrackValidatorBase::label
std::vector< edm::InputTag > label
Definition: MuonTrackValidatorBase.h:338
MuonTrackValidator::h_ptrmsh
std::vector< MonitorElement * > h_ptrmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::h_efficPt_Quality075
std::vector< MonitorElement * > h_efficPt_Quality075
Definition: MuonTrackValidatorBase.h:393
MuonTrackValidatorBase::totASS2eta
std::vector< std::vector< int > > totASS2eta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::totSIMpT
std::vector< std::vector< int > > totSIMpT
Definition: MuonTrackValidatorBase.h:420
MuonTrackValidatorBase::totREC_dxy
std::vector< std::vector< int > > totREC_dxy
Definition: MuonTrackValidatorBase.h:423
MuonTrackValidator::h_phirmshPhi
std::vector< MonitorElement * > h_phirmshPhi
Definition: MuonTrackValidator.h:179
MuonTrackValidatorBase::h_assoc2dz
std::vector< MonitorElement * > h_assoc2dz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidatorBase::totSIM_hit
std::vector< std::vector< int > > totSIM_hit
Definition: MuonTrackValidatorBase.h:421
MuonTrackValidator::ptpull_vs_eta
std::vector< MonitorElement * > ptpull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidator::h_dxyrmshPt
std::vector< MonitorElement * > h_dxyrmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidatorBase::h_effic_vs_phi_Quality05
std::vector< MonitorElement * > h_effic_vs_phi_Quality05
Definition: MuonTrackValidatorBase.h:394
MuonTrackValidatorBase::h_recophi
std::vector< MonitorElement * > h_recophi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidatorBase::totASSpT
std::vector< std::vector< int > > totASSpT
Definition: MuonTrackValidatorBase.h:420
MuonTrackValidator::dxypull_vs_eta
std::vector< MonitorElement * > dxypull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidator::h_ptrmshPhi
std::vector< MonitorElement * > h_ptrmshPhi
Definition: MuonTrackValidator.h:179
MuonTrackValidatorBase::totASS_phi_Quality075
std::vector< std::vector< int > > totASS_phi_Quality075
Definition: MuonTrackValidatorBase.h:431
MuonTrackValidatorBase::h_effic_Quality05
std::vector< MonitorElement * > h_effic_Quality05
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidatorBase::h_effic
std::vector< MonitorElement * > h_effic
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidatorBase::h_effic_vs_phi_Quality075
std::vector< MonitorElement * > h_effic_vs_phi_Quality075
Definition: MuonTrackValidatorBase.h:394
MuonTrackValidatorBase::totSIM_vertpos
std::vector< std::vector< int > > totSIM_vertpos
Definition: MuonTrackValidatorBase.h:425
MuonTrackValidatorBase::totASSpT_Quality05
std::vector< std::vector< int > > totASSpT_Quality05
Definition: MuonTrackValidatorBase.h:430
MuonTrackValidatorBase::totREC_hit
std::vector< std::vector< int > > totREC_hit
Definition: MuonTrackValidatorBase.h:421
MuonTrackValidatorBase::h_simulhit
std::vector< MonitorElement * > h_simulhit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidator::h_cotThetarmshPt
std::vector< MonitorElement * > h_cotThetarmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidator::phipull_vs_eta
std::vector< MonitorElement * > phipull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidatorBase::h_effic_vs_dxy
std::vector< MonitorElement * > h_effic_vs_dxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidatorBase::totSIMeta
std::vector< std::vector< int > > totSIMeta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::totASS2_dxy
std::vector< std::vector< int > > totASS2_dxy
Definition: MuonTrackValidatorBase.h:423
MuonTrackValidator::thetapull_vs_eta
std::vector< MonitorElement * > thetapull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidatorBase::h_effic_vs_phi
std::vector< MonitorElement * > h_effic_vs_phi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidator::h_phirmshPt
std::vector< MonitorElement * > h_phirmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidatorBase::h_assocphi
std::vector< MonitorElement * > h_assocphi
Definition: MuonTrackValidatorBase.h:386
MuonTrackValidator::h_dzrmshPt
std::vector< MonitorElement * > h_dzrmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidatorBase::totASS2_phi
std::vector< std::vector< int > > totASS2_phi
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidator::h_thetapullphimean
std::vector< MonitorElement * > h_thetapullphimean
Definition: MuonTrackValidator.h:192
MuonTrackValidatorBase::totASSpT_Quality075
std::vector< std::vector< int > > totASSpT_Quality075
Definition: MuonTrackValidatorBase.h:430
MuonTrackValidatorBase::h_assocpT
std::vector< MonitorElement * > h_assocpT
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::h_simulpT
std::vector< MonitorElement * > h_simulpT
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::h_assocphi_Quality05
std::vector< MonitorElement * > h_assocphi_Quality05
Definition: MuonTrackValidatorBase.h:394
MuonTrackValidator::h_dzmeanhPt
std::vector< MonitorElement * > h_dzmeanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidatorBase::h_simulzpos
std::vector< MonitorElement * > h_simulzpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidatorBase::h_recodz
std::vector< MonitorElement * > h_recodz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidatorBase::h_effic_Quality075
std::vector< MonitorElement * > h_effic_Quality075
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidator::chi2_vs_eta
std::vector< MonitorElement * > chi2_vs_eta
Definition: MuonTrackValidator.h:168
MuonTrackValidatorBase::h_assocphi_Quality075
std::vector< MonitorElement * > h_assocphi_Quality075
Definition: MuonTrackValidatorBase.h:394
MuonTrackValidator::h_phirmsh
std::vector< MonitorElement * > h_phirmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::totSIM_phi
std::vector< std::vector< int > > totSIM_phi
Definition: MuonTrackValidatorBase.h:422
MuonTrackValidatorBase::fillPlotFromVector
void fillPlotFromVector(MonitorElement *h, std::vector< int > &vec)
Definition: MuonTrackValidatorBase.h:130
MuonTrackValidatorBase::h_hits_eta
std::vector< MonitorElement * > h_hits_eta
Definition: MuonTrackValidatorBase.h:408
MuonTrackValidator::ptpull_vs_phi
std::vector< MonitorElement * > ptpull_vs_phi
Definition: MuonTrackValidator.h:188
MuonTrackValidatorBase::h_fake_vs_dz
std::vector< MonitorElement * > h_fake_vs_dz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidator::h_ptrmshPt
std::vector< MonitorElement * > h_ptrmshPt
Definition: MuonTrackValidator.h:178
MuonTrackValidatorBase::h_effic_vs_vertpos
std::vector< MonitorElement * > h_effic_vs_vertpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::h_thetapulletamean
std::vector< MonitorElement * > h_thetapulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::totASSeta
std::vector< std::vector< int > > totASSeta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidatorBase::h_simulvertpos
std::vector< MonitorElement * > h_simulvertpos
Definition: MuonTrackValidatorBase.h:389
MuonTrackValidator::ptres_vs_pt
std::vector< MonitorElement * > ptres_vs_pt
Definition: MuonTrackValidator.h:175
MuonTrackValidatorBase::totRECeta
std::vector< std::vector< int > > totRECeta
Definition: MuonTrackValidatorBase.h:419
MuonTrackValidator::h_chi2meanhitsh
std::vector< MonitorElement * > h_chi2meanhitsh
Definition: MuonTrackValidator.h:171
MuonTrackValidatorBase::h_recodxy
std::vector< MonitorElement * > h_recodxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidator::phires_vs_eta
std::vector< MonitorElement * > phires_vs_eta
Definition: MuonTrackValidator.h:174
MuonTrackValidator::h_ptpulletamean
std::vector< MonitorElement * > h_ptpulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidator::h_dxymeanh
std::vector< MonitorElement * > h_dxymeanh
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::h_assocdz
std::vector< MonitorElement * > h_assocdz
Definition: MuonTrackValidatorBase.h:388
MuonTrackValidator::h_hits_phi
std::vector< MonitorElement * > h_hits_phi
Definition: MuonTrackValidator.h:170
MuonTrackValidatorBase::h_assocpT_Quality05
std::vector< MonitorElement * > h_assocpT_Quality05
Definition: MuonTrackValidatorBase.h:393
MuonTrackValidatorBase::h_RPChits_eta
std::vector< MonitorElement * > h_RPChits_eta
Definition: MuonTrackValidatorBase.h:408
MuonTrackValidator::h_ptmeanhPt
std::vector< MonitorElement * > h_ptmeanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidatorBase::h_assoc2pT
std::vector< MonitorElement * > h_assoc2pT
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::h_fake_vs_dxy
std::vector< MonitorElement * > h_fake_vs_dxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidator::h_phimeanhPt
std::vector< MonitorElement * > h_phimeanhPt
Definition: MuonTrackValidator.h:183
MuonTrackValidator::h_ptpullphimean
std::vector< MonitorElement * > h_ptpullphimean
Definition: MuonTrackValidator.h:192
MuonTrackValidator::h_phimeanh
std::vector< MonitorElement * > h_phimeanh
Definition: MuonTrackValidator.h:182
MuonTrackValidatorBase::totSIM_dz
std::vector< std::vector< int > > totSIM_dz
Definition: MuonTrackValidatorBase.h:424
MuonTrackValidatorBase::associators
std::vector< std::string > associators
Definition: MuonTrackValidatorBase.h:344
MuonTrackValidator::nlosthits_vs_eta
std::vector< MonitorElement * > nlosthits_vs_eta
Definition: MuonTrackValidator.h:168
MuonTrackValidatorBase::h_assoceta_Quality05
std::vector< MonitorElement * > h_assoceta_Quality05
Definition: MuonTrackValidatorBase.h:392
MuonTrackValidator::dzpull_vs_eta
std::vector< MonitorElement * > dzpull_vs_eta
Definition: MuonTrackValidator.h:187
MuonTrackValidatorBase::h_assocdxy
std::vector< MonitorElement * > h_assocdxy
Definition: MuonTrackValidatorBase.h:387
MuonTrackValidator::h_phimeanhPhi
std::vector< MonitorElement * > h_phimeanhPhi
Definition: MuonTrackValidator.h:184
MuonTrackValidatorBase::nRPChits_vs_eta
std::vector< MonitorElement * > nRPChits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidatorBase::totASS_vertpos
std::vector< std::vector< int > > totASS_vertpos
Definition: MuonTrackValidatorBase.h:425
MuonTrackValidator::h_phipullphimean
std::vector< MonitorElement * > h_phipullphimean
Definition: MuonTrackValidator.h:192
MuonTrackValidatorBase::h_assochit
std::vector< MonitorElement * > h_assochit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::totASS_hit
std::vector< std::vector< int > > totASS_hit
Definition: MuonTrackValidatorBase.h:421
MuonTrackValidator::h_chi2meanh
std::vector< MonitorElement * > h_chi2meanh
Definition: MuonTrackValidator.h:169
MuonTrackValidatorBase::h_fakeratePt
std::vector< MonitorElement * > h_fakeratePt
Definition: MuonTrackValidatorBase.h:384
MuonTrackValidatorBase::nCSChits_vs_eta
std::vector< MonitorElement * > nCSChits_vs_eta
Definition: MuonTrackValidatorBase.h:405
MuonTrackValidator::h_phipulletamean
std::vector< MonitorElement * > h_phipulletamean
Definition: MuonTrackValidator.h:191
MuonTrackValidatorBase::h_fake_vs_hit
std::vector< MonitorElement * > h_fake_vs_hit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::h_effic_vs_hit
std::vector< MonitorElement * > h_effic_vs_hit
Definition: MuonTrackValidatorBase.h:385
MuonTrackValidatorBase::h_simuleta
std::vector< MonitorElement * > h_simuleta
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidator::h_ptmeanhPhi
std::vector< MonitorElement * > h_ptmeanhPhi
Definition: MuonTrackValidator.h:184
MuonTrackValidator::h_dzrmsh
std::vector< MonitorElement * > h_dzrmsh
Definition: MuonTrackValidator.h:177
MuonTrackValidatorBase::h_assoc2eta
std::vector< MonitorElement * > h_assoc2eta
Definition: MuonTrackValidatorBase.h:383
MuonTrackValidator::phipull_vs_phi
std::vector< MonitorElement * > phipull_vs_phi
Definition: MuonTrackValidator.h:188
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 1124 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().

1125  {
1126  pt = track.pt();
1127  ptError = track.ptError();
1128  qoverp = track.qoverp();
1129  qoverpError = track.qoverpError();
1130  lambda = track.lambda();
1131  lambdaError = track.lambdaError();
1132  phi = track.phi();
1133  phiError = track.phiError();
1134 
1135 }
reco::TrackBase::qoverp
double qoverp() const
q/p
Definition: TrackBase.h:114
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:138
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:130
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:187
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:200
reco::TrackBase::lambda
double lambda() const
Lambda angle.
Definition: TrackBase.h:118
ExpressReco_HICollisions_FallBack.pt
tuple pt
Definition: ExpressReco_HICollisions_FallBack.py:1070
reco::TrackBase::qoverpError
double qoverpError() const
error on signed transverse curvature
Definition: TrackBase.h:185
reco::TrackBase::lambdaError
double lambdaError() const
error on lambda
Definition: TrackBase.h:196
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 1138 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().

1139  {
1140 
1141  pt = gsfTrack.ptMode();
1142  ptError = gsfTrack.ptModeError();
1143  qoverp = gsfTrack.qoverpMode();
1144  qoverpError = gsfTrack.qoverpModeError();
1145  lambda = gsfTrack.lambdaMode();
1146  lambdaError = gsfTrack.lambdaModeError();
1147  phi = gsfTrack.phiMode();
1148  phiError = gsfTrack.phiModeError();
1149 
1150 }
reco::GsfTrack::phiModeError
double phiModeError() const
error on phi from mode
Definition: GsfTrack.h:93
reco::GsfTrack::lambdaMode
double lambdaMode() const
Lambda angle from mode.
Definition: GsfTrack.h:44
reco::GsfTrack::qoverpMode
double qoverpMode() const
q/p from mode
Definition: GsfTrack.h:40
reco::GsfTrack::ptModeError
double ptModeError() const
error on Pt (set to 1000 TeV if charge==0 for safety) from mode
Definition: GsfTrack.h:80
ExpressReco_HICollisions_FallBack.pt
tuple pt
Definition: ExpressReco_HICollisions_FallBack.py:1070
reco::GsfTrack::qoverpModeError
double qoverpModeError() const
error on signed transverse curvature from mode
Definition: GsfTrack.h:78
reco::GsfTrack::phiMode
double phiMode() const
azimuthal angle of momentum vector from mode
Definition: GsfTrack.h:56
reco::GsfTrack::lambdaModeError
double lambdaModeError() const
error on lambda from mode
Definition: GsfTrack.h:89
reco::GsfTrack::ptMode
double ptMode() const
track transverse momentum from mode
Definition: GsfTrack.h:48
phi
Definition: DDAxes.h:10

Member Data Documentation

edm::InputTag MuonTrackValidator::associatormap
private

Definition at line 138 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::BiDirectional_RecoToSim_association
private

Definition at line 149 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

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

Definition at line 168 of file MuonTrackValidator.h.

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

Definition at line 157 of file MuonTrackValidator.h.

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

Definition at line 159 of file MuonTrackValidator.h.

CosmicTrackingParticleSelector MuonTrackValidator::cosmictpSelector
private

Definition at line 146 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

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

Definition at line 174 of file MuonTrackValidator.h.

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

Definition at line 175 of file MuonTrackValidator.h.

std::string MuonTrackValidator::dirName_
private

Definition at line 137 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

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

Definition at line 187 of file MuonTrackValidator.h.

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

Definition at line 174 of file MuonTrackValidator.h.

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

Definition at line 175 of file MuonTrackValidator.h.

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

Definition at line 187 of file MuonTrackValidator.h.

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

Definition at line 174 of file MuonTrackValidator.h.

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

Definition at line 175 of file MuonTrackValidator.h.

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

Definition at line 157 of file MuonTrackValidator.h.

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

Definition at line 165 of file MuonTrackValidator.h.

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

Definition at line 165 of file MuonTrackValidator.h.

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

Definition at line 171 of file MuonTrackValidator.h.

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

Definition at line 169 of file MuonTrackValidator.h.

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

Definition at line 171 of file MuonTrackValidator.h.

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

Definition at line 182 of file MuonTrackValidator.h.

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

Definition at line 183 of file MuonTrackValidator.h.

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

Definition at line 177 of file MuonTrackValidator.h.

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

Definition at line 178 of file MuonTrackValidator.h.

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

Definition at line 182 of file MuonTrackValidator.h.

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

Definition at line 183 of file MuonTrackValidator.h.

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

Definition at line 189 of file MuonTrackValidator.h.

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

Definition at line 191 of file MuonTrackValidator.h.

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

Definition at line 177 of file MuonTrackValidator.h.

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

Definition at line 178 of file MuonTrackValidator.h.

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

Definition at line 182 of file MuonTrackValidator.h.

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

Definition at line 183 of file MuonTrackValidator.h.

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

Definition at line 189 of file MuonTrackValidator.h.

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

Definition at line 191 of file MuonTrackValidator.h.

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

Definition at line 177 of file MuonTrackValidator.h.

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

Definition at line 178 of file MuonTrackValidator.h.

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

Definition at line 170 of file MuonTrackValidator.h.

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

Definition at line 154 of file MuonTrackValidator.h.

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

Definition at line 169 of file MuonTrackValidator.h.

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

Definition at line 154 of file MuonTrackValidator.h.

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

Definition at line 154 of file MuonTrackValidator.h.

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

Definition at line 182 of file MuonTrackValidator.h.

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

Definition at line 184 of file MuonTrackValidator.h.

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

Definition at line 183 of file MuonTrackValidator.h.

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

Definition at line 189 of file MuonTrackValidator.h.

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

Definition at line 191 of file MuonTrackValidator.h.

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

Definition at line 190 of file MuonTrackValidator.h.

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

Definition at line 192 of file MuonTrackValidator.h.

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

Definition at line 177 of file MuonTrackValidator.h.

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

Definition at line 179 of file MuonTrackValidator.h.

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

Definition at line 178 of file MuonTrackValidator.h.

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

Definition at line 182 of file MuonTrackValidator.h.

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

Definition at line 184 of file MuonTrackValidator.h.

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

Definition at line 183 of file MuonTrackValidator.h.

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

Definition at line 189 of file MuonTrackValidator.h.

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

Definition at line 191 of file MuonTrackValidator.h.

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

Definition at line 190 of file MuonTrackValidator.h.

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

Definition at line 192 of file MuonTrackValidator.h.

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

Definition at line 177 of file MuonTrackValidator.h.

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

Definition at line 179 of file MuonTrackValidator.h.

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

Definition at line 178 of file MuonTrackValidator.h.

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

Definition at line 158 of file MuonTrackValidator.h.

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

Definition at line 189 of file MuonTrackValidator.h.

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

Definition at line 191 of file MuonTrackValidator.h.

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

Definition at line 190 of file MuonTrackValidator.h.

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

Definition at line 192 of file MuonTrackValidator.h.

bool MuonTrackValidator::MABH
private

Definition at line 151 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

double MuonTrackValidator::maxPhi
private

Definition at line 140 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

double MuonTrackValidator::minPhi
private

Definition at line 140 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

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

Definition at line 159 of file MuonTrackValidator.h.

int MuonTrackValidator::nintPhi
private

Definition at line 141 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

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

Definition at line 168 of file MuonTrackValidator.h.

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

Definition at line 162 of file MuonTrackValidator.h.

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

Definition at line 187 of file MuonTrackValidator.h.

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

Definition at line 188 of file MuonTrackValidator.h.

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

Definition at line 174 of file MuonTrackValidator.h.

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

Definition at line 159 of file MuonTrackValidator.h.

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

Definition at line 175 of file MuonTrackValidator.h.

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

Definition at line 162 of file MuonTrackValidator.h.

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

Definition at line 187 of file MuonTrackValidator.h.

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

Definition at line 188 of file MuonTrackValidator.h.

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

Definition at line 174 of file MuonTrackValidator.h.

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

Definition at line 159 of file MuonTrackValidator.h.

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

Definition at line 175 of file MuonTrackValidator.h.

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

Definition at line 187 of file MuonTrackValidator.h.

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

Definition at line 188 of file MuonTrackValidator.h.

TrackingParticleSelector MuonTrackValidator::tpSelector
private

Definition at line 145 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::UseAssociators
private

Definition at line 139 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::useGsf
private

Definition at line 142 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

Generated for CMSSW Reference Manual by doxygen 1.8.5