--- PamCut/TrkCuts/TrkDedxHCut/TrkDedxHCut.h 2009/05/27 13:30:08 1.1 +++ PamCut/TrkCuts/TrkDedxHCut/TrkDedxHCut.h 2010/03/18 14:46:11 1.2 @@ -1,8 +1,8 @@ /* * TrkDedxHCut.h * - * Created on: 26-mar-2009 - * Author: Nico De Simone, S. Ricciarini + * Created on: 15-mar-2010 + * Author: S. Ricciarini */ /*! @file TrkDedxHCut.h The TrkDedxHCut class definition file */ @@ -12,11 +12,9 @@ #include "../../PamCutBase/PamCutBase.h" -#include - /*! @brief The dE/dX cut to select H and anti-H. - * This cut discards all the events whose rigidity modulus and dE/dX (both given by TRK) are outside a pre-defined band (by Nico De Simone). - * Here rigidity is defined as p/Z (GV) where Z is the particle charge (WITH SIGN) and p the momentum modulus: therefore rigidity can be positive or negative. + * This cut discards all the events whose rigidity modulus and dE/dX (both given by TRK) are outside a pre-defined band (lower band: by Nico De Simone; upper band: by Nicola Mori). + * Here the rigidity is defined as pc/Ze (GV) where Z is the particle charge (WITH SIGN) and p the momentum modulus: therefore the rigidity can be positive or negative. * CUT DEPENDENCIES: TrkPhysSin for object trkTrack */ class TrkDedxHCut: public PamCut { @@ -30,38 +28,11 @@ TrkDedxHCut(const char *cutName) : PamCut(cutName) { -#ifdef DEBUGPAMCUT - - TString hId; - TString hTitle; - - for (UInt_t j = 0; j < 2; j++) { - hId.Form("h_trk_he_dedx_vs_rigmod_%i", j); - hTitle.Form("TRK HE dE/dX vs |rig| (%i)", j); - h_trk_he_dedx_rigmod[j] = new TH2F(hId.Data(), hTitle.Data(), 50, 0, 20, 50, 0, 10); - } - for (UInt_t j = 0; j < 2; j++) { - hId.Form("h_trk_le_dedx_vs_rigmod_%i", j); - hTitle.Form("TRK LE dE/dX vs |rig| (%i)", j); - h_trk_le_dedx_rigmod[j] = new TH2F(hId.Data(), hTitle.Data(), 60, 0, 3, 50, 0, 10); - } - -#endif - } /*! @brief Destructor. */ ~TrkDedxHCut() { -#ifdef DEBUGPAMCUT - - for (UInt_t j = 0; j < 2; j++) { - h_trk_he_dedx_rigmod[j]->Write(); - h_trk_le_dedx_rigmod[j]->Write(); - } - -#endif - } /*! @brief The rigidity check. @@ -74,22 +45,10 @@ private: -#ifdef DEBUGPAMCUT - - TH2F* h_trk_he_dedx_rigmod[2]; - TH2F* h_trk_le_dedx_rigmod[2]; - -#endif - // lower cut line static const int _nLow = 23; static const Double_t _xLow[]; static const Double_t _yLow[]; - // higher cut line - static const int _nHigh = 21; - static const Double_t _xHigh[]; - static const Double_t _yHigh[]; - }; -#endif +#endif /* TRKDEDXHCUT_H_ */