ThreeBodyCalculators.h

Go to the documentation of this file.

#ifndef AMPGEN_THREEBODYCALCULATORS_H
#define AMPGEN_THREEBODYCALCULATORS_H
 
#include "AmpGen/DalitzIntegrator.h"
#include "AmpGen/Expression.h"
#include "AmpGen/CoherentSum.h"
#include "AmpGen/Tensor.h"
 
class TGraph;
 
namespace AmpGen
{
  class MinuitParameterSet;
 
  class ThreeBodyCalculator
  {
  private:
    struct PartialWidth {
      CoherentSum fcs;
      DalitzIntegrator integrator;
      CompiledExpression<complex_v(const real_t*, const real_v*) > totalWidth;
      EventType type;
      std::vector<CompiledExpression< complex_v(const real_t*, const real_v*)>> partialWidths; 
      double getWidth( const double& m );
      PartialWidth( const EventType& type, MinuitParameterSet& mps );
      Expression spinAverageMatrixElement( const std::vector<std::pair<Particle, TotalCoupling> >& elements,
                                           DebugSymbols* msym );
    };
    Expression calculateSAME( const std::string& particle );
 
    double       m_min;
    double       m_max;
    double       m_norm;
    double       m_step;
    size_t       m_nKnots;
    std::string  m_name;
    std::vector<PartialWidth> m_widths;
    MinuitParameterSet* m_mps; 
  public:
    ThreeBodyCalculator( const std::string& head, MinuitParameterSet& mps, const size_t& nKnots=999, const double& min=-1, const double& max=-1 );
 
    TGraph* widthGraph( const double& mNorm=-1 );
    TGraph* widthGraph( const size_t& steps, const double& min, const double& max );
    TGraph* runningMass( const double& mass, const double& min, const double& max, const size_t& nSteps, const size_t& nSubtractions=2 );
    TGraph* fastRunningMass( const double& mass, const double& min, const double& max, const size_t& nSteps, const size_t& nSubtractions=2 );
 
    double getWidth( const double& );
 
    void updateRunningWidth( MinuitParameterSet& mps, const double& mNorm = 0 );    
    void setNorm( const double& mNorm );
    void setAxis( const size_t& nKnots, const double& min, const double& max );
    void prepare();
    void makePlots(const double& mass=-1, const size_t& x=0, const size_t& y=0);
    void debug( const double& m, const double& theta );
  };
} // namespace AmpGen
#endif /* end of include guard: AMPGEN_THREEBODYCALCULATORS_H */