Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/externals/g4tools/include/tools/histo/axes

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

  1 // Copyright (C) 2010, Guy Barrand. All rights reserved.
  2 // See the file tools.license for terms.
  3 
  4 #ifndef tools_histo_axes
  5 #define tools_histo_axes
  6 
  7 #include "axis"
  8 
  9 namespace tools {
 10 namespace histo {
 11 
 12 //TC is for a coordinate.
 13 //TO is for an offset used to identify a bin.
 14 
 15 template <class TC,class TO>
 16 inline bool is_out(const std::vector< axis<TC,TO> >& a_axes,TO a_offset) {
 17   TO offset = a_offset;
 18   int index;
 19   typename std::vector< axis<TC,TO> >::size_type dimension = a_axes.size();
 20   for(int iaxis=int(dimension)-1;iaxis>=0;iaxis--) {
 21     index = int(offset/a_axes[iaxis].m_offset);
 22     if(index==0) return true;
 23     if(index==(int(a_axes[iaxis].m_number_of_bins)+1)) return true;
 24     offset -= index * a_axes[iaxis].m_offset;
 25   }
 26   return false;
 27 }
 28 
 29 template <class TC,class TO>
 30 inline void get_indices(const std::vector< axis<TC,TO> >& a_axes,TO a_offset,std::vector<int>& a_is) {
 31   TO offset = a_offset;
 32   typename std::vector< axis<TC,TO> >::size_type dimension = a_axes.size();
 33  {for(int iaxis=int(dimension)-1;iaxis>=0;iaxis--) {
 34     a_is[iaxis] = int(offset/a_axes[iaxis].m_offset);
 35     offset -= a_is[iaxis] * a_axes[iaxis].m_offset;
 36   }}
 37   typedef unsigned int dim_t;
 38   for(dim_t iaxis=0;iaxis<dimension;iaxis++) {
 39     if(a_is[iaxis]==0) {
 40       a_is[iaxis] = axis_UNDERFLOW_BIN;
 41     } else if(a_is[iaxis]==int(a_axes[iaxis].m_number_of_bins)+1) {
 42       a_is[iaxis] = axis_OVERFLOW_BIN;
 43     } else {
 44       a_is[iaxis]--;
 45     }
 46   }
 47 }
 48 
 49 template <class TC,class TO>
 50 inline bool get_offset(const std::vector< axis<TC,TO> >& a_axes,const std::vector<int>& a_is,TO& a_offset) {
 51   // a_is[iaxis] is given in in-range indexing :
 52   //  - [0,n[iaxis]-1] for in-range bins
 53   //  - UNDERFLOW_BIN for the iaxis underflow bin
 54   //  - OVERFLOW_BIN for the iaxis overflow bin
 55   a_offset = 0;
 56   if(a_axes.empty()) return false;
 57   typename std::vector< axis<TC,TO> >::size_type dimension = a_axes.size();
 58   typename axis<TC,TO>::bn_t ibin;
 59   typedef unsigned int dim_t;
 60   for(dim_t iaxis=0;iaxis<dimension;iaxis++) {
 61     if(!a_axes[iaxis].in_range_to_absolute_index(a_is[iaxis],ibin)) {
 62       a_offset = 0;
 63       return false;
 64     }
 65     a_offset += ibin * a_axes[iaxis].m_offset;
 66   }
 67   return true;
 68 }
 69 
 70 }}
 71 
 72 #endif
 73 
 74 
 75 
 76