|
FlyoDoc_2011 Pisa 2011 by GmP --- 011
|
flyopunta/src/src_ana/anak6g.h
00001 /*************************************************************************** 00002 anak6g.h - description 00003 ------------------- 00004 begin : Sat Jan 18 2003 00005 copyright : (C) 2003 by Giuseppe Pierazzini 00006 email : peppe@unipi.it 00007 ***************************************************************************/ 00008 00009 #ifndef ANAK6G_H 00010 #define ANAK6G_H 00011 00012 #include "anak4g.h" 00013 #include "gmatrix.h" 00014 00020 class AnaK6g : public AnaK4g 00021 { 00022 public: 00023 00024 AnaK6g(); 00025 ~AnaK6g(); 00026 //--- root/paw interface-see also K4gam ------- 00027 int ngm,dead6, permut_fit,iter_fit; 00028 float eg_fit[6],dist_fit,cog_fit,egamt_fit,bmax; 00029 float ws6[3], ss6[3],vss6[3],ss6_sort,vss6_sort,dist_sort; 00030 // second choice 00031 double last_chiq,dchi2; 00032 //----------------------- 00033 // stuff for fit. 00034 00035 gmatrix Bm64,Az41; 00036 gmatrix F41,R41; 00037 gmatrix De61,Lmd41; 00038 gmatrix Zm11; 00039 gmatrix Var66, InV66; 00040 gmatrix BVB; 00041 gmatrix GB; 00042 gmatrix AGBA; 00043 gmatrix Gx; 00044 gmatrix CHQ; 00045 00046 00047 // for Alan fit routine.. 00048 float Ass[3], Ae_fit[6],Ax_fit[6],Ay_fit[6],Ase_fit[6],Asx_fit[6],Asy_fit[6],Azeta,Avita, 00049 Vertex[3],Vsigsq[3],Ener,Esigsq,ChiKin; 00050 00051 //=================================== 00052 // 00053 int iter_flag,permut_flag; 00054 float xg_fit[6],yg_fit[6]; 00055 double em[6],ec[6], egtk6; 00056 double sq[6],xg[6],yg[6]; 00057 double dist_vertex; 00058 double dq[6][6]; 00059 double xg_sav[12],yg_sav[12],zg_sav[12],eg_sav[12]; 00060 int idpr_fit[12],idpr_sav[12]; 00061 // double dq_fit[6][6]; 00062 00063 // 00064 // calcolo dei parametri cinematici ss.--------- 00065 qvet Gq[6], Pio01, Pio23, Pio45; 00066 gvet Versg,Pgam; 00067 //-------------------------- 00068 00069 void Reset_parm(); 00070 int Fisica(); 00071 void Permutazione ( int ); 00072 int Get_Data_from_saved ( int ); 00073 int * Get_Permuta ( int ); 00074 int Fitta(); 00075 double Calcolo_Vincoli (); 00076 double F_Vinc ( int i4,double *ec ,double dist,int id,double delta ); 00077 void Reset_Dummy_Data(); 00078 void Save_good_fit(); 00079 }; 00080 #endif
