367 std::cout <<
"Inside analyze(). " << std::endl;
377 if (!aodTriggerEvent.isValid())
381 unordered_map<string, bool> firedMap =
acceptMap;
384 if (firedMap[
plot->pathName])
386 bool triggerAccept =
false;
392 pathName =
plot->pathName;
393 triggerAccept =
true;
396 if (firedMap[
plot->pathNameOR])
399 pathName =
plot->pathNameOR;
400 triggerAccept =
true;
404 unsigned int moduleFilterIndex = aodTriggerEvent->filterIndex(moduleFilter);
406 if (moduleFilterIndex + 1 > aodTriggerEvent->sizeFilters())
408 const Keys&
keys = aodTriggerEvent->filterKeys(moduleFilterIndex);
418 for (
const auto&
key : keys)
424 for (
const auto&
key : keys)
430 for (
const auto&
key : keys) {
439 for (
const auto&
key : keys) {
448 for (
const auto&
key : keys) {
457 for (
const auto&
key : keys) {
466 for (
const auto&
key : keys) {
475 double onlineMR = 0, onlineRsq = 0;
476 for (
const auto&
key : keys) {
477 if (objects[
key].
id() == 0) {
478 onlineMR = objects[
key].px();
479 onlineRsq = objects[
key].py();
488 std::vector<ROOT::Math::LorentzVector<ROOT::Math::PtEtaPhiM4D<double>>> alphaT_jets;
489 for (
const auto&
key : keys) {
490 ROOT::Math::LorentzVector<ROOT::Math::PtEtaPhiM4D<double>> JetLVec(
492 alphaT_jets.push_back(JetLVec);
495 float alphaT = AlphaT(alphaT_jets,
false).value();
501 for (
const auto&
key : keys)
507 for (
const auto&
key : keys) {
515 for (
const auto&
key : keys) {
516 if (objects[
key].
id() == 89)
523 for (
const auto&
key : keys) {
531 for (
const auto&
key : keys) {
539 for (
const auto&
key : keys) {
553 for (
auto iter = csvPfTags->begin(); iter != csvPfTags->end(); iter++)
564 for (
auto iter = csvCaloTags->begin(); iter != csvCaloTags->end(); iter++)
578 for (
const auto&
key : keys) {
580 if (muon_dxy != -99.)
595 unsigned int kCnt0 = 0;
596 for (
const auto& key0 : keys) {
597 unsigned int kCnt1 = 0;
598 for (
const auto& key1 : keys) {
602 if (
abs(objects[key0].
id()) == 13 &&
603 (objects[key0].
id() + objects[key1].
id() == 0))
605 TLorentzVector mu1, mu2, dimu;
606 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(),
mu_mass);
607 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(),
mu_mass);
620 unsigned int kCnt0 = 0;
621 for (
const auto& key0 : keys) {
622 unsigned int kCnt1 = 0;
623 for (
const auto& key1 : keys) {
627 if (
abs(objects[key0].
id()) == 13 &&
628 (objects[key0].
id() + objects[key1].
id() == 0))
630 TLorentzVector mu1, mu2, dimu;
631 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(),
mu_mass);
632 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(),
mu_mass);
645 unsigned int kCnt0 = 0;
646 for (
const auto& key0 : keys) {
647 unsigned int kCnt1 = 0;
648 for (
const auto& key1 : keys) {
654 TLorentzVector mu1, mu2, dimu;
655 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(),
mu_mass);
656 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(),
mu_mass);
670 unsigned int kCnt0 = 0;
671 for (
const auto& key0 : keys) {
672 unsigned int kCnt1 = 0;
673 for (
const auto& key1 : keys) {
677 if (
abs(objects[key0].
id()) == 13 &&
678 (objects[key0].
id() + objects[key1].
id() == 0))
680 TLorentzVector mu1, mu2, dimu;
681 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(),
mu_mass);
682 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(),
mu_mass);
696 unsigned int kCnt0 = 0;
697 for (
const auto& key0 : keys) {
698 unsigned int kCnt1 = 0;
699 for (
const auto& key1 : keys) {
703 if (
abs(objects[key0].
id()) == 13 &&
704 (objects[key0].
id() + objects[key1].
id() == 0))
706 TLorentzVector mu1, mu2, dimu;
707 mu1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(),
mu_mass);
708 mu2.SetPtEtaPhiM(objects[key1].
pt(), objects[key1].
eta(), objects[key1].
phi(),
mu_mass);
721 unsigned int kCnt0 = 0;
722 for (
const auto& key0 : keys) {
723 unsigned int kCnt1 = 0;
724 for (
const auto& key1 : keys) {
730 TLorentzVector el1, el2, diEl;
731 el1.SetPtEtaPhiM(objects[key0].
pt(), objects[key0].
eta(), objects[key0].
phi(), 0);
732 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)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken_
edm::EDGetTokenT< vector< reco::PFJet > > csvPfJetsToken_
hltPlot pAL2DoubleMuZMass_
edm::EDGetTokenT< trigger::TriggerEvent > aodTriggerToken_
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > chargedCandToken_
tuple key
prepare the HTCondor submission files and eventually submit them
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
static constexpr float mu_mass
unordered_map< string, bool > acceptMap