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