CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HeavyFlavorValidation.cc
Go to the documentation of this file.
1 
8 // Original Author: Zoltan Gecse
9 
10 #include <memory>
11 #include "FWCore/Framework/interface/Frameworkfwd.h"
12 #include "FWCore/Framework/interface/EDAnalyzer.h"
13 #include "FWCore/Framework/interface/Event.h"
14 #include "FWCore/Framework/interface/MakerMacros.h"
15 #include "FWCore/MessageLogger/interface/MessageLogger.h"
16 #include "FWCore/ServiceRegistry/interface/Service.h"
17 #include "FWCore/ParameterSet/interface/ParameterSet.h"
18 #include "FWCore/Utilities/interface/InputTag.h"
19 #include "FWCore/Common/interface/TriggerNames.h"
20 
21 #include "DataFormats/Common/interface/RefToBase.h"
22 #include "DataFormats/TrackReco/interface/Track.h"
23 #include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"
24 #include "DataFormats/RecoCandidate/interface/RecoChargedCandidateFwd.h"
25 #include "DataFormats/L1Trigger/interface/L1MuonParticle.h"
26 #include "DataFormats/L1Trigger/interface/L1MuonParticleFwd.h"
27 #include "DataFormats/Candidate/interface/Candidate.h"
28 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
29 #include "DataFormats/MuonReco/interface/MuonFwd.h"
30 #include "DataFormats/MuonReco/interface/Muon.h"
31 #include "DataFormats/Math/interface/deltaR.h"
32 #include "DataFormats/Common/interface/TriggerResults.h"
33 #include "DataFormats/Candidate/interface/Particle.h"
34 #include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
35 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
36 
37 #include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
38 
39 #include "DQMServices/Core/interface/DQMStore.h"
40 #include "DQMServices/Core/interface/MonitorElement.h"
41 
42 #include "CommonTools/Utils/interface/PtComparator.h"
43 
44 #include "TLorentzVector.h"
45 
46 using namespace std;
47 using namespace edm;
48 using namespace reco;
49 using namespace l1extra;
50 using namespace trigger;
51 
52 class HeavyFlavorValidation : public edm::EDAnalyzer {
53  public:
54  explicit HeavyFlavorValidation(const edm::ParameterSet&);
55  ~HeavyFlavorValidation();
56  private:
57  virtual void beginRun(const Run & iRun, const EventSetup & iSetup) override;
58  virtual void analyze(const edm::Event&, const edm::EventSetup&) override;
59  virtual void endJob() override;
60  int getMotherId( const Candidate * p );
61  void match( MonitorElement * me, vector<LeafCandidate> & from, vector<LeafCandidate> & to, double deltaRMatchingCut, vector<int> & map );
62  void myBook2D( TString name, vector<double> &xBins, TString xLabel, vector<double> &yBins, TString yLabel, TString title);
63  void myBook2D( TString name, vector<double> &xBins, TString xLabel, vector<double> &yBins, TString yLabel){
64  myBook2D( name, xBins, xLabel, yBins, yLabel, name);
65  }
66  void myBookProfile2D( TString name, vector<double> &xBins, TString xLabel, vector<double> &yBins, TString yLabel, TString title);
67  void myBookProfile2D( TString name, vector<double> &xBins, TString xLabel, vector<double> &yBins, TString yLabel){
68  myBookProfile2D( name, xBins, xLabel, yBins, yLabel, name);
69  }
70  void myBook1D( TString name, vector<double> &xBins, TString label, TString title );
71  void myBook1D( TString name, vector<double> &xBins, TString label ){
72  myBook1D( name, xBins, label, name );
73  }
74  string dqmFolder;
75  string triggerProcessName;
76  string triggerPathName;
77 
78  EDGetTokenT<TriggerEventWithRefs> triggerSummaryRAWTag;
79  EDGetTokenT<TriggerEvent> triggerSummaryAODTag;
80  EDGetTokenT<TriggerResults> triggerResultsTag;
81  EDGetTokenT<MuonCollection> recoMuonsTag;
82  EDGetTokenT<GenParticleCollection> genParticlesTag;
83 
84  vector<int> motherIDs;
85  double genGlobDeltaRMatchingCut;
86  double globL1DeltaRMatchingCut;
87  double globL2DeltaRMatchingCut;
88  double globL3DeltaRMatchingCut;
89  vector<double> deltaEtaBins;
90  vector<double> deltaPhiBins;
91  vector<double> muonPtBins;
92  vector<double> muonEtaBins;
93  vector<double> muonPhiBins;
94  vector<double> dimuonPtBins;
95  vector<double> dimuonEtaBins;
96  vector<double> dimuonDRBins;
97  DQMStore* dqmStore;
98  map<TString, MonitorElement *> ME;
99  vector<pair<string,int> > filterNamesLevels;
100  const double muonMass;
101 };
102 
103 HeavyFlavorValidation::HeavyFlavorValidation(const ParameterSet& pset):
104 //get parameters
105  dqmFolder(pset.getUntrackedParameter<string>("DQMFolder")),
106  triggerProcessName(pset.getUntrackedParameter<string>("TriggerProcessName")),
107  triggerPathName(pset.getUntrackedParameter<string>("TriggerPathName")),
108  motherIDs(pset.getUntrackedParameter<vector<int> >("MotherIDs")),
109  genGlobDeltaRMatchingCut(pset.getUntrackedParameter<double>("GenGlobDeltaRMatchingCut")),
110  globL1DeltaRMatchingCut(pset.getUntrackedParameter<double>("GlobL1DeltaRMatchingCut")),
111  globL2DeltaRMatchingCut(pset.getUntrackedParameter<double>("GlobL2DeltaRMatchingCut")),
112  globL3DeltaRMatchingCut(pset.getUntrackedParameter<double>("GlobL3DeltaRMatchingCut")),
113  deltaEtaBins(pset.getUntrackedParameter<vector<double> >("DeltaEtaBins")),
114  deltaPhiBins(pset.getUntrackedParameter<vector<double> >("DeltaPhiBins")),
115  muonPtBins(pset.getUntrackedParameter<vector<double> >("MuonPtBins")),
116  muonEtaBins(pset.getUntrackedParameter<vector<double> >("MuonEtaBins")),
117  muonPhiBins(pset.getUntrackedParameter<vector<double> >("MuonPhiBins")),
118  dimuonPtBins(pset.getUntrackedParameter<vector<double> >("DimuonPtBins")),
119  dimuonEtaBins(pset.getUntrackedParameter<vector<double> >("DimuonEtaBins")),
120  dimuonDRBins(pset.getUntrackedParameter<vector<double> >("DimuonDRBins")),
121  muonMass(0.106)
122 {
123  triggerSummaryRAWTag = consumes<TriggerEventWithRefs>(InputTag( pset.getUntrackedParameter<string>("TriggerSummaryRAW"), "", triggerProcessName));
124  triggerSummaryAODTag = consumes<TriggerEvent>(InputTag( pset.getUntrackedParameter<string>("TriggerSummaryAOD"), "", triggerProcessName));
125  triggerResultsTag = consumes<TriggerResults>(InputTag( pset.getUntrackedParameter<string>("TriggerResults"), "", triggerProcessName));
126  recoMuonsTag = consumes<MuonCollection>(pset.getParameter<InputTag>("RecoMuons"));
127  genParticlesTag = consumes<GenParticleCollection>(pset.getParameter<InputTag>("GenParticles"));
128 }
129 
130 void HeavyFlavorValidation::beginRun(const Run & iRun, const EventSetup & iSetup){
131 //discover HLT configuration
132  HLTConfigProvider hltConfig;
133  bool isChanged;
134  if(hltConfig.init(iRun, iSetup, triggerProcessName, isChanged)){
135  LogDebug("HLTriggerOfflineHeavyFlavor") << "Successfully initialized HLTConfigProvider with process name: "<<triggerProcessName<<endl;
136  }else{
137  LogWarning("HLTriggerOfflineHeavyFlavor") << "Could not initialize HLTConfigProvider with process name: "<<triggerProcessName<<endl;
138  return;
139  }
140  stringstream os;
141  vector<string> triggerNames = hltConfig.triggerNames();
142  for( size_t i = 0; i < triggerNames.size(); i++) {
143  TString triggerName = triggerNames[i];
144  if (triggerName.Contains(triggerPathName)){
145  vector<string> moduleNames = hltConfig.moduleLabels( triggerNames[i] );
146  for( size_t j = 0; j < moduleNames.size(); j++) {
147  TString name = moduleNames[j];
148  if(name.Contains("Filter")){
149  int level = 0;
150  if(name.Contains("L1"))
151  level = 1;
152  else if(name.Contains("L2"))
153  level = 2;
154  else if(name.Contains("L3"))
155  level = 3;
156  else if(name.Contains("mumuFilter") || name.Contains("JpsiTrackMass"))
157  level = 4;
158  filterNamesLevels.push_back( pair<string,int>(moduleNames[j],level) );
159  os<<" "<<moduleNames[j];
160  }
161  }
162  break;
163  }
164  }
165  if(filterNamesLevels.size()==0){
166  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Bad Trigger Path: "<<triggerPathName<<endl;
167  return;
168  }else{
169  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Trigger Path: "<<triggerPathName<<" has filters:"<<os.str();
170  }
171 
172 //create Monitor Elements
173  dqmStore = Service<DQMStore>().operator->();
174  if( !dqmStore ){
175  LogError("HLTriggerOfflineHeavyFlavor") << "Could not find DQMStore service\n";
176  return;
177  }
178  dqmStore->setVerbose(0);
179  dqmStore->setCurrentFolder((dqmFolder+"/")+triggerProcessName+"/"+triggerPathName);
180 // Eta Pt Single
181  myBook2D( "genMuon_genEtaPt", muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)");
182  myBook2D( "globMuon_genEtaPt", muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)");
183  myBook2D( "globMuon_recoEtaPt", muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)");
184  for(size_t i=0; i<filterNamesLevels.size(); i++){
185  myBook2D( TString::Format("filt%dMuon_recoEtaPt",int(i+1)), muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)", filterNamesLevels[i].first);
186  }
187  myBook2D( "pathMuon_recoEtaPt", muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)", triggerPathName);
188  myBook2D( "resultMuon_recoEtaPt", muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)");
189 // Eta Pt Single Resolution
190  myBookProfile2D( "resGlobGen_genEtaPt", muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)");
191  for(size_t i=0; i<filterNamesLevels.size(); i++){
192  myBookProfile2D( TString::Format("resFilt%dGlob_recoEtaPt",int(i+1)), muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)", filterNamesLevels[i].first);
193  }
194  myBookProfile2D( "resPathGlob_recoEtaPt", muonEtaBins, "#mu eta", muonPtBins, " #mu pT (GeV)", triggerPathName);
195 // Eta Pt Double
196  myBook2D( "genDimuon_genEtaPt", dimuonEtaBins, "#mu#mu eta", dimuonPtBins, " #mu#mu pT (GeV)");
197  myBook2D( "globDimuon_genEtaPt", dimuonEtaBins, "#mu#mu eta", dimuonPtBins, " #mu#mu pT (GeV)");
198  myBook2D( "globDimuon_recoEtaPt", dimuonEtaBins, "#mu#mu eta", dimuonPtBins, " #mu#mu pT (GeV)");
199  for(size_t i=0; i<filterNamesLevels.size(); i++){
200  myBook2D( TString::Format("filt%dDimuon_recoEtaPt",int(i+1)), dimuonEtaBins, "#mu#mu eta", dimuonPtBins, " #mu#mu pT (GeV)", filterNamesLevels[i].first);
201  }
202  myBook2D( "pathDimuon_recoEtaPt", dimuonEtaBins, "#mu#mu eta", dimuonPtBins, " #mu#mu pT (GeV)", triggerPathName);
203  myBook2D( "resultDimuon_recoEtaPt", dimuonEtaBins, "#mu#mu eta", dimuonPtBins, " #mu#mu pT (GeV)");
204  for(size_t i=0; i<filterNamesLevels.size(); i++){
205  myBook2D( TString::Format("diFilt%dDimuon_recoEtaPt",int(i+1)), dimuonEtaBins, "#mu#mu eta", dimuonPtBins, " #mu#mu pT (GeV)", filterNamesLevels[i].first);
206  }
207  myBook2D( "diPathDimuon_recoEtaPt", dimuonEtaBins, "#mu#mu eta", dimuonPtBins, " #mu#mu pT (GeV)", triggerPathName);
208 // Eta Phi Single
209  myBook2D( "genMuon_genEtaPhi", muonEtaBins, "#mu eta", muonPhiBins, "#mu phi");
210  myBook2D( "globMuon_genEtaPhi", muonEtaBins, "#mu eta", muonPhiBins, "#mu phi");
211  myBook2D( "globMuon_recoEtaPhi", muonEtaBins, "#mu eta", muonPhiBins, "#mu phi");
212  for(size_t i=0; i<filterNamesLevels.size(); i++){
213  myBook2D( TString::Format("filt%dMuon_recoEtaPhi",int(i+1)), muonEtaBins, "#mu eta", muonPhiBins, "#mu phi", filterNamesLevels[i].first);
214  }
215  myBook2D( "pathMuon_recoEtaPhi", muonEtaBins, "#mu eta", muonPhiBins, "#mu phi", triggerPathName);
216  myBook2D( "resultMuon_recoEtaPhi", muonEtaBins, "#mu eta", muonPhiBins, "#mu phi");
217 // Rap Pt Double
218  myBook2D( "genDimuon_genRapPt", dimuonEtaBins, "#mu#mu rapidity", dimuonPtBins, " #mu#mu pT (GeV)");
219  myBook2D( "globDimuon_genRapPt", dimuonEtaBins, "#mu#mu rapidity", dimuonPtBins, " #mu#mu pT (GeV)");
220  myBook2D( "globDimuon_recoRapPt", dimuonEtaBins, "#mu#mu rapidity", dimuonPtBins, " #mu#mu pT (GeV)");
221  for(size_t i=0; i<filterNamesLevels.size(); i++){
222  myBook2D( TString::Format("filt%dDimuon_recoRapPt",int(i+1)), dimuonEtaBins, "#mu#mu rapidity", dimuonPtBins, " #mu#mu pT (GeV)", filterNamesLevels[i].first);
223  }
224  myBook2D( "pathDimuon_recoRapPt", dimuonEtaBins, "#mu#mu rapidity", dimuonPtBins, " #mu#mu pT (GeV)", triggerPathName);
225  myBook2D( "resultDimuon_recoRapPt", dimuonEtaBins, "#mu#mu rapidity", dimuonPtBins, " #mu#mu pT (GeV)");
226  for(size_t i=0; i<filterNamesLevels.size(); i++){
227  myBook2D( TString::Format("diFilt%dDimuon_recoRapPt",int(i+1)), dimuonEtaBins, "#mu#mu rapidity", dimuonPtBins, " #mu#mu pT (GeV)", filterNamesLevels[i].first);
228  }
229  myBook2D( "diPathDimuon_recoRapPt", dimuonEtaBins, "#mu#mu rapidity", dimuonPtBins, " #mu#mu pT (GeV)", triggerPathName);
230 // Pt DR Double
231  myBook2D( "genDimuon_genPtDR", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R at IP");
232  myBook2D( "globDimuon_genPtDR", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R at IP");
233  myBook2D( "globDimuon_recoPtDR", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R at IP");
234  for(size_t i=0; i<filterNamesLevels.size(); i++){
235  myBook2D( TString::Format("filt%dDimuon_recoPtDR",int(i+1)), dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R at IP", filterNamesLevels[i].first);
236  }
237  myBook2D( "pathDimuon_recoPtDR", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R at IP", triggerPathName);
238  for(size_t i=0; i<filterNamesLevels.size(); i++){
239  myBook2D( TString::Format("diFilt%dDimuon_recoPtDR",int(i+1)), dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R at IP", filterNamesLevels[i].first);
240  }
241  myBook2D( "diPathDimuon_recoPtDR", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R at IP", triggerPathName);
242  myBook2D( "resultDimuon_recoPtDR", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R at IP");
243 // Pt DRpos Double
244 // myBook2D( "genDimuon_genPtDRpos", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R in MS");
245 // myBook2D( "globDimuon_genPtDRpos", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R in MS");
246  myBook2D( "globDimuon_recoPtDRpos", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R in MS");
247  for(size_t i=0; i<filterNamesLevels.size(); i++){
248  myBook2D( TString::Format("filt%dDimuon_recoPtDRpos",int(i+1)), dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R in MS", filterNamesLevels[i].first);
249  }
250  myBook2D( "pathDimuon_recoPtDRpos", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R in MS", triggerPathName);
251  for(size_t i=0; i<filterNamesLevels.size(); i++){
252  myBook2D( TString::Format("diFilt%dDimuon_recoPtDRpos",int(i+1)), dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R in MS", filterNamesLevels[i].first);
253  }
254  myBook2D( "diPathDimuon_recoPtDRpos", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R in MS", triggerPathName);
255  myBook2D( "resultDimuon_recoPtDRpos", dimuonPtBins, " #mu#mu pT (GeV)", dimuonDRBins, "#mu#mu #Delta R in MS");
256 
257 /* ME["genDimuon_genPt+Pt-"] = myBook2D( "genDimuon_genPt+Pt-", muonPtBins, "#mu+ pt (GeV)", muonPtBins, "#mu- pT (GeV)");
258  ME["genGlobDimuon_genPt+Pt-"] = myBook2D( "genGlobDimuon_genPt+Pt-", muonPtBins, "#mu+ pt (GeV)", muonPtBins, "#mu- pT (GeV)");
259  ME["genGlobDimuon_recoPt+Pt-"] = myBook2D( "genGlobDimuon_recoPt+Pt-", muonPtBins, "#mu+ pt (GeV)", muonPtBins, "#mu- pT (GeV)");
260  ME["genGlobL1Dimuon_recoPt+Pt-"] = myBook2D( "genGlobL1Dimuon_recoPt+Pt-", muonPtBins, "#mu+ pt (GeV)", muonPtBins, "#mu- pT (GeV)");
261  ME["genGlobL1L2Dimuon_recoPt+Pt-"] = myBook2D( "genGlobL1L2Dimuon_recoPt+Pt-", muonPtBins, "#mu+ pt (GeV)", muonPtBins, "#mu- pT (GeV)");
262  ME["genGlobL1L2L3Dimuon_recoPt+Pt-"] = myBook2D( "genGlobL1L2L3Dimuon_recoPt+Pt-", muonPtBins, "#mu+ pt (GeV)", muonPtBins, "#mu- pT (GeV)");*/
263 
264 // Matching
265  myBook2D( "globGen_deltaEtaDeltaPhi", deltaEtaBins, "#Delta eta", deltaPhiBins, "#Delta phi");
266  for(size_t i=0; i<filterNamesLevels.size(); i++){
267  myBook2D( TString::Format("filt%dGlob_deltaEtaDeltaPhi",int(i+1)), deltaEtaBins, "#Delta eta", deltaPhiBins, "#Delta phi", filterNamesLevels[i].first);
268  }
269  myBook2D( "pathGlob_deltaEtaDeltaPhi", deltaEtaBins, "#Delta eta", deltaPhiBins, "#Delta phi", triggerPathName);
270 // Size of containers
271  vector<double> sizeBins; sizeBins.push_back(10); sizeBins.push_back(0); sizeBins.push_back(10);
272  myBook1D( "genMuon_size", sizeBins, "container size" );
273  myBook1D( "globMuon_size", sizeBins, "container size" );
274  for(size_t i=0; i<filterNamesLevels.size(); i++){
275  myBook1D( TString::Format("filt%dMuon_size",int(i+1)), sizeBins, "container size", filterNamesLevels[i].first);
276  }
277  myBook1D( "pathMuon_size", sizeBins, "container size", triggerPathName );
278 }
279 
280 void HeavyFlavorValidation::analyze(const Event& iEvent, const EventSetup& iSetup){
281  if( !dqmStore ){
282  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Could not access DQM Store service"<<endl;
283  return;
284  }
285  if(filterNamesLevels.size()==0){
286  return;
287  }
288 //access the containers and create LeafCandidate copies
289  vector<LeafCandidate> genMuons;
290  Handle<GenParticleCollection> genParticles;
291  iEvent.getByToken(genParticlesTag, genParticles);
292  if(genParticles.isValid()){
293  for(GenParticleCollection::const_iterator p=genParticles->begin(); p!= genParticles->end(); ++p){
294  if( p->status() == 1 && std::abs(p->pdgId())==13 &&
295  ( find( motherIDs.begin(), motherIDs.end(), -1 )!=motherIDs.end() || find( motherIDs.begin(), motherIDs.end(), getMotherId(&(*p)) )!=motherIDs.end() ) ){
296  genMuons.push_back( *p );
297  }
298  }
299  }else{
300  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Could not access GenParticleCollection"<<endl;
301  }
302  sort(genMuons.begin(), genMuons.end(), GreaterByPt<LeafCandidate>());
303  ME["genMuon_size"]->Fill(genMuons.size());
304  LogDebug("HLTriggerOfflineHeavyFlavor")<<"GenParticleCollection from "<<genParticlesTag<<" has size: "<<genMuons.size()<<endl;
305 
306  vector<LeafCandidate> globMuons;
307  vector<LeafCandidate> globMuons_position;
308  Handle<MuonCollection> recoMuonsHandle;
309  iEvent.getByToken(recoMuonsTag, recoMuonsHandle);
310  if(recoMuonsHandle.isValid()){
311  for(MuonCollection::const_iterator p=recoMuonsHandle->begin(); p!= recoMuonsHandle->end(); ++p){
312  if(p->isGlobalMuon()){
313  globMuons.push_back(*p);
314  globMuons_position.push_back( LeafCandidate( p->charge(), math::XYZTLorentzVector(p->outerTrack()->innerPosition().x(), p->outerTrack()->innerPosition().y(), p->outerTrack()->innerPosition().z(), 0.) ) );
315  }
316  }
317  }else{
318  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Could not access reco Muons"<<endl;
319  }
320  ME["globMuon_size"]->Fill(globMuons.size());
321  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Global Muons from "<<recoMuonsTag<<" has size: "<<globMuons.size()<<endl;
322 
323 // access RAW trigger event
324  vector<vector<LeafCandidate> > muonsAtFilter;
325  vector<vector<LeafCandidate> > muonPositionsAtFilter;
326  for(size_t i=0; i<filterNamesLevels.size(); i++){
327  muonsAtFilter.push_back(vector<LeafCandidate>());
328  muonPositionsAtFilter.push_back(vector<LeafCandidate>());
329  }
330  Handle<TriggerEventWithRefs> rawTriggerEvent;
331  iEvent.getByToken( triggerSummaryRAWTag, rawTriggerEvent );
332  if( rawTriggerEvent.isValid() ){
333  for(size_t i=0; i<filterNamesLevels.size(); i++){
334  size_t index = rawTriggerEvent->filterIndex(InputTag( filterNamesLevels[i].first, "", triggerProcessName ));
335  if ( index < rawTriggerEvent->size() ){
336  if( filterNamesLevels[i].second==1 ){
337  vector<L1MuonParticleRef> l1Cands;
338  rawTriggerEvent->getObjects( index, TriggerL1Mu, l1Cands );
339  for(size_t j=0; j<l1Cands.size(); j++){
340  muonsAtFilter[i].push_back(*l1Cands[j]);
341  }
342  }else{
343  vector<RecoChargedCandidateRef> hltCands;
344  rawTriggerEvent->getObjects( index, TriggerMuon, hltCands );
345  for(size_t j=0; j<hltCands.size(); j++){
346  muonsAtFilter[i].push_back(*hltCands[j]);
347  if( filterNamesLevels[i].second==2 ){
348  muonPositionsAtFilter[i].push_back( LeafCandidate( hltCands[j]->charge(), math::XYZTLorentzVector(hltCands[j]->track()->innerPosition().x(), hltCands[j]->track()->innerPosition().y(), hltCands[j]->track()->innerPosition().z(), 0.) ) );
349  }
350  }
351  }
352  }
353  ME[TString::Format("filt%dMuon_size",int(i+1))]->Fill(muonsAtFilter[i].size());
354  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Filter \""<<filterNamesLevels[i].first<<"\" has "<<muonsAtFilter[i].size()<<" muons"<<endl;
355  }
356  }else{
357  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Could not access RAWTriggerEvent"<<endl;
358  }
359 
360 // access AOD trigger event
361  vector<LeafCandidate> pathMuons;
362  Handle<TriggerEvent> aodTriggerEvent;
363  iEvent.getByToken(triggerSummaryAODTag,aodTriggerEvent);
364  if(aodTriggerEvent.isValid()){
365  TriggerObjectCollection allObjects = aodTriggerEvent->getObjects();
366  for(int i=0; i<aodTriggerEvent->sizeFilters(); i++){
367  if(aodTriggerEvent->filterTag(i)==InputTag((filterNamesLevels.end()-1)->first,"",triggerProcessName)){
368  Keys keys = aodTriggerEvent->filterKeys(i);
369  for(size_t j=0; j<keys.size(); j++){
370  pathMuons.push_back( LeafCandidate( allObjects[keys[j]].id()>0 ? 1 : -1, math::PtEtaPhiMLorentzVector( allObjects[keys[j]].pt(), allObjects[keys[j]].eta(), allObjects[keys[j]].phi(), muonMass ) ) );
371  }
372  }
373  }
374  ME["pathMuon_size"]->Fill(pathMuons.size());
375  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Path \""<<triggerPathName<<"\" has "<<pathMuons.size()<<" muons at last filter \""<<(filterNamesLevels.end()-1)->first<<"\""<<endl;
376  }else{
377  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Could not access AODTriggerEvent"<<endl;
378  }
379 
380 // access Trigger Results
381  bool triggerFired = false;
382  Handle<TriggerResults> triggerResults;
383  iEvent.getByToken(triggerResultsTag,triggerResults);
384  if(triggerResults.isValid()){
385  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Successfully initialized "<<triggerResultsTag<<endl;
386  const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
387  bool hlt_exists = false;
388  for ( unsigned int i=0; i!=triggerNames.size(); i++) {
389  TString hlt_name = triggerNames.triggerName(i);
390  if (hlt_name.Contains(triggerPathName)) {
391  triggerFired = triggerResults->accept( i );
392  hlt_exists = true;
393  break;
394  }
395  }
396  if (!hlt_exists) {
397  LogDebug("HLTriggerOfflineHeavyFlavor")<<triggerResultsTag<<" has no trigger: "<<triggerPathName<<endl;
398  }
399  }else{
400  LogDebug("HLTriggerOfflineHeavyFlavor")<<"Could not initialize "<<triggerResultsTag<<endl;
401  }
402 
403 //create matching maps
404  vector<int> glob_gen(genMuons.size(),-1);
405  match( ME["globGen_deltaEtaDeltaPhi"], genMuons, globMuons, genGlobDeltaRMatchingCut, glob_gen );
406  vector<vector<int> > filt_glob;
407  for(size_t i=0; i<filterNamesLevels.size(); i++){
408  filt_glob.push_back( vector<int>(globMuons.size(),-1) );
409  if( filterNamesLevels[i].second == 1 ){
410  match( ME[TString::Format("filt%dGlob_deltaEtaDeltaPhi",int(i+1))], globMuons_position, muonsAtFilter[i] ,globL1DeltaRMatchingCut, filt_glob[i] );
411  }else if( filterNamesLevels[i].second == 2 ){
412  match( ME[TString::Format("filt%dGlob_deltaEtaDeltaPhi",int(i+1))], globMuons_position, muonPositionsAtFilter[i] ,globL2DeltaRMatchingCut, filt_glob[i] );
413  }else if( filterNamesLevels[i].second == 3 || filterNamesLevels[i].second == 4){
414  match( ME[TString::Format("filt%dGlob_deltaEtaDeltaPhi",int(i+1))], globMuons, muonsAtFilter[i] ,globL3DeltaRMatchingCut, filt_glob[i] );
415  }
416  }
417  vector<int> path_glob(globMuons.size(),-1);
418  if( (filterNamesLevels.end()-1)->second == 1 ){
419  match( ME["pathGlob_deltaEtaDeltaPhi"], globMuons_position, pathMuons ,globL1DeltaRMatchingCut, path_glob );
420  }else if( (filterNamesLevels.end()-1)->second == 2 ){
421  match( ME["pathGlob_deltaEtaDeltaPhi"], globMuons, pathMuons ,globL2DeltaRMatchingCut, path_glob );
422  }else if( (filterNamesLevels.end()-1)->second == 3 || (filterNamesLevels.end()-1)->second == 4){
423  match( ME["pathGlob_deltaEtaDeltaPhi"], globMuons, pathMuons ,globL3DeltaRMatchingCut, path_glob );
424  }
425 
426 //fill histos
427  bool first = true;
428  for(size_t i=0; i<genMuons.size(); i++){
429  ME["genMuon_genEtaPt"]->Fill(genMuons[i].eta(), genMuons[i].pt());
430  ME["genMuon_genEtaPhi"]->Fill(genMuons[i].eta(), genMuons[i].phi());
431  if(glob_gen[i] != -1){
432  ME["resGlobGen_genEtaPt"]->Fill(genMuons[i].eta(), genMuons[i].pt(), (globMuons[glob_gen[i]].pt()-genMuons[i].pt())/genMuons[i].pt() );
433  ME["globMuon_genEtaPt"]->Fill(genMuons[i].eta(), genMuons[i].pt());
434  ME["globMuon_genEtaPhi"]->Fill(genMuons[i].eta(), genMuons[i].phi());
435  ME["globMuon_recoEtaPt"]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].pt());
436  ME["globMuon_recoEtaPhi"]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].phi());
437  for(size_t f=0; f<filterNamesLevels.size(); f++){
438  if(filt_glob[f][glob_gen[i]] != -1){
439  ME[TString::Format("resFilt%dGlob_recoEtaPt",int(f+1))]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].pt(), (muonsAtFilter[f][filt_glob[f][glob_gen[i]]].pt()-globMuons[glob_gen[i]].pt())/globMuons[glob_gen[i]].pt() );
440  ME[TString::Format("filt%dMuon_recoEtaPt",int(f+1))]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].pt());
441  ME[TString::Format("filt%dMuon_recoEtaPhi",int(f+1))]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].phi());
442  }else{
443  break;
444  }
445  }
446  if(path_glob[glob_gen[i]] != -1){
447  ME["resPathGlob_recoEtaPt"]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].pt(), (pathMuons[path_glob[glob_gen[i]]].pt()-globMuons[glob_gen[i]].pt())/globMuons[glob_gen[i]].pt() );
448  ME["pathMuon_recoEtaPt"]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].pt());
449  ME["pathMuon_recoEtaPhi"]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].phi());
450  }
451 //highest pt muon
452  if( first ){
453  first = false;
454  if( triggerFired ){
455  ME["resultMuon_recoEtaPt"]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].pt());
456  ME["resultMuon_recoEtaPhi"]->Fill(globMuons[glob_gen[i]].eta(), globMuons[glob_gen[i]].phi());
457  }
458  }
459  }
460  }
461 
462 //fill dimuon histograms (highest pT, opposite charge)
463  int secondMuon = 0;
464  for(size_t j=1; j<genMuons.size(); j++){
465  if(genMuons[0].charge()*genMuons[j].charge()==-1){
466  secondMuon = j;
467  break;
468  }
469  }
470  if(secondMuon > 0){
471 // int pos = genMuons[0].charge()>0 ? 0 : secondMuon ;
472 // int neg = genMuons[0].charge()<0 ? 0 : secondMuon ;
473 //two generated
474  double genDimuonPt = (genMuons[0].p4()+genMuons[secondMuon].p4()).pt();
475  double genDimuonEta = (genMuons[0].p4()+genMuons[secondMuon].p4()).eta();
476  double genDimuonRap = (genMuons[0].p4()+genMuons[secondMuon].p4()).Rapidity();
477  double genDimuonDR = deltaR<LeafCandidate,LeafCandidate>( genMuons[0], genMuons[secondMuon] );
478  bool highPt = genMuons[0].pt()>7. && genMuons[secondMuon].pt()>7;
479  ME["genDimuon_genEtaPt"]->Fill( genDimuonEta, genDimuonPt );
480  ME["genDimuon_genRapPt"]->Fill( genDimuonRap, genDimuonPt );
481  if(highPt) ME["genDimuon_genPtDR"]->Fill( genDimuonPt, genDimuonDR );
482 //two global
483  if(glob_gen[0]!=-1 && glob_gen[secondMuon]!=-1){
484  ME["globDimuon_genEtaPt"]->Fill( genDimuonEta, genDimuonPt );
485  ME["globDimuon_genRapPt"]->Fill( genDimuonRap, genDimuonPt );
486  if(highPt) ME["globDimuon_genPtDR"]->Fill( genDimuonPt, genDimuonDR );
487  double globDimuonPt = (globMuons[glob_gen[0]].p4()+globMuons[glob_gen[secondMuon]].p4()).pt();
488  double globDimuonEta = (globMuons[glob_gen[0]].p4()+globMuons[glob_gen[secondMuon]].p4()).eta();
489  double globDimuonRap = (globMuons[glob_gen[0]].p4()+globMuons[glob_gen[secondMuon]].p4()).Rapidity();
490  double globDimuonDR = deltaR<LeafCandidate,LeafCandidate>( globMuons[glob_gen[0]], globMuons[glob_gen[secondMuon]] );
491  double globDimuonDRpos = deltaR<LeafCandidate,LeafCandidate>( globMuons_position[glob_gen[0]], globMuons_position[glob_gen[secondMuon]] );
492  ME["globDimuon_recoEtaPt"]->Fill( globDimuonEta, globDimuonPt );
493  ME["globDimuon_recoRapPt"]->Fill( globDimuonRap, globDimuonPt );
494  if(highPt) ME["globDimuon_recoPtDR"]->Fill( globDimuonPt, globDimuonDR );
495  if(highPt) ME["globDimuon_recoPtDRpos"]->Fill( globDimuonPt, globDimuonDRpos );
496 //two filter objects
497  for(size_t f=0; f<filterNamesLevels.size(); f++){
498  if(filt_glob[f][glob_gen[0]] != -1 && filt_glob[f][glob_gen[secondMuon]] != -1){
499  ME[TString::Format("diFilt%dDimuon_recoEtaPt",int(f+1))]->Fill( globDimuonEta, globDimuonPt );
500  ME[TString::Format("diFilt%dDimuon_recoRapPt",int(f+1))]->Fill( globDimuonRap, globDimuonPt );
501  if(highPt) ME[TString::Format("diFilt%dDimuon_recoPtDR",int(f+1))]->Fill( globDimuonPt, globDimuonDR );
502  if(highPt) ME[TString::Format("diFilt%dDimuon_recoPtDRpos",int(f+1))]->Fill( globDimuonPt, globDimuonDRpos );
503  }else{
504  break;
505  }
506  }
507 //one filter object
508  for(size_t f=0; f<filterNamesLevels.size(); f++){
509  if(filt_glob[f][glob_gen[0]] != -1 || filt_glob[f][glob_gen[secondMuon]] != -1){
510  ME[TString::Format("filt%dDimuon_recoEtaPt",int(f+1))]->Fill( globDimuonEta, globDimuonPt );
511  ME[TString::Format("filt%dDimuon_recoRapPt",int(f+1))]->Fill( globDimuonRap, globDimuonPt );
512  if(highPt) ME[TString::Format("filt%dDimuon_recoPtDR",int(f+1))]->Fill( globDimuonPt, globDimuonDR );
513  if(highPt) ME[TString::Format("filt%dDimuon_recoPtDRpos",int(f+1))]->Fill( globDimuonPt, globDimuonDRpos );
514  }else{
515  break;
516  }
517  }
518 //two path objects
519  if(path_glob[glob_gen[0]] != -1 && path_glob[glob_gen[secondMuon]] != -1){
520  ME["diPathDimuon_recoEtaPt"]->Fill( globDimuonEta, globDimuonPt );
521  ME["diPathDimuon_recoRapPt"]->Fill( globDimuonRap, globDimuonPt );
522  if(highPt) ME["diPathDimuon_recoPtDR"]->Fill( globDimuonPt, globDimuonDR );
523  if(highPt) ME["diPathDimuon_recoPtDRpos"]->Fill( globDimuonPt, globDimuonDRpos );
524  }
525 //one path object
526  if(path_glob[glob_gen[0]] != -1 || path_glob[glob_gen[secondMuon]] != -1){
527  ME["pathDimuon_recoEtaPt"]->Fill( globDimuonEta, globDimuonPt );
528  ME["pathDimuon_recoRapPt"]->Fill( globDimuonRap, globDimuonPt );
529  if(highPt) ME["pathDimuon_recoPtDR"]->Fill( globDimuonPt, globDimuonDR );
530  if(highPt) ME["pathDimuon_recoPtDRpos"]->Fill( globDimuonPt, globDimuonDRpos );
531  }
532 //trigger result
533  if( triggerFired ){
534  ME["resultDimuon_recoEtaPt"]->Fill( globDimuonEta, globDimuonPt );
535  ME["resultDimuon_recoRapPt"]->Fill( globDimuonRap, globDimuonPt );
536  if(highPt) ME["resultDimuon_recoPtDR"]->Fill( globDimuonPt, globDimuonDR );
537  if(highPt) ME["resultDimuon_recoPtDRpos"]->Fill( globDimuonPt, globDimuonDRpos );
538  }
539  }
540  }
541 }
542 
543 void HeavyFlavorValidation::endJob(){
544 }
545 
546 int HeavyFlavorValidation::getMotherId( const Candidate * p ){
547  const Candidate* mother = p->mother();
548  if( mother ){
549  if( mother->pdgId() == p->pdgId() ){
550  return getMotherId(mother);
551  }else{
552  return mother->pdgId();
553  }
554  }else{
555  return 0;
556  }
557 }
558 
559 void HeavyFlavorValidation::match( MonitorElement * me, vector<LeafCandidate> & from, vector<LeafCandidate> & to, double dRMatchingCut, vector<int> & map ){
560  vector<double> dR(from.size());
561  for(size_t i=0; i<from.size(); i++){
562  map[i] = -1;
563  dR[i] = 10.;
564  //find closest
565  for(size_t j=0; j<to.size(); j++){
566  double dRtmp = deltaR<double>(from[i].eta(), from[i].phi(), to[j].eta(), to[j].phi());
567  if( dRtmp < dR[i] ){
568  dR[i] = dRtmp;
569  map[i] = j;
570  }
571  }
572  //fill matching histo
573  if( map[i] != -1 ){
574  me->Fill( to[map[i]].eta()-from[i].eta(), deltaPhi<double>(to[map[i]].phi(), from[i].phi()) );
575  }
576  //apply matching cut
577  if( dR[i] > dRMatchingCut ){
578  map[i] = -1;
579  }
580  //remove duplication
581  if( map[i] != -1 ){
582  for(size_t k=0; k<i; k++){
583  if( map[k] != -1 && map[i] == map[k] ){
584  if( dR[i] < dR[k] ){
585  map[k] = -1;
586  }else{
587  map[i] = -1;
588  }
589  break;
590  }
591  }
592  }
593  }
594 }
595 
596 void HeavyFlavorValidation::myBook2D( TString name, vector<double> &ptBins, TString ptLabel, vector<double> &etaBins, TString etaLabel, TString title )
597 {
598 // dqmStore->setCurrentFolder(dqmFolder+"/"+folder);
599  int ptN = ptBins.size()==3 ? (int)ptBins[0]+1 : ptBins.size();
600  Double_t *pt = new Double_t[ ptN ];
601  for(int i=0; i<ptN; i++){
602  pt[i] = ptBins.size()==3 ? ptBins[1] + i*(ptBins[2]-ptBins[1])/ptBins[0] : ptBins[i] ;
603  }
604  int etaN = etaBins.size()==3 ? (int)etaBins[0]+1 : etaBins.size();
605  Double_t *eta = new Double_t[ etaN ];
606  for(int i=0; i<etaN; i++){
607  eta[i] = etaBins.size()==3 ? etaBins[1] + i*(etaBins[2]-etaBins[1])/etaBins[0] : etaBins[i] ;
608  }
609  TH2F *h = new TH2F( name, name, ptN-1, pt, etaN-1, eta );
610  h->SetXTitle(ptLabel);
611  h->SetYTitle(etaLabel);
612  h->SetTitle(title);
613  ME[name] = dqmStore->book2D( name.Data(), h );
614  delete h;
615 }
616 
617 void HeavyFlavorValidation::myBookProfile2D( TString name, vector<double> &ptBins, TString ptLabel, vector<double> &etaBins, TString etaLabel, TString title )
618 {
619 // dqmStore->setCurrentFolder(dqmFolder+"/"+folder);
620  int ptN = ptBins.size()==3 ? (int)ptBins[0]+1 : ptBins.size();
621  Double_t *pt = new Double_t[ ptN ];
622  for(int i=0; i<ptN; i++){
623  pt[i] = ptBins.size()==3 ? ptBins[1] + i*(ptBins[2]-ptBins[1])/ptBins[0] : ptBins[i] ;
624  }
625  int etaN = etaBins.size()==3 ? (int)etaBins[0]+1 : etaBins.size();
626  Double_t *eta = new Double_t[ etaN ];
627  for(int i=0; i<etaN; i++){
628  eta[i] = etaBins.size()==3 ? etaBins[1] + i*(etaBins[2]-etaBins[1])/etaBins[0] : etaBins[i] ;
629  }
630  TProfile2D *h = new TProfile2D( name, name, ptN-1, pt, etaN-1, eta );
631  h->SetXTitle(ptLabel);
632  h->SetYTitle(etaLabel);
633  h->SetTitle(title);
634  ME[name] = dqmStore->bookProfile2D( name.Data(), h );
635  delete h;
636 }
637 
638 void HeavyFlavorValidation::myBook1D( TString name, vector<double> &bins, TString label, TString title )
639 {
640 // dqmStore->setCurrentFolder(dqmFolder+"/"+folder);
641  int binsN = bins.size()==3 ? (int)bins[0]+1 : bins.size();
642  Double_t *myBins = new Double_t[ binsN ];
643  for(int i=0; i<binsN; i++){
644  myBins[i] = bins.size()==3 ? bins[1] + i*(bins[2]-bins[1])/bins[0] : bins[i] ;
645  }
646  TH1F *h = new TH1F( name, name, binsN-1, myBins );
647  h->SetXTitle(label);
648  h->SetTitle(title);
649  ME[name] = dqmStore->book1D( name.Data(), h );
650  delete h;
651 }
652 
653 HeavyFlavorValidation::~HeavyFlavorValidation(){
654 }
655 
656 //define this as a plug-in
657 DEFINE_FWK_MODULE(HeavyFlavorValidation);
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
jptDQMConfig_cff.ptBins
tuple ptBins
Definition: jptDQMConfig_cff.py:62
edm::Event::triggerNames
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
Definition: Event.cc:204
reco::Candidate::mother
virtual const Candidate * mother(size_type i=0) const =0
return pointer to mother
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:872
HeavyFlavorValidation::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: HeavyFlavorValidation.cc:280
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
jptDQMConfig_cff.deltaPhiBins
tuple deltaPhiBins
Definition: jptDQMConfig_cff.py:68
trigger::TriggerL1Mu
enum start value shifted to 81 so as to avoid clashes with PDG codes
Definition: TriggerTypeDefs.h:31
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
HLTConfigProvider::triggerNames
const std::vector< std::string > & triggerNames() const
names of trigger paths
Definition: HLTConfigProvider.h:100
HeavyFlavorValidation::filterNamesLevels
vector< pair< string, int > > filterNamesLevels
Definition: HeavyFlavorValidation.cc:99
HeavyFlavorValidation::myBook1D
void myBook1D(TString name, vector< double > &xBins, TString label, TString title)
Definition: HeavyFlavorValidation.cc:638
jptDQMConfig_cff.deltaEtaBins
tuple deltaEtaBins
Definition: jptDQMConfig_cff.py:66
muonMass
const double muonMass
Definition: ParticleReplacerZtautau.cc:33
edm::TriggerNames::size
Strings::size_type size() const
Definition: TriggerNames.cc:39
genParticleCandidates2GenParticles_cfi.genParticles
tuple genParticles
Definition: genParticleCandidates2GenParticles_cfi.py:3
eta
T eta() const
Definition: Basic3DVectorLD.h:179
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:105
HeavyFlavorValidation::dimuonEtaBins
vector< double > dimuonEtaBins
Definition: HeavyFlavorValidation.cc:95
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
DeDxDiscriminatorTools::charge
double charge(const std::vector< uint8_t > &Ampls)
Definition: DeDxDiscriminatorTools.cc:43
detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:15
HeavyFlavorValidation::triggerSummaryAODTag
EDGetTokenT< TriggerEvent > triggerSummaryAODTag
Definition: HeavyFlavorValidation.cc:79
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HeavyFlavorValidation::genParticlesTag
EDGetTokenT< GenParticleCollection > genParticlesTag
Definition: HeavyFlavorValidation.cc:82
HeavyFlavorValidation::myBook2D
void myBook2D(TString name, vector< double > &xBins, TString xLabel, vector< double > &yBins, TString yLabel)
Definition: HeavyFlavorValidation.cc:63
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:161
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:250
math::PtEtaPhiMLorentzVector
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
ME
Definition: ME.h:11
HeavyFlavorValidation::endJob
virtual void endJob() override
Definition: HeavyFlavorValidation.cc:543
HeavyFlavorValidation::getMotherId
int getMotherId(const Candidate *p)
Definition: HeavyFlavorValidation.cc:546
iEvent
int iEvent
Definition: GenABIO.cc:243
HeavyFlavorValidation::match
void match(MonitorElement *me, vector< LeafCandidate > &from, vector< LeafCandidate > &to, double deltaRMatchingCut, vector< int > &map)
Definition: HeavyFlavorValidation.cc:559
HeavyFlavorValidation::triggerSummaryRAWTag
EDGetTokenT< TriggerEventWithRefs > triggerSummaryRAWTag
Definition: HeavyFlavorValidation.cc:78
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
j
int j
Definition: DBlmapReader.cc:9
f
double f[11][100]
Definition: MuScleFitUtils.cc:77
DQMStore::setVerbose
void setVerbose(unsigned level)
Definition: DQMStore.cc:548
PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi.dR
tuple dR
Definition: PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi.py:19
triggerResults
static std::string const triggerResults
Definition: EdmProvDump.cc:41
first
bool first
Definition: L1TdeRCT.cc:79
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
HeavyFlavorValidation::HeavyFlavorValidation
HeavyFlavorValidation(const edm::ParameterSet &)
Definition: HeavyFlavorValidation.cc:103
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:81
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:113
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:81
HeavyFlavorValidation::triggerResultsTag
EDGetTokenT< TriggerResults > triggerResultsTag
Definition: HeavyFlavorValidation.cc:80
h
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
Definition: Activities.doc:4
edm::TriggerNames::triggerName
std::string const & triggerName(unsigned int index) const
Definition: TriggerNames.cc:27
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:20
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
Definition: HLTConfigProvider.cc:65
HeavyFlavorValidation::myBookProfile2D
void myBookProfile2D(TString name, vector< double > &xBins, TString xLabel, vector< double > &yBins, TString yLabel, TString title)
Definition: HeavyFlavorValidation.cc:617
HeavyFlavorValidation::myBook1D
void myBook1D(TString name, vector< double > &xBins, TString label)
Definition: HeavyFlavorValidation.cc:71
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
HeavyFlavorValidation::myBook2D
void myBook2D(TString name, vector< double > &xBins, TString xLabel, vector< double > &yBins, TString yLabel, TString title)
Definition: HeavyFlavorValidation.cc:596
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:6
testEve_cfg.level
tuple level
Definition: testEve_cfg.py:34
x
Definition: DDAxes.h:10
jptDQMConfig_cff.etaBins
tuple etaBins
Definition: jptDQMConfig_cff.py:53
HeavyFlavorValidation::ME
map< TString, MonitorElement * > ME
Definition: HeavyFlavorValidation.cc:98
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:1000
HeavyFlavorValidation::beginRun
virtual void beginRun(const Run &iRun, const EventSetup &iSetup) override
Definition: HeavyFlavorValidation.cc:130
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
HeavyFlavorValidation::myBookProfile2D
void myBookProfile2D(TString name, vector< double > &xBins, TString xLabel, vector< double > &yBins, TString yLabel)
Definition: HeavyFlavorValidation.cc:67
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:584
HeavyFlavorValidation::recoMuonsTag
EDGetTokenT< MuonCollection > recoMuonsTag
Definition: HeavyFlavorValidation.cc:81
edm::Run
Definition: Run.h:41
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:1330
phi
Definition: DDAxes.h:10
Generated for CMSSW Reference Manual by doxygen 1.8.5