added more output

Domenico Tiziani
1 parent 71062026
Showing 3 changed files with 26 additions and 12 deletions
inst/cta/src/GCTACubePsf.cpp
inst/cta/src/GCTAResponseTable.cpp
src/support/GNodeArray.cpp
@@ -845,7 +845,18 @@ void GCTACubePsf::fill_cube(const GCTAObservation&amp; obs,
                     int imap = offset(idelta, iebin);
  
                     // supress output
-                    //std::cout.setstate(std::ios_base::failbit);
+                    std::cout.setstate(std::ios_base::failbit);
+
+                    if( (iebin == m_energies.size()-2 || iebin == m_energies.size()-3) && idelta == 0){
+                        if(pixel > 1521 && pixel < 1525){
+                            std::cout.clear();
+                            //std::cout << "delta = " << delta * gammalib::rad2deg << std::endl;
+                            //std::cout << "theta  = " << theta  * gammalib::rad2deg << std::endl;
+                            //std::cout << "psf value " << rsp->psf(delta, theta, 0.0, 0.0, 0.0, logE) <<" pixel " <<  pixel << std::endl;
+                            //std::cout << "--------------------------" << std::endl;
+                            //std::cout.setstate(std::ios_base::failbit);
+                        }
+                    }
  
                     // Add on PSF cube
                     m_cube(pixel, imap) +=
@@ -853,16 +864,7 @@ void GCTACubePsf::fill_cube(const GCTAObservation&amp; obs,
                     std::cout << "energy bin " << iebin << " of " << m_energies.size() << std::endl;
                     std::cout << "psf value " << rsp->psf(delta, theta, 0.0, 0.0, 0.0, logE) <<" pixel " <<  pixel << std::endl;
  
-                    if( (iebin == m_energies.size()-2 || iebin == m_energies.size()-3) && idelta == 0){
-                    	if(pixel > 1521 && pixel < 1525){
-                    		//std::cout.clear();
-                    		//std::cout << "delta = " << delta * gammalib::rad2deg << std::endl;
-                    		//std::cout << "theta  = " << theta  * gammalib::rad2deg << std::endl;
-                    		//std::cout << "psf value " << rsp->psf(delta, theta, 0.0, 0.0, 0.0, logE) <<" pixel " <<  pixel << std::endl;
-                    		//std::cout << "--------------------------" << std::endl;
-                    		//std::cout.setstate(std::ios_base::failbit);
-                    	}
-                    }
+
                     //std::cout << "Energy: " << iebin << " == " << m_energies.size()-1 << std::endl;
                 } // endfor: looped over delta bins
  
@@ -288,6 +288,13 @@ std::vector&lt;double&gt; GCTAResponseTable::operator()(const double&amp; arg1,
                     m_wgt2 * m_tables[i][m_inx2] +
                     m_wgt3 * m_tables[i][m_inx3] +
                     m_wgt4 * m_tables[i][m_inx4];
+        if (i == 4) {
+            std::cout << "Values  for entry 4:" << std::endl;
+            std::cout << m_wgt1 << " * " << m_tables[i][m_inx1] << std::endl;
+            std::cout << m_wgt2 << " * " << m_tables[i][m_inx2] << std::endl;
+            std::cout << m_wgt3 << " * " << m_tables[i][m_inx3] << std::endl;
+            std::cout << m_wgt4 << " * " << m_tables[i][m_inx4] << std::endl;
+        }
     }
  
     // Return result vector
@@ -1984,7 +1991,8 @@ void GCTAResponseTable::update(const double&amp; arg1, const double&amp; arg2) const
     m_wgt2 = nodes1->wgt_left()  * nodes2->wgt_right();
     m_wgt3 = nodes1->wgt_right() * nodes2->wgt_left();
     m_wgt4 = nodes1->wgt_right() * nodes2->wgt_right();
-
+    std::cout << "Energy weights: " << nodes1->wgt_left() << ", " << nodes1->wgt_right() << std::endl;
+    //std::cout <<
     // Return
     return;
 }
@@ -675,6 +675,10 @@ void GNodeArray::set_value(const double&amp; value) const
             m_inx_right = m_inx_left + 1;
  
             // Set weighting factors
+            std::cout << "Calculating weights: (" << value << " - " << m_node[m_inx_left] <<
+                    ") / " << m_step[m_inx_left] << std::endl;
+            std::cout << "Last: " << m_node[nodes-1] << std::endl;
+
             m_wgt_right = (value - m_node[m_inx_left]) / m_step[m_inx_left];
             m_wgt_left  = 1.0 - m_wgt_right;