CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Functions
pat::helper::ParametrizationHelper Namespace Reference

Functions

AlgebraicVector4 diffToParameters (pat::CandKinResolution::Parametrization parametrization, const math::XYZTLorentzVector &p4ini, const math::XYZTLorentzVector &p4fin)
 Expresses the difference between two 4-momentum vectors as a shift in coordinates in a given frame. More...
 
AlgebraicVector4 diffToParameters (pat::CandKinResolution::Parametrization parametrization, const math::PtEtaPhiMLorentzVector &p4ini, const math::PtEtaPhiMLorentzVector &p4fin)
 Expresses the difference between two 4-momentum vectors as a shift in coordinates in a given frame. More...
 
uint32_t dimension (pat::CandKinResolution::Parametrization parametrization)
 Returns the number of free parameters in a parametrization (3 or 4) More...
 
pat::CandKinResolution::Parametrization fromString (const std::string &name)
 Convert a name into a parametrization code. More...
 
bool isAlwaysMassive (pat::CandKinResolution::Parametrization parametrization)
 Is this parametrization usable only with massive objects? More...
 
bool isAlwaysMassless (pat::CandKinResolution::Parametrization parametrization)
 Is this parametrization usable only with massless objects? More...
 
bool isMassConstrained (pat::CandKinResolution::Parametrization parametrization)
 If this parametrization has a mass constraint (including the 'isAlwaysMassless' case) More...
 
bool isPhysical (pat::CandKinResolution::Parametrization parametrization, const AlgebraicVector4 &v4, const math::PtEtaPhiMLorentzVector &initialP4)
 
const char * name (pat::CandKinResolution::Parametrization param)
 Convert a number into a string. More...
 
math::XYZTLorentzVector p4fromParameters (pat::CandKinResolution::Parametrization parametrization, const AlgebraicVector4 &parameters, const math::XYZTLorentzVector &initialP4)
 
AlgebraicVector4 parametersFromP4 (pat::CandKinResolution::Parametrization parametrization, const math::XYZTLorentzVector &p4)
 Returns a vector of coordinates values given a coordinate frame and a 4-vector. More...
 
AlgebraicVector4 parametersFromP4 (pat::CandKinResolution::Parametrization parametrization, const math::PtEtaPhiMLorentzVector &p4)
 Returns a vector of coordinates values given a coordinate frame and a 4-vector. More...
 
math::PtEtaPhiMLorentzVector polarP4fromParameters (pat::CandKinResolution::Parametrization parametrization, const AlgebraicVector4 &parameters, const math::PtEtaPhiMLorentzVector &initialP4)
 
template<typename T >
void setParametersFromAnyVector (pat::CandKinResolution::Parametrization parametrization, AlgebraicVector4 &pars, const T &p4)
 For internal use only, so we provide only the interface. Use the 'setParametersFromP4'. More...
 
void setParametersFromP4 (pat::CandKinResolution::Parametrization parametrization, AlgebraicVector4 &pars, const math::XYZTLorentzVector &p4)
 Set the values of the parameters for a given 4-momentum. More...
 
void setParametersFromP4 (pat::CandKinResolution::Parametrization parametrization, AlgebraicVector4 &pars, const math::PtEtaPhiMLorentzVector &p4)
 Set the values of the parameters for a given 4-momentum. More...
 

Function Documentation

AlgebraicVector4 pat::helper::ParametrizationHelper::diffToParameters ( pat::CandKinResolution::Parametrization  parametrization,
const math::XYZTLorentzVector &  p4ini,
const math::XYZTLorentzVector &  p4fin 
)

Expresses the difference between two 4-momentum vectors as a shift in coordinates in a given frame.

Basically, if you do: pars = parametersFromP4(param, simp4); diff = diffToParameters(param, simP4, recP4); then up to roundoff errors recP4 == p4fromParameters(param, pars+diff, simP4);

Definition at line 323 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::Invalid, M_PI, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, p1, p2, parametersFromP4(), run_regression::ret, and pat::CandKinResolution::Spher.

Referenced by diffToParameters().

