00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046 #include <iostream>
00047 #include "CLHEP/Matrix/SymMatrix.h"
00048
00049 using namespace CLHEP;
00050
00051 static const double Large( 1.0e13 );
00052
00053
00054
00055
00056
00057
00058 bool Ext_err_valid( bool msg, HepSymMatrix &error, const int dimension )
00059 {
00060 bool valid( 1 );
00061 double trace( 0 );
00062
00063 for( int i = 1; i<=dimension; i++ ){
00064 double elem( error( i, i ) );
00065 trace += elem;
00066 if( elem < 0.0 ){
00067 valid = 0;
00068 if( msg ){
00069 std::cout << "%ERROR detected at Ext_err_valid: error("
00070 << i << "," << i << ") = " << elem << " < 0.0. "
00071 << "Force to 0.0." << std::endl;
00072 }
00073 error( i, i ) = 0.0;
00074 } else if( elem > Large ){
00075 valid = 0;
00076 if( msg ){
00077 std::cout << "%ERROR detected at Ext_err_valid: error("
00078 << i << "," << i << ") = " << elem << " > " << Large
00079 << ". Force to " << Large << std::endl;
00080 }
00081 error( i, i ) = Large;
00082 }
00083 }
00084 if( !trace ) valid = 0;
00085 return( valid );
00086 }
00087
00088
00089
00090
00091
00092 bool Ext_err_valid( bool msg, const HepSymMatrix &error, const int dimension )
00093 {
00094 bool valid( 1 );
00095 double trace( 0 );
00096
00097 for( int i = 1; i<=dimension; i++ ){
00098 double elem( error( i, i ) );
00099 trace += elem;
00100 if( elem < 0.0 ){
00101 valid = 0;
00102 if( msg ){
00103 std::cout << "%ERROR detected at Ext_err_valid: error matrix: error("
00104 << i << "," << i << ")= " << elem << " < 0.0." << std::endl;
00105 }
00106 } else if( elem > Large ){
00107 valid = 0;
00108 if( msg ){
00109 std::cout << "%ERROR detected at Ext_err_valid: error matrix: error("
00110 << i << "," << i << ")= " << elem << " > " << Large << std::endl;
00111 }
00112 }
00113 }
00114 if( !trace ) valid = 0;
00115 return( valid );
00116 }