26 defineParameter(
"InputWidth", 3);
27 defineParameter(
"MaxTob", 0);
28 defineParameter(
"NumResultBits", 3);
29 defineParameter(
"MinET1",1,0);
30 defineParameter(
"MinET2",1,0);
31 defineParameter(
"MinDeltaEta", 0, 0);
32 defineParameter(
"MaxDeltaEta", 49, 0);
33 defineParameter(
"MinDeltaPhi", 0, 0);
34 defineParameter(
"MaxDeltaPhi", 63, 0);
35 defineParameter(
"MinET1",1,1);
36 defineParameter(
"MinET2",1,1);
37 defineParameter(
"MinDeltaEta", 0, 1);
38 defineParameter(
"MaxDeltaEta", 49, 1);
39 defineParameter(
"MinDeltaPhi", 0, 1);
40 defineParameter(
"MaxDeltaPhi", 63, 1);
41 defineParameter(
"MinET1",1,2);
42 defineParameter(
"MinET2",1,2);
43 defineParameter(
"MinDeltaEta", 0, 2);
44 defineParameter(
"MaxDeltaEta", 49, 2);
45 defineParameter(
"MinDeltaPhi", 0, 2);
46 defineParameter(
"MaxDeltaPhi", 63, 2);
47 setNumberOutputBits(3);
56 if(parameter(
"MaxTob").
value() > 0) {
57 p_NumberLeading1 = parameter(
"MaxTob").value();
58 p_NumberLeading2 = parameter(
"MaxTob").value();
60 p_NumberLeading1 = parameter(
"InputWidth").value();
61 p_NumberLeading2 = parameter(
"InputWidth").value();
64 for(
unsigned int i=0;
i<numberOutputBits(); ++
i) {
65 p_DeltaPhiMin[
i] = parameter(
"MinDeltaPhi",
i).value();
66 p_DeltaPhiMax[
i] = parameter(
"MaxDeltaPhi",
i).value();
67 p_DeltaEtaMin[
i] = parameter(
"MinDeltaEta",
i).value();
68 p_DeltaEtaMax[
i] = parameter(
"MaxDeltaEta",
i).value();
69 p_MinET1[
i] = parameter(
"MinET1",
i).value();
70 p_MinET2[
i] = parameter(
"MinET2",
i).value();
85 for(
unsigned int i=0;
i<numberOutputBits(); ++
i) {
86 std::string hname_accept =
"hDeltaEtaPhiIncl1_accept_bit"+
std::to_string((
int)
i);
87 std::string hname_reject =
"hDeltaEtaPhiIncl1_reject_bit"+
std::to_string((
int)
i);
89 bookHist(m_histAccept, hname_accept,
"DETA vs DPHI", 100, p_DeltaEtaMin[
i], p_DeltaEtaMax[
i], 100, p_DeltaPhiMin[
i], p_DeltaPhiMax[
i]);
90 bookHist(m_histReject, hname_reject,
"DETA vs DPHI", 100, p_DeltaEtaMin[
i], p_DeltaEtaMax[
i], 100, p_DeltaPhiMin[
i], p_DeltaPhiMax[
i]);
100 const std::vector<TCS::TOBArray *> & output,
103 if(input.size() == 1) {
105 unsigned int nLeading = p_NumberLeading1;
107 tob1 != input[0]->end() &&
distance( input[0]->
begin(), tob1) < nLeading;
112 tob2 != input[0]->end() &&
distance( input[0]->
begin(), tob2) < p_NumberLeading2;
115 unsigned int deltaPhi = calcDeltaPhiBW( *tob1, *tob2 );
117 unsigned int deltaEta = calcDeltaEtaBW( *tob1, *tob2 );
120 std::stringstream msgss;
121 msgss <<
" phi1=" << (*tob1)->phi() <<
" , phi2=" << (*tob2)->phi()
123 for(
unsigned int i=0;
i<numberOutputBits(); ++
i) {
129 const bool fillAccept = fillHistos() and (fillHistosBasedOnHardware() ? getDecisionHardwareBit(
i) :
accept);
130 const bool fillReject = fillHistos() and not fillAccept;
131 const bool alreadyFilled = decision.
bit(
i);
136 if(fillAccept and not alreadyFilled) {
138 }
else if(fillReject) {
141 msgss <<
"DeltaEtaPhiIncl1 alg bit" <<
i << (
accept?
" pass":
" fail") <<
"|";
146 for (
unsigned int i=0;
i < numberOutputBits(); ++
i) {
149 output[
i]->setAmbiguityFlag(hasAmbiguousInputs);
152 TCS_EXCEPTION(
"DeltaEtaPhiIncl1 alg must have 1 input, but got " << input.size());
159 const std::vector<TCS::TOBArray *> & output,
162 if(input.size() == 1) {
164 unsigned int nLeading = p_NumberLeading1;
166 tob1 != input[0]->end() &&
distance( input[0]->
begin(), tob1) < nLeading;
171 tob2 != input[0]->end() &&
distance( input[0]->
begin(), tob2) < p_NumberLeading2;
174 unsigned int deltaPhi = calcDeltaPhi( *tob1, *tob2 );
176 unsigned int deltaEta = calcDeltaEta( *tob1, *tob2 );
179 std::stringstream msgss;
180 msgss <<
" Combination : " <<
distance( input[0]->
begin(), tob1)
182 <<
" phi1=" << (*tob1)->phi()
183 <<
" , phi2=" << (*tob2)->phi()
185 <<
", DeltaEta = " <<
deltaEta <<
" -> ";
186 for(
unsigned int i=0;
i<numberOutputBits(); ++
i) {
192 const bool fillAccept = fillHistos() and (fillHistosBasedOnHardware() ? getDecisionHardwareBit(
i) :
accept);
193 const bool fillReject = fillHistos() and not fillAccept;
194 const bool alreadyFilled = decision.
bit(
i);
199 if(fillAccept and not alreadyFilled) {
201 }
else if(fillReject) {
204 msgss <<
"DeltaEtaPhiIncl1 alg bit" <<
i << (
accept?
" pass":
" fail") <<
"|";
210 TCS_EXCEPTION(
"DeltaEtaPhiIncl1 alg must have 1 input, but got " << input.size());