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 <iomanip>
#include <iostream>
#include <string>
#include <vector>
// #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);