9{
11 BesEMCRecGeo& geo=BesEMCRecGeo::GetInstance();
12 BesEMCRecCrystal acry;
13
15
16 id=EMCalID::crystal_id(2,22,30);
17 acry=geo.GetCrystal(id);
18 aPoint=geo.GetCCenter(id);
19 acry.Dump();
20 cout<<aPoint<<endl;
21 acry.Checkout();
22 cout<<endl;
23
24 id=EMCalID::crystal_id(2,22,31);
25 acry=geo.GetCrystal(id);
26 aPoint=geo.GetCCenter(id);
27 acry.Dump();
28 cout<<aPoint<<endl;
29 acry.Checkout();
30 cout<<endl;
31
32 id=EMCalID::crystal_id(2,22,1);
33 acry=geo.GetCrystal(id);
34 aPoint=geo.GetCCenter(id);
35 acry.Dump();
36 cout<<aPoint<<endl;
37 acry.Checkout();
38 cout<<endl;
39
40 id=EMCalID::crystal_id(2,22,60);
41 acry=geo.GetCrystal(id);
42 aPoint=geo.GetCCenter(id);
43 acry.Dump();
44 cout<<aPoint<<endl;
45 acry.Checkout();
46 cout<<endl;
47
48}
HepGeom::Point3D< double > HepPoint3D