diff --git a/core/CMakeLists.txt b/core/CMakeLists.txt
index eb7c93b..5d144fb 100644
--- a/core/CMakeLists.txt
+++ b/core/CMakeLists.txt
@@ -84,9 +84,7 @@ target_link_libraries(${ELSPECTRO} ${ROOT_LIBRARIES} )
 install(TARGETS ${ELSPECTRO} 
   LIBRARY DESTINATION "${CMAKE_INSTALL_LIBDIR}")
 
-FILE(GLOB ROOT_ELSPECTRO_PCM  ${CMAKE_BINARY_DIR}/${DIRNAME}/libelSpectro_rdict.pcm)
-
-install (FILES  ${ROOT_ELSPECTRO_PCM}    
+install (FILES  ${CMAKE_CURRENT_BINARY_DIR}/libelSpectro_rdict.pcm    
   DESTINATION "${CMAKE_INSTALL_LIBDIR}")
 
 #include_directories(${CMAKE_CURRENT_SOURCE_DIR}/src)
diff --git a/core/DecayModelst.cpp b/core/DecayModelst.cpp
index 6a6e404..b21afbd 100644
--- a/core/DecayModelst.cpp
+++ b/core/DecayModelst.cpp
@@ -152,7 +152,7 @@ namespace elSpectro{
 	for(int it=0;it<Npoints;it++){
 	  WtVals[1]=tming-it*trange/Npoints;
 	  auto val = Fmax(WtVals);
-	  if(it==0)std::cout<<WtVals[0]<<" "<<WtVals[1]<<" val "<<val<<" "<<PhaseSpaceFactor()<<" s "<<_s<<" pgam "<<PgammaCMsq()<<" at Q2 0  = "<<kine::PDK2(_W,0,_target->M())<<std::endl;
+	  //if(it==0)std::cout<<WtVals[0]<<" "<<WtVals[1]<<" val "<<val<<" "<<PhaseSpaceFactor()<<" s "<<_s<<" pgam "<<PgammaCMsq()<<" at Q2 0  = "<<kine::PDK2(_W,0,_target->M())<<std::endl;
 	  if(val<gridMin) {
 	    gridMin=val;
 	    gridW=WtVals[0];
diff --git a/core/DistVirtPhotFlux_xy.h b/core/DistVirtPhotFlux_xy.h
index 0d0ecc5..42a53fc 100644
--- a/core/DistVirtPhotFlux_xy.h
+++ b/core/DistVirtPhotFlux_xy.h
@@ -138,6 +138,7 @@ namespace elSpectro{
     return escat::flux_dlnxdlny(_ebeam,lnx,lny);
     }*/
   inline double DistVirtPhotFlux_xy::Eval(const double *x) const{
+    //given x and y evaluate photon flux
 
     double lnx=TMath::Log(x[0]);
     double lny=TMath::Log(x[1]);
diff --git a/core/ElectronScattering.cpp b/core/ElectronScattering.cpp
index 4c61de3..1f99504 100644
--- a/core/ElectronScattering.cpp
+++ b/core/ElectronScattering.cpp
@@ -397,8 +397,8 @@ namespace elSpectro{
     gBenchmark->Stop("RooFitIntegral");
     gBenchmark->Print("RooFitIntegral");
     
-    std::cout<<" ElectronScattering::IntegrateCrossSection() ROOFIT "<<RFintegral<<" "<<" giving a photon flux weighted average photoproduction cross section of "<<RFintegral/photonFlux->Dist().Integral()<<std::endl;
-    std::cout<<" W range "<<_Wmin<<" "<< collision.M() <<" "<< ( collision.M()- _Wmin)<<std::endl;
+    std::cout<<" ElectronScattering::IntegrateCrossSection()  "<<RFintegral<<" nb "<<std::endl<<" giving a photon flux weighted average photoproduction cross section of "<<RFintegral/photonFlux->Dist().Integral()<<" nb"<<std::endl;
+    std::cout<<" W range "<<_Wmin<<" - "<< collision.M() <<" =  "<< ( collision.M()- _Wmin)<<std::endl;
     xvar.Print();
     yvar.Print();
     cthvar.Print();
diff --git a/examples/EIC_JPAC_nZc_Hists.C b/examples/EIC_JPAC_nZc_Hists.C
index 2d060ae..39c6e71 100644
--- a/examples/EIC_JPAC_nZc_Hists.C
+++ b/examples/EIC_JPAC_nZc_Hists.C
@@ -125,7 +125,6 @@ void EIC_JPAC_nZc_Hists(string ampPar="high",double ebeamE = 5, double pbeamE =
   production->SetCombinedBranchingFraction(0.06); //Just Jpsi->e+e-
   generator().SetNEvents_via_LuminosityTime(nLumi,24*60*60*nDays);
 
-  //or can just do generator().SetNEvents(1E6);
   auto fastIntegral=production->IntegrateCrossSectionFast();
   std::cout<<"       check fast cross section "<<fastIntegral<<std::endl;