/data/refman/pasoursint/CMSSW_6_1_2_SLHC2/src/Alignment/TwoBodyDecay/interface/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( AlgebraicVector curv, double zMagField );

  AlgebraicVector convertCurvilinearToCartesian( 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