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