PamVMC_update/examples/geom.C

{
  gSystem->Load("libGeom.so");
  new TGeoManager();
  gGeoManager->Import("geometry.root");
  gGeoManager->CloseGeometry();
//Print out global positions of the TSPA sensors
 for (Int_t plNo = 1; plNo<=6; plNo++){
   for(Int_t padNo = 1; padNo<=6; padNo++){
     Int_t i = (plNo-1)*6+padNo;
     TString path = "/PAME_1/SPEB";
     if(plNo<=5){ 
       path+="_1/MGF";
       path += plNo;
     }
     path +="_1/TRP";
     path += plNo;
     path +="_1/TPA";
     path += plNo;
     path +="_";
     path += plNo;
     path +="/TRSL_";
     path += i;
     //cout<<"changing path to "<<path<<endl;
     gGeoManager->cd(path);
     Double_t master[3];
     Double_t local[3];

     for (Int_t j=0; j<3; ++j) local[j] = (Double_t)0.;

     Double_t * m = &master[0];
     Double_t * l = &local[0];
     
     gGeoManager->LocalToMaster(l,m);
     
     cout<<"TSPA_"<<i<<" x,y,z,"<<master[0]<<","<<master[1]<<","<<master[2]<<","<<endl;
     gGeoManager->GetCurrentNode()->GetMatrix()->Print();
     
   }
 }
}
1	{
2	gSystem->Load("libGeom.so");
3	new TGeoManager();
4	gGeoManager->Import("geometry.root");
5	gGeoManager->CloseGeometry();
6	//Print out global positions of the TSPA sensors
7	for (Int_t plNo = 1; plNo<=6; plNo++){
8	for(Int_t padNo = 1; padNo<=6; padNo++){
9	Int_t i = (plNo-1)*6+padNo;
10	TString path = "/PAME_1/SPEB";
11	if(plNo<=5){
12	path+="_1/MGF";
13	path += plNo;
14	}
15	path +="_1/TRP";
16	path += plNo;
17	path +="_1/TPA";
18	path += plNo;
19	path +="_";
20	path += plNo;
21	path +="/TRSL_";
22	path += i;
23	//cout<<"changing path to "<<path<<endl;
24	gGeoManager->cd(path);
25	Double_t master[3];
26	Double_t local[3];
27
28	for (Int_t j=0; j<3; ++j) local[j] = (Double_t)0.;
29
30	Double_t * m = &master[0];
31	Double_t * l = &local[0];
32
33	gGeoManager->LocalToMaster(l,m);
34
35	cout<<"TSPA_"<<i<<" x,y,z,"<<master[0]<<","<<master[1]<<","<<master[2]<<","<<endl;
36	gGeoManager->GetCurrentNode()->GetMatrix()->Print();
37
38	}
39	}
40	}