325 {
326  AlgebraicVector4 ret;
327  switch (parametrization) {
328  case pat::CandKinResolution::Cart:
329  case pat::CandKinResolution::ECart:
330  case pat::CandKinResolution::MCCart:
331  case pat::CandKinResolution::Spher:
332  ret = parametersFromP4(parametrization,p4fin) - parametersFromP4(parametrization,p4ini);
333  break;
334  case pat::CandKinResolution::ESpher:
335  case pat::CandKinResolution::MCSpher:
336  case pat::CandKinResolution::MCPInvSpher:
337  case pat::CandKinResolution::EtThetaPhi:
338  case pat::CandKinResolution::EtEtaPhi:
339  ret = parametersFromP4(parametrization,p4fin) - parametersFromP4(parametrization,p4ini);
340  while(ret[2] > +M_PI) ret[2] -= (2*M_PI);
341  while(ret[2] < -M_PI) ret[2] += (2*M_PI);
342  break;
343  case pat::CandKinResolution::MCMomDev:
344  case pat::CandKinResolution::MomDev:
345  case pat::CandKinResolution::EMomDev:
346  case pat::CandKinResolution::EScaledMomDev:
347  {
348  typedef ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D<double> > V3Cart;
349  V3Cart p1 = p4ini.Vect(), p2 = p4fin.Vect();
350  V3Cart ur = p1.Unit();
351  V3Cart uz(0,0,1);
352  V3Cart uph = uz.Cross(ur).Unit();
353  V3Cart uth = ur.Cross(uph).Unit();
354  ret[0] = p2.Dot(ur)/p1.R() - 1.0;
355  ret[1] = (p2 - p1).Dot(uth);
356  ret[2] = (p2 - p1).Dot(uph);
357  if (parametrization == pat::CandKinResolution::MomDev) {
358  ret[3] = p4fin.mass()/p4ini.mass() - 1.0;
359  } else if (parametrization == pat::CandKinResolution::EMomDev) {
360  ret[3] = p4fin.energy()/p4ini.energy() - 1.0;
361  }
362  }
363  break;
364  case pat::CandKinResolution::Invalid:
365  throw cms::Exception("Invalid parametrization") << parametrization;
366  default:
367  throw cms::Exception("Not Implemented") << "diffToParameters not yet implemented for parametrization " << parametrization ;
368  }
369  return ret;
370 }
pat::helper::ParametrizationHelper::parametersFromP4
AlgebraicVector4 parametersFromP4(pat::CandKinResolution::Parametrization parametrization, const math::XYZTLorentzVector &p4)
Returns a vector of coordinates values given a coordinate frame and a 4-vector.
Definition: ParametrizationHelper.cc:279
run_regression.ret
int ret
Definition: run_regression.py:388
p2
double p2[4]
Definition: TauolaWrapper.h:90
M_PI
#define M_PI
Definition: readTestVector.cc:25
p1
double p1[4]
Definition: TauolaWrapper.h:89
AlgebraicVector4
ROOT::Math::SVector< double, 4 > AlgebraicVector4
Definition: AlgebraicROOTObjects.h:21
AlgebraicVector4 pat::helper::ParametrizationHelper::diffToParameters ( pat::CandKinResolution::Parametrization  parametrization,
const math::PtEtaPhiMLorentzVector &  p4ini,
const math::PtEtaPhiMLorentzVector &  p4fin 
)

Expresses the difference between two 4-momentum vectors as a shift in coordinates in a given frame.

Basically, if you do: pars = parametersFromP4(param, simp4); diff = diffToParameters(param, simP4, recP4); then up to roundoff errors recP4 == p4fromParameters(param, pars+diff, simP4);

Definition at line 286 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, diffToParameters(), pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::Invalid, M_PI, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, parametersFromP4(), run_regression::ret, and pat::CandKinResolution::Spher.

288 {
289  AlgebraicVector4 ret;
290  switch (parametrization) {
291  case pat::CandKinResolution::Cart:
292  case pat::CandKinResolution::ECart:
293  case pat::CandKinResolution::MCCart:
294  ret = parametersFromP4(parametrization,p4fin) - parametersFromP4(parametrization,p4ini);
295  break;
296  case pat::CandKinResolution::Spher:
297  case pat::CandKinResolution::ESpher:
298  case pat::CandKinResolution::MCSpher:
299  case pat::CandKinResolution::MCPInvSpher:
300  case pat::CandKinResolution::EtThetaPhi:
301  case pat::CandKinResolution::EtEtaPhi:
302  ret = parametersFromP4(parametrization,p4fin) - parametersFromP4(parametrization,p4ini);
303  while(ret[2] > +M_PI) ret[2] -= (2*M_PI);
304  while(ret[2] < -M_PI) ret[2] += (2*M_PI);
305  break;
306  case pat::CandKinResolution::MCMomDev:
307  case pat::CandKinResolution::MomDev:
308  case pat::CandKinResolution::EMomDev:
309  case pat::CandKinResolution::EScaledMomDev:
310  return diffToParameters(parametrization,
311  math::XYZTLorentzVector(p4ini),
312  math::XYZTLorentzVector(p4fin));
313  case pat::CandKinResolution::Invalid:
314  throw cms::Exception("Invalid parametrization") << parametrization;
315  default:
316  throw cms::Exception("Not Implemented") << "diffToParameters not yet implemented for parametrization " << parametrization ;
317  }
318  return ret;
319 }
pat::helper::ParametrizationHelper::parametersFromP4
AlgebraicVector4 parametersFromP4(pat::CandKinResolution::Parametrization parametrization, const math::XYZTLorentzVector &p4)
Returns a vector of coordinates values given a coordinate frame and a 4-vector.
Definition: ParametrizationHelper.cc:279
run_regression.ret
int ret
Definition: run_regression.py:388
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
M_PI
#define M_PI
Definition: readTestVector.cc:25
pat::helper::ParametrizationHelper::diffToParameters
AlgebraicVector4 diffToParameters(pat::CandKinResolution::Parametrization parametrization, const math::XYZTLorentzVector &p4ini, const math::XYZTLorentzVector &p4fin)
Expresses the difference between two 4-momentum vectors as a shift in coordinates in a given frame...
Definition: ParametrizationHelper.cc:323
AlgebraicVector4
ROOT::Math::SVector< double, 4 > AlgebraicVector4
Definition: AlgebraicROOTObjects.h:21
uint32_t pat::helper::ParametrizationHelper::dimension ( pat::CandKinResolution::Parametrization  parametrization)
inline

