////////////////////////////////////////////////////////////////////////
//
//  JSFLTKCLTrack
//
//  A class to save linked track-cluster track
//  Algorithm and performance of Link-Track-Cluster routines are described
//  in  New Combined Track Bank(PS 329KB)  .
// 
//  Data member in this classes are those created by the subroutine CMBTRK. 
//  Bank data can be obtained by member functions of this class.
// 
//(Warning) 
//  When this class is created when reading the SIMDST data, some information
//  such as fEMC is not given properly as no corresponding information is 
//  available in SIMDST format.  Please use them with care.
//$Id: JSFLTKCLTrack.cxx,v 1.16 2003/05/23 07:04:35 miyamoto Exp $
//
////////////////////////////////////////////////////////////////////////

#include "JSFConfig.h"
#include "JSFSteer.h"
#include "JSFLTKCLTrack.h"


ClassImp(JSFLTKCLTrack)

  const Char_t *gTypeName[14]={"illegal track type",
               "pure gamma", "gamma in mixed EMC cluster",
	       "pure nutral hadron", "hadron in mixed HDC cluster",
	       "pure charged hadron", "unmatched track",
	       "illegal track type",  "illegal track type",
	       "illegal track type",  "illegal track type",
	       "electron candidate",  "illegal track type",
	       "muon candidate"};

//_____________________________________________________________________________
 JSFLTKCLTrack::JSFLTKCLTrack()
{
  fCDC=NULL;
  for(Int_t i=0;i<4;i++){ fP[i]=0.0; }
  fEcl=0.0; fNEMC=0; fCharge=0 ; fType=0;
  fSource=0; fNCDC=0 ; f1stCDC=0 ;
  fCDCs=NULL;
  fEMGen=NULL;
  fIDCDC=NULL;
  fIDEMGen=NULL;
  fNEMGen=0;

  fVTXDR=0.0;
  fVTXDZ=0.0;
  fVTXDDR=1.E10;
  fVTXDDZ=1.E10;
  fVTXNSig=0.0;

}
//_____________________________________________________________________________
 JSFLTKCLTrack::~JSFLTKCLTrack()
{
  // Delete opbjects

  if( fCDCs ) { 
    fCDCs->Clear();
    delete fCDCs;
  }
  if( fEMGen ) { 
    fEMGen->Clear();
    delete fEMGen;
  }
  if( fIDCDC ) { delete fIDCDC; }
  if( fIDEMGen ) { delete fIDEMGen; }
}


//_____________________________________________________________________________
 JSFLTKCLTrack::JSFLTKCLTrack(EJSFLTKCLTrackBank bank, Float_t data[])
{
    fBank=bank;
    fP[1]=data[0] ;  fP[2]=data[1] ;  fP[3]=data[2] ;  fP[0]=data[3] ; 
    fEcl=data[4]  ;  
    fNEMC=(Int_t)data[6]; 
    fCharge=(Int_t)data[8]; fType=(Int_t)data[9];
    fSource=(Int_t)data[10]; fNCDC=(Int_t)data[11];
    if( fNCDC > 0 ) { 
      f1stCDC=(Int_t)data[12];
      fIDCDC=new Int_t[fNCDC];
      for(Int_t i=0;i<fNCDC;i++){ fIDCDC[i]=(Int_t)data[12+i]; }
    }
    else {
      f1stCDC=-1 ;
      fIDCDC=NULL;
    }
    fIDEMGen=NULL;
    fNEMGen=0;
    fCDC=NULL;
    fCDCs=new TObjArray();
    fEMGen=new TObjArray();

    fVTXDR=0.0;
    fVTXDZ=0.0;
    fVTXDDR=1.E10;
    fVTXDDZ=1.E10;
    fVTXNSig=0.0;
}

