00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015 #ifndef EvtHis2F_HH
00016 #define EvtHis2F_HH
00017
00018 #include "TH1.h"
00019 #include "TAxis.h"
00020 #include "TH2.h"
00021 #include "TH3.h"
00022 #include "TFile.h"
00023 #include "TApplication.h"
00024 #include "TROOT.h"
00025 #include "TChain.h"
00026 #include "TDirectory.h"
00027 #include "EvtGenBase/EvtRandom.hh"
00028
00029 class EvtHis2F {
00030
00031 public:
00032
00033 EvtHis2F(){}
00034 virtual ~EvtHis2F();
00035 void init();
00036 void setFile(const char* dtfile);
00037 void setHTitle(const char* htitle);
00038 const char* getFile();
00039 const char* getHTitle();
00040 void HFill(double xmass2, double ymass2);
00041 double getZvalue(double m12_square, double m13_square);
00042 void HReweight();
00043 bool AR(double xmass2,double ymass2);
00044 bool AR(double xmass2,double ymass2, double zmax, TH2F *h2);
00045 void setZmax( );
00046 void setZmax(TH2F *H2);
00047 double getZmax();
00048 void setHmc(TH2F *H2);
00049 void setHdt(TH2F *H2);
00050 void setHwt(TH2F *H2){ HWT=(TH2F*) H2->Clone("HWT");}
00051
00052 TH2F* getHmc( );
00053 TH2F* getHdt( );
00054 TH2F* getHwt( );
00055
00056 void setBINSx(int bx){BINSx=bx;}
00057 void setBINSy(int by){BINSy=by;}
00058 void setXlow (double xl){xlow=xl;}
00059 void setXup(double xu){xup=xu;}
00060 void setYlow (double yl){ylow=yl;}
00061 void setYup(double yu){yup=yu;}
00062
00063 void setBins(TH2F *h2);
00064 int getBINSx(){return BINSx;}
00065 int getBINSy(){return BINSy;}
00066 double getXlow(){return xlow;}
00067 double getYlow(){return ylow;}
00068 double getXup(){return xup;}
00069 double getYup(){return yup;}
00070
00071 void show(TH2F *h2);
00072 void showFrame(TH2F *h2);
00073
00074 void init(TH2F *hmc, TH2F *hdt, TH2F *hwt);
00075
00076 private:
00077 TH2F *HDATA, *HMC, *HWT;
00078 TFile *dataf;
00079 double zmax;
00080 int BINSx, BINSy;
00081 double xlow,xup,ylow,yup;
00082 const char *datafile, *datatitle;
00083 TDirectory dir;
00084 int icount;
00085 double max1,max2,max3;
00086 TH2F *WT1,*WT2,*WT3;
00087 };
00088
00089
00090 #endif
00091