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

#include <HLTInclusiveVBFSource.h>

Inheritance diagram for HLTInclusiveVBFSource:
DQMEDAnalyzer edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >

Classes

class  PathInfo
 
class  PathInfoCollection
 

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 HLTInclusiveVBFSource (const edm::ParameterSet &)
 
 ~HLTInclusiveVBFSource () override
 
- Public Member Functions inherited from DQMEDAnalyzer
void accumulate (edm::Event const &event, edm::EventSetup const &setup) final
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &setup) final
 
void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer ()
 
void endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void endRun (edm::Run const &run, edm::EventSetup const &setup) final
 
virtual bool getCanSaveByLumi ()
 
- Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProducer & operator= (const EDProducer &)=delete
 

Private Member Functions

virtual bool isBarrel (double eta)
 
virtual bool isEndCap (double eta)
 
virtual bool isForward (double eta)
 
virtual bool isHLTPathAccepted (std::string pathName)
 
virtual bool isTriggerObjectFound (std::string objectName)
 
virtual bool validPathHLT (std::string path)
 

Private Attributes

reco::CaloJetCollection calojet
 
edm::Handle
< reco::CaloJetCollection > 		calojetColl_
 
edm::EDGetTokenT
< reco::CaloJetCollection > 		caloJetsToken
 
edm::Handle
< reco::CaloMETCollection > 		calometColl_
 
edm::EDGetTokenT
< reco::CaloMETCollection > 		caloMetToken
 
bool check_mjj650_Pt35_DEta3p5
 
bool check_mjj650_Pt40_DEta3p5
 
bool check_mjj700_Pt35_DEta3p5
 
bool check_mjj700_Pt40_DEta3p5
 
bool check_mjj750_Pt35_DEta3p5
 
bool check_mjj750_Pt40_DEta3p5
 
bool check_mjj800_Pt35_DEta3p5
 
bool check_mjj800_Pt40_DEta3p5
 
bool checkdR_crossOrder
 
bool checkdR_sameOrder
 
bool checkHLT
 
bool checkHLTIndex
 
bool checkOffline
 
bool debug_
 
double deltaRMatch_
 
std::string dirname_
 
float dR_HLT_RECO_11
 
float dR_HLT_RECO_12
 
float dR_HLT_RECO_21
 
float dR_HLT_RECO_22
 
bool etaOpposite_
 
std::string filtername
 
double hlt_deltaetajet
 
double hlt_deltaphijet
 
double hlt_ejet1
 
double hlt_ejet2
 
double hlt_etajet1
 
double hlt_etajet2
 
double hlt_invmassjet
 
double hlt_phijet1
 
double hlt_phijet2
 
double hlt_ptjet1
 
double hlt_ptjet2
 
double hlt_pxjet1
 
double hlt_pxjet2
 
double hlt_pyjet1
 
double hlt_pyjet2
 
double hlt_pzjet1
 
double hlt_pzjet2
 
HLTConfigProvider hltConfig_
 
PathInfoCollection hltPathsAll_
 
double minDeltaEta_
 
double minInvMass_
 
double minPtHigh_
 
double minPtLow_
 
int nCount_
 
std::vector< std::string > path_
 
std::string pathname
 
reco::PFJetCollection pfjet
 
edm::Handle
< reco::PFJetCollection > 		pfjetColl_
 
edm::EDGetTokenT< edm::View
< reco::PFJet > > 		pfJetsToken
 
edm::Handle
< reco::PFMETCollection > 		pfmetColl_
 
edm::EDGetTokenT< edm::View
< reco::PFMET > > 		pfMetToken
 
std::vector< int > prescUsed_
 
std::string processname_
 
double reco_deltaetajet
 
double reco_deltaphijet
 
double reco_ejet1
 
double reco_ejet2
 
double reco_etajet1
 
double reco_etajet2
 
double reco_invmassjet
 
double reco_phijet1
 
double reco_phijet2
 
double reco_ptjet1
 
double reco_ptjet2
 
double reco_pxjet1
 
double reco_pxjet2
 
double reco_pyjet1
 
double reco_pyjet2
 
double reco_pzjet1
 
double reco_pzjet2
 
edm::TriggerNames triggerNames_
 
edm::Handle
< trigger::TriggerEvent > 		triggerObj_
 
edm::Handle< edm::TriggerResults > triggerResults_
 
edm::EDGetTokenT
< edm::TriggerResults > 		triggerResultsFUToken
 
edm::InputTag triggerResultsLabel_
 
edm::EDGetTokenT
< edm::TriggerResults > 		triggerResultsToken
 
edm::EDGetTokenT
< trigger::TriggerEvent > 		triggerSummaryFUToken
 
edm::InputTag triggerSummaryLabel_
 
edm::EDGetTokenT
< trigger::TriggerEvent > 		triggerSummaryToken
 

Additional Inherited Members

- Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore DQMStore
 
typedef dqm::reco::MonitorElement MonitorElement
 
- Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using CacheTypes = CacheContexts< T...>
 
using GlobalCache = typename CacheTypes::GlobalCache
 
using HasAbility = AbilityChecker< T...>
 
using InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
 
using LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
 
using LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
 
using LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
 
using RunCache = typename CacheTypes::RunCache
 
using RunContext = RunContextT< RunCache, GlobalCache >
 
using RunSummaryCache = typename CacheTypes::RunSummaryCache
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static void globalEndJob (DQMEDAnalyzerGlobalCache const *)
 
static void globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
 
static void globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
 
static std::unique_ptr
< DQMEDAnalyzerGlobalCache > 		initializeGlobalCache (edm::ParameterSet const &)
 
- Protected Member Functions inherited from DQMEDAnalyzer
uint64_t meId () const
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken > lumiToken_
 
edm::EDPutTokenT< DQMToken > runToken_
 
unsigned int streamId_
 

Detailed Description

Definition at line 50 of file HLTInclusiveVBFSource.h.

Constructor & Destructor Documentation

HLTInclusiveVBFSource::HLTInclusiveVBFSource ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 45 of file HLTInclusiveVBFSource.cc.

References debug_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_FULL_cff::InputTag, edm::InputTag::instance(), edm::InputTag::label(), LogDebug, and AlCaHLTBitMon_QueryRunRegistry::string.

