Get the analysis.
411 int myLuminosityBlock;
415 if (myLuminosityBlock<
_LSBegin)
return;
422 bool bPhysicsDeclared =
false;
428 if( triggerNames.
triggerIndex(
"HLT_PhysicsDeclared") != triggerNames.
size() )
429 if (triggerResults->accept(triggerNames.
triggerIndex(
"HLT_PhysicsDeclared"))) {
437 if (triggerResults->accept(
j)){
482 std::cout <<
">>> Trigger Lo = " << JetLoPass
483 <<
" Hi = " << JetHiPass
490 bool bPrimaryVertex =
true;
492 bPrimaryVertex =
false;
498 LogInfo(
"JetMETAnalyzer") <<
"JetMETAnalyzer: Could not find vertex collection" << std::endl;
499 if (
DEBUG)
std::cout <<
"JetMETAnalyzer: Could not find vertex collection" << std::endl;
504 int vertex_number = vertexCollection.size();
505 reco::VertexCollection::const_iterator
v = vertexCollection.begin();
506 double vertex_chi2 = v->normalizedChi2();
507 double vertex_ndof = v->ndof();
508 bool fakeVtx = v->isFake();
509 double vertex_Z = v->z();
515 && fabs(vertex_Z)<
_vtxz_max ) bPrimaryVertex =
true;
523 if (!gtReadoutRecord.
isValid()) {
524 LogInfo(
"JetMETAnalyzer") <<
"JetMETAnalyzer: Could not find GT readout record" << std::endl;
525 if (
DEBUG)
std::cout <<
"JetMETAnalyzer: Could not find GT readout record product" << std::endl;
528 bool bTechTriggers =
true;
529 bool bTechTriggersAND =
true;
530 bool bTechTriggersOR =
false;
531 bool bTechTriggersNOT =
false;
533 if (gtReadoutRecord.
isValid()) {
534 const TechnicalTriggerWord& technicalTriggerWordBeforeMask = gtReadoutRecord->technicalTriggerWord();
537 bTechTriggersAND =
true;
540 bTechTriggersAND = bTechTriggersAND && technicalTriggerWordBeforeMask.at(
_techTrigsAND.at(ttr));
544 bTechTriggersOR =
true;
546 for (
unsigned ttr = 0; ttr !=
_techTrigsOR.size(); ttr++) {
547 bTechTriggersOR = bTechTriggersOR || technicalTriggerWordBeforeMask.at(
_techTrigsOR.at(ttr));
550 bTechTriggersNOT =
false;
553 bTechTriggersNOT = bTechTriggersNOT || technicalTriggerWordBeforeMask.at(
_techTrigsNOT.at(ttr));
558 bTechTriggersAND =
true;
559 bTechTriggersOR =
true;
560 bTechTriggersNOT =
false;
564 bTechTriggersAND =
true;
566 bTechTriggersOR =
true;
568 bTechTriggersNOT =
false;
570 bTechTriggers = bTechTriggersAND && bTechTriggersOR && !bTechTriggersNOT;
572 bool bJetCleanup = bTechTriggers && bPrimaryVertex && bPhysicsDeclared;
589 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Jet Pt anti-Kt analyzer";
595 if(caloJets.
isValid() && bJetCleanup) {
603 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned Jet Pt anti-Kt analyzer";
609 if(caloJets.
isValid() && bJetCleanup){
620 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Jet Pt anti-Kt analyzer";
625 if(caloJets.
isValid() && bJetCleanup){
628 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned Jet Pt anti-Kt analyzer";
644 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Jet Pt SisCone analyzer";
649 if(caloJets.
isValid() && bJetCleanup){
657 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned Jet Pt SisCone analyzer";
673 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Jet Pt ICone analyzer";
678 if(caloJets.
isValid() && bJetCleanup){
686 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned Jet Pt ICone analyzer";
725 LogTrace(
metname)<<
"[JetMETAnalyzer] Call to the Cleaned PFJet analyzer";
unsigned int size() const
number of trigger paths in trigger table
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
void analyze(const edm::Event &, const edm::EventSetup &, const reco::PFJetCollection &pfJets)
Get 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 analyze(const edm::Event &event, const edm::EventSetup &eventSetup, const reco::JPTJet &jptJet, double &pt1, double &pt2, double &pt3)
Do the analysis.
void setJetLoPass(int pass)
unsigned int triggerIndex(std::string const &name) const
void setJetLoPass(int pass)
bool accept(const edm::Event &event, const edm::EventSetup &setup)
To be called from analyze/filter() methods.
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)
T const * product() const
void analyze(const edm::Event &, const edm::EventSetup &, const reco::CaloJetCollection &caloJets)
Get the analysis.