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