|
|
|
#include <L1RCTElectronIsolationCard.h>
Public Member Functions | |
| int | cardNumber () |
| int | crateNumber () |
| void | fillElectronCandidates () |
| unsigned short | getIsoElectrons (int i) |
| unsigned short | getNonIsoElectrons (int i) |
| L1RCTElectronIsolationCard (int crateNumber, int cardNumber, const L1RCTLookupTables *rctLookupTables) | |
| void | print () |
| void | printEdges () |
| void | setRegion (int i, const L1RCTRegion ®ion) |
| ~L1RCTElectronIsolationCard () | |
Private Member Functions | |
| std::vector< unsigned short > | calcElectronCandidates (const L1RCTRegion ®ion, int regionNum) |
| unsigned short | calcMaxSum (unsigned short primaryEt, unsigned short northEt, unsigned short southEt, unsigned short eastEt, unsigned short westEt) |
| L1RCTElectronIsolationCard () | |
Private Attributes | |
| unsigned short | cardNo |
| unsigned short | crtNo |
| L1RCTRegion | empty |
| std::vector< unsigned short > | isoElectrons |
| std::vector< unsigned short > | nonIsoElectrons |
| const L1RCTLookupTables * | rctLookupTables_ |
| std::vector< L1RCTRegion > | regions |
Definition at line 19 of file L1RCTElectronIsolationCard.h.
| L1RCTElectronIsolationCard::L1RCTElectronIsolationCard | ( | int | crateNumber, |
| int | cardNumber, | ||
| const L1RCTLookupTables * | rctLookupTables | ||
| ) |
Definition at line 8 of file L1RCTElectronIsolationCard.cc.
References regions.
: crtNo(crateNumber),cardNo(cardNumber), rctLookupTables_(rctLookupTables), isoElectrons(2),nonIsoElectrons(2), regions(2) { regions.push_back(L1RCTRegion()); regions.push_back(L1RCTRegion()); }
| L1RCTElectronIsolationCard::~L1RCTElectronIsolationCard | ( | ) |
Definition at line 19 of file L1RCTElectronIsolationCard.cc.
References regions.
{
regions.clear();
}
| L1RCTElectronIsolationCard::L1RCTElectronIsolationCard | ( | ) | [private] |
| std::vector< unsigned short > L1RCTElectronIsolationCard::calcElectronCandidates | ( | const L1RCTRegion & | region, |
| int | regionNum | ||
| ) | [private] |
Definition at line 44 of file L1RCTElectronIsolationCard.cc.
References calcMaxSum(), cardNo, cardNumber(), benchmark_cfg::cerr, crateNumber(), L1RCTParameters::eicIsolationThreshold(), L1RCTRegion::getEtIn7Bits(), L1RCTRegion::getHE_FGBit(), i, j, n, python::Vispa::Plugins::EdmBrowser::EdmDataAccessor::ne(), rctLookupTables_, L1RCTLookupTables::rctParameters(), and asciidump::s.
Referenced by fillElectronCandidates().
{
unsigned short nonIsoElectron = 0;
unsigned short isoElectron = 0;
//i is row and j is column
for(int i = 0; i<4; i++){
for(int j = 0; j<4; j++){
unsigned short primaryEt = region.getEtIn7Bits(i,j);
unsigned short primaryHE_FG = region.getHE_FGBit(i,j);
unsigned short northEt = region.getEtIn7Bits(i-1, j);
unsigned short southEt = region.getEtIn7Bits(i+1, j);
unsigned short westEt = region.getEtIn7Bits( i,j-1);
unsigned short eastEt = region.getEtIn7Bits( i,j+1);
unsigned short neEt = region.getEtIn7Bits(i-1,j+1);
unsigned short nwEt = region.getEtIn7Bits(i-1,j-1);
unsigned short seEt = region.getEtIn7Bits(i+1,j+1);
unsigned short swEt = region.getEtIn7Bits(i+1,j-1);
unsigned short northHE_FG = region.getHE_FGBit(i-1, j);
unsigned short southHE_FG = region.getHE_FGBit(i+1, j);
unsigned short westHE_FG = region.getHE_FGBit( i,j-1);
unsigned short eastHE_FG = region.getHE_FGBit( i,j+1);
unsigned short neHE_FG = region.getHE_FGBit(i-1,j+1);
unsigned short nwHE_FG = region.getHE_FGBit(i-1,j-1);
unsigned short seHE_FG = region.getHE_FGBit(i+1,j+1);
unsigned short swHE_FG = region.getHE_FGBit(i+1,j-1);
bool top = false;
int nCrate = crateNumber();
int nCard = cardNumber();
int nRegion = regionNum;
// top row of crate
if (nCard == 0 || nCard == 2 || nCard == 4 || (nCard == 6 && nRegion == 0))
{
top = true;
}
// bottom row of crate
else if (nCard == 1 || nCard == 3 || nCard == 5 || (nCard == 6 && nRegion == 1))
{} // top already false
else
{
std::cerr << "Error! EIC top assignment" << std::endl; // this shouldn't happen!
}
// The following values are used for zeroing and determining whether or
// not a tower is a "candidate". The original primaryEt, northEt, neEt, etc.
// are used to calculate vetoes and must not be zeroed.
unsigned short primaryTowerEt = primaryEt;
unsigned short northTowerEt = northEt;
unsigned short southTowerEt = southEt;
unsigned short eastTowerEt = eastEt;
unsigned short westTowerEt = westEt;
// In order to ensure proper selection of candidate tower and neighbor,
// if two neighbor energies are equal, one is set to zero (in the
// appropriate regions, those for which tp_lf bit is set to 0 in
// Pam's JCCTest/EGWithShare.cc).
if (primaryEt > 0) // this value should maybe be customizable?
{
if (nCard != 6) // all cards except 6
{
if(top && nCrate >= 9) // top row of regions in positive-eta crate
{
if(westTowerEt == eastTowerEt) westTowerEt=0;
if(southTowerEt == northTowerEt) southTowerEt=0;
if(southTowerEt == eastTowerEt) southTowerEt=0;
if(westTowerEt == northTowerEt) westTowerEt=0;
}
else if((!top) && nCrate < 9) // bottom row of regions in negative-eta crate
{
if(eastTowerEt == westTowerEt) eastTowerEt=0;
if(northTowerEt == southTowerEt) northTowerEt=0;
if(northTowerEt == westTowerEt) northTowerEt=0;
if(eastTowerEt == southTowerEt) eastTowerEt=0;
}
}
else // card 6
{
// only +eta card 6 needs to have zeroing. Pam sez.
// -eta card 6 does what it's supposed to even w/o zeroing.
if(nRegion == 0 && nCrate >=9)
{
if(westTowerEt == eastTowerEt) westTowerEt=0;
if(southTowerEt == northTowerEt) southTowerEt=0;
if(southTowerEt == eastTowerEt) southTowerEt=0;
if(westTowerEt == northTowerEt) westTowerEt=0;
}
}
}
// This section compares the energies in the primary tower with the
// surrounding towers to determine whether or not the primary tower
// should be considered a "candidate".
bool candidate = false;
// for case where primary tower et greater than all neighbors -> candidate
if (primaryEt > northEt && primaryEt > southEt && primaryEt > eastEt
&& primaryEt > westEt && !primaryHE_FG)
{
candidate = true;
}
// if primary et less than any neighbors (or HE_FG veto set) NOT a candidate!
else if (primaryEt < northEt || primaryEt < southEt || primaryEt < eastEt
|| primaryEt < westEt || primaryHE_FG)
{} // candidate already false
else // Case of primary tower et being equal to any of its neighbors.
// This section determines which tower gets the candidate.
// See AboutTP.pdf document, figure on p. 4, for clarification.
// Zeroed values are used in this calculation.
{
if (primaryEt > 0)
{
if (nCrate >=9) // positive eta
{
if (top) // top row of regions in crate. tp_lf == 0
// priority order: east < south < north < west
{
if (westTowerEt == primaryTowerEt)
candidate = true;
else if (northTowerEt == primaryTowerEt)
candidate = false;
else if (southTowerEt == primaryTowerEt)
candidate = true;
else if (eastTowerEt == primaryTowerEt)
candidate = false;
}
else // bottom row of regions in crate. tp_lf == 1
// priority order: west < north < south < east
{
if (eastTowerEt == primaryTowerEt)
candidate = true;
else if (southTowerEt == primaryTowerEt)
candidate = true;
else if (northTowerEt == primaryTowerEt)
candidate = false;
else if (westTowerEt == primaryTowerEt)
candidate = false;
if (nCard == 6) // card 6. tp_lf == 1
{
// priority order: east < north < south < west
if (westTowerEt == primaryTowerEt)
candidate = true;
else if (southTowerEt == primaryTowerEt)
candidate = true;
else if (northTowerEt == primaryTowerEt)
candidate = false;
else if (eastTowerEt == primaryTowerEt)
candidate = false;
}
}
}
else // negative eta
{
if (top) // top row of regions in crate. tp_lf == 1
// priority order: east < south < north < west
{
if (westTowerEt == primaryTowerEt)
candidate = true;
else if (northTowerEt == primaryTowerEt)
candidate = true;
else if (southTowerEt == primaryTowerEt)
candidate = false;
else if (eastTowerEt == primaryTowerEt)
candidate = false;
if (nCard == 6) // card 6. tp_lf == 0
// east < south < north < west
{
if (westTowerEt == primaryTowerEt)
candidate = false;
else if (northTowerEt == primaryTowerEt)
candidate = false;
else if (southTowerEt == primaryTowerEt)
candidate = true;
else if (eastTowerEt == primaryTowerEt)
candidate = true;
}
}
else // bottom row of regions. tp_lf == 0
// west < north < south < east
{
if (eastTowerEt == primaryTowerEt)
candidate = true;
else if (southTowerEt == primaryTowerEt)
candidate = false;
else if (northTowerEt == primaryTowerEt)
candidate = true;
else if (westTowerEt == primaryTowerEt)
candidate = false;
if (nCard == 6) // card 6. tp_lf == 1
// west < north < south < east
{
if (eastTowerEt == primaryTowerEt)
candidate = true;
else if (southTowerEt == primaryTowerEt)
candidate = true;
else if (northTowerEt == primaryTowerEt)
candidate = false;
else if (westTowerEt == primaryTowerEt)
candidate = false;
}
}
}
}
} // end of if (primary == neighbors)
if (candidate) {
// Either zeroed or non-zeroed set of values can be used here --
// neighbor tower only zeroed if another neighbor tower of same
// energy. Max sum calculated from primary and only one neighbor
// tower, and always one neighbor tower left over, so value of sum
// is not affected. Currently using non-zeroed.
unsigned short candidateEt = calcMaxSum(primaryEt,northEt,southEt,
eastEt,westEt);
// neighbor HE_FG veto true if neighbor has HE_FG set
bool neighborVeto = (nwHE_FG || northHE_FG || neHE_FG || westHE_FG ||
eastHE_FG || swHE_FG || southHE_FG || seHE_FG);
// threshold for five-tower corner quiet veto
//int quietThreshold = 3; // 3 - loose isolation 0 - very tight isolation
//int quietThreshold = 7; // ECALGREN
//int quietThreshold = 0; // HCALGREN
unsigned quietThreshold = rctLookupTables_->rctParameters()->eicIsolationThreshold();
bool nw = false;
bool ne = false;
bool sw = false;
bool se = false;
bool n = false;
bool w = false;
bool s = false;
bool e = false;
// individual neighbor vetoes set if neighbor is over threshold
if (nwEt >= quietThreshold) nw = true;
if (neEt >= quietThreshold) ne = true;
if (swEt >= quietThreshold) sw = true;
if (seEt >= quietThreshold) se = true;
if (northEt >= quietThreshold) n = true;
if (southEt >= quietThreshold) s = true;
if (westEt >= quietThreshold) w = true;
if (eastEt >= quietThreshold) e = true;
// veto TRUE for each corner set if any individual tower in each set is over threshold
bool nwC = (sw || w || nw || n || ne);
bool neC = (nw || n || ne || e || se);
bool seC = (ne || e || se || s || sw);
bool swC = (se || s || sw || w || nw);
// overall quiet veto TRUE only if NO corner sets are quiet
// (all are "loud") -> non-isolated
bool quietVeto = (nwC && neC && seC && swC);
// only isolated if both vetoes are false
// Note: quietThreshold = 0 forces all candidates to be non-iso
if(!(quietVeto || neighborVeto)){
if(candidateEt > isoElectron)
isoElectron = candidateEt;
}
// otherwise, non-isolated
else if(candidateEt > nonIsoElectron)
nonIsoElectron = candidateEt;
}
}
}
std::vector<unsigned short> candidates;
unsigned short fullIsoElectron = isoElectron*16 + cardNo*2; // leaves room for last bit -- region number, added in Crate.cc
candidates.push_back(fullIsoElectron);
unsigned short fullNonIsoElectron = nonIsoElectron*16 + cardNo*2; // leaves room for region info in last bit
candidates.push_back(fullNonIsoElectron);
return candidates;
}
| unsigned short L1RCTElectronIsolationCard::calcMaxSum | ( | unsigned short | primaryEt, |
| unsigned short | northEt, | ||
| unsigned short | southEt, | ||
| unsigned short | eastEt, | ||
| unsigned short | westEt | ||
| ) | [private] |
Definition at line 335 of file L1RCTElectronIsolationCard.cc.
References i, max(), and MultipleCompare::test.
Referenced by calcElectronCandidates().
| int L1RCTElectronIsolationCard::cardNumber | ( | ) | [inline] |
Definition at line 29 of file L1RCTElectronIsolationCard.h.
References cardNo.
Referenced by calcElectronCandidates().
{return cardNo;}
| int L1RCTElectronIsolationCard::crateNumber | ( | ) | [inline] |
Definition at line 28 of file L1RCTElectronIsolationCard.h.
References crtNo.
Referenced by calcElectronCandidates().
{return crtNo;}
| void L1RCTElectronIsolationCard::fillElectronCandidates | ( | ) |
Definition at line 25 of file L1RCTElectronIsolationCard.cc.
References calcElectronCandidates(), isoElectrons, nonIsoElectrons, and regions.
{
std::vector<unsigned short> region0Electrons = calcElectronCandidates(regions.at(0),0);
std::vector<unsigned short> region1Electrons = calcElectronCandidates(regions.at(1),1);
isoElectrons.at(0) = region0Electrons.at(0);
isoElectrons.at(1) = region1Electrons.at(0);
nonIsoElectrons.at(0) = region0Electrons.at(1);
nonIsoElectrons.at(1) = region1Electrons.at(1);
}
| unsigned short L1RCTElectronIsolationCard::getIsoElectrons | ( | int | i | ) | [inline] |
Definition at line 37 of file L1RCTElectronIsolationCard.h.
References isoElectrons.
{
return isoElectrons.at(i);
}
| unsigned short L1RCTElectronIsolationCard::getNonIsoElectrons | ( | int | i | ) | [inline] |
Definition at line 41 of file L1RCTElectronIsolationCard.h.
References nonIsoElectrons.
{
return nonIsoElectrons.at(i);
}
| void L1RCTElectronIsolationCard::print | ( | void | ) |
Definition at line 348 of file L1RCTElectronIsolationCard.cc.
References cardNo, gather_cfg::cout, isoElectrons, nonIsoElectrons, and regions.
{
std::cout << "Electron isolation card " << cardNo << std::endl;
std::cout << "Region 0 Information" << std::endl;
regions.at(0).print();
std::cout << "IsoElectron Candidate " << isoElectrons.at(0) << std::endl;
std::cout << "NonIsoElectron Candidate " << nonIsoElectrons.at(0) << std::endl << std::endl;
std::cout << "Region 1 Information" << std::endl;
regions.at(1).print();
std::cout << "IsoElectron Candidate " << isoElectrons.at(1) << std::endl;
std::cout << "NonIsoElectron Candidate " << nonIsoElectrons.at(1) << std::endl;
}
| void L1RCTElectronIsolationCard::printEdges | ( | ) | [inline] |
Definition at line 45 of file L1RCTElectronIsolationCard.h.
References regions.
| void L1RCTElectronIsolationCard::setRegion | ( | int | i, |
| const L1RCTRegion & | region | ||
| ) | [inline] |
Definition at line 32 of file L1RCTElectronIsolationCard.h.
References regions.
unsigned short L1RCTElectronIsolationCard::cardNo [private] |
Definition at line 57 of file L1RCTElectronIsolationCard.h.
Referenced by calcElectronCandidates(), cardNumber(), and print().
unsigned short L1RCTElectronIsolationCard::crtNo [private] |
Definition at line 56 of file L1RCTElectronIsolationCard.h.
Referenced by crateNumber().
L1RCTRegion L1RCTElectronIsolationCard::empty [private] |
Definition at line 61 of file L1RCTElectronIsolationCard.h.
std::vector<unsigned short> L1RCTElectronIsolationCard::isoElectrons [private] |
Definition at line 63 of file L1RCTElectronIsolationCard.h.
Referenced by fillElectronCandidates(), getIsoElectrons(), and print().
std::vector<unsigned short> L1RCTElectronIsolationCard::nonIsoElectrons [private] |
Definition at line 64 of file L1RCTElectronIsolationCard.h.
Referenced by fillElectronCandidates(), getNonIsoElectrons(), and print().
const L1RCTLookupTables* L1RCTElectronIsolationCard::rctLookupTables_ [private] |
Definition at line 59 of file L1RCTElectronIsolationCard.h.
Referenced by calcElectronCandidates().
std::vector<L1RCTRegion> L1RCTElectronIsolationCard::regions [private] |
Definition at line 65 of file L1RCTElectronIsolationCard.h.
Referenced by fillElectronCandidates(), L1RCTElectronIsolationCard(), print(), printEdges(), setRegion(), and ~L1RCTElectronIsolationCard().
1.7.3