/[PAMELA software]/PamCut/TrkCuts/TrkDedxHCut/TrkDedxHCut.h
ViewVC logotype

Diff of /PamCut/TrkCuts/TrkDedxHCut/TrkDedxHCut.h

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

revision 1.1 by pam-fi, Wed May 27 13:30:08 2009 UTC revision 1.2 by pam-fi, Thu Mar 18 14:46:11 2010 UTC
# Line 1  Line 1 
1  /*  /*
2   * TrkDedxHCut.h   * TrkDedxHCut.h
3   *   *
4   *  Created on: 26-mar-2009   *  Created on: 15-mar-2010
5   *      Author: Nico De Simone, S. Ricciarini   *      Author: S. Ricciarini
6   */   */
7    
8  /*! @file TrkDedxHCut.h The TrkDedxHCut class definition file */  /*! @file TrkDedxHCut.h The TrkDedxHCut class definition file */
# Line 12  Line 12 
12    
13  #include "../../PamCutBase/PamCutBase.h"  #include "../../PamCutBase/PamCutBase.h"
14    
 #include <TH2F.h>  
   
15  /*! @brief The dE/dX cut to select H and anti-H.  /*! @brief The dE/dX cut to select H and anti-H.
16   * 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).   * 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).
17   * 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.   * 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.
18   * CUT DEPENDENCIES: TrkPhysSin for object trkTrack   * CUT DEPENDENCIES: TrkPhysSin for object trkTrack
19   */   */
20  class TrkDedxHCut: public PamCut {  class TrkDedxHCut: public PamCut {
# Line 30  public: Line 28  public:
28    TrkDedxHCut(const char *cutName) :    TrkDedxHCut(const char *cutName) :
29      PamCut(cutName) {      PamCut(cutName) {
30    
 #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  
   
31    }    }
32    
33    /*! @brief Destructor. */    /*! @brief Destructor. */
34    ~TrkDedxHCut() {    ~TrkDedxHCut() {
35    
 #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  
   
36    }    }
37    
38    /*! @brief The rigidity check.    /*! @brief The rigidity check.
# Line 74  public: Line 45  public:
45    
46  private:  private:
47    
 #ifdef DEBUGPAMCUT  
   
   TH2F* h_trk_he_dedx_rigmod[2];  
   TH2F* h_trk_le_dedx_rigmod[2];  
   
 #endif  
   
48    // lower cut line    // lower cut line
49    static const int _nLow = 23;    static const int _nLow = 23;
50    static const Double_t _xLow[];    static const Double_t _xLow[];
51    static const Double_t _yLow[];    static const Double_t _yLow[];
52    
   // higher cut line  
   static const int _nHigh = 21;  
   static const Double_t _xHigh[];  
   static const Double_t _yHigh[];  
   
53  };  };
54  #endif  #endif /* TRKDEDXHCUT_H_ */
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.23