Returns the number of free parameters in a parametrization (3 or 4)

Definition at line 11 of file ParametrizationHelper.h.

Referenced by MillePedeAlignmentAlgorithm::callMille2D(), CosmicMuonSeedGenerator::createSeed(), SETSeedFinder::estimateMomentum(), pat::CandKinResolution::fillMatrix(), pat::CandKinResolution::fillVector(), and SETSeedFinder::pre_prune().

11  {
12  return (static_cast<uint32_t>(parametrization) & 0x0F);
13  }
pat::CandKinResolution::Parametrization pat::helper::ParametrizationHelper::fromString ( const std::string &  name)

Convert a name into a parametrization code.

Definition at line 15 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, and pat::CandKinResolution::Spher.

Referenced by StringResolutionProvider::StringResolutionProvider(), and ora::Version::thisSchemaVersion().

15  {
16  using namespace std;
17 
18  typedef pat::CandKinResolution::Parametrization Parametrization;
19 
20  static map<string,Parametrization> parMaps;
21  if (parMaps.empty()) {
22  parMaps["Cart"] = pat::CandKinResolution::Cart;
23  parMaps["ECart"] = pat::CandKinResolution::ECart;
24  parMaps["MCCart"] = pat::CandKinResolution::MCCart;
25  parMaps["Spher"] = pat::CandKinResolution::Spher;
26  parMaps["ESpher"] = pat::CandKinResolution::ESpher;
27  parMaps["MCSpher"] = pat::CandKinResolution::MCSpher;
28  parMaps["MCPInvSpher"] = pat::CandKinResolution::MCPInvSpher;
29  parMaps["EtEtaPhi"] = pat::CandKinResolution::EtEtaPhi;
30  parMaps["EtThetaPhi"] = pat::CandKinResolution::EtThetaPhi;
31  parMaps["MomDev"] = pat::CandKinResolution::MomDev;
32  parMaps["EMomDev"] = pat::CandKinResolution::EMomDev;
33  parMaps["MCMomDev"] = pat::CandKinResolution::MCMomDev;
34  parMaps["EScaledMomDev"] = pat::CandKinResolution::EScaledMomDev;
35  }
36  map<string,Parametrization>::const_iterator itP = parMaps.find(name);
37  if (itP == parMaps.end()) {
38  throw cms::Exception("StringResolutionProvider") << "Bad parametrization '" << name.c_str() << "'";
39  }
40  return itP->second;
41 }
bool pat::helper::ParametrizationHelper::isAlwaysMassive ( pat::CandKinResolution::Parametrization  parametrization)

Is this parametrization usable only with massive objects?

Definition at line 399 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::Invalid, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, and pat::CandKinResolution::Spher.

400 {
401  switch (parametrization) {
402  case pat::CandKinResolution::Cart:
403  case pat::CandKinResolution::ECart:
404  case pat::CandKinResolution::MCCart:
405  case pat::CandKinResolution::Spher:
406  case pat::CandKinResolution::ESpher:
407  case pat::CandKinResolution::MCSpher:
408  case pat::CandKinResolution::MCPInvSpher:
409  case pat::CandKinResolution::MCMomDev:
410  case pat::CandKinResolution::EMomDev:
411  case pat::CandKinResolution::EScaledMomDev:
412  case pat::CandKinResolution::EtThetaPhi:
413  case pat::CandKinResolution::EtEtaPhi:
414  return false;
415  case pat::CandKinResolution::MomDev:
416  return true;
417  case pat::CandKinResolution::Invalid:
418  throw cms::Exception("Invalid parametrization") << parametrization;
419  default:
420  throw cms::Exception("Not Implemented") << "isAlwaysMassless not yet implemented for parametrization " << parametrization ;
421  }
422 }
bool pat::helper::ParametrizationHelper::isAlwaysMassless ( pat::CandKinResolution::Parametrization  parametrization)

Is this parametrization usable only with massless objects?

Definition at line 373 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::Invalid, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, and pat::CandKinResolution::Spher.

374 {
375  switch (parametrization) {
376  case pat::CandKinResolution::Cart:
377  case pat::CandKinResolution::ECart:
378  case pat::CandKinResolution::MCCart:
379  case pat::CandKinResolution::Spher:
380  case pat::CandKinResolution::ESpher:
381  case pat::CandKinResolution::MCSpher:
382  case pat::CandKinResolution::MCPInvSpher:
383  case pat::CandKinResolution::MCMomDev:
384  case pat::CandKinResolution::MomDev:
385  case pat::CandKinResolution::EMomDev:
386  case pat::CandKinResolution::EScaledMomDev:
387  return false;
388  case pat::CandKinResolution::EtThetaPhi:
389  case pat::CandKinResolution::EtEtaPhi:
390  return true;
391  case pat::CandKinResolution::Invalid:
392  throw cms::Exception("Invalid parametrization") << parametrization;
393  default:
394  throw cms::Exception("Not Implemented") << "isAlwaysMassless not yet implemented for parametrization " << parametrization ;
395  }
396 }
bool pat::helper::ParametrizationHelper::isMassConstrained ( pat::CandKinResolution::Parametrization  parametrization)

