for(i=0; i<XEN1->GetNdaughters(); i++) {
TGeoNodeMatrix *sm = (TGeoNodeMatrix*)XEN1->GetDaughter(i);
fMatrixOfSM[i] = sm->GetMatrix();
for(i=0; i<XEN1->GetNdaughters(); i++) {
TGeoNodeMatrix *sm = (TGeoNodeMatrix*)XEN1->GetDaughter(i);
fMatrixOfSM[i] = sm->GetMatrix();