For a given FanF the function computes the subfan F1 of all smooth cones.
Let's consider the fan consisting of the following three dimensional cone and all of its faces:
i1 : C = coneFromVData matrix {{1,-1,0},{1,1,0},{1,1,1}}
o1 = C
o1 : Cone
i2 : F = fan C
o2 = F
o2 : Fan
This cone is not smooth, therefore also the fan is not. But if we remove the interior and one of the two dimensional faces the resulting subfan is smooth.