// SimpleCFTDetector.cpp

#include "SimpleCFTDetector.h"
#include "trfcyl/LayerCylinder.h"
#include "trfcyl/ClusterFindCyl.h"
#include "trfbase/Miss.h"
#include "trfutil/Tuple.h"

using std::string;

using namespace trf;

//**********************************************************************

// constructor

SimpleCFTDetector::SimpleCFTDetector() {

  // Name layers.
  string pname[] = { "CFTX1", "CFTX2", "CFTX3", "CFTX4",
                     "CFTX5", "CFTX6", "CFTX7", "CFTX8" };
  string sname[] = { "CFTU1", "CFTV1", "CFTU2", "CFTV2",
                     "CFTU3", "CFTV3", "CFTU4", "CFTV4" };

  // define radii for axial measurements
  int nbarrel = 8;
  double rcyl[]={19.50, 23.41, 28.09, 32.77, 37.46, 42.14, 48.75, 51.50};
  double zlen[]={166.0, 166.0, 252.0, 252.0, 252.0, 252.0, 252.0, 252.0};

  // phi msmts
  int ib;
  for ( ib=0; ib<nbarrel; ++ib ) {
    double r = rcyl[ib];
    double zmax = zlen[ib]/2.0;
    double zmin = -zmax;
    MutableClusterFindManagerPtr pfind( 
      new ClusterFindCyl( r, 0.0, 11+ib ) );
    LayerPtr plyr;
    plyr = new LayerCylinder( r, zmin, zmax, pfind );
    add_layer( pname[ib], plyr );
  }

  // stereo msmts
  double stereo = -0.001;
  for ( ib=0; ib<nbarrel; ++ib ) {
    stereo *= -1.0;
    double r = rcyl[ib] + 0.2;
    double zmax = zlen[ib]/2.0;
    double zmin = -zmax;
    MutableClusterFindManagerPtr pfind(
      new ClusterFindCyl( r, stereo, 21+ib ) );
    LayerPtr plyr;
    plyr = new LayerCylinder( r, zmin, zmax, pfind );
    add_layer( sname[ib], plyr );
  }

}
  
//**********************************************************************

// destructor

SimpleCFTDetector::~SimpleCFTDetector() { }

//**********************************************************************