/data/refman/pasoursint/CMSSW_5_3_0/src/Alignment/TwoBodyDecay/interface/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