If this parametrization has a mass constraint (including the 'isAlwaysMassless' case)

Definition at line 425 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::Invalid, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, and pat::CandKinResolution::Spher.

426 {
427  switch (parametrization) {
428  case pat::CandKinResolution::Cart:
429  case pat::CandKinResolution::ECart:
430  case pat::CandKinResolution::Spher:
431  case pat::CandKinResolution::ESpher:
432  case pat::CandKinResolution::EMomDev:
433  case pat::CandKinResolution::EScaledMomDev:
434  case pat::CandKinResolution::EtThetaPhi:
435  case pat::CandKinResolution::EtEtaPhi:
436  case pat::CandKinResolution::MomDev:
437  return false;
438  case pat::CandKinResolution::MCCart:
439  case pat::CandKinResolution::MCSpher:
440  case pat::CandKinResolution::MCPInvSpher:
441  case pat::CandKinResolution::MCMomDev:
442  return true;
443  case pat::CandKinResolution::Invalid:
444  throw cms::Exception("Invalid parametrization") << parametrization;
445  default:
446  throw cms::Exception("Not Implemented") << "isAlwaysMassless not yet implemented for parametrization " << parametrization ;
447  }
448 }
bool pat::helper::ParametrizationHelper::isPhysical ( pat::CandKinResolution::Parametrization  parametrization,
const AlgebraicVector4 &  v4,
const math::PtEtaPhiMLorentzVector &  initialP4 
)

If this value of the parameters is meaningful in this parametrization. It can be used e.g. when doing random smearing to check you're still within the physical region This DOES check inequalities (e.g. E >= P, M >= 0, theta in [0,PI], ..) This DOES NOT check strict equalities (e.g. M == 0) This DOES NOT check that your parameters comply with your constraints (e.g. fixed mass constraint)

Definition at line 451 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::Invalid, M_PI, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, AlCaHLTBitMon_ParallelJobs::p, and pat::CandKinResolution::Spher.

453  {
454  switch (parametrization) {
455  // ======= CARTESIAN ==========
456  case pat::CandKinResolution::Cart:
457  return parameters[3] >= 0; // M >= 0
458  case pat::CandKinResolution::MCCart:
459  return true;
460  case pat::CandKinResolution::ECart:
461  return (parameters[0]*parameters[0] + parameters[1]*parameters[1] + parameters[2]*parameters[2] <= parameters[3]*parameters[3]); // E >= P
462  case pat::CandKinResolution::Spher:
463  return (parameters[0] >= 0 ) && // P >= 0
464  (parameters[3] >= 0 ) && // M >= 0
465  (parameters[1] >= 0 ) && // theta >= 0
466  (parameters[1] <= M_PI); // theta <= pi
467  case pat::CandKinResolution::MCSpher:
468  return (parameters[0] >= 0 ) && // P >= 0
469  (parameters[1] >= 0 ) && // theta >= 0
470  (parameters[1] <= M_PI); // theta <= pi
471  case pat::CandKinResolution::ESpher: //
472  return (parameters[0] >= 0 ) && // P >= 0
473  (parameters[3] >= parameters[0]) && // E >= P
474  (parameters[1] >= 0 ) && // theta >= 0
475  (parameters[1] <= M_PI ) ; // theta <= PI
476  case pat::CandKinResolution::MCPInvSpher: //
477  return (parameters[0] > 0 ) && // 1/P > 0
478  (parameters[1] >= 0 ) && // theta >= 0
479  (parameters[1] <= M_PI); // theta <= pi
480  // ======= HEP CYLINDRICAL ==========
481  case pat::CandKinResolution::EtThetaPhi:
482  return (parameters[0] > 0 ) && // Et >= 0
483  (parameters[1] >= 0 ) && // theta >= 0
484  (parameters[1] <= M_PI); // theta <= pi
485  case pat::CandKinResolution::EtEtaPhi: // as simple as that
486  return (parameters[0] > 0); // Et >= 0
487  // ======= MomentumDeviates ==========
488  case pat::CandKinResolution::MomDev:
489  return (parameters[0] >= 0) && // P >= 0
490  (parameters[3] >= 0); // m/M0 >= 0
491  case pat::CandKinResolution::MCMomDev:
492  return (parameters[0] >= 0); // P >= 0
493  case pat::CandKinResolution::EMomDev:
494  case pat::CandKinResolution::EScaledMomDev:
495  {
496  if (parameters[0] <= 0) return false;
497  ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D<double> > p = initialP4.Vect(), uz(0,0,1), uph, uth;
498  uph = uz.Cross(p).Unit();
499  uth = p.Cross(uph).Unit();
500  p *= parameters[0];
501  p += uth * parameters[1] + uph * parameters[2];
502  if (parametrization == pat::CandKinResolution::EMomDev) {
503  if (parameters[3] < 0) return false;
504  double m2 = ROOT::Math::Square(parameters[3] * initialP4.energy()) - p.Mag2();
505  if (m2 < 0) return false;
506  } else if (parametrization == pat::CandKinResolution::EScaledMomDev) {
507  if (parameters[3] < 0) return false;
508  double m2 = ROOT::Math::Square(p.R()*initialP4.E()/initialP4.P()) - p.Mag2();
509  if (m2 < 0) return false;
510  }
511  return true;
512  }
513  // ======= OTHER ==========
514  case pat::CandKinResolution::Invalid:
515  throw cms::Exception("Invalid parametrization") << parametrization;
516  default:
517  throw cms::Exception("Not Implemented") << "getResolEta not yet implemented for parametrization " << parametrization ;
518  }
519 }
Parameters.parameters
dictionary parameters
Definition: Parameters.py:2
M_PI
#define M_PI
Definition: readTestVector.cc:25
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
const char * pat::helper::ParametrizationHelper::name ( pat::CandKinResolution::Parametrization  param)