//_____________________________________________________________________________
 JSFLTKCLTrack::JSFLTKCLTrack(Float_t data[])
{
  // Load data obtained from SIMDST data.
    fP[1]=data[0] ;  fP[2]=data[1] ;  fP[3]=data[2] ;  fP[0]=data[3] ; 
    fCharge=(Int_t)data[4]; fType=(Int_t)data[5];
    fCDC=NULL;
    fSource=0;
    fNCDC=0;
    f1stCDC=-1;
    fEcl=0;
    fNEMC=0;
    if( fType !=1 || fType != 3 ) f1stCDC=(Int_t)data[7];
    fCDCs=new TObjArray();
    fEMGen=new TObjArray();
    fIDCDC=NULL;
    fIDEMGen=NULL;
    fNEMGen=0;

    fVTXDR=0.0;
    fVTXDZ=0.0;
    fVTXDDR=1.E10;
    fVTXDDZ=1.E10;
    fVTXNSig=0.0;
}
//_____________________________________________________________________________
 JSFLTKCLTrack::JSFLTKCLTrack(JSFLTKCLTrack& t)
{
  // Load data obtained from SIMDST data.

  fBank=t.fBank;
  for(Int_t i=0;i<4;i++){ fP[i]=t.fP[i]; }
  fEcl=t.fEcl;  fNEMC=t.fNEMC;
  fCharge=t.fCharge;  fType=t.fType;
  fSource=t.fSource;   fNCDC=t.fNCDC;   f1stCDC=t.f1stCDC;
  fCDC=t.fCDC;
  if( fNCDC > 0 ) {
    fIDCDC=new Int_t[fNCDC];
    for(Int_t i=0;i<t.fNCDC;i++){ fIDCDC[i]=t.fIDCDC[i]; }
    fCDCs=new TObjArray();
    for(Int_t i=0;i<t.fCDCs->GetEntries();i++){
      fCDCs->Add(t.fCDCs->UncheckedAt(i));
    }
  }
  else {
    fIDCDC=NULL;
    fCDCs=NULL;
  }
  fEMGen=NULL;
  if( t.fEMGen != NULL ) {
    fEMGen=new TObjArray();
    fEMGen->SetOwner(kFALSE);
    for(Int_t i=0;i<t.fEMGen->GetEntries();i++){
      JSFGeneratorParticle *gp=(JSFGeneratorParticle*)t.fEMGen->UncheckedAt(i);
      fEMGen->Add(gp);
    }
  }
  fNEMGen=t.fNEMGen;
  if( fNEMGen > 0 ) {
    fIDEMGen=new Int_t[fNEMGen];
    for(Int_t i=0;i<fNEMGen;i++){ fIDEMGen[i]=t.fIDEMGen[i]; }
  }
  else {
    fIDEMGen=NULL;
    fNEMGen=0;
  }
  fVTXDR=t.fVTXDR;;
  fVTXDZ=t.fVTXDZ;
  fVTXDDR=t.fVTXDDR;
  fVTXDDZ=t.fVTXDDZ;
  fVTXNSig=t.fVTXNSig;
}

//_____________________________________________________________________________
 const Char_t *JSFLTKCLTrack::GetTypeName()
{
  // Return Track type by name
  return gTypeName[fType];
}

#if ROOT_VERSION_CODE >= ROOT_VERSION(3,0,0) && JSF_VERSION_CODE >= JSF_VERSION(1,17,26) && JSF_VERSION_CODE < JSF_VERSION(1,17,30)

   
//______________________________________________________________________________
 void JSFLTKCLTrack::Streamer(TBuffer &R__b)
{
   // Stream an object of class JSFLTKCLTrack.

   if (R__b.IsReading()) {
     UInt_t R__s, R__c;
     Version_t R__v = R__b.ReadVersion(&R__s, &R__c);
     if( R__v > 1 ) {
       JSFLTKCLTrack::Class()->ReadBuffer(R__b, this, R__v, R__s, R__c);
       return;
     }
     TObject::Streamer(R__b);
     R__b >> (Int_t&)fBank;
     R__b.ReadStaticArray(fP);
     R__b >> fEcl;
     R__b >> fNEMC;
     R__b >> fCharge;
     R__b >> fType;
     R__b >> fSource;
     R__b >> fNCDC;
     R__b >> f1stCDC;
     R__b.CheckByteCount(R__s, R__c, JSFLTKCLTrack::IsA());

   } else {
     JSFLTKCLTrack::Class()->WriteBuffer(R__b, this);
   }
}

#endif


#if ROOT_VERSION_CODE < ROOT_VERSION(3,0,0) || JSF_VERSION_CODE < JSF_VERSION(1,17,26)

//______________________________________________________________________________
 void JSFLTKCLTrack::Streamer(TBuffer &R__b)
{
   // Stream an object of class JSFLTKCLTrack.

   if (R__b.IsReading()) {
      Version_t R__v = R__b.ReadVersion(); if (R__v) { }
      TObject::Streamer(R__b);
      R__b >> (Int_t&)fBank;
      R__b.ReadStaticArray(fP);
      R__b >> fEcl;
      R__b >> fNEMC;
      R__b >> fCharge;
      R__b >> fType;
      R__b >> fSource;
      R__b >> fNCDC;
      R__b >> f1stCDC;
   } else {
      R__b.WriteVersion(JSFLTKCLTrack::IsA());
      TObject::Streamer(R__b);
      R__b << (Int_t)fBank;
      R__b.WriteArray(fP, 4);
      R__b << fEcl;
      R__b << fNEMC;
      R__b << fCharge;
      R__b << fType;
      R__b << fSource;
      R__b << fNCDC;
      R__b << f1stCDC;
   }
}

#endif