Changeset 7829 for trunk/MagicSoft/Mars/mcalib/MCalibrationQEPix.cc
- Timestamp:
- 08/01/06 14:13:34 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mcalib/MCalibrationQEPix.cc
r5412 r7829 87 87 // by one of the three methods. 88 88 // 89 // 11) A call to GetQECascadesBlindPixel( zenith), GetQECascadesFFactor(zenith) or90 // GetQECascadesPINDiode( zenith) returns then the normalization multiplied with an average QE89 // 11) A call to GetQECascadesBlindPixel(), GetQECascadesFFactor() or 90 // GetQECascadesPINDiode() returns then the normalization multiplied with an average QE 91 91 // folded into a cascades spectrum. This number should be dependent on zenith angle, but 92 92 // this feature is not yet implemented, instead a fixed number gkDefaultAverageQE is used. … … 210 210 const Float_t MCalibrationQEPix::gkLightGuidesEffCT1 = 0.94 ; 211 211 const Float_t MCalibrationQEPix::gkLightGuidesEffCT1Err = 0.03 ; 212 212 213 // -------------------------------------------------------------------------- 213 214 // … … 447 448 // ----------------------------------------------------------------- 448 449 // 449 // Return the average Default QE (depending on zenith angle) 450 // FIXME: The zenith angle dependency is not yet implemented 451 // 452 const Float_t MCalibrationQEPix::GetAverageQE( const Float_t zenith ) const 453 { 454 // return gkDefaultAverageQE ; 450 // Return the average Default QE 451 // 452 const Float_t MCalibrationQEPix::GetAverageQE() const 453 { 455 454 return fAverageQE; 456 455 } … … 458 457 // ----------------------------------------------------------------- 459 458 // 460 // Return the relative variance of the average Default QE (depending on zenith angle) 461 // FIXME: The zenith angle dependency is not yet implemented 462 // 463 const Float_t MCalibrationQEPix::GetAverageQERelVar( const Float_t zenith ) const 459 // Return the relative variance of the average Default QE 460 // 461 const Float_t MCalibrationQEPix::GetAverageQERelVar() const 464 462 { 465 463 return gkDefaultAverageQEErr * gkDefaultAverageQEErr / (gkDefaultAverageQE * gkDefaultAverageQE ); … … 469 467 // 470 468 // Return the relative variance of the average normalization (Blind Pixel Method) 471 // FIXME: The zenith angle dependency is not yet implemented472 469 // 473 470 const Float_t MCalibrationQEPix::GetAvNormBlindPixelRelVar( ) const … … 479 476 // 480 477 // Return the relative variance of the average normalization (Combined Method) 481 // FIXME: The zenith angle dependency is not yet implemented482 478 // 483 479 const Float_t MCalibrationQEPix::GetAvNormCombinedRelVar( ) const … … 498 494 // 499 495 // Return the relative variance of the average normalization (PIN Diode Method) 500 // FIXME: The zenith angle dependency is not yet implemented501 496 // 502 497 const Float_t MCalibrationQEPix::GetAvNormPINDiodeRelVar( ) const … … 624 619 // ------------------------------------------------------------------------------ 625 620 // 626 // Get the light guides efficiency for Cherenkov spectra, 627 // depending on the zenith angle of the telescope 621 // Get the light guides efficiency for Cherenkov spectra. 628 622 // FIXME: Lacking detailed measurement, these number are not yet available 629 623 // for the individual colours and therefore, only one same number is 630 624 // returned, namely gkLightGuidesEffBlue 631 625 // 632 Float_t MCalibrationQEPix::GetLightGuidesEff( const Float_t zenith) const626 Float_t MCalibrationQEPix::GetLightGuidesEff() const 633 627 { 634 628 return gkLightGuidesEffBlue; … … 638 632 // ------------------------------------------------------------------------------ 639 633 // 640 // Get the relative variance of the light guides efficiency for Cherenkov spectra 641 // depending on the zenith angle of the telescope 634 // Get the relative variance of the light guides efficiency for Cherenkov spectra. 642 635 // FIXME: Lacking detailed measurement, these number are not yet available 643 636 // for the individual colours and therefore, only one same number is 644 637 // returned, namely gkLightGuidesEffBlueErr^2 / gkLightGuidesBlueEff^2 645 638 // 646 Float_t MCalibrationQEPix::GetLightGuidesEffRelVar( const Float_t zenith) const639 Float_t MCalibrationQEPix::GetLightGuidesEffRelVar() const 647 640 { 648 641 return gkLightGuidesEffBlueErr * gkLightGuidesEffBlueErr / gkLightGuidesEffBlue / gkLightGuidesEffBlue; … … 846 839 // with the blind pixel method and averaged over the results from the different colours. 847 840 // 848 Float_t MCalibrationQEPix::GetQECascadesBlindPixel( const Float_t zenith) const849 { 850 return fAvNormBlindPixel * GetAverageQE ( zenith);841 Float_t MCalibrationQEPix::GetQECascadesBlindPixel() const 842 { 843 return fAvNormBlindPixel * GetAverageQE(); 851 844 } 852 845 … … 857 850 // different colours. 858 851 // 859 Float_t MCalibrationQEPix::GetQECascadesBlindPixelVar( const Float_t zenith) const860 { 861 return ( GetAvNormBlindPixelRelVar() + GetAverageQERelVar( zenith))862 * GetQECascadesBlindPixel( zenith) * GetQECascadesBlindPixel(zenith);852 Float_t MCalibrationQEPix::GetQECascadesBlindPixelVar() const 853 { 854 return ( GetAvNormBlindPixelRelVar() + GetAverageQERelVar()) 855 * GetQECascadesBlindPixel() * GetQECascadesBlindPixel(); 863 856 } 864 857 … … 869 862 // different colours. 870 863 // 871 Float_t MCalibrationQEPix::GetQECascadesBlindPixelErr( const Float_t zenith) const872 { 873 const Float_t var = GetQECascadesBlindPixelVar( zenith);864 Float_t MCalibrationQEPix::GetQECascadesBlindPixelErr() const 865 { 866 const Float_t var = GetQECascadesBlindPixelVar(); 874 867 875 868 if (var < 0.) … … 885 878 // from the different colours. 886 879 // 887 Float_t MCalibrationQEPix::GetQECascadesCombined( const Float_t zenith )const888 { 889 return fAvNormCombined * GetAverageQE ( zenith);880 Float_t MCalibrationQEPix::GetQECascadesCombined() const 881 { 882 return fAvNormCombined * GetAverageQE(); 890 883 } 891 884 … … 896 889 // different colours. 897 890 // 898 Float_t MCalibrationQEPix::GetQECascadesCombinedErr( const Float_t zenith )const899 { 900 const Float_t var = GetQECascadesCombinedVar( zenith);891 Float_t MCalibrationQEPix::GetQECascadesCombinedErr() const 892 { 893 const Float_t var = GetQECascadesCombinedVar(); 901 894 902 895 if (var < 0.) … … 912 905 // different colours. 913 906 // 914 Float_t MCalibrationQEPix::GetQECascadesCombinedVar( const Float_t zenith )const915 { 916 return ( GetAvNormCombinedRelVar() + GetAverageQERelVar( zenith))917 * GetQECascadesCombined( zenith) * GetQECascadesCombined(zenith);907 Float_t MCalibrationQEPix::GetQECascadesCombinedVar() const 908 { 909 return ( GetAvNormCombinedRelVar() + GetAverageQERelVar()) 910 * GetQECascadesCombined() * GetQECascadesCombined(); 918 911 } 919 912 … … 923 916 // with the F-Factor method and averaged over the results from the different colours. 924 917 // 925 Float_t MCalibrationQEPix::GetQECascadesFFactor( const Float_t zenith) const926 { 927 return fAvNormFFactor * GetAverageQE ( zenith);918 Float_t MCalibrationQEPix::GetQECascadesFFactor() const 919 { 920 return fAvNormFFactor * GetAverageQE(); 928 921 } 929 922 … … 934 927 // different colours. 935 928 // 936 Float_t MCalibrationQEPix::GetQECascadesFFactorErr( const Float_t zenith) const937 { 938 const Float_t var = GetQECascadesFFactorVar( zenith);929 Float_t MCalibrationQEPix::GetQECascadesFFactorErr() const 930 { 931 const Float_t var = GetQECascadesFFactorVar(); 939 932 940 933 if (var < 0.) … … 950 943 // different colours. 951 944 // 952 Float_t MCalibrationQEPix::GetQECascadesFFactorVar( const Float_t zenith) const953 { 954 return ( GetAvNormFFactorRelVar() + GetAverageQERelVar( zenith))955 * GetQECascadesFFactor( zenith) * GetQECascadesFFactor(zenith);945 Float_t MCalibrationQEPix::GetQECascadesFFactorVar() const 946 { 947 return ( GetAvNormFFactorRelVar() + GetAverageQERelVar()) 948 * GetQECascadesFFactor() * GetQECascadesFFactor(); 956 949 } 957 950 … … 961 954 // with the PIN Diode method and averaged over the results from the different colours. 962 955 // 963 Float_t MCalibrationQEPix::GetQECascadesPINDiode( const Float_t zenith) const964 { 965 return fAvNormPINDiode * GetAverageQE ( zenith);956 Float_t MCalibrationQEPix::GetQECascadesPINDiode() const 957 { 958 return fAvNormPINDiode * GetAverageQE(); 966 959 } 967 960 … … 972 965 // different colours. 973 966 // 974 Float_t MCalibrationQEPix::GetQECascadesPINDiodeErr( const Float_t zenith )const975 { 976 const Float_t var = GetQECascadesPINDiodeVar( zenith);967 Float_t MCalibrationQEPix::GetQECascadesPINDiodeErr() const 968 { 969 const Float_t var = GetQECascadesPINDiodeVar(); 977 970 978 971 if (var < 0.) … … 988 981 // different colours. 989 982 // 990 Float_t MCalibrationQEPix::GetQECascadesPINDiodeVar( const Float_t zenith )const991 { 992 return ( GetAvNormPINDiodeRelVar() + GetAverageQERelVar( zenith))993 * GetQECascadesPINDiode( zenith) * GetQECascadesPINDiode(zenith);983 Float_t MCalibrationQEPix::GetQECascadesPINDiodeVar() const 984 { 985 return ( GetAvNormPINDiodeRelVar() + GetAverageQERelVar()) 986 * GetQECascadesPINDiode() * GetQECascadesPINDiode(); 994 987 } 995 988
Note:
See TracChangeset
for help on using the changeset viewer.