Convert a number into a string.

Definition at line 44 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, pat::CandKinResolution::Invalid, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, and pat::CandKinResolution::Spher.

44  {
45  switch (param) {
46  case pat::CandKinResolution::Cart: return "Cart";
47  case pat::CandKinResolution::ECart: return "ECart";
48  case pat::CandKinResolution::MCCart: return "MCCart";
49  case pat::CandKinResolution::Spher: return "Spher";
50  case pat::CandKinResolution::ESpher: return "ESpher";
51  case pat::CandKinResolution::MCSpher: return "MCSpher";
52  case pat::CandKinResolution::MCPInvSpher: return "MCPInvSpher";
53  case pat::CandKinResolution::EtEtaPhi: return "EtEtaPhi";
54  case pat::CandKinResolution::EtThetaPhi: return "EtThetaPhi";
55  case pat::CandKinResolution::MomDev: return "MomDev";
56  case pat::CandKinResolution::EMomDev: return "EMomDev";
57  case pat::CandKinResolution::MCMomDev: return "MCMomDev";
58  case pat::CandKinResolution::EScaledMomDev: return "EScaledMomDev";
59  case pat::CandKinResolution::Invalid: return "Invalid";
60  default: return "UNKNOWN";
61  }
62 }
math::XYZTLorentzVector pat::helper::ParametrizationHelper::p4fromParameters ( pat::CandKinResolution::Parametrization  parametrization,
const AlgebraicVector4 &  parameters,
const math::XYZTLorentzVector &  initialP4 
)

Given a choice of coordinate frame, a vector of coordinates and a reference 4-vector, produce a new 4 vector with the specified parameters. The new 4-vector is not guaranteed to satisfy the constraints of these parametrization if the initial 4-vector does not satisfy them. In the future this method will throw an exception if you go in an unphysical point of the coordinate system (e.g. E^2 < P^2)

Definition at line 152 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MomDev, AlCaHLTBitMon_ParallelJobs::p, polarP4fromParameters(), run_regression::ret, and mathSSE::sqrt().

154  {
155  math::XYZTLorentzVector ret;
156  switch (parametrization) {
157  // ======= CARTESIAN ==========
158  case pat::CandKinResolution::Cart:
159  ret.SetCoordinates(parameters[0], parameters[1], parameters[2],
160  sqrt(parameters[0]*parameters[0] +
161  parameters[1]*parameters[1] +
162  parameters[2]*parameters[2] +
163  parameters[3]*parameters[3]) );
164  break;
165  case pat::CandKinResolution::MCCart: // same as above
166  ret.SetCoordinates(parameters[0], parameters[1], parameters[2],
167  sqrt(parameters[0]*parameters[0] +
168  parameters[1]*parameters[1] +
169  parameters[2]*parameters[2] +
170  initialP4.mass()*initialP4.mass()) );
171  break;
172  case pat::CandKinResolution::ECart:
173  ret.SetCoordinates(parameters[0], parameters[1], parameters[2], parameters[3]);
174  break;
175  // ======= MomentumDeviates ==========
176  case pat::CandKinResolution::MomDev:
177  case pat::CandKinResolution::EMomDev:
178  case pat::CandKinResolution::MCMomDev:
179  case pat::CandKinResolution::EScaledMomDev:
180  {
181  ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D<double> > p = initialP4.Vect(), uz(0,0,1), uph, uth;
182  uph = uz.Cross(p).Unit();
183  uth = p.Cross(uph).Unit();
184  p *= parameters[0];
185  p += uth * parameters[1] + uph * parameters[2];
186  if (parametrization == pat::CandKinResolution::MomDev) {
187  ret.SetCoordinates(p.X(), p.Y(), p.Z(), sqrt(p.Mag2() + ROOT::Math::Square(initialP4.mass() * parameters[3])) );
188  } else if (parametrization == pat::CandKinResolution::EMomDev) {
189  ret.SetCoordinates(p.X(), p.Y(), p.Z(), parameters[3] * initialP4.energy());
190  } else if (parametrization == pat::CandKinResolution::EMomDev) {
191  ret.SetCoordinates(p.X(), p.Y(), p.Z(), p.R() * initialP4.E()/initialP4.P());
192  } else {
193  ret.SetCoordinates(p.X(), p.Y(), p.Z(), sqrt(p.Mag2() + initialP4.mass()*initialP4.mass()));
194  }
195  break;
196  }
197  // ======= ALL OTHERS ==========
198  default:
199  ret = polarP4fromParameters(parametrization, parameters, math::PtEtaPhiMLorentzVector(initialP4));
200  }
201  return ret;
202 }
Parameters.parameters
dictionary parameters
Definition: Parameters.py:2
run_regression.ret
int ret
Definition: run_regression.py:388
math::PtEtaPhiMLorentzVector
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
pat::helper::ParametrizationHelper::polarP4fromParameters
math::PtEtaPhiMLorentzVector polarP4fromParameters(pat::CandKinResolution::Parametrization parametrization, const AlgebraicVector4 &parameters, const math::PtEtaPhiMLorentzVector &initialP4)
Definition: ParametrizationHelper.cc:65
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
AlgebraicVector4 pat::helper::ParametrizationHelper::parametersFromP4 ( pat::CandKinResolution::Parametrization  parametrization,
const math::XYZTLorentzVector &  p4 
)

