| 254 |  |  | 
| 255 |  | free(working_line); | 
| 256 |  | return(ntokens); | 
| 257 | + | } | 
| 258 | + |  | 
| 259 | + |  | 
| 260 | + | double SHAPE::getSigmaAt(double costheta, double phi) { | 
| 261 | + |  | 
| 262 | + | vector<SHFunc*>::iterator sigmaIter; | 
| 263 | + | double sigma; | 
| 264 | + |  | 
| 265 | + | sigma = 0.0; | 
| 266 | + |  | 
| 267 | + | for(sigmaIter = sigmaFuncs.begin();  sigmaIter != sigmaFuncs.end(); | 
| 268 | + | ++sigmaIter) | 
| 269 | + | sigma += (*sigmaIter)->getValueAt(costheta, phi); | 
| 270 | + |  | 
| 271 | + | return sigma; | 
| 272 |  | } | 
| 273 | + |  | 
| 274 | + | double SHAPE::getSAt(double costheta, double phi) { | 
| 275 | + |  | 
| 276 | + | vector<SHFunc*>::iterator sIter; | 
| 277 | + | double s; | 
| 278 | + |  | 
| 279 | + | s = 0.0; | 
| 280 | + |  | 
| 281 | + | for(sIter = sFuncs.begin(); sIter != sFuncs.end();  ++sIter) | 
| 282 | + | s += (*sIter)->getValueAt(costheta, phi); | 
| 283 | + |  | 
| 284 | + | return s; | 
| 285 | + | } | 
| 286 | + |  | 
| 287 | + | double SHAPE::getEpsAt(double costheta, double phi) { | 
| 288 | + |  | 
| 289 | + | vector<SHFunc*>::iterator epsIter; | 
| 290 | + | double eps; | 
| 291 | + |  | 
| 292 | + | eps = 0.0; | 
| 293 | + |  | 
| 294 | + | for(epsIter = epsFuncs.begin(); epsIter != epsFuncs.end(); ++epsIter) | 
| 295 | + | eps += (*epsIter)->getValueAt(costheta, phi); | 
| 296 | + |  | 
| 297 | + | return eps; | 
| 298 | + | } |