Changeset 8917 for trunk/MagicSoft/Mars/mcalib
- Timestamp:
- 06/03/08 12:29:27 (17 years ago)
- Location:
- trunk/MagicSoft/Mars/mcalib
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.cc
r8147 r8917 18 18 ! Author(s): Markus Gaug 02/2004 <mailto:markus@ifae.es> 19 19 ! 20 ! Copyright: MAGIC Software Development, 2000-200 420 ! Copyright: MAGIC Software Development, 2000-2008 21 21 ! 22 22 \* ======================================================================== */ … … 65 65 // FADC slices extraction results, but can be retrieved converted to High-Gain amplifications 66 66 // by calls to: GetConvertedMean() or GetConvertedSigma() 67 // 68 // Inline Functions: 69 // ----------------- 70 // 71 // GetPedRms(): Get the pedestals RMS: Test bit kHiGainSaturation: 72 // If yes, return square root of fLoGainPedRmsSquare (if greater than 0, 73 // otherwise -1.), If no, return fPedRms 74 // 75 // GetConvertedMean(): Get the Low Gain Mean Charge converted to High Gain 76 // amplification: Returns fLoGainMean multiplied with fConversionHiLo if 77 // IsHiGainSaturation(), else return fHiGainMean 78 // 79 // GetConvertedSigma(): Get the Low Gain Sigma converted to High Gain 80 // amplification: Returns fLoGainSigma multiplied with fConversionHiLo 81 // if IsHiGainSaturation() else return fHiGainSigma 82 // 83 // GetConvertedRSigmaSquare(): Get the reduced Sigma Square: 84 // If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), return -1. 85 // Test bit kHiGainSaturation: 86 // If yes, return fRSigmaSquare, multiplied with fConversionHiLo^2, 87 // If no , return fRSigmaSquare 88 // 89 // GetPheFFactorMethodRelVar(): Get the relative variance on the number of 90 // photo-electrons (F-Factor Method): If fPheFFactorMethodVar is smaller 91 // than 0 (i.e. has not yet been set), return -1. If fPheFFactorMethod 92 // is 0, return -1. Else returns fPheFFactorMethodVar / fPheFFactorMethod^2 93 // 94 // GetMeanConvFADC2PheErr(): Get the error on the mean conversion factor 95 // (FFactor Method): If fMeanConvFADC2PheVar is smaller than 0 (i.e. has 96 // not yet been set), return -1. Else returns the square root of 97 // fMeanConvFADC2PheVar 98 // 67 99 // 68 100 // See also: MCalibrationChargeCam, MCalibrationChargeCalc, … … 88 120 ///////////////////////////////////////////////////////////////////////////// 89 121 #include "MCalibrationChargePix.h" 122 123 #include <TMath.h> 90 124 91 125 #include "MLog.h" … … 264 298 } 265 299 266 267 // Inline Functions: 268 // ----------------- 269 // 270 // GetConversionHiLoErr: 271 // Get the conversion Error Hi-Gain to Low-Gain: 272 // If fConversionHiLoVar is smaller than 0 (i.e. has not yet 273 // been set), return -1. 300 // -------------------------------------------------------------------------- 301 // 302 // Get the conversion Error Hi-Gain to Low-Gain: 303 // If fConversionHiLoVar is smaller than 0 (i.e. has not yet 304 // been set), return -1. 274 305 // 275 // GetPedRms(): Get the pedestals RMS: Test bit kHiGainSaturation: 276 // If yes, return square root of fLoGainPedRmsSquare (if greater than 0, 277 // otherwise -1.), If no, return fPedRms 278 // 279 // GetConvertedMean(): Get the Low Gain Mean Charge converted to High Gain 280 // amplification: Returns fLoGainMean multiplied with fConversionHiLo if 281 // IsHiGainSaturation(), else return fHiGainMean 282 // 283 // GetConvertedMeanErr(): Get the Error of the converted Low Gain Mean: 284 // Returns -1 if the variable fLoGainMean or fLoGainMeanVar are smaller than 0. 285 // Returns the square root of the quadratic sum of the relative variances of 286 // the fLoGainMean and fConversionHiLo, mulitplied with GetConvertedMean() 287 // in case of HiGain Saturation, 288 // else return GetMeanErr() 289 // 290 // GetConvertedSigma(): Get the Low Gain Sigma converted to High Gain 291 // amplification: Returns fLoGainSigma multiplied with fConversionHiLo 292 // if IsHiGainSaturation() else return fHiGainSigma 293 // 294 // GetConvertedSigmaErr(): Get the Error of the converted Sigma: 295 // Returns -1 if the variable fLoGainSigma or fLoGainSigmaVar are smaller than 0. 296 // if IsHiGainSaturatio() 297 // returns the square root of the quadratic sum of the relative variances of 298 // the fLoGainSigma and fConversionHiLo, mulitplied with GetConvertedSigma() 299 // else returns GetSigmaErr() 300 // 301 // GetConvertedRSigma(): Get the converted reduced Sigma: 302 // If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), return -1. 303 // Test bit kHiGainSaturation: 304 // If yes, return square root of fRSigmaSquare, multiplied with fConversionHiLo, 305 // If no , return square root of fRSigmaSquare 306 // 307 // GetRSigma(): Get the reduced Sigma: 308 // If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), return -1. 309 // 310 // GetConvertedRSigmaSquare(): Get the reduced Sigma Square: 311 // If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), return -1. 312 // Test bit kHiGainSaturation: 313 // If yes, return fRSigmaSquare, multiplied with fConversionHiLo^2, 314 // If no , return fRSigmaSquare 306 Float_t MCalibrationChargePix::GetConversionHiLoErr() const 307 { 308 return fConversionHiLoVar<0 ? -1 : TMath::Sqrt(fConversionHiLoVar); 309 } 310 311 // -------------------------------------------------------------------------- 312 // 313 // Get the Error of the converted Low Gain Mean: Returns -1 if the 314 // variable fLoGainMean or fLoGainMeanVar are smaller than 0. Returns the 315 // square root of the quadratic sum of the relative variances of the 316 // fLoGainMean and fConversionHiLo, mulitplied with GetConvertedMean() in 317 // case of HiGain Saturation, else return GetMeanErr() 318 // 319 Float_t MCalibrationChargePix::GetConvertedMeanErr() const 320 { 321 if (!IsHiGainSaturation()) 322 return GetMeanErr(); 323 324 const Float_t logainrelvar = GetLoGainMeanRelVar(); 325 return logainrelvar<0 ? -1 : TMath::Sqrt(logainrelvar + GetConversionHiLoRelVar()) * GetConvertedMean(); 326 } 327 328 // -------------------------------------------------------------------------- 329 // 330 // Get the Error of the converted Sigma: Returns -1 if the variable 331 // fLoGainSigma or fLoGainSigmaVar are smaller than 0. if 332 // IsHiGainSaturatio() returns the square root of the quadratic sum of the 333 // relative variances of the fLoGainSigma and fConversionHiLo, mulitplied 334 // with GetConvertedSigma() else returns GetSigmaErr() 335 // 336 Float_t MCalibrationChargePix::GetConvertedSigmaErr() const 337 { 338 if (!IsHiGainSaturation()) 339 return GetSigmaErr(); 340 341 if (fLoGainSigmaVar<0 || fLoGainSigma<0) 342 return -1.; 343 344 const Float_t sigmaRelVar = fLoGainSigmaVar/(fLoGainSigma*fLoGainSigma); 345 return TMath::Sqrt(sigmaRelVar+GetConversionHiLoRelVar()) * GetConvertedSigma(); 346 } 347 348 // -------------------------------------------------------------------------- 349 // 350 // Get the converted reduced Sigma: If fRSigmaSquare is smaller than 0 351 // (i.e. has not yet been set), return -1. Test bit kHiGainSaturation: If 352 // yes, return square root of fRSigmaSquare, multiplied with 353 // fConversionHiLo, If no , return square root of fRSigmaSquare 354 // 355 Float_t MCalibrationChargePix::GetConvertedRSigma() const 356 { 357 if (fRSigmaSquare < 0) 358 return -1; 359 360 const Float_t rsigma = TMath::Sqrt(fRSigmaSquare); 361 return IsHiGainSaturation() ? rsigma*fConversionHiLo : rsigma ; 362 } 363 364 // -------------------------------------------------------------------------- 365 // 366 Float_t MCalibrationChargePix::GetMeanConvFADC2PheErr() const 367 { 368 return fMeanConvFADC2PheStatVar<0 ? -1 : TMath::Sqrt(fMeanConvFADC2PheStatVar); 369 } 370 371 // -------------------------------------------------------------------------- 372 // 373 Float_t MCalibrationChargePix::GetMeanConvFADC2PheSystErr() const 374 { 375 return fMeanConvFADC2PheSystVar<0 ? -1 : TMath::Sqrt(fMeanConvFADC2PheSystVar); 376 } 377 378 // -------------------------------------------------------------------------- 379 // 380 Float_t MCalibrationChargePix::GetMeanConvFADC2PheTotErr() const 381 { 382 if (fMeanConvFADC2PheSystVar<0 || fMeanConvFADC2PheStatVar<0) 383 return -1.; 384 385 return TMath::Sqrt(fMeanConvFADC2PheSystVar+fMeanConvFADC2PheStatVar); 386 } 387 388 // -------------------------------------------------------------------------- 389 // 390 // Get the error on the mean total F-Factor of the signal readout 391 // (F-Factor Method): If fMeanFFactorFADC2PhotVar is smaller than 0 392 // (i.e. has not yet been set), return -1. Else returns the square 393 // root of fMeanFFactorFADC2PhotVar 394 // 395 Float_t MCalibrationChargePix::GetMeanFFactorFADC2PhotErr() const 396 { 397 return fMeanFFactorFADC2PhotVar<0 ? -1. : TMath::Sqrt(fMeanFFactorFADC2PhotVar); 398 } 399 400 // -------------------------------------------------------------------------- 401 // 402 // Get the relative variance of the conversion factor between higain and logain: 403 // - If fConversionHiLo is 0, return -1. 404 // - If fConversionHiLoVar is smaller than 0, return -1. 405 // - Else returns: fConversionHiLoVar / fConversionHiLo^2 406 // 407 const Float_t MCalibrationChargePix::GetConversionHiLoRelVar() const 408 { 409 410 if (fConversionHiLoVar < 0.) 411 return -1.; 412 413 if (fConversionHiLo == 0.) 414 return -1.; 415 416 return fConversionHiLoVar / (fConversionHiLo * fConversionHiLo); 417 } 418 419 // -------------------------------------------------------------------------- 420 // 421 // Get the relative variance of the electronics pedestal RMS 422 // - returns rel. variance of gkElectronicPedRms 423 // 424 const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar() const 425 { 426 return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRms / gkElectronicPedRms; 427 } 428 429 430 // -------------------------------------------------------------------------- 431 // 432 // Get the relative variance of the conversion factor between higain and logain: 433 // - If gkFFactor is 0, return -1. 434 // - If gkFFactorErr is smaller than 0, return -1. 435 // - Else returns: gkFFactorErr^2 / gkFFactor*^2 436 // 437 const Float_t MCalibrationChargePix::GetFFactorRelVar() const 438 { 439 440 if (gkFFactorErr < 0.) 441 return -1.; 442 443 if (gkFFactor == 0.) 444 return -1.; 445 446 return gkFFactorErr * gkFFactorErr / (gkFFactor * gkFFactor); 447 } 448 449 // -------------------------------------------------------------------------- 450 // 451 Float_t MCalibrationChargePix::GetPedErr() const 452 { 453 return TMath::Sqrt(fPedVar); 454 } 455 456 // -------------------------------------------------------------------------- 457 // 458 Float_t MCalibrationChargePix::GetPedRms() const 459 { 460 if (!IsHiGainSaturation()) 461 return fPedRms; 462 463 return fLoGainPedRmsSquare<0 ? -1 : TMath::Sqrt(fLoGainPedRmsSquare); 464 } 465 466 // -------------------------------------------------------------------------- 467 // 468 // Get the Error of the pedestals RMS: 469 // - Test bit kHiGainSaturation: 470 // If yes, return square root of (0.25*fLoGainPedRmsSquareVar/ fLoGainPedRmsSquare) (if greater than 0, otherwise -1.) 471 // If no , return square root of (fPedVar) (if greater than 0, otherwise -1.), divided by 2. 472 // 473 Float_t MCalibrationChargePix::GetPedRmsErr() const 474 { 475 if (IsHiGainSaturation()) 476 return fLoGainPedRmsSquareVar < 0. ? -1. :TMath::Sqrt(0.25*fLoGainPedRmsSquareVar/fLoGainPedRmsSquare); 477 else 478 return fPedVar < 0. ? -1. : TMath::Sqrt(fPedVar)/2.; 479 } 480 481 // -------------------------------------------------------------------------- 315 482 // 316 483 // GetPheFFactorMethodErr(): Get the error on the number of photo-electrons … … 319 486 // fPheFFactorMethodVar 320 487 // 321 // GetMeanFFactorFADC2PhotErr(): Get the error on the mean total F-Factor 322 // of the signal readout (F-Factor Method): If fMeanFFactorFADC2PhotVar 323 // is smaller than 0 (i.e. has not yet been set), return -1. Else returns 324 // the square root of fMeanFFactorFADC2PhotVar 325 // 326 // GetPheFFactorMethodRelVar(): Get the relative variance on the number of 327 // photo-electrons (F-Factor Method): If fPheFFactorMethodVar is smaller 328 // than 0 (i.e. has not yet been set), return -1. If fPheFFactorMethod 329 // is 0, return -1. Else returns fPheFFactorMethodVar / fPheFFactorMethod^2 330 // 331 // GetMeanConvFADC2PheErr(): Get the error on the mean conversion factor 332 // (FFactor Method): If fMeanConvFADC2PheVar is smaller than 0 (i.e. has 333 // not yet been set), return -1. Else returns the square root of 334 // fMeanConvFADC2PheVar 335 336 // -------------------------------------------------------------------------- 337 // 338 // Get the relative variance of the conversion factor between higain and logain: 339 // - If fConversionHiLo is 0, return -1. 340 // - If fConversionHiLoVar is smaller than 0, return -1. 341 // - Else returns: fConversionHiLoVar / fConversionHiLo^2 342 // 343 const Float_t MCalibrationChargePix::GetConversionHiLoRelVar() const 344 { 345 346 if (fConversionHiLoVar < 0.) 347 return -1.; 348 349 if (fConversionHiLo == 0.) 350 return -1.; 351 352 return fConversionHiLoVar / (fConversionHiLo * fConversionHiLo); 353 } 354 355 // -------------------------------------------------------------------------- 356 // 357 // Get the relative variance of the electronics pedestal RMS 358 // - returns rel. variance of gkElectronicPedRms 359 // 360 const Float_t MCalibrationChargePix::GetElectronicPedRmsRelVar() const 361 { 362 return gkElectronicPedRmsErr * gkElectronicPedRmsErr / gkElectronicPedRms / gkElectronicPedRms; 363 } 364 365 366 // -------------------------------------------------------------------------- 367 // 368 // Get the relative variance of the conversion factor between higain and logain: 369 // - If gkFFactor is 0, return -1. 370 // - If gkFFactorErr is smaller than 0, return -1. 371 // - Else returns: gkFFactorErr^2 / gkFFactor*^2 372 // 373 const Float_t MCalibrationChargePix::GetFFactorRelVar() const 374 { 375 376 if (gkFFactorErr < 0.) 377 return -1.; 378 379 if (gkFFactor == 0.) 380 return -1.; 381 382 return gkFFactorErr * gkFFactorErr / (gkFFactor * gkFFactor); 383 } 384 385 386 387 // -------------------------------------------------------------------------- 388 // 389 // Get the Error of the pedestals RMS: 390 // - Test bit kHiGainSaturation: 391 // If yes, return square root of (0.25*fLoGainPedRmsSquareVar/ fLoGainPedRmsSquare) (if greater than 0, otherwise -1.) 392 // If no , return square root of (fPedVar) (if greater than 0, otherwise -1.), divided by 2. 393 // 394 Float_t MCalibrationChargePix::GetPedRmsErr() const 395 { 396 if (IsHiGainSaturation()) 397 if (fLoGainPedRmsSquareVar < 0.) 398 return -1.; 399 else 400 return TMath::Sqrt(0.25*fLoGainPedRmsSquareVar/fLoGainPedRmsSquare); 401 else 402 if (fPedVar < 0.) 403 return -1.; 404 else 405 return TMath::Sqrt(fPedVar)/2.; 406 } 407 488 Float_t MCalibrationChargePix::GetPheFFactorMethodErr() const 489 { 490 return fPheFFactorMethodStatVar<0 ? -1 : TMath::Sqrt(fPheFFactorMethodStatVar); 491 } 492 493 // -------------------------------------------------------------------------- 494 // 495 Float_t MCalibrationChargePix::GetPheFFactorMethodSystErr() const 496 { 497 return fPheFFactorMethodSystVar<0 ? -1 : TMath::Sqrt(fPheFFactorMethodSystVar); 498 } 499 500 // -------------------------------------------------------------------------- 501 // 502 Float_t MCalibrationChargePix::GetPheFFactorMethodTotErr() const 503 { 504 if (fPheFFactorMethodStatVar<0 || fPheFFactorMethodSystVar<0) 505 return -1.; 506 507 return TMath::Sqrt(fPheFFactorMethodStatVar+fPheFFactorMethodSystVar); 508 } 408 509 409 510 // -------------------------------------------------------------------------- … … 439 540 // -------------------------------------------------------------------------- 440 541 // 542 // Get the reduced Sigma: 543 // If fRSigmaSquare is smaller than 0 (i.e. has not yet been set), 544 // return -1. 545 // 546 Float_t MCalibrationChargePix::GetRSigma() const 547 { 548 return fRSigmaSquare<0 ? -1 : TMath::Sqrt(fRSigmaSquare); 549 } 550 551 // -------------------------------------------------------------------------- 552 // 441 553 // Get the error of the reduced Sigma: 442 554 // - If fRSigmaSquareVar is smaller than 0 (i.e. has not yet been set), return -1. -
trunk/MagicSoft/Mars/mcalib/MCalibrationChargePix.h
r7189 r8917 96 96 Float_t GetAbsTimeMean () const { return fAbsTimeMean; } 97 97 Float_t GetAbsTimeRms () const { return fAbsTimeRms; } 98 98 99 Float_t GetConversionHiLo () const { return fConversionHiLo; } 99 Float_t GetConversionHiLoErr () const { return fConversionHiLoVar<0 ? -1 : TMath::Sqrt(fConversionHiLoVar); }100 Float_t GetConversionHiLoErr () const; 100 101 Float_t GetConversionHiLoSigma() const { return fConversionHiLoSigma; } 102 101 103 Float_t GetConvertedMean () const { return IsHiGainSaturation() ? fLoGainMean * fConversionHiLo : fHiGainMean; } 102 Float_t GetConvertedMeanErr () const 103 { 104 if (!IsHiGainSaturation()) 105 return GetMeanErr(); 106 const Float_t logainrelvar = GetLoGainMeanRelVar(); 107 return logainrelvar<0 ? -1 : TMath::Sqrt(logainrelvar + GetConversionHiLoRelVar()) * GetConvertedMean(); 108 } 104 Float_t GetConvertedMeanErr () const; 105 109 106 Float_t GetConvertedSigma() const { return IsHiGainSaturation() ? fLoGainSigma * fConversionHiLo : fHiGainSigma; } 110 Float_t GetConvertedSigmaErr() const 111 { 112 if (!IsHiGainSaturation()) 113 return GetSigmaErr(); 107 Float_t GetConvertedSigmaErr() const; 114 108 115 if (fLoGainSigmaVar<0 || fLoGainSigma<0) 116 return -1.; 117 118 const Float_t sigmaRelVar = fLoGainSigmaVar/(fLoGainSigma*fLoGainSigma); 119 return TMath::Sqrt(sigmaRelVar+GetConversionHiLoRelVar()) * GetConvertedSigma(); 120 } 121 Float_t GetConvertedRSigma() const 122 { 123 if (fRSigmaSquare < 0) 124 return -1; 125 const Float_t rsigma = TMath::Sqrt(fRSigmaSquare); 126 return IsHiGainSaturation() ? rsigma*fConversionHiLo : rsigma ; 127 } 109 Float_t GetConvertedRSigma() const; 128 110 Float_t GetConvertedRSigmaErr() const; 129 111 Float_t GetConvertedRSigmaSquare() const … … 133 115 return IsHiGainSaturation() ? fRSigmaSquare*fConversionHiLo*fConversionHiLo : fRSigmaSquare ; 134 116 } 117 135 118 Float_t GetMeanConvFADC2Phe() const { return fMeanConvFADC2Phe; } 136 Float_t GetMeanConvFADC2PheErr () const { return fMeanConvFADC2PheStatVar<0 ? -1 : TMath::Sqrt(fMeanConvFADC2PheStatVar); } 137 Float_t GetMeanConvFADC2PheSystErr() const { return fMeanConvFADC2PheSystVar<0 ? -1 : TMath::Sqrt(fMeanConvFADC2PheSystVar); } 138 Float_t GetMeanConvFADC2PheTotErr() const 139 { 140 if (fMeanConvFADC2PheSystVar<0 || fMeanConvFADC2PheStatVar<0) 141 return -1.; 142 return TMath::Sqrt(fMeanConvFADC2PheSystVar+fMeanConvFADC2PheStatVar); 143 } 119 Float_t GetMeanConvFADC2PheErr() const; 120 Float_t GetMeanConvFADC2PheSystErr() const; 121 Float_t GetMeanConvFADC2PheTotErr() const; 122 144 123 Float_t GetFFactorFADC2Phe () const { return gkFFactor; } 145 124 Float_t GetMeanConvFADC2PheVar () const { return fMeanConvFADC2PheStatVar; } … … 147 126 148 127 Float_t GetMeanFFactorFADC2Phot () const { return fMeanFFactorFADC2Phot; } 149 Float_t GetMeanFFactorFADC2PhotErr() const { return fMeanFFactorFADC2PhotVar<0 ? -1. : TMath::Sqrt(fMeanFFactorFADC2PhotVar); }128 Float_t GetMeanFFactorFADC2PhotErr() const; 150 129 Float_t GetMeanFFactorFADC2PhotVar() const { return fMeanFFactorFADC2PhotVar; } 151 Int_t GetNumSaturated () const { return fNumSaturated; } 152 Float_t GetPed () const { return fPed; } 153 Float_t GetPedErr () const { return TMath::Sqrt(fPedVar); } 154 Float_t GetPedRms () const 155 { 156 if (!IsHiGainSaturation()) 157 return fPedRms; 158 return fLoGainPedRmsSquare<0 ? -1 : TMath::Sqrt(fLoGainPedRmsSquare); 159 } 160 Float_t GetPedRmsErr () const; 161 Float_t GetPheFFactorMethod () const { return fPheFFactorMethod; } 162 Float_t GetPheFFactorMethodErr () const { return fPheFFactorMethodStatVar<0 ? -1 : TMath::Sqrt(fPheFFactorMethodStatVar); } 163 Float_t GetPheFFactorMethodSystErr() const { return fPheFFactorMethodSystVar<0 ? -1 : TMath::Sqrt(fPheFFactorMethodSystVar); } 164 Float_t GetPheFFactorMethodTotErr () const 165 { 166 if (fPheFFactorMethodStatVar<0 || fPheFFactorMethodSystVar<0) 167 return -1.; 168 return TMath::Sqrt(fPheFFactorMethodStatVar+fPheFFactorMethodSystVar); 169 } 130 131 Int_t GetNumSaturated () const { return fNumSaturated; } 132 133 Float_t GetPed() const { return fPed; } 134 Float_t GetPedErr() const; 135 Float_t GetPedRms() const; 136 Float_t GetPedRmsErr() const; 137 138 Float_t GetPheFFactorMethod() const { return fPheFFactorMethod; } 139 Float_t GetPheFFactorMethodErr() const; 140 Float_t GetPheFFactorMethodSystErr() const; 141 Float_t GetPheFFactorMethodTotErr() const; 142 170 143 Float_t GetPheFFactorMethodVar () const { return fPheFFactorMethodStatVar; } 171 144 Float_t GetPheFFactorMethodSystVar() const { return fPheFFactorMethodSystVar; } 172 145 Float_t GetPheFFactorMethodRelVar () const { return fPheFFactorMethodStatVar<=0 ? -1 : fPheFFactorMethodStatVar / (fPheFFactorMethod * fPheFFactorMethod); } 173 146 Float_t GetPheFFactorMethodRelSystVar() const { return fPheFFactorMethodSystVar<=0 ? -1. : fPheFFactorMethodSystVar / (fPheFFactorMethod * fPheFFactorMethod); } 174 Float_t GetRSigma () const { return fRSigmaSquare<0 ? -1 : TMath::Sqrt(fRSigmaSquare); }175 147 176 Float_t GetRSigmaErr () const; 177 Float_t GetRSigmaRelVar () const; 178 Float_t GetRSigmaPerCharge () const; 148 Float_t GetRSigma() const; 149 Float_t GetRSigmaErr() const; 150 Float_t GetRSigmaRelVar() const; 151 Float_t GetRSigmaPerCharge() const; 179 152 Float_t GetRSigmaPerChargeErr() const; 180 153
Note:
See TracChangeset
for help on using the changeset viewer.