387 if (!triggerResults.
isValid())
return;
391 if (!aodTriggerEvent.
isValid())
return;
394 unordered_map<string, bool> firedMap =
acceptMap;
397 if (firedMap[
plot->pathName])
continue;
398 bool triggerAccept =
false;
405 pathName =
plot->pathName;
406 triggerAccept =
true;
410 if (firedMap[
plot->pathNameOR])
continue;
412 pathName =
plot->pathNameOR;
413 triggerAccept =
true;
418 unsigned int moduleFilterIndex = aodTriggerEvent->
filterIndex(moduleFilter);
420 if (moduleFilterIndex+1 > aodTriggerEvent->
sizeFilters())
return;
444 for (
const auto &
key : keys)
455 for (
const auto &
key : keys)
466 for (
const auto &
key : keys)
477 for (
const auto &
key : keys)
488 for (
const auto &
key : keys)
499 double onlineMR = 0, onlineRsq = 0;
500 for (
const auto &
key : keys)
502 if(objects[
key].
id() == 0){
503 onlineMR = objects[
key].px();
504 onlineRsq = objects[
key].py();
514 std::vector<ROOT::Math::LorentzVector<ROOT::Math::PtEtaPhiM4D<double>>> alphaT_jets;
515 for (
const auto &
key : keys)
517 ROOT::Math::LorentzVector<ROOT::Math::PtEtaPhiM4D<double>> JetLVec(objects[
key].
pt(),objects[
key].
eta(),objects[
key].
phi(),objects[
key].
mass());
518 alphaT_jets.push_back(JetLVec);
521 float alphaT = AlphaT(alphaT_jets,
false).value();
534 for (
const auto &
key : keys)
544 for (
const auto &
key : keys)
553 for (
const auto &
key : keys)
563 for (
const auto &
key : keys)
573 for (
const auto &
key : keys)
617 for (
const auto &
key : keys)
634 const double mu_mass(.105658);
635 unsigned int kCnt0 = 0;
636 for (
const auto & key0: keys)
638 unsigned int kCnt1 = 0;
639 for (
const auto & key1: keys)
641 if (key0 != key1 && kCnt1 > kCnt0)
643 if (
abs(objects[key0].
id()) == 13 && (objects[key0].
id()+objects[key1].
id()==0))
645 TLorentzVector mu1, mu2, dimu;
646 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
647 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
660 const double mu_mass(.105658);
661 unsigned int kCnt0 = 0;
662 for (
const auto & key0: keys)
664 unsigned int kCnt1 = 0;
665 for (
const auto & key1: keys)
667 if (key0 != key1 && kCnt1 > kCnt0)
669 if (
abs(objects[key0].
id()) == 13 && (objects[key0].
id()+objects[key1].
id()==0))
671 TLorentzVector mu1, mu2, dimu;
672 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
673 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
686 const double mu_mass(.105658);
687 unsigned int kCnt0 = 0;
688 for (
const auto & key0: keys)
690 unsigned int kCnt1 = 0;
691 for (
const auto & key1: keys)
693 if (key0 != key1 && kCnt1 > kCnt0)
697 TLorentzVector mu1, mu2, dimu;
698 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
699 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
712 const double mu_mass(.105658);
713 unsigned int kCnt0 = 0;
714 for (
const auto & key0: keys)
716 unsigned int kCnt1 = 0;
717 for (
const auto & key1: keys)
719 if (key0 != key1 && kCnt1 > kCnt0)
721 if (
abs(objects[key0].
id()) == 13 && (objects[key0].
id()+objects[key1].
id()==0))
723 TLorentzVector mu1, mu2, dimu;
724 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
725 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
738 const double mu_mass(.105658);
739 unsigned int kCnt0 = 0;
740 for (
const auto & key0: keys)
742 unsigned int kCnt1 = 0;
743 for (
const auto & key1: keys)
745 if (key0 != key1 && kCnt1 > kCnt0)
747 if (
abs(objects[key0].
id()) == 13 && (objects[key0].
id()+objects[key1].
id()==0))
749 TLorentzVector mu1, mu2, dimu;
750 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
751 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
764 unsigned int kCnt0 = 0;
765 for (
const auto & key0: keys)
767 unsigned int kCnt1 = 0;
768 for (
const auto & key1: keys)
770 if (key0 != key1 && kCnt1 > kCnt0)
774 TLorentzVector el1, el2, diEl;
775 el1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), 0);
776 el2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), 0);
787 firedMap[pathName] =
true;
793 double wallTime = end -
start;
bool saveTags(const std::string &module) const
Is module an L3 filter (ie, tracked saveTags=true)
trigger::size_type sizeFilters() const
bool getByToken(EDGetToken token, Handle< PROD > &result) const
bool accept() const
Has at least one path accepted the event?
const Keys & filterKeys(trigger::size_type index) const
trigger::size_type filterIndex(const edm::InputTag &filterTag) const
find index of filter in data-member vector from filter tag
const_iterator end() const
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken_
edm::EDGetTokenT< vector< reco::PFJet > > csvPfJetsToken_
hltPlot pAL2DoubleMuZMass_
edm::EDGetTokenT< trigger::TriggerEvent > aodTriggerToken_
const TriggerObjectCollection & getObjects() const
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > chargedCandToken_
Abs< T >::type abs(const T &t)
static std::string const triggerResults
hltPlot pAL1DoubleMuZMass_
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
hltPlot pAL3DoubleMuZMass_
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
HLTConfigProvider hltConfig_
double dxyFinder(double, double, edm::Handle< reco::RecoChargedCandidateCollection >, edm::Handle< reco::BeamSpot >)
std::vector< size_type > Keys
edm::EDGetTokenT< reco::JetTagCollection > csvCaloTagsToken_
edm::EDGetTokenT< reco::JetTagCollection > csvPfTagsToken_
double get_wall_time(void)
edm::EDGetTokenT< vector< reco::CaloJet > > csvCaloJetsToken_
vector< hltPlot * > plotList
const_iterator begin() const
unordered_map< string, bool > acceptMap