Reconstruction/MdcPatRec/MdcGeom/MdcGeom-00-01-17/test/Test.cxx

Go to the documentation of this file.

#include <iostream>
#include "MdcGeomSvc/MdcGeomSvc.h"
#include "MdcGeom/MdcSagTraj.h"
#include "CLHEP/Geometry/Point3D.h"
#include "MdcGeom/MdcSWire.h"
#include "MdcGeom/MdcDetector.h"
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;
 
}