21 #include "GaudiKernel/MsgStream.h"
27 static const double DELTA = 0.001;
36 if (
this != &mfeature) {
61 int a =
static_cast<int>(
m_flag/10000);
67 int a =
static_cast<int>(
m_flag/10000);
68 int b =
static_cast<int>(
m_flag-
a*10000);
79 m_sumet01 ( feat->getsumet01() ),
80 m_sumet02 ( feat->getsumet02() ),
81 m_sumet03 ( feat->getsumet03() ),
82 m_sumet04 ( feat->getsumet04() ),
83 m_sumpt01 ( feat->getsumpt01() ),
84 m_sumpt02 ( feat->getsumpt02() ),
85 m_sumpt03 ( feat->getsumpt03() ),
86 m_sumpt04 ( feat->getsumpt04() ),
87 m_PtMuID ( feat->getPtMuID() ),
88 m_MaxPtID ( feat->getMaxPtID() ),
90 m_RoiIdMu ( feat->getRoiIdMu() ),
91 m_PtMu ( feat->
pt() ),
92 m_QMu ( feat->getQMu() ),
93 m_EtaMu ( feat->getEtaMu() ),
94 m_PhiMu ( feat->getPhiMu() ) {}
103 m_sumet01 ( feat.getsumet01() ),
104 m_sumet02 ( feat.getsumet02() ),
105 m_sumet03 ( feat.getsumet03() ),
106 m_sumet04 ( feat.getsumet04() ),
107 m_sumpt01 ( feat.getsumpt01() ),
108 m_sumpt02 ( feat.getsumpt02() ),
109 m_sumpt03 ( feat.getsumpt03() ),
110 m_sumpt04 ( feat.getsumpt04() ),
111 m_PtMuID ( feat.getPtMuID() ),
112 m_MaxPtID ( feat.getMaxPtID() ),
114 m_RoiIdMu ( feat.getRoiIdMu() ),
115 m_PtMu ( feat.
pt() ),
116 m_QMu ( feat.getQMu() ),
117 m_EtaMu ( feat.getEtaMu() ),
118 m_PhiMu ( feat.getPhiMu() ) {}
121 std::cout << std::endl;
122 std::cout <<
"m_sumet01: " <<
m_sumet01 << std::endl;
123 std::cout <<
"m_sumet02: " <<
m_sumet02 << std::endl;
124 std::cout <<
"m_sumet03: " <<
m_sumet03 << std::endl;
125 std::cout <<
"m_sumet04: " <<
m_sumet04 << std::endl;
126 std::cout <<
"m_sumpt01: " <<
m_sumpt01 << std::endl;
127 std::cout <<
"m_sumpt02: " <<
m_sumpt02 << std::endl;
128 std::cout <<
"m_sumpt03: " <<
m_sumpt03 << std::endl;
129 std::cout <<
"m_sumpt04: " <<
m_sumpt04 << std::endl;
130 std::cout <<
"m_PtMuID: " <<
m_PtMuID << std::endl;
131 std::cout <<
"m_MaxPtID: " <<
m_MaxPtID << std::endl;
132 std::cout <<
"m_flag: " <<
m_flag << std::endl;
133 std::cout <<
"m_RoiIdMu: " <<
m_RoiIdMu << std::endl;
134 std::cout <<
"m_PtMu: " <<
m_PtMu << std::endl;
135 std::cout <<
"m_QMu: " <<
m_QMu << std::endl;
136 std::cout <<
"m_EtaMu: " <<
m_EtaMu << std::endl;
137 std::cout <<
"m_PhiMu: " <<
m_PhiMu << std::endl;
138 std::cout << std::endl;
165 std::stringstream
ss;
166 ss <<
"sumet01: " <<
d.getsumet01()
167 <<
"; sumet02: " <<
d.getsumet02()
168 <<
"; sumet03: " <<
d.getsumet03()
169 <<
"; sumet04: " <<
d.getsumet04()
170 <<
"; sumpt01: " <<
d.getsumpt01()
171 <<
"; sumpt02: " <<
d.getsumpt02()
172 <<
"; sumpt03: " <<
d.getsumpt03()
173 <<
"; sumpt04: " <<
d.getsumpt04()
174 <<
"; PtMuID: " <<
d.getPtMuID()
175 <<
"; MaxPtID: " <<
d.getMaxPtID()
176 <<
"; flag: " <<
d.getFlag()
177 <<
"; RoiIdMu: " <<
d.getRoiIdMu()
178 <<
"; RoiIdMu: " <<
d.getRoiIdMu()
179 <<
"; PtMu: " <<
d.pt()
180 <<
"; QMu: " <<
d.getQMu()
181 <<
"; EtaMu: " <<
d.getEtaMu()
182 <<
"; PhiMu: " <<
d.getPhiMu();
189 return (
m <<
str(
d ) );
194 if( std::abs(
a.getsumet01() -
b.getsumet01() ) > DELTA )
196 if( std::abs(
a.getsumet02() -
b.getsumet02() ) > DELTA )
198 if( std::abs(
a.getsumet03() -
b.getsumet03() ) > DELTA )
200 if( std::abs(
a.getsumet04() -
b.getsumet04() ) > DELTA )
202 if( std::abs(
a.getsumpt01() -
b.getsumpt01() ) > DELTA )
204 if( std::abs(
a.getsumpt02() -
b.getsumpt02() ) > DELTA )
206 if( std::abs(
a.getsumpt03() -
b.getsumpt03() ) > DELTA )
208 if( std::abs(
a.getsumpt04() -
b.getsumpt04() ) > DELTA )
210 if( std::abs(
a.getPtMuID() -
b.getPtMuID() ) > DELTA )
212 if( std::abs(
a.getMaxPtID() -
b.getMaxPtID() ) > DELTA )
214 if(
a.getFlag() !=
b.getFlag() )
216 if(
a.getRoiIdMu() !=
b.getRoiIdMu() )
218 if( std::abs(
a.pt() -
b.pt() ) > DELTA )
220 if( std::abs(
a.getQMu() -
b.getQMu() ) > DELTA )
222 if( std::abs(
a.getEtaMu() -
b.getEtaMu() ) > DELTA )
224 if( std::abs(
a.getPhiMu() -
b.getPhiMu() ) > DELTA )
240 std::map< std::string, double >& variableChange ) {
242 if( std::abs(
a.getsumet01() -
b.getsumet01() ) > DELTA ) {
243 variableChange[
"sumet01" ] =
a.getsumet01() -
b.getsumet01();
245 if( std::abs(
a.getsumet02() -
b.getsumet02() ) > DELTA ) {
246 variableChange[
"sumet02" ] =
a.getsumet02() -
b.getsumet02();
248 if( std::abs(
a.getsumet03() -
b.getsumet03() ) > DELTA ) {
249 variableChange[
"sumet03" ] =
a.getsumet03() -
b.getsumet03();
251 if( std::abs(
a.getsumet04() -
b.getsumet04() ) > DELTA ) {
252 variableChange[
"sumet04" ] =
a.getsumet04() -
b.getsumet04();
254 if( std::abs(
a.getsumpt01() -
b.getsumpt01() ) > DELTA ) {
255 variableChange[
"sumpt01" ] =
a.getsumpt01() -
b.getsumpt01();
257 if( std::abs(
a.getsumet02() -
b.getsumet02() ) > DELTA ) {
258 variableChange[
"sumet02" ] =
a.getsumet02() -
b.getsumpt02();
260 if( std::abs(
a.getsumpt03() -
b.getsumpt03() ) > DELTA ) {
261 variableChange[
"sumpt03" ] =
a.getsumpt03() -
b.getsumpt03();
263 if( std::abs(
a.getsumpt04() -
b.getsumpt04() ) > DELTA ) {
264 variableChange[
"sumpt04" ] =
a.getsumpt04() -
b.getsumpt04();
266 if(
a.getFlag() !=
b.getFlag() ) {
267 variableChange[
"flag" ] =
static_cast< double >(
a.getFlag() -
270 if(
a.getRoiIdMu() !=
b.getRoiIdMu() ) {
271 variableChange[
"RoiIdMu" ] =
static_cast< double >(
a.getRoiIdMu() -
274 if( std::abs(
a.pt() -
b.pt() ) > DELTA ) {
275 variableChange[
"PtMu" ] =
a.pt() -
b.pt();
277 if( std::abs(
a.getQMu() -
b.getQMu() ) > DELTA ) {
278 variableChange[
"QMu" ] =
a.getQMu() -
b.getQMu();
280 if( std::abs(
a.getEtaMu() -
b.getEtaMu() ) > DELTA ) {
281 variableChange[
"EtaMu" ] =
a.getEtaMu() -
b.getEtaMu();
283 if( std::abs(
a.getPhiMu() -
b.getPhiMu() ) > DELTA ) {
284 variableChange[
"PhiMu" ] =
a.getPhiMu() -
b.getPhiMu();