Returns a vector of coordinates values given a coordinate frame and a 4-vector.

Definition at line 279 of file ParametrizationHelper.cc.

References run_regression::ret, and setParametersFromP4().

Referenced by diffToParameters().

279  {
280  AlgebraicVector4 ret;
281  setParametersFromP4(parametrization, ret, p4);
282  return ret;
283 }
run_regression.ret
int ret
Definition: run_regression.py:388
p4
double p4[4]
Definition: TauolaWrapper.h:92
pat::helper::ParametrizationHelper::setParametersFromP4
void setParametersFromP4(pat::CandKinResolution::Parametrization parametrization, AlgebraicVector4 &pars, const math::XYZTLorentzVector &p4)
Set the values of the parameters for a given 4-momentum.
Definition: ParametrizationHelper.cc:266
AlgebraicVector4
ROOT::Math::SVector< double, 4 > AlgebraicVector4
Definition: AlgebraicROOTObjects.h:21
AlgebraicVector4 pat::helper::ParametrizationHelper::parametersFromP4 ( pat::CandKinResolution::Parametrization  parametrization,
const math::PtEtaPhiMLorentzVector &  p4 
)

Returns a vector of coordinates values given a coordinate frame and a 4-vector.

Definition at line 272 of file ParametrizationHelper.cc.

References run_regression::ret, and setParametersFromP4().

272  {
273  AlgebraicVector4 ret;
274  setParametersFromP4(parametrization, ret, p4);
275  return ret;
276 }
run_regression.ret
int ret
Definition: run_regression.py:388
p4
double p4[4]
Definition: TauolaWrapper.h:92
pat::helper::ParametrizationHelper::setParametersFromP4
void setParametersFromP4(pat::CandKinResolution::Parametrization parametrization, AlgebraicVector4 &pars, const math::XYZTLorentzVector &p4)
Set the values of the parameters for a given 4-momentum.
Definition: ParametrizationHelper.cc:266
AlgebraicVector4
ROOT::Math::SVector< double, 4 > AlgebraicVector4
Definition: AlgebraicROOTObjects.h:21
math::PtEtaPhiMLorentzVector pat::helper::ParametrizationHelper::polarP4fromParameters ( pat::CandKinResolution::Parametrization  parametrization,
const AlgebraicVector4 &  parameters,
const math::PtEtaPhiMLorentzVector &  initialP4 
)

Given a choice of coordinate frame, a vector of coordinates and a reference 4-vector, produce a new 4 vector with the specified parameters. The new 4-vector is not guaranteed to satisfy the constraints of these parametrization if the initial 4-vector does not satisfy them. In the future this method will throw an exception if you go in an unphysical point of the coordinate system (e.g. E^2 < P^2)

Definition at line 65 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::Invalid, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, AlCaHLTBitMon_ParallelJobs::p, run_regression::ret, pat::CandKinResolution::Spher, and mathSSE::sqrt().

Referenced by p4fromParameters().

