#ifndef _RVERTEX_H_
#define _RVERTEX_H_

class Rvertex {
 public:
  float rver[3];
  float mass2;
  float prob;  
  float p[2];
  int hsi[2];
  int hin[2];
  int hou[2];
  int hri[2];
  int hca[2];
  int htr[2];
  int hmu[2];
  int fit[2];
  float chi[2];
  float lre[2];
  float lrmu[2];
  float lrpi[2];
  float lrk[2];
  float lrp[2];
  int flr[2];
  int otopp[2];
  float lte[2];
  float lth[2];
  int flt[2];
  float lee[2];
  float lemi[2];
  float leh[2];
  int fle[2];
  float lmmu[2];
  float lmh[2];
  int flm[2];   
  Rvertex(){ };       
  Rvertex( float *x){ for (int i=0;i<3;i++) rver[i]=x[i]; };
  Rvertex( RecoVertex *v,TrackBase* t1, TrackBase* t2);
  ~Rvertex(){};
  int Write( FILE *fp);
};     


inline Rvertex::Rvertex( RecoVertex *v,TrackBase* t1, TrackBase* t2 ){
  rver[0]= v->point().X();
  rver[1]= v->point().Y();
  rver[2]= v->point().Z();
  mass2= ( t1->GetLorentzVector() + t2->GetLorentzVector() ).Mag();;
  prob   = v->prob(); 
  p[0]   = t1->getRecoTrack()->p();
  p[1]   = t2->getRecoTrack()->p();;
  hsi[0] = t1->getRecoTrack()->GetHsi();
  hin[0] = t1->getRecoTrack()->GetHin();
  hou[0] = t1->getRecoTrack()->GetHou();
  hri[0] = t1->getRecoTrack()->GetHri();
  hca[0] = t1->getRecoTrack()->GetHca();
  htr[0] = t1->getRecoTrack()->GetHtr();
  hmu[0] = t1->getRecoTrack()->GetHmu();
  fit[0] = t1->getRecoTrack()->GetFit();
  chi[0] = t1->getRecoTrack()->GetChi2();
  lre[0] = t1->getRecoTrack()->GetLre();
  lrmu[0] = t1->getRecoTrack()->GetLrmu();
  lrpi[0] = t1->getRecoTrack()->GetLrpi();
  lrk[0] = t1->getRecoTrack()->GetLrk();
  lrp[0] = t1->getRecoTrack()->GetLrp();
  flr[0] = t1->getRecoTrack()->GetFlr();
  lte[0] = t1->getRecoTrack()->GetLte();
  lth[0] = t1->getRecoTrack()->GetLth();
  flt[0] = t1->getRecoTrack()->GetFlt();
  lee[0] = t1->getRecoTrack()->GetLee();
  lemi[0] = t1->getRecoTrack()->GetLemi();
  leh[0] = t1->getRecoTrack()->GetLeh();
  fle[0] = t1->getRecoTrack()->GetFle();
  lmmu[0] = t1->getRecoTrack()->GetLmmu();
  lmh[0] = t1->getRecoTrack()->GetLmh();
  flm[0] = t1->getRecoTrack()->GetFlm();       


  hsi[1] = t2->getRecoTrack()->GetHsi();
  hin[1] = t2->getRecoTrack()->GetHin();
  hou[1] = t2->getRecoTrack()->GetHou();
  hri[1] = t2->getRecoTrack()->GetHri();
  hca[1] = t2->getRecoTrack()->GetHca();
  htr[1] = t2->getRecoTrack()->GetHtr();
  hmu[1] = t2->getRecoTrack()->GetHmu();
  fit[1] = t2->getRecoTrack()->GetFit();
  chi[1] = t2->getRecoTrack()->GetChi2();
  lre[1] = t2->getRecoTrack()->GetLre();
  lrmu[1] = t2->getRecoTrack()->GetLrmu();
  lrpi[1] = t2->getRecoTrack()->GetLrpi();
  lrk[1] = t2->getRecoTrack()->GetLrk();
  lrp[1] = t2->getRecoTrack()->GetLrp();
  flr[1] = t2->getRecoTrack()->GetFlr();
  lte[1] = t2->getRecoTrack()->GetLte();
  lth[1] = t2->getRecoTrack()->GetLth();
  flt[1] = t2->getRecoTrack()->GetFlt();
  lee[1] = t2->getRecoTrack()->GetLee();
  lemi[1] = t2->getRecoTrack()->GetLemi();
  leh[1] = t2->getRecoTrack()->GetLeh();
  fle[1] = t2->getRecoTrack()->GetFle();
  lmmu[1] = t2->getRecoTrack()->GetLmmu();
  lmh[1] = t2->getRecoTrack()->GetLmh();
  flm[1] = t2->getRecoTrack()->GetFlm(); 

};

inline int WriteRvertex(Rvertex *r,FILE *fp){ return fwrite(r,sizeof(Rvertex),1,fp); }
inline int ReadRvertex(Rvertex *r,FILE *fp){ return fread(r,sizeof(Rvertex)-4,1,fp); }

#endif