CMS 3D CMS Logo
List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
metsig::SignAlgoResolutions Class Reference

#include <SignAlgoResolutions.h>

Public Member Functions

void addResolutions (const edm::ParameterSet &iConfig)
 
double eval (const resolutionType &type, const resolutionFunc &func, const double &et, const double &phi, const double &eta, const double &p) const
 
double eval (const resolutionType &type, const resolutionFunc &func, const double &et, const double &phi, const double &eta) const
 
metsig::SigInputObj evalPF (const reco::PFCandidate *candidate) const
 
metsig::SigInputObj evalPFJet (const reco::Jet *jet) const
 
bool isFilled () const
 
 SignAlgoResolutions ()
 
 SignAlgoResolutions (const edm::ParameterSet &iConfig)
 

Private Types

typedef std::pair< metsig::resolutionType, metsig::resolutionFunc > functionCombo
 
typedef std::vector< double > functionPars
 

Private Member Functions

void addfunction (const resolutionType type, const resolutionFunc func, const std::vector< double > &parameters)
 
double ElectronPtResolution (const reco::PFCandidate *c) const
 
double EtFunction (const functionPars &x, const functionPars &par) const
 
double getfunc (const resolutionType &type, const resolutionFunc &func, std::vector< double > &x) const
 
void initializeJetResolutions (const edm::ParameterSet &iConfig)
 
double PFunction (const functionPars &x, const functionPars &par) const
 
double PhiConstFunction (const functionPars &x, const functionPars &par) const
 
double PhiFunction (const functionPars &x, const functionPars &par) const
 

Private Attributes

std::map< functionCombo, functionPars > functionmap_
 
std::vector< double > jdphi [10]
 
std::vector< double > jdpt [10]
 
PFEnergyResolution * pfresol_
 
JetResolution * phiResol_
 
JetResolution * ptResol_
 
double ptResolThreshold_
 

Detailed Description

Definition at line 61 of file SignAlgoResolutions.h.

Member Typedef Documentation

◆ functionCombo

typedef std::pair<metsig::resolutionType, metsig::resolutionFunc> metsig::SignAlgoResolutions::functionCombo
private

Definition at line 87 of file SignAlgoResolutions.h.

◆ functionPars

typedef std::vector<double> metsig::SignAlgoResolutions::functionPars
private

Definition at line 88 of file SignAlgoResolutions.h.

Constructor & Destructor Documentation

◆ SignAlgoResolutions() [1/2]

metsig::SignAlgoResolutions::SignAlgoResolutions ( )
inline

Definition at line 63 of file SignAlgoResolutions.h.

63 : functionmap_() { ; }
metsig::SignAlgoResolutions::functionmap_
std::map< functionCombo, functionPars > functionmap_
Definition: SignAlgoResolutions.h:89

◆ SignAlgoResolutions() [2/2]

metsig::SignAlgoResolutions::SignAlgoResolutions ( const edm::ParameterSet &  iConfig)

Definition at line 34 of file SignAlgoResolutions.cc.

References addResolutions().

35  : functionmap_(), ptResol_(nullptr), phiResol_(nullptr) {
36  addResolutions(iConfig);
37 }
metsig::SignAlgoResolutions::functionmap_
std::map< functionCombo, functionPars > functionmap_
Definition: SignAlgoResolutions.h:89
metsig::SignAlgoResolutions::addResolutions
void addResolutions(const edm::ParameterSet &iConfig)
Definition: SignAlgoResolutions.cc:152

Member Function Documentation

◆ addfunction()

void metsig::SignAlgoResolutions::addfunction ( const resolutionType  type,
const resolutionFunc  func,
const std::vector< double > &  parameters 
)
private

Definition at line 311 of file SignAlgoResolutions.cc.

References EcalMonitorTask_cff::func.

