/[PAMELA software]/PamCut/CollectionActions/RigFillAction/RigFillAction.h
ViewVC logotype

Diff of /PamCut/CollectionActions/RigFillAction/RigFillAction.h

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

revision 1.2 by pam-fi, Wed Aug 5 17:05:01 2009 UTC revision 1.4 by pam-fi, Tue Oct 27 10:18:21 2009 UTC
# Line 5  Line 5 
5   *      Author: Nicola Mori   *      Author: Nicola Mori
6   */   */
7    
8  /*! @file RigFillAction.h The RigFillAction class declaration file */  /*! @file RigFillAction.h The RigFillAction class declaration file. */
9    
10  #ifndef RIGFILLACTION_H_  #ifndef RIGFILLACTION_H_
11  #define RIGFILLACTION_H_  #define RIGFILLACTION_H_
# Line 16  Line 16 
16    
17  /*! @brief The rigidity vs threshold rigidity histogram filling.  /*! @brief The rigidity vs threshold rigidity histogram filling.
18   *   *
19   * This class builds a 2D histogram binned in rigidity modulus   * This class builds a 2D histogram binned in event rigidity modulus
20   * and threshold rigidity. Each 2D bin will contain the number of events whose   * and threshold rigidity. Each 2D bin will contain the number of events whose
21   * rigidity modulus and threshold rigidity (eg., Stoermer cutoff rigidity times a threshold coefficient)   * event rigidity modulus and threshold rigidity (eg., Stoermer cutoff rigidity times a
22   * lie in that bin. Note that the meaning of this threshold coefficient is the same as in   * threshold coefficient) lie in that bin. Note that the meaning of this threshold coefficient
23   * TrkRigGeoCut, so it must have the same value used for TrkRigGeoCut (consider   * is the same as in TrkRigGeoCut, so it must have the same value used for TrkRigGeoCut
24   * using TrkRigGeoCut::GetThresholdCoeff() to retrieve its value).   * (consider using TrkRigGeoCut::GetThresholdCoeff() to retrieve its value).
25   * Events whose critical rigidity is below the lower limit of the axis are recorded in a vector of   * Events whose threshold rigidity is below the lower limit of the bins are recorded in a vector
26   * "zero bins", which range from zero to the the lowest critical rigidity.   * of event rigidities with threshold rigidity in the "inferior threshold bin", which ranges
27   * Events are discarded only if their rigidity modulus or cutoff rigidity lies   * from zero to the the lowest threshold rigidity of the "normal" bins.
28   * outside the histogram bounds.   * Events are discarded only if their rigidity modulus or threshold rigidity lies
29     * outside the histogram boundaries.
30   *   *
31   * CUT DEPENDECIES: TrkPhSinCut for single physical track, TrkRigGeoCut for galactic event.   * CUT DEPENDENCIES: TrkPhSinCut for single physical track, TrkRigGeoCut for galactic event.
32   *   *
33   */   */
34  class RigFillAction: public CollectionAction {  class RigFillAction: public CollectionAction {
# Line 72  public: Line 73  public:
73    
74    /*! @brief Fills histogram with the selected event.    /*! @brief Fills histogram with the selected event.
75     *     *
76     * The current event will be added to the bin corresponding to its rigidity and critical rigidity multiplied     * The current event will be added to the event rigidity bin corresponding to its event rigidity R and to the threshold bin corresponding to the Stoermer cutoff rigidity (S) multiplied by the threshold coefficient.
    * by the threshold coefficient, eg., to the bin (Rc*threshold, R).  
77     *     *
78     * @param event The selected event.     * @param event The selected event.
79     */     */
# Line 89  public: Line 89  public:
89    
90    /*! @brief Returns the histogram.    /*! @brief Returns the histogram.
91     *     *
92     * This method returns a SimpleMatrix. Its [i][j] element contain the number of analyzed events     * This method returns a SimpleMatrix. Its [i][j] element contains the number of analyzed events
93     * whose rigidity modulus falls in the i-th rigidity bin and whose threshold rigidity multiplied by the     * whose event rigidity modulus falls in the i-th normal rigidity bin and whose threshold rigidity falls in the j-th normal bin.
    * threshold falls in the j-th rigidity bin.  
94     *     *
95     * @return The rigidity modulus - threshold rigidity 2D histogram.     * @return The rigidity modulus - threshold rigidity 2D histogram.
96     */     */
97    SimpleMatrix<UInt_t>& GetHisto() {    SimpleMatrix<UInt_t>& GetHistoThreshNormBins() {
98      return _textHisto;      return _textHisto;
99    }    }
100    
101    /*! @brief Returns the vector of cutoff rigidity zero bins */    /*! @brief whose [i] element contains the number of analyzed events whose event rigidity modulus falls in the i-th normal rigidity bin, and whose threshold rigidity falls in any normal bin */
102    vector<UInt_t>& GetCutoffZeroBins() {    vector<UInt_t>& GetTotalHistoThreshNormBins() {
103        return _totalTextHisto;
104      }
105    
106      /*! @brief Returns the vector whose [i] element contains the number of analyzed events whose event rigidity modulus falls in the i-th normal rigidity bin, and whose threshold rigidity falls in the inferior bin */
107      vector<UInt_t>& GetHistoThreshInfBin() {
108      return _zeroCutoffBins;      return _zeroCutoffBins;
109    }    }
110    
# Line 111  private: Line 115  private:
115    TH2I _rootHisto;    TH2I _rootHisto;
116    SimpleMatrix<UInt_t> _textHisto;    SimpleMatrix<UInt_t> _textHisto;
117    vector<UInt_t> _zeroCutoffBins;    vector<UInt_t> _zeroCutoffBins;
118      vector<UInt_t> _totalTextHisto;
119    float _thresholdCoeff;    float _thresholdCoeff;
120    
121    void _InitHistos(vector<float> &bins);    void _InitHistos(vector<float> &bins);
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.23