|
|
|
#include <L1GtEtaPhiConversions.h>
Public Member Functions | |
| const bool | convertEtaIndex (const L1GtObject &, const unsigned int initialIndex, unsigned int &convertedIndex) const |
| void | convertL1Scales (const L1CaloGeometry *, const L1MuTriggerScales *, const int, const int) |
| perform all scale conversions | |
| const bool | convertPhiIndex (const unsigned int pairIndex, const unsigned int positionPair, const unsigned int initialIndex, unsigned int &convertedIndex) const |
| const unsigned int | gtObjectNrBinsPhi (const unsigned int) const |
| const unsigned int | gtObjectNrBinsPhi (const L1GtObject &) const |
| return the number of phi bins for a GT object | |
| const unsigned int | gtObjectNrBinsPhi (const L1GtObject &, const L1GtObject &) const |
| return the number of phi bins for a pair of GT objects, according to conversion rules | |
| const unsigned int | gtObjectPairIndex (const L1GtObject &, const L1GtObject &) const |
| L1GtEtaPhiConversions () | |
| constructor | |
| virtual void | print (std::ostream &myCout) const |
| print all the performed conversions | |
| void | setVerbosity (const int verbosity) |
| virtual | ~L1GtEtaPhiConversions () |
| destructor | |
Private Member Functions | |
| const double | rad2deg (const double &) const |
| convert phi from rad (-pi, pi] to deg (0, 360) | |
Private Attributes | |
| std::vector< std::pair < L1GtObject, L1GtObject > > | m_gtObjectPairVec |
| vector of all L1GtObject pairs | |
| bool | m_isDebugEnabled |
| cached edm::isDebugEnabled() | |
| const L1CaloGeometry * | m_l1CaloGeometry |
| pointer to calorimetry scales - updated in convertl1Scales method | |
| const L1MuTriggerScales * | m_l1MuTriggerScales |
| pointer to muon scales - updated in convertl1Scales method | |
| std::vector< unsigned int > | m_lutEtaCentralToCommonCalo |
| std::vector< unsigned int > | m_lutEtaForJetToCommonCalo |
| eta conversion of ForJet to the common calorimeter eta scale defined before | |
| std::vector< unsigned int > | m_lutEtaMuToCommonCalo |
| eta conversion of Mu to the common calorimeter eta scale defined before | |
| std::vector< unsigned int > | m_lutPhiEtmToHtm |
| phi conversion for ETM to HTM | |
| std::vector< unsigned int > | m_lutPhiEtmToJetEg |
| phi conversion for ETM to (*Jet, EG) | |
| std::vector< unsigned int > | m_lutPhiHtmToJetEg |
| phi conversion for HTM to (*Jet, EG) | |
| std::vector< unsigned int > | m_lutPhiJetEgToJetEg |
| std::vector< unsigned int > | m_lutPhiMuToEtm |
| phi conversion for Mu to ETM | |
| std::vector< unsigned int > | m_lutPhiMuToHtm |
| phi conversion for Mu to HTM | |
| std::vector< unsigned int > | m_lutPhiMuToJetEg |
| phi conversion for Mu to (*Jet, EG) | |
| unsigned int | m_nrBinsEtaCommon |
| number of eta bins for common scale | |
| unsigned int | m_nrBinsPhiEtm |
| number of phi bins for ETM | |
| unsigned int | m_nrBinsPhiHtm |
| number of phi bins for HTM | |
| unsigned int | m_nrBinsPhiJetEg |
| number of phi bins for calorimeter objects (*Jet, *EG) | |
| unsigned int | m_nrBinsPhiMu |
| number of phi bins for muons | |
| std::vector< bool > | m_pairConvertPhiFirstGtObject |
| std::vector< const unsigned int * > | m_pairNrPhiBinsVec |
| std::vector< const std::vector < unsigned int > * > | m_pairPhiConvVec |
| int | m_verbosity |
| verbosity level | |
Static Private Attributes | |
| static const unsigned int | badIndex = 999999 |
| static const double | PiConversion = 180. / acos(-1.) |
Description: convert eta and phi between various L1 trigger objects.
Implementation: <TODO: enter implementation details>
$Date$ $Revision$
Definition at line 38 of file L1GtEtaPhiConversions.h.
| L1GtEtaPhiConversions::L1GtEtaPhiConversions | ( | ) |
constructor
Definition at line 40 of file L1GtEtaPhiConversions.cc.
References CenJet, ETM, ForJet, HTM, IsoEG, l1GtObjectEnumToString(), LogTrace, m_gtObjectPairVec, m_isDebugEnabled, m_lutPhiEtmToHtm, m_lutPhiEtmToJetEg, m_lutPhiHtmToJetEg, m_lutPhiJetEgToJetEg, m_lutPhiMuToEtm, m_lutPhiMuToHtm, m_lutPhiMuToJetEg, m_nrBinsPhiEtm, m_nrBinsPhiHtm, m_nrBinsPhiJetEg, m_pairConvertPhiFirstGtObject, m_pairNrPhiBinsVec, m_pairPhiConvVec, Mu, NoIsoEG, and TauJet.
:
m_nrBinsPhiMu(0), m_nrBinsPhiJetEg(0), m_nrBinsPhiEtm(0),
m_nrBinsPhiHtm(0), m_nrBinsEtaCommon(0), m_verbosity(0),
m_isDebugEnabled(edm::isDebugEnabled()) {
// prepare the pairs of L1GtObjects, reserve (by hand) memory for vectors
// the index of the pair in m_gtObjectPairVec is used to extract
// the information from the other vectors, so the push_back must be done
// coherently
std::pair < L1GtObject, L1GtObject > gtObjPair;
m_gtObjectPairVec.reserve(56);
m_pairConvertPhiFirstGtObject.reserve(56);
m_pairNrPhiBinsVec.reserve(56);
m_pairPhiConvVec.reserve(56);
// Mu -> Jet & EG & ETM & HTM
//
gtObjPair = std::make_pair(Mu, CenJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
gtObjPair = std::make_pair(CenJet, Mu);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
//
gtObjPair = std::make_pair(Mu, ForJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
gtObjPair = std::make_pair(ForJet, Mu);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
//
gtObjPair = std::make_pair(Mu, TauJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
gtObjPair = std::make_pair(TauJet, Mu);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
//
gtObjPair = std::make_pair(Mu, NoIsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
gtObjPair = std::make_pair(NoIsoEG, Mu);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
//
gtObjPair = std::make_pair(Mu, IsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
gtObjPair = std::make_pair(IsoEG, Mu);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiMuToJetEg);
//
gtObjPair = std::make_pair(Mu, ETM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiEtm);
m_pairPhiConvVec.push_back(&m_lutPhiMuToEtm);
gtObjPair = std::make_pair(ETM, Mu);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiEtm);
m_pairPhiConvVec.push_back(&m_lutPhiMuToEtm);
//
gtObjPair = std::make_pair(Mu, HTM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiHtm);
m_pairPhiConvVec.push_back(&m_lutPhiMuToHtm);
gtObjPair = std::make_pair(HTM, Mu);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiHtm);
m_pairPhiConvVec.push_back(&m_lutPhiMuToHtm);
// ETM -> Jet & EG
//
gtObjPair = std::make_pair(ETM, CenJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
gtObjPair = std::make_pair(CenJet, ETM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
//
gtObjPair = std::make_pair(ETM, ForJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
gtObjPair = std::make_pair(ForJet, ETM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
//
gtObjPair = std::make_pair(ETM, TauJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
gtObjPair = std::make_pair(TauJet, ETM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
//
gtObjPair = std::make_pair(ETM, NoIsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
gtObjPair = std::make_pair(NoIsoEG, ETM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
//
gtObjPair = std::make_pair(ETM, IsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
gtObjPair = std::make_pair(IsoEG, ETM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToJetEg);
// HTM -> Jet & EG
//
gtObjPair = std::make_pair(HTM, CenJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
gtObjPair = std::make_pair(CenJet, HTM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
//
gtObjPair = std::make_pair(HTM, ForJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
gtObjPair = std::make_pair(ForJet, HTM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
//
gtObjPair = std::make_pair(HTM, TauJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
gtObjPair = std::make_pair(TauJet, HTM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
//
gtObjPair = std::make_pair(HTM, NoIsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
gtObjPair = std::make_pair(NoIsoEG, HTM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
//
gtObjPair = std::make_pair(HTM, IsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
gtObjPair = std::make_pair(IsoEG, HTM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiHtmToJetEg);
// ETM -> HTM
//
gtObjPair = std::make_pair(ETM, HTM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiHtm);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToHtm);
gtObjPair = std::make_pair(HTM, ETM);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiHtm);
m_pairPhiConvVec.push_back(&m_lutPhiEtmToHtm);
// Jet & EG -> Jet & EG
//
gtObjPair = std::make_pair(CenJet, ForJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(ForJet, CenJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(CenJet, TauJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(TauJet, CenJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(CenJet, NoIsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(NoIsoEG, CenJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(CenJet, IsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(IsoEG, CenJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(ForJet, TauJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(TauJet, ForJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(ForJet, NoIsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(NoIsoEG, ForJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(ForJet, IsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(IsoEG, ForJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(TauJet, NoIsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(NoIsoEG, TauJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(TauJet, IsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(IsoEG, TauJet);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
//
gtObjPair = std::make_pair(NoIsoEG, IsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(true);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
gtObjPair = std::make_pair(IsoEG, NoIsoEG);
m_gtObjectPairVec.push_back(gtObjPair);
m_pairConvertPhiFirstGtObject.push_back(false);
m_pairNrPhiBinsVec.push_back(&m_nrBinsPhiJetEg);
m_pairPhiConvVec.push_back(&m_lutPhiJetEgToJetEg);
// m_verbosity can not be used here, as L1GtEtaPhiConversions is called
// in L1GlobalTriggerGTL constructor, where m_verbosity is not yet set
if (m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << "\nm_gtObjectPairVec size: "
<< (m_gtObjectPairVec.size()) << std::endl;
unsigned int iPair = 0;
for (std::vector<std::pair<L1GtObject, L1GtObject> >::const_iterator
cIter = m_gtObjectPairVec.begin();
cIter != m_gtObjectPairVec.end(); ++cIter) {
LogTrace("L1GlobalTrigger") << "m_gtObjectPairVec vector element ["
<< l1GtObjectEnumToString((*cIter).first) << ", "
<< l1GtObjectEnumToString((*cIter).second)
<< "], \t\tpair index = " << iPair << std::endl;
iPair++;
}
}
}
| L1GtEtaPhiConversions::~L1GtEtaPhiConversions | ( | ) | [virtual] |
| const bool L1GtEtaPhiConversions::convertEtaIndex | ( | const L1GtObject & | gtObject, |
| const unsigned int | initialIndex, | ||
| unsigned int & | convertedIndex | ||
| ) | const |
convert the eta index initialIndex for a L1GtObject object to common scale converted index returned by reference method return true, if initial index within scale size otherwise (hardware error), return false
Definition at line 640 of file L1GtEtaPhiConversions.cc.
References badIndex, BPTX, CenJet, ETM, ETT, ForJet, GtExternal, HfBitCounts, HfRingEtSums, HTM, HTT, IsoEG, JetCounts, l1GtObjectEnumToString(), LogTrace, m_isDebugEnabled, m_lutEtaCentralToCommonCalo, m_lutEtaForJetToCommonCalo, m_lutEtaMuToCommonCalo, m_verbosity, Mu, NoIsoEG, ObjNull, TauJet, and TechTrig.
Referenced by L1GtCorrelationCondition::evaluateCondition().
{
unsigned int newIndex = badIndex;
bool conversionStatus = false;
switch (gtObject) {
case Mu: {
// check if initial index is within the scale size
// could be outside the scale size if there are hardware errors
// or wrong scale conversions
if (initialIndex >= m_lutEtaMuToCommonCalo.size()) {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " has initial eta index " << initialIndex
<< " >= " << (m_lutEtaMuToCommonCalo.size())
<< " scale size. Conversion failed." << std::endl;
}
} else {
// convert the index
newIndex = m_lutEtaMuToCommonCalo[initialIndex];
if (newIndex != badIndex) {
conversionStatus = true;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " initial eta index " << initialIndex
<< " (within scale size "
<< (m_lutEtaMuToCommonCalo.size())
<< ") converted to " << newIndex << std::endl;
}
} else {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " initial eta index " << initialIndex
<< " (within scale size "
<< (m_lutEtaMuToCommonCalo.size())
<< ") converted to badIndex" << newIndex
<< " Conversion failed." << std::endl;
}
}
}
}
break;
case NoIsoEG:
case IsoEG:
case CenJet:
case TauJet: {
// check if initial index is within the scale size
// could be outside the scale size if there are hardware errors
// or wrong scale conversions
if (initialIndex >= m_lutEtaCentralToCommonCalo.size()) {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " has initial eta index " << initialIndex
<< " >= " << (m_lutEtaCentralToCommonCalo.size())
<< " scale size. Conversion failed." << std::endl;
}
} else {
// convert the index
newIndex = m_lutEtaCentralToCommonCalo[initialIndex];
if (newIndex != badIndex) {
conversionStatus = true;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " initial eta index " << initialIndex
<< " (within scale size "
<< (m_lutEtaMuToCommonCalo.size())
<< ") converted to " << newIndex << std::endl;
}
} else {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " initial eta index " << initialIndex
<< " (within scale size "
<< (m_lutEtaCentralToCommonCalo.size())
<< ") converted to badIndex" << newIndex
<< " Conversion failed." << std::endl;
}
}
}
}
break;
case ForJet: {
// check if initial index is within the scale size
// could be outside the scale size if there are hardware errors
// or wrong scale conversions
if (initialIndex >= m_lutEtaForJetToCommonCalo.size()) {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " has initial eta index " << initialIndex
<< " >= " << (m_lutEtaForJetToCommonCalo.size())
<< " scale size. Conversion failed." << std::endl;
}
} else {
// convert the index
newIndex = m_lutEtaForJetToCommonCalo[initialIndex];
if (newIndex != badIndex) {
conversionStatus = true;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " initial eta index " << initialIndex
<< " (within scale size "
<< (m_lutEtaMuToCommonCalo.size())
<< ") converted to " << newIndex << std::endl;
}
} else {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << " L1 GT object "
<< (l1GtObjectEnumToString(gtObject))
<< " initial eta index " << initialIndex
<< " (within scale size "
<< (m_lutEtaForJetToCommonCalo.size())
<< ") converted to badIndex" << newIndex
<< " Conversion failed." << std::endl;
}
}
}
}
break;
case ETM:
case ETT:
case HTT:
case HTM:
case JetCounts:
case HfBitCounts:
case HfRingEtSums:
case TechTrig:
case Castor:
case BPTX:
case GtExternal:
case ObjNull: {
//no conversions needed, there is no eta quantity for these objects
conversionStatus = false;
}
break;
default: {
edm::LogInfo("L1GtObject") << "\n '" << (l1GtObjectEnumToString(
gtObject)) << "' is not a recognized L1GtObject. "
<< "\n Conversion failed. " << std::endl;
conversionStatus = false;
}
break;
}
//
convertedIndex = newIndex;
return conversionStatus;
}
| void L1GtEtaPhiConversions::convertL1Scales | ( | const L1CaloGeometry * | l1CaloGeometry, |
| const L1MuTriggerScales * | l1MuTriggerScales, | ||
| const int | ifCaloEtaNumberBits, | ||
| const int | ifMuEtaNumberBits | ||
| ) |
perform all scale conversions
Definition at line 950 of file L1GtEtaPhiConversions.cc.
References badIndex, L1CaloGeometry::emJetPhiBinHighEdge(), L1CaloGeometry::emJetPhiBinLowEdge(), L1CaloGeometry::etaBinCenter(), L1CaloGeometry::etSumPhiBinHighEdge(), L1CaloGeometry::etSumPhiBinLowEdge(), L1MuTriggerScales::getGMTEtaScale(), L1MuScale::getHighEdge(), L1MuScale::getLowEdge(), L1MuScale::getNBins(), L1MuTriggerScales::getPhiScale(), L1MuScale::getValue(), L1CaloGeometry::globalEtaBinLowEdge(), L1CaloGeometry::globalEtaIndex(), L1CaloGeometry::htSumPhiBinHighEdge(), L1CaloGeometry::htSumPhiBinLowEdge(), l1CaloGeometry_cfi::l1CaloGeometry, LogDebug, LogTrace, m_isDebugEnabled, m_l1CaloGeometry, m_l1MuTriggerScales, m_lutEtaCentralToCommonCalo, m_lutEtaForJetToCommonCalo, m_lutEtaMuToCommonCalo, m_lutPhiEtmToHtm, m_lutPhiEtmToJetEg, m_lutPhiHtmToJetEg, m_lutPhiJetEgToJetEg, m_lutPhiMuToEtm, m_lutPhiMuToHtm, m_lutPhiMuToJetEg, m_nrBinsEtaCommon, m_nrBinsPhiEtm, m_nrBinsPhiHtm, m_nrBinsPhiJetEg, m_nrBinsPhiMu, m_verbosity, L1CaloGeometry::numberGctCentralEtaBinsPerHalf(), L1CaloGeometry::numberGctEmJetPhiBins(), L1CaloGeometry::numberGctEtSumPhiBins(), L1CaloGeometry::numberGctForwardEtaBinsPerHalf(), and L1CaloGeometry::numberGctHtSumPhiBins().
Referenced by L1GlobalTriggerGTL::run().
{
// no bullet-proof method, depends on binning ...
// decide "by hand / by eyes" which object converts to which object
// update the scales used
m_l1CaloGeometry = l1CaloGeometry;
m_l1MuTriggerScales = l1MuTriggerScales;
// number of bins for all phi scales used
m_nrBinsPhiMu = 144; // FIXME ask Ivan for size() ...
//m_nrBinsPhiMu = m_l1MuTriggerScales->getPhiScale()->size();
m_nrBinsPhiJetEg = m_l1CaloGeometry->numberGctEmJetPhiBins();
m_nrBinsPhiEtm = m_l1CaloGeometry->numberGctEtSumPhiBins();
m_nrBinsPhiHtm = m_l1CaloGeometry->numberGctHtSumPhiBins();
//
// convert phi scale for muon (finer) to phi scale for (*Jet, EG) / ETM / HTM (coarser)
//
m_lutPhiMuToJetEg.clear();
m_lutPhiMuToJetEg.resize(m_nrBinsPhiMu, badIndex);
m_lutPhiMuToEtm.clear();
m_lutPhiMuToEtm.resize(m_nrBinsPhiMu, badIndex);
m_lutPhiMuToHtm.clear();
m_lutPhiMuToHtm.resize(m_nrBinsPhiMu, badIndex);
for (unsigned int phiMuInd = 0; phiMuInd < m_nrBinsPhiMu; ++phiMuInd) {
double phiMuLowEdge = m_l1MuTriggerScales->getPhiScale()->getLowEdge(
phiMuInd);
double phiMuHighEdge = m_l1MuTriggerScales->getPhiScale()->getHighEdge(
phiMuInd);
// to avoid precision problems, add a small quantity to phiMuLowEdge
double phiMuLowEdgeSmallShiftRight = phiMuLowEdge + (phiMuHighEdge
- phiMuLowEdge) / 100.;
// phi Mu -> (*Jet, EG)
unsigned int nrBins = m_nrBinsPhiJetEg;
for (unsigned int iBin = nrBins;; --iBin) {
double phiLowEdge = m_l1CaloGeometry->emJetPhiBinLowEdge(iBin);
double phiHighEdge = m_l1CaloGeometry->emJetPhiBinHighEdge(iBin);
if (phiMuLowEdgeSmallShiftRight >= phiLowEdge) {
m_lutPhiMuToJetEg[phiMuInd] = iBin % nrBins;
LogTrace("L1GlobalTrigger") << " phiMuIndex \t" << phiMuInd
<< " [ " << phiMuLowEdge << " \t, "
<< phiMuHighEdge << "] \t==>\t phiMuJetEG \t"
<< m_lutPhiMuToJetEg[phiMuInd] << " [ "
<< phiLowEdge << "\t, " << phiHighEdge << " ]"
<< std::endl;
break;
}
}
// phi Mu -> ETM
nrBins = m_nrBinsPhiEtm;
for (unsigned int iBin = nrBins;; --iBin) {
double phiLowEdge = m_l1CaloGeometry->etSumPhiBinLowEdge(iBin);
double phiHighEdge = m_l1CaloGeometry->etSumPhiBinHighEdge(iBin);
if (phiMuLowEdgeSmallShiftRight >= phiLowEdge) {
m_lutPhiMuToEtm[phiMuInd] = iBin % nrBins;
LogTrace("L1GlobalTrigger") << " phiMuIndex \t" << phiMuInd
<< " [ " << phiMuLowEdge << " \t, "
<< phiMuHighEdge << "] \t==>\t phiMuToEtm \t"
<< m_lutPhiMuToEtm[phiMuInd] << " [ " << phiLowEdge
<< "\t, " << phiHighEdge << " ]" << std::endl;
break;
}
}
// phi Mu -> HTM
nrBins = m_nrBinsPhiHtm;
for (unsigned int iBin = nrBins;; --iBin) {
double phiLowEdge = m_l1CaloGeometry->htSumPhiBinLowEdge(iBin);
double phiHighEdge = m_l1CaloGeometry->htSumPhiBinHighEdge(iBin);
if (phiMuLowEdgeSmallShiftRight >= phiLowEdge) {
m_lutPhiMuToHtm[phiMuInd] = iBin % nrBins;
LogTrace("L1GlobalTrigger") << " phiMuIndex \t" << phiMuInd
<< " [ " << phiMuLowEdge << " \t, "
<< phiMuHighEdge << "] \t==>\t phiMuToHtm \t"
<< m_lutPhiMuToHtm[phiMuInd] << " [ " << phiLowEdge
<< "\t, " << phiHighEdge << " ]" << std::endl;
break;
}
}
LogTrace("L1GlobalTrigger") << std::endl;
}
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << "Mu phi conversions" << std::endl;
for (unsigned int iBin = 0; iBin < m_nrBinsPhiMu; ++iBin) {
LogTrace("L1GlobalTrigger") << " Mu phiIndex \t" << iBin
<< "\t converted to index:"
<< "\t Jet-EG \t" << m_lutPhiMuToJetEg.at(iBin)
<< "\t ETM \t" << m_lutPhiMuToEtm.at(iBin)
<< "\t HTM \t" << m_lutPhiMuToHtm.at(iBin)
<< std::endl;
}
LogTrace("L1GlobalTrigger") << std::endl;
}
//
// convert phi scale for ETM to phi scale for (*Jet, EG) / HTM (coarser)
//
m_lutPhiEtmToJetEg.clear();
m_lutPhiEtmToJetEg.resize(m_nrBinsPhiEtm, badIndex);
m_lutPhiEtmToHtm.clear();
m_lutPhiEtmToHtm.resize(m_nrBinsPhiEtm, badIndex);
for (unsigned int phiEtmInd = 0; phiEtmInd < m_nrBinsPhiEtm; ++phiEtmInd) {
double phiEtmLowEdge = m_l1CaloGeometry->etSumPhiBinLowEdge(phiEtmInd);
double phiEtmHighEdge = m_l1CaloGeometry->etSumPhiBinHighEdge(phiEtmInd);
// to avoid precision problems, add a small quantity to phiEtmLowEdge
double phiEtmLowEdgeSmallShiftRight = phiEtmLowEdge + (phiEtmHighEdge
- phiEtmLowEdge) / 100.;
// phi ETM -> (*Jet, EG)
unsigned int nrBins = m_nrBinsPhiJetEg;
for (unsigned int iBin = nrBins;; --iBin) {
double phiLowEdge = m_l1CaloGeometry->emJetPhiBinLowEdge(iBin);
double phiHighEdge = m_l1CaloGeometry->emJetPhiBinHighEdge(iBin);
if (phiEtmLowEdgeSmallShiftRight >= phiLowEdge) {
m_lutPhiEtmToJetEg[phiEtmInd] = iBin % nrBins;
LogTrace("L1GlobalTrigger") << " phiEtmIndex \t" << phiEtmInd
<< " [ " << phiEtmLowEdge << " \t, "
<< phiEtmHighEdge << "] \t==>\t phiEtmJetEG \t"
<< m_lutPhiEtmToJetEg[phiEtmInd] << " [ "
<< phiLowEdge << "\t, " << phiHighEdge << " ]"
<< std::endl;
break;
}
}
// phi ETM -> HTM
nrBins = m_nrBinsPhiHtm;
for (unsigned int iBin = nrBins;; --iBin) {
double phiLowEdge = m_l1CaloGeometry->htSumPhiBinLowEdge(iBin);
double phiHighEdge = m_l1CaloGeometry->htSumPhiBinHighEdge(iBin);
if (phiEtmLowEdgeSmallShiftRight >= phiLowEdge) {
m_lutPhiEtmToHtm[phiEtmInd] = iBin % nrBins;
LogTrace("L1GlobalTrigger") << " phiEtmIndex \t" << phiEtmInd
<< " [ " << phiEtmLowEdge << " \t, "
<< phiEtmHighEdge << "] \t==>\t phiEtmToHtm \t"
<< m_lutPhiEtmToHtm[phiEtmInd] << " [ " << phiLowEdge
<< "\t, " << phiHighEdge << " ]" << std::endl;
break;
}
}
LogTrace("L1GlobalTrigger") << std::endl;
}
//
// convert phi scale for HTM to phi scale for (*Jet, EG)
//
m_lutPhiHtmToJetEg.clear();
m_lutPhiHtmToJetEg.resize(m_nrBinsPhiHtm, badIndex);
for (unsigned int phiHtmInd = 0; phiHtmInd < m_nrBinsPhiHtm; ++phiHtmInd) {
double phiHtmLowEdge = m_l1CaloGeometry->htSumPhiBinLowEdge(phiHtmInd);
double phiHtmHighEdge = m_l1CaloGeometry->htSumPhiBinHighEdge(phiHtmInd);
// to avoid precision problems, add a small quantity to phiHtmLowEdge
double phiHtmLowEdgeSmallShiftRight = phiHtmLowEdge + (phiHtmHighEdge
- phiHtmLowEdge) / 100.;
unsigned int nrBins = m_nrBinsPhiJetEg;
for (unsigned int iBin = nrBins;; --iBin) {
double phiLowEdge = m_l1CaloGeometry->emJetPhiBinLowEdge(iBin);
double phiHighEdge = m_l1CaloGeometry->emJetPhiBinHighEdge(iBin);
if (phiHtmLowEdgeSmallShiftRight >= phiLowEdge) {
m_lutPhiHtmToJetEg[phiHtmInd] = iBin % nrBins;
LogTrace("L1GlobalTrigger") << " phiHtmIndex \t" << phiHtmInd
<< " [ " << phiHtmLowEdge << " \t, "
<< phiHtmHighEdge << "] \t==>\t phiHtmJetEG \t"
<< m_lutPhiHtmToJetEg[phiHtmInd] << " [ "
<< phiLowEdge << "\t, " << phiHighEdge << " ]"
<< std::endl;
break;
}
}
}
//
// convert phi scale for (*Jet, EG) to (*Jet, EG)
// dummy - return the same index as the input index
m_lutPhiJetEgToJetEg.clear();
m_lutPhiJetEgToJetEg.resize(m_nrBinsPhiJetEg, badIndex);
for (unsigned int phiInd = 0; phiInd < m_nrBinsPhiJetEg; ++phiInd) {
m_lutPhiJetEgToJetEg[phiInd] = phiInd;
}
//
// eta conversions
//
// all objects are converted to a common central / forward calorimeter eta scale,
// built by setting together the forward scale and the central scale
//
// eta is signed, MSB is the sign for all objects - must be taken into account
// in conversion - the common scale is, from 0 to m_nrBinsEtaCommon:
//
// [ForJet negative bins][Central Jet/IsoEG/NoIsoEG negative bins][Central Jet/IsoEG/NoIsoEG positive bins][ForJet positive bins]
//
unsigned int nrGctCentralEtaBinsPerHalf =
m_l1CaloGeometry->numberGctCentralEtaBinsPerHalf();
unsigned int nrGctForwardEtaBinsPerHalf =
m_l1CaloGeometry->numberGctForwardEtaBinsPerHalf();
unsigned int nrGctTotalEtaBinsPerHalf = nrGctCentralEtaBinsPerHalf
+ nrGctForwardEtaBinsPerHalf;
m_nrBinsEtaCommon = 2*nrGctTotalEtaBinsPerHalf;
//
// convert eta scale for CenJet/TauJet & IsoEG/NoIsoEG to a common
// central / forward calorimeter eta scale
//
// get the sign and the index absolute value
LogTrace("L1GlobalTrigger")
<< " \nEta conversion: CenJet/TauJet & IsoEG/NoIsoEG to a common calorimeter scale\n"
<< std::endl;
m_lutEtaCentralToCommonCalo.clear();
m_lutEtaCentralToCommonCalo.resize(
(nrGctCentralEtaBinsPerHalf | (1 << (ifCaloEtaNumberBits - 1))),
badIndex);
for (unsigned int etaInd = 0; etaInd < nrGctCentralEtaBinsPerHalf; ++etaInd) {
// for positive values, the index is etaInd
unsigned int globalIndex = m_l1CaloGeometry->globalEtaIndex(
m_l1CaloGeometry->etaBinCenter(etaInd, true));
m_lutEtaCentralToCommonCalo[etaInd] = globalIndex;
LogTrace("L1GlobalTrigger") << " etaIndex " << etaInd << "\t [hex: "
<< std::hex << etaInd << "] " << std::dec
<< " ==> etaIndexGlobal " << globalIndex << std::endl;
// for negative values, one adds (binary) 1 as MSB to the index
unsigned int etaIndNeg = etaInd | (1 << (ifCaloEtaNumberBits - 1));
globalIndex = m_l1CaloGeometry->globalEtaIndex(
m_l1CaloGeometry->etaBinCenter(etaIndNeg, true));
m_lutEtaCentralToCommonCalo[etaIndNeg] = globalIndex;
LogTrace("L1GlobalTrigger") << " etaIndex " << etaIndNeg
<< "\t [hex: " << std::hex << etaIndNeg << "] " << std::dec
<< " ==> etaIndexGlobal " << globalIndex << std::endl;
}
//
// convert eta scale for ForJet to a common
// central / forward calorimeter eta scale
//
LogTrace("L1GlobalTrigger")
<< " \nEta conversion: ForJet to a common calorimeter scale\n"
<< std::endl;
m_lutEtaForJetToCommonCalo.clear();
m_lutEtaForJetToCommonCalo.resize(
(nrGctForwardEtaBinsPerHalf | (1 << (ifCaloEtaNumberBits - 1))),
badIndex);
for (unsigned int etaInd = 0; etaInd < nrGctForwardEtaBinsPerHalf; ++etaInd) {
// for positive values, the index is etaInd
unsigned int globalIndex = m_l1CaloGeometry->globalEtaIndex(
m_l1CaloGeometry->etaBinCenter(etaInd, false));
m_lutEtaForJetToCommonCalo[etaInd] = globalIndex;
LogTrace("L1GlobalTrigger") << " etaIndex " << etaInd << "\t [hex: "
<< std::hex << etaInd << "] " << std::dec
<< " ==> etaIndexGlobal " << globalIndex << std::endl;
// for negative values, one adds (binary) 1 as MSB to the index
unsigned int etaIndNeg = etaInd | (1 << (ifCaloEtaNumberBits - 1));
globalIndex = m_l1CaloGeometry->globalEtaIndex(
m_l1CaloGeometry->etaBinCenter(etaIndNeg, false));
m_lutEtaForJetToCommonCalo[etaIndNeg] = globalIndex;
LogTrace("L1GlobalTrigger") << " etaIndex " << etaIndNeg << "\t [hex: "
<< std::hex << etaIndNeg << "] " << std::dec
<< " ==> etaIndexGlobal " << globalIndex << std::endl;
}
//
// convert eta scale for Mu to a common
// central / forward calorimeter eta scale
//
LogDebug("L1GlobalTrigger")
<< " \nEta conversion: Mu to a common calorimeter scale\n"
<< std::endl;
// eta scale defined for positive values - need to be symmetrized
unsigned int nrBinsEtaMuPerHalf =
m_l1MuTriggerScales->getGMTEtaScale()->getNBins();
LogTrace("L1GlobalTrigger") << " \nnrBinsEtaMuPerHalf = "
<< nrBinsEtaMuPerHalf << "\n" << std::endl;
m_lutEtaMuToCommonCalo.clear();
m_lutEtaMuToCommonCalo.resize(
(nrBinsEtaMuPerHalf | (1 << (ifMuEtaNumberBits - 1))), badIndex);
for (unsigned int etaMuInd = 0; etaMuInd < nrBinsEtaMuPerHalf; ++etaMuInd) {
double etaMuLowEdge = m_l1MuTriggerScales->getGMTEtaScale()->getValue(
etaMuInd);
double etaMuHighEdge = m_l1MuTriggerScales->getGMTEtaScale()->getValue(
etaMuInd + 1);
// to avoid precision problems, add a small quantity to etaMuLowEdge
double etaMuLowEdgeSmallShiftRight = etaMuLowEdge + (etaMuHighEdge
- etaMuLowEdge) / 100.;
// positive values
for (unsigned int iBin = m_nrBinsEtaCommon;; --iBin) {
double etaLowEdge = m_l1CaloGeometry->globalEtaBinLowEdge(iBin);
double etaHighEdge = 0.0;
if (iBin == m_nrBinsEtaCommon) {
etaHighEdge = etaLowEdge;
} else {
etaHighEdge = m_l1CaloGeometry->globalEtaBinLowEdge(iBin + 1);
}
if (etaMuLowEdgeSmallShiftRight >= etaLowEdge) {
m_lutEtaMuToCommonCalo[etaMuInd] = iBin % m_nrBinsEtaCommon;
LogTrace("L1GlobalTrigger") << " etaMuIndex \t" << etaMuInd
<< "\t [ " << etaMuLowEdge << ", \t" << etaMuHighEdge
<< "] ==> etaMuJetEG \t"
<< m_lutEtaMuToCommonCalo[etaMuInd] << "\t [ "
<< etaLowEdge << ", \t" << etaHighEdge << " ]"
<< std::endl;
break;
}
}
// for negative values, one adds (binary) 1 as MSB to the index
unsigned int etaMuIndNeg = etaMuInd | (1 << (ifMuEtaNumberBits - 1));
m_lutEtaMuToCommonCalo[etaMuIndNeg] = m_lutEtaMuToCommonCalo[0]
- (m_lutEtaMuToCommonCalo[etaMuInd] - m_lutEtaMuToCommonCalo[0]
+ 1);
LogTrace("L1GlobalTrigger") << " etaMuIndexNeg \t" << etaMuIndNeg
<< "\t [ " << (-1.0 * etaMuLowEdge) << ", \t" << (-1.0
* etaMuHighEdge) << "] ==> etaMuJetEG \t"
<< m_lutEtaMuToCommonCalo[etaMuIndNeg] << "\t [ "
<< m_l1CaloGeometry->globalEtaBinLowEdge(
m_lutEtaMuToCommonCalo[etaMuIndNeg]) << ", \t"
<< m_l1CaloGeometry->globalEtaBinLowEdge(
m_lutEtaMuToCommonCalo[etaMuIndNeg] + 1) << " ]"
<< std::endl;
}
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << std::endl;
LogTrace("L1GlobalTrigger") << std::endl;
}
}
| const bool L1GtEtaPhiConversions::convertPhiIndex | ( | const unsigned int | pairIndex, |
| const unsigned int | positionPair, | ||
| const unsigned int | initialIndex, | ||
| unsigned int & | convertedIndex | ||
| ) | const |
convert the phi index initialIndex for an object from pair pairIndex, with position of object in pair positionPair to common scale for the L1GtObject pair converted index returned by reference method return true, if initial index within scale size otherwise (hardware error), return false
Definition at line 493 of file L1GtEtaPhiConversions.cc.
References badIndex, Exception, LogTrace, m_isDebugEnabled, m_pairConvertPhiFirstGtObject, m_pairPhiConvVec, and m_verbosity.
Referenced by L1GtCorrelationCondition::evaluateCondition().
{
unsigned int newIndex = badIndex;
bool conversionStatus = false;
// check if initial index is within the scale size
// could be outside the scale size if there are hardware errors
// or wrong scale conversions
if (initialIndex >= (*(m_pairPhiConvVec.at(pairIndex))).size()) {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << (positionPair ? " Second"
: "\n First") << " object from pair " << pairIndex
<< ": initial phi index " << initialIndex << " >= "
<< ((*(m_pairPhiConvVec.at(pairIndex))).size())
<< " Conversion failed." << std::endl;
}
} else {
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << (positionPair ? " Second"
: "\n First") << " object from pair " << pairIndex
<< ": initial phi index " << initialIndex
<< " within scale size " <<
((*(m_pairPhiConvVec.at(pairIndex))).size())
<< std::endl;
}
}
// convert the index
switch (positionPair) {
case 0: {
if (m_pairConvertPhiFirstGtObject.at(pairIndex)) {
newIndex = (*(m_pairPhiConvVec.at(pairIndex))).at(initialIndex);
if (newIndex != badIndex) {
conversionStatus = true;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger")
<< (positionPair ? " Second" : "\n First")
<< " object from pair " << pairIndex
<< ": initial phi index " << initialIndex
<< " converted to " << newIndex << std::endl;
}
} else {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger")
<< (positionPair ? " Second" : "\n First")
<< " object from pair " << pairIndex
<< ": converted phi index " << newIndex
<< "is equal to badIndex " << badIndex
<< " Conversion failed." << std::endl;
}
}
} else {
newIndex = initialIndex;
conversionStatus = true;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << (positionPair ? " Second"
: "\n First") << " object from pair " << pairIndex
<< ": initial phi index " << initialIndex
<< " not requested to be converted, return index "
<< newIndex << std::endl;
}
}
}
break;
case 1: {
if (m_pairConvertPhiFirstGtObject.at(pairIndex)) {
newIndex = initialIndex;
conversionStatus = true;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger") << (positionPair ? " Second"
: "\n First") << " object from pair " << pairIndex
<< ": initial phi index " << initialIndex
<< " not requested to be converted, return index, return index "
<< newIndex << std::endl;
}
} else {
newIndex = (*(m_pairPhiConvVec.at(pairIndex))).at(initialIndex);
if (newIndex != badIndex) {
conversionStatus = true;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger")
<< (positionPair ? " Second" : "\n First")
<< " object from pair " << pairIndex
<< ": initial phi index " << initialIndex
<< " converted to " << newIndex << std::endl;
}
} else {
conversionStatus = false;
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger")
<< (positionPair ? " Second" : "\n First")
<< " object from pair " << pairIndex
<< ": converted phi index " << newIndex
<< "is equal to badIndex " << badIndex
<< " Conversion failed." << std::endl;
}
}
}
}
break;
default: {
// should not happen (programming error)
throw cms::Exception("FailModule")
<< "\n Wrong position in the object pair " << positionPair
<< "\n Programming error - position must be either 0 or 1..."
<< std::endl;
}
break;
}
//
convertedIndex = newIndex;
return conversionStatus;
}
| const unsigned int L1GtEtaPhiConversions::gtObjectNrBinsPhi | ( | const L1GtObject & | gtObject | ) | const |
return the number of phi bins for a GT object
Definition at line 846 of file L1GtEtaPhiConversions.cc.
References BPTX, CenJet, ETM, ETT, ForJet, GtExternal, HfBitCounts, HfRingEtSums, HTM, HTT, IsoEG, JetCounts, l1GtObjectEnumToString(), m_nrBinsPhiEtm, m_nrBinsPhiHtm, m_nrBinsPhiJetEg, m_nrBinsPhiMu, Mu, NoIsoEG, ObjNull, TauJet, and TechTrig.
Referenced by L1GtCorrelationCondition::evaluateCondition(), and L1GlobalTriggerGTL::run().
{
switch (gtObject) {
case Mu: {
return m_nrBinsPhiMu;
}
break;
case NoIsoEG:
case IsoEG:
case CenJet:
case ForJet:
case TauJet: {
return m_nrBinsPhiJetEg;
}
break;
case ETM: {
return m_nrBinsPhiEtm;
}
break;
case ETT:
case HTT: {
return 0;
}
break;
case HTM: {
return m_nrBinsPhiHtm;
}
break;
case JetCounts:
case HfBitCounts:
case HfRingEtSums:
case TechTrig:
case Castor:
case BPTX:
case GtExternal:
case ObjNull: {
return 0;
}
break;
default: {
edm::LogInfo("L1GtObject") << "\n '"
<< (l1GtObjectEnumToString(gtObject))
<< "' is not a recognized L1GtObject. "
<< "\n Return 0 bins.";
return 0;
}
break;
}
}
| const unsigned int L1GtEtaPhiConversions::gtObjectNrBinsPhi | ( | const L1GtObject & | obj0, |
| const L1GtObject & | obj1 | ||
| ) | const |
return the number of phi bins for a pair of GT objects, according to conversion rules
Definition at line 905 of file L1GtEtaPhiConversions.cc.
References l1GtObjectEnumToString(), LogTrace, m_gtObjectPairVec, and m_pairNrPhiBinsVec.
{
std::pair < L1GtObject, L1GtObject > gtObjPair;
gtObjPair = std::make_pair(obj0, obj1);
//LogTrace("L1GlobalTrigger") << "\nCompute gtObjectNrBinsPhi ["
// << (l1GtObjectEnumToString(obj0)) << ", "
// << (l1GtObjectEnumToString(obj1)) << "]\n" << std::endl;
int iPair = 0;
for (std::vector<std::pair<L1GtObject, L1GtObject> >::const_iterator cIter =
m_gtObjectPairVec.begin(); cIter != m_gtObjectPairVec.end(); ++cIter) {
if (*cIter == gtObjPair) {
LogTrace("L1GlobalTrigger") << "\n gtObjectNrBinsPhi ["
<< l1GtObjectEnumToString(obj0) << ", "
<< l1GtObjectEnumToString(obj1) << "] = "
<< (*(m_pairNrPhiBinsVec.at(iPair))) << std::endl;
return *(m_pairNrPhiBinsVec.at(iPair));
}
iPair++;
}
return 0;
}
| const unsigned int L1GtEtaPhiConversions::gtObjectNrBinsPhi | ( | const unsigned int | pairIndex | ) | const |
return the number of phi bins for a pair of GT objects, according to conversion rules, when the index of the pair is used
Definition at line 934 of file L1GtEtaPhiConversions.cc.
References LogTrace, m_isDebugEnabled, m_pairNrPhiBinsVec, and m_verbosity.
{
if (m_verbosity && m_isDebugEnabled) {
LogTrace("L1GlobalTrigger")
<< "\n gtObjectNrBinsPhi for L1 GT object pair index "
<< pairIndex << " = " << (*(m_pairNrPhiBinsVec.at(pairIndex)))
<< std::endl;
}
return *(m_pairNrPhiBinsVec.at(pairIndex));
}
| const unsigned int L1GtEtaPhiConversions::gtObjectPairIndex | ( | const L1GtObject & | obj0, |
| const L1GtObject & | obj1 | ||
| ) | const |
return the index of a pair in the vector m_gtObjectPairVec, to be used to extract the number of phi bins, the conversion of the indices, etc
Definition at line 460 of file L1GtEtaPhiConversions.cc.
References l1GtObjectEnumToString(), LogTrace, and m_gtObjectPairVec.
Referenced by L1GtCorrelationCondition::evaluateCondition().
{
std::pair < L1GtObject, L1GtObject > gtObjPair;
gtObjPair = std::make_pair(obj0, obj1);
//LogTrace("L1GlobalTrigger") << "\nCompute index for pair ["
// << (l1GtObjectEnumToString(obj0)) << ", "
// << (l1GtObjectEnumToString(obj1)) << "]\n" << std::endl;
unsigned int iPair = 0;
for (std::vector<std::pair<L1GtObject, L1GtObject> >::const_iterator cIter =
m_gtObjectPairVec.begin(); cIter != m_gtObjectPairVec.end(); ++cIter) {
if (*cIter == gtObjPair) {
LogTrace("L1GlobalTrigger") << "\n Index for pair ["
<< l1GtObjectEnumToString(obj0) << ", "
<< l1GtObjectEnumToString(obj1) << "] = "
<< iPair << std::endl;
return iPair;
}
iPair++;
}
// if the pair was not found, return index outside vector size
// it should never happen, except due to programming error
// by using .at one gets an exception when using index outside vector size,
// due to the programming error...
return m_gtObjectPairVec.size();
}
| void L1GtEtaPhiConversions::print | ( | std::ostream & | myCout | ) | const [virtual] |
print all the performed conversions
Definition at line 1382 of file L1GtEtaPhiConversions.cc.
References badIndex, L1CaloGeometry::emJetPhiBinHighEdge(), L1CaloGeometry::emJetPhiBinLowEdge(), L1CaloGeometry::etaBinCenter(), L1CaloGeometry::etSumPhiBinHighEdge(), L1CaloGeometry::etSumPhiBinLowEdge(), L1MuTriggerScales::getGMTEtaScale(), L1MuScale::getHighEdge(), L1MuScale::getLowEdge(), L1MuScale::getNBins(), L1MuTriggerScales::getPhiScale(), L1MuScale::getValue(), L1CaloGeometry::globalEtaBinLowEdge(), L1CaloGeometry::globalEtaIndex(), L1CaloGeometry::htSumPhiBinHighEdge(), L1CaloGeometry::htSumPhiBinLowEdge(), m_l1CaloGeometry, m_l1MuTriggerScales, m_lutEtaCentralToCommonCalo, m_lutEtaForJetToCommonCalo, m_lutEtaMuToCommonCalo, m_lutPhiEtmToHtm, m_lutPhiEtmToJetEg, m_lutPhiHtmToJetEg, m_lutPhiMuToEtm, m_lutPhiMuToHtm, m_lutPhiMuToJetEg, and rad2deg().
Referenced by L1GlobalTriggerGTL::run().
{
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout
<< "\n---++Conversion tables for phi and eta variables of the trigger objects used in correlation conditions \n"
<< std::endl;
//
// phi conversions
//
// phi Mu -> (*Jet, EG)
myCout
<< "\n---+++Phi conversion for muons to jets and e-gamma common phi scale \n"
<< std::endl;
size_t lutPhiMuToJetEgSize = m_lutPhiMuToJetEg.size();
myCout << "Size of look-up table = " << lutPhiMuToJetEgSize << "\n"
<< std::endl;
myCout << "| *Initial Phi Hardware Index* "
<< "|| *Initial Phi Range* ||"
<< " *Converted Phi Hardware Index* "
<< "|| *Converted Phi Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv < lutPhiMuToJetEgSize; ++indexToConv) {
double lowEdgeToConv = m_l1MuTriggerScales->getPhiScale()->getLowEdge(
indexToConv);
double highEdgeToConv =
m_l1MuTriggerScales->getPhiScale()->getHighEdge(indexToConv);
unsigned int convIndex = m_lutPhiMuToJetEg[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->emJetPhiBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->emJetPhiBinHighEdge(convIndex);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << (rad2deg(lowEdgeToConv)) << ", |"
<< std::setw(10) << std::left << (rad2deg(highEdgeToConv))
<< " ) | 0x" << std::setw(6) << std::hex << std::left
<< convIndex << " | " << std::dec << std::setw(6) << convIndex
<< " |[ " << std::setw(10) << std::left << (rad2deg(
convLowEdge)) << ", |" << std::setw(10) << std::left
<< (rad2deg(convHighEdge)) << " ) |" << std::right
<< std::endl;
}
// phi Mu -> ETM
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout << "\n---+++Phi conversion for muons to ETM phi scale \n"
<< std::endl;
size_t lutPhiMuToEtmSize = m_lutPhiMuToEtm.size();
myCout << "Size of look-up table = " << lutPhiMuToEtmSize << "\n"
<< std::endl;
myCout << "| *Initial Phi Hardware Index* "
<< "|| *Initial Phi Range* ||"
<< " *Converted Phi Hardware Index* "
<< "|| *Converted Phi Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv < lutPhiMuToEtmSize; ++indexToConv) {
double lowEdgeToConv = m_l1MuTriggerScales->getPhiScale()->getLowEdge(
indexToConv);
double highEdgeToConv =
m_l1MuTriggerScales->getPhiScale()->getHighEdge(indexToConv);
unsigned int convIndex = m_lutPhiMuToEtm[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->etSumPhiBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->etSumPhiBinHighEdge(convIndex);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << (rad2deg(lowEdgeToConv)) << ", |"
<< std::setw(10) << std::left << (rad2deg(highEdgeToConv))
<< " ) | 0x" << std::setw(6) << std::hex << std::left
<< convIndex << " | " << std::dec << std::setw(6) << convIndex
<< " |[ " << std::setw(10) << std::left << (rad2deg(
convLowEdge)) << ", |" << std::setw(10) << std::left
<< (rad2deg(convHighEdge)) << " ) |" << std::right
<< std::endl;
}
// phi Mu -> HTM
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout << "\n---+++Phi conversion for muons to HTM phi scale \n"
<< std::endl;
size_t lutPhiMuToHtmSize = m_lutPhiMuToHtm.size();
myCout << "Size of look-up table = " << lutPhiMuToHtmSize << "\n"
<< std::endl;
myCout << "| *Initial Phi Hardware Index* "
<< "|| *Initial Phi Range* ||"
<< " *Converted Phi Hardware Index* "
<< "|| *Converted Phi Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv < lutPhiMuToHtmSize; ++indexToConv) {
double lowEdgeToConv = m_l1MuTriggerScales->getPhiScale()->getLowEdge(
indexToConv);
double highEdgeToConv =
m_l1MuTriggerScales->getPhiScale()->getHighEdge(indexToConv);
unsigned int convIndex = m_lutPhiMuToHtm[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->htSumPhiBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->htSumPhiBinHighEdge(convIndex);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << (rad2deg(lowEdgeToConv)) << ", |"
<< std::setw(10) << std::left << (rad2deg(highEdgeToConv))
<< " ) | 0x" << std::setw(6) << std::hex << std::left
<< convIndex << " | " << std::dec << std::setw(6) << convIndex
<< " |[ " << std::setw(10) << std::left << (rad2deg(
convLowEdge)) << ", |" << std::setw(10) << std::left
<< (rad2deg(convHighEdge)) << " ) |" << std::right
<< std::endl;
}
// phi ETM -> (*Jet, EG)
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout
<< "\n---+++Phi conversion for ETM to jets and e-gamma scale common phi scale \n"
<< std::endl;
size_t lutPhiEtmToJetEgSize = m_lutPhiEtmToJetEg.size();
myCout << "Size of look-up table = " << lutPhiEtmToJetEgSize << "\n"
<< std::endl;
myCout << "| *Initial Phi Hardware Index* "
<< "|| *Initial Phi Range* ||"
<< " *Converted Phi Hardware Index* "
<< "|| *Converted Phi Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv < lutPhiEtmToJetEgSize; ++indexToConv) {
double lowEdgeToConv =
m_l1CaloGeometry->etSumPhiBinLowEdge(indexToConv);
double highEdgeToConv = m_l1CaloGeometry->etSumPhiBinHighEdge(
indexToConv);
unsigned int convIndex = m_lutPhiEtmToJetEg[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->emJetPhiBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->emJetPhiBinHighEdge(convIndex);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << (rad2deg(lowEdgeToConv)) << ", |"
<< std::setw(10) << std::left << (rad2deg(highEdgeToConv))
<< " ) | 0x" << std::setw(6) << std::hex << std::left
<< convIndex << " | " << std::dec << std::setw(6) << convIndex
<< " |[ " << std::setw(10) << std::left << (rad2deg(
convLowEdge)) << ", |" << std::setw(10) << std::left
<< (rad2deg(convHighEdge)) << " ) |" << std::right
<< std::endl;
}
// phi ETM -> HTM
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout << "\n---+++Phi conversion for ETM to HTM phi scale \n" << std::endl;
size_t lutPhiEtmToHtmSize = m_lutPhiEtmToHtm.size();
myCout << "Size of look-up table = " << lutPhiEtmToHtmSize << "\n"
<< std::endl;
myCout << "| *Initial Phi Hardware Index* "
<< "|| *Initial Phi Range* ||"
<< " *Converted Phi Hardware Index* "
<< "|| *Converted Phi Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv < lutPhiEtmToHtmSize; ++indexToConv) {
double lowEdgeToConv =
m_l1CaloGeometry->etSumPhiBinLowEdge(indexToConv);
double highEdgeToConv = m_l1CaloGeometry->etSumPhiBinHighEdge(
indexToConv);
unsigned int convIndex = m_lutPhiEtmToHtm[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->htSumPhiBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->htSumPhiBinHighEdge(convIndex);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << (rad2deg(lowEdgeToConv)) << ", |"
<< std::setw(10) << std::left << (rad2deg(highEdgeToConv))
<< " ) | 0x" << std::setw(6) << std::hex << std::left
<< convIndex << " | " << std::dec << std::setw(6) << convIndex
<< " |[ " << std::setw(10) << std::left << (rad2deg(
convLowEdge)) << ", |" << std::setw(10) << std::left
<< (rad2deg(convHighEdge)) << " ) |" << std::right
<< std::endl;
}
// phi HTM -> (*Jet, EG)
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout
<< "\n---+++Phi conversion for HTM to jets and e-gamma scale common phi scale \n"
<< std::endl;
size_t lutPhiHtmToJetEgSize = m_lutPhiHtmToJetEg.size();
myCout << "Size of look-up table = " << lutPhiHtmToJetEgSize << "\n"
<< std::endl;
myCout << "| *Initial Phi Hardware Index* "
<< "|| *Initial Phi Range* ||"
<< " *Converted Phi Hardware Index* "
<< "|| *Converted Phi Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv < lutPhiHtmToJetEgSize; ++indexToConv) {
double lowEdgeToConv =
m_l1CaloGeometry->htSumPhiBinLowEdge(indexToConv);
double highEdgeToConv = m_l1CaloGeometry->htSumPhiBinHighEdge(
indexToConv);
unsigned int convIndex = m_lutPhiHtmToJetEg[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->emJetPhiBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->emJetPhiBinHighEdge(convIndex);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << (rad2deg(lowEdgeToConv)) << ", |"
<< std::setw(10) << std::left << (rad2deg(highEdgeToConv))
<< " ) | 0x" << std::setw(6) << std::hex << std::left
<< convIndex << " | " << std::dec << std::setw(6) << convIndex
<< " |[ " << std::setw(10) << std::left << (rad2deg(
convLowEdge)) << ", |" << std::setw(10) << std::left
<< (rad2deg(convHighEdge)) << " ) |" << std::right
<< std::endl;
}
//
// eta conversions
//
// CenJet/TauJet & IsoEG/NoIsoEG to a common central / forward calorimeter eta scale
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout
<< "\n---+++Eta conversion for central and tau jets and e-gamma objects to a common central and forward calorimeter eta scale \n"
<< std::endl;
size_t lutEtaCentralToCommonCaloSize = m_lutEtaCentralToCommonCalo.size();
myCout << "Size of look-up table = " << lutEtaCentralToCommonCaloSize
<< "\n" << std::endl;
myCout << "| *Initial Eta Hardware Index* "
<< "|| *Initial Eta Range* ||"
<< " *Converted Eta Hardware Index* "
<< "|| *Converted Eta Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv
< lutEtaCentralToCommonCaloSize; ++indexToConv) {
double lowEdgeToConv = m_l1CaloGeometry->globalEtaBinLowEdge(
m_l1CaloGeometry->globalEtaIndex(
m_l1CaloGeometry->etaBinCenter(indexToConv, true)));
double highEdgeToConv = m_l1CaloGeometry->globalEtaBinLowEdge(
m_l1CaloGeometry->globalEtaIndex(
m_l1CaloGeometry->etaBinCenter(indexToConv, true)) + 1);
unsigned int convIndex = m_lutEtaCentralToCommonCalo[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->globalEtaBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->globalEtaBinLowEdge(convIndex + 1);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << lowEdgeToConv << ", |" << std::setw(10)
<< std::left << highEdgeToConv << " ) | 0x" << std::setw(6)
<< std::hex << std::left << convIndex << " | " << std::dec
<< std::setw(6) << convIndex << " |[ " << std::setw(10)
<< std::left << convLowEdge << ", |" << std::setw(10)
<< std::left << convHighEdge << " ) |" << std::right
<< std::endl;
}
// ForJet to a common central / forward calorimeter eta scale
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout
<< "\n---+++Eta conversion for forward jets to a common central and forward calorimeter eta scale \n"
<< std::endl;
size_t lutEtaForJetToCommonCaloSize = m_lutEtaForJetToCommonCalo.size();
myCout << "Size of look-up table = " << lutEtaForJetToCommonCaloSize
<< "\n" << std::endl;
myCout << "| *Initial Eta Hardware Index* "
<< "|| *Initial Eta Range* ||"
<< " *Converted Eta Hardware Index* "
<< "|| *Converted Eta Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv
< lutEtaForJetToCommonCaloSize; ++indexToConv) {
double lowEdgeToConv = m_l1CaloGeometry->globalEtaBinLowEdge(
m_l1CaloGeometry->globalEtaIndex(
m_l1CaloGeometry->etaBinCenter(indexToConv, false)));
double highEdgeToConv =
m_l1CaloGeometry->globalEtaBinLowEdge(
m_l1CaloGeometry->globalEtaIndex(
m_l1CaloGeometry->etaBinCenter(indexToConv,
false)) + 1);
unsigned int convIndex = m_lutEtaForJetToCommonCalo[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->globalEtaBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->globalEtaBinLowEdge(convIndex + 1);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << lowEdgeToConv << ", |" << std::setw(10)
<< std::left << highEdgeToConv << " ) | 0x" << std::setw(6)
<< std::hex << std::left << convIndex << " | " << std::dec
<< std::setw(6) << convIndex << " |[ " << std::setw(10)
<< std::left << convLowEdge << ", |" << std::setw(10)
<< std::left << convHighEdge << " ) |" << std::right
<< std::endl;
}
// Mu to a common central / forward calorimeter eta scale
// force a page break before each group
myCout << "<p style=\"page-break-before: always\"> </p>";
myCout
<< "\n---+++Eta conversion for muons to a common central and forward calorimeter eta scale \n"
<< std::endl;
size_t lutEtaMuToCommonCaloSize = m_lutEtaMuToCommonCalo.size();
myCout << "Size of look-up table = " << lutEtaMuToCommonCaloSize << "\n"
<< std::endl;
unsigned int nrBinsEtaMuPerHalf =
m_l1MuTriggerScales->getGMTEtaScale()->getNBins();
myCout << "| *Initial Eta Hardware Index* "
<< "|| *Initial Eta Range* ||"
<< " *Converted Eta Hardware Index* "
<< "|| *Converted Eta Range* ||" << "\n"
<< "| *hex* | *dec* | ^|^| *hex* | *dec* |^|^|"
<< std::endl;
for (unsigned int indexToConv = 0; indexToConv < lutEtaMuToCommonCaloSize; ++indexToConv) {
// Mu scale defined for positive values only, need to be symmetrized
unsigned int iBinOffset = 0;
double etaSign = 1.;
if (indexToConv > nrBinsEtaMuPerHalf) {
iBinOffset = nrBinsEtaMuPerHalf + 1;
etaSign = -1.;
}
double lowEdgeToConv = etaSign
* m_l1MuTriggerScales->getGMTEtaScale()->getValue(
indexToConv - iBinOffset);
double highEdgeToConv = etaSign
* m_l1MuTriggerScales->getGMTEtaScale()->getValue(
indexToConv + 1 - iBinOffset);
unsigned int convIndex = m_lutEtaMuToCommonCalo[indexToConv];
double convLowEdge = 0.;
double convHighEdge = 0.;
if (convIndex != badIndex) {
convLowEdge = m_l1CaloGeometry->globalEtaBinLowEdge(convIndex);
convHighEdge = m_l1CaloGeometry->globalEtaBinLowEdge(convIndex + 1);
} else {
// badIndex means a bad initialIndex
lowEdgeToConv = 0.;
highEdgeToConv = 0.;
}
myCout << "| 0x" << std::setw(3) << std::hex << std::left
<< indexToConv << " | " << std::dec << std::setw(3)
<< std::left << indexToConv << " |[ " << std::setw(10)
<< std::left << lowEdgeToConv << ", |" << std::setw(10)
<< std::left << highEdgeToConv << " ) | 0x" << std::setw(6)
<< std::hex << std::left << convIndex << " | " << std::dec
<< std::setw(6) << convIndex << " |[ " << std::setw(10)
<< std::left << convLowEdge << ", |" << std::setw(10)
<< std::left << convHighEdge << " ) |" << std::right
<< std::endl;
}
}
| const double L1GtEtaPhiConversions::rad2deg | ( | const double & | phiRad | ) | const [private] |
convert phi from rad (-pi, pi] to deg (0, 360)
Definition at line 1910 of file L1GtEtaPhiConversions.cc.
References PiConversion.
Referenced by print().
{
if (phiRad < 0.) {
return (phiRad * PiConversion) + 360.;
} else {
return (phiRad * PiConversion);
}
}
| void L1GtEtaPhiConversions::setVerbosity | ( | const int | verbosity | ) | [inline] |
Definition at line 86 of file L1GtEtaPhiConversions.h.
References m_verbosity, and verbosity.
Referenced by L1GlobalTriggerGTL::L1GlobalTriggerGTL(), and L1GlobalTriggerGTL::run().
{
m_verbosity = verbosity;
}
const unsigned int L1GtEtaPhiConversions::badIndex = 999999 [static, private] |
a bad index value, treated specially when performing the conversions or printing the conversion vectors
Definition at line 100 of file L1GtEtaPhiConversions.h.
Referenced by convertEtaIndex(), convertL1Scales(), convertPhiIndex(), and print().
std::vector<std::pair<L1GtObject, L1GtObject> > L1GtEtaPhiConversions::m_gtObjectPairVec [private] |
vector of all L1GtObject pairs
Definition at line 112 of file L1GtEtaPhiConversions.h.
Referenced by gtObjectNrBinsPhi(), gtObjectPairIndex(), and L1GtEtaPhiConversions().
bool L1GtEtaPhiConversions::m_isDebugEnabled [private] |
cached edm::isDebugEnabled()
Definition at line 199 of file L1GtEtaPhiConversions.h.
Referenced by convertEtaIndex(), convertL1Scales(), convertPhiIndex(), gtObjectNrBinsPhi(), and L1GtEtaPhiConversions().
const L1CaloGeometry* L1GtEtaPhiConversions::m_l1CaloGeometry [private] |
pointer to calorimetry scales - updated in convertl1Scales method
Definition at line 134 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), and print().
const L1MuTriggerScales* L1GtEtaPhiConversions::m_l1MuTriggerScales [private] |
pointer to muon scales - updated in convertl1Scales method
Definition at line 137 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutEtaCentralToCommonCalo [private] |
eta conversion of CenJet/TauJet & IsoEG/NoIsoEG to a common calorimeter eta scale
Definition at line 185 of file L1GtEtaPhiConversions.h.
Referenced by convertEtaIndex(), convertL1Scales(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutEtaForJetToCommonCalo [private] |
eta conversion of ForJet to the common calorimeter eta scale defined before
Definition at line 188 of file L1GtEtaPhiConversions.h.
Referenced by convertEtaIndex(), convertL1Scales(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutEtaMuToCommonCalo [private] |
eta conversion of Mu to the common calorimeter eta scale defined before
Definition at line 191 of file L1GtEtaPhiConversions.h.
Referenced by convertEtaIndex(), convertL1Scales(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutPhiEtmToHtm [private] |
phi conversion for ETM to HTM
Definition at line 171 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), L1GtEtaPhiConversions(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutPhiEtmToJetEg [private] |
phi conversion for ETM to (*Jet, EG)
Definition at line 168 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), L1GtEtaPhiConversions(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutPhiHtmToJetEg [private] |
phi conversion for HTM to (*Jet, EG)
Definition at line 174 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), L1GtEtaPhiConversions(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutPhiJetEgToJetEg [private] |
phi conversion for (*Jet, EG) to (*Jet, EG) return the same index as the input index, introduced only to simplify convertPhiIndex
Definition at line 179 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), and L1GtEtaPhiConversions().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutPhiMuToEtm [private] |
phi conversion for Mu to ETM
Definition at line 162 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), L1GtEtaPhiConversions(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutPhiMuToHtm [private] |
phi conversion for Mu to HTM
Definition at line 165 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), L1GtEtaPhiConversions(), and print().
std::vector<unsigned int> L1GtEtaPhiConversions::m_lutPhiMuToJetEg [private] |
phi conversion for Mu to (*Jet, EG)
Definition at line 159 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), L1GtEtaPhiConversions(), and print().
unsigned int L1GtEtaPhiConversions::m_nrBinsEtaCommon [private] |
number of eta bins for common scale
Definition at line 154 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales().
unsigned int L1GtEtaPhiConversions::m_nrBinsPhiEtm [private] |
number of phi bins for ETM
Definition at line 148 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), gtObjectNrBinsPhi(), and L1GtEtaPhiConversions().
unsigned int L1GtEtaPhiConversions::m_nrBinsPhiHtm [private] |
number of phi bins for HTM
Definition at line 151 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), gtObjectNrBinsPhi(), and L1GtEtaPhiConversions().
unsigned int L1GtEtaPhiConversions::m_nrBinsPhiJetEg [private] |
number of phi bins for calorimeter objects (*Jet, *EG)
Definition at line 145 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), gtObjectNrBinsPhi(), and L1GtEtaPhiConversions().
unsigned int L1GtEtaPhiConversions::m_nrBinsPhiMu [private] |
number of phi bins for muons
Definition at line 142 of file L1GtEtaPhiConversions.h.
Referenced by convertL1Scales(), and gtObjectNrBinsPhi().
std::vector<bool> L1GtEtaPhiConversions::m_pairConvertPhiFirstGtObject [private] |
decide which object to convert: if m_pairConvertPhiFirstGtObject true, convert pair.first and do not convert pair.second if m_pairConvertPhiFirstGtObject false, do not convert pair.first and convert pair.second
Definition at line 117 of file L1GtEtaPhiConversions.h.
Referenced by convertPhiIndex(), and L1GtEtaPhiConversions().
std::vector<const unsigned int*> L1GtEtaPhiConversions::m_pairNrPhiBinsVec [private] |
number of phi bins for each L1GtObject pair in the scale used for that pair it is filled correlated with m_gtObjectPairVec, so the index of the pair in m_gtObjectPairVec is the index of the m_pairNrPhiBinsVec element containing the number of phi bins
Definition at line 123 of file L1GtEtaPhiConversions.h.
Referenced by gtObjectNrBinsPhi(), and L1GtEtaPhiConversions().
std::vector<const std::vector<unsigned int>*> L1GtEtaPhiConversions::m_pairPhiConvVec [private] |
constant references to conversion LUT for a given L1GtObject pair it is filled correlated with m_gtObjectPairVec, so the index of the pair in m_gtObjectPairVec is the index of the m_pairPhiConvVec element containing the reference
Definition at line 129 of file L1GtEtaPhiConversions.h.
Referenced by convertPhiIndex(), and L1GtEtaPhiConversions().
int L1GtEtaPhiConversions::m_verbosity [private] |
verbosity level
Definition at line 196 of file L1GtEtaPhiConversions.h.
Referenced by convertEtaIndex(), convertL1Scales(), convertPhiIndex(), gtObjectNrBinsPhi(), and setVerbosity().
const double L1GtEtaPhiConversions::PiConversion = 180. / acos(-1.) [static, private] |
Definition at line 103 of file L1GtEtaPhiConversions.h.
Referenced by rad2deg().
1.7.3