Develop

examples/phys/FlatteLineShape.inl

@@ -252,7 +252,7 @@ public:
     Flatte() = delete;
 
     Flatte(hydra::Parameter const& c_re, hydra::Parameter const& c_im, hydra::Parameter const& mean,
-           std::vector<std::vector<double>> const& params,double mother_mass,	double daugther1_mass, double daugther2_mass,
+           std::array<std::array<double,3>,2> const &params,double mother_mass,	double daugther1_mass, double daugther2_mass,
            double daugther3_mass, double radi):
             hydra::BaseFunctor<Flatte<CHANNEL,L>, hydra::complex<double>, 3>{c_re, c_im, mean},
             fLineShape(mean,params,mother_mass,daugther1_mass,daugther2_mass,daugther3_mass,radi)
@@ -360,7 +360,8 @@ int main(int argv, char** argc)
 
     //Flatté
     double pi_MASS = 0.13957018;
-    std::vector<std::vector<double>> params = {{pi_MASS,pi_MASS,f0_rho1},{Kplus_MASS,Kplus_MASS,f0_rho2}};
+    std::array<std::array<double,3>,2> params = {{{pi_MASS,pi_MASS,f0_rho1},{Kplus_MASS,Kplus_MASS,f0_rho2}}};
+    //std::vector<std::vector<double>> params = {{pi_MASS,pi_MASS,f0_rho1},{Kplus_MASS,Kplus_MASS,f0_rho2}};
 
 	//======================================================
 	//Phi

hydra/functions/FlatteLineShape.h

@@ -80,7 +80,7 @@ namespace hydra {
  * @param radi decay vertex radio.
  */
         FlatteLineShape(hydra::Parameter const& mean,
-                        std::vector<std::vector<double>> params,
+                        std::array<std::array<double,3>,2> params,
                         double mother_mass,
                         double daugther1_mass, double daugther2_mass, double daugther3_mass,
                         double radi) :
@@ -121,12 +121,12 @@ namespace hydra {
         }
 
         __hydra_host__  __hydra_device__ inline
-        std::vector<std::vector<double>>GetParams() const {
+        std::array<std::array<double,3>,2> GetParams() const {
             return fParams;
         }
 
         __hydra_host__  __hydra_device__ inline
-        void SetDaughter1Mass(std::vector<std::vector<double>> _params) {
+        void SetDaughter1Mass(std::array<std::array<double,3>,2> _params) {
             fParams = _params;
         }
 
@@ -204,32 +204,34 @@ namespace hydra {
         }
 
     private:
+
         __hydra_host__ __hydra_device__ inline
         hydra::complex<double>
         sqrtCplx(const double in) const { return (in > 0) ? hydra::complex<double>(::sqrt(in), 0.0) : hydra::complex<double>(0.0, ::sqrt(-in)); }
 
-        __hydra_host__ __hydra_device__ inline
-        hydra::complex<double>
-        LineShape(const double s, const double resonance_mass) const {
 
-		hydra::complex<double> w;
+        __hydra_host__ __hydra_device__ inline hydra::complex<double> LineShape(const double s, const double resonance_mass) const  {
 
-for(size_t i = 0; i < fParams.size() ; i++) {
+            hydra::complex<double> w;
 
-        double m1a = fParams[i][0];
-        double m1b = fParams[i][1];
-        double g   = fParams[i][2];
 
-        w += g * g * sqrtCplx((1 - (m1a - m1b) * (m1a - m1b) / s*s ) *
-                          (1 - (m1a + m1b) * (m1a + m1b) / s*s ));
 
-}
+            for(size_t i = 0; i < fParams.size() ; i++) {
+
+                double m1a = fParams[i][0];
+                double m1b = fParams[i][1];
+                double g   = fParams[i][2];
+
+                w += g * g * sqrtCplx((1 - (m1a - m1b) * (m1a - m1b) / s*s ) *
+                                      (1 - (m1a + m1b) * (m1a + m1b) / s*s ));
+
+            }
 
-		hydra::complex<double> denom =   resonance_mass - s*s - hydra::complex<double>(0.0,1.0)*w;
+            hydra::complex<double> denom =   resonance_mass - s*s - hydra::complex<double>(0.0,1.0)*w;
 
-		hydra::complex<double> ampl = hydra::complex<double>(1.0,0.0)/denom;
+            hydra::complex<double> ampl = hydra::complex<double>(1.0,0.0)/denom;
 
-		return ampl;
+            return ampl;
 
         }
 
@@ -238,7 +240,7 @@ for(size_t i = 0; i < fParams.size() ; i++) {
         double fBachelorMass;
         double fMotherMass;
         double fRadi;
-        std::vector<std::vector<double>> fParams;
+        std::array<std::array<double,3>,2> fParams;
 
     };