|
|
#include <PATMHTProducer.h>
Definition at line 60 of file PATMHTProducer.h.
◆ PATMHTProducer()
Definition at line 8 of file PATMHTProducer.cc.
51 produces<pat::MHTCollection>();
References CaloTowerTag_, controlledUncertainty_, eleEtaMax_, eleEtUncertaintyParameter0_, elePhiUncertaintyParameter0_, elePtMin_, eleToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), jetEMfracMax_, jetEtaMax_, jetEtUncertaintyParameter0_, jetEtUncertaintyParameter1_, jetEtUncertaintyParameter2_, jetPhiUncertaintyParameter0_, jetPhiUncertaintyParameter1_, jetPhiUncertaintyParameter2_, jetPtMin_, jetToken_, muonEtaMax_, muonEtUncertaintyParameter0_, muonPhiUncertaintyParameter0_, muonPtMin_, muoToken_, noHF_, phoToken_, setUncertaintyParameters(), tauToken_, towerEtThreshold_, uncertaintyScaleFactor_, useHO_, and verbose_.
◆ ~PATMHTProducer()
PATMHTProducer::~PATMHTProducer |
( |
| ) |
|
|
override |
◆ getElectrons()
Definition at line 197 of file PATMHTProducer.cc.
200 double number_of_electrons_ = 0.0;
218 std::cout <<
"electron pt = " << electron_iter->pt() <<
" eta : " << electron_iter->eta() << std::endl;
221 double electron_et = electron_iter->et();
222 double electron_phi = electron_iter->phi();
224 double sigma_et, sigma_phi;
235 std::cout <<
"electron sigma_et : " << sigma_et <<
", electron sigma_phi : " << sigma_phi << std::endl;
238 if (sigma_et < 0 || sigma_phi < 0)
239 edm::LogWarning(
"PATMHTProducer") <<
" uncertainties for " << objectname <<
" are (et, phi): " << sigma_et <<
","
240 << sigma_phi <<
" (et,phi): " << electron_et <<
"," << electron_phi;
247 metsig::SigInputObj tmp_electron(objectname, electron_et, electron_phi, sigma_et, sigma_phi);
249 number_of_electrons_++;
274 return number_of_electrons_;
References Abs(), gather_cfg::cout, pwdgSkimBPark_cfi::electrons, iEvent, and AlCaHLTBitMon_QueryRunRegistry::string.
◆ getJets()
Definition at line 117 of file PATMHTProducer.cc.
120 double number_of_jets_ = 0.0;
132 double jet_et = jet_iter->et();
133 double jet_phi = jet_iter->phi();
136 std::cout <<
"jet pt : " << jet_iter->pt() <<
" eta : " << jet_iter->eta()
137 <<
" EMF: " << jet_iter->emEnergyFraction() << std::endl;
140 double sigma_et, sigma_phi;
151 std::cout <<
"jet sigma_et : " << sigma_et <<
", jet sigma_phi : " << sigma_phi << std::endl;
154 if (sigma_et <= 0 || sigma_phi <= 0)
155 edm::LogWarning(
"PATMHTProducer") <<
" uncertainties for " << objectname <<
" are (et, phi): " << sigma_et <<
","
156 << sigma_phi <<
" (et,phi): " << jet_et <<
"," << jet_phi;
169 std::cout <<
"v101> " << number_of_jets_ <<
" " << jet_et <<
" " << sigma_et <<
" " << jet_phi <<
" "
170 << sigma_phi << std::endl;
178 std::vector<CaloTowerPtr> v_towers = jet_iter->getCaloConstituents();
181 for (
unsigned int ii = 0;
ii < v_towers.size();
ii++) {
188 std::cout <<
"v101> --------------------------------------------" << std::endl;
191 return number_of_jets_;
References Abs(), gather_cfg::cout, iEvent, cuy::ii, singleTopDQM_cfi::jets, and AlCaHLTBitMon_QueryRunRegistry::string.
◆ getMuons()
Definition at line 281 of file PATMHTProducer.cc.
288 std::cout <<
">>> PATMHTSelector not valid muon Handle!" << std::endl;
292 double number_of_muons_ = 0.0;
299 std::cout <<
"muon pt = " << muon_iter->pt() <<
" eta : " << muon_iter->eta() << std::endl;
302 double muon_pt = muon_iter->pt();
303 double muon_phi = muon_iter->phi();
305 double sigma_et, sigma_phi;
316 std::cout <<
"muon sigma_et : " << sigma_et <<
", muon sigma_phi : " << sigma_phi << std::endl;
319 if (sigma_et < 0 || sigma_phi < 0)
320 edm::LogWarning(
"PATMHTProducer") <<
" uncertainties for " << objectname <<
" are (et, phi): " << sigma_et <<
","
321 << sigma_phi <<
" (pt,phi): " << muon_pt <<
"," << muon_phi;
334 return number_of_muons_;
References Abs(), gather_cfg::cout, iEvent, edm::HandleBase::isValid(), PDWG_BPHSkim_cff::muons, and AlCaHLTBitMon_QueryRunRegistry::string.
◆ getTowers()
◆ produce()
Definition at line 56 of file PATMHTProducer.cc.
72 std::cout <<
">>>---> Number of jets: " << number_of_jets << std::endl;
73 std::cout <<
">>>---> Number of electrons: " << number_of_jets << std::endl;
74 std::cout <<
">>>---> Number of muons: " << number_of_muons << std::endl;
83 auto themetsigcoll = std::make_unique<pat::MHTCollection>();
93 met_x = met_et *
cos(met_phi);
94 met_y = met_et *
sin(met_phi);
103 themetsigobj.setNumberOfJets(number_of_jets);
104 themetsigobj.setNumberOfElectrons(number_of_electrons);
105 themetsigobj.setNumberOfMuons(number_of_muons);
107 themetsigcoll->push_back(themetsigobj);
References metsig::significanceAlgo::addObjects(), funct::cos(), gather_cfg::cout, iEvent, eostools::move(), pat::MHT::setNumberOfElectrons(), pat::MHT::setNumberOfJets(), pat::MHT::setNumberOfMuons(), met_cff::significance, metsig::significanceAlgo::significance(), and funct::sin().
◆ setUncertaintyParameters()
void PATMHTProducer::setUncertaintyParameters |
( |
| ) |
|
|
private |
Definition at line 338 of file PATMHTProducer.cc.
346 std::make_unique<TF1>(
"ecalEBEtFunc",
"x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
356 std::make_unique<TF1>(
"ecalEEEtFunc",
"x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
368 std::make_unique<TF1>(
"hcalHBEtFunc",
"x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
378 std::make_unique<TF1>(
"hcalHEEtFunc",
"x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
388 std::make_unique<TF1>(
"hcalHOEtFunc",
"x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
398 std::make_unique<TF1>(
"hcalHFEtFunc",
"x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
420 std::make_unique<TF1>(
"jetPhiFunc",
"x*sqrt(([0]*[0]/(x*x))+([1]*[1]/x)+([2]*[2]))", 3);
Referenced by PATMHTProducer().
◆ CaloJetAlgorithmTag_
◆ CaloTowerTag_
◆ controlledUncertainty_
bool pat::PATMHTProducer::controlledUncertainty_ |
|
private |
◆ CorJetAlgorithmTag_
◆ ecalEBUncertainty
◆ ecalEEUncertainty
◆ ElectronTag_
◆ eleEtaMax_
double pat::PATMHTProducer::eleEtaMax_ |
|
private |
◆ eleEtUncertaintyParameter0_
double pat::PATMHTProducer::eleEtUncertaintyParameter0_ |
|
private |
◆ elePhiUncertaintyParameter0_
double pat::PATMHTProducer::elePhiUncertaintyParameter0_ |
|
private |
◆ elePtMin_
double pat::PATMHTProducer::elePtMin_ |
|
private |
◆ eleToken_
◆ eleUncertainty
◆ hcalHBUncertainty
◆ hcalHEUncertainty
◆ hcalHFUncertainty
◆ hcalHOUncertainty
◆ JetCorrectionService_
std::string pat::PATMHTProducer::JetCorrectionService_ |
|
private |
◆ jetCorrUncertainty
◆ jetEMfracMax_
double pat::PATMHTProducer::jetEMfracMax_ |
|
private |
◆ jetEtaMax_
double pat::PATMHTProducer::jetEtaMax_ |
|
private |
◆ jetEtUncertaintyParameter0_
double pat::PATMHTProducer::jetEtUncertaintyParameter0_ |
|
private |
◆ jetEtUncertaintyParameter1_
double pat::PATMHTProducer::jetEtUncertaintyParameter1_ |
|
private |
◆ jetEtUncertaintyParameter2_
double pat::PATMHTProducer::jetEtUncertaintyParameter2_ |
|
private |
◆ jetPhiUncertaintyParameter0_
double pat::PATMHTProducer::jetPhiUncertaintyParameter0_ |
|
private |
◆ jetPhiUncertaintyParameter1_
double pat::PATMHTProducer::jetPhiUncertaintyParameter1_ |
|
private |
◆ jetPhiUncertaintyParameter2_
double pat::PATMHTProducer::jetPhiUncertaintyParameter2_ |
|
private |
◆ jetPtMin_
double pat::PATMHTProducer::jetPtMin_ |
|
private |
◆ jetToken_
◆ jetUncertainty
◆ metCollectionLabel_
std::string pat::PATMHTProducer::metCollectionLabel_ |
|
private |
◆ mhtLabel_
◆ muonChiSqMax_
double pat::PATMHTProducer::muonChiSqMax_ |
|
private |
◆ muonCorrUncertainty
◆ muonDPtMax_
double pat::PATMHTProducer::muonDPtMax_ |
|
private |
◆ muonEtaMax_
double pat::PATMHTProducer::muonEtaMax_ |
|
private |
◆ muonEtUncertaintyParameter0_
double pat::PATMHTProducer::muonEtUncertaintyParameter0_ |
|
private |
◆ muonNHitsMin_
int pat::PATMHTProducer::muonNHitsMin_ |
|
private |
◆ muonPhiUncertaintyParameter0_
double pat::PATMHTProducer::muonPhiUncertaintyParameter0_ |
|
private |
◆ muonPtMin_
double pat::PATMHTProducer::muonPtMin_ |
|
private |
◆ MuonTag_
◆ muonTrackD0Max_
double pat::PATMHTProducer::muonTrackD0Max_ |
|
private |
◆ muonTrackDzMax_
double pat::PATMHTProducer::muonTrackDzMax_ |
|
private |
◆ muonUncertainty
◆ muoToken_
◆ noHF_
bool pat::PATMHTProducer::noHF_ |
|
private |
◆ phoToken_
◆ physobjvector_
◆ s_clusteredTowers
◆ significanceLabel_
std::string pat::PATMHTProducer::significanceLabel_ |
|
private |
◆ tauToken_
◆ towerEtThreshold_
double pat::PATMHTProducer::towerEtThreshold_ |
|
private |
◆ uncertaintyScaleFactor_
double pat::PATMHTProducer::uncertaintyScaleFactor_ |
|
private |
◆ useCaloTowers_
bool pat::PATMHTProducer::useCaloTowers_ |
|
private |
◆ useElectrons_
bool pat::PATMHTProducer::useElectrons_ |
|
private |
◆ useHO_
bool pat::PATMHTProducer::useHO_ |
|
private |
◆ useJets_
bool pat::PATMHTProducer::useJets_ |
|
private |
◆ useMuons_
bool pat::PATMHTProducer::useMuons_ |
|
private |
◆ verbose_
double pat::PATMHTProducer::verbose_ |
|
private |
uncertaintyFunctions hcalHFUncertainty
edm::EDGetTokenT< edm::View< pat::Electron > > eleToken_
bool controlledUncertainty_
double getMuons(edm::Event &, const edm::EventSetup &)
edm::EDGetTokenT< edm::View< pat::Jet > > jetToken_
uncertaintyFunctions muonUncertainty
double jetEtUncertaintyParameter0_
uncertaintyFunctions jetUncertainty
edm::EDGetTokenT< edm::View< pat::Photon > > phoToken_
std::vector< metsig::SigInputObj > physobjvector_
double jetPhiUncertaintyParameter0_
T getUntrackedParameter(std::string const &, T const &) const
uncertaintyFunctions ecalEEUncertainty
Log< level::Warning, false > LogWarning
double jetEtUncertaintyParameter1_
double elePhiUncertaintyParameter0_
math::XYZTLorentzVector LorentzVector
Lorentz vector.
Sin< T >::type sin(const T &t)
const void addObjects(const std::vector< metsig::SigInputObj > &EventVec)
Cos< T >::type cos(const T &t)
std::unique_ptr< TF1 > phiUncertainty
double jetPhiUncertaintyParameter2_
uncertaintyFunctions eleUncertainty
double getElectrons(edm::Event &, const edm::EventSetup &)
edm::EDGetTokenT< edm::View< pat::Muon > > muoToken_
double getJets(edm::Event &, const edm::EventSetup &)
uncertaintyFunctions hcalHEUncertainty
edm::EDGetTokenT< edm::View< pat::Tau > > tauToken_
double eleEtUncertaintyParameter0_
const double significance(double &met_r, double &met_phi, double &met_set)
std::set< CaloTowerDetId > s_clusteredTowers
edm::InputTag CaloTowerTag_
std::unique_ptr< TF1 > etUncertainty
uncertaintyFunctions hcalHBUncertainty
double jetEtUncertaintyParameter2_
uncertaintyFunctions ecalEBUncertainty
uncertaintyFunctions hcalHOUncertainty
T getParameter(std::string const &) const
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
void setUncertaintyParameters()
double jetPhiUncertaintyParameter1_
double muonPhiUncertaintyParameter0_
double muonEtUncertaintyParameter0_
double uncertaintyScaleFactor_