L1TNtuples/src/L1AnalysisRecoTau.cc

0001 #include "L1Trigger/L1TNtuples/interface/L1AnalysisRecoTau.h"
0002
0003 //#include "JetMETCorrections/Objects/interface/JetCorrector.h"
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 // void L1Analysis::L1AnalysisRecoTau::SetPFJet(const edm::Event& event,
0053 //                  const edm::EventSetup& setup,
0054 //                  edm::Handle<reco::PFJetCollection> taus,
0055 //                  unsigned maxTau)
0056 // {
0057
0058 //   recoPFJet_.nTaus=0;
0059
0060 //   for(reco::CaloJetCollection::const_iterator it=taus->begin();
0061 //       it!=taus->end() && recoTau_.nTaus < maxTau;
0062 //       ++it) {
0063
0064 //     recoTau_.et.push_back(it->et());
0065 //     // recoTau_.etCorr.push_back(it->et());// * scale);
0066 //     // recoTau_.corrFactor.push_back(1.);//scale);
0067 //     recoTau_.eta.push_back(it->eta());
0068 //     recoTau_.phi.push_back(it->phi());
0069 //     recoTau_.e.push_back(it->energy());
0070 //     // recoTau_.eEMF.push_back(it->emEnergyFraction());
0071 //     // recoTau_.eEmEB.push_back(it->emEnergyInEB());
0072 //     // recoTau_.eEmEE.push_back(it->emEnergyInEE());
0073 //     // recoTau_.eEmHF.push_back(it->emEnergyInHF());
0074 //     // recoTau_.eHadHB.push_back(it->hadEnergyInHB());
0075 //     // recoTau_.eHadHE.push_back(it->hadEnergyInHE());
0076 //     // recoTau_.eHadHO.push_back(it->hadEnergyInHO());
0077 //     // recoTau_.eHadHF.push_back(it->hadEnergyInHF());
0078 //     // recoTau_.eMaxEcalTow.push_back(it->maxEInEmTowers());
0079 //     // recoTau_.eMaxHcalTow.push_back(it->maxEInHadTowers());
0080 //     // recoTau_.towerArea.push_back(it->towersArea());
0081 //     //    recoTau_.towerSize.push_back(static_cast<int>(it->getCaloConstituents().size()));
0082 //     // recoTau_.n60.push_back(it->n60());
0083 //     // recoTau_.n90.push_back(it->n90());
0084
0085 //     // recoTau_.n90hits.push_back(1.); //int((*jetsID)[jetRef].n90Hits));
0086 //     // recoTau_.fHPD.push_back(1.); //(*jetsID)[jetRef].fHPD);
0087 //     // recoTau_.fRBX.push_back(1.); //(*jetsID)[jetRef].fRBX);
0088
0089 //     recoTau_.nTaus++;
0090
0091 //   }
0092 // }