00001 00002 00003 00004 #include "MdcGeomSvc/MdcGeoWire.h" 00005 #include "MdcGeomSvc/MdcGeomSvc.h" 00006 00007 00008 // according to the mathematicial formula: 00009 // Sg = L^2*g*rhol/8*T 00010 // rhol is the linear density 00011 // rhol is fixed at 9.47E-3 g/m 00012 const double MdcGeoWire::Sag(void) const{ 00013 00014 if(MdcGeomSvc::getSagFlag()){ 00015 // unit of wire length is mm 00016 const double L = (fForward - fBackward).mag(); 00017 // sag unit is mm 00018 return (L*L*1.184E-6/fTension); 00019 } 00020 // do not consider wire sagita 00021 else return 0.; 00022 } 00023 00024 00025 // calculate the sag at length z 00026 // attention, the z is not the z-coordinate of the point 00027 // but the length from this point to the backward point 00028 const double MdcGeoWire::Sagz(const double z) const{ 00029 return z*0.; //FIXME 00030 } 00031 00032 00033 00034