|
|
|
#include <FastL1GlobalAlgo.h>
Description: Global algorithm.
Implementation: <Notes on="" implementation>="">
Definition at line 70 of file FastL1GlobalAlgo.h.
| FastL1GlobalAlgo::FastL1GlobalAlgo | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 28 of file FastL1GlobalAlgo.cc.
References L1Config::CrystalEBThreshold, L1Config::CrystalEEThreshold, L1Config::DoEMCorr, L1Config::DoJetCorr, L1Config::EcalTPInput, L1Config::EMActiveLevel, L1Config::EmInputs, L1Config::EMLSB, L1Config::EMNoiseLevel, L1Config::EMSeedEnThreshold, L1Config::FGEBThreshold, L1Config::FGEEThreshold, edm::FileInPath::fullPath(), FastL1RegionMap::getFastL1RegionMap(), edm::ParameterSet::getParameter(), L1Config::HadActiveLevel, L1Config::HadNoiseLevel, L1Config::HcalLUT, L1Config::HcalTPInput, L1Config::hOeThreshold, i, j, L1Config::JetLSB, L1Config::JetSeedEtThreshold, m_DoBitInfo, m_GctIso, m_hcaluncomp, m_IsolationEt, m_L1Config, m_Regions, m_RMap, L1Config::MuonNoiseLevel, L1Config::noFGThreshold, L1Config::noTauVetoLevel, L1Config::QuietRegionThreshold, L1Config::TowerEBScale, L1Config::TowerEBThreshold, L1Config::TowerEEScale, L1Config::TowerEEThreshold, L1Config::TowerEMLSB, L1Config::TowerHadLSB, L1Config::TowerHBScale, L1Config::TowerHBThreshold, L1Config::TowerHEScale, L1Config::TowerHEThreshold, and L1Config::TowerInput.
{
m_RMap = FastL1RegionMap::getFastL1RegionMap();
// Get L1 config
m_L1Config.DoEMCorr = iConfig.getParameter<bool>("DoEMCorr");
m_L1Config.DoJetCorr = iConfig.getParameter<bool>("DoJetCorr");
m_DoBitInfo = iConfig.getParameter<bool>("DoBitInfo");
m_GctIso = iConfig.getParameter<bool>("GctIso");
m_IsolationEt = iConfig.getParameter<double>("IsolationEt");
// get uncompressed hcal et
m_L1Config.HcalLUT = iConfig.getParameter<edm::FileInPath>("HcalLUT");
m_L1Config.EMSeedEnThreshold = iConfig.getParameter<double>("EMSeedEnThreshold");
m_L1Config.EMActiveLevel = iConfig.getParameter<double>("EMActiveLevel");
m_L1Config.HadActiveLevel = iConfig.getParameter<double>("HadActiveLevel");
m_L1Config.noTauVetoLevel = iConfig.getParameter<double>("noTauVetoLevel");
m_L1Config.hOeThreshold = iConfig.getParameter<double>("hOeThreshold");
m_L1Config.FGEBThreshold = iConfig.getParameter<double>("FGEBThreshold");
m_L1Config.FGEEThreshold = iConfig.getParameter<double>("FGEEThreshold");
m_L1Config.MuonNoiseLevel = iConfig.getParameter<double>("MuonNoiseLevel");
m_L1Config.EMNoiseLevel = iConfig.getParameter<double>("EMNoiseLevel");
m_L1Config.HadNoiseLevel = iConfig.getParameter<double>("HadNoiseLevel");
m_L1Config.QuietRegionThreshold = iConfig.getParameter<double>("QuietRegionThreshold");
m_L1Config.JetSeedEtThreshold = iConfig.getParameter<double>("JetSeedEtThreshold");
m_L1Config.TowerEMLSB = iConfig.getParameter<double>("TowerEMLSB");
m_L1Config.TowerHadLSB = iConfig.getParameter<double>("TowerHadLSB");
m_L1Config.EMLSB = iConfig.getParameter<double>("EMLSB");
m_L1Config.JetLSB = iConfig.getParameter<double>("JetLSB");
m_L1Config.CrystalEBThreshold = iConfig.getParameter<double>("CrystalEBThreshold");
m_L1Config.CrystalEEThreshold = iConfig.getParameter<double>("CrystalEEThreshold");
m_L1Config.TowerEBThreshold = iConfig.getParameter<double>("TowerEBThreshold");
m_L1Config.TowerEEThreshold = iConfig.getParameter<double>("TowerEEThreshold");
m_L1Config.TowerHBThreshold = iConfig.getParameter<double>("TowerHBThreshold");
m_L1Config.TowerHEThreshold = iConfig.getParameter<double>("TowerHEThreshold");
m_L1Config.TowerEBScale = iConfig.getParameter<double>("TowerEBScale");
m_L1Config.TowerEEScale = iConfig.getParameter<double>("TowerEEScale");
m_L1Config.TowerHBScale = iConfig.getParameter<double>("TowerHBScale");
m_L1Config.TowerHEScale = iConfig.getParameter<double>("TowerHEScale");
m_L1Config.noFGThreshold = iConfig.getParameter<double>("noFGThreshold");
m_L1Config.EmInputs = iConfig.getParameter <std::vector<edm::InputTag> >("EmInputs");
m_L1Config.TowerInput = iConfig.getParameter<edm::InputTag>("TowerInput");
m_L1Config.EcalTPInput = iConfig.getParameter<edm::InputTag>("EcalTPInput");
m_L1Config.HcalTPInput = iConfig.getParameter<edm::InputTag>("HcalTPInput");
//Load Hcal LUT file
std::ifstream userfile;
const std::string userfileName = m_L1Config.HcalLUT.fullPath();
userfile.open(userfileName.c_str());
static const int etabound = 32;
static const int tpgmax = 256;
for (int i=0; i<tpgmax; i++) {
for(int j = 1; j <=etabound; j++) {
userfile >> m_hcaluncomp[j][i];
}
}
userfile.close();
FastL1Region tr;
for (int i=0;i<396;i++)
m_Regions.push_back(tr);
}
| FastL1GlobalAlgo::~FastL1GlobalAlgo | ( | ) |
Definition at line 103 of file FastL1GlobalAlgo.cc.
{
}
| void FastL1GlobalAlgo::addJet | ( | int | rgnId, |
| bool | taubit | ||
| ) | [private] |
Definition at line 208 of file FastL1GlobalAlgo.cc.
References abs, corrJetEt(), funct::cos(), L1Config::DoJetCorr, ExpressReco_HICollisions_FallBack::e, ExpressReco_HICollisions_FallBack::et, eta(), funct::exp(), GCTEnergyTrunc(), FastL1RegionMap::getRegionCenterEtaPhi(), greaterEt(), L1Config::JetLSB, m_CenJets, m_DoBitInfo, m_ForJets, m_L1Config, m_Regions, m_RMap, m_TauJets, L1Config::noTauVetoLevel, L1TEmulatorMonitor_cff::p, phi, funct::sin(), python::multivaluedict::sort(), and theta().
Referenced by findJets().
{
std::pair<double, double> p = m_RMap->getRegionCenterEtaPhi(iRgn);
//double e = m_Regions.at(iRgn).GetJetE();
//double et = GCTEnergyTrunc(m_Regions.at(iRgn).GetJetEt(), m_L1Config.JetLSB, false);
double et = m_Regions.at(iRgn).GetJetEt();
double eta = p.first;
double phi = p.second;
if (m_L1Config.DoJetCorr) {
et = GCTEnergyTrunc(corrJetEt(et,eta), m_L1Config.JetLSB, false);
} else {
et = GCTEnergyTrunc(et, m_L1Config.JetLSB, false);
}
double theta = 2.*atan(exp(-eta));
double ex = et*cos(phi);
double ey = et*sin(phi);
//double ex = e*sin(theta)*cos(phi);
//double ey = e*sin(theta)*sin(phi);
double e = ex/sin(theta)/cos(phi);
double ez = e*cos(theta);
if (m_DoBitInfo){
m_Regions[iRgn].BitInfo.setEt(et);
m_Regions[iRgn].BitInfo.setEnergy(e);
}
reco::Particle::LorentzVector rp4(ex,ey,ez,e);
l1extra::L1JetParticle tjet(rp4);
if (et>=5.) {
if (m_DoBitInfo) m_Regions[iRgn].BitInfo.setHard(true);
if (m_DoBitInfo) m_Regions[iRgn].BitInfo.setSoft(false);
if ((taubit || et>m_L1Config.noTauVetoLevel) && (std::abs(eta)<3.0) ) {
m_TauJets.push_back(tjet);
// sort by et
std::sort(m_TauJets.begin(),m_TauJets.end(), myspace::greaterEt);
} else {
if (std::abs(eta)<3.0) {
m_CenJets.push_back(tjet);
std::sort(m_CenJets.begin(),m_CenJets.end(), myspace::greaterEt);
} else {
m_ForJets.push_back(tjet);
std::sort(m_ForJets.begin(),m_ForJets.end(), myspace::greaterEt);
}
}
}
else{
if (m_DoBitInfo) m_Regions[iRgn].BitInfo.setSoft(true);
if (m_DoBitInfo) m_Regions[iRgn].BitInfo.setHard(false);
}
}
| void FastL1GlobalAlgo::CaloTowersDump | ( | edm::Event const & | e | ) |
Definition at line 111 of file FastL1GlobalAlgo.cc.
References edm::Event::getByLabel(), collect_tpl::input, j, m_L1Config, and L1Config::TowerInput.
{
/*
std::vector<edm::Handle<CaloTowerCollection> > prods;
edm::LogInfo("FastL1GlobalAlgo::CaloTowersDump") << "Start!";
e.getManyByType(prods);
std::vector<edm::Handle<CaloTowerCollection> >::iterator i;
for (i=prods.begin(); i!=prods.end(); i++) {
const CaloTowerCollection& c=*(*i);
for (CaloTowerCollection::const_iterator j=c.begin(); j!=c.end(); j++) {
edm::LogInfo("FastL1GlobalAlgo::CaloTowersDump") << *j;
}
}
edm::LogInfo("FastL1GlobalAlgo::CaloTowersDump") << "End!";
*/
edm::Handle<CaloTowerCollection> input;
edm::LogInfo("FastL1GlobalAlgo::CaloTowersDump") << "Start!";
e.getByLabel(m_L1Config.TowerInput,input);
for (CaloTowerCollection::const_iterator j=input->begin(); j!=input->end(); j++) {
edm::LogInfo("FastL1GlobalAlgo::CaloTowersDump") << *j;
}
edm::LogInfo("FastL1GlobalAlgo::CaloTowersDump") << "End!";
}
| void FastL1GlobalAlgo::checkMapping | ( | ) | [private] |
Definition at line 1341 of file FastL1GlobalAlgo.cc.
References FastL1RegionMap::convertFromECal_to_HCal_iphi(), gather_cfg::cout, FastL1RegionMap::getRegionCenterEtaPhi(), FastL1RegionMap::getRegionEtaPhiIndex(), FastL1RegionMap::getRegionIndex(), FastL1RegionMap::getRegionTowerIndex(), i, j, and m_RMap.
{
// loop over towers ieta,iphi
for (int j=1;j<=72;j++) {
for (int i=-28; i<=28; i++) {
if (i==0) continue;
int iRgn = m_RMap->getRegionIndex(i,j);
std::pair<double, double> RgnEtaPhi = m_RMap->getRegionCenterEtaPhi(iRgn);
//int iTwr = m_RMap->getRegionTowerIndex(i,j);
std::pair<int, int> iRgnEtaPhi = m_RMap->getRegionEtaPhiIndex(iRgn);
std::pair<int, int> iRgnEtaPhi2 = m_RMap->getRegionEtaPhiIndex(std::pair<int, int>(i,j));
std::cout<<"---------------------------------------------------------------------------"<<std::endl;
std::cout<<"Region: "<<iRgn<<" | "<<RgnEtaPhi.first<<", "<<RgnEtaPhi.second*180./3.141<<std::endl;
std::cout<<" - "<<iRgnEtaPhi.first<<", "<<iRgnEtaPhi.second<<std::endl;
std::cout<<" - "<<iRgnEtaPhi2.first<<", "<<iRgnEtaPhi2.second<<std::endl;
std::cout<<" Tower: "<<i<<", "<<m_RMap->convertFromECal_to_HCal_iphi(j)<<std::endl;
std::cout<<" TowerId: "<<m_RMap->getRegionTowerIndex(i,j)<<std::endl;
}
}
}
| void FastL1GlobalAlgo::FillBitInfos | ( | ) |
Definition at line 832 of file FastL1GlobalAlgo.cc.
References i, m_BitInfos, m_DoBitInfo, and m_Regions.
Referenced by FastL1CaloSim::produce().
{
if (m_DoBitInfo){
m_BitInfos.clear();
for (int i=0; i<396; i++) {
m_BitInfos.push_back(m_Regions[i].getBitInfo());
}
}
}
| void FastL1GlobalAlgo::FillEgammas | ( | edm::Event const & | e | ) |
Definition at line 304 of file FastL1GlobalAlgo.cc.
References edm::Event::getByLabel(), greaterEt(), collect_tpl::input, isEMCand(), m_Egammas, m_isoEgammas, m_L1Config, python::multivaluedict::sort(), and L1Config::TowerInput.
Referenced by FastL1CaloSim::produce().
{
m_Egammas.clear();
m_isoEgammas.clear();
//std::vector< edm::Handle<CaloTowerCollection> > input;
//e.getManyByType(input);
edm::Handle<CaloTowerCollection> input;
e.getByLabel(m_L1Config.TowerInput,input);
//std::vector< edm::Handle<CaloTowerCollection> >::iterator j;
//for (j=input.begin(); j!=input.end(); j++) {
// const CaloTowerCollection& c=*(*j);
l1extra::L1EmParticle* ph = new l1extra::L1EmParticle();
//l1extra::L1EmParticle ph;
//for (CaloTowerCollection::const_iterator cnd=c.begin(); cnd!=c.end(); cnd++) {
for (CaloTowerCollection::const_iterator cnd=input->begin(); cnd!=input->end(); cnd++) {
reco::Particle::LorentzVector rp4(0.,0.,0.,0.);
//l1extra::L1EmParticle* ph = new l1extra::L1EmParticle(rp4);
*ph = l1extra::L1EmParticle(rp4);
CaloTowerDetId cid = cnd->id();
int emTag = isEMCand(cid,ph,e);
// 1 = non-iso EM, 2 = iso EM
if (emTag==1) {
m_Egammas.push_back(*ph);
} else if (emTag==2) {
m_isoEgammas.push_back(*ph);
}
}
//}
std::sort(m_Egammas.begin(),m_Egammas.end(), myspace::greaterEt);
std::sort(m_isoEgammas.begin(),m_isoEgammas.end(), myspace::greaterEt);
delete ph;
}
| void FastL1GlobalAlgo::FillEgammasTP | ( | edm::Event const & | e | ) |
Definition at line 266 of file FastL1GlobalAlgo.cc.
References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), greaterEt(), i, isEMCand(), m_Egammas, m_isoEgammas, m_Regions, and python::multivaluedict::sort().
Referenced by FastL1CaloSim::produce().
{
m_Egammas.clear();
m_isoEgammas.clear();
l1extra::L1EmParticle* ph = new l1extra::L1EmParticle();
//l1extra::L1EmParticle ph;
for (int i=0; i<396; i++) {
CaloTowerCollection towers = m_Regions[i].GetCaloTowers();
for (CaloTowerCollection::const_iterator cnd=towers.begin(); cnd!=towers.end(); cnd++) {
if (cnd->emEt()<0.01 && cnd->hadEt()<0.01) continue;
//if (cnd->emEt()<0.01) continue;
reco::Particle::LorentzVector rp4(0.,0.,0.,0.);
// l1extra::L1EmParticle* ph = new l1extra::L1EmParticle(rp4);
*ph = l1extra::L1EmParticle(rp4);
CaloTowerDetId cid = cnd->id();
int emTag = isEMCand(cid,ph,e);
// 1 = non-iso EM, 2 = iso EM
if (emTag==1) {
m_Egammas.push_back(*ph);
} else if (emTag==2) {
m_isoEgammas.push_back(*ph);
}
}
std::sort(m_Egammas.begin(),m_Egammas.end(), myspace::greaterEt);
std::sort(m_isoEgammas.begin(),m_isoEgammas.end(), myspace::greaterEt);
}
delete ph;
}
| void FastL1GlobalAlgo::FillJets | ( | const edm::EventSetup & | e | ) | [inline] |
Definition at line 93 of file FastL1GlobalAlgo.h.
References findJets().
Referenced by FastL1CaloSim::produce().
{ findJets(); };
| void FastL1GlobalAlgo::FillL1Regions | ( | edm::Event const & | e, |
| const edm::EventSetup & | c | ||
| ) |
Definition at line 750 of file FastL1GlobalAlgo.cc.
References abs, L1Config::EmInputs, edm::EventSetup::get(), edm::Event::getByLabel(), FastL1RegionMap::getRegionEtaPhiIndex(), FastL1RegionMap::getRegionTowerIndex(), i, InitL1Regions(), collect_tpl::input, m_L1Config, m_Regions, m_RMap, edm::ESHandle< T >::product(), and L1Config::TowerInput.
Referenced by FastL1CaloSim::produce().
{
InitL1Regions();
edm::Handle<CaloTowerCollection> input;
e.getByLabel(m_L1Config.TowerInput,input);
edm::ESHandle<CaloTowerConstituentsMap> cttopo;
c.get<IdealGeometryRecord>().get(cttopo);
const CaloTowerConstituentsMap* theTowerConstituentsMap = cttopo.product();
edm::ESHandle<CaloTopology> calotopo;
c.get<CaloTopologyRecord>().get(calotopo);
edm::ESHandle<CaloGeometry> cGeom;
//c.get<IdealGeometryRecord>().get(cGeom);
c.get<CaloGeometryRecord>().get(cGeom);
edm::Handle<EcalRecHitCollection> ec1;
e.getByLabel(m_L1Config.EmInputs.at(1),ec1);
edm::Handle<EcalRecHitCollection> ec0;
e.getByLabel(m_L1Config.EmInputs.at(0),ec0);
// ascii visualisation of mapping
//m_RMap->display();
// works for barrel/endcap region only right now!
//std::vector< edm::Handle<CaloTowerCollection> > input;
//e.getManyByType(input);
//edm::Handle<CaloTowerCollection> input;
//e.getByLabel(m_L1Config.TowerInput,input);
//std::vector< edm::Handle<CaloTowerCollection> >::iterator j;
//for (j=input.begin(); j!=input.end(); j++) {
// const CaloTowerCollection& c=*(*j);
// for (CaloTowerCollection::const_iterator cnd=c.begin(); cnd!=c.end(); cnd++) {
for (CaloTowerCollection::const_iterator cnd=input->begin(); cnd!=input->end(); cnd++) {
CaloTowerDetId cid = cnd->id();
std::pair<int, int> pep = m_RMap->getRegionEtaPhiIndex(cid);
int rgnid = 999;
int twrid = 999;
if (std::abs(pep.first)<=22) {
rgnid = pep.second*22 + pep.first;
twrid = m_RMap->getRegionTowerIndex(cid);
//std::cout << rgnid << " " << twrid << std::endl;
//std::cout << cnd->emEt() << " " << cnd->hadEt() << std::endl;
}
if (rgnid<396 && twrid<16) {
m_Regions[rgnid].FillTower_Scaled(*cnd,twrid,true,cGeom);
m_Regions[rgnid].SetRegionBits(e);
}
}
//}
// Fill EM Crystals
for (int i=0; i<396; i++) {
//m_Regions[i].FillEMCrystals(e,iConfig,m_RMap);
m_Regions[i].FillEMCrystals(theTowerConstituentsMap,
&(*calotopo),
&(*cGeom),
&(*ec0), &(*ec1),
m_RMap);
}
//checkMapping();
}
| void FastL1GlobalAlgo::FillL1RegionsTP | ( | edm::Event const & | e, |
| const edm::EventSetup & | c | ||
| ) |
Definition at line 531 of file FastL1GlobalAlgo.cc.
References abs, corrEmEt(), L1Config::DoEMCorr, L1Config::EcalTPInput, L1Config::EMLSB, ExpressReco_HICollisions_FallBack::et, eta(), fastmath::etaphi(), edm::EventSetup::get(), edm::Event::getByLabel(), FastL1RegionMap::getRegionCenterEtaPhi(), FastL1RegionMap::getRegionIndex(), FastL1RegionMap::getRegionTowerIndex(), hcaletValue(), L1Config::HcalTPInput, L1Config::hOeThreshold, i, InitL1Regions(), j, m_L1Config, m_Regions, m_RMap, L1Config::noFGThreshold, phi, and matplotRender::t.
Referenced by FastL1CaloSim::produce().
{
//edm::ESHandle<CaloTowerConstituentsMap> cttopo;
//s.get<IdealGeometryRecord>().get(cttopo);
//const CaloTowerConstituentsMap* theTowerConstituentsMap = cttopo.product();
InitL1Regions();
edm::ESHandle<CaloGeometry> cGeom;
//c.get<IdealGeometryRecord>().get(cGeom);
s.get<CaloGeometryRecord>().get(cGeom);
edm::Handle<EcalTrigPrimDigiCollection> ETPinput;
e.getByLabel(m_L1Config.EcalTPInput,ETPinput);
edm::Handle<HcalTrigPrimDigiCollection> HTPinput;
e.getByLabel(m_L1Config.HcalTPInput,HTPinput);
//CaloTowerCollection towers;
int hEtV [396][16] = {{0}};
int hFGV [396][16] = {{0}};
int hiEtaV[396][16] = {{0}};
int hiPhiV[396][16] = {{0}};
for (HcalTrigPrimDigiCollection::const_iterator hTP=HTPinput->begin();
hTP!=HTPinput->end(); hTP++) {
int rgnid = 999;
int twrid = 999;
int hiphi = hTP->id().iphi();
int hieta = hTP->id().ieta();
/*
if(abs(hieta) > 20 && hiphi%2 == 0) hiphi--;
std::pair<int, int> prim_tower_feed; // prim tower indeces iEta +/- 1~32, iPhi (+) 1~72
prim_tower_feed.first = hieta;
prim_tower_feed.second = hiphi;
std::pair<int, int> rct_index = m_RMap-> getRegionEtaPhiIndex(prim_tower_feed); // convert prim tower indeces into RCT indeces ieta 0~21, iphi 0~17
rgnid = rct_index.second*22 + rct_index.first; // converting fastL1 obsolete RCT numbering
*/
rgnid = m_RMap->getRegionIndex(hieta,hiphi);
twrid = m_RMap->getRegionTowerIndex(hieta,hiphi);
/*
if (std::abs(htwrid.ieta())>28) {
std::cout<<htwrid.ieta()<<" "<<htwrid.iphi()<<" "<<rgnid<<" "<<twrid<<" "<<std::endl;
}
*/
/*
if (hTP->SOI_compressedEt()>0) {
std::cout<<"hcalTP >>> "<<hTP->SOI_compressedEt()<<" "<<hTP->SOI_fineGrain()<<" "
<<rgnid<<" "<<twrid<<std::endl;
}
*/
if(rgnid < 396 && twrid < 16){
hEtV[rgnid][twrid] = (int)hTP->SOI_compressedEt();
hFGV[rgnid][twrid] = (int)hTP->SOI_fineGrain();
hiEtaV[rgnid][twrid] = hieta;
hiPhiV[rgnid][twrid] = hiphi;
}
}
int emEtV [396][16] = {{0}};
int emFGV [396][16] = {{0}};
int emiEtaV[396][16] = {{0}};
int emiPhiV[396][16] = {{0}};
//double emEtaV[396][16] = {0.};
for (EcalTrigPrimDigiCollection::const_iterator eTP=ETPinput->begin();
eTP!=ETPinput->end(); eTP++) {
int eieta = eTP->id().ieta();
if(abs(eieta)> 28) continue; // no crystal in HF
else{
int eiphi = eTP->id().iphi();
//int teiphi = eiphi;
int rgnid = 999;
int twrid = 999;
//if(abs(eieta) > 20 && eiphi%2 == 0) teiphi=eiphi-1;
/*
if(abs(eieta) > 20 && eiphi%2 == 0) eiphi--;
std::pair<int, int> prim_tower_feed; // prim tower indeces iEta +/- 1~28, iPhi (+) 1~72
prim_tower_feed.first = eieta;
prim_tower_feed.second = eiphi;
std::pair<int, int> rct_index = m_RMap-> getRegionEtaPhiIndex(prim_tower_feed); // convert prim tower indeces into RCT indeces ieta 0~21, iphi 0~17
rgnid = rct_index.second*22 + rct_index.first; // converting fastL1 obsolete RCT numbering
*/
/*
rgnid = m_RMap->getRegionIndex(eieta,eiphi);
twrid = m_RMap->getRegionTowerIndex(eieta,eiphi);
if (eTP->compressedEt()>0) {
std::cout<<"ecalTP *** "<<eTP->compressedEt()<<" "<<eTP->fineGrain()<<" "
<<rgnid<<" "<<twrid<<std::endl;
}
*/
if(rgnid < 396 && twrid < 16){
emEtV[rgnid][twrid] = (int)eTP->compressedEt();
emFGV[rgnid][twrid] = (int)eTP->fineGrain();
emiEtaV[rgnid][twrid] = eieta;
emiPhiV[rgnid][twrid] = eiphi;
//edm::ESHandle<CaloGeometry> cGeom;
//s.get<IdealGeometryRecord>().get(cGeom);
//CaloTowerDetId towerDetId = CaloTowerDetId( eieta, teiphi);
//CaloTowerDetId towerDetId2 = CaloTowerDetId( eieta, eiphi);
//const GlobalPoint gP1 = cGeom->getPosition(towerDetId);
//const GlobalPoint gP12 = cGeom->getPosition(towerDetId2);
//double eta = gP1.eta();
//if(abs(eieta) > 20)
//std::cout<<eiphi<<" "<<gP1.phi()<<" "<<eieta<<" "<<gP1.eta()<<" "<<std::endl;
//double phi = gP1.phi();
//emEtaV[rgnid][twrid] = eta;
}
} // else
} // ecalTP
for (int i=0; i<396; i++) {
for (int j=0; j<16; j++) {
if (emEtV[i][j]>0 || hEtV[i][j]>0) {
std::pair<double, double> etaphi
= m_RMap->getRegionCenterEtaPhi(i);
double eta = etaphi.first;
double phi = etaphi.second;
double emEt = ((double) emEtV[i][j]) * m_L1Config.EMLSB;
//double hadEt = ((double )hEtV[i][j]) * m_L1Config.JetLSB * cos(2.*atan(exp(-hiEtaV[i][j])));
int iAbsTwrEta = std::abs(hiEtaV[i][j]);
//double hadEt = ((double )hcaletValue(iAbsTwrEta, hEtV[i][j])) * m_L1Config.JetLSB;
double hadEt = ((double )hcaletValue(iAbsTwrEta, hEtV[i][j]));
//double hadEt = hEtV[i][j] * m_L1Config.JetLSB;
if (m_L1Config.DoEMCorr) {
//emEt = corrEmEt(emEt,emEtaV[i][j]);
emEt = corrEmEt(emEt,std::abs(emiEtaV[i][j]));
}
double et = emEt + hadEt;
edm::ESHandle<CaloGeometry> cGeom;
//s.get<IdealGeometryRecord>().get(cGeom);
s.get<CaloGeometryRecord>().get(cGeom);
//math::RhoEtaPhiVector lvec(et,eta,phi);
math::XYZTLorentzVector lvec(et,eta,phi,0.);
//math::PtEtaPhiMLorentzVector lvec(et,eta,phi,0.);
CaloTowerDetId towerDetId;
if (emEtV[i][j]>0)
towerDetId = CaloTowerDetId(emiEtaV[i][j],emiPhiV[i][j]);
else
towerDetId = CaloTowerDetId(hiEtaV[i][j],hiPhiV[i][j]);
GlobalPoint gP = cGeom->getPosition(towerDetId);
CaloTower t = CaloTower(towerDetId,
emEt, hadEt, 0.,
0, 0, lvec,
gP, gP);
//m_Regions[i].FillTower_Scaled(t,j,false);
m_Regions[i].FillTower_Scaled(t,j,true,cGeom);
/*
if (et>0) {
std::cout<<"+++ "<<emEt<<" "<<hadEt<<" "
<<i<<" "<<j<<" "
<<std::endl<<"-- "
<<t.emEt()<<" "<<t.hadEt()<<" "
<<t.eta()<<" "<<t.phi()<<" "
<<std::endl;
}
*/
// Explicitely take care of integer boolean conversion
if (emEt > 3.0 && emEt < m_L1Config.noFGThreshold && (int)emFGV[i][j]!=0)
m_Regions[i].SetFGBit(j,true);
else
m_Regions[i].SetFGBit(j,false);
if (emEt > 3.0){
if (emEt < 60. && hadEt/emEt > m_L1Config.hOeThreshold)
m_Regions[i].SetHOEBit(j,true);
else
m_Regions[i].SetHOEBit(j,false);
}
else
m_Regions[i].SetHOEBit(j,false);
m_Regions[i].SetRegionBits(e);
}
}
}
}
| void FastL1GlobalAlgo::FillMET | ( | ) |
Definition at line 463 of file FastL1GlobalAlgo.cc.
References funct::cos(), ExpressReco_HICollisions_FallBack::et, fastmath::etaphi(), FastL1RegionMap::getRegionCenterEtaPhi(), i, l1extra::L1EtMissParticle::kMET, m_METs, m_Regions, m_RMap, phi, funct::sin(), and mathSSE::sqrt().
Referenced by FastL1CaloSim::produce().
{
m_METs.clear();
//double sum_e = 0.0;
double sum_et = 0.0;
double sum_ex = 0.0;
double sum_ey = 0.0;
//double sum_ez = 0.0;
for (int i=0; i<396; i++) {
std::pair<double, double> etaphi = m_RMap->getRegionCenterEtaPhi(i);
double phi = etaphi.second;
//double eta = etaphi.first;
//double theta = 2.*atan(exp(-eta));
double et = m_Regions[i].SumEt();
//double e = m_Regions[i].SumE();
//sum_et += et;
//sum_et += RCTEnergyTrunc(et,0.5,2048);
//sum_et += et;
//sum_ex += RCTEnergyTrunc(et*cos(phi),0.5,2048);
//sum_ey += RCTEnergyTrunc(et*sin(phi),0.5,2048);
sum_ex += et*cos(phi);
sum_ey += et*sin(phi);
//sum_ex += e*sin(theta)*cos(phi);
//sum_ey += e*sin(theta)*sin(phi);
//sum_e += e;
//sum_e += RCTEnergyTrunc(et/sin(theta),0.5,2048);
//sum_ez += RCTEnergyTrunc(et*cos(theta)/sin(theta),0.5,2048);
//sum_e += et/sin(theta);
//sum_ez += et*cos(theta)/sin(theta);
//sum_ez += e*cos(theta);
//sum_et += e*sin(theta);
}
//sum_et = RCTEnergyTrunc(std::sqrt(sum_ex*sum_ex + sum_ey*sum_ey),0.5,2048);
sum_et = std::sqrt(sum_ex*sum_ex + sum_ey*sum_ey);
//reco::Particle::LorentzVector rp4(-sum_ex,-sum_ey,-sum_ez,sum_e);
reco::Particle::LorentzVector rp4(-sum_ex,-sum_ey,0.,std::sqrt(sum_ex*sum_ex + sum_ey*sum_ey));
//edm::LogInfo("********** FastL1GlobalAlgo::FillMET()")<<rp4.mass()<<std::endl;
//m_MET = l1extra::L1EtMissParticle(rp4,sum_et,0.);
m_METs.push_back(l1extra::L1EtMissParticle(rp4,l1extra::L1EtMissParticle::kMET,sum_et));
}
| void FastL1GlobalAlgo::FillMET | ( | edm::Event const & | e | ) |
Definition at line 347 of file FastL1GlobalAlgo.cc.
References abs, funct::cos(), ExpressReco_HICollisions_FallBack::et, eta(), funct::exp(), edm::Event::getByLabel(), collect_tpl::input, l1extra::L1EtMissParticle::kMET, m_L1Config, m_METs, phi, RCTEnergyTrunc(), funct::sin(), mathSSE::sqrt(), theta(), L1Config::TowerEBScale, L1Config::TowerEBThreshold, L1Config::TowerEEScale, L1Config::TowerEEThreshold, L1Config::TowerHBScale, L1Config::TowerHBThreshold, L1Config::TowerHEScale, and L1Config::TowerInput.
{
m_METs.clear();
//std::vector< edm::Handle<CaloTowerCollection> > input;
//e.getManyByType(input);
edm::Handle<CaloTowerCollection> input;
e.getByLabel(m_L1Config.TowerInput,input);
double sum_hade = 0.0;
double sum_hadet = 0.0;
double sum_hadex = 0.0;
double sum_hadey = 0.0;
double sum_hadez = 0.0;
double sum_e = 0.0;
double sum_et = 0.0;
double sum_ex = 0.0;
double sum_ey = 0.0;
double sum_ez = 0.0;
//std::vector< edm::Handle<CaloTowerCollection> >::iterator i;
//for (i=input.begin(); i!=input.end(); i++) {
// const CaloTowerCollection& c=*(*i);
// for (CaloTowerCollection::const_iterator candidate=c.begin(); candidate!=c.end(); candidate++) {
for (CaloTowerCollection::const_iterator candidate=input->begin(); candidate!=input->end(); candidate++) {
//double eme = candidate->emEnergy();
//double hade = candidate->hadEnergy();
double eme = candidate->emEt();
double hade = candidate->hadEt();
double EThres = 0.;
double HThres = 0.;
double EBthres = m_L1Config.TowerEBThreshold;
double HBthres = m_L1Config.TowerHBThreshold;
double EEthres = m_L1Config.TowerEBThreshold;
double HEthres = m_L1Config.TowerEEThreshold;
//if(std::abs(candidate->eta())<1.479) {
if(std::abs(candidate->eta())<2.322) {
EThres = EBthres;
} else {
EThres = EEthres;
}
//if(std::abs(candidate->eta())<1.305) {
if(std::abs(candidate->eta())<2.322) {
HThres = HBthres;
} else {
HThres = HEthres;
}
// rescale energies
double emScale = 1.0;
double hadScale = 1.0;
if (std::abs(candidate->eta()>1.3050) && std::abs(candidate->eta())<3.0) {
hadScale = m_L1Config.TowerHEScale;
emScale = m_L1Config.TowerEEScale;
}
if (std::abs(candidate->eta()<1.3050)) {
hadScale = m_L1Config.TowerHBScale;
emScale = m_L1Config.TowerEBScale;
}
eme *= emScale;
hade *= hadScale;
if (eme>=EThres || hade>=HThres) {
double phi = candidate->phi();
double eta = candidate->eta();
//double et = candidate->et();
//double e = candidate->energy();
double theta = 2.*atan(exp(-eta));
double et = 0.;
double e = 0.;
double had_et = 0.;
double had_e = 0.;
if (eme>=EThres) {
et += candidate->emEt();
e += candidate->emEnergy();
}
if (hade>=HThres) {
et += candidate->hadEt();
e += candidate->hadEnergy();
had_et += candidate->hadEt();
had_e += candidate->hadEnergy();
}
//sum_et += et;
sum_et += RCTEnergyTrunc(et,1.0,1024);
sum_ex += et*cos(phi);
sum_ey += et*sin(phi);
//sum_ex += e*sin(theta)*cos(phi);
//sum_ey += e*sin(theta)*sin(phi);
//sum_e += e;
sum_e += et/sin(theta);
sum_ez += et*cos(theta)/sin(theta);
sum_hadet += had_et;
sum_hadex += had_et*cos(phi);
sum_hadey += had_et*sin(phi);
//sum_hadex += had_e*sin(theta)*cos(phi);
//sum_hadey += had_e*sin(theta)*sin(phi);
//sum_hade += had_e;
sum_hade += had_et/sin(theta);
sum_hadez += had_et*cos(theta)/sin(theta);
}
}
//}
reco::Particle::LorentzVector rp4(-sum_ex,-sum_ey,0.,std::sqrt(sum_ex*sum_ex + sum_ey*sum_ey));
//m_MET = l1extra::L1EtMissParticle(rp4,sum_et,0.);
m_METs.push_back(l1extra::L1EtMissParticle(rp4,l1extra::L1EtMissParticle::kMET,sum_et));
}
| void FastL1GlobalAlgo::findJets | ( | ) | [private] |
Definition at line 155 of file FastL1GlobalAlgo.cc.
References addJet(), eta(), FastL1RegionMap::getRegionCenterEtaPhi(), i, isMaxEtRgn_Window33(), isTauJet(), L1Config::JetSeedEtThreshold, m_CenJets, m_DoBitInfo, m_ForJets, m_GctIso, m_L1Config, m_Regions, m_RMap, m_TauJets, L1TEmulatorMonitor_cff::p, phi, and TauIsolation().
Referenced by FillJets().
{
m_TauJets.clear();
m_CenJets.clear();
m_ForJets.clear();
for (int i=0; i<396; i++) {
// barrel/endcap part only:
//if ((i%22)>3 && (i%22)<18) {
std::pair<double, double> p = m_RMap->getRegionCenterEtaPhi(i);
//if (m_Regions.at(i).SumEt()>0.)
// std::cout<<"******** TEST "<<i<<": "<<m_Regions.at(i).SumEt()<<std::endl;
double eta = p.first;
double phi = p.second;
if (m_DoBitInfo){
m_Regions[i].BitInfo.setEta(eta);
m_Regions[i].BitInfo.setPhi(phi);
}
if (m_Regions.at(i).SumEt()>m_L1Config.JetSeedEtThreshold) {
if (isMaxEtRgn_Window33(i)) {
if (m_GctIso == true) {
if (TauIsolation(i) && (i%22)>3 && (i%22)<18 ) {
addJet(i,true);
} else {
addJet(i,false);
}
} else {
if (isTauJet(i) && (i%22)>3 && (i%22)<18 ) {
addJet(i,true);
} else {
addJet(i,false);
}
}
if (m_DoBitInfo) m_Regions[i].BitInfo.setMaxEt(true);
}
else {
if (m_DoBitInfo) m_Regions[i].BitInfo.setMaxEt(false);
}
if (m_DoBitInfo) m_Regions[i].BitInfo.setSumEtBelowThres (false);
} else {
if (m_DoBitInfo) m_Regions[i].BitInfo.setSumEtBelowThres (true);
}
//}
}
}
| FastL1BitInfoCollection FastL1GlobalAlgo::getBitInfos | ( | ) | [inline] |
Definition at line 84 of file FastL1GlobalAlgo.h.
References m_BitInfos.
Referenced by FastL1CaloSim::produce().
{ return m_BitInfos; }
| std::vector<FastL1Region> FastL1GlobalAlgo::GetCaloRegions | ( | ) | [inline] |
| l1extra::L1JetParticleCollection FastL1GlobalAlgo::getCenJets | ( | ) | const [inline] |
Definition at line 80 of file FastL1GlobalAlgo.h.
References m_CenJets.
Referenced by FastL1CaloSim::produce().
{ return m_CenJets; }
| l1extra::L1EmParticleCollection FastL1GlobalAlgo::getEgammas | ( | ) | const [inline] |
Definition at line 82 of file FastL1GlobalAlgo.h.
References m_Egammas.
Referenced by FastL1CaloSim::produce().
{ return m_Egammas; }
| l1extra::L1JetParticleCollection FastL1GlobalAlgo::getForJets | ( | ) | const [inline] |
Definition at line 81 of file FastL1GlobalAlgo.h.
References m_ForJets.
Referenced by FastL1CaloSim::produce().
{ return m_ForJets; }
| l1extra::L1EmParticleCollection FastL1GlobalAlgo::getisoEgammas | ( | ) | const [inline] |
Definition at line 83 of file FastL1GlobalAlgo.h.
References m_isoEgammas.
Referenced by FastL1CaloSim::produce().
{ return m_isoEgammas; }
| l1extra::L1EtMissParticleCollection FastL1GlobalAlgo::getMET | ( | ) | const [inline] |
Definition at line 78 of file FastL1GlobalAlgo.h.
References m_METs.
Referenced by FastL1CaloSim::produce().
{ return m_METs; }
| l1extra::L1JetParticleCollection FastL1GlobalAlgo::getTauJets | ( | ) | const [inline] |
Definition at line 79 of file FastL1GlobalAlgo.h.
References m_TauJets.
Referenced by FastL1CaloSim::produce().
{ return m_TauJets; }
| bool FastL1GlobalAlgo::greaterEt | ( | const reco::Candidate & | a, |
| const reco::Candidate & | b | ||
| ) | [private] |
Referenced by addJet(), FillEgammas(), and FillEgammasTP().
| double FastL1GlobalAlgo::hcaletValue | ( | const int | ieta, |
| const int | compET | ||
| ) | [private] |
Definition at line 1366 of file FastL1GlobalAlgo.cc.
References m_hcaluncomp.
Referenced by FillL1RegionsTP().
{
double etvalue = m_hcaluncomp[ieta][compET];//*cos(eta_ave);
return etvalue;
}
| void FastL1GlobalAlgo::InitL1Regions | ( | ) |
Definition at line 511 of file FastL1GlobalAlgo.cc.
References trackerHits::c, FastL1RegionMap::getRegionEtaPhiIndex(), i, m_DoBitInfo, m_L1Config, m_Regions, m_RMap, and L1TEmulatorMonitor_cff::p.
Referenced by FillL1Regions(), and FillL1RegionsTP().
{
m_Regions.clear();
m_Regions = std::vector<FastL1Region>(396); // ieta: 1-22, iphi: 1-18
// init regions
for (int i=0; i<396; i++) {
m_Regions[i].SetParameters(m_L1Config);
m_Regions[i].SetDoBitInfo(m_DoBitInfo);
std::pair<int, int> p = m_RMap->getRegionEtaPhiIndex(i);
m_Regions[i].SetEtaPhiIndex(p.first,p.second,i);
CaloTower c;
for (int twrid=0; twrid<16; twrid++) {
m_Regions[i].FillTowerZero(c,twrid);
}
}
}
| int FastL1GlobalAlgo::isEMCand | ( | CaloTowerDetId | cid, |
| l1extra::L1EmParticle * | p, | ||
| const edm::Event & | e | ||
| ) | [private] |
Definition at line 932 of file FastL1GlobalAlgo.cc.
References trackerHits::c, funct::cos(), L1Config::EMLSB, L1Config::EMSeedEnThreshold, funct::exp(), GCTEnergyTrunc(), FastL1RegionMap::getRegionCenterEtaPhi(), FastL1RegionMap::getRegionIndex(), FastL1RegionMap::getRegionTowerIndex(), FastL1RegionMap::GetTowerEastEtaPhi(), FastL1RegionMap::GetTowerNEEtaPhi(), FastL1RegionMap::GetTowerNorthEtaPhi(), FastL1RegionMap::GetTowerNWEtaPhi(), FastL1RegionMap::GetTowerSEEtaPhi(), FastL1RegionMap::GetTowerSouthEtaPhi(), FastL1RegionMap::GetTowerSWEtaPhi(), FastL1RegionMap::GetTowerWestEtaPhi(), CaloTowerDetId::ieta(), CaloTowerDetId::iphi(), m_L1Config, m_Regions, m_RMap, max(), python::Vispa::Plugins::EdmBrowser::EdmDataAccessor::ne(), L1Config::QuietRegionThreshold, funct::sin(), and L1Config::TowerEBThreshold.
Referenced by FillEgammas(), and FillEgammasTP().
{
// center tower
int crgn = m_RMap->getRegionIndex(cid);
int ctwr = m_RMap->getRegionTowerIndex(cid);
// crystals only in barrel/endcap part
if ((crgn%22)<4 || (crgn%22)>17) return 0;
if (crgn>395 || crgn < 0 || ctwr > 15 || ctwr < 0) return 0;
CaloTowerCollection c = m_Regions.at(crgn).GetCaloTowers();
double cenEt = c[ctwr].et();
//double cenEt = RCTEnergyTrunc(c[ctwr].et(),0.5,64);
//double cenE = c[ctwr].emEnergy();
// Using region position rather than tower position
std::pair<double, double> crpos = m_RMap->getRegionCenterEtaPhi(crgn);
//double cenEta = c[ctwr].eta();
//double cenPhi = c[ctwr].phi();
double cenEta = crpos.first;
double cenPhi = crpos.second;
double cenFGbit = m_Regions.at(crgn).GetFGBit(ctwr);
double cenHOEbit = m_Regions.at(crgn).GetHOEBit(ctwr);
if (cenEt<m_L1Config.TowerEBThreshold) return 0;
// check fine grain bit
if (cenFGbit) return 0;
// check H/E bit
if (cenHOEbit) return 0;
// check neighbours
std::pair<int, int> no = m_RMap->GetTowerNorthEtaPhi(cid.ieta(),cid.iphi());
std::pair<int, int> so = m_RMap->GetTowerSouthEtaPhi(cid.ieta(),cid.iphi());
std::pair<int, int> we = m_RMap->GetTowerWestEtaPhi(cid.ieta(),cid.iphi());
std::pair<int, int> ea = m_RMap->GetTowerEastEtaPhi(cid.ieta(),cid.iphi());
std::pair<int, int> nw = m_RMap->GetTowerNWEtaPhi(cid.ieta(),cid.iphi());
std::pair<int, int> ne = m_RMap->GetTowerNEEtaPhi(cid.ieta(),cid.iphi());
std::pair<int, int> sw = m_RMap->GetTowerSWEtaPhi(cid.ieta(),cid.iphi());
std::pair<int, int> se = m_RMap->GetTowerSEEtaPhi(cid.ieta(),cid.iphi());
if (no.first>29 || no.first<-29 || no.second>72 || no.second<0) return 0;
if (so.first>29 || so.first<-29 || so.second>72 || so.second<0) return 0;
if (we.first>29 || we.first<-29 || we.second>72 || we.second<0) return 0;
if (ea.first>29 || ea.first<-29 || ea.second>72 || ea.second<0) return 0;
if (nw.first>29 || nw.first<-29 || nw.second>72 || nw.second<0) return 0;
if (ne.first>29 || ne.first<-29 || ne.second>72 || ne.second<0) return 0;
if (sw.first>29 || sw.first<-29 || sw.second>72 || sw.second<0) return 0;
if (se.first>29 || se.first<-29 || se.second>72 || se.second<0) return 0;
int notwr = m_RMap->getRegionTowerIndex(no);
int norgn = m_RMap->getRegionIndex(no.first,no.second);
int sotwr = m_RMap->getRegionTowerIndex(so);
int sorgn = m_RMap->getRegionIndex(so.first,so.second);
int wetwr = m_RMap->getRegionTowerIndex(we);
int wergn = m_RMap->getRegionIndex(we.first,we.second);
int eatwr = m_RMap->getRegionTowerIndex(ea);
int eargn = m_RMap->getRegionIndex(ea.first,ea.second);
int setwr = m_RMap->getRegionTowerIndex(se);
int sergn = m_RMap->getRegionIndex(se.first,sw.second);
int swtwr = m_RMap->getRegionTowerIndex(sw);
int swrgn = m_RMap->getRegionIndex(sw.first,sw.second);
int netwr = m_RMap->getRegionTowerIndex(ne);
int nergn = m_RMap->getRegionIndex(ne.first,ne.second);
int nwtwr = m_RMap->getRegionTowerIndex(nw);
int nwrgn = m_RMap->getRegionIndex(nw.first,nw.second);
//
if (norgn>395 || norgn < 0 || notwr > 15 || notwr < 0) return 0;
c = m_Regions[norgn].GetCaloTowers();
double noEt = c[notwr].et();
//double noEt = RCTEnergyTrunc(c[notwr].et(),0.5,64);
//double noE = c[notwr].emEnergy();
// check fine grain bit
bool noFGbit = m_Regions[norgn].GetFGBit(notwr);
// check H/E bit
bool noHOEbit = m_Regions[norgn].GetHOEBit(notwr);
//
if (sorgn>395 || sorgn < 0 || sotwr > 15 || sotwr < 0) return 0;
c = m_Regions[sorgn].GetCaloTowers();
double soEt = c[sotwr].et();
//double soEt = RCTEnergyTrunc(c[sotwr].et(),0.5,64);
//double soE = c[sotwr].emEnergy();
// check fine grain bit
bool soFGbit = m_Regions[sorgn].GetFGBit(sotwr);
// check H/E bit
bool soHOEbit = m_Regions[sorgn].GetHOEBit(sotwr);
//
if (wergn>395 || wergn < 0 || wetwr > 15 || wetwr < 0) return 0;
c = m_Regions[wergn].GetCaloTowers();
double weEt = c[wetwr].et();
//double weEt = RCTEnergyTrunc(c[wetwr].et(),0.5,64);
//double weE = c[wetwr].emEnergy();
// check fine grain bit
bool weFGbit = m_Regions[wergn].GetFGBit(wetwr);
// check H/E bit
bool weHOEbit = m_Regions[wergn].GetHOEBit(wetwr);
//
if (eargn>395 || eargn < 0 || eatwr > 15 || eatwr < 0) return 0;
c = m_Regions[eargn].GetCaloTowers();
double eaEt = c[eatwr].et();
//double eaEt = RCTEnergyTrunc(c[eatwr].et(),0.5,64);
//double eaE = c[eatwr].emEnergy();
// check fine grain bit
bool eaFGbit = m_Regions[eargn].GetFGBit(eatwr);
// check H/E bit
bool eaHOEbit = m_Regions[eargn].GetHOEBit(eatwr);
//
if (nwrgn>395 || nwrgn < 0 || nwtwr > 15 || nwtwr < 0) return 0;
c = m_Regions[nwrgn].GetCaloTowers();
double nwEt = c[nwtwr].et();
//double nwEt = RCTEnergyTrunc(c[nwtwr].et(),0.5,64);
//double nwE = c[nwtwr].emEnergy();
// check fine grain bit
bool nwFGbit = m_Regions[nwrgn].GetFGBit(nwtwr);
// check H/E bit
bool nwHOEbit = m_Regions[nwrgn].GetHOEBit(nwtwr);
//
if (nergn>395 || nergn < 0 || netwr > 15 || netwr < 0) return 0;
c = m_Regions[nergn].GetCaloTowers();
double neEt = c[netwr].et();
//double neEt = RCTEnergyTrunc(c[netwr].et(),0.5,64);
//double neE = c[netwr].emEnergy();
// check fine grain bit
bool neFGbit = m_Regions[nergn].GetFGBit(netwr);
// check H/E bit
bool neHOEbit = m_Regions[nergn].GetHOEBit(netwr);
//
if (swrgn>395 || swrgn < 0 || swtwr > 15 || swtwr < 0) return 0;
c = m_Regions[swrgn].GetCaloTowers();
double swEt = c[swtwr].et();
//double swEt = RCTEnergyTrunc(c[swtwr].et(),0.5,64);
//double swE = c[swtwr].emEnergy();
// check fine grain bit
bool swFGbit = m_Regions[swrgn].GetFGBit(swtwr);
// check H/E bit
bool swHOEbit = m_Regions[swrgn].GetHOEBit(swtwr);
//
if (sergn>395 || sergn < 0 || setwr > 15 || setwr < 0) return 0;
c = m_Regions[sergn].GetCaloTowers();
double seEt = c[setwr].et();
//double seEt = RCTEnergyTrunc(c[setwr].et(),0.5,64);
//double seE = c[setwr].emEnergy();
// check fine grain bit
bool seFGbit = m_Regions[sergn].GetFGBit(setwr);
// check H/E bit
bool seHOEbit = m_Regions[sergn].GetHOEBit(setwr);
// check if highest et tower
bool isHit = false;
if ( cenEt > noEt && cenEt >= soEt && cenEt > weEt &&
cenEt >= eaEt && cenEt > nwEt && cenEt > neEt &&
cenEt >= swEt && cenEt >= seEt ) isHit = true;
else
return 0;
// find highest neighbour
double hitEt = cenEt;
//double hitE = cenE;
double maxEt = std::max(noEt,std::max(soEt,std::max(weEt,eaEt)));
//double maxE = std::max(noE,std::max(soE,std::max(weE,eaE)));
// check 2 tower Et
float emEtThres = m_L1Config.EMSeedEnThreshold;
// at this point candidate is at least non-iso Egamma
//double eme = (hitE+maxE);
//double eme = (hitE+maxE);
double emet = (hitEt+maxEt);
/*
if (m_L1Config.DoEMCorr) {
emet = GCTEnergyTrunc(corrEmEt(emet,cenEta),m_L1Config.EMLSB, true);
//emet = GCTEnergyTrunc(emet,m_L1Config.EMLSB, true);
} else {
emet = GCTEnergyTrunc(emet,m_L1Config.EMLSB, true);
}
*/
emet = GCTEnergyTrunc(emet,m_L1Config.EMLSB, true);
if ((emet)<emEtThres) return 0;
double emtheta = 2.*atan(exp(-cenEta));
//double emet = eme*sin(emtheta);
double emex = emet*cos(cenPhi);
double emey = emet*sin(cenPhi);
//double emex = eme*sin(emtheta)*cos(cenPhi);
//double emey = eme*sin(emtheta)*sin(cenPhi);
double eme = emex/sin(emtheta)/cos(cenPhi);
double emez = eme*cos(emtheta);
reco::Particle::LorentzVector rp4(emex,emey,emez,eme);
//reco::Particle::Point rp3(0.,0.,0.);
//reco::Particle::Charge q(0);
//*ph = reco::Photon(q,rp4,rp3);
*ph = l1extra::L1EmParticle(rp4);
//ph = l1extra::L1EmParticle(rp4);
//std::cout<<"EM eme : "<<eme<<std::endl;
//std::cout<<"EM rp4.et(): "<<rp4.Et()<<std::endl;
//std::cout<<"EM ph->et() : "<<ph->et()<<std::endl;
//if (emet>0.) {
// std::cout << "em region et, eta, phi: "<< emet<<" "<< cenEta<<" "<< cenPhi<<" " << std::endl;
// std::cout << "em lv et, eta, phi: "<< ph->et()<<" "<< ph->eta()<<" "<< ph->phi()<<" " << std::endl;
//}
// check isolation FG bits
if (noFGbit || soFGbit || weFGbit || eaFGbit ||
nwFGbit || neFGbit || swFGbit || seFGbit ||
noHOEbit || soHOEbit || weHOEbit || eaHOEbit ||
nwHOEbit || neHOEbit || swHOEbit || seHOEbit)
return 1;
// check isolation corners
//double corThres = 0.4;
//double quietThres = m_L1Config.EMNoiseLevel;
double quietThres = m_L1Config.QuietRegionThreshold;
bool isoVeto1 = false,isoVeto2 = false,isoVeto3 = false,isoVeto4 = false;
if (swEt>quietThres || weEt>quietThres || nwEt>quietThres || noEt>quietThres || neEt>quietThres ) {
//if ((swEt + weEt + nwEt + noEt + neEt)/cenEt > corThres)
isoVeto1 = true;
}
if (neEt>quietThres || eaEt>quietThres || seEt>quietThres || soEt>quietThres || swEt>quietThres ) {
//if ((neEt + eaEt + seEt + soEt + swEt)/cenEt > corThres)
isoVeto2 = true;
}
if (nwEt>quietThres || noEt>quietThres || neEt>quietThres || eaEt>quietThres || seEt>quietThres ) {
//if ((nwEt + noEt + neEt + eaEt + seEt)/cenEt > corThres)
isoVeto3 = true;
}
if (seEt>quietThres || soEt>quietThres || swEt>quietThres || weEt>quietThres || nwEt>quietThres ) {
//if ((seEt + soEt + swEt + weEt + nwEt)/cenEt > corThres)
isoVeto4 = true;
}
if (isoVeto1 && isoVeto2 && isoVeto3 && isoVeto4)
return 1;
return 2;
}
| bool FastL1GlobalAlgo::isMaxEtRgn_Window33 | ( | int | rgnid | ) | [private] |
Definition at line 1186 of file FastL1GlobalAlgo.cc.
References m_Regions, and evf::utils::sid.
Referenced by findJets().
{
int nwid = m_Regions.at(crgn).GetNWId();
int nid = m_Regions.at(crgn).GetNorthId();
int neid = m_Regions.at(crgn).GetNEId();
int wid = m_Regions.at(crgn).GetWestId();
int eid = m_Regions.at(crgn).GetEastId();
int swid = m_Regions.at(crgn).GetSWId();
int sid = m_Regions.at(crgn).GetSouthId();
int seid = m_Regions.at(crgn).GetSEId();
//Use 3x2 window at eta borders!
// east border:
if ((crgn%22)==21) {
if (nwid==999 || nid==999 || swid==999 || sid==999 || wid==999 ) {
return false;
}
double cenet = m_Regions.at(crgn).SumEt();
double nwet = m_Regions[nwid].SumEt();
double noet = m_Regions[nid].SumEt();
double weet = m_Regions[wid].SumEt();
double swet = m_Regions[swid].SumEt();
double soet = m_Regions[sid].SumEt();
if ( cenet > nwet && cenet > noet &&
cenet >= weet && cenet >= soet &&
cenet >= swet )
{
double cene = m_Regions.at(crgn).SumE();
double nwe = m_Regions[nwid].SumE();
double noe = m_Regions[nid].SumE();
double wee = m_Regions[wid].SumE();
double swe = m_Regions[swid].SumE();
double soe = m_Regions[sid].SumE();
// if region is central: jet energy is sum of 3x3 region
// surrounding the central region
double jE = cene + nwe + noe + wee + swe + soe;
double jEt = cenet + nwet + noet + weet + swet + soet;
m_Regions.at(crgn).SetJetE(jE);
m_Regions.at(crgn).SetJetEt(jEt);
m_Regions.at(crgn).SetJetE3x3(cene);
m_Regions.at(crgn).SetJetEt3x3(cenet);
return true;
} else { return false; }
}
// west border:
if ((crgn%22)==0) {
if (neid==999 || nid==999 || seid==999 || sid==999 || eid==999 ) {
return false;
}
double cenet = m_Regions.at(crgn).SumEt();
double neet = m_Regions[neid].SumEt();
double noet = m_Regions[nid].SumEt();
double eaet = m_Regions[eid].SumEt();
double seet = m_Regions[seid].SumEt();
double soet = m_Regions[sid].SumEt();
if ( cenet > neet && cenet > noet &&
cenet >= eaet && cenet >= soet &&
cenet >= seet )
{
double cene = m_Regions.at(crgn).SumE();
double nee = m_Regions[neid].SumE();
double noe = m_Regions[nid].SumE();
double eae = m_Regions[eid].SumE();
double see = m_Regions[seid].SumE();
double soe = m_Regions[sid].SumE();
// if region is central: jet energy is sum of 3x3 region
// surrounding the central region
double jE = cene + nee + noe + eae + see + soe;
double jEt = cenet + neet + noet + eaet + seet + soet;
m_Regions.at(crgn).SetJetE(jE);
m_Regions.at(crgn).SetJetEt(jEt);
m_Regions.at(crgn).SetJetE3x3(cene);
m_Regions.at(crgn).SetJetEt3x3(cenet);
return true;
} else { return false; }
}
if (nwid==999 || neid==999 || nid==999 || swid==999 || seid==999 || sid==999 || wid==999 ||
eid==999 ) {
//std::cerr << "FastL1GlobalAlgo::isMaxEtRgn_Window33(): RegionId out of bounds: " << std::endl
// << nwid << " " << nid << " " << neid << " " << std::endl
// << wid << " " << crgn << " " << eid << " " << std::endl
// << swid << " " << sid << " " << seid << " " << std::endl;
return false;
}
double cenet = m_Regions.at(crgn).SumEt();
double nwet = m_Regions[nwid].SumEt();
double noet = m_Regions[nid].SumEt();
double neet = m_Regions[neid].SumEt();
double weet = m_Regions[wid].SumEt();
double eaet = m_Regions[eid].SumEt();
double swet = m_Regions[swid].SumEt();
double soet = m_Regions[sid].SumEt();
double seet = m_Regions[seid].SumEt();
if ( cenet > nwet && cenet > noet &&
cenet > neet && cenet >= eaet &&
cenet > weet && cenet >= soet &&
cenet >= swet && cenet >= seet )
{
double cene = m_Regions.at(crgn).SumE();
double nwe = m_Regions[nwid].SumE();
double noe = m_Regions[nid].SumE();
double nee = m_Regions[neid].SumE();
double wee = m_Regions[wid].SumE();
double eae = m_Regions[eid].SumE();
double swe = m_Regions[swid].SumE();
double soe = m_Regions[sid].SumE();
double see = m_Regions[seid].SumE();
// if region is central: jet energy is sum of 3x3 region
// surrounding the central region
double jE = cene + nwe + noe + nee + wee + eae + swe + soe + see;
double jEt = cenet + nwet + noet + neet + weet + eaet + swet + soet + seet;
m_Regions.at(crgn).SetJetE(jE);
m_Regions.at(crgn).SetJetEt(jEt);
m_Regions.at(crgn).SetJetE3x3(cene);
m_Regions.at(crgn).SetJetEt3x3(cenet);
return true;
} else { return false; }
}
| bool FastL1GlobalAlgo::isTauJet | ( | int | rgnid | ) | [private] |
Definition at line 843 of file FastL1GlobalAlgo.cc.
References m_DoBitInfo, m_Regions, and evf::utils::sid.
Referenced by findJets().
{
// Barrel and Endcap only
if ((cRgn%22)<4 || (cRgn%22)>17)
return false;
int shower_shape = 0;
int et_isolation = 0;
if (m_Regions[cRgn].GetTauBit()) {
shower_shape = 1;
//if (m_DoBitInfo) m_Regions[cRgn].BitInfo.setTauVeto(true);
} else {
shower_shape = 0;
//if (m_DoBitInfo) m_Regions[cRgn].BitInfo.setTauVeto(false);
}
int nwid = m_Regions[cRgn].GetNWId();
int nid = m_Regions[cRgn].GetNorthId();
int neid = m_Regions[cRgn].GetNEId();
int wid = m_Regions[cRgn].GetWestId();
int eid = m_Regions[cRgn].GetEastId();
int swid = m_Regions[cRgn].GetSWId();
int sid = m_Regions[cRgn].GetSouthId();
int seid = m_Regions[cRgn].GetSEId();
//Use 3x2 window at eta borders!
// west border:
if ((cRgn%22)==4) {
if (m_Regions[nid].GetTauBit() ||
m_Regions[neid].GetTauBit() ||
m_Regions[eid].GetTauBit() ||
m_Regions[seid].GetTauBit() ||
m_Regions[sid].GetTauBit() ) {
et_isolation = 1;
}
}
// east border:
if ((cRgn%22)==17) {
if (m_Regions[nid].GetTauBit() ||
m_Regions[nwid].GetTauBit() ||
m_Regions[wid].GetTauBit() ||
m_Regions[swid].GetTauBit() ||
m_Regions[sid].GetTauBit() ) {
et_isolation = 1;
}
}
if ( (cRgn%22)>4 && (cRgn%22)<17){ // non-border
if (nwid==999 || neid==999 || nid==999 || swid==999 || seid==999 || sid==999 || wid==999 ||
eid==999 ) {
return false;
}
if (m_Regions[nwid].GetTauBit() ||
m_Regions[nid].GetTauBit() ||
m_Regions[neid].GetTauBit() ||
m_Regions[wid].GetTauBit() ||
m_Regions[eid].GetTauBit() ||
m_Regions[swid].GetTauBit() ||
m_Regions[seid].GetTauBit() ||
m_Regions[sid].GetTauBit() ) {
et_isolation = 1;
}
} // non-border
if (m_DoBitInfo) {
if (et_isolation == 1) {
//std::cout<<"*********************************"<<std::endl;
//std::cout<<"Rgn iso veto: "<<et_isolation<<std::endl;
m_Regions[cRgn].BitInfo.setIsolationVeto(true);
} else {
m_Regions[cRgn].BitInfo.setIsolationVeto(false);
}
// if (shower_shape == 1) {
//std::cout<<"###############################"<<std::endl;
//std::cout<<"Rgn shower veto: "<<shower_shape<<std::endl;
//}
}
if (et_isolation == 1 || shower_shape == 1) return false;
else return true;
}
| bool FastL1GlobalAlgo::TauIsolation | ( | int | rgnid | ) | [private] |
Definition at line 1373 of file FastL1GlobalAlgo.cc.
References m_DoBitInfo, m_IsolationEt, m_Regions, and evf::utils::sid.
Referenced by findJets().
{
if ((cRgn%22)<4 || (cRgn%22)>17) return false;
double iso_threshold = m_IsolationEt; // arbitrarily set
int shower_shape = 0;
int et_isolation = 0;
unsigned int iso_count = 0;
int nwid = m_Regions[cRgn].GetNWId();
int nid = m_Regions[cRgn].GetNorthId();
int neid = m_Regions[cRgn].GetNEId();
int wid = m_Regions[cRgn].GetWestId();
int eid = m_Regions[cRgn].GetEastId();
int swid = m_Regions[cRgn].GetSWId();
int sid = m_Regions[cRgn].GetSouthId();
int seid = m_Regions[cRgn].GetSEId();
if (m_Regions[cRgn].GetTauBit()) { // check center
shower_shape = 1;
//if (m_DoBitInfo) m_Regions[cRgn].BitInfo.setTauVeto(true);
}
if((cRgn%22)==4 || (cRgn%22)==17 ) {
// west border
if ((cRgn%22)==4) {
if( m_Regions[neid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[neid].GetTauBit()) iso_count++;
}
if( m_Regions[nid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[nid].GetTauBit()) iso_count++;
}
if( m_Regions[eid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[eid].GetTauBit()) iso_count++;
}
if( m_Regions[seid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[seid].GetTauBit()) iso_count++;
}
if( m_Regions[sid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[sid].GetTauBit()) iso_count++;
}
} // west bd
// east border:
if ((cRgn%22)==17) {
if( m_Regions[nwid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[nwid].GetTauBit()) iso_count++;
}
if( m_Regions[nid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[nid].GetTauBit()) iso_count++;
}
if( m_Regions[wid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[wid].GetTauBit()) iso_count++;
}
if( m_Regions[swid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[swid].GetTauBit()) iso_count++;
}
if( m_Regions[sid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[sid].GetTauBit()) iso_count++;
}
} // east bd
}
if ( (cRgn%22)>4 && (cRgn%22)<17){ // non-border
if (nwid==999 || neid==999 || nid==999 || swid==999 || seid==999 || sid==999 || wid==999 ||
eid==999 ) {
return false;
}
if( m_Regions[neid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[neid].GetTauBit()) iso_count++;
}
if( m_Regions[nid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[nid].GetTauBit()) iso_count++;
}
if( m_Regions[eid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[eid].GetTauBit()) iso_count++;
}
if( m_Regions[seid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[seid].GetTauBit()) iso_count++;
}
if( m_Regions[sid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[sid].GetTauBit()) iso_count++;
}
if( m_Regions[nwid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[nwid].GetTauBit()) iso_count++;
}
if( m_Regions[wid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[wid].GetTauBit()) iso_count++;
}
if( m_Regions[swid].SumEt() > iso_threshold){
iso_count ++;
if (m_Regions[swid].GetTauBit()) iso_count++;
}
}// non-border
if (iso_count >= 2 ){
et_isolation = 1;
}
else {
et_isolation = 0;
}
if (m_DoBitInfo){
if (et_isolation == 1)
m_Regions[cRgn].BitInfo.setIsolationVeto (true);
else
m_Regions[cRgn].BitInfo.setIsolationVeto (false);
}
if (et_isolation == 1 || shower_shape == 1) return false;
else return true;
}//
Definition at line 125 of file FastL1GlobalAlgo.h.
Referenced by FillBitInfos(), and getBitInfos().
Definition at line 121 of file FastL1GlobalAlgo.h.
Referenced by addJet(), findJets(), and getCenJets().
bool FastL1GlobalAlgo::m_DoBitInfo [private] |
Definition at line 129 of file FastL1GlobalAlgo.h.
Referenced by addJet(), FastL1GlobalAlgo(), FillBitInfos(), findJets(), InitL1Regions(), isTauJet(), and TauIsolation().
Definition at line 123 of file FastL1GlobalAlgo.h.
Referenced by FillEgammas(), FillEgammasTP(), and getEgammas().
Definition at line 122 of file FastL1GlobalAlgo.h.
Referenced by addJet(), findJets(), and getForJets().
bool FastL1GlobalAlgo::m_GctIso [private] |
Definition at line 131 of file FastL1GlobalAlgo.h.
Referenced by FastL1GlobalAlgo(), and findJets().
double FastL1GlobalAlgo::m_hcaluncomp[33][256] [private] |
Definition at line 135 of file FastL1GlobalAlgo.h.
Referenced by FastL1GlobalAlgo(), and hcaletValue().
Definition at line 124 of file FastL1GlobalAlgo.h.
Referenced by FillEgammas(), FillEgammasTP(), and getisoEgammas().
double FastL1GlobalAlgo::m_IsolationEt [private] |
Definition at line 132 of file FastL1GlobalAlgo.h.
Referenced by FastL1GlobalAlgo(), and TauIsolation().
L1Config FastL1GlobalAlgo::m_L1Config [private] |
Definition at line 134 of file FastL1GlobalAlgo.h.
Referenced by addJet(), CaloTowersDump(), FastL1GlobalAlgo(), FillEgammas(), FillL1Regions(), FillL1RegionsTP(), FillMET(), findJets(), InitL1Regions(), and isEMCand().
Definition at line 119 of file FastL1GlobalAlgo.h.
std::vector<FastL1Region> FastL1GlobalAlgo::m_Regions [private] |
Definition at line 127 of file FastL1GlobalAlgo.h.
Referenced by addJet(), FastL1GlobalAlgo(), FillBitInfos(), FillEgammasTP(), FillL1Regions(), FillL1RegionsTP(), FillMET(), findJets(), GetCaloRegions(), InitL1Regions(), isEMCand(), isMaxEtRgn_Window33(), isTauJet(), and TauIsolation().
FastL1RegionMap* FastL1GlobalAlgo::m_RMap [private] |
Definition at line 128 of file FastL1GlobalAlgo.h.
Referenced by addJet(), checkMapping(), FastL1GlobalAlgo(), FillL1Regions(), FillL1RegionsTP(), FillMET(), findJets(), InitL1Regions(), and isEMCand().
Definition at line 120 of file FastL1GlobalAlgo.h.
Referenced by addJet(), findJets(), and getTauJets().
1.7.3