Get the analysis.
437 int myLuminosityBlock;
441 if (myLuminosityBlock<
_LSBegin)
return;
448 bool bPhysicsDeclared =
false;
457 if( triggerNames.
triggerIndex(
"HLT_PhysicsDeclared") != triggerNames.
size() )
458 if (triggerResults->accept(triggerNames.
triggerIndex(
"HLT_PhysicsDeclared"))) {
473 const unsigned int nTrig(triggerNames.
size());
474 for (
unsigned int i=0;
i<nTrig;++
i)
495 std::cout <<
">>> Trigger Lo = " << JetLoPass
496 <<
" Hi = " << JetHiPass
504 bool bPrimaryVertex =
true;
506 bPrimaryVertex =
false;
512 LogInfo(
"JetMETAnalyzer") <<
"JetMETAnalyzer: Could not find vertex collection" << std::endl;
513 if (
DEBUG)
std::cout <<
"JetMETAnalyzer: Could not find vertex collection" << std::endl;
518 int vertex_number = vertexCollection.size();
519 reco::VertexCollection::const_iterator
v = vertexCollection.begin();
520 for ( ; v != vertexCollection.end(); ++
v) {
521 double vertex_chi2 = v->normalizedChi2();
522 double vertex_ndof = v->ndof();
523 bool fakeVtx = v->isFake();
524 double vertex_Z = v->z();
531 bPrimaryVertex =
true;
543 if (!gtReadoutRecord.
isValid()) {
544 LogInfo(
"JetMETAnalyzer") <<
"JetMETAnalyzer: Could not find GT readout record" << std::endl;
545 if (
DEBUG)
std::cout <<
"JetMETAnalyzer: Could not find GT readout record product" << std::endl;
548 bool bTechTriggers =
true;
549 bool bTechTriggersAND =
true;
550 bool bTechTriggersOR =
false;
551 bool bTechTriggersNOT =
false;
553 if (gtReadoutRecord.
isValid()) {
554 const TechnicalTriggerWord& technicalTriggerWordBeforeMask = gtReadoutRecord->technicalTriggerWord();
557 bTechTriggersAND =
true;
560 bTechTriggersAND = bTechTriggersAND && technicalTriggerWordBeforeMask.at(
_techTrigsAND.at(ttr));
564 bTechTriggersOR =
true;
566 for (
unsigned ttr = 0; ttr !=
_techTrigsOR.size(); ttr++) {
567 bTechTriggersOR = bTechTriggersOR || technicalTriggerWordBeforeMask.at(
_techTrigsOR.at(ttr));
570 bTechTriggersNOT =
false;
573 bTechTriggersNOT = bTechTriggersNOT || technicalTriggerWordBeforeMask.at(
_techTrigsNOT.at(ttr));
578 bTechTriggersAND =
true;
579 bTechTriggersOR =
true;
580 bTechTriggersNOT =
false;
584 bTechTriggersAND =
true;
586 bTechTriggersOR =
true;
588 bTechTriggersNOT =
false;
590 bTechTriggers = bTechTriggersAND && bTechTriggersOR && !bTechTriggersNOT;
592 bool bJetCleanup = bTechTriggers && bPrimaryVertex && bPhysicsDeclared;
619 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Jet Pt anti-Kt analyzer";
625 if(caloJets.
isValid() && bJetCleanup) {
633 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned Jet Pt anti-Kt analyzer";
639 if(caloJets.
isValid() && bJetCleanup){
650 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Jet Pt anti-Kt analyzer";
655 if(caloJets.
isValid() && bJetCleanup){
658 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned Jet Pt anti-Kt analyzer";
674 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Jet Pt SisCone analyzer";
679 if(caloJets.
isValid() && bJetCleanup){
687 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned Jet Pt SisCone analyzer";
703 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Jet Pt ICone analyzer";
708 if(caloJets.
isValid() && bJetCleanup){
716 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned Jet Pt ICone analyzer";
755 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned PFJet analyzer";
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
void analyze(const edm::Event &event, const edm::EventSetup &eventSetup, const reco::JPTJet &jptJet, double &pt1, double &pt2, double &pt3, const int numPV)
Do the analysis.
edm::LuminosityBlockNumber_t luminosityBlock() const
void analyze(const edm::Event &, const edm::EventSetup &, const edm::TriggerResults &)
Get the analysis.
std::vector< Vertex > VertexCollection
collection of Vertex objects
Strings::size_type size() const
void setJetLoPass(int pass)
unsigned int triggerIndex(std::string const &name) const
void setJetLoPass(int pass)
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
void analyze(const edm::Event &, const edm::EventSetup &, const reco::CaloJetCollection &caloJets)
Get the analysis.
void setJetHiPass(int pass)
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
bool filter(const edm::Event &evt, const edm::EventSetup &es)
void setJetHiPass(int pass)
std::string const & triggerName(unsigned int index) const
T const * product() const
void analyze(const edm::Event &, const edm::EventSetup &, const reco::PFJetCollection &pfJets, const int numPV)
Get the analysis.
void analyze(const edm::Event &, const edm::EventSetup &, const reco::CaloJetCollection &caloJets, const int numPV)
Get the analysis.