flyopunta/src/src_ana/analisi.h

/***************************************************************************
                          analisi.h  -  description
                             -------------------
    begin                : Sat Jan 18 2003
    copyright            : (C) 2003 by Giuseppe Pierazzini
    email                : peppe@unipi.it
 ***************************************************************************/

#ifndef ANALISI_H
#define ANALISI_H
#include "flyoh.h"
#include "device.h"
#include "gigabeam.h"
using namespace std;
struct Vtxstruct
{
  int Ntr;
  gvet *P,*V;
  gvet Vtx;
  double chiq;  
  gmatrix *Proj;
  Vtxstruct(int n);
  void Vertice();
  
};



 class Analisi
  {
   protected:

   const char  *nome;
   int nev;
   int tipo;
   double ana_tempo,ana_dtim;
   double msAvo,msqAvo,ctauAvo;
   double reckik;
   GigaBeam *Pk_ric;
   
   public:
   const char *titol;
   int count_call,count_decay,getta_gg,getta_en,getta_cog,getta_tau,getta_gf;
  
   int getta_trk,getta_ch ; // counter for wrong charge 
   int getta_fit,count_wnt,getta_char;
   int getta_kabes;
   int counter[30];
   void Set_tempo_parms();
   
   
//------------------------------------------
//dummy variable for calculation...

   int idpr[12],mishit,nvar,nfit,nextra,Kmax;
   double x[12],y[12],e[12],xgood[12],ygood[12],egood[12],ek;
   double distq,dist,dsumq, diffq;
   double smaq,smbq,smcq;
   double egij[15],exij[15],eyij[15],dxij[15],dyij[15],ddij[15],mij[15];
   double dersq[24];
   double sv[24];
   gvet Pa,Pb;
   qvet q0,q1,q2,qs;
   int g0,g1,g2,g3,g4,g5,g6,flag[12];   // gi gamma index in the lkr
//--------------------------------------------

   Analisi *up,*next;
   Analisi();
   virtual ~Analisi();
   virtual int Fisica(){return 0;};
   virtual void print_scale();
   virtual int kappa(){return 0;};
   virtual void Reset_parm(){};
   virtual int MassFit(){return 0;};
   virtual int MassChiq(){return 0;};
   virtual double Distvtx(){return 0.;};
   virtual double Degij(int,int) {return 0.;};
   virtual void Ddistq(){};
   virtual void Simq(int, int){};
   virtual float Get_chi(){return 0.;};
   void Write_titol();
   int Get_Tipo() {return tipo;};
   const char *Nome(){return nome;};
   double Vertice2(gvet &Xvert,gvet P1,gvet V1,gvet P2,gvet V2);
   
   
};



#if defined(MAIN)
  Analisi    *Primofit=0;
#endif
#if !defined(MAIN)
 extern Analisi *Primofit;
#endif

#if defined(MAIN)
int Ausr=0;
int K4g=0,K6g=0,K2p=0,K3p=0,K3l=0,Ka2pi=0,Pi2g=0,E2g=0,E6g=0,K2p2g=0,Bea=0;
int Kl2=0,Ke2g=0,Kpc2g=0,Kpcpi0=0,Kpcp0g=0,Kpc2p0=0,K2p0enu=0,K2pcenu=0,Kpc3g=0,Kpivv=0,Kmunu=0,Kpovv=0;
int Kpip0=0,Kpip0ee=0;
int Klambda=0,Kevg=0,Kpievg=0,Kmvg=0;
#endif


#if !defined(MAIN)
extern int Ausr;
extern int K4g,K6g,K2p,K3p,K3l,Ka2pi,Pi2g,E2g,E6g,K2p2g,Bea;
extern int Kl2,Ke2g,Kpc2g,Kpcpi0,Kpcp0g,Kpc2p0,K2p0enu,K2pcenu,Kpc3g,Kpivv,Kmunu,Kpovv,Kpip0,Kpip0ee;
extern int Klambda,Kevg,Kpievg,Kmvg;
#endif
#endif