24 {
25
26
27
28 SG::WriteDecorHandle< xAOD::JetContainer, float > decoratorD2(
m_decorationName, ctx);
29
30
32
33 if(!largeRjets.isValid()) {
35 return StatusCode::FAILURE;
36 }
37
38
39 for( const auto *jet: *largeRjets){
40
41 float ecf1 = jet->getAttribute<float>("ECF1");
42 float ecf2 = jet->getAttribute<float>("ECF2");
43 float ecf3 = jet->getAttribute<float>("ECF3");
44
45
46 float D2=-999;
47 if(std::abs(ecf2)>1e-8)
48 D2=ecf3 * std::pow(ecf1, 3.0) / std::pow(ecf2, 3.0);
49 decoratorD2(*jet) = D2;
50 }
51
52 return StatusCode::SUCCESS;
53 }
SG::ReadHandleKey< xAOD::JetContainer > m_jetContainerKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decorationName