/[PAMELA software]/PamCut/CollectionActions/Histo2DActions/Histo2DAction/Histo2DAction.h
ViewVC logotype

Diff of /PamCut/CollectionActions/Histo2DActions/Histo2DAction/Histo2DAction.h

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by pam-fi, Fri Sep 25 15:36:45 2009 UTC revision 1.5 by pam-fi, Tue Dec 8 17:38:39 2009 UTC
# Line 106  public: Line 106  public:
106     */     */
107    void SetYAxis(TString label, unsigned int nBins, float min, float max, bool logBinning = false);    void SetYAxis(TString label, unsigned int nBins, float min, float max, bool logBinning = false);
108    
109      /*! @brief Sets the ROOT histogram's title.
110       *
111       * @param title The histogram title as it will appear on the histogram itself.
112       */
113      void SetTitle(TString &title) {
114        _title = title;
115      }
116    /*! @brief Sets up the histogram    /*! @brief Sets up the histogram
117     *     *
118     * This routine effectively prepares the histogram, after the desired parameters has been set by #SetXAxis() and #SetYAxis().     * This routine effectively prepares the histogram, after the desired parameters has been set by #SetXAxis() and #SetYAxis().
119     *     *
120     * @param events Pointer to PamLevel2 events (unused).     * @param events Pointer to PamLevel2 events (unused).
121     */     */
   
   /*! @brief Sets the ROOT histogram's title. */  
   void SetTitle(TString &title){  
     _title = title;  
   }  