313  {
314  functionCombo mypair(type, func);
315  functionmap_[mypair] = parameters;
316 }
EcalMonitorTask_cff.func
func
Definition: EcalMonitorTask_cff.py:10
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
metsig::SignAlgoResolutions::functionmap_
std::map< functionCombo, functionPars > functionmap_
Definition: SignAlgoResolutions.h:89
metsig::SignAlgoResolutions::functionCombo
std::pair< metsig::resolutionType, metsig::resolutionFunc > functionCombo
Definition: SignAlgoResolutions.h:87

◆ addResolutions()

void metsig::SignAlgoResolutions::addResolutions ( const edm::ParameterSet &  iConfig)

Definition at line 152 of file SignAlgoResolutions.cc.

References metsig::caloEB, metsig::caloEE, metsig::caloHB, metsig::caloHE, metsig::caloHF, metsig::caloHO, metsig::ET, edm::ParameterSet::getParameter(), hcalRecHitTable_cff::ieta, initializeJetResolutions(), metsig::PFtype1, metsig::PFtype2, metsig::PFtype3, metsig::PFtype4, metsig::PFtype5, metsig::PFtype6, metsig::PFtype7, and metsig::PHI.

Referenced by SignAlgoResolutions().

152  {
153  using namespace std;
154 
155  // Jet Resolutions - for now load from the files. Migrate to EventSetup asap.
156  metsig::SignAlgoResolutions::initializeJetResolutions(iConfig);
157 
158  ptResolThreshold_ = iConfig.getParameter<double>("ptresolthreshold");
159 
160  //get temporary low pT pfjet resolutions
161  for (int ieta = 0; ieta < 10; ieta++) {
162  jdpt[ieta] = iConfig.getParameter<std::vector<double> >(Form("jdpt%d", ieta));
163  jdphi[ieta] = iConfig.getParameter<std::vector<double> >(Form("jdphi%d", ieta));
164  }
165 
166  // for now: do this by hand - this can obviously also be done via ESSource etc.
167  functionPars etparameters(3, 0);
168  functionPars phiparameters(1, 0);
169  // set the parameters per function:
170  // ECAL, BARREL:
171  std::vector<double> ebet = iConfig.getParameter<std::vector<double> >("EB_EtResPar");
172  std::vector<double> ebphi = iConfig.getParameter<std::vector<double> >("EB_PhiResPar");
173 
174  etparameters[0] = ebet[0];
175  etparameters[1] = ebet[1];
176  etparameters[2] = ebet[2];
177  phiparameters[0] = ebphi[0];
178  addfunction(caloEB, ET, etparameters);
179  addfunction(caloEB, PHI, phiparameters);
180  // ECAL, ENDCAP:
181  std::vector<double> eeet = iConfig.getParameter<std::vector<double> >("EE_EtResPar");
182  std::vector<double> eephi = iConfig.getParameter<std::vector<double> >("EE_PhiResPar");
183 
184  etparameters[0] = eeet[0];
185  etparameters[1] = eeet[1];
186  etparameters[2] = eeet[2];
187  phiparameters[0] = eephi[0];
188  addfunction(caloEE, ET, etparameters);
189  addfunction(caloEE, PHI, phiparameters);
190  // HCAL, BARREL:
191  std::vector<double> hbet = iConfig.getParameter<std::vector<double> >("HB_EtResPar");
192  std::vector<double> hbphi = iConfig.getParameter<std::vector<double> >("HB_PhiResPar");
193 
194  etparameters[0] = hbet[0];
195  etparameters[1] = hbet[1];
196  etparameters[2] = hbet[2];
197  phiparameters[0] = hbphi[0];
198  addfunction(caloHB, ET, etparameters);
199  addfunction(caloHB, PHI, phiparameters);
200  // HCAL, ENDCAP:
201  std::vector<double> heet = iConfig.getParameter<std::vector<double> >("HE_EtResPar");
202  std::vector<double> hephi = iConfig.getParameter<std::vector<double> >("HE_PhiResPar");
203 
204  etparameters[0] = heet[0];
205  etparameters[1] = heet[1];
206  etparameters[2] = heet[2];
207  phiparameters[0] = hephi[0];
208  addfunction(caloHE, ET, etparameters);
209  addfunction(caloHE, PHI, phiparameters);
210  // HCAL, Outer
211  std::vector<double> hoet = iConfig.getParameter<std::vector<double> >("HO_EtResPar");
212  std::vector<double> hophi = iConfig.getParameter<std::vector<double> >("HO_PhiResPar");
213 
214  etparameters[0] = hoet[0];
215  etparameters[1] = hoet[1];
216  etparameters[2] = hoet[2];
217  phiparameters[0] = hophi[0];
218  addfunction(caloHO, ET, etparameters);
219  addfunction(caloHO, PHI, phiparameters);
220  // HCAL, Forward
221  std::vector<double> hfet = iConfig.getParameter<std::vector<double> >("HF_EtResPar");
222  std::vector<double> hfphi = iConfig.getParameter<std::vector<double> >("HF_PhiResPar");
223 
224  etparameters[0] = hfet[0];
225  etparameters[1] = hfet[1];
226  etparameters[2] = hfet[2];
227  phiparameters[0] = hfphi[0];
228  addfunction(caloHF, ET, etparameters);
229  addfunction(caloHF, PHI, phiparameters);
230 
231  // PF objects:
232  // type 1:
233  std::vector<double> pf1et = iConfig.getParameter<std::vector<double> >("PF_EtResType1");
234  std::vector<double> pf1phi = iConfig.getParameter<std::vector<double> >("PF_PhiResType1");
235  etparameters[0] = pf1et[0];
236  etparameters[1] = pf1et[1];
237  etparameters[2] = pf1et[2];
238  phiparameters[0] = pf1phi[0];
239  addfunction(PFtype1, ET, etparameters);
240  addfunction(PFtype1, PHI, phiparameters);
241 
242  // PF objects:
243  // type 2:
244  std::vector<double> pf2et = iConfig.getParameter<std::vector<double> >("PF_EtResType2");
245  std::vector<double> pf2phi = iConfig.getParameter<std::vector<double> >("PF_PhiResType2");
246  etparameters[0] = pf2et[0];
247  etparameters[1] = pf2et[1];
248  etparameters[2] = pf2et[2];
249  phiparameters[0] = pf2phi[0];
250  addfunction(PFtype2, ET, etparameters);
251  addfunction(PFtype2, PHI, phiparameters);
252 
253  // PF objects:
254  // type 3:
255  std::vector<double> pf3et = iConfig.getParameter<std::vector<double> >("PF_EtResType3");
256  std::vector<double> pf3phi = iConfig.getParameter<std::vector<double> >("PF_PhiResType3");
257  etparameters[0] = pf3et[0];
258  etparameters[1] = pf3et[1];
259  etparameters[2] = pf3et[2];
260  phiparameters[0] = pf3phi[0];
261  addfunction(PFtype3, ET, etparameters);
262  addfunction(PFtype3, PHI, phiparameters);
263 
264  // PF objects:
265  // type 4:
266  std::vector<double> pf4et = iConfig.getParameter<std::vector<double> >("PF_EtResType4");
267  std::vector<double> pf4phi = iConfig.getParameter<std::vector<double> >("PF_PhiResType4");
268  etparameters[0] = pf4et[0];
269  etparameters[1] = pf4et[1];
270  etparameters[2] = pf4et[2];
271  //phiparameters[0]=pf4phi[0];
272  addfunction(PFtype4, ET, etparameters);
273  addfunction(PFtype4, PHI, pf4phi); //use the same functional form for photon phi error as for pT, pass whole vector
274 
275  // PF objects:
276  // type 5:
277  std::vector<double> pf5et = iConfig.getParameter<std::vector<double> >("PF_EtResType5");
278  std::vector<double> pf5phi = iConfig.getParameter<std::vector<double> >("PF_PhiResType5");
279  etparameters[0] = pf5et[0];
280  etparameters[1] = pf5et[1];
281  etparameters[2] = pf5et[2];
282  phiparameters[0] = pf5phi[0];
283  addfunction(PFtype5, ET, etparameters);
284  addfunction(PFtype5, PHI, pf5phi);
285 
286  // PF objects:
287  // type 6:
288  std::vector<double> pf6et = iConfig.getParameter<std::vector<double> >("PF_EtResType6");
289  std::vector<double> pf6phi = iConfig.getParameter<std::vector<double> >("PF_PhiResType6");
290  etparameters[0] = pf6et[0];
291  etparameters[1] = pf6et[1];
292  etparameters[2] = pf6et[2];
293  phiparameters[0] = pf6phi[0];
294  addfunction(PFtype6, ET, etparameters);
295  addfunction(PFtype6, PHI, phiparameters);
296 
297  // PF objects:
298  // type 7:
299  std::vector<double> pf7et = iConfig.getParameter<std::vector<double> >("PF_EtResType7");
300  std::vector<double> pf7phi = iConfig.getParameter<std::vector<double> >("PF_PhiResType7");
301  etparameters[0] = pf7et[0];
302  etparameters[1] = pf7et[1];
303  etparameters[2] = pf7et[2];
304  phiparameters[0] = pf7phi[0];
305  addfunction(PFtype7, ET, etparameters);
306  addfunction(PFtype7, PHI, phiparameters);
307 
308  return;
309 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
metsig::SignAlgoResolutions::addfunction
void addfunction(const resolutionType type, const resolutionFunc func, const std::vector< double > &parameters)
Definition: SignAlgoResolutions.cc:311
metsig::SignAlgoResolutions::jdpt
std::vector< double > jdpt[10]
Definition: SignAlgoResolutions.h:101
hcalRecHitTable_cff.ieta
ieta
Definition: hcalRecHitTable_cff.py:15
metsig::SignAlgoResolutions::functionPars
std::vector< double > functionPars
Definition: SignAlgoResolutions.h:88
metsig::SignAlgoResolutions::initializeJetResolutions
void initializeJetResolutions(const edm::ParameterSet &iConfig)
Definition: SignAlgoResolutions.cc:381
metsig::SignAlgoResolutions::jdphi
std::vector< double > jdphi[10]
Definition: SignAlgoResolutions.h:102

◆ ElectronPtResolution()

double metsig::SignAlgoResolutions::ElectronPtResolution ( const reco::PFCandidate *  c) const
private

Definition at line 399 of file SignAlgoResolutions.cc.

References DummyCfis::c, HBHEDarkening_cff::energy, and PVValHelper::eta.

399  {
400  double eta = c->eta();
401  double energy = c->energy();
402  double dEnergy = pfresol_->getEnergyResolutionEm(energy, eta);
403 
404  return dEnergy / cosh(eta);
405 }
PFEnergyResolution::getEnergyResolutionEm
double getEnergyResolutionEm(double CorrectedEnergy, double eta) const
Definition: PFEnergyResolution.cc:18
HBHEDarkening_cff.energy
energy
Definition: HBHEDarkening_cff.py:10

◆ EtFunction()

double metsig::SignAlgoResolutions::EtFunction ( const functionPars &  x,
const functionPars &  par 
) const
private

Definition at line 345 of file SignAlgoResolutions.cc.

References l1tnanotables_cff::et, mps_fire::result, mathSSE::sqrt(), and x.

345  {
346  if (par.size() < 3)
347  return 0.;
348  if (x.empty())
349  return 0.;
350  double et = x[0];
351  if (et <= 0.)
352  return 0.;
353  double result = et * sqrt((par[2] * par[2]) + (par[1] * par[1] / et) + (par[0] * par[0] / (et * et)));
354  return result;
355 }
mps_fire.result
result
Definition: mps_fire.py:311
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:23
x
float x
Definition: beamSpotDipStandalone.cc:55
l1tnanotables_cff.et
et
Definition: l1tnanotables_cff.py:32

◆ eval() [1/2]

double metsig::SignAlgoResolutions::eval ( const resolutionType &  type,
const resolutionFunc &  func,
const double &  et,
const double &  phi,
const double &  eta,
const double &  p 
) const

Definition at line 50 of file SignAlgoResolutions.cc.

References l1tnanotables_cff::et, PVValHelper::eta, EcalMonitorTask_cff::func, class-composition::getfunc, AlCaHLTBitMon_ParallelJobs::p, and x.

Referenced by PFMEtSignInterfaceBase::compResolution().

55  {
56  functionPars x(4);
57  x[0] = et;
58  x[1] = phi;
59  x[2] = eta;
60  x[3] = p;
61  // std::cout << "getting function of type " << type << " " << func << " " << x[0] << " " << x[1] << " " << x[2] << " " << x[3] << std::endl;
62  return getfunc(type, func, x);
63 }
EcalMonitorTask_cff.func
func
Definition: EcalMonitorTask_cff.py:10
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
metsig::SignAlgoResolutions::functionPars
std::vector< double > functionPars
Definition: SignAlgoResolutions.h:88
x
float x
Definition: beamSpotDipStandalone.cc:55
l1tnanotables_cff.et
et
Definition: l1tnanotables_cff.py:32
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
metsig::SignAlgoResolutions::getfunc
double getfunc(const resolutionType &type, const resolutionFunc &func, std::vector< double > &x) const
Definition: SignAlgoResolutions.cc:318

◆ eval() [2/2]

double metsig::SignAlgoResolutions::eval ( const resolutionType &  type,
const resolutionFunc &  func,
const double &  et,
const double &  phi,
const double &  eta 
) const

Definition at line 39 of file SignAlgoResolutions.cc.

References l1tnanotables_cff::et, PVValHelper::eta, JetChargeProducer_cfi::exp, EcalMonitorTask_cff::func, AlCaHLTBitMon_ParallelJobs::p, funct::sin(), and tauSpinnerTable_cfi::theta.

43  {
44  // derive p from et and eta;
45  double theta = 2 * atan(exp(-eta));
46  double p = et / sin(theta); // rough assumption: take e and p equivalent...
47  return eval(type, func, et, phi, eta, p);
48 }
EcalMonitorTask_cff.func
func
Definition: EcalMonitorTask_cff.py:10
tauSpinnerTable_cfi.theta
theta
Definition: tauSpinnerTable_cfi.py:7
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
metsig::SignAlgoResolutions::eval
double eval(const resolutionType &type, const resolutionFunc &func, const double &et, const double &phi, const double &eta, const double &p) const
Definition: SignAlgoResolutions.cc:50
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
l1tnanotables_cff.et
et
Definition: l1tnanotables_cff.py:32
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153

◆ evalPF()

metsig::SigInputObj metsig::SignAlgoResolutions::evalPF ( const reco::PFCandidate *  candidate) const

Definition at line 65 of file SignAlgoResolutions.cc.

References reco::LeafCandidate::energy(), l1tnanotables_cff::et, metsig::ET, PVValHelper::eta, reco::LeafCandidate::eta(), edm::Ref< C, T, F >::isNull(), mergeVDriftHistosByStation::name, reco::PFCandidate::particleId(), metsig::PFtype1, metsig::PFtype2, metsig::PFtype3, metsig::PFtype4, metsig::PFtype5, metsig::PFtype6, metsig::PFtype7, metsig::PHI, reco::LeafCandidate::phi(), funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, reco::LeafCandidate::theta(), and reco::PFCandidate::trackRef().

Referenced by PFMEtSignInterfaceBase::compResolution().

65  {
66  double eta = candidate->eta();
67  double phi = candidate->phi();
68  double et = candidate->energy() * sin(candidate->theta());
69  resolutionType thetype;
70  std::string name;
71  int type = candidate->particleId();
72  switch (type) {
73  case 1:
74  thetype = PFtype1;
75  name = "PFChargedHadron";
76  break;
77  case 2:
78  thetype = PFtype2;
79  name = "PFChargedEM";
80  break;
81  case 3:
82  thetype = PFtype3;
83  name = "PFMuon";
84  break;
85  case 4:
86  thetype = PFtype4;
87  name = "PFNeutralEM";
88  break;
89  case 5:
90  thetype = PFtype5;
91  name = "PFNeutralHadron";
92  break;
93  case 6:
94  thetype = PFtype6;
95  name = "PFtype6";
96  break;
97  case 7:
98  thetype = PFtype7;
99  name = "PFtype7";
100  break;
101  default:
102  thetype = PFtype7;
103  name = "PFunknown";
104  break;
105  }
106 
107  double d_et = 0, d_phi = 0; //d_phi here is the error on phi component of the et
108  reco::TrackRef trackRef = candidate->trackRef();
109  if (!trackRef.isNull()) {
110  d_phi = et * trackRef->phiError();
111  d_et = (type == 2) ? ElectronPtResolution(candidate) : trackRef->ptError();
112  //if(type==2) std::cout << eval(thetype,ET,et,phi,eta) << " " << trackRef->ptError() << " "<< ElectronPtResolution(candidate) << std::endl;
113  } else {
114  d_et = eval(thetype, ET, et, phi, eta);
115  d_phi = eval(thetype, PHI, et, phi, eta);
116  }
117 
118  metsig::SigInputObj resultingobj(name, et, phi, d_et, d_phi);
119  return resultingobj;
120 }
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
reco::LeafCandidate::theta
double theta() const final
momentum polar angle
Definition: LeafCandidate.h:150
edm::Ref::isNull
bool isNull() const
Checks for null.
Definition: Ref.h:229
metsig::SignAlgoResolutions::eval
double eval(const resolutionType &type, const resolutionFunc &func, const double &et, const double &phi, const double &eta, const double &p) const
Definition: SignAlgoResolutions.cc:50
metsig::SignAlgoResolutions::ElectronPtResolution
double ElectronPtResolution(const reco::PFCandidate *c) const
Definition: SignAlgoResolutions.cc:399
reco::PFCandidate::trackRef
reco::TrackRef trackRef() const
Definition: PFCandidate.cc:437
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:148
l1tnanotables_cff.et
et
Definition: l1tnanotables_cff.py:32
reco::LeafCandidate::energy
double energy() const final
energy
Definition: LeafCandidate.h:125
reco::PFCandidate::particleId
virtual ParticleType particleId() const
Definition: PFCandidate.h:392
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:152

◆ evalPFJet()

metsig::SigInputObj metsig::SignAlgoResolutions::evalPFJet ( const reco::Jet *  jet) const

Definition at line 122 of file SignAlgoResolutions.cc.

References JetMETHLTOfflineSource_cfi::feta, hcalRecHitTable_cff::ieta, createfilelist::int, metsig::jet, METSignificanceParams_cfi::jeta, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by PFMEtSignInterfaceBase::compResolution().

122  {
123  double jpt = jet->pt();
124  double jphi = jet->phi();
125  double jeta = jet->eta();
126  double jdeltapt = 999.;
127  double jdeltapphi = 999.;
128 
129  if (jpt < ptResolThreshold_ && jpt < 20.) { //use temporary fix for low pT jets
130  double feta = TMath::Abs(jeta);
131  int ieta = feta < 5. ? int(feta / 0.5) : 9; //bin size = 0.5
132  int ipt = jpt > 3. ? int(jpt - 3. / 2) : 0; //bin size =2, starting from ptmin=3GeV
133  jdeltapt = jdpt[ieta][ipt];
134  jdeltapphi = jpt * jdphi[ieta][ipt];
135  } else {
136  //use the resolution functions at |eta|=5 to avoid crash for jets with large eta.
137  if (jeta > 5)
138  jeta = 5;
139  if (jeta < -5)
140  jeta = -5;
141  double jptForEval = jpt > ptResolThreshold_ ? jpt : ptResolThreshold_;
142  jdeltapt = jpt * ptResol_->parameterEtaEval("sigma", jeta, jptForEval);
143  jdeltapphi = jpt * phiResol_->parameterEtaEval("sigma", jeta, jptForEval);
144  }
145 
146  std::string inputtype = "jet";
147  metsig::SigInputObj obj_jet(inputtype, jpt, jphi, jdeltapt, jdeltapphi);
148  //std::cout << "RESOLUTIONS JET: " << jpt << " " << jphi<< " " <<jdeltapt << " " << jdeltapphi << std::endl;
149  return obj_jet;
150 }
metsig::SignAlgoResolutions::jdpt
std::vector< double > jdpt[10]
Definition: SignAlgoResolutions.h:101
METSignificanceParams_cfi.jeta
jeta
Definition: METSignificanceParams_cfi.py:12
JetResolution::parameterEtaEval
double parameterEtaEval(const std::string &parameterName, float eta, float pt)
Definition: JetResolution.cc:170
JetMETHLTOfflineSource_cfi.feta
feta
Definition: JetMETHLTOfflineSource_cfi.py:30
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
hcalRecHitTable_cff.ieta
ieta
Definition: hcalRecHitTable_cff.py:15
metsig::SignAlgoResolutions::jdphi
std::vector< double > jdphi[10]
Definition: SignAlgoResolutions.h:102

◆ getfunc()

double metsig::SignAlgoResolutions::getfunc ( const resolutionType &  type,
const resolutionFunc &  func,
std::vector< double > &  x 
) const
private

Definition at line 318 of file SignAlgoResolutions.cc.

References metsig::CONSTPHI, metsig::ET, EcalMonitorTask_cff::func, metsig::PHI, mps_fire::result, edm::second(), metsig::TRACKP, contentValuesCheck::values, and x.

320  {
321  double result = 0;
322  functionCombo mypair(type, func);
323 
324  if (functionmap_.count(mypair) == 0) {
325  return result;
326  }
327 
328  functionPars values = (functionmap_.find(mypair))->second;
329  switch (func) {
330  case metsig::ET:
331  return EtFunction(x, values);
332  case metsig::PHI:
333  return PhiFunction(x, values);
334  case metsig::TRACKP:
335  return PFunction(x, values);
336  case metsig::CONSTPHI:
337  return PhiConstFunction(x, values);
338  }
339 
340  // std::cout << "returning function " << type << " " << func << " " << result << " " << x[0] << std::endl;
341 
342  return result;
343 }
EcalMonitorTask_cff.func
func
Definition: EcalMonitorTask_cff.py:10
metsig::SignAlgoResolutions::EtFunction
double EtFunction(const functionPars &x, const functionPars &par) const
Definition: SignAlgoResolutions.cc:345
mps_fire.result
result
Definition: mps_fire.py:311
metsig::SignAlgoResolutions::PFunction
double PFunction(const functionPars &x, const functionPars &par) const
Definition: SignAlgoResolutions.cc:372
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
metsig::SignAlgoResolutions::PhiFunction
double PhiFunction(const functionPars &x, const functionPars &par) const
Definition: SignAlgoResolutions.cc:357
metsig::SignAlgoResolutions::functionmap_
std::map< functionCombo, functionPars > functionmap_
Definition: SignAlgoResolutions.h:89
metsig::SignAlgoResolutions::functionPars
std::vector< double > functionPars
Definition: SignAlgoResolutions.h:88
metsig::SignAlgoResolutions::PhiConstFunction
double PhiConstFunction(const functionPars &x, const functionPars &par) const
Definition: SignAlgoResolutions.cc:377
metsig::SignAlgoResolutions::functionCombo
std::pair< metsig::resolutionType, metsig::resolutionFunc > functionCombo
Definition: SignAlgoResolutions.h:87

◆ initializeJetResolutions()

void metsig::SignAlgoResolutions::initializeJetResolutions ( const edm::ParameterSet &  iConfig)
private

Definition at line 381 of file SignAlgoResolutions.cc.

References edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), METSigParams_cfi::resolutionsAlgo, METSigParams_cfi::resolutionsEra, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by addResolutions().

