MucGeoTest1.cxx

Go to the documentation of this file.

//$id$
//
//$log$
 
/*
 *    2003/09/06   Zhengyun You      Peking University
 * 
 *    2004/09/09   Zhengyun You     Peking University
 *                 transplanted to Gaudi framework
 */
 
using namespace std;
 
#include <vector>
#include <iostream>
#include <iomanip>
#include <string>
//#include "TRint.h"
//#include <TGeoManager.h>
 
#include "MucGeomSvc/MucGeoGeneral.h"
 
int
main()
{
  //new TRint("ROOT GDML converter", 0, 0);
  //if (gGeoManager) delete gGeoManager;
  //new TGeoManager("BesGeo", "Bes geometry");
 
  cout << endl << "****** Initialize MucGeoGeneral! ******" << endl << endl;
  MucGeoGeneral::Instance()->Init();
 
    
  std::cout<<"================================ begin event " << std::endl;
 
  HepPoint3D gPoint(0,0,0);
  Hep3Vector gVect(-1,0.8,-1.4);
 
  cout << "********************" << endl << endl ;
  cout << " Origin    of the track : " << gPoint << endl;
  cout << " Direction of the track : " << gVect  << endl;
 
  vector<Identifier> idGap, idGapTemp;
  for(unsigned int part = 0; part < MucID::getPartNum(); part++) {
    for(unsigned int gap = 0; gap < MucID::getGapNum(part); gap++) {
      idGapTemp = MucGeoGeneral::Instance()->FindIntersectGaps(part, gap, gPoint, gVect);
      
      for( vector<Identifier>::iterator iter  = idGapTemp.begin(); 
       iter != idGapTemp.end(); ++iter) {
    idGap.push_back(*iter);
      }
    }
  }
  
  cout << endl << "********************" << endl << endl;
  cout << " Number of intersect gaps : " << idGap.size() << endl << endl;
  
  for( unsigned int i = 0; i< idGap.size(); i++){
    cout << idGap[i] << endl;
  }
  
  vector<Identifier> idStrip, idStripTemp;
  for(unsigned int part = 0; part < MucID::getPartNum(); part++) {
    for(unsigned int gap = 0; gap < MucID::getGapNum(part); gap++) {
      idStripTemp = MucGeoGeneral::Instance()->FindIntersectStrips(part, gap, gPoint, gVect);
      
      for( vector<Identifier>::iterator iter  = idStripTemp.begin(); 
       iter != idStripTemp.end(); ++iter) {
    idStrip.push_back(*iter);
      }
    }
  }
 
  cout << endl << "********************" << endl << endl;
  cout << " Number of intersect strips : " << idStrip.size() << endl << endl;
  
  float stripCenterX, stripCenterY, stripCenterZ;
  HepPoint3D localStripCenter, stripCenter;
  MucGeoStrip* pStrip;
  MucGeoGap*   pGap;
 
 
  for( unsigned int i = 0; i < idStrip.size(); i++) {
    pStrip = MucGeoGeneral::Instance()->GetStrip(MucID::part(idStrip[i]),
                         MucID::seg(idStrip[i]),
                         MucID::gap(idStrip[i]),
                         MucID::strip(idStrip[i]));
    pGap = pStrip->GetGap();
 
    pStrip->GetCenterPos(stripCenterX, 
             stripCenterY, 
             stripCenterZ);
    
    localStripCenter.setX(stripCenterX);
    localStripCenter.setY(stripCenterY);
    localStripCenter.setZ(stripCenterZ);
 
    stripCenter = pGap->TransformToGlobal(localStripCenter);
 
    cout << idStrip[i] << "   center : "  << stripCenter << endl;
  }
 
  vector<HepPoint3D> intersection, intersectionTemp;
  HepPoint3D origin(0,0,0);
 
  for(unsigned int part = 0; part < MucID::getPartNum(); part++) {
    for(unsigned int gap = 0; gap < MucID::getGapNum(part); gap++) {
      intersectionTemp =  MucGeoGeneral::Instance()->FindIntersections(part, gap, gPoint, gVect);
      
      for( vector<HepPoint3D>::iterator iter  = intersectionTemp.begin(); 
       iter != intersectionTemp.end(); ++iter) {
    intersection.push_back(*iter);
      }
    }
  }
 
  cout << endl << "********************" << endl << endl;
  cout << " Number of intersections : " << intersection.size() << endl << endl;
 
 
  for( unsigned int i = 0; i < intersection.size(); i++) {
    cout << intersection[i] << endl;
  }
 
  cout << endl << "********************" << endl << endl;
 
  cout << " MucGeoTest1:: Execute  Successfully!!"
       << endl;
}
 
 
 
 
  
 
 
//HepPoint3D intersection[9];
//  intersection[i] = pmuc->GetIntersection(0, i, gPoint, gVect);