/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Muc/MucGeomSvc/MucGeomSvc-00-02-25/test/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);

---