#ifndef _KinLine_H
#define _KinLine_H
// System Headers:
#include <iostream>
#include <string>
#include <vector>
#include <cmath>
#include <fstream>
// ROOT Headers:
#include "TLorentzVector.h"
#include "TVector3.h"
#include "TMatrixD.h"
// Local Headers:
//#include "ClasRuns.h"
//_____________________________________________________________________________


typedef std::string String; 
// Full covaraiance matrix from tracking one 
TMatrixD CorrectCLAS_V(const TMatrixD &__covTrack,
		      const std::vector<String> &__p_names,
		      const std::vector<TLorentzVector> &__p4,
		      const std::vector<TVector3> &__vert,
		      bool __multi,bool __isMC, String &__xid);

// functions to get tracking parameters from 4-momentum
double LambdaTrack(const TLorentzVector &__p4);
double PhiTrack(const TLorentzVector &__p4);
double AlphaTrack(const TLorentzVector &__p4);

// reads in resolution parameters
bool ReadInResParams(double __pPars[3][6][15][3],double __lamPars[3][6][15],
		     double __phiPars[3][6][15],String &__xid);
		     

//_____________________________________________________________________________

// inline functions:

/// Get CLAS sector number from 4-momentum
inline int GetSectorFromP4(const TLorentzVector &__p4){

  int sector = 0;
  double pi = 3.14159;
  double phi_lab = __p4.Phi();
  double phi = (180./pi)*phi_lab;

  if(std::abs(phi) <= 30.) sector = 1;
  else if(phi > 0.){
    if(phi <= 90.) sector = 2;
    else if(phi <= 150) sector = 3;
    else sector = 4;
  }
  else {
    // phi < 0
    if(std::abs(phi) <= 90.) sector = 6;
    else if(std::abs(phi) <= 150.) sector = 5;
    else sector = 4;
  }
  return sector;
}
//_____________________________________________________________________________

/// Returns the \f$ \theta \f$ bin (for resolution parameters)
inline int GetThetaBin(double __theta){

  double theta = __theta*180./3.14159;
  if(theta < 10.) return 0;
  else if(theta < 15) return 1;
  else if(theta < 20) return 2;
  else if(theta < 25) return 3;
  else if(theta < 30) return 4;
  else if(theta < 35) return 5;
  else if(theta < 40) return 6;
  else if(theta < 45) return 7;
  else if(theta < 50) return 8;
  else if(theta < 60) return 9;
  else if(theta < 70) return 10;
  else if(theta < 80) return 11;
  else if(theta < 90) return 12;
  else if(theta < 110) return 13;
  else return 14;
}
//_____________________________________________________________________________

/// Checks to see if resolution parameter file was read correctly
inline bool CheckRowsRead(const String &__fileName,int __rowsRead){

  if(__rowsRead != 90){ // should be 6 sectors X 15 rows/sector = 90 rows
    std::cout << "Error! <ReadInResParams> File read error: " << __fileName
	      << " Read in incorrect number of rows (" << __rowsRead 
	      << " instead of 90)" << std::endl;
    return false;
  }
  return true;
}
//_____________________________________________________________________________

#endif