//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      Implementation of class TrackFitStat
//      see TrackFitStat.hh for details
//
// Environment:
//      Software developed for the PANDA Detector at FAIR.
//
// Author List:
//      Martin Berger    TUM            (original author)
//
//
//-----------------------------------------------------------

// Panda Headers ----------------------

// This Class' Header ------------------
#include "TpcCosmics.h"

// C/C++ Headers ----------------------
#include <iostream>


// Collaborating Class Headers --------


// Class Member definitions -----------_TpcDigiPrfWeights

ClassImp(TpcCosmics)

TpcCosmics::TpcCosmics()
    : _p(-1000),
      _pmc(-1000),
      _pstart(-1000),
      _theta(-1000),
      _phi(-1000),
      _length(-1),
      _clmeanamp(-1),
      _clmediamp(-1),
      _clsumamp(-1),
      _clsizemed2d(-1),
      _clsizemed3d(-1),
      _mom(0,0,0),
      _momerr(0,0,0),
      _pos(0,0,-1),
      _poserr(0,0,0),
      _MCpos(0,0,-1),
      _MCmom(0,0,0),
      _zAxisPoca(0,0,0),
      _zAxisDir(0,0,0),
      _firstMCPoca(0,0,0),
      _firstMCDir(0,0,0),
      _charge(0),
      _mccharge(0),
      _sigp(0),
      _nPDG(0),
      _motherID(0),
      _statflag(-100),
      _meanz(-1),
      _event(-1),
      _trackChi(-1),
      _failedHits(-1),
      _hasMc(false),
      _ChiSqu(-1),
      _NDF(-1),
      _TpcResX(),
      _TpcResY(),
      _TpcResZ(),
      _TpcResU(),
      _TpcResV(),
      _TpcResMCU(),
      _TpcResMCV(),
      _TpcErrU(),
      _TpcErrV(),
      _TpcResR(),
      _TpcResP(),
      _TpcResA(),
      _TpcResAMC(),
      _TpcSigX(),
      _TpcSigY(),
      _TpcSigZ(),
      _TpcSigC(),
      _TpcSigA(),
      _TpcSigT(),
      _TpcSigResT(),
      _TpcSigResMC(),
      _TpcAxis1(),
      _TpcAxis2(),
      _TpcAxis3(),
      _TpcAxis1MCErr(),
      _TpcAxis2MCErr(),
      _TpcAxis3MCErr(),
      _TpcAxis1Diff(),
      _TpcAxis2Diff(),
      _TpcAxis3Diff(),
      _PlaneU(),
      _PlaneV(),
      _PlaneW(),
      _PlaneO(),
      _MidPlaneU(),
      _MidPlaneV(),
      _MidPlaneW(),
      _MidPlaneO(),
      _TpcResMcX(),
      _TpcResMcY(),
      _TpcResMcZ(),
      _TpcResMcP(),
      _TpcResMcR(),
      _TpcMcRefPosX(),
      _TpcMcRefPosY(),
      _TpcMcRefPosZ(),
      _TpcMcRefPosR(),
      _TpcMcPosMom(),
      _TpcGTrack(),
      _TpcChi2X(),
      _TpcChi2Y(),
      _TpcChi2Z(),
      _TpcClusterSize(),
      _TpcClusterSizeU(),
      _Tpc2DClusterSize(),
      _TpcFullClusterSize(),
      _TpcFullClusterSizeTrack(),
      _TpcFullClusterSizeAxis(),
      _TpcClusterAmp(),
      _TpcClusterPhi(),
      _TpcProj(),
      _TpcProjErr(),
      _TpcProjMom(),
      _TpcProjMomErr(),
      _TpcPosX(),
      _TpcPosY(),
      _TpcPosZ(),
      _TpcPosR(),
      _TpcHitID(),
      _TpcClusterDist(),
      _TpcMaxDigiAmp(),
	  _TpcMaxSampleAmp(),
      _TpcDigisOnPlane(),
      _TpcDigiMCRes(),
      _TpcDigiRes(),
      _TpcDigisOnPlaneAmp(),
      _TpcDigiLen(),
      _TpcSaturated(),
      _TpcDSNR(),
      _TpcCSNR(),
	  _TpcSSNR(),
      _TpcIsGeomErr(),
      _TpcClusterFitGood(),
      _TpcEigenValOnPlane(),
      _TpcMinSigma(),
      _TpcResMCTrack(),
      _TpcResMCTrackUV(),
      _TpcResMCTrackA(),
      _TpcCogPos(),
      _TpcMinNll(),
      _TpcClusterCov(),
      _TpcOuterAmp(),
      _TpcInnerAmp(),
      _TpcDAFWeights(),
      _TpcCovOnPlane(),
      _TpcDigiId()
{}

