//////////////////////////////////////////////////////////////////////////
//
// File: TMBIsoTrks.cpp
//                                                                      //
// TMBTree class for isolated charged particles
//
// Created: 4-APR-2002     S. Protopopescu
//                                                                      //
//////////////////////////////////////////////////////////////////////////

#include "tmb_tree/TMBIsoTrks.hpp"
#include "tmb_tree/TMBTrks.hpp"

ClassImp(TMBIsoTrks)


//_____________________________________________________________________________

 TMBIsoTrks::TMBIsoTrks():
      _pshwrdphi(999.9), _pshwrdeta(999.9), 
      _iso_e(999.9), _iso_n(999.9), _svx(999.9), _svy(999.9), _svz(999.9),
      _PSphi(999.9), _PSz(999.9), _PSR(999.9), _e_mtc_tot(0.), _i_mtc(0), _tref(),
      _CPS(), _FPS(), _bestcps_match_chi2(999.9), _cps_match_chi2()
{    
    for(int j=0; j<14;j++)
      _e_mtc[j] = 0;
    for(int j=0; j<3; j++)
      _e33[j] = 0;
    for(int j=0; j<3; j++)
      _e55[j] = 0;
}



 TMBIsoTrks::TMBIsoTrks(Float_t pshwrdphi, Float_t pshwrdeta,
             Float_t iso_e, Float_t iso_n, Float_t svx, Float_t svy, Float_t svz,
             Float_t PSphi, Float_t PSz, Float_t PSR, Float_t e_mtc_tot,
		       Float_t* e_mtc, Float_t* e33, Float_t* e55, Int_t i_mtc, TRef& tref,
		       TRefArray* CPS, TRefArray* FPS, 
		       Float_t bestcps_match_chi2, TArrayF* cps_match_chi2)
{
       _pshwrdphi          = pshwrdphi;
      _pshwrdeta           = pshwrdeta;
       _iso_e              = iso_e; 
       _iso_n              = iso_n;
       _svx                = svx; 
       _svy                = svy; 
       _svz                = svz;
       _PSphi              = PSphi; 
       _PSz                = PSz; 
       _PSR                = PSR;
       _e_mtc_tot          = e_mtc_tot;

       for (Int_t i=0; i<14; i++)
           _e_mtc[i]      = e_mtc[i];
       for (Int_t i=0; i<3; i++)
           _e33[i]      = e33[i];
       for (Int_t i=0; i<3; i++)
           _e55[i]      = e55[i];

       _i_mtc=i_mtc;
       _tref=tref;
       _CPS.Clear();
       if (CPS) {
	 for (int i=0; i<CPS->GetLast()+1; i++) {
	   _CPS.Add( CPS->At(i) );
	 }
       }
       _FPS.Clear();
       if (FPS) {
	 for (int i=0; i<FPS->GetLast()+1; i++) {
	   _FPS.Add( FPS->At(i) );
	 }
       }
       _bestcps_match_chi2 = bestcps_match_chi2;
       if ( cps_match_chi2 ) _cps_match_chi2 = *cps_match_chi2;
}



//______________________________________________________________________________