/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Emc/EmcRecGeoSvc/EmcRecGeoSvc-01-01-07/test/main.cxx

Go to the documentation of this file.

#include "Identifier/Identifier.h"
#include "Identifier/EmcID.h"
#include "EmcRecGeoSvc/EmcRecGeoSvc.h"
#include "EmcRecGeoSvc/EmcRecROOTGeo.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();
}

---