TpcCosmics::TpcCosmics(const TpcCosmics& _c)

  : _p(_c._p),
    _pmc(_c._pmc),
    _pstart(_c._pstart),
    _theta(_c._theta),
    _phi(_c._phi),
    _length(_c._length),
    _clmeanamp(_c._clmeanamp),
    _clmediamp(_c._clmediamp),
    _clsumamp(_c._clsumamp),
    _clsizemed2d(_c._clsizemed2d),
    _clsizemed3d(_c._clsizemed3d),
    _mom(_c._mom),
    _momerr(_c._momerr),
    _pos(_c._pos),
    _poserr(_c._poserr),
    _MCpos(_c._MCpos),
    _MCmom(_c._MCmom),
    _zAxisPoca(_c._zAxisPoca),
    _zAxisDir(_c._zAxisDir),
    _firstMCPoca(_c._firstMCPoca),
    _firstMCDir(_c._firstMCDir),
    _charge(_c._charge),
    _mccharge(_c._mccharge),
    _sigp(_c._sigp),
    _nPDG(_c._nPDG),
    _motherID(_c._motherID),
    _statflag(_c._statflag),
    _meanz(_c._meanz),
    _event(_c._event),
    _trackChi(_c._trackChi),
    _failedHits(_c._failedHits),
    _hasMc(_c._hasMc),
    _ChiSqu(_c._ChiSqu),
    _NDF(_c._NDF),
    _TpcResX(_c._TpcResX),
    _TpcResY(_c._TpcResY),
    _TpcResZ(_c._TpcResZ),
    _TpcResU(_c._TpcResU),
    _TpcResV(_c._TpcResV),
    _TpcResMCU(_c._TpcResMCU),
    _TpcResMCV(_c._TpcResMCV),
    _TpcErrU(_c._TpcErrU),
    _TpcErrV(_c._TpcErrV),
    _TpcResR(_c._TpcResR),
    _TpcResP(_c._TpcResP),
    _TpcResA(_c._TpcResA),
    _TpcResAMC(_c._TpcResAMC),
    _TpcSigX(_c._TpcSigX),
    _TpcSigY(_c._TpcSigY),
    _TpcSigZ(_c._TpcSigZ),
    _TpcSigC(_c._TpcSigC),
    _TpcSigA(_c._TpcSigA),
    _TpcSigT(_c._TpcSigT),
    _TpcSigResT(_c._TpcSigResT),
    _TpcSigResMC(_c._TpcSigResMC),
    _TpcAxis1(_c._TpcAxis1),
    _TpcAxis2(_c._TpcAxis2),
    _TpcAxis3(_c._TpcAxis3),
    _TpcAxis1MCErr(_c._TpcAxis1MCErr),
    _TpcAxis2MCErr(_c._TpcAxis2MCErr),
    _TpcAxis3MCErr(_c._TpcAxis3MCErr),
    _TpcAxis1Diff(_c._TpcAxis1Diff),
    _TpcAxis2Diff(_c._TpcAxis2Diff),
    _TpcAxis3Diff(_c._TpcAxis3Diff),
    _PlaneU(_c._PlaneU),
    _PlaneV(_c._PlaneV),
    _PlaneW(_c._PlaneW),
    _PlaneO(_c._PlaneO),
    _MidPlaneU(_c._MidPlaneU),
    _MidPlaneV(_c._MidPlaneV),
    _MidPlaneW(_c._MidPlaneW),
    _MidPlaneO(_c._MidPlaneO),
    _TpcResMcX(_c._TpcResMcX),
    _TpcResMcY(_c._TpcResMcY),
    _TpcResMcZ(_c._TpcResMcZ),
    _TpcResMcP(_c._TpcResMcP),
    _TpcResMcR(_c._TpcResMcR),
    _TpcMcRefPosX(_c._TpcMcRefPosX),
    _TpcMcRefPosY(_c._TpcMcRefPosY),
    _TpcMcRefPosZ(_c._TpcMcRefPosZ),
    _TpcMcRefPosR(_c._TpcMcRefPosR),
    _TpcMcPosMom(_c._TpcMcPosMom),
    _TpcGTrack(_c._TpcGTrack),
    _TpcChi2X(_c._TpcChi2X),
    _TpcChi2Y(_c._TpcChi2Y),
    _TpcChi2Z(_c._TpcChi2Z),
    _TpcClusterSize(_c._TpcClusterSize),
    _TpcClusterSizeU(_c._TpcClusterSizeU),
    _Tpc2DClusterSize(_c._Tpc2DClusterSize),
    _TpcFullClusterSize(_c._TpcFullClusterSize),
    _TpcFullClusterSizeTrack(_c._TpcFullClusterSizeTrack),
    _TpcFullClusterSizeAxis(_c._TpcFullClusterSizeAxis),
    _TpcClusterAmp(_c._TpcClusterAmp),
    _TpcClusterPhi(_c._TpcClusterPhi),
    _TpcProj(_c._TpcProj),
    _TpcProjErr(_c._TpcProjErr),
    _TpcProjMom(_c._TpcProjMom),
    _TpcProjMomErr(_c._TpcProjMomErr),
    _TpcPosX(_c._TpcPosX),
    _TpcPosY(_c._TpcPosY),
    _TpcPosZ(_c._TpcPosZ),
    _TpcPosR(_c._TpcPosR),
    _TpcHitID(_c._TpcHitID),
    _TpcClusterDist(_c._TpcClusterDist),
    _TpcMaxDigiAmp(_c._TpcMaxDigiAmp),
	_TpcMaxSampleAmp(_c._TpcMaxSampleAmp),
    _TpcDigisOnPlane(_c._TpcDigisOnPlane),
    _TpcDigiPos(_c._TpcDigiPos),
    _TpcDigiMCPos(_c._TpcDigiMCPos),
    _TpcDigisOnPlaneAmp(_c._TpcDigisOnPlaneAmp),
    _TpcSaturated(_c._TpcSaturated),
    _TpcDSNR(_c._TpcDSNR),
    _TpcCSNR(_c._TpcCSNR),
	_TpcSSNR(_c._TpcSSNR),
    _TpcIsGeomErr(_c._TpcIsGeomErr),
    _TpcClusterFitGood(_c._TpcClusterFitGood),
    _TpcEigenValOnPlane(_c._TpcEigenValOnPlane),
    _TpcMinSigma(_c._TpcMinSigma),
    _TpcResMCTrack(_c._TpcResMCTrack),
    _TpcResMCTrackUV(_c._TpcResMCTrackUV),
    _TpcCogPos(_c._TpcCogPos),
    _TpcMinNll(_c._TpcMinNll),
    _TpcClusterCov(_c._TpcClusterCov),
    _TpcOuterAmp(_c._TpcOuterAmp),
    _TpcInnerAmp(_c._TpcInnerAmp),
	_TpcDAFWeights(_c._TpcDAFWeights),
	_TpcCovOnPlane(_c._TpcCovOnPlane),
	_TpcUnwCovOnPlane(_c._TpcUnwCovOnPlane),
	_TpcDigiMCRes(_c._TpcDigiMCRes),
	_TpcDigiTrackRes(_c._TpcDigiTrackRes),
    _TpcDigiRes(_c._TpcDigiRes),
	_TpcDigiLen(_c._TpcDigiLen),
	_TpcDigiMom(_c._TpcDigiMom),
	_TpcDigiId(_c._TpcDigiId),
	_TpcShapeCov(_c._TpcShapeCov),
	_TpcUnwShapeCov(_c._TpcUnwShapeCov),
	_TpcShapeCovOnPlane(_c._TpcShapeCovOnPlane),
	_TpcUnwShapeCovOnPlane(_c._TpcUnwShapeCovOnPlane),
	_TpcRMS(_c._TpcRMS),
	_TpcUnwRMS(_c._TpcUnwRMS),
	_TpcRMS_UV(_c._TpcRMS_UV),
	_TpcUnwRMS_UV(_c._TpcUnwRMS_UV),
	_TpcRMS_A(_c._TpcRMS_A),
	_TpcUnwRMS_A(_c._TpcUnwRMS_A),
	_TpcRMS_A_UV(_c._TpcRMS_A_UV),
	_TpcUnwRMS_A_UV(_c._TpcUnwRMS_A_UV),
	_TpcDigiPrfWeights(_c._TpcDigiPrfWeights),
	_TpcFieldCorr(_c._TpcFieldCorr)
{
//	std::cout<<"copied cosmic \n";
//	for(int i=0;i<_TpcDigiPrfWeights.size();i++)
//		std::cout<<_TpcDigiPrfWeights[i].size()<<"\n";
}

TpcCosmics::~TpcCosmics()
{
  ;
}

const std::vector<TVector3> TpcCosmics::GetResMC()
{
  std::vector<TVector3> vec;
  for(int i=0;i<_TpcResMcX.size();i++)
    {
      vec.push_back(TVector3(_TpcResMcX[i],_TpcResMcY[i],_TpcResMcZ[i]));
    }
  return vec;
}

const std::vector<TVector3> TpcCosmics::GetRes()
{
  std::vector<TVector3> vec;
  for(int i=0;i<_TpcResX.size();i++)
    vec.push_back(TVector3(_TpcResX[i],_TpcResY[i],_TpcResZ[i]));
  return vec;
}