120 pfJetsMissing_(
false),
121 pfMetMissing_(
false),
122 pfJetCorrMissing_(
false),
123 caloMetMissing_(
false),
124 caloMetBEMissing_(
false)
145 tree_->Branch(
"Jet",
"L1Analysis::L1AnalysisRecoJetDataFormat", &jet_data, 32000, 3);
146 tree_->Branch(
"Sums",
"L1Analysis::L1AnalysisRecoMetDataFormat", &
met_data, 32000, 3);
253 for(
auto it=caloJets->begin();
287 for(
auto it=pfJets->begin();
296 jet_data->
chef.push_back(it->chargedHadronEnergyFraction());
297 jet_data->
nhef.push_back(it->neutralHadronEnergyFraction());
298 jet_data->
pef.push_back(it->photonEnergyFraction());
299 jet_data->
eef.push_back(it->electronEnergyFraction());
312 jet_data->
cmef.push_back(it->chargedMuEnergyFraction());
328 float corrFactor = 1.;
339 for(
auto it=pfJets->begin();
340 it!=pfJets->end() && nJets <
maxJet_;
343 corrFactor = pfJetCorr.
product()->correction(*it);
351 mHx += -1.*it->px()*corrFactor;
352 mHy += -1.*it->py()*corrFactor;
358 TVector2 *tv2 =
new TVector2(mHx,mHy);
433 if (fabs(jet.
eta()) < 2.4) {
438 if (fabs(jet.
eta()) > 3.0) {
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
void doPFMet(edm::Handle< reco::PFMETCollection > pfMet)
L1Analysis::L1AnalysisRecoMetDataFormat * met_data
tuple caloMet
____________________________________________________________________________||
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::Service< TFileService > fs_
#define DEFINE_FWK_MODULE(type)
float chargedHadronEnergyFraction() const
chargedHadronEnergyFraction
void doPFJetCorr(edm::Handle< reco::PFJetCollection > pfJets, edm::Handle< reco::JetCorrector > pfJetCorr)
edm::EDGetTokenT< reco::CaloMETCollection > caloMetToken_
virtual double phi() const final
momentum azimuthal angle
void doPFJets(edm::Handle< reco::PFJetCollection > pfJets)
T * make(const Args &...args) const
make new ROOT object
void doCaloMetBE(edm::Handle< reco::CaloMETCollection > caloMetBE)
int chargedMultiplicity() const
chargedMultiplicity
Jets made from PFObjects.
L1Analysis::L1AnalysisRecoJetDataFormat * jet_data
float neutralEmEnergyFraction() const
neutralEmEnergyFraction
edm::EDGetTokenT< reco::JetCorrector > jecToken_
tuple pfMet
____________________________________________________________________________||
float neutralHadronEnergyFraction() const
neutralHadronEnergyFraction
edm::EDGetTokenT< reco::PFJetCollection > pfJetToken_
virtual void analyze(const edm::Event &, const edm::EventSetup &)
void doCaloMet(edm::Handle< reco::CaloMETCollection > caloMet)
int neutralMultiplicity() const
neutralMultiplicity
edm::EDGetTokenT< reco::PFMETCollection > pfMetToken_
float chargedEmEnergyFraction() const
chargedEmEnergyFraction
std::vector< reco::CaloMET > CaloMETCollection
collection of CaloMET objects
T const * product() const
L1JetRecoTreeProducer(const edm::ParameterSet &)
void doCaloJets(edm::Handle< reco::CaloJetCollection > caloJets)
int muonMultiplicity() const
muonMultiplicity
std::vector< std::vector< double > > tmp
virtual void beginJob(void)
std::vector< reco::PFMET > PFMETCollection
collection of PFMET objects
bool jetId(const reco::PFJet &jet)
virtual double et() const final
transverse energy
virtual double eta() const final
momentum pseudorapidity
volatile std::atomic< bool > shutdown_flag false
float muonEnergyFraction() const
muonEnergyFraction
int electronMultiplicity() const
electronMultiplicity
edm::EDGetTokenT< reco::CaloMETCollection > caloMetBEToken_
tuple caloMetBE
____________________________________________________________________________||