Commit 14a4610799a811dfad5f75b84e62b693770ea02f

Authored by Jürgen Knödlseder
1 parent fba08ce0

Remove the event cube bin weighting in ctmodel::fill_cube (#2694)

Showing 1 changed file with 2 additions and 1 deletions
src/ctmodel/ctmodel.cpp
@@ -891,6 +891,7 @@ void ctmodel::fill_cube(const GCTAObservation* obs, GModels& models) @@ -891,6 +891,7 @@ void ctmodel::fill_cube(const GCTAObservation* obs, GModels& models)
891 // Set event bin attributes that are constant for a counts cube 891 // Set event bin attributes that are constant for a counts cube
892 bin.time(m_time); 892 bin.time(m_time);
893 bin.ontime(obs->ontime()); 893 bin.ontime(obs->ontime());
  894 + bin.weight(1.0);
894 895
895 // Loop over all the spatial bins 896 // Loop over all the spatial bins
896 for (int i = 0; i < npix; ++i) { 897 for (int i = 0; i < npix; ++i) {
@@ -926,7 +927,7 @@ void ctmodel::fill_cube(const GCTAObservation* obs, GModels&amp; models) @@ -926,7 +927,7 @@ void ctmodel::fill_cube(const GCTAObservation* obs, GModels&amp; models)
926 bin.energy(m_energy[iebin]); 927 bin.energy(m_energy[iebin]);
927 bin.ewidth(m_ewidth[iebin]); 928 bin.ewidth(m_ewidth[iebin]);
928 bin.ieng(iebin); 929 bin.ieng(iebin);
929 - bin.weight(m_cube[ibin]->weight()); 930 + //bin.weight(m_cube[ibin]->weight()); // Bad residuals for H.E.S.S. Crab
930 931
931 // Compute model value for cube bin 932 // Compute model value for cube bin
932 double model = models.eval(bin, *obs) * bin.size(); 933 double model = models.eval(bin, *obs) * bin.size();