Reconstruction/MdcPatRec/MdcGeom/test/Test.cxx

Go to the documentation of this file.

#include "CLHEP/Geometry/Point3D.h"
#include "MdcGeom/MdcDetector.h"
#include "MdcGeom/MdcSWire.h"
#include "MdcGeom/MdcSagTraj.h"
#include "MdcGeomSvc/MdcGeomSvc.h"
#include <iostream>
using namespace std;
 
int main() {
  /*   const std::string& name = " ";
     ISvcLocator* svcloc = NULL ;
     MdcGeomSvc fGeom(name, svcloc);
     fGeom.initialize();
     cout <<"dump"<<endl;
     fGeom.Dump();
     cout<<"Hi, event routine is running"<<endl;
     cout<<"1st wire id:"<<fGeom.Wire(0)->Id()<<endl;
     cout<<"2nd wire lyr id:"<<fGeom.Wire(1)->Lyr()->Id()<<endl;
     cout<<"6860th wire lyr id:"<<fGeom.Wire(6859)->Lyr()->Id()<<endl;
     fGeom.finalize();
  */
  MdcDetector det();
  //  HepPoint3D p1(78.906,6.21,393.);
  //  HepPoint3D p2(76.9631,-18.4772,-393.);
  //  HepPoint3D p2(76.906,6.21,-393.);
  HepPoint3D        p1( 0.0, 0.0, 393. );
  HepPoint3D        p2( 0.0, 0.0, -393. );
  MdcSWire          wire( p1, p2, 0.5 );
  const MdcSagTraj* trr = wire.getTraj();
  MdcSagTraj        tr( 0.5, p1, p2 );
  cout << tr.sag() << endl;
  cout << tr.rawDirection() << endl;
  cout << tr.position( 300.0 ) << endl;
  cout << trr->sag() << endl;
  cout << trr->rawDirection() << endl;
  cout << trr->position( 300.0 ) << endl;
}