Emc/EmcRecGeoSvc/test/main.cxx

Go to the documentation of this file.

#include "EmcRecGeoSvc/EmcRecGeoSvc.h"
#include "EmcRecGeoSvc/EmcRecROOTGeo.h"
#include "Identifier/EmcID.h"
#include "Identifier/Identifier.h"
#include <fstream>
 
using namespace std;
 
int main() {
  Identifier      id;
  EmcRecBarrelGeo geo;
  // EmcRecEndCapGeo geo;
  EmcRecROOTGeo aROOTGeo;
  aROOTGeo.InitFromXML();
 
  int nphi[6];
  nphi[0] = 64;
  nphi[1] = 64;
  nphi[2] = 80;
  nphi[3] = 80;
  nphi[4] = 96;
  nphi[5] = 96;
 
  // int layer=0;
  // ofstream out;
  // out.open("endcap.txt");
  /*for(int part=0;part<2;part++) {
    for(int layer=0;layer<6;layer++) {
      for(int i=0;i<nphi[layer];i++) {
        id = EmcID::crystal_id(part*2,layer,i);
        if(part==0&&(layer==1)&&(i==16||i==17||i==18)) {
          cout<<part*2<<"\t"<<layer<<"\t"<<i<<endl;
          cout<<geo.GetCrystal(id)
            <<"\n-------------\n"<<aROOTGeo.GetCrystal(id)<<"\n"<<std::endl;
        }
          //cout<<id<<"\t"<<geo.GetCFrontCenter(id)
          //<<"\t"<<aROOTGeo.GetCFrontCenter(id)<<endl;
          //out<<part*2<<"\t"<<layer<<"\t"<<i<<"\t"
            //<<(aROOTGeo.GetCCenter(id)-geo.GetCCenter(id)).mag()<<"\t"
            //<<(aROOTGeo.GetCFrontCenter(id)-geo.GetCFrontCenter(id)).mag()<<endl;
      }
    }
  }*/
  // id=EmcID::crystal_id(0,1,);
 
  ofstream out;
  out.open( "barrel.txt" );
  for ( int part = 1; part < 2; part++ )
  {
    for ( int phi = 0; phi < 120; phi++ )
    {
      for ( int theta = 0; theta < 44; theta++ )
      {
        // cout<<part*2<<"\t"<<phi<<"\t"<<theta<<endl;
        id = EmcID::crystal_id( part, theta, phi );
        if ( ( phi == 31 || phi == 30 ) && ( theta == 21 || theta == 22 ) )
        {
          std::cout << id << "\t" << phi << "\t" << theta << endl;
          cout << geo.GetCrystal( id ) << "\n-------------\n"
               << aROOTGeo.GetCrystal( id ) << "\n"
               << std::endl;
        }
        out << part << "\t" << phi << "\t" << theta << "\t"
            << ( aROOTGeo.GetCCenter( id ) - geo.GetCCenter( id ) ).mag() << "\t"
            << ( aROOTGeo.GetCFrontCenter( id ) - geo.GetCFrontCenter( id ) ).mag() << endl;
      }
    }
  }
  out.close();
}