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