368 std::cout <<
"Inside analyze(). " << std::endl;
378 if (!aodTriggerEvent.
isValid())
382 unordered_map<string, bool> firedMap =
acceptMap;
385 if (firedMap[
plot->pathName])
393 pathName =
plot->pathName;
394 triggerAccept =
true;
395 }
else if (
plot->pathIndexOR > 0 && triggerResults->
accept(
plot->pathIndexOR) &&
397 if (firedMap[
plot->pathNameOR])
400 pathName =
plot->pathNameOR;
401 triggerAccept =
true;
405 unsigned int moduleFilterIndex = aodTriggerEvent->
filterIndex(moduleFilter);
407 if (moduleFilterIndex + 1 > aodTriggerEvent->
sizeFilters())
419 for (
const auto&
key : keys)
425 for (
const auto&
key : keys)
431 for (
const auto&
key : keys) {
440 for (
const auto&
key : keys) {
449 for (
const auto&
key : keys) {
458 for (
const auto&
key : keys) {
467 for (
const auto&
key : keys) {
476 double onlineMR = 0, onlineRsq = 0;
477 for (
const auto&
key : keys) {
478 if (objects[
key].
id() == 0) {
479 onlineMR = objects[
key].px();
480 onlineRsq = objects[
key].py();
489 std::vector<ROOT::Math::LorentzVector<ROOT::Math::PtEtaPhiM4D<double>>> alphaT_jets;
490 for (
const auto&
key : keys) {
491 ROOT::Math::LorentzVector<ROOT::Math::PtEtaPhiM4D<double>> JetLVec(
493 alphaT_jets.push_back(JetLVec);
496 float alphaT = AlphaT(alphaT_jets,
false).value();
502 for (
const auto&
key : keys)
508 for (
const auto&
key : keys) {
516 for (
const auto&
key : keys) {
517 if (objects[
key].
id() == 89)
524 for (
const auto&
key : keys) {
532 for (
const auto&
key : keys) {
540 for (
const auto&
key : keys) {
554 for (
auto iter = csvPfTags->
begin(); iter != csvPfTags->
end(); iter++)
565 for (
auto iter = csvCaloTags->
begin(); iter != csvCaloTags->
end(); iter++)
579 for (
const auto&
key : keys) {
581 if (muon_dxy != -99.)
595 const double mu_mass(.105658);
596 unsigned int kCnt0 = 0;
597 for (
const auto& key0 : keys) {
598 unsigned int kCnt1 = 0;
599 for (
const auto& key1 : keys) {
603 if (
abs(objects[key0].
id()) == 13 &&
604 (objects[key0].
id() + objects[key1].
id() == 0))
606 TLorentzVector mu1, mu2, dimu;
607 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
608 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
620 const double mu_mass(.105658);
621 unsigned int kCnt0 = 0;
622 for (
const auto& key0 : keys) {
623 unsigned int kCnt1 = 0;
624 for (
const auto& key1 : keys) {
628 if (
abs(objects[key0].
id()) == 13 &&
629 (objects[key0].
id() + objects[key1].
id() == 0))
631 TLorentzVector mu1, mu2, dimu;
632 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
633 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
645 const double mu_mass(.105658);
646 unsigned int kCnt0 = 0;
647 for (
const auto& key0 : keys) {
648 unsigned int kCnt1 = 0;
649 for (
const auto& key1 : keys) {
655 TLorentzVector mu1, mu2, dimu;
656 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
657 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
670 const double mu_mass(.105658);
671 unsigned int kCnt0 = 0;
672 for (
const auto& key0 : keys) {
673 unsigned int kCnt1 = 0;
674 for (
const auto& key1 : keys) {
678 if (
abs(objects[key0].
id()) == 13 &&
679 (objects[key0].
id() + objects[key1].
id() == 0))
681 TLorentzVector mu1, mu2, dimu;
682 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
683 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
696 const double mu_mass(.105658);
697 unsigned int kCnt0 = 0;
698 for (
const auto& key0 : keys) {
699 unsigned int kCnt1 = 0;
700 for (
const auto& key1 : keys) {
704 if (
abs(objects[key0].
id()) == 13 &&
705 (objects[key0].
id() + objects[key1].
id() == 0))
707 TLorentzVector mu1, mu2, dimu;
708 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), mu_mass);
709 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), mu_mass);
722 unsigned int kCnt0 = 0;
723 for (
const auto& key0 : keys) {
724 unsigned int kCnt1 = 0;
725 for (
const auto& key1 : keys) {
731 TLorentzVector el1, el2, diEl;
732 el1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), 0);
733 el2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(), 0);
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
1.8.11
