Merge pull request #89 from dbrundu/master

Solving small bugs - issues #87 and #88

Author: AAAlvesJr

hydra/SPlot.h

@@ -118,7 +118,7 @@ class SPlot: public detail::AddPdfBase<PDF1,PDF2,PDFs...>
 		//fCovMatrix << 0.0, 0.0, 0.0, 0.0;
 
 
-		Eigen::Matrix<double, npdfs, npdfs>  init{};
+		Eigen::Matrix<double, npdfs, npdfs>  init = Eigen::Matrix<double, npdfs, npdfs>::Zero();
 		//init << 0.0, 0.0, 0.0, 0.0;
 
 		fCovMatrix = hydra_thrust::transform_reduce(system_type(), first, last,

hydra/detail/functors/DecayMother.h

@@ -180,7 +180,7 @@ struct DecayMother
 	GReal_t process(size_t evt, Vector4R (&daugters)[N])
 	{
 
-		GRND randEng;//(  );
+		GRND randEng(fSeed);
 		randEng.discard(evt+3*N);
 		hydra_thrust::uniform_real_distribution<GReal_t> uniDist(0.0, 1.0);
 

hydra/detail/functors/ProcessSPlot.h

@@ -203,7 +203,7 @@ struct SWeights
 		Eigen::Matrix<double, nfunctors,1> values_vector(Eigen::Map<Eigen::Matrix<double, nfunctors,1> >(values).eval());
 		Eigen::Matrix<double, nfunctors,1> sweights(fICovMatrix*values_vector.eval());
 		auto wfvalues = detail::multiply_array_tuple(fCoeficients, fvalues);
-		GReal_t denominator   = 1;
+		GReal_t denominator   = 0.0;
 		detail::add_tuple_values(denominator, wfvalues);
 		sweights /=denominator;