added FVector PVector etc. II

PwaDynamics/KMatrixPiPiS.cc

+#include "PwaDynamics/KMatrixPiPiS.hh"
+#include "PwaDynamics/KPole.hh"
+#include "PwaDynamics/AbsPhaseSpace.hh"
+#include "PwaDynamics/PhaseSpaceIsobar.hh"
+#include "PwaDynamics/PhaseSpace4Pi.hh"
+#include "qft++/relativistic-quantum-mechanics/Utils.hh"
+#include "qft++/matrix/IdentityMatrix.hh"
+#include <boost/multi_array.hpp>
+
+
+KMatrixPiPiS::KMatrixPiPiS() :
+  KMatrixSlowAdlerCorRel(5,5)
+{
+  initASParam1900();
+}
+
+
+KMatrixPiPiS::~KMatrixPiPiS(){
+}
+
+void  KMatrixPiPiS::initASParam1900(){
+  const double piMass=0.1349766;
+  const double KplusMass=0.493677;
+  const double K0Mass= 0.497614;
+  const double etaMass= 0.547853;
+  const double etaprimeMass= 0.95778;
+
+  std::map<int, std::vector<double> > gFactorsMap;
+  std::vector<double> poleMasses;
+  poleMasses.resize(5);
+  poleMasses[0]=0.651;
+  poleMasses[1]=1.2036;
+  poleMasses[2]=1.55817;
+  poleMasses[3]=1.21;
+  poleMasses[4]=1.82206;
+
+  gFactorsMap[0].push_back(0.22889); //pipi
+  gFactorsMap[0].push_back(-0.55377); //KK
+  gFactorsMap[0].push_back(0.); //4pi
+  gFactorsMap[0].push_back(-0.39899); //eta eta
+  gFactorsMap[0].push_back(-0.34639); //eta eta'
+
+  gFactorsMap[1].push_back(0.94128); //pipi
+  gFactorsMap[1].push_back(0.55095); //KK
+  gFactorsMap[1].push_back(0.); //4pi
+  gFactorsMap[1].push_back(0.39065); //eta eta
+  gFactorsMap[1].push_back(0.31503); //eta eta'
+
+  gFactorsMap[2].push_back(0.36856); //pipi
+  gFactorsMap[2].push_back(0.23888); //KK
+  gFactorsMap[2].push_back(0.55639); //4pi
+  gFactorsMap[2].push_back(0.18340); //eta eta
+  gFactorsMap[2].push_back(0.18681); //eta eta'
+
+  gFactorsMap[3].push_back(0.3365); //pipi
+  gFactorsMap[3].push_back(0.40907); //KK
+  gFactorsMap[3].push_back(0.85679); 
+  gFactorsMap[3].push_back(0.19906); //eta eta
+  gFactorsMap[3].push_back(-0.00984); //eta eta'
+
+  gFactorsMap[4].push_back(0.18171); //pipi
+  gFactorsMap[4].push_back(-0.17558); //KK
+  gFactorsMap[4].push_back(-0.79658); //4pi
+  gFactorsMap[4].push_back(-0.00355); //eta eta
+  gFactorsMap[4].push_back(0.22358); //eta eta'  
+
+  boost::shared_ptr<AbsPhaseSpace> pipiPhp(new PhaseSpaceIsobar(piMass, piMass));
+  boost::shared_ptr<AbsPhaseSpace> kkPhp(new PhaseSpaceIsobar(KplusMass, K0Mass));
+  boost::shared_ptr<AbsPhaseSpace> pipipipiPhp(new PhaseSpace4Pi());
+  boost::shared_ptr<AbsPhaseSpace> etaetaPhp(new PhaseSpaceIsobar(etaMass, etaMass));
+  boost::shared_ptr<AbsPhaseSpace> etaetapPhp(new PhaseSpaceIsobar(etaMass, etaprimeMass));
+
+  _phpVecs.push_back(pipiPhp);
+  _phpVecs.push_back(kkPhp);
+  _phpVecs.push_back(pipipipiPhp);
+  _phpVecs.push_back(etaetaPhp);
+  _phpVecs.push_back(etaetapPhp);
+
+  //  boost::shared_ptr< boost::multi_array<double, 2> > fScatMatr(new boost::multi_array<double, 2>(boost::extents[gFactors.size()][gFactors.size()]));
+
+  _fScatPtr=boost::shared_ptr< boost::multi_array<double, 2> >(new boost::multi_array<double, 2>(boost::extents[NumRows()][NumCols()]));
+
+  for (int i=0; i<NumRows(); ++i){
+    for (int j=0; j<NumCols(); ++j){
+      (*_fScatPtr)[i][j]=0.;
+    }
+  }
+  (*_fScatPtr)[0][0]=0.23399;
+  (*_fScatPtr)[0][1]=0.15044;
+  (*_fScatPtr)[0][2]=-0.20545;
+  (*_fScatPtr)[0][3]=0.32825;
+  (*_fScatPtr)[0][4]=0.35412;
+
+  (*_fScatPtr)[1][0]=(*_fScatPtr)[0][1];
+  (*_fScatPtr)[2][0]=(*_fScatPtr)[0][2];
+  (*_fScatPtr)[3][0]=(*_fScatPtr)[0][3];
+  (*_fScatPtr)[4][0]=(*_fScatPtr)[0][4];
+
+  std::vector< boost::shared_ptr<KPole> > kPoles;
+
+  std::map<int, std::vector<double> >::iterator itgFac;
+  for (itgFac=gFactorsMap.begin(); itgFac!=gFactorsMap.end(); ++itgFac){
+    boost::shared_ptr<KPole> currentPole(new KPole(itgFac->second, poleMasses[itgFac->first]));
+    _KPoles.push_back(currentPole);
+  }
+
+  _s0Scat=-3.92637;
+  _sAdler0=-0.15;
+  _sAdler=1.;
+
+
+}

PwaDynamics/KMatrixPiPiS.hh

+//KMatrixPiPiS class definition file. -*- C++ -*-
+// Copyright 2010 Bertram Kopf
+
+#pragma once 
+
+//_____________________________________________________________________________
+// @file KMatrixPiPiS.h
+//_____________________________________________________________________________
+
+#include "PwaDynamics/KMatrixSlowAdlerCorRel.hh"
+#include <iostream>
+#include <vector>
+#include <boost/shared_ptr.hpp>
+#include <boost/multi_array.hpp>
+
+
+typedef boost::multi_array< double, 2> array_type_2d;
+
+using namespace std;
+//_____________________________________________________________________________
+//_____________________________________________________________________________
+
+class KMatrixPiPiS : public KMatrixSlowAdlerCorRel {
+
+public:
+
+  /// Constructor 
+ KMatrixPiPiS(); 
+
+  /// Destructor
+  virtual ~KMatrixPiPiS();
+
+protected:
+  void initASParam1900();
+private:
+};
+//_____________________________________________________________________________
+
+