|
|
|
#include <CSCSectorReceiverMiniLUT.h>
Static Public Member Functions | |
| static global_eta_data | calcGlobalEtaMEMini (unsigned short endcap, unsigned short sector, unsigned short station, unsigned short subsector, unsigned theadd) |
| static global_phi_data | calcGlobalPhiMBMini (unsigned short endcap, unsigned short sector, unsigned short subsector, unsigned theadd) |
| static global_phi_data | calcGlobalPhiMEMini (unsigned short endcap, unsigned short sector, unsigned short station, unsigned short subsector, unsigned theadd) |
| static lclphidat | calcLocalPhiMini (unsigned theadd) |
Static Private Attributes | |
| static const unsigned short int | gbl_eta_bounds [2][6][4][2][4][9][2] |
| static const float | gbl_eta_params [2][6][4][2][4][9][3] |
| static const float | gbl_phi_mb_params [2][6][2][9][2] |
| static const float | gbl_phi_me_params [2][6][4][2][9][2] |
| static const float | lcl_phi_param0 [1<< 4] |
| static const float | lcl_phi_param1 = 6.4 |
Provides a new way of defining the Lookup tables used by the core. Defines the lookup tables as parameterized functions in order to save on memory usage when compared to the current standard definitions of LUTs
Definition at line 15 of file CSCSectorReceiverMiniLUT.h.
| global_eta_data CSCSectorReceiverMiniLUT::calcGlobalEtaMEMini | ( | unsigned short | endcap, |
| unsigned short | sector, | ||
| unsigned short | station, | ||
| unsigned short | subsector, | ||
| unsigned | theadd | ||
| ) | [static] |
Definition at line 51 of file CSCSectorReceiverMiniLUT.cc.
References data, gbl_eta_bounds, gbl_eta_params, create_public_lumi_plots::log, and CSCConstants::MAX_NUM_WIRES.
Referenced by CSCSectorReceiverLUT::globalEtaME().
{
if(endcap < 1 || endcap > 2)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of endcap, " << endcap
<< ", is out of bounds, [1, 2] +++\n";
if(sector < 1 || sector > 6)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of sector, " << sector
<< ", is out of bounds, [1, 6] +++\n";
if(station < 1 || station > 4)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of station, " << station
<< ", is out of bounds, [1, 4] +++\n";
gbletadat data(0);
unsigned short int cscid = ((theadd >> 15) & 0xf);
unsigned short int lclPhi = ((theadd >> 6) & 0x3);
unsigned short int WG = ((theadd >> 8) & 0x7f);
unsigned short int bend = ((theadd) & 0x3f);
int eta_temp=999, eta_min=999, eta_max=999;
if((cscid > 0) && (cscid <= 9) && (WG < CSCConstants::MAX_NUM_WIRES))
{
if(station == 1)
{
eta_temp = (gbl_eta_params[endcap-1][sector-1][station-1][subsector-1][lclPhi][cscid-1][0] +
gbl_eta_params[endcap-1][sector-1][station-1][subsector-1][lclPhi][cscid-1][1] *
log(gbl_eta_params[endcap-1][sector-1][station-1][subsector-1][lclPhi][cscid-1][2] + WG));
eta_min = gbl_eta_bounds[endcap-1][sector-1][station-1][subsector-1][lclPhi][cscid-1][0];
eta_max = gbl_eta_bounds[endcap-1][sector-1][station-1][subsector-1][lclPhi][cscid-1][1];
// add offset to ME+11a, subtract for ME-11a (wire tilt and strip direction)1
if ( (cscid < 4)&&(lclPhi == 3) ) {
if (endcap == 1)
eta_temp += 3;
else
eta_temp -= 3;
}
}
else
{
eta_temp = (gbl_eta_params[endcap-1][sector-1][station-1][0][lclPhi][cscid-1][0] +
gbl_eta_params[endcap-1][sector-1][station-1][0][lclPhi][cscid-1][1] *
log(gbl_eta_params[endcap-1][sector-1][station-1][0][lclPhi][cscid-1][2] + WG));
eta_min = gbl_eta_bounds[endcap-1][sector-1][station-1][0][lclPhi][cscid-1][0];
eta_max = gbl_eta_bounds[endcap-1][sector-1][station-1][0][lclPhi][cscid-1][1];
}
}
else
{
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of cscid, " << cscid
<< ", is out of bounds, [1, 9] -- or --"
<< " Value of wire group, " << WG
<< ", exceeds max allowed, " << CSCConstants::MAX_NUM_WIRES << " +++\n";
}
// protect from negative numbers. If the value of eta_temp is <0, set global eta to the minimum value
if((eta_temp >= eta_min) &&
(eta_temp <= eta_max))
data.global_eta = eta_temp;
else if(eta_temp < eta_min)
data.global_eta = eta_min;
else
data.global_eta = eta_max;
// data.global_bend = 0;
// Just pass through lowest 5 bits of local bend (drop 1 MSB)
data.global_bend = bend & 0x1F;
return data;
}
| global_phi_data CSCSectorReceiverMiniLUT::calcGlobalPhiMBMini | ( | unsigned short | endcap, |
| unsigned short | sector, | ||
| unsigned short | subsector, | ||
| unsigned | theadd | ||
| ) | [static] |
Definition at line 212 of file CSCSectorReceiverMiniLUT.cc.
References data, gbl_phi_mb_params, CSCBitWidths::kGlobalPhiDataBitWidth, and CSCBitWidths::kLocalPhiDataBitWidth.
{
if(endcap < 1 || endcap > 2)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of endcap, " << endcap
<< ", is out of bounds, [1, 2] +++\n";
if(sector < 1 || sector > 6)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of sector, " << sector
<< ", is out of bounds, [1, 6] +++\n";
gblphidat data(0);
unsigned short int maxPhiL = 1<<CSCBitWidths::kLocalPhiDataBitWidth;
unsigned short int maxPhiG = 1<<CSCBitWidths::kGlobalPhiDataBitWidth;
unsigned short int cscid = ((theadd >> 15)&0xf);
unsigned short int lclPhi = (theadd & 0x3ff);
if((cscid <= 3) || (cscid >= 7))
maxPhiL = maxPhiL*(64./80); // currently a hack that is in place to handle the different number of strips in ME1/1 and ME1/3
if((cscid > 0) && (cscid <= 9))
{
if(lclPhi < maxPhiL)
data.global_phi = (gbl_phi_mb_params[endcap-1][sector-1][subsector-1][cscid-1][0] +
gbl_phi_mb_params[endcap-1][sector-1][subsector-1][cscid-1][1]*lclPhi);
else
data.global_phi = (gbl_phi_mb_params[endcap-1][sector-1][subsector-1][cscid-1][0] +
gbl_phi_mb_params[endcap-1][sector-1][subsector-1][cscid-1][1]*(maxPhiL-1));
}
else
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of cscid, " << cscid
<< ", is out of bounds, [1, 9] +++\n";
if(data.global_phi >= maxPhiG)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of global_phi, " << data.global_phi
<< ", exceeds max allowed, " << maxPhiG-1 << " +++\n";
if(data.global_phi >= maxPhiG)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of global_phi, " << data.global_phi
<< ", exceeds max allowed, " << maxPhiG-1 << " +++\n";
return data;
}
| global_phi_data CSCSectorReceiverMiniLUT::calcGlobalPhiMEMini | ( | unsigned short | endcap, |
| unsigned short | sector, | ||
| unsigned short | station, | ||
| unsigned short | subsector, | ||
| unsigned | theadd | ||
| ) | [static] |
Definition at line 131 of file CSCSectorReceiverMiniLUT.cc.
References data, gbl_phi_me_params, CSCBitWidths::kGlobalPhiDataBitWidth, CSCBitWidths::kLocalPhiDataBitWidth, CSCConstants::MAX_NUM_STRIPS, and strip().
Referenced by CSCSectorReceiverLUT::globalPhiME().
{
if(endcap < 1 || endcap > 2)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of endcap, " << endcap
<< ", is out of bounds, [1, 2] +++\n";
if(sector < 1 || sector > 6)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of sector, " << sector
<< ", is out of bounds, [1, 6] +++\n";
if(station < 1 || station > 4)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of station, " << station
<< ", is out of bounds, [1, 4] +++\n";
gblphidat data(0);
unsigned short int maxPhiL = 1<<CSCBitWidths::kLocalPhiDataBitWidth;
unsigned short int maxPhiG = 1<<CSCBitWidths::kGlobalPhiDataBitWidth;
unsigned short int cscid = ((theadd >> 15)&0xf);
unsigned short int lclPhi = (theadd & 0x3ff);
// 12/11/09
// GP et DA: how to identify the strip number and isolate and shift the localPhi value
const double binPhiL = static_cast<double>(maxPhiL)/(2*CSCConstants::MAX_NUM_STRIPS);
int strip = static_cast<int>(lclPhi/binPhiL);
// edm::LogWarning("GP Input") << " CSCID " << cscid << " strip:" << strip << " lclPhi: " << lclPhi << "theadd: " << theadd;
if (station == 1 && (cscid <= 3) && (strip >= 127 && strip < 160)){
//edm::LogWarning("GP Input") << " -> Matched Selection Criteria";
// in this case need to redefine lclPhi in order to
// place local phi in the middle of the 5th CFEB
// and not on the first third of the CFEB as default
lclPhi = (strip-127+31)*(4*binPhiL/3);
}
// end GP et DA
if(station == 1 && ((cscid <= 3) || (cscid >= 7))) {
//if ( (strip >= 127 && strip < 160) || (cscid >= 10) ) // VK: the || (cscid >= 10) for unganged ME1a
// maxPhiL = maxPhiL*(48./80); // GP et DA: currently a hack that is in place to handle the different number of strips in ME1/1a and ME1/3
//else
maxPhiL = maxPhiL*(64./80); // currently a hack that is in place to handle the different number of strips in ME1/1 and ME1/3
}
// VK: The the unganged ME1a hack
if (station == 1 && (cscid >= 10)) {
lclPhi = strip*(4*binPhiL/3);
cscid = cscid - 9; // back to normal 1-9 range
}
// end VK
if((cscid > 0) && (cscid <= 9))
{
if((station == 1) && (lclPhi < maxPhiL))
data.global_phi = (gbl_phi_me_params[endcap-1][sector-1][station-1][subsector-1][cscid-1][0] +
gbl_phi_me_params[endcap-1][sector-1][station-1][subsector-1][cscid-1][1]*lclPhi);
else if((station == 1) && (lclPhi >= maxPhiL))
data.global_phi = (gbl_phi_me_params[endcap-1][sector-1][station-1][subsector-1][cscid-1][0] +
gbl_phi_me_params[endcap-1][sector-1][station-1][subsector-1][cscid-1][1]*(maxPhiL-1));
else
data.global_phi = (gbl_phi_me_params[endcap-1][sector-1][station-1][0][cscid-1][0] +
gbl_phi_me_params[endcap-1][sector-1][station-1][0][cscid-1][1]*lclPhi);
}
else
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of cscid, " << cscid
<< ", is out of bounds, [1, 9] +++\n";
if(data.global_phi >= maxPhiG)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of global_phi, " << data.global_phi
<< ", exceeds max allowed, " << maxPhiG-1 << " +++\n";
return data;
}
| lclphidat CSCSectorReceiverMiniLUT::calcLocalPhiMini | ( | unsigned | theadd | ) | [static] |
Definition at line 21 of file CSCSectorReceiverMiniLUT.cc.
References data, CSCBitWidths::kLocalPhiDataBitWidth, lcl_phi_param0, lcl_phi_param1, CSCConstants::MAX_NUM_STRIPS, CSCConstants::NUM_CLCT_PATTERNS, listBenchmarks::pattern, and strip().
Referenced by CSCSectorReceiverLUT::localPhi().
{
// This method is ripped from CSCSectorReceverLUT.cc with minor changes
lclphidat data;
static int maxPhiL = 1<<CSCBitWidths::kLocalPhiDataBitWidth;
unsigned short int pattern = ((theadd >> 8) & 0xf);
unsigned short int strip = (theadd & 0xff);
if(strip < 2*CSCConstants::MAX_NUM_STRIPS && pattern < CSCConstants::NUM_CLCT_PATTERNS)
data.phi_local = static_cast<unsigned>((lcl_phi_param0[pattern] + strip)*lcl_phi_param1);
else
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of strip, " << strip
<< ", exceeds max allowed, " << 2*CSCConstants::MAX_NUM_STRIPS-1
<< " +++\n";
if(data.phi_local >= maxPhiL)
edm::LogWarning("CSCSectorReceiverMiniLUT")
<< "+++ Value of phi_local, " << data.phi_local
<< ", exceeds max allowed, " << CSCConstants::NUM_CLCT_PATTERNS-1 << " +++\n";
// data.phi_bend_local = 0;
// Just pass through all bits of pattern as bend angle (so 2 MSB unfilled)
data.phi_bend_local = pattern & 0x3F;
return data;
}
const unsigned short int CSCSectorReceiverMiniLUT::gbl_eta_bounds [static, private] |
Definition at line 28 of file CSCSectorReceiverMiniLUT.h.
Referenced by calcGlobalEtaMEMini().
const float CSCSectorReceiverMiniLUT::gbl_eta_params [static, private] |
Definition at line 27 of file CSCSectorReceiverMiniLUT.h.
Referenced by calcGlobalEtaMEMini().
const float CSCSectorReceiverMiniLUT::gbl_phi_mb_params [static, private] |
Definition at line 32 of file CSCSectorReceiverMiniLUT.h.
Referenced by calcGlobalPhiMBMini().
const float CSCSectorReceiverMiniLUT::gbl_phi_me_params [static, private] |
Definition at line 30 of file CSCSectorReceiverMiniLUT.h.
Referenced by calcGlobalPhiMEMini().
const float CSCSectorReceiverMiniLUT::lcl_phi_param0 [static, private] |
{
0.50, 0.50, -0.10, 1.10, -0.14, 1.14, 0.27, 0.73, 0.29, 0.71, 0.50, 0, 0, 0, 0, 0
}
Definition at line 24 of file CSCSectorReceiverMiniLUT.h.
Referenced by calcLocalPhiMini().
const float CSCSectorReceiverMiniLUT::lcl_phi_param1 = 6.4 [static, private] |
Definition at line 25 of file CSCSectorReceiverMiniLUT.h.
Referenced by calcLocalPhiMini().
1.7.3