67  {
68  math::PtEtaPhiMLorentzVector ret;
69  ROOT::Math::CylindricalEta3D<double> converter;
70  double m2;
71  switch (parametrization) {
72  // ======= CARTESIAN ==========
73  case pat::CandKinResolution::Cart:
74  ret = math::XYZTLorentzVector(parameters[0], parameters[1], parameters[2],
75  sqrt(parameters[0]*parameters[0] +
76  parameters[1]*parameters[1] +
77  parameters[2]*parameters[2] +
78  parameters[3]*parameters[3]) );
79  ret.SetM(parameters[3]); // to be sure about roundoffs
80  break;
81  case pat::CandKinResolution::MCCart:
82  ret = math::XYZTLorentzVector(parameters[0], parameters[1], parameters[2],
83  sqrt(parameters[0]*parameters[0] +
84  parameters[1]*parameters[1] +
85  parameters[2]*parameters[2] +
86  initialP4.mass()*initialP4.mass()) );
87  ret.SetM(initialP4.mass()); // to be sure about roundoffs
88  break;
89  case pat::CandKinResolution::ECart:
90  ret = math::XYZTLorentzVector(parameters[0], parameters[1], parameters[2], parameters[3]);
91  break;
92  // ======= SPHERICAL ==========
93  case pat::CandKinResolution::Spher:
94  converter = ROOT::Math::Polar3D<double>(parameters[0], parameters[1], 0);
95  ret.SetCoordinates(converter.Rho(),converter.Eta(),parameters[2],parameters[3]);
96  break;
97  case pat::CandKinResolution::MCSpher: // same as above
98  converter = ROOT::Math::Polar3D<double>(parameters[0], parameters[1], 0);
99  ret.SetCoordinates(converter.Rho(),converter.Eta(),parameters[2],initialP4.mass());
100  break;
101  case pat::CandKinResolution::ESpher: //
102  converter = ROOT::Math::Polar3D<double>(parameters[0], parameters[1], 0);
103  m2 = - parameters[0]*parameters[0] + parameters[3]*parameters[3];
104  ret.SetCoordinates(converter.Rho(),converter.Eta(),parameters[2],(m2 > 0 ? sqrt(m2) : 0.0));
105  break;
106  case pat::CandKinResolution::MCPInvSpher: //
107  converter = ROOT::Math::Polar3D<double>(1.0/parameters[0], parameters[1], 0);
108  ret.SetCoordinates(converter.Rho(),converter.Eta(),parameters[2],initialP4.mass());
109  break;
110  // ======= HEP CYLINDRICAL ==========
111  case pat::CandKinResolution::EtThetaPhi:
112  converter = ROOT::Math::Polar3D<double>(1.0, parameters[1], 0);
113  ret.SetCoordinates(parameters[0],converter.Eta(),parameters[2],0);
114  break;
115  case pat::CandKinResolution::EtEtaPhi: // as simple as that
116  ret.SetCoordinates(parameters[0], parameters[1], parameters[2], 0);
117  break;
118  // ======= MomentumDeviates ==========
119  case pat::CandKinResolution::MomDev:
120  case pat::CandKinResolution::EMomDev:
121  case pat::CandKinResolution::MCMomDev:
122  case pat::CandKinResolution::EScaledMomDev:
123  {
124  ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D<double> > p = initialP4.Vect(), uz(0,0,1), uph, uth;
125  uph = uz.Cross(p).Unit();
126  uth = p.Cross(uph).Unit();
127  p *= parameters[0];
128  p += uth * parameters[1] + uph * parameters[2];
129  if (parametrization == pat::CandKinResolution::MomDev) {
130  ret.SetCoordinates(p.Rho(), p.Eta(), p.Phi(), initialP4.mass() * parameters[3]);
131  } else if (parametrization == pat::CandKinResolution::EMomDev) {
132  double m2 = ROOT::Math::Square(parameters[3] * initialP4.energy()) - p.Mag2();
133  ret.SetCoordinates(p.Rho(), p.Eta(), p.Phi(), (m2 > 0 ? sqrt(m2) : 0.0));
134  } else if (parametrization == pat::CandKinResolution::EScaledMomDev) {
135  double m2 = ROOT::Math::Square(p.R()*initialP4.E()/initialP4.P()) - p.Mag2();
136  ret.SetCoordinates(p.Rho(), p.Eta(), p.Phi(), (m2 > 0 ? sqrt(m2) : 0.0));
137  } else if (parametrization == pat::CandKinResolution::MCMomDev) {
138  ret.SetCoordinates(p.Rho(), p.Eta(), p.Phi(), initialP4.mass());
139  }
140  break;
141  }
142  // ======= OTHER ==========
143  case pat::CandKinResolution::Invalid:
144  throw cms::Exception("Invalid parametrization") << parametrization;
145  default:
146  throw cms::Exception("Not Implemented") << "getResolEta not yet implemented for parametrization " << parametrization ;
147  }
148  return ret;
149 }
Parameters.parameters
dictionary parameters
Definition: Parameters.py:2
run_regression.ret
int ret
Definition: run_regression.py:388
math::PtEtaPhiMLorentzVector
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
template<typename T >
void pat::helper::ParametrizationHelper::setParametersFromAnyVector ( pat::CandKinResolution::Parametrization  parametrization,
AlgebraicVector4 &  pars,
const T &  p4 
)

For internal use only, so we provide only the interface. Use the 'setParametersFromP4'.

Definition at line 206 of file ParametrizationHelper.cc.

References pat::CandKinResolution::Cart, pat::CandKinResolution::ECart, pat::CandKinResolution::EMomDev, pat::CandKinResolution::EScaledMomDev, pat::CandKinResolution::ESpher, pat::CandKinResolution::EtEtaPhi, pat::CandKinResolution::EtThetaPhi, Exception, pat::CandKinResolution::Invalid, pat::CandKinResolution::MCCart, pat::CandKinResolution::MCMomDev, pat::CandKinResolution::MCPInvSpher, pat::CandKinResolution::MCSpher, pat::CandKinResolution::MomDev, and pat::CandKinResolution::Spher.

