TwoBodyDecayModel.h

Go to the documentation of this file.

#ifndef Alignment_TwoBodyDecay_TwoBodyDecayModel_h
#define Alignment_TwoBodyDecay_TwoBodyDecayModel_h

#include "Alignment/TwoBodyDecay/interface/TwoBodyDecay.h"

class TwoBodyDecayModel {
public:
  TwoBodyDecayModel(double mPrimary = 91.1876, double mSecondary = 0.105658);
  ~TwoBodyDecayModel();

  AlgebraicMatrix rotationMatrix(double px, double py, double pz);

  AlgebraicMatrix curvilinearToCartesianJacobian(double rho, double theta, double phi, double zMagField);

  AlgebraicMatrix curvilinearToCartesianJacobian(const AlgebraicVector &curv, double zMagField);

  AlgebraicVector convertCurvilinearToCartesian(const AlgebraicVector &curv, double zMagField);

  const std::pair<AlgebraicVector, AlgebraicVector> cartesianSecondaryMomenta(const AlgebraicVector &param);

  const std::pair<AlgebraicVector, AlgebraicVector> cartesianSecondaryMomenta(const TwoBodyDecay &tbd);

  const std::pair<AlgebraicVector, AlgebraicVector> cartesianSecondaryMomenta(const TwoBodyDecayParameters &tbdparam);

private:
  double thePrimaryMass;
  double theSecondaryMass;
};

#endif