00001 #ifndef RAD_TKinemCut 00002 #define RAD_TKinemCut 00003 00004 #include "TRadGlobal.h" 00005 #include "TConstants.h" 00006 #include <iostream> 00007 00008 class TKinemCut{ 00009 protected: 00010 double fDeltaTheta; 00011 double fDeltaPhi; 00012 double fPAverage; 00013 double fPCross; 00014 double fEMin; 00015 double fPAverageAbs; 00016 double fPCrossAbs; 00017 double fEMinAbs; 00018 double fCosPsi; 00019 double fAverageThetaMin; 00020 double fAverageThetaMax; 00021 double fThetaMin; 00022 double fThetaMax; 00023 double fThetaMinMinus; 00024 double fThetaMaxMinus; 00025 double fThetaMinPlus; 00026 double fThetaMaxPlus; 00027 public: 00028 TKinemCut(); 00029 ~TKinemCut(){}; 00030 double DeltaTheta(){return fDeltaTheta;} 00031 double DeltaPhi(){return fDeltaPhi;} 00032 double PAverage(){return fPAverage;} 00033 double PCross(){return fPCross;} 00034 double EMin(){return fEMin;} 00035 double CosPsi(){return fCosPsi;} 00036 double AverageThetaMin(){return fAverageThetaMin;} 00037 double AverageThetaMax(){return fAverageThetaMax;} 00038 double ThetaMin(){return fThetaMin;} 00039 double ThetaMax(){return fThetaMax;} 00040 double ThetaMinM(){return fThetaMinMinus;} 00041 double ThetaMaxM(){return fThetaMaxMinus;} 00042 double ThetaMinP(){return fThetaMinPlus;} 00043 double ThetaMaxP(){return fThetaMaxPlus;} 00044 00045 void DeltaTheta(const double& x){ fDeltaTheta = x;} 00046 void DeltaPhi(const double& x){ fDeltaPhi = x;} 00047 void PAverage(const double& x){ fPAverageAbs = x;} 00048 void PCross(const double& x){ fPCrossAbs = x;} 00049 void EMin(const double& x){ fEMinAbs = x;} 00050 void CosPsi(const double& x){ fCosPsi = x;} 00051 void AverageTheta(const double& x){ 00052 fAverageThetaMin = x; 00053 fAverageThetaMax = gConst->Pi()-x; 00054 } 00055 void AverageThetaMin(const double& x){ fAverageThetaMin = x;} 00056 void AverageThetaMax(const double& x){ fAverageThetaMax = x;} 00057 void ThetaMinM(const double &x){fThetaMinMinus = x; fThetaMaxMinus = gConst->Pi()-x;} 00058 void ThetaMinP(const double &x){fThetaMinPlus = x; fThetaMaxPlus = gConst->Pi()-x;} 00059 void ThetaRangeM(const double &th_min, const double &th_max){ 00060 fThetaMinMinus = th_min; 00061 fThetaMaxMinus = th_max; 00062 } 00063 void ThetaRangeP(const double &th_min, const double &th_max){ 00064 fThetaMinPlus = th_min; 00065 fThetaMaxPlus = th_max; 00066 } 00067 00068 void Init(); 00069 void Default(); 00070 void Print(); 00071 }; 00072 00073 extern TKinemCut *gCut; 00074 00075 00076 #endif // ifndef RAD_TKinemCut