#include <PionBinData.h>
Inheritance diagram for Atlfast::PionBinData:
Public Member Functions | |
| PionBinData (BinID &, vector< ParameterResolutions * >, vector< ParameterResolutions * >, vector< ParameterResolutions * >, vector< ParameterResolutions * >, int randSeed) | |
| virtual | ~PionBinData () |
| HepSymMatrix | getMatrix (const TrackTrajectory &track) const |
Private Attributes | |
| BinID | m_id |
| vector< ParameterResolutions * > | m_cores |
| vector< ParameterResolutions * > | m_tails |
| vector< ParameterResolutions * > | m_fractions |
| vector< ParameterResolutions * > | m_correlations |
| HepRandomEngine * | m_randomEngine |
Definition at line 34 of file PionBinData.h.
| Atlfast::PionBinData::PionBinData | ( | BinID & | , | |
| vector< ParameterResolutions * > | , | |||
| vector< ParameterResolutions * > | , | |||
| vector< ParameterResolutions * > | , | |||
| vector< ParameterResolutions * > | , | |||
| int | randSeed | |||
| ) |
Constructor from flat file contents
Definition at line 9 of file PionBinData.cxx.
00015 : 00016 m_id(id), 00017 m_cores(core), 00018 m_tails(tails), 00019 m_fractions(fractions), 00020 m_correlations(correlations) 00021 { 00022 m_randomEngine = new HepJamesRandom(randSeed); 00023 }
| Atlfast::PionBinData::~PionBinData | ( | ) | [virtual] |
returns the correlation matrix depending on an internally generated random number and track trajectory supplied as argument
Definition at line 26 of file PionBinData.cxx.
00027 { 00028 delete m_randomEngine; 00029 for (size_t i=0; i < m_fractions.size(); i++) 00030 delete m_fractions[i]; 00031 for (size_t i=0; i < m_correlations.size(); i++) 00032 delete m_correlations[i]; 00033 for (size_t i=0; i < m_cores.size(); i++) 00034 delete m_cores[i]; 00035 for (size_t i=0; i < m_tails.size(); i++) 00036 delete m_tails[i]; 00037 }
| HepSymMatrix Atlfast::PionBinData::getMatrix | ( | const TrackTrajectory & | track | ) | const [virtual] |
returns the correlation matrix depending on a random number supplied as argument
Implements Atlfast::IBinData.
Definition at line 58 of file PionBinData.cxx.
00059 { 00060 HepSymMatrix Sigma(5,0); 00061 00062 double fraction, random[5]; 00063 00064 // order of parameters: d0, z0, phi0, cot(theta0), q/pT 00065 // introduce "correlation" between core <--> tails 00066 random[0] = m_randomEngine->flat(); 00067 random[1] = m_randomEngine->flat(); 00068 random[2] = random[0]; 00069 random[3] = random[1]; 00070 random[4] = random[0]; 00071 00072 // diagonals 00073 for ( int param = 0; param < 5; param++ ) 00074 { 00075 fraction = m_fractions[param]->resolution(traj); 00076 if ( fraction > 1.0 ) fraction = 1.0; 00077 Sigma[param][param] = ( random[param] < fraction ) ? 00078 std::pow( m_cores[param]->resolution(traj), 2 ) : 00079 std::pow( m_tails[param]->resolution(traj), 2 ); 00080 } 00081 00082 // off-diagonals 00083 // NOTE: m_correlations[] holds correlation coefficients, need covariances 00084 00085 // (1,3) ... cov(d0,phi0) 00086 // (1,5) ... cov(d0,q/pT) 00087 // (3,5) ... cov(phi0,q/pT) 00088 double rho13 = m_correlations[0]->resolution(traj); 00089 double rho15 = m_correlations[1]->resolution(traj); 00090 double rho35 = m_correlations[2]->resolution(traj); 00091 00092 // covariance sub-matrix of transverse parameters needs to be positive definite 00093 // in order that its square root (cf. PionMatrixManager) exists 00094 double det3 = 1 - rho13 * rho13 - rho15 * rho15 - rho35 * rho35 - 2 * rho13 * rho15 * rho35; 00095 if ( det3 < 0 ) rho13 = rho15 = rho35 = 0; 00096 00097 // make sure that correlation coefficients stay within [-1,+1] 00098 if ( std::abs(rho13) > 1 ) rho13 *= 0.99 / std::abs(rho13); 00099 if ( std::abs(rho15) > 1 ) rho15 *= 0.99 / std::abs(rho15); 00100 if ( std::abs(rho35) > 1 ) rho35 *= 0.99 / std::abs(rho35); 00101 00102 Sigma(1,3) = Sigma(3,1) = rho13 * std::sqrt( Sigma(1,1) * Sigma(3,3) ); 00103 Sigma(1,5) = Sigma(5,1) = rho15 * std::sqrt( Sigma(1,1) * Sigma(5,5) ); 00104 Sigma(3,5) = Sigma(5,3) = rho35 * std::sqrt( Sigma(3,3) * Sigma(5,5) ); 00105 00106 00107 // (2,4) ... cov(z0,cot(theta0)) 00108 double rho24 = m_correlations[3]->resolution(traj); 00109 // make sure that correlation coefficient stays within [-1,+1] 00110 if ( std::abs(rho24) > 1 ) rho24 *= 0.99 / std::abs(rho24); 00111 Sigma(2,4) = Sigma(4,2) = rho24 * std::sqrt( Sigma(2,2) * Sigma(4,4) ); 00112 00113 // DONE! 00114 return Sigma; 00115 00116 }
BinID Atlfast::PionBinData::m_id [private] |
Definition at line 58 of file PionBinData.h.
vector< ParameterResolutions* > Atlfast::PionBinData::m_cores [private] |
Definition at line 60 of file PionBinData.h.
vector< ParameterResolutions* > Atlfast::PionBinData::m_tails [private] |
Definition at line 61 of file PionBinData.h.
vector< ParameterResolutions* > Atlfast::PionBinData::m_fractions [private] |
Definition at line 62 of file PionBinData.h.
vector< ParameterResolutions* > Atlfast::PionBinData::m_correlations [private] |
Definition at line 63 of file PionBinData.h.
HepRandomEngine* Atlfast::PionBinData::m_randomEngine [private] |
Definition at line 65 of file PionBinData.h.
1.5.1