30 if( alphaFix < 0.0 ) {
44 return StatusCode::SUCCESS;
52 std::string
suffix = moment.second.suffix;
54 if( !moment.second.acc_Tau1->isAvailable(
jet) ||
55 !moment.second.acc_Tau2->isAvailable(
jet) ||
56 !moment.second.acc_Tau3->isAvailable(
jet) ||
57 !moment.second.acc_Tau4->isAvailable(
jet) ||
58 !moment.second.acc_Tau1_wta->isAvailable(
jet) ||
59 !moment.second.acc_Tau2_wta->isAvailable(
jet) ||
60 !moment.second.acc_Tau3_wta->isAvailable(
jet) ||
61 !moment.second.acc_Tau4_wta->isAvailable(
jet)) {
70 if( !moment.second.acc_Tau2_ungroomed->isAvailable(
jet) ||
71 !moment.second.acc_Tau3_ungroomed->isAvailable(
jet) ||
72 !moment.second.acc_Tau4_ungroomed->isAvailable(
jet) ||
73 !moment.second.acc_Tau2_wta_ungroomed->isAvailable(
jet) ||
74 !moment.second.acc_Tau3_wta_ungroomed->isAvailable(
jet) ||
75 !moment.second.acc_Tau4_wta_ungroomed->isAvailable(
jet)) {
85 float Tau1 = (*moment.second.acc_Tau1)(
jet);
86 float Tau2 = (*moment.second.acc_Tau2)(
jet);
87 float Tau3 = (*moment.second.acc_Tau3)(
jet);
88 float Tau4 = (*moment.second.acc_Tau4)(
jet);
90 float Tau2_ungroomed = -999.0;
91 float Tau3_ungroomed = -999.0;
92 float Tau4_ungroomed = -999.0;
95 Tau2_ungroomed = (*moment.second.acc_Tau2_ungroomed)(
jet);
96 Tau3_ungroomed = (*moment.second.acc_Tau3_ungroomed)(
jet);
97 Tau4_ungroomed = (*moment.second.acc_Tau4_ungroomed)(
jet);
100 float Tau21 = -999.0;
101 float Tau32 = -999.0;
102 float Tau42 = -999.0;
104 float Tau21_dichroic = -999.0;
105 float Tau32_dichroic = -999.0;
106 float Tau42_dichroic = -999.0;
113 if( Tau2_ungroomed > 1
e-8 ) {
114 Tau21_dichroic = Tau2_ungroomed/Tau1;
125 if(Tau3_ungroomed > 1
e-8) {
126 Tau32_dichroic = Tau3_ungroomed / Tau2;
129 if(Tau4_ungroomed > 1
e-8) {
130 Tau42_dichroic = Tau4_ungroomed / Tau2;
136 float Tau1_wta = (*moment.second.acc_Tau1_wta)(
jet);
137 float Tau2_wta = (*moment.second.acc_Tau2_wta)(
jet);
138 float Tau3_wta = (*moment.second.acc_Tau3_wta)(
jet);
139 float Tau4_wta = (*moment.second.acc_Tau4_wta)(
jet);
141 float Tau2_wta_ungroomed = -999.0;
142 float Tau3_wta_ungroomed = -999.0;
143 float Tau4_wta_ungroomed = -999.0;
146 Tau2_wta_ungroomed = (*moment.second.acc_Tau2_wta_ungroomed)(
jet);
147 Tau3_wta_ungroomed = (*moment.second.acc_Tau3_wta_ungroomed)(
jet);
148 Tau4_wta_ungroomed = (*moment.second.acc_Tau4_wta_ungroomed)(
jet);
151 float Tau21_wta = -999.0;
152 float Tau32_wta = -999.0;
153 float Tau42_wta = -999.0;
155 float Tau21_wta_dichroic = -999.0;
156 float Tau32_wta_dichroic = -999.0;
157 float Tau42_wta_dichroic = -999.0;
160 if( Tau1_wta > 1
e-8 ) {
162 Tau21_wta = Tau2_wta / Tau1_wta;
164 if( Tau2_wta_ungroomed > 1
e-8 ) {
165 Tau21_wta_dichroic = Tau2_wta_ungroomed / Tau1_wta;
171 if( Tau2_wta > 1
e-8 ) {
173 Tau32_wta = Tau3_wta / Tau2_wta;
174 Tau42_wta = Tau4_wta / Tau2_wta;
176 if( Tau3_wta_ungroomed > 1
e-8 ) {
177 Tau32_wta_dichroic = Tau3_wta_ungroomed / Tau2_wta;
180 if(Tau4_wta_ungroomed > 1
e-8) {
181 Tau42_wta_dichroic = Tau4_wta_ungroomed / Tau2_wta;
186 (*moment.second.dec_Tau21)(
jet) =
Tau21;
187 (*moment.second.dec_Tau32)(
jet) =
Tau32;
188 (*moment.second.dec_Tau42)(
jet) = Tau42;
190 (*moment.second.dec_Tau21_dichroic)(
jet) = Tau21_dichroic;
191 (*moment.second.dec_Tau32_dichroic)(
jet) = Tau32_dichroic;
192 (*moment.second.dec_Tau42_dichroic)(
jet) = Tau42_dichroic;
194 (*moment.second.dec_Tau21_wta)(
jet) = Tau21_wta;
195 (*moment.second.dec_Tau32_wta)(
jet) = Tau32_wta;
196 (*moment.second.dec_Tau42_wta)(
jet) = Tau42_wta;
198 (*moment.second.dec_Tau21_wta_dichroic)(
jet) = Tau21_wta_dichroic;
199 (*moment.second.dec_Tau32_wta_dichroic)(
jet) = Tau32_wta_dichroic;
200 (*moment.second.dec_Tau42_wta_dichroic)(
jet) = Tau42_wta_dichroic;