Referenced by setParametersFromP4().

208  {
209  switch (parametrization) {
210  // ======= CARTESIAN ==========
211  case pat::CandKinResolution::Cart:
212  ret[0] = p4.px(); ret[1] = p4.py(); ret[2] = p4.pz(); ret[3] = p4.mass();
213  break;
214  case pat::CandKinResolution::MCCart:
215  ret[0] = p4.px(); ret[1] = p4.py(); ret[2] = p4.pz(); ret[3] = p4.mass();
216  break;
217  case pat::CandKinResolution::ECart:
218  ret[0] = p4.px(); ret[1] = p4.py(); ret[2] = p4.pz(); ret[3] = p4.energy();
219  break;
220  // ======= SPHERICAL ==========
221  case pat::CandKinResolution::Spher:
222  ret[0] = p4.P(); ret[1] = p4.theta(); ret[2] = p4.phi(); ret[3] = p4.mass();
223  break;
224  case pat::CandKinResolution::MCSpher:
225  ret[0] = p4.P(); ret[1] = p4.theta(); ret[2] = p4.phi(); ret[3] = p4.mass();
226  break;
227  case pat::CandKinResolution::ESpher:
228  ret[0] = p4.P(); ret[1] = p4.theta(); ret[2] = p4.phi(); ret[3] = p4.energy();
229  break;
230  case pat::CandKinResolution::MCPInvSpher:
231  ret[0] = 1.0/p4.P(); ret[1] = p4.theta(); ret[2] = p4.phi(); ret[3] = p4.mass();
232  break;
233  // ======= HEP CYLINDRICAL ==========
234  case pat::CandKinResolution::EtThetaPhi:
235  ret[0] = p4.pt(); ret[1] = p4.theta(); ret[2] = p4.phi(); ret[3] = 0;
236  break;
237  case pat::CandKinResolution::EtEtaPhi:
238  ret[0] = p4.pt(); ret[1] = p4.eta(); ret[2] = p4.phi(); ret[3] = 0;
239  break;
240  // ======= DEVIATES ==========
241  case pat::CandKinResolution::MomDev:
242  case pat::CandKinResolution::EMomDev:
243  ret[0] = 1.0; ret[1] = 0.0; ret[2] = 0.0; ret[3] = 1.0;
244  break;
245  case pat::CandKinResolution::MCMomDev:
246  ret[0] = 1.0; ret[1] = 0.0; ret[2] = 0.0; ret[3] = p4.mass();
247  break;
248  case pat::CandKinResolution::EScaledMomDev:
249  ret[0] = 1.0; ret[1] = 0.0; ret[2] = 0.0; ret[3] = p4.E()/p4.P();
250  break;
251  // ======= OTHER ==========
252  case pat::CandKinResolution::Invalid:
253  throw cms::Exception("Invalid parametrization") << parametrization;
254  default:
255  throw cms::Exception("Not Implemented") << "getResolEta not yet implemented for parametrization " << parametrization ;
256  }
257 }
run_regression.ret
int ret
Definition: run_regression.py:388
p4
double p4[4]
Definition: TauolaWrapper.h:92
void pat::helper::ParametrizationHelper::setParametersFromP4 ( pat::CandKinResolution::Parametrization  parametrization,
AlgebraicVector4 &  pars,
const math::XYZTLorentzVector &  p4 
)

Set the values of the parameters for a given 4-momentum.

Definition at line 266 of file ParametrizationHelper.cc.

References setParametersFromAnyVector().

Referenced by parametersFromP4().

267  {
268  setParametersFromAnyVector(parametrization, v, p4);
269 }
p4
double p4[4]
Definition: TauolaWrapper.h:92
pat::helper::ParametrizationHelper::setParametersFromAnyVector
void setParametersFromAnyVector(pat::CandKinResolution::Parametrization parametrization, AlgebraicVector4 &pars, const T &p4)
For internal use only, so we provide only the interface. Use the &#39;setParametersFromP4&#39;.
Definition: ParametrizationHelper.cc:206
void pat::helper::ParametrizationHelper::setParametersFromP4 ( pat::CandKinResolution::Parametrization  parametrization,
AlgebraicVector4 &  pars,
const math::PtEtaPhiMLorentzVector &  p4 
)

Set the values of the parameters for a given 4-momentum.

Definition at line 260 of file ParametrizationHelper.cc.

References setParametersFromAnyVector().

261  {
262  setParametersFromAnyVector(parametrization, v, p4);
263 }
p4
double p4[4]
Definition: TauolaWrapper.h:92
pat::helper::ParametrizationHelper::setParametersFromAnyVector
void setParametersFromAnyVector(pat::CandKinResolution::Parametrization parametrization, AlgebraicVector4 &pars, const T &p4)
For internal use only, so we provide only the interface. Use the &#39;setParametersFromP4&#39;.
Definition: ParametrizationHelper.cc:206
Generated for CMSSW Reference Manual by doxygen 1.8.5