381  {
382  using namespace std;
383 
384  // only reinitialize the resolutsion if the pointers are zero
385  if (ptResol_ == nullptr) {
386  string resolutionsAlgo = iConfig.getParameter<std::string>("resolutionsAlgo");
387  string resolutionsEra = iConfig.getParameter<std::string>("resolutionsEra");
388  string ptFileName = edm::FileInPath("CondFormats/JetMETObjects/data/" + resolutionsEra + "_PtResolution_" +
389  resolutionsAlgo + ".txt")
390  .fullPath();
391  string phiFileName = edm::FileInPath("CondFormats/JetMETObjects/data/" + resolutionsEra + "_PhiResolution_" +
392  resolutionsAlgo + ".txt")
393  .fullPath();
394  ptResol_ = new JetResolution(ptFileName, false);
395  phiResol_ = new JetResolution(phiFileName, false);
396  }
397 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
METSigParams_cfi.resolutionsEra
resolutionsEra
Definition: METSigParams_cfi.py:50
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
METSigParams_cfi.resolutionsAlgo
resolutionsAlgo
Definition: METSigParams_cfi.py:51
edm::FileInPath::fullPath
const std::string & fullPath() const
Definition: FileInPath.cc:144

◆ isFilled()

bool metsig::SignAlgoResolutions::isFilled ( ) const
inline

Definition at line 80 of file SignAlgoResolutions.h.

References functionmap_.

80 { return !functionmap_.empty(); }
metsig::SignAlgoResolutions::functionmap_
std::map< functionCombo, functionPars > functionmap_
Definition: SignAlgoResolutions.h:89

◆ PFunction()

double metsig::SignAlgoResolutions::PFunction ( const functionPars &  x,
const functionPars &  par 
) const
private

Definition at line 372 of file SignAlgoResolutions.cc.

372  {
373  // not currently implemented
374  return 0;
375 }

◆ PhiConstFunction()

double metsig::SignAlgoResolutions::PhiConstFunction ( const functionPars &  x,
const functionPars &  par 
) const
private

Definition at line 377 of file SignAlgoResolutions.cc.

377  {
378  return par[0];
379 }

◆ PhiFunction()

double metsig::SignAlgoResolutions::PhiFunction ( const functionPars &  x,
const functionPars &  par 
) const
private

Definition at line 357 of file SignAlgoResolutions.cc.

References l1tnanotables_cff::et, mathSSE::sqrt(), and x.

357  {
358  double et = x[0];
359  if (et <= 0.) {
360  return 0.;
361  }
362 
363  //if 1 parameter is C provided, returns C*pT, if three parameters N, S, C are provided, it returns the usual resolution value, as for sigmaPt
364  if (par.size() != 1 && par.size() != 3) { //only 1 or 3 parameters supported for phi function
365  return 0.;
366  } else if (par.size() == 1) {
367  return par[0] * et;
368  } else {
369  return et * sqrt((par[2] * par[2]) + (par[1] * par[1] / et) + (par[0] * par[0] / (et * et)));
370  }
371 }
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:23
x
float x
Definition: beamSpotDipStandalone.cc:55
l1tnanotables_cff.et
et
Definition: l1tnanotables_cff.py:32

Member Data Documentation

◆ functionmap_

std::map<functionCombo, functionPars> metsig::SignAlgoResolutions::functionmap_
private

Definition at line 89 of file SignAlgoResolutions.h.

Referenced by isFilled().

◆ jdphi

std::vector<double> metsig::SignAlgoResolutions::jdphi[10]
private

Definition at line 102 of file SignAlgoResolutions.h.

◆ jdpt

std::vector<double> metsig::SignAlgoResolutions::jdpt[10]
private

Definition at line 101 of file SignAlgoResolutions.h.

◆ pfresol_

PFEnergyResolution* metsig::SignAlgoResolutions::pfresol_
private

Definition at line 106 of file SignAlgoResolutions.h.

◆ phiResol_

JetResolution* metsig::SignAlgoResolutions::phiResol_
private

Definition at line 105 of file SignAlgoResolutions.h.

◆ ptResol_

JetResolution* metsig::SignAlgoResolutions::ptResol_
private

Definition at line 104 of file SignAlgoResolutions.h.

◆ ptResolThreshold_

double metsig::SignAlgoResolutions::ptResolThreshold_
private

Definition at line 97 of file SignAlgoResolutions.h.

Generated for CMSSW Reference Manual by doxygen 1.8.14