File indexing completed on 2024-04-06 12:21:16
0001 #include "L1Trigger/L1TNtuples/interface/L1AnalysisRecoTau.h"
0002
0003
0004 #include "JetMETCorrections/JetCorrector/interface/JetCorrector.h"
0005
0006 using namespace std;
0007
0008 L1Analysis::L1AnalysisRecoTau::L1AnalysisRecoTau() {}
0009
0010 L1Analysis::L1AnalysisRecoTau::~L1AnalysisRecoTau() {}
0011
0012 void L1Analysis::L1AnalysisRecoTau::SetTau(const edm::Event& event,
0013 const edm::EventSetup& setup,
0014 edm::Handle<reco::PFTauCollection> taus,
0015 edm::Handle<reco::PFTauDiscriminator> DMFindingOldTaus,
0016 edm::Handle<reco::PFTauDiscriminator> DMFindingTaus,
0017 edm::Handle<reco::PFTauDiscriminator> TightIsoTaus,
0018 edm::Handle<reco::PFTauDiscriminator> TightRawIsoTaus,
0019 edm::Handle<reco::PFTauDiscriminator> LooseIsoTaus,
0020 edm::Handle<reco::PFTauDiscriminator> LooseAntiMuon,
0021 edm::Handle<reco::PFTauDiscriminator> TightAntiMuon,
0022 edm::Handle<reco::PFTauDiscriminator> VLooseAntiElectron,
0023 edm::Handle<reco::PFTauDiscriminator> LooseAntiElectron,
0024 edm::Handle<reco::PFTauDiscriminator> TightAntiElectron,
0025 unsigned maxTau) {
0026 recoTau_.nTaus = 0;
0027
0028 for (reco::PFTauCollection::const_iterator it = taus->begin(); it != taus->end() && recoTau_.nTaus < maxTau; ++it) {
0029 recoTau_.e.push_back(it->energy());
0030 recoTau_.pt.push_back(it->pt());
0031 recoTau_.et.push_back(it->et());
0032 recoTau_.eta.push_back(it->eta());
0033 recoTau_.phi.push_back(it->phi());
0034
0035 edm::Ref<reco::PFTauCollection> tauEdmRef(taus, recoTau_.nTaus);
0036 recoTau_.TightIsoFlag.push_back((*TightIsoTaus)[tauEdmRef]);
0037 recoTau_.RawIso.push_back((*TightRawIsoTaus)[tauEdmRef]);
0038 recoTau_.LooseIsoFlag.push_back((*LooseIsoTaus)[tauEdmRef]);
0039 recoTau_.LooseAntiMuonFlag.push_back((*LooseAntiMuon)[tauEdmRef]);
0040 recoTau_.TightAntiMuonFlag.push_back((*TightAntiMuon)[tauEdmRef]);
0041 recoTau_.VLooseAntiElectronFlag.push_back((*VLooseAntiElectron)[tauEdmRef]);
0042 recoTau_.LooseAntiElectronFlag.push_back((*LooseAntiElectron)[tauEdmRef]);
0043 recoTau_.TightAntiElectronFlag.push_back((*TightAntiElectron)[tauEdmRef]);
0044 recoTau_.DMFindingOldDMs.push_back((*DMFindingOldTaus)[tauEdmRef]);
0045 recoTau_.DMFindingNewDMs.push_back((*DMFindingTaus)[tauEdmRef]);
0046 recoTau_.charge.push_back(it->charge());
0047
0048 recoTau_.nTaus++;
0049 }
0050 }
0051
0052
0053
0054
0055
0056
0057
0058
0059
0060
0061
0062
0063
0064
0065
0066
0067
0068
0069
0070
0071
0072
0073
0074
0075
0076
0077
0078
0079
0080
0081
0082
0083
0084
0085
0086
0087
0088
0089
0090
0091
0092