ATLAS Offline Software
Trigger
TrigT1
L1Topo
L1TopoAlgorithms
Root
gJetSort.cxx
Go to the documentation of this file.
1
/*
2
Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
3
*/
4
// gJetSort.cxx
5
// TopoCore
6
// algorithm to make sorted gJets lists
7
//
8
#include "
L1TopoAlgorithms/gJetSort.h
"
9
#include "
L1TopoEvent/TOBArray.h
"
10
#include "
L1TopoEvent/gJetTOBArray.h
"
11
#include "
L1TopoEvent/GenericTOB.h
"
12
#include <algorithm>
13
14
REGISTER_ALG_TCS
(gJetSort)
15
16
bool
SortByEtLargestgJet
(
TCS
::GenericTOB* tob1,
TCS
::GenericTOB* tob2)
17
{
18
return
tob1->Et() > tob2->Et();
19
}
20
21
22
// constructor
23
TCS::gJetSort::gJetSort
(
const
std::string &
name
) :
24
SortingAlg
(
name
)
25
{
26
defineParameter
(
"InputWidth"
, 64 );
// for FW
27
defineParameter
(
"InputWidth1stStage"
, 16 );
// for FW
28
defineParameter
(
"OutputWidth"
, 10 );
29
defineParameter
(
"MinEta"
, 0 );
30
defineParameter
(
"MaxEta"
, 196 );
31
}
32
33
34
TCS::gJetSort::~gJetSort
()
35
{}
36
37
38
39
TCS::StatusCode
40
TCS::gJetSort::initialize
() {
41
m_numberOfJets = parameter(
"OutputWidth"
).value();
42
m_minEta = parameter(
"MinEta"
).value();
43
m_maxEta = parameter(
"MaxEta"
).value();
44
return
TCS::StatusCode::SUCCESS
;
45
}
46
47
48
TCS::StatusCode
49
TCS::gJetSort::sort
(
const
InputTOBArray
&
input
,
TOBArray
&
output
) {
50
const
gJetTOBArray
&
jets
=
dynamic_cast<
const
gJetTOBArray
&
>
(
input
);
51
52
// fill output array with GenericTOBs builds from jets
53
for
(
gJetTOBArray::const_iterator
jet
=
jets
.begin();
jet
!=
jets
.end(); ++
jet
) {
54
if
(
parType_t
(std::abs((*jet)->
eta
())) < m_minEta)
continue
;
55
if
(
parType_t
(std::abs((*jet)->
eta
())) > m_maxEta)
continue
;
56
output
.push_back(
GenericTOB
(**
jet
) );
57
}
58
59
// sort
60
output
.sort(
SortByEtLargestgJet
);
61
62
63
// keep only max number of jets
64
int
par
= m_numberOfJets;
65
unsigned
int
maxNumberOfJets = std::clamp(
par
, 0, std::abs(
par
));
66
if
(maxNumberOfJets>0) {
67
while
(
output
.size()> maxNumberOfJets ) {
68
if
(
output
.size() == (maxNumberOfJets+1)) {
69
bool
isAmbiguous =
output
[maxNumberOfJets-1].EtDouble() ==
output
[maxNumberOfJets].EtDouble();
70
if
(isAmbiguous) {
output
.setAmbiguityFlag(
true
); }
71
}
72
output
.pop_back();
73
}
74
}
75
return
TCS::StatusCode::SUCCESS
;
76
}
77
TCS::gJetSort::~gJetSort
virtual ~gJetSort()
Definition:
gJetSort.cxx:34
TCS::StatusCode::SUCCESS
@ SUCCESS
Definition:
Trigger/TrigT1/L1Topo/L1TopoCommon/L1TopoCommon/StatusCode.h:17
TCS::parType_t
uint32_t parType_t
Definition:
Parameter.h:22
TCS::DataArrayImpl< gJetTOB >::const_iterator
data_t::const_iterator const_iterator
Definition:
DataArrayImpl.h:18
TCS::gJetSort::sort
virtual TCS::StatusCode sort(const InputTOBArray &input, TOBArray &output) override final
Definition:
gJetSort.cxx:49
eta
Scalar eta() const
pseudorapidity method
Definition:
AmgMatrixBasePlugin.h:79
jet
Definition:
JetCalibTools_PlotJESFactors.cxx:23
TCS::ConfigurableAlg::defineParameter
void defineParameter(const std::string &name, TCS::parType_t value)
Definition:
ConfigurableAlg.cxx:201
TCS::InputTOBArray
Definition:
InputTOBArray.h:15
PlotPulseshapeFromCool.input
input
Definition:
PlotPulseshapeFromCool.py:106
gJetSort.h
TCS::gJetSort::gJetSort
gJetSort(const std::string &name)
Definition:
gJetSort.cxx:23
TOBArray.h
TCS::TOBArray
Definition:
TOBArray.h:24
TCS::gJetSort::initialize
virtual TCS::StatusCode initialize() override
Definition:
gJetSort.cxx:40
merge.output
output
Definition:
merge.py:17
TCS::GenericTOB
Definition:
GenericTOB.h:35
REGISTER_ALG_TCS
#define REGISTER_ALG_TCS(CLASS)
Definition:
AlgFactory.h:62
name
std::string name
Definition:
Control/AthContainers/Root/debug.cxx:192
createCoolChannelIdFile.par
par
Definition:
createCoolChannelIdFile.py:29
TCS::SortingAlg
Definition:
SortingAlg.h:18
GenericTOB.h
TCS
Definition:
Global/GlobalSimulation/src/IO/Decision.h:18
TCS::gJetTOBArray
Definition:
gJetTOBArray.h:19
defineDB.jets
list jets
Definition:
JetTagCalibration/share/defineDB.py:24
gJetTOBArray.h
xAOD::bool
setBGCode setTAP setLVL2ErrorBits bool
Definition:
TrigDecision_v1.cxx:60
TCS::StatusCode
Definition:
Trigger/TrigT1/L1Topo/L1TopoCommon/L1TopoCommon/StatusCode.h:15
SortByEtLargestgJet
bool SortByEtLargestgJet(TCS::GenericTOB *tob1, TCS::GenericTOB *tob2)
Definition:
gJetSort.cxx:16
Generated on Sun Jun 30 2024 21:16:27 for ATLAS Offline Software by
1.8.18