45  {
46  LogDebug("HLTInclusiveVBFSource") << "constructor....";
47  nCount_ = 0;
48 
49  dirname_ = iConfig.getUntrackedParameter("dirname", std::string("HLT/InclusiveVBF"));
50  processname_ = iConfig.getParameter<std::string>("processname");
51  triggerSummaryLabel_ = iConfig.getParameter<edm::InputTag>("triggerSummaryLabel");
52  triggerResultsLabel_ = iConfig.getParameter<edm::InputTag>("triggerResultsLabel");
53  triggerSummaryToken = consumes<trigger::TriggerEvent>(triggerSummaryLabel_);
54  triggerResultsToken = consumes<edm::TriggerResults>(triggerResultsLabel_);
55  triggerSummaryFUToken = consumes<trigger::TriggerEvent>(
56  edm::InputTag(triggerSummaryLabel_.label(), triggerSummaryLabel_.instance(), std::string("FU")));
57  triggerResultsFUToken = consumes<edm::TriggerResults>(
58  edm::InputTag(triggerResultsLabel_.label(), triggerResultsLabel_.instance(), std::string("FU")));
59 
60  //path_ = iConfig.getUntrackedParameter<std::vector<std::string> >("paths");
61  //l1path_ = iConfig.getUntrackedParameter<std::vector<std::string> >("l1paths");
62  debug_ = iConfig.getUntrackedParameter<bool>("debug", false);
63 
64  caloJetsToken = consumes<reco::CaloJetCollection>(iConfig.getParameter<edm::InputTag>("CaloJetCollectionLabel"));
65  caloMetToken = consumes<reco::CaloMETCollection>(iConfig.getParameter<edm::InputTag>("CaloMETCollectionLabel"));
66  pfJetsToken = consumes<edm::View<reco::PFJet> >(iConfig.getParameter<edm::InputTag>("PFJetCollectionLabel"));
67  pfMetToken = consumes<edm::View<reco::PFMET> >(iConfig.getParameter<edm::InputTag>("PFMETCollectionLabel"));
68  //jetID = new reco::helper::JetIDHelper(iConfig.getParameter<ParameterSet>("JetIDParams"));
69 
70  minPtHigh_ = iConfig.getUntrackedParameter<double>("minPtHigh", 40.);
71  minPtLow_ = iConfig.getUntrackedParameter<double>("minPtLow", 40.);
72  minDeltaEta_ = iConfig.getUntrackedParameter<double>("minDeltaEta", 3.5);
73  deltaRMatch_ = iConfig.getUntrackedParameter<double>("deltaRMatch", 0.1);
74  minInvMass_ = iConfig.getUntrackedParameter<double>("minInvMass", 1000.0);
75  etaOpposite_ = iConfig.getUntrackedParameter<bool>("etaOpposite", true);
76 
77  check_mjj650_Pt35_DEta3p5 = false;
78  check_mjj700_Pt35_DEta3p5 = false;
79  check_mjj750_Pt35_DEta3p5 = false;
80  check_mjj800_Pt35_DEta3p5 = false;
81  check_mjj650_Pt40_DEta3p5 = false;
82  check_mjj700_Pt40_DEta3p5 = false;
83  check_mjj750_Pt40_DEta3p5 = false;
84  check_mjj800_Pt40_DEta3p5 = false;
85 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HLTInclusiveVBFSource::triggerResultsFUToken
edm::EDGetTokenT< edm::TriggerResults > triggerResultsFUToken
Definition: HLTInclusiveVBFSource.h:93
HLTInclusiveVBFSource::caloMetToken
edm::EDGetTokenT< reco::CaloMETCollection > caloMetToken
Definition: HLTInclusiveVBFSource.h:100
HLTInclusiveVBFSource::triggerSummaryToken
edm::EDGetTokenT< trigger::TriggerEvent > triggerSummaryToken
Definition: HLTInclusiveVBFSource.h:94
HLTInclusiveVBFSource::triggerSummaryLabel_
edm::InputTag triggerSummaryLabel_
Definition: HLTInclusiveVBFSource.h:86
HLTInclusiveVBFSource::triggerResultsToken
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken
Definition: HLTInclusiveVBFSource.h:92
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HLTInclusiveVBFSource::caloJetsToken
edm::EDGetTokenT< reco::CaloJetCollection > caloJetsToken
Definition: HLTInclusiveVBFSource.h:99
HLTInclusiveVBFSource::pfJetsToken
edm::EDGetTokenT< edm::View< reco::PFJet > > pfJetsToken
Definition: HLTInclusiveVBFSource.h:97
HLTInclusiveVBFSource::pfMetToken
edm::EDGetTokenT< edm::View< reco::PFMET > > pfMetToken
Definition: HLTInclusiveVBFSource.h:98
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
HLT_FULL_cff.InputTag
tuple InputTag
Definition: HLT_FULL_cff.py:71427
HLTInclusiveVBFSource::triggerResultsLabel_
edm::InputTag triggerResultsLabel_
Definition: HLTInclusiveVBFSource.h:88
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
HLTInclusiveVBFSource::triggerSummaryFUToken
edm::EDGetTokenT< trigger::TriggerEvent > triggerSummaryFUToken
Definition: HLTInclusiveVBFSource.h:95
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233
HLTInclusiveVBFSource::~HLTInclusiveVBFSource ( )
override

Definition at line 87 of file HLTInclusiveVBFSource.cc.

87  {
88  //
89  // do anything here that needs to be done at desctruction time
90  // (e.g. close files, deallocate resources etc.)
91 }

Member Function Documentation

void HLTInclusiveVBFSource::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 93 of file HLTInclusiveVBFSource.cc.

References funct::abs(), gather_cfg::cout, debug_, reco::deltaPhi(), reco::deltaR(), relval_parameters_module::energy, PVValHelper::eta, edm::Event::getByToken(), fwrapper::jets, HLT_FULL_cff::jetSrc, PatBasicFWLiteJetAnalyzer_Selector_cfg::metSrc, reco::btau::neutralHadronEnergyFraction, timingPdfMaker::pathname, DiDispStaMuonMonitor_cfi::pt, dt_dqm_sourceclient_common_cff::reco, mathSSE::sqrt(), edm::Event::triggerNames(), and findQualityFiles::v.

93  {
94  using namespace std;
95  using namespace edm;
96  using namespace trigger;
97  using namespace reco;
98 
99  if (debug_)
100  cout << "DEBUG-0: Start to analyze" << endl;
101 
102  //****************************************************
103  // Get trigger information.
104  //****************************************************
105  //
106  //---------- triggerResults ----------
107  iEvent.getByToken(triggerResultsToken, triggerResults_);
108  if (!triggerResults_.isValid()) {
109  iEvent.getByToken(triggerResultsFUToken, triggerResults_);
110  if (!triggerResults_.isValid()) {
111  edm::LogInfo("HLTInclusiveVBFSource") << "TriggerResults not found, "
112  "skipping event";
113  return;
114  }
115  }
116 
117  // Check how many HLT triggers are in triggerResults
118  triggerNames_ = iEvent.triggerNames(*triggerResults_);
119 
120  //---------- triggerSummary ----------
121  iEvent.getByToken(triggerSummaryToken, triggerObj_);
122  if (!triggerObj_.isValid()) {
123  iEvent.getByToken(triggerSummaryFUToken, triggerObj_);
124  if (!triggerObj_.isValid()) {
125  edm::LogInfo("HLTInclusiveVBFSource") << "TriggerEvent not found, "
126  "skipping event";
127  return;
128  }
129  }
130 
131  if (debug_)
132  cout << "DEBUG-1: Trigger information" << endl;
133 
134  //****************************************************
135  // Get AOD information
136  //****************************************************
137  //
138  edm::Handle<edm::View<reco::PFMET> > metSrc;
139  bool ValidPFMET_ = iEvent.getByToken(pfMetToken, metSrc);
140  if (!ValidPFMET_)
141  return;
142 
143  edm::Handle<edm::View<reco::PFJet> > jetSrc;
144  bool ValidPFJet_ = iEvent.getByToken(pfJetsToken, jetSrc);
145  if (!ValidPFJet_)
146  return;
147 
148  if (!metSrc.isValid())
149  return;
150  if (!jetSrc.isValid())
151  return;
152  const edm::View<reco::PFMET>& mets = *metSrc;
153  const edm::View<reco::PFJet>& jets = *jetSrc;
154  if (jets.empty())
155  return;
156  if (mets.empty())
157  return;
158 
159  if (debug_)
160  cout << "DEBUG-2: AOD Information" << endl;
161 
162  //****************************************************
163  // Variable setting
164  //****************************************************
165  //
166  pathname = "dummy";
167  filtername = "dummy";
168 
169  //
170  reco_ejet1 = 0.;
171  //reco_etjet1 = 0.;
172  reco_pxjet1 = 0.;
173  reco_pyjet1 = 0.;
174  reco_pzjet1 = 0.;
175  reco_ptjet1 = 0.;
176  reco_etajet1 = 0.;
177  reco_phijet1 = 0.;
178 
179  //
180  reco_ejet2 = 0.;
181  //reco_etjet2 = 0.;
182  reco_pxjet2 = 0.;
183  reco_pyjet2 = 0.;
184  reco_pzjet2 = 0.;
185  reco_ptjet2 = 0.;
186  reco_etajet2 = 0.;
187  reco_phijet2 = 0.;
188 
189  //
190  hlt_ejet1 = 0.;
191  //hlt_etjet1 = 0.;
192  hlt_pxjet1 = 0.;
193  hlt_pyjet1 = 0.;
194  hlt_pzjet1 = 0.;
195  hlt_ptjet1 = 0.;
196  hlt_etajet1 = 0.;
197  hlt_phijet1 = 0.;
198 
199  //
200  hlt_ejet2 = 0.;
201  //hlt_etjet2 = 0.;
202  hlt_pxjet2 = 0.;
203  hlt_pyjet2 = 0.;
204  hlt_pzjet2 = 0.;
205  hlt_ptjet2 = 0.;
206  hlt_etajet2 = 0.;
207  hlt_phijet2 = 0.;
208 
209  //
210  checkOffline = false;
211  checkHLT = false;
212  checkHLTIndex = false;
213 
214  //
215  dR_HLT_RECO_11 = 0.;
216  dR_HLT_RECO_22 = 0.;
217  dR_HLT_RECO_12 = 0.;
218  dR_HLT_RECO_21 = 0.;
219 
220  //
221  checkdR_sameOrder = false;
222  checkdR_crossOrder = false;
223 
224  //
225  reco_deltaetajet = 0.;
226  reco_deltaphijet = 0.;
227  reco_invmassjet = 0.;
228  hlt_deltaetajet = 0.;
229  hlt_deltaphijet = 0.;
230  hlt_invmassjet = 0.;
231 
232  //****************************************************
233  // Offline analysis
234  //****************************************************
235  //
236  checkOffline = false;
237  for (unsigned int ijet1 = 0; ijet1 < jets.size(); ijet1++) {
238  if (jets[ijet1].neutralHadronEnergyFraction() > 0.99)
239  continue;
240  if (jets[ijet1].neutralEmEnergyFraction() > 0.99)
241  continue;
242  for (unsigned int ijet2 = ijet1 + 1; ijet2 < jets.size(); ijet2++) {
243  if (jets[ijet2].neutralHadronEnergyFraction() > 0.99)
244  continue;
245  if (jets[ijet2].neutralEmEnergyFraction() > 0.99)
246  continue;
247  //
248  reco_ejet1 = jets[ijet1].energy();
249  //reco_etjet1 = jets[ijet1].et();
250  reco_pxjet1 = jets[ijet1].momentum().X();
251  reco_pyjet1 = jets[ijet1].momentum().Y();
252  reco_pzjet1 = jets[ijet1].momentum().Z();
253  reco_ptjet1 = jets[ijet1].pt();
254  reco_etajet1 = jets[ijet1].eta();
255  reco_phijet1 = jets[ijet1].phi();
256  //
257  reco_ejet2 = jets[ijet2].energy();
258  //reco_etjet2 = jets[ijet2].et();
259  reco_pxjet2 = jets[ijet2].momentum().X();
260  reco_pyjet2 = jets[ijet2].momentum().Y();
261  reco_pzjet2 = jets[ijet2].momentum().Z();
262  reco_ptjet2 = jets[ijet2].pt();
263  reco_etajet2 = jets[ijet2].eta();
264  reco_phijet2 = jets[ijet2].phi();
265  //
266  reco_deltaetajet = reco_etajet1 - reco_etajet2;
267  reco_deltaphijet = reco::deltaPhi(reco_phijet1, reco_phijet2);
268  reco_invmassjet = sqrt((reco_ejet1 + reco_ejet2) * (reco_ejet1 + reco_ejet2) -
269  (reco_pxjet1 + reco_pxjet2) * (reco_pxjet1 + reco_pxjet2) -
270  (reco_pyjet1 + reco_pyjet2) * (reco_pyjet1 + reco_pyjet2) -
271  (reco_pzjet1 + reco_pzjet2) * (reco_pzjet1 + reco_pzjet2));
272 
273  //
274  if (reco_ptjet1 < minPtHigh_)
275  continue;
276  if (reco_ptjet2 < minPtLow_)
277  continue;
278  if (etaOpposite_ == true && reco_etajet1 * reco_etajet2 > 0)
279  continue;
280  if (std::abs(reco_deltaetajet) < minDeltaEta_)
281  continue;
282  if (std::abs(reco_invmassjet) < minInvMass_)
283  continue;
284 
285  //
286  if (debug_)
287  cout << "DEBUG-3" << endl;
288  checkOffline = true;
289  break;
290  }
291  if (checkOffline == true)
292  break;
293  }
294  if (checkOffline == false)
295  return;
296 
297  //****************************************************
298  // Trigger efficiency: Loop for all VBF paths
299  //****************************************************
300  //const unsigned int numberOfPaths(hltConfig_.size());
301  const trigger::TriggerObjectCollection& toc(triggerObj_->getObjects());
302  for (auto& v : hltPathsAll_) {
303  checkHLT = false;
304  checkHLTIndex = false;
305 
306  //
307  v.getMEhisto_RECO_deltaEta_DiJet()->Fill(reco_deltaetajet);
308  v.getMEhisto_RECO_deltaPhi_DiJet()->Fill(reco_deltaphijet);
309  v.getMEhisto_RECO_invMass_DiJet()->Fill(reco_invmassjet);
310 
311  //
312  if (debug_)
313  cout << "DEBUG-4-0: Path loops" << endl;
314 
315  //
316  if (isHLTPathAccepted(v.getPath()) == false)
317  continue;
318  checkHLT = true;
319 
320  //
321  if (debug_)
322  cout << "DEBUG-4-1: Path is accepted. Now we are looking for " << v.getLabel() << " module." << endl;
323 
324  //
325  edm::InputTag hltTag(v.getLabel(), "", processname_);
326  const int hltIndex = triggerObj_->filterIndex(hltTag);
327  if (hltIndex >= triggerObj_->sizeFilters())
328  continue;
329  checkHLT = true;
330  if (debug_)
331  cout << "DEBUG-4-2: HLT module " << v.getLabel() << " exists" << endl;
332  const trigger::Keys& khlt = triggerObj_->filterKeys(hltIndex);
333  auto kj = khlt.begin();
334  for (; kj != khlt.end(); kj += 2) {
335  if (debug_)
336  cout << "DEBUG-5" << endl;
337  checkdR_sameOrder = false;
338  checkdR_crossOrder = false; //
339  hlt_ejet1 = toc[*kj].energy();
340  //hlt_etjet1 = toc[*kj].et();
341  hlt_pxjet1 = toc[*kj].px();
342  hlt_pyjet1 = toc[*kj].py();
343  hlt_pzjet1 = toc[*kj].pz();
344  hlt_ptjet1 = toc[*kj].pt();
345  hlt_etajet1 = toc[*kj].eta();
346  hlt_phijet1 = toc[*kj].phi();
347  //
348  hlt_ejet2 = toc[*(kj + 1)].energy();
349  //hlt_etjet2 = toc[*(kj+1)].et();
350  hlt_pxjet2 = toc[*(kj + 1)].px();
351  hlt_pyjet2 = toc[*(kj + 1)].py();
352  hlt_pzjet2 = toc[*(kj + 1)].pz();
353  hlt_ptjet2 = toc[*(kj + 1)].pt();
354  hlt_etajet2 = toc[*(kj + 1)].eta();
355  hlt_phijet2 = toc[*(kj + 1)].phi();
356  //
357  dR_HLT_RECO_11 = reco::deltaR(hlt_etajet1, hlt_phijet1, reco_etajet1, reco_phijet1);
358  dR_HLT_RECO_22 = reco::deltaR(hlt_etajet2, hlt_phijet2, reco_etajet2, reco_phijet2);
359  dR_HLT_RECO_12 = reco::deltaR(hlt_etajet1, hlt_phijet1, reco_etajet2, reco_phijet2);
360  dR_HLT_RECO_21 = reco::deltaR(hlt_etajet2, hlt_phijet2, reco_etajet1, reco_phijet1);
361  if (dR_HLT_RECO_11 < deltaRMatch_ && dR_HLT_RECO_22 < deltaRMatch_)
362  checkdR_sameOrder = true;
363  if (dR_HLT_RECO_12 < deltaRMatch_ && dR_HLT_RECO_21 < deltaRMatch_)
364  checkdR_crossOrder = true;
365  if (checkdR_sameOrder == false && checkdR_crossOrder == false)
366  continue;
367  checkHLTIndex = true;
368  //
369  if (debug_)
370  cout << "DEBUG-6: Match" << endl;
371  hlt_deltaetajet = hlt_etajet1 - hlt_etajet2;
372  hlt_deltaphijet = reco::deltaPhi(hlt_phijet1, hlt_phijet2);
373  if (checkdR_crossOrder) {
374  hlt_deltaetajet = (-1) * hlt_deltaetajet;
375  hlt_deltaphijet = reco::deltaPhi(hlt_phijet2, hlt_phijet1);
376  }
377  hlt_invmassjet = sqrt((hlt_ejet1 + hlt_ejet2) * (hlt_ejet1 + hlt_ejet2) -
378  (hlt_pxjet1 + hlt_pxjet2) * (hlt_pxjet1 + hlt_pxjet2) -
379  (hlt_pyjet1 + hlt_pyjet2) * (hlt_pyjet1 + hlt_pyjet2) -
380  (hlt_pzjet1 + hlt_pzjet2) * (hlt_pzjet1 + hlt_pzjet2));
381  v.getMEhisto_HLT_deltaEta_DiJet()->Fill(hlt_deltaetajet);
382  v.getMEhisto_HLT_deltaPhi_DiJet()->Fill(hlt_deltaphijet);
383  v.getMEhisto_HLT_invMass_DiJet()->Fill(hlt_invmassjet);
384  //
385  v.getMEhisto_RECO_deltaEta_DiJet_Match()->Fill(reco_deltaetajet);
386  v.getMEhisto_RECO_deltaPhi_DiJet_Match()->Fill(reco_deltaphijet);
387  v.getMEhisto_RECO_invMass_DiJet_Match()->Fill(reco_invmassjet);
388  //
389  v.getMEhisto_RECOHLT_deltaEta()->Fill(reco_deltaetajet, hlt_deltaetajet);
390  v.getMEhisto_RECOHLT_deltaPhi()->Fill(reco_deltaphijet, hlt_deltaphijet);
391  v.getMEhisto_RECOHLT_invMass()->Fill(reco_invmassjet, hlt_invmassjet);
392  //
393  if (checkHLTIndex == true)
394  break;
395  }
396 
397  //****************************************************
398  // Match information
399  //****************************************************
400  if (checkHLT == true && checkHLTIndex == true) {
401  if (debug_)
402  cout << "DEBUG-7: Match" << endl;
403  v.getMEhisto_NumberOfMatches()->Fill(1);
404  } else {
405  if (debug_)
406  cout << "DEBUG-8: Not match" << endl;
407  v.getMEhisto_NumberOfMatches()->Fill(0);
408  }
409  }
410 
411  //****************************************************
412  //
413  //****************************************************
414  for (auto& v : hltPathsAll_) {
415  if (isHLTPathAccepted(v.getPath()) == false)
416  continue;
417  if (debug_)
418  cout << "DEBUG-9: Loop for rate approximation: " << v.getPath() << endl;
419  check_mjj650_Pt35_DEta3p5 = false;
420  check_mjj700_Pt35_DEta3p5 = false;
421  check_mjj750_Pt35_DEta3p5 = false;
422  check_mjj800_Pt35_DEta3p5 = false;
423  check_mjj650_Pt40_DEta3p5 = false;
424  check_mjj700_Pt40_DEta3p5 = false;
425  check_mjj750_Pt40_DEta3p5 = false;
426  check_mjj800_Pt40_DEta3p5 = false;
427  edm::InputTag hltTag(v.getLabel(), "", processname_);
428  const int hltIndex = triggerObj_->filterIndex(hltTag);
429  if (hltIndex >= triggerObj_->sizeFilters())
430  continue;
431  const trigger::Keys& khlt = triggerObj_->filterKeys(hltIndex);
432  auto kj = khlt.begin();
433  for (; kj != khlt.end(); kj += 2) {
434  checkdR_sameOrder = false;
435  checkdR_crossOrder = false;
436  //
437  hlt_ejet1 = toc[*kj].energy();
438  //hlt_etjet1 = toc[*kj].et();
439  hlt_pxjet1 = toc[*kj].px();
440  hlt_pyjet1 = toc[*kj].py();
441  hlt_pzjet1 = toc[*kj].pz();
442  hlt_ptjet1 = toc[*kj].pt();
443  hlt_etajet1 = toc[*kj].eta();
444  hlt_phijet1 = toc[*kj].phi();
445  //
446  hlt_ejet2 = toc[*(kj + 1)].energy();
447  //hlt_etjet2 = toc[*(kj+1)].et();
448  hlt_pxjet2 = toc[*(kj + 1)].px();
449  hlt_pyjet2 = toc[*(kj + 1)].py();
450  hlt_pzjet2 = toc[*(kj + 1)].pz();
451  hlt_ptjet2 = toc[*(kj + 1)].pt();
452  hlt_etajet2 = toc[*(kj + 1)].eta();
453  hlt_phijet2 = toc[*(kj + 1)].phi();
454  //
455  hlt_deltaetajet = hlt_etajet1 - hlt_etajet2;
456  hlt_deltaphijet = reco::deltaPhi(hlt_phijet1, hlt_phijet2);
457  hlt_invmassjet = sqrt((hlt_ejet1 + hlt_ejet2) * (hlt_ejet1 + hlt_ejet2) -
458  (hlt_pxjet1 + hlt_pxjet2) * (hlt_pxjet1 + hlt_pxjet2) -
459  (hlt_pyjet1 + hlt_pyjet2) * (hlt_pyjet1 + hlt_pyjet2) -
460  (hlt_pzjet1 + hlt_pzjet2) * (hlt_pzjet1 + hlt_pzjet2));
461  //
462  if (check_mjj650_Pt35_DEta3p5 == false && hlt_ptjet1 > 35. && hlt_ptjet2 >= 35. && hlt_invmassjet > 650 &&
463  std::abs(hlt_deltaetajet) > 3.5) {
464  check_mjj650_Pt35_DEta3p5 = true;
465  }
466  if (check_mjj700_Pt35_DEta3p5 == false && hlt_ptjet1 > 35. && hlt_ptjet2 >= 35. && hlt_invmassjet > 700 &&
467  std::abs(hlt_deltaetajet) > 3.5) {
468  check_mjj700_Pt35_DEta3p5 = true;
469  }
470  if (check_mjj750_Pt35_DEta3p5 == false && hlt_ptjet1 > 35. && hlt_ptjet2 >= 35. && hlt_invmassjet > 750 &&
471  std::abs(hlt_deltaetajet) > 3.5) {
472  check_mjj750_Pt35_DEta3p5 = true;
473  }
474  if (check_mjj800_Pt35_DEta3p5 == false && hlt_ptjet1 > 35. && hlt_ptjet2 >= 35. && hlt_invmassjet > 800 &&
475  std::abs(hlt_deltaetajet) > 3.5) {
476  check_mjj800_Pt35_DEta3p5 = true;
477  }
478  if (check_mjj650_Pt40_DEta3p5 == false && hlt_ptjet1 > 40. && hlt_ptjet2 >= 40. && hlt_invmassjet > 650 &&
479  std::abs(hlt_deltaetajet) > 3.5) {
480  check_mjj650_Pt40_DEta3p5 = true;
481  }
482  if (check_mjj700_Pt40_DEta3p5 == false && hlt_ptjet1 > 40. && hlt_ptjet2 >= 40. && hlt_invmassjet > 700 &&
483  std::abs(hlt_deltaetajet) > 3.5) {
484  check_mjj700_Pt40_DEta3p5 = true;
485  }
486  if (check_mjj750_Pt40_DEta3p5 == false && hlt_ptjet1 > 40. && hlt_ptjet2 >= 40. && hlt_invmassjet > 750 &&
487  std::abs(hlt_deltaetajet) > 3.5) {
488  check_mjj750_Pt40_DEta3p5 = true;
489  }
490  if (check_mjj800_Pt40_DEta3p5 == false && hlt_ptjet1 > 40. && hlt_ptjet2 >= 40. && hlt_invmassjet > 800 &&
491  std::abs(hlt_deltaetajet) > 3.5) {
492  check_mjj800_Pt40_DEta3p5 = true;
493  }
494  }
495  if (check_mjj650_Pt35_DEta3p5 == true)
496  v.getMEhisto_NumberOfEvents()->Fill(0);
497  if (check_mjj700_Pt35_DEta3p5 == true)
498  v.getMEhisto_NumberOfEvents()->Fill(1);
499  if (check_mjj750_Pt35_DEta3p5 == true)
500  v.getMEhisto_NumberOfEvents()->Fill(2);
501  if (check_mjj800_Pt35_DEta3p5 == true)
502  v.getMEhisto_NumberOfEvents()->Fill(3);
503  if (check_mjj650_Pt40_DEta3p5 == true)
504  v.getMEhisto_NumberOfEvents()->Fill(4);
505  if (check_mjj700_Pt40_DEta3p5 == true)
506  v.getMEhisto_NumberOfEvents()->Fill(5);
507  if (check_mjj750_Pt40_DEta3p5 == true)
508  v.getMEhisto_NumberOfEvents()->Fill(6);
509  if (check_mjj800_Pt40_DEta3p5 == true)
510  v.getMEhisto_NumberOfEvents()->Fill(7);
511  }
512 }
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
HLTInclusiveVBFSource::triggerResultsFUToken
edm::EDGetTokenT< edm::TriggerResults > triggerResultsFUToken
Definition: HLTInclusiveVBFSource.h:93
DiDispStaMuonMonitor_cfi.pt
tuple pt
Definition: DiDispStaMuonMonitor_cfi.py:39
HLTInclusiveVBFSource::hltPathsAll_
PathInfoCollection hltPathsAll_
Definition: HLTInclusiveVBFSource.h:295
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
HLTInclusiveVBFSource::triggerSummaryToken
edm::EDGetTokenT< trigger::TriggerEvent > triggerSummaryToken
Definition: HLTInclusiveVBFSource.h:94
HLT_FULL_cff.jetSrc
tuple jetSrc
Definition: HLT_FULL_cff.py:24898
HLTInclusiveVBFSource::triggerObj_
edm::Handle< trigger::TriggerEvent > triggerObj_
Definition: HLTInclusiveVBFSource.h:87
HLTInclusiveVBFSource::triggerResultsToken
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken
Definition: HLTInclusiveVBFSource.h:92
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:111
edm::Event::triggerNames
edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const override
Definition: Event.cc:275
HLTInclusiveVBFSource::pfJetsToken
edm::EDGetTokenT< edm::View< reco::PFJet > > pfJetsToken
Definition: HLTInclusiveVBFSource.h:97
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
fwrapper::jets
vector< PseudoJet > jets
Definition: fastjetfortran_madfks.cc:44
HLTInclusiveVBFSource::pfMetToken
edm::EDGetTokenT< edm::View< reco::PFMET > > pfMetToken
Definition: HLTInclusiveVBFSource.h:98
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
HLTInclusiveVBFSource::triggerNames_
edm::TriggerNames triggerNames_
Definition: HLTInclusiveVBFSource.h:90
reco::deltaR
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:75
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:19
PatBasicFWLiteJetAnalyzer_Selector_cfg.metSrc
tuple metSrc
Definition: PatBasicFWLiteJetAnalyzer_Selector_cfg.py:22
HLTInclusiveVBFSource::isHLTPathAccepted
virtual bool isHLTPathAccepted(std::string pathName)
Definition: HLTInclusiveVBFSource.cc:765
relval_parameters_module.energy
string energy
Definition: relval_parameters_module.py:29
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:144
HLTInclusiveVBFSource::triggerSummaryFUToken
edm::EDGetTokenT< trigger::TriggerEvent > triggerSummaryFUToken
Definition: HLTInclusiveVBFSource.h:95
HLTInclusiveVBFSource::triggerResults_
edm::Handle< edm::TriggerResults > triggerResults_
Definition: HLTInclusiveVBFSource.h:89
void HLTInclusiveVBFSource::bookHistograms ( DQMStore::IBooker &  iBooker,
edm::Run const &  run,
edm::EventSetup const &  c 
)
overridevirtual

Implements DQMEDAnalyzer.

Definition at line 515 of file HLTInclusiveVBFSource.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookFloat(), gather_cfg::cout, debug_, TrackerOfflineValidation_Dqm_cff::dirName, dqm::impl::MonitorElement::getTH1F(), dqm::impl::MonitorElement::getTH2F(), mps_fire::i, LogDebug, timingPdfMaker::pathname, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, trigger::TriggerJet, HLT_FULL_cff::triggerType, and findQualityFiles::v.

515  {
516  iBooker.setCurrentFolder(dirname_);
517 
518  //--- htlConfig_
519  bool changed(true);
520  if (!hltConfig_.init(run, c, processname_, changed)) {
521  LogDebug("HLTInclusiveVBFSource") << "HLTConfigProvider failed to initialize.";
522  }
523 
524  const unsigned int numberOfPaths(hltConfig_.size());
525  for (unsigned int i = 0; i != numberOfPaths; ++i) {
526  bool numFound = false;
527  pathname = hltConfig_.triggerName(i);
528  filtername = "dummy";
529  unsigned int usedPrescale = 1;
530  unsigned int objectType = 0;
531  std::string triggerType = "";
532 
533  if (pathname.find("HLT_Di") == std::string::npos)
534  continue;
535  if (pathname.find("Jet") == std::string::npos)
536  continue;
537  if (pathname.find("MJJ") == std::string::npos)
538  continue;
539  if (pathname.find("VBF_v") == std::string::npos)
540  continue;
541 
542  if (debug_) {
543  cout << " - Startup:Path = " << pathname << endl;
544  //cout<<" - Startup:PS = "<<hltConfig_.prescaleSize()<<endl;
545  }
546 
547  triggerType = "DiJet_Trigger";
548  objectType = trigger::TriggerJet;
549 
550  // Checking if the trigger exist in HLT table or not
551  for (unsigned int i = 0; i != numberOfPaths; ++i) {
552  std::string HLTname = hltConfig_.triggerName(i);
553  if (HLTname == pathname)
554  numFound = true;
555  }
556 
557  if (numFound == false)
558  continue;
559  std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
560  auto numpathmodule = numpathmodules.begin();
561  for (; numpathmodule != numpathmodules.end(); ++numpathmodule) {
562  edm::InputTag testTag(*numpathmodule, "", processname_);
563  if (hltConfig_.moduleType(*numpathmodule) == "HLTCaloJetVBFFilter" ||
564  hltConfig_.moduleType(*numpathmodule) == "HLTPFJetVBFFilter") {
565  filtername = *numpathmodule;
566  if (debug_)
567  cout << " - Startup:Module = " << hltConfig_.moduleType(*numpathmodule) << ", FilterName = " << filtername
568  << endl;
569  }
570  }
571  if (debug_)
572  cout << " - Startup:Final filter = " << filtername << endl;
573 
574  if (objectType == 0 || numFound == false)
575  continue;
576  //if(debug_){
577  //cout<<"Pathname = "<<pathname
578  // <<", Filtername = "<<filtername
579  // <<", ObjectType = "<<objectType<<endl;
580  //}
581  hltPathsAll_.push_back(PathInfo(usedPrescale, pathname, filtername, processname_, objectType, triggerType));
582  } //Loop over paths
583 
584  //if(debug_) cout<<"== end hltPathsEff_.push_back ======" << endl;
585 
586  std::string dirName = dirname_ + "/MonitorInclusiveVBFTrigger/";
587  for (auto& v : hltPathsAll_) {
588  if (debug_)
589  cout << "Storing: " << v.getPath() << ", Prescale = " << v.getprescaleUsed() << endl;
590  //if(v->getprescaleUsed()!=1) continue;
591 
592  std::string subdirName = dirName + v.getPath();
593  std::string trigPath = "(" + v.getPath() + ")";
594  iBooker.setCurrentFolder(subdirName);
595 
596  MonitorElement* RECO_deltaEta_DiJet;
597  MonitorElement* RECO_deltaPhi_DiJet;
598  MonitorElement* RECO_invMass_DiJet;
599  MonitorElement* HLT_deltaEta_DiJet;
600  MonitorElement* HLT_deltaPhi_DiJet;
601  MonitorElement* HLT_invMass_DiJet;
602  MonitorElement* RECO_deltaEta_DiJet_Match;
603  MonitorElement* RECO_deltaPhi_DiJet_Match;
604  MonitorElement* RECO_invMass_DiJet_Match;
605  MonitorElement* RECOHLT_deltaEta;
606  MonitorElement* RECOHLT_deltaPhi;
607  MonitorElement* RECOHLT_invMass;
608  MonitorElement* NumberOfMatches;
609  MonitorElement* NumberOfEvents;
610 
611  //dummy = iBooker.bookFloat("dummy");
612  RECO_deltaEta_DiJet = iBooker.bookFloat("RECO_deltaEta_DiJet");
613  RECO_deltaPhi_DiJet = iBooker.bookFloat("RECO_deltaPhi_DiJet");
614  RECO_invMass_DiJet = iBooker.bookFloat("RECO_invMass_DiJet");
615  HLT_deltaEta_DiJet = iBooker.bookFloat("HLT_deltaEta_DiJet");
616  HLT_deltaPhi_DiJet = iBooker.bookFloat("HLT_deltaPhi_DiJet ");
617  HLT_invMass_DiJet = iBooker.bookFloat("HLT_invMass_DiJet");
618  RECO_deltaEta_DiJet_Match = iBooker.bookFloat("RECO_deltaEta_DiJet_Match");
619  RECO_deltaPhi_DiJet_Match = iBooker.bookFloat("RECO_deltaPhi_DiJet_Match");
620  RECO_invMass_DiJet_Match = iBooker.bookFloat("RECO_invMass_DiJet_Match");
621  RECOHLT_deltaEta = iBooker.bookFloat("RECOHLT_deltaEta");
622  RECOHLT_deltaPhi = iBooker.bookFloat("RECOHLT_deltaPhi ");
623  RECOHLT_invMass = iBooker.bookFloat("RECOHLT_invMass");
624  NumberOfMatches = iBooker.bookFloat("NumberOfMatches");
625  NumberOfEvents = iBooker.bookFloat("NumberOfEvents");
626 
627  std::string labelname("ME");
628  std::string histoname(labelname + "");
629  std::string title(labelname + "");
630 
631  //RECO_deltaEta_DiJet
632  histoname = labelname + "_RECO_deltaEta_DiJet";
633  title = labelname + "_RECO_deltaEta_DiJet " + trigPath;
634  RECO_deltaEta_DiJet = iBooker.book1D(histoname.c_str(), title.c_str(), 50, -10., 10.);
635  RECO_deltaEta_DiJet->getTH1F();
636 
637  //RECO_deltaPhi_DiJet
638  histoname = labelname + "_RECO_deltaPhi_DiJet";
639  title = labelname + "_RECO_deltaPhi_DiJet " + trigPath;
640  RECO_deltaPhi_DiJet = iBooker.book1D(histoname.c_str(), title.c_str(), 35, -3.5, 3.5);
641  RECO_deltaPhi_DiJet->getTH1F();
642 
643  //RECO_invMass_DiJet
644  histoname = labelname + "_RECO_invMass_DiJet";
645  title = labelname + "_RECO_invMass_DiJet " + trigPath;
646  RECO_invMass_DiJet = iBooker.book1D(histoname.c_str(), title.c_str(), 100, 500., 2000.);
647  RECO_invMass_DiJet->getTH1F();
648 
649  //HLT_deltaEta_DiJet
650  histoname = labelname + "_HLT_deltaEta_DiJet";
651  title = labelname + "_HLT_deltaEta_DiJet " + trigPath;
652  HLT_deltaEta_DiJet = iBooker.book1D(histoname.c_str(), title.c_str(), 50, -10., 10.);
653  HLT_deltaEta_DiJet->getTH1F();
654 
655  //HLT_deltaPhi_DiJet
656  histoname = labelname + "_HLT_deltaPhi_DiJet";
657  title = labelname + "_HLT_deltaPhi_DiJet " + trigPath;
658  HLT_deltaPhi_DiJet = iBooker.book1D(histoname.c_str(), title.c_str(), 35, -3.5, 3.5);
659  HLT_deltaPhi_DiJet->getTH1F();
660 
661  //HLT_invMass_DiJet
662  histoname = labelname + "_HLT_invMass_DiJet";
663  title = labelname + "_HLT_invMass_DiJet " + trigPath;
664  HLT_invMass_DiJet = iBooker.book1D(histoname.c_str(), title.c_str(), 100, 500., 2000.);
665  HLT_invMass_DiJet->getTH1F();
666 
667  //RECO_deltaEta_DiJet_Match
668  histoname = labelname + "_RECO_deltaEta_DiJet_Match";
669  title = labelname + "_RECO_deltaEta_DiJet_Match " + trigPath;
670  RECO_deltaEta_DiJet_Match = iBooker.book1D(histoname.c_str(), title.c_str(), 50, -10., 10.);
671  RECO_deltaEta_DiJet_Match->getTH1F();
672 
673  //RECO_deltaPhi_DiJet_Match
674  histoname = labelname + "_RECO_deltaPhi_DiJet_Match";
675  title = labelname + "_RECO_deltaPhi_DiJet_Match " + trigPath;
676  RECO_deltaPhi_DiJet_Match = iBooker.book1D(histoname.c_str(), title.c_str(), 35, -3.5, 3.5);
677  RECO_deltaPhi_DiJet_Match->getTH1F();
678 
679  //RECO_invMass_DiJet_Match
680  histoname = labelname + "_RECO_invMass_DiJet_Match";
681  title = labelname + "_RECO_invMass_DiJet_Match " + trigPath;
682  RECO_invMass_DiJet_Match = iBooker.book1D(histoname.c_str(), title.c_str(), 100, 500., 2000.);
683  RECO_invMass_DiJet_Match->getTH1F();
684 
685  //RECOHLT_deltaEta
686  histoname = labelname + "_RECOHLT_deltaEta";
687  title = labelname + "_RECOHLT_deltaEta " + trigPath;
688  RECOHLT_deltaEta = iBooker.book2D(histoname.c_str(), title.c_str(), 50, -10., 10., 50, -10., 10.);
689  RECOHLT_deltaEta->getTH2F();
690 
691  //RECOHLT_deltaPhi
692  histoname = labelname + "_RECOHLT_deltaPhi";
693  title = labelname + "_RECOHLT_deltaPhi " + trigPath;
694  RECOHLT_deltaPhi = iBooker.book2D(histoname.c_str(), title.c_str(), 35, -3.5, 3.5, 35, -3.5, 3.5);
695  RECOHLT_deltaPhi->getTH2F();
696 
697  //RECOHLT_invMass
698  histoname = labelname + "_RECOHLT_invMass";
699  title = labelname + "_RECOHLT_invMass " + trigPath;
700  RECOHLT_invMass = iBooker.book2D(histoname.c_str(), title.c_str(), 100, 500., 2000., 100, 500., 2000.);
701  RECOHLT_invMass->getTH2F();
702 
703  //NumberOfMatches
704  histoname = labelname + "_NumberOfMatches ";
705  title = labelname + "_NumberOfMatches " + trigPath;
706  NumberOfMatches = iBooker.book1D(histoname.c_str(), title.c_str(), 2, 0., 2.);
707  NumberOfMatches->getTH1F();
708 
709  //NumberOfEvents
710  histoname = labelname + "_NumberOfEvents";
711  title = labelname + "_NumberOfEvents " + trigPath;
712  NumberOfEvents = iBooker.book1D(histoname.c_str(), title.c_str(), 10, 0., 10.);
713  NumberOfEvents->getTH1F();
714 
715  //}
716  v.setHistos(RECO_deltaEta_DiJet,
717  RECO_deltaPhi_DiJet,
718  RECO_invMass_DiJet,
719  HLT_deltaEta_DiJet,
720  HLT_deltaPhi_DiJet,
721  HLT_invMass_DiJet,
722  RECO_deltaEta_DiJet_Match,
723  RECO_deltaPhi_DiJet_Match,
724  RECO_invMass_DiJet_Match,
725  RECOHLT_deltaEta,
726  RECOHLT_deltaPhi,
727  RECOHLT_invMass,
728  NumberOfMatches,
729  NumberOfEvents);
730  //break;//We need only the first unprescale paths
731  }
732 }
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:63
dqm::legacy::MonitorElement::getTH2F
virtual TH2F * getTH2F() const
Definition: MonitorElement.h:491
dqm::implementation::IBooker::bookFloat
MonitorElement * bookFloat(TString const &name, FUNC onbooking=NOOP())
Definition: DQMStore.h:80
c
const edm::EventSetup & c
Definition: SiStripLAProfileBooker.cc:66
mps_fire.i
i
Definition: mps_fire.py:428
HLTConfigProvider::moduleType
const std::string moduleType(const std::string &module) const
C++ class name of module.
Definition: HLTConfigProvider.h:104
HLTConfigProvider::triggerName
const std::string & triggerName(unsigned int triggerIndex) const
Definition: HLTConfigProvider.h:70
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
HLTInclusiveVBFSource::hltPathsAll_
PathInfoCollection hltPathsAll_
Definition: HLTInclusiveVBFSource.h:295
dqm::legacy::MonitorElement::getTH1F
virtual TH1F * getTH1F() const
Definition: MonitorElement.h:479
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
TrackerOfflineValidation_Dqm_cff.dirName
tuple dirName
Definition: TrackerOfflineValidation_Dqm_cff.py:55
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:76
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
Definition: HLTConfigProvider.cc:36
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:177
HLTInclusiveVBFSource::hltConfig_
HLTConfigProvider hltConfig_
Definition: HLTInclusiveVBFSource.h:109
HLT_FULL_cff.triggerType
tuple triggerType
Definition: HLT_FULL_cff.py:7695
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:144
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233
bool HLTInclusiveVBFSource::isBarrel ( double  eta)
privatevirtual

Definition at line 734 of file HLTInclusiveVBFSource.cc.

References convertSQLitetoXML_cfg::output.

734  {
735  bool output = false;
736  if (fabs(eta) <= 1.3)
737  output = true;
738  return output;
739 }
convertSQLitetoXML_cfg.output
tuple output
Definition: convertSQLitetoXML_cfg.py:72
bool HLTInclusiveVBFSource::isEndCap ( double  eta)
privatevirtual

Definition at line 741 of file HLTInclusiveVBFSource.cc.

References convertSQLitetoXML_cfg::output.

741  {
742  bool output = false;
743  if (fabs(eta) <= 3.0 && fabs(eta) > 1.3)
744  output = true;
745  return output;
746 }
convertSQLitetoXML_cfg.output
tuple output
Definition: convertSQLitetoXML_cfg.py:72
bool HLTInclusiveVBFSource::isForward ( double  eta)
privatevirtual

Definition at line 748 of file HLTInclusiveVBFSource.cc.

References convertSQLitetoXML_cfg::output.

748  {
749  bool output = false;
750  if (fabs(eta) > 3.0)
751  output = true;
752  return output;
753 }
convertSQLitetoXML_cfg.output
tuple output
Definition: convertSQLitetoXML_cfg.py:72
bool HLTInclusiveVBFSource::isHLTPathAccepted ( std::string  pathName)
privatevirtual

Definition at line 765 of file HLTInclusiveVBFSource.cc.

References convertSQLitetoXML_cfg::output.

765  {
766  // triggerResults_, triggerNames_ has to be defined first before calling this method
767  bool output = false;
768  if (triggerResults_.isValid()) {
769  unsigned index = triggerNames_.triggerIndex(pathName);
770  if (index < triggerNames_.size() && triggerResults_->accept(index))
771  output = true;
772  }
773  return output;
774 }
edm::TriggerNames::size
std::size_t size() const
Definition: TriggerNames.cc:59
edm::TriggerNames::triggerIndex
unsigned int triggerIndex(std::string_view name) const
Definition: TriggerNames.cc:52
EgammaValidation_cff.pathName
tuple pathName
Definition: EgammaValidation_cff.py:269
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
HLTInclusiveVBFSource::triggerNames_
edm::TriggerNames triggerNames_
Definition: HLTInclusiveVBFSource.h:90
convertSQLitetoXML_cfg.output
tuple output
Definition: convertSQLitetoXML_cfg.py:72
HLTInclusiveVBFSource::triggerResults_
edm::Handle< edm::TriggerResults > triggerResults_
Definition: HLTInclusiveVBFSource.h:89
bool HLTInclusiveVBFSource::isTriggerObjectFound ( std::string  objectName)
privatevirtual

Definition at line 776 of file HLTInclusiveVBFSource.cc.

References isotrackApplyRegressor::k, and convertSQLitetoXML_cfg::output.

776  {
777  // processname_, triggerObj_ has to be defined before calling this method
778  bool output = false;
779  edm::InputTag testTag(objectName, "", processname_);
780  const int index = triggerObj_->filterIndex(testTag);
781  if (index >= triggerObj_->sizeFilters()) {
782  edm::LogInfo("HLTInclusiveVBFSource") << "no index " << index << " of that name ";
783  } else {
784  const trigger::Keys& k = triggerObj_->filterKeys(index);
785  if (!k.empty())
786  output = true;
787  }
788  return output;
789 }
HLTInclusiveVBFSource::triggerObj_
edm::Handle< trigger::TriggerEvent > triggerObj_
Definition: HLTInclusiveVBFSource.h:87
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
convertSQLitetoXML_cfg.output
tuple output
Definition: convertSQLitetoXML_cfg.py:72
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:19
bool HLTInclusiveVBFSource::validPathHLT ( std::string  path)
privatevirtual

Definition at line 755 of file HLTInclusiveVBFSource.cc.

References dqmiolumiharvest::j, convertSQLitetoXML_cfg::output, and timingPdfMaker::pathname.

755  {
756  // hltConfig_ has to be defined first before calling this method
757  bool output = false;
758  for (unsigned int j = 0; j != hltConfig_.size(); ++j) {
759  if (hltConfig_.triggerName(j) == pathname)
760  output = true;
761  }
762  return output;
763 }
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:63
HLTConfigProvider::triggerName
const std::string & triggerName(unsigned int triggerIndex) const
Definition: HLTConfigProvider.h:70
convertSQLitetoXML_cfg.output
tuple output
Definition: convertSQLitetoXML_cfg.py:72
HLTInclusiveVBFSource::hltConfig_
HLTConfigProvider hltConfig_
Definition: HLTInclusiveVBFSource.h:109

Member Data Documentation

reco::CaloJetCollection HLTInclusiveVBFSource::calojet
private

Definition at line 107 of file HLTInclusiveVBFSource.h.

edm::Handle<reco::CaloJetCollection> HLTInclusiveVBFSource::calojetColl_
private

Definition at line 102 of file HLTInclusiveVBFSource.h.

edm::EDGetTokenT<reco::CaloJetCollection> HLTInclusiveVBFSource::caloJetsToken
private

Definition at line 99 of file HLTInclusiveVBFSource.h.

edm::Handle<reco::CaloMETCollection> HLTInclusiveVBFSource::calometColl_
private

Definition at line 103 of file HLTInclusiveVBFSource.h.

edm::EDGetTokenT<reco::CaloMETCollection> HLTInclusiveVBFSource::caloMetToken
private

Definition at line 100 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::check_mjj650_Pt35_DEta3p5
private

Definition at line 111 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::check_mjj650_Pt40_DEta3p5
private

Definition at line 115 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::check_mjj700_Pt35_DEta3p5
private

Definition at line 112 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::check_mjj700_Pt40_DEta3p5
private

Definition at line 116 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::check_mjj750_Pt35_DEta3p5
private

Definition at line 113 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::check_mjj750_Pt40_DEta3p5
private

Definition at line 117 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::check_mjj800_Pt35_DEta3p5
private

Definition at line 114 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::check_mjj800_Pt40_DEta3p5
private

Definition at line 118 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::checkdR_crossOrder
private

Definition at line 168 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::checkdR_sameOrder
private

Definition at line 167 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::checkHLT
private

Definition at line 160 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::checkHLTIndex
private

Definition at line 161 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::checkOffline
private

Definition at line 159 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::debug_
private

Definition at line 77 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::deltaRMatch_
private

Definition at line 83 of file HLTInclusiveVBFSource.h.

std::string HLTInclusiveVBFSource::dirname_
private

Definition at line 72 of file HLTInclusiveVBFSource.h.

float HLTInclusiveVBFSource::dR_HLT_RECO_11
private

Definition at line 163 of file HLTInclusiveVBFSource.h.

float HLTInclusiveVBFSource::dR_HLT_RECO_12
private

Definition at line 165 of file HLTInclusiveVBFSource.h.

float HLTInclusiveVBFSource::dR_HLT_RECO_21
private

Definition at line 166 of file HLTInclusiveVBFSource.h.

float HLTInclusiveVBFSource::dR_HLT_RECO_22
private

Definition at line 164 of file HLTInclusiveVBFSource.h.

bool HLTInclusiveVBFSource::etaOpposite_
private

Definition at line 84 of file HLTInclusiveVBFSource.h.

std::string HLTInclusiveVBFSource::filtername
private

Definition at line 121 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_deltaetajet
private

Definition at line 173 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_deltaphijet
private

Definition at line 174 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_ejet1
private

Definition at line 141 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_ejet2
private

Definition at line 150 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_etajet1
private

Definition at line 147 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_etajet2
private

Definition at line 156 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_invmassjet
private

Definition at line 175 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_phijet1
private

Definition at line 148 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_phijet2
private

Definition at line 157 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_ptjet1
private

Definition at line 146 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_ptjet2
private

Definition at line 155 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_pxjet1
private

Definition at line 143 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_pxjet2
private

Definition at line 152 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_pyjet1
private

Definition at line 144 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_pyjet2
private

Definition at line 153 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_pzjet1
private

Definition at line 145 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::hlt_pzjet2
private

Definition at line 154 of file HLTInclusiveVBFSource.h.

HLTConfigProvider HLTInclusiveVBFSource::hltConfig_
private

Definition at line 109 of file HLTInclusiveVBFSource.h.

PathInfoCollection HLTInclusiveVBFSource::hltPathsAll_
private

Definition at line 295 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::minDeltaEta_
private

Definition at line 81 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::minInvMass_
private

Definition at line 82 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::minPtHigh_
private

Definition at line 79 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::minPtLow_
private

Definition at line 80 of file HLTInclusiveVBFSource.h.

int HLTInclusiveVBFSource::nCount_
private

Definition at line 68 of file HLTInclusiveVBFSource.h.

std::vector<std::string> HLTInclusiveVBFSource::path_
private

Definition at line 75 of file HLTInclusiveVBFSource.h.

std::string HLTInclusiveVBFSource::pathname
private

Definition at line 120 of file HLTInclusiveVBFSource.h.

reco::PFJetCollection HLTInclusiveVBFSource::pfjet
private

Definition at line 108 of file HLTInclusiveVBFSource.h.

edm::Handle<reco::PFJetCollection> HLTInclusiveVBFSource::pfjetColl_
private

Definition at line 104 of file HLTInclusiveVBFSource.h.

edm::EDGetTokenT<edm::View<reco::PFJet> > HLTInclusiveVBFSource::pfJetsToken
private

Definition at line 97 of file HLTInclusiveVBFSource.h.

edm::Handle<reco::PFMETCollection> HLTInclusiveVBFSource::pfmetColl_
private

Definition at line 105 of file HLTInclusiveVBFSource.h.

edm::EDGetTokenT<edm::View<reco::PFMET> > HLTInclusiveVBFSource::pfMetToken
private

Definition at line 98 of file HLTInclusiveVBFSource.h.

std::vector<int> HLTInclusiveVBFSource::prescUsed_
private

Definition at line 70 of file HLTInclusiveVBFSource.h.

std::string HLTInclusiveVBFSource::processname_
private

Definition at line 73 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_deltaetajet
private

Definition at line 170 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_deltaphijet
private

Definition at line 171 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_ejet1
private

Definition at line 123 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_ejet2
private

Definition at line 132 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_etajet1
private

Definition at line 129 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_etajet2
private

Definition at line 138 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_invmassjet
private

Definition at line 172 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_phijet1
private

Definition at line 130 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_phijet2
private

Definition at line 139 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_ptjet1
private

Definition at line 128 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_ptjet2
private

Definition at line 137 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_pxjet1
private

Definition at line 125 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_pxjet2
private

Definition at line 134 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_pyjet1
private

Definition at line 126 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_pyjet2
private

Definition at line 135 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_pzjet1
private

Definition at line 127 of file HLTInclusiveVBFSource.h.

double HLTInclusiveVBFSource::reco_pzjet2
private

Definition at line 136 of file HLTInclusiveVBFSource.h.

edm::TriggerNames HLTInclusiveVBFSource::triggerNames_
private

Definition at line 90 of file HLTInclusiveVBFSource.h.

edm::Handle<trigger::TriggerEvent> HLTInclusiveVBFSource::triggerObj_
private

Definition at line 87 of file HLTInclusiveVBFSource.h.

edm::Handle<edm::TriggerResults> HLTInclusiveVBFSource::triggerResults_
private

Definition at line 89 of file HLTInclusiveVBFSource.h.

edm::EDGetTokenT<edm::TriggerResults> HLTInclusiveVBFSource::triggerResultsFUToken
private

Definition at line 93 of file HLTInclusiveVBFSource.h.

edm::InputTag HLTInclusiveVBFSource::triggerResultsLabel_
private

Definition at line 88 of file HLTInclusiveVBFSource.h.

edm::EDGetTokenT<edm::TriggerResults> HLTInclusiveVBFSource::triggerResultsToken
private

Definition at line 92 of file HLTInclusiveVBFSource.h.

edm::EDGetTokenT<trigger::TriggerEvent> HLTInclusiveVBFSource::triggerSummaryFUToken
private

Definition at line 95 of file HLTInclusiveVBFSource.h.

edm::InputTag HLTInclusiveVBFSource::triggerSummaryLabel_
private

Definition at line 86 of file HLTInclusiveVBFSource.h.

edm::EDGetTokenT<trigger::TriggerEvent> HLTInclusiveVBFSource::triggerSummaryToken
private

Definition at line 94 of file HLTInclusiveVBFSource.h.

Generated for CMSSW Reference Manual by doxygen 1.8.5