122    void Setup(PamLevel2 *events) {    void Setup(PamLevel2 *events) {
123      CollectionAction::Setup(events);      CollectionAction::Setup(events);
124      _InitHistos();      _InitHistos();
# Line 141  public: Line 143  public:
143      return _histo;      return _histo;
144    }    }
145    
146    /*! Fills the ROOT and the vector histogram. */    /*! @brief Returns a pointer to the ROOT histogram.
147       *
148       * @return A pointer to the root histogram
149       */
150      TH2 *GetRootHisto() {
151        return _rootHisto;
152      }
153    
154      /*! Fills the ROOT and the vector histogram.
155       *
156       * @param xValue The value of the X coordinate associated to the event.
157       * @param yValue The value of the Y coordinate associated to the event.
158       * @param weight The weight which will be applied to the event.
159       */
160    void Fill(double xValue, double yValue, double weight = 1.);    void Fill(double xValue, double yValue, double weight = 1.);
161    
162    /*! @brief Gets the X overflow histogram.    /*! @brief Gets the X overflow histogram.
# Line 236  public: Line 251  public:
251    
252  protected:  protected:
253    
254      /*! @brief The vector containing the limits of the X bins(from lower to higher). */
255    std::vector<float> _xBins;    std::vector<float> _xBins;
256      /*! @brief The vector containing the limits of the Y bins(from lower to higher). */
257    std::vector<float> _yBins;    std::vector<float> _yBins;
258      /*! @brief A matrix containing the value of the histogram for each X-Y bin. */
259    SimpleMatrix<HistoType> _histo;    SimpleMatrix<HistoType> _histo;
260      /*! @brief The ROOT histogram. */
261    TH2 *_rootHisto;    TH2 *_rootHisto;
262    
263  private:  private:
# Line 262  void Histo2DAction<HistoType>::_CreateHi Line 281  void Histo2DAction<HistoType>::_CreateHi
281    _rootHisto = NULL;    _rootHisto = NULL;
282  }  }
283    
284    // Specializations for _CreateHistos(). See Histo2DAction.cpp
285    template<>
286    void Histo2DAction<Int_t>::_CreateHisto();
287    
288    template<>
289    void Histo2DAction<Float_t>::_CreateHisto();
290    
291    template<>
292    void Histo2DAction<Double_t>::_CreateHisto();
293    
294  template<class HistoType>  template<class HistoType>
295  void Histo2DAction<HistoType>::_InitHistos() {  void Histo2DAction<HistoType>::_InitHistos() {
296    
297    _CreateHisto();    _CreateHisto();
298      if (_xBins.size() < 2) // SetXAxis not called by the main program, or wrongly filled (only 1 bin limit)
299    if (_xBins.size() < 2){      SetXAxis("Default X", 10, 0., 1.);
300      _xBins.resize(2);    if (_yBins.size() < 2) // SetYAxis not called by the main program, or wrongly filled (only 1 bin limit)
301      _xBins[0] = 0.;      SetYAxis("Default Y", 10, 0., 1.);
     _xBins[1] = 1.;  
   }  
   
   if (_yBins.size() < 2){  
     _yBins.resize(2);  
     _yBins[0] = 0.;  
     _yBins[1] = 1.;  
   }  
302    
303    if (_rootHisto) {    if (_rootHisto) {
304      Double_t *auxXArray = new Double_t[_xBins.size()];      Double_t *auxXArray = new Double_t[_xBins.size()];
# Line 314  void Histo2DAction<HistoType>::_InitHist Line 335  void Histo2DAction<HistoType>::_InitHist
335    
336  template<class HistoType>  template<class HistoType>
337  Histo2DAction<HistoType>::~Histo2DAction() {  Histo2DAction<HistoType>::~Histo2DAction() {
338    
339      delete _rootHisto;
340      _rootHisto = NULL;
341  }  }
342    
343  template<class HistoType>  template<class HistoType>
# Line 444  void Histo2DAction<HistoType>::SetYAxis( Line 468  void Histo2DAction<HistoType>::SetYAxis(
468  }  }
469    
470  template<class HistoType>  template<class HistoType>
471  inline void Histo2DAction<HistoType>::Fill(double xValue, double yValue, double weight) {  void Histo2DAction<HistoType>::Fill(double xValue, double yValue, double weight) {
472    
473    _rootHisto->Fill(xValue, yValue, weight);    _rootHisto->Fill(xValue, yValue, weight);
474    
# Line 457  inline void Histo2DAction<HistoType>::Fi Line 481  inline void Histo2DAction<HistoType>::Fi
481      UOflow = true;      UOflow = true;
482    }    }
483    
484    if (xValue > _xBins.back()) {    if (xValue >= _xBins.back()) {
485      xBin = _xBins.size();      xBin = _xBins.size();
486      UOflow = true;      UOflow = true;
487    }    }
# Line 467  inline void Histo2DAction<HistoType>::Fi Line 491  inline void Histo2DAction<HistoType>::Fi
491      UOflow = true;      UOflow = true;
492    }    }
493    
494    if (yValue > _yBins.back()) {    if (yValue >= _yBins.back()) {
495      yBin = _yBins.size();      yBin = _yBins.size();
496      UOflow = true;      UOflow = true;
497    }    }
# Line 501  inline void Histo2DAction<HistoType>::Fi Line 525  inline void Histo2DAction<HistoType>::Fi
525          return;          return;
526        }        }
527        else {        else {
528          if (yBin == (int)_yBins.size()) {          if (yBin == (int) _yBins.size()) {
529            _xUnderYOverflow += (HistoType) weight;            _xUnderYOverflow += (HistoType) weight;
530            return;            return;
531          }          }
# Line 512  inline void Histo2DAction<HistoType>::Fi Line 536  inline void Histo2DAction<HistoType>::Fi
536        }        }
537      }      }
538    
539      if (xBin == (int)_xBins.size()) {      if (xBin == (int) _xBins.size()) {
540        if (yBin == -1) {        if (yBin == -1) {
541          _xOverYUnderflow += (HistoType) weight;          _xOverYUnderflow += (HistoType) weight;
542          return;          return;
543        }        }
544        else {        else {
545          if (yBin ==(int) _yBins.size()) {          if (yBin == (int) _yBins.size()) {
546            _xOverYOverflow += (HistoType) weight;            _xOverYOverflow += (HistoType) weight;
547            return;            return;
548          }          }
# Line 534  inline void Histo2DAction<HistoType>::Fi Line 558  inline void Histo2DAction<HistoType>::Fi
558        return;        return;
559      }      }
560    
561      if (yBin == (int)_yBins.size()) {      if (yBin == (int) _yBins.size()) {
562        _yOverflow[xBin] += (HistoType) weight;        _yOverflow[xBin] += (HistoType) weight;
563        return;        return;
564      }      }
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23