- Timestamp:
- 02/18/05 11:17:24 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SDContainer.cc
r465 r473 189 189 // spectrum_ is [,,,nChan] 190 190 // How annoying. 191 // nPol==2192 191 { 193 192 … … 230 229 // 231 230 // flags is [nChan,nPol] 232 // flags_ is [,,,nChan] 233 // How annoying. 234 // there are no separate flags for XY so make 235 // them up from X and Y 231 // flags_ is [nBeam,nIF,nPol,nChan] 232 // 233 // Note that there are no separate flags for the XY cross polarization so we make 234 // them up from the X and Y which is all the silly code below. This means 235 // that nPol==2 on input but 4 on output 236 236 // 237 237 { … … 271 271 } 272 272 } 273 //274 273 inIt.next(); 275 274 } else if (pos(asap::PolAxis)==2) { 276 outIt.vector() = maskPol01; 275 if (hasXPol) { 276 outIt.vector() = maskPol01; 277 } else { 278 outIt.vector() = inIt.vector(); 279 inIt.next(); 280 } 281 277 282 } else if (pos(asap::PolAxis)==3) { 278 outIt.vector() = maskPol01; 283 if (hasXPol) { 284 outIt.vector() = maskPol01; 285 } else { 286 outIt.vector() = inIt.vector(); 287 inIt.next(); 288 } 279 289 } 280 //281 290 outIt.next(); 282 291 } … … 290 299 Bool hasXpol) 291 300 // 301 // TSys shape is [nPol] 292 302 // Tsys does not depend upon channel but is replicated 293 303 // for simplicity of use. … … 317 327 outIt.vector() = tsys(i++); 318 328 } else { 319 outIt.vector() = sqrt(tsys[0]*tsys[1]); 329 if (hasXpol) { 330 outIt.vector() = sqrt(tsys[0]*tsys[1]); 331 } else { 332 outIt.vector() = tsys(i++); 333 } 320 334 } 321 335 // … … 385 399 // 386 400 // Input [nBeam,nIF,nPol,nChan] 387 // Output [nPol] (drop channel dependency and select first value )401 // Output [nPol] (drop channel dependency and select first value only) 388 402 // 389 403 {
Note:
See TracChangeset
for help on using the changeset viewer.