TwoBodyDecayDerivatives.h

Go to the documentation of this file.

#ifndef Alignment_TwoBodyDecay_TwoBodyDecayDerivatives_h
#define Alignment_TwoBodyDecay_TwoBodyDecayDerivatives_h
 
#include "Alignment/TwoBodyDecay/interface/TwoBodyDecay.h"
 
class TwoBodyDecayDerivatives {
public:
  enum { dimension = 6 };
 
  enum DerivativeParameterName { px = 1, py = 2, pz = 3, theta = 4, phi = 5, mass = 6 };
 
  TwoBodyDecayDerivatives(double mPrimary = 91.1876, double mSecondary = 0.105658);
  ~TwoBodyDecayDerivatives();
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> derivatives(const TwoBodyDecay &tbd) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> derivatives(const TwoBodyDecayParameters &param) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> selectedDerivatives(const TwoBodyDecay &tbd,
                                                                        const std::vector<bool> &selector) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> selectedDerivatives(const TwoBodyDecayParameters &param,
                                                                        const std::vector<bool> &selector) const;
 
private:
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> dqsdpx(const TwoBodyDecayParameters &param) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> dqsdpy(const TwoBodyDecayParameters &param) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> dqsdpz(const TwoBodyDecayParameters &param) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> dqsdtheta(const TwoBodyDecayParameters &param) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> dqsdphi(const TwoBodyDecayParameters &param) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> dqsdm(const TwoBodyDecayParameters &param) const;
 
  const std::pair<AlgebraicMatrix, AlgebraicMatrix> dqsdzi(const TwoBodyDecayParameters &param,
                                                           const DerivativeParameterName &i) const;
 
  double thePrimaryMass;
  double theSecondaryMass;
};
 
#endif