--- PamCut/CollectionActions/RigFillAction/RigFillAction.h 2009/08/05 17:05:01 1.2 +++ PamCut/CollectionActions/RigFillAction/RigFillAction.h 2009/08/26 16:34:15 1.3 @@ -16,16 +16,15 @@ /*! @brief The rigidity vs threshold rigidity histogram filling. * - * This class builds a 2D histogram binned in rigidity modulus + * This class builds a 2D histogram binned in event rigidity modulus * and threshold rigidity. Each 2D bin will contain the number of events whose - * rigidity modulus and threshold rigidity (eg., Stoermer cutoff rigidity times a threshold coefficient) - * lie in that bin. Note that the meaning of this threshold coefficient is the same as in + * event rigidity modulus and threshold rigidity (eg., Stoermer cutoff rigidity times a threshold coefficient) lie in that bin. Note that the meaning of this threshold coefficient is the same as in * TrkRigGeoCut, so it must have the same value used for TrkRigGeoCut (consider * using TrkRigGeoCut::GetThresholdCoeff() to retrieve its value). - * Events whose critical rigidity is below the lower limit of the axis are recorded in a vector of - * "zero bins", which range from zero to the the lowest critical rigidity. - * Events are discarded only if their rigidity modulus or cutoff rigidity lies - * outside the histogram bounds. + * Events whose threshold rigidity is below the lower limit of the bins are recorded in a vector of + * event rigidities with threshold rigidity in the "inferior threshold bin", which ranges from zero to the the lowest threshold rigidity of the "normal" bins. + * Events are discarded only if their rigidity modulus or threshold rigidity lies + * outside the histogram boundaries. * * CUT DEPENDECIES: TrkPhSinCut for single physical track, TrkRigGeoCut for galactic event. * @@ -72,8 +71,7 @@ /*! @brief Fills histogram with the selected event. * - * The current event will be added to the bin corresponding to its rigidity and critical rigidity multiplied - * by the threshold coefficient, eg., to the bin (Rc*threshold, R). + * 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. * * @param event The selected event. */ @@ -89,18 +87,22 @@ /*! @brief Returns the histogram. * - * This method returns a SimpleMatrix. Its [i][j] element contain the number of analyzed events - * whose rigidity modulus falls in the i-th rigidity bin and whose threshold rigidity multiplied by the - * threshold falls in the j-th rigidity bin. + * This method returns a SimpleMatrix. Its [i][j] 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 j-th normal bin. * * @return The rigidity modulus - threshold rigidity 2D histogram. */ - SimpleMatrix& GetHisto() { + SimpleMatrix& GetHistoThreshNormBins() { return _textHisto; } - /*! @brief Returns the vector of cutoff rigidity zero bins */ - vector& GetCutoffZeroBins() { + /*! @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 */ + vector& GetTotalHistoThreshNormBins() { + return _totalTextHisto; + } + + /*! @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 */ + vector& GetHistoThreshInfBin() { return _zeroCutoffBins; } @@ -111,6 +113,7 @@ TH2I _rootHisto; SimpleMatrix _textHisto; vector _zeroCutoffBins; + vector _totalTextHisto; float _thresholdCoeff; void _InitHistos(vector &bins);