CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
MuonTrackValidator.h
Go to the documentation of this file.
1 #ifndef MuonTrackValidator_h
2 #define MuonTrackValidator_h
3 
12 
17 
19  public:
22  dirName_ = pset.getParameter<std::string>("dirName");
23  associatormap = pset.getParameter< edm::InputTag >("associatormap");
24  UseAssociators = pset.getParameter< bool >("UseAssociators");
25  tpSelector = TrackingParticleSelector(pset.getParameter<double>("ptMinTP"),
26  pset.getParameter<double>("minRapidityTP"),
27  pset.getParameter<double>("maxRapidityTP"),
28  pset.getParameter<double>("tipTP"),
29  pset.getParameter<double>("lipTP"),
30  pset.getParameter<int>("minHitTP"),
31  pset.getParameter<bool>("signalOnlyTP"),
32  pset.getParameter<bool>("chargedOnlyTP"),
33  pset.getParameter<bool>("stableOnlyTP"),
34  pset.getParameter<std::vector<int> >("pdgIdTP"));
36  pset.getParameter<double>("minRapidityTP"),
37  pset.getParameter<double>("maxRapidityTP"),
38  pset.getParameter<double>("tipTP"),
39  pset.getParameter<double>("lipTP"),
40  pset.getParameter<int>("minHitTP"),
41  pset.getParameter<bool>("chargedOnlyTP"),
42  pset.getParameter<std::vector<int> >("pdgIdTP"));
43 
44  minPhi = pset.getParameter<double>("minPhi");
45  maxPhi = pset.getParameter<double>("maxPhi");
46  nintPhi = pset.getParameter<int>("nintPhi");
47  useGsf = pset.getParameter<bool>("useGsf");
48  BiDirectional_RecoToSim_association = pset.getParameter<bool>("BiDirectional_RecoToSim_association");
49 
50  // dump cfg parameters
51  edm::LogVerbatim("MuonTrackValidator") << "constructing MuonTrackValidator: " << pset.dump();
52 
53  // Declare consumes (also for the base class)
54  bsSrc_Token = consumes<reco::BeamSpot>(bsSrc);
55  tp_effic_Token = consumes<TrackingParticleCollection>(label_tp_effic);
56  tp_fake_Token = consumes<TrackingParticleCollection>(label_tp_fake);
57  for (unsigned int www=0;www<label.size();www++){
59  }
60  simToRecoCollection_Token = consumes<reco::SimToRecoCollection>(associatormap);
61  recoToSimCollection_Token = consumes<reco::RecoToSimCollection>(associatormap);
62 
63  MABH = false;
64  if (!UseAssociators) {
65  // flag MuonAssociatorByHits
66  if (associators[0] == "MuonAssociationByHits") MABH = true;
67  // reset string associators to the map label
68  associators.clear();
69  associators.push_back(associatormap.label());
70  edm::LogVerbatim("MuonTrackValidator") << "--> associators reset to: " <<associators[0];
71  } else {
72  for (auto const& associator :associators) {
73  consumes<reco::TrackToTrackingParticleAssociator>(edm::InputTag(associator));
74  }
75  }
76 
77  // inform on which SimHits will be counted
78  if (usetracker) edm::LogVerbatim("MuonTrackValidator")
79  <<"\n usetracker = TRUE : Tracker SimHits WILL be counted";
80  else edm::LogVerbatim("MuonTrackValidator")
81  <<"\n usetracker = FALSE : Tracker SimHits WILL NOT be counted";
82  if (usemuon) edm::LogVerbatim("MuonTrackValidator")
83  <<" usemuon = TRUE : Muon SimHits WILL be counted";
84  else edm::LogVerbatim("MuonTrackValidator")
85  <<" usemuon = FALSE : Muon SimHits WILL NOT be counted"<<std::endl;
86 
87  // loop over the reco::Track collections to validate: check for inconsistent input settings
88  for (unsigned int www=0;www<label.size();www++) {
89  std::string recoTracksLabel = label[www].label();
90  std::string recoTracksInstance = label[www].instance();
91 
92  // tracks with hits only on tracker
93  if (recoTracksLabel=="generalTracks" ||
94  (recoTracksLabel.find("cutsRecoTracks") != std::string::npos) ||
95  recoTracksLabel=="ctfWithMaterialTracksP5LHCNavigation" ||
96  recoTracksLabel=="hltL3TkTracksFromL2" ||
97  (recoTracksLabel=="hltL3Muons" && recoTracksInstance=="L2Seeded"))
98  {
99  if (usemuon) {
100  edm::LogWarning("MuonTrackValidator")
101  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
102  <<"\n with usemuon == true"<<"\n ---> please change to usemuon == false ";
103  }
104  if (!usetracker) {
105  edm::LogWarning("MuonTrackValidator")
106  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
107  <<"\n with usetracker == false"<<"\n ---> please change to usetracker == true ";
108  }
109  }
110 
111  // tracks with hits only on muon detectors
112  else if (recoTracksLabel=="standAloneMuons" ||
113  recoTracksLabel=="standAloneSETMuons" ||
114  recoTracksLabel=="cosmicMuons" ||
115  recoTracksLabel=="hltL2Muons")
116  {
117  if (usetracker) {
118  edm::LogWarning("MuonTrackValidator")
119  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
120  <<"\n with usetracker == true"<<"\n ---> please change to usetracker == false ";
121  }
122  if (!usemuon) {
123  edm::LogWarning("MuonTrackValidator")
124  <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
125  <<"\n with usemuon == false"<<"\n ---> please change to usemuon == true ";
126  }
127  }
128 
129  } // for (unsigned int www=0;www<label.size();www++)
130  }
131 
133  virtual ~MuonTrackValidator(){ }
134 
136  // void beginRun(edm::Run const&, edm::EventSetup const&);
138  void analyze(const edm::Event&, const edm::EventSetup& );
140  void endRun(edm::Run const&, edm::EventSetup const&);
142 
143 private:
145  void getRecoMomentum (const reco::Track& track, double& pt, double& ptError,
146  double& qoverp, double& qoverpError, double& lambda, double& lambdaError,
147  double& phi, double& phiError ) const;
149  void getRecoMomentum (const reco::GsfTrack& gsfTrack, double& pt, double& ptError,
150  double& qoverp, double& qoverpError, double& lambda, double& lambdaError,
151  double& phi, double& phiError) const;
152 
153  private:
159  double minPhi, maxPhi;
160  int nintPhi;
161  bool useGsf;
162  // select tracking particles
163  //(i.e. "denominator" of the efficiency ratio)
166 
167  // flag new validation logic (bidirectional RecoToSim association)
169  // flag MuonAssociatorByHits
170  bool MABH;
171 
172  //1D
173  std::vector<MonitorElement*> h_nchi2, h_nchi2_prob, h_losthits;
174 
175  //2D
176  std::vector<MonitorElement*> chi2_vs_nhits, etares_vs_eta;
177  std::vector<MonitorElement*> h_ptshifteta;
178  std::vector<MonitorElement*> ptres_vs_phi, chi2_vs_phi, nhits_vs_phi, phires_vs_phi;
179 
180  //Profile2D
181  std::vector<MonitorElement*> ptmean_vs_eta_phi, phimean_vs_eta_phi;
182 
183  //assoc chi2
184  std::vector<MonitorElement*> h_assochi2, h_assochi2_prob;
185 
186  //chi2 and # lost hits vs eta: to be used with doProfileX
187  std::vector<MonitorElement*> chi2_vs_eta, nlosthits_vs_eta;
188  std::vector<MonitorElement*> h_chi2meanh, h_losthits_eta;
189  std::vector<MonitorElement*> h_hits_phi;
190  std::vector<MonitorElement*> h_chi2meanhitsh, h_chi2mean_vs_phi;
191 
192  //resolution of track params: to be used with fitslicesytool
195 
196  //pulls of track params vs eta: to be used with fitslicesytool
198  std::vector<MonitorElement*> ptpull_vs_phi, phipull_vs_phi, thetapull_vs_phi;
200  std::vector<MonitorElement*> h_ptpullphi, h_phipullphi, h_thetapullphi;
201 
202 };
203 
204 
205 #endif
std::vector< MonitorElement * > h_dxypulleta
T getParameter(std::string const &) const
std::vector< MonitorElement * > h_ptpulleta
std::vector< MonitorElement * > etares_vs_eta
std::vector< MonitorElement * > ptres_vs_phi
std::vector< edm::EDGetTokenT< edm::View< reco::Track > > > track_Collection_Token
std::vector< MonitorElement * > h_losthits_eta
std::string dump(unsigned int indent=0) const
std::vector< MonitorElement * > h_thetapullphi
std::vector< MonitorElement * > h_ptpullphi
std::vector< MonitorElement * > ptres_vs_eta
std::vector< MonitorElement * > cotThetares_vs_eta
std::vector< MonitorElement * > h_dzpulleta
std::vector< MonitorElement * > h_thetapulleta
std::vector< MonitorElement * > phimean_vs_eta_phi
virtual ~MuonTrackValidator()
Destructor.
std::vector< MonitorElement * > phires_vs_pt
std::vector< MonitorElement * > h_assochi2_prob
std::vector< MonitorElement * > cotThetares_vs_pt
std::vector< MonitorElement * > h_ptshifteta
std::vector< MonitorElement * > dxyres_vs_eta
std::vector< MonitorElement * > dzres_vs_pt
std::vector< MonitorElement * > dzres_vs_eta
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) ...
std::vector< MonitorElement * > thetapull_vs_phi
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
std::vector< MonitorElement * > h_assochi2
std::vector< MonitorElement * > chi2_vs_phi
std::vector< MonitorElement * > h_phipulleta
std::vector< MonitorElement * > chi2_vs_nhits
std::vector< MonitorElement * > phires_vs_phi
std::vector< MonitorElement * > h_chi2mean_vs_phi
std::vector< MonitorElement * > h_phipullphi
SingleObjectSelector< TrackingParticleCollection,::TrackingParticleSelector > TrackingParticleSelector
CosmicTrackingParticleSelector cosmictpSelector
std::vector< MonitorElement * > nhits_vs_phi
edm::EDGetTokenT< TrackingParticleCollection > tp_fake_Token
std::vector< MonitorElement * > dxyres_vs_pt
std::vector< edm::InputTag > label
std::vector< MonitorElement * > h_nchi2
std::vector< MonitorElement * > h_nchi2_prob
std::vector< MonitorElement * > ptpull_vs_eta
std::vector< MonitorElement * > dxypull_vs_eta
void analyze(const edm::Event &, const edm::EventSetup &)
Method called before the event loop.
std::vector< MonitorElement * > phipull_vs_eta
std::vector< MonitorElement * > thetapull_vs_eta
MuonTrackValidator(const edm::ParameterSet &pset)
Constructor.
ObjectSelector< CosmicTrackingParticleSelector > CosmicTrackingParticleSelector
std::vector< MonitorElement * > ptmean_vs_eta_phi
edm::EDGetTokenT< reco::RecoToSimCollection > recoToSimCollection_Token
edm::EDGetTokenT< reco::SimToRecoCollection > simToRecoCollection_Token
std::vector< MonitorElement * > chi2_vs_eta
TrackingParticleSelector tpSelector
std::vector< MonitorElement * > h_losthits
std::vector< MonitorElement * > ptpull_vs_phi
edm::InputTag associatormap
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
std::vector< MonitorElement * > ptres_vs_pt
std::vector< MonitorElement * > h_chi2meanhitsh
std::vector< MonitorElement * > phires_vs_eta
std::vector< MonitorElement * > h_hits_phi
std::vector< std::string > associators
std::vector< MonitorElement * > nlosthits_vs_eta
std::vector< MonitorElement * > dzpull_vs_eta
edm::EDGetTokenT< TrackingParticleCollection > tp_effic_Token
Definition: Run.h:41
std::vector< MonitorElement * > h_chi2meanh
void endRun(edm::Run const &, edm::EventSetup const &)
Method called at the end of the event loop.
std::vector< MonitorElement * > phipull_vs_phi
edm::EDGetTokenT< reco::BeamSpot > bsSrc_Token
Definition: DDAxes.h:10