FastL1Region Class Reference

#include <FastL1Region.h>

Public Member Functions
double	CalcSumE ()
double	CalcSumEmE ()
double	CalcSumEmEt ()
double	CalcSumEt ()
double	CalcSumHadE ()
double	CalcSumHadEt ()
void	Dump ()
	FastL1Region ()
void	FillEMCrystals (const CaloTowerConstituentsMap *theTowerConstituentsMap, const CaloTopology *calotopo, const CaloGeometry *cGeom, const EcalRecHitCollection *ec0, const EcalRecHitCollection *ec1, FastL1RegionMap *m_RMap)
void	FillTower (const CaloTower &t, int &tid, edm::ESHandle< CaloGeometry > &cGeom)
void	FillTower_Scaled (const CaloTower &t, int &tid, bool doRCTTrunc, edm::ESHandle< CaloGeometry > &cGeom)
void	FillTowerZero (const CaloTower &t, int &tid)
FastL1BitInfo	getBitInfo ()
CaloTowerCollection	GetCaloTowers ()
int	GetEastId ()
double	GetEMCrystalEnergy (int itwr, int icell)
bool	GetFGBit (int i)
bool	GetHCFGBit (int i)
bool	GetHOEBit (int i)
int	GetId ()
int	GetiEta ()
int	GetiPhi ()
double	GetJetE ()
double	GetJetE3x3 ()
double	GetJetEt ()
double	GetJetEt3x3 ()
bool	GetMIPBit ()
int	GetNEId ()
int	GetNorthId ()
int	GetNWId ()
bool	GetQuietBit ()
std::pair< double, double >	getRegionCenterEtaPhi (const edm::EventSetup &c)
int	GetSEId ()
int	GetSouthId ()
int	GetSWId ()
bool	GetTauBit ()
std::pair< int, int >	GetTowerNorthEtaPhi (int ieta, int iphi)
int	GetWestId ()
int	HighestEmEtTowerID ()
int	HighestEtTowerID ()
int	HighestHadEtTowerID ()
void	SetDoBitInfo (bool doIt)
void	SetEMCrystalEnergy (int itwr, int icell, double en)
void	SetEtaPhiIndex (int eta, int phi, int ind)
void	SetFGBit (int twrid, bool FGBIT)
void	SetHCFGBit (int twrid, bool FGBIT)
void	SetHOEBit (int twrid, bool FGBIT)
void	SetJetE (double jE)
void	SetJetE3x3 (double jE)
void	SetJetEt (double jEt)
void	SetJetEt3x3 (double jEt)
void	SetParameters (L1Config)
void	SetRegionBits (edm::Event const &e)
void	SetRegionEnergy ()
void	SetTowerBits ()
double	SumE ()
double	SumEmE ()
double	SumEmEt ()
double	SumEt ()
double	SumHadE ()
double	SumHadEt ()
	~FastL1Region ()

Public Attributes
FastL1BitInfo	BitInfo
bool	doBitInfo

Private Member Functions
void	SetFGBit ()
void	SetHCFGBit ()
void	SetHOEBit ()
void	SetMIPBit ()
void	SetQuietBit ()
void	SetTauBit (edm::Event const &e)

Private Attributes
L1Config	Config
double	EMCrystalEnergy [16][25]
bool	fgBit [16]
bool	hcfgBit [16]
bool	hOeBit [16]
int	id
int	ieta
int	iphi
double	jetE
double	jetE3x3
double	jetEt
double	jetEt3x3
bool	mipBit
bool	quietBit
double	sumE
double	sumEt
bool	tauBit
CaloTowerCollection	Towers

Detailed Description

Description: Container class for L1 regions.

Implementation: <Notes on="" implementation>="">

Definition at line 110 of file FastL1Region.h.

Constructor & Destructor Documentation

FastL1Region::FastL1Region

(

)

Definition at line 22 of file FastL1Region.cc.

References Config, L1Config::CrystalEBThreshold, L1Config::CrystalEEThreshold, L1Config::EMActiveLevel, EMCrystalEnergy, L1Config::EMLSB, L1Config::EMNoiseLevel, L1Config::EMSeedEnThreshold, fgBit, L1Config::FGEBThreshold, L1Config::FGEEThreshold, L1Config::HadActiveLevel, L1Config::HadNoiseLevel, hcfgBit, hOeBit, L1Config::hOeThreshold, i, ieta, iphi, j, jetE, jetEt, L1Config::JetLSB, L1Config::JetSeedEtThreshold, mipBit, L1Config::MuonNoiseLevel, L1Config::noFGThreshold, L1Config::noTauVetoLevel, quietBit, L1Config::QuietRegionThreshold, tauBit, L1Config::TowerEBScale, L1Config::TowerEBThreshold, L1Config::TowerEEScale, L1Config::TowerEEThreshold, L1Config::TowerEMLSB, L1Config::TowerHadLSB, L1Config::TowerHBScale, L1Config::TowerHBThreshold, L1Config::TowerHEScale, L1Config::TowerHEThreshold, and Towers.

{
  Towers = CaloTowerCollection(16);

  jetE = 0.;
  jetEt = 0.;

  id = 999;
  ieta = 999;
  iphi = 999;

  tauBit = false;
  quietBit = false;
  mipBit = false;
  for(int i=0;i<16;i++) {
    hcfgBit[i] = false;
    fgBit[i] = false;
    hOeBit[i] = false;
    for (int j=0;j<25;j++) {
      EMCrystalEnergy[i][j] = 0. ; // 16x25 Crystals
    }
  }

  // default values
  Config.EMSeedEnThreshold = 2.;
  Config.EMActiveLevel = 3.;
  Config.HadActiveLevel = 3.;
  Config.noTauVetoLevel = 10000.;   
  Config.hOeThreshold = 0.05;
  Config.FGEBThreshold = 0.8;
  Config.noFGThreshold = 50.;
  Config.FGEEThreshold = 0.8;
  Config.MuonNoiseLevel = 2.;
  Config.EMNoiseLevel = 2.;
  Config.HadNoiseLevel = 2.;
  Config.QuietRegionThreshold = 2.;  
  Config.JetSeedEtThreshold = 2.;  
  Config.CrystalEBThreshold = 0.09;
  Config.CrystalEEThreshold = 0.45;

  Config.TowerEMLSB = 1.;
  Config.TowerHadLSB = 1.;
  Config.EMLSB = 1.;
  Config.JetLSB = 1.;

  Config.TowerEBThreshold = 0.2;
  Config.TowerEEThreshold = 0.45;
  Config.TowerHBThreshold = 0.9;
  Config.TowerHEThreshold = 1.4;

  Config.TowerEBScale = 1.0;
  Config.TowerEEScale = 1.0;
  Config.TowerHBScale = 1.0;
  Config.TowerHEScale = 1.0;


  //Config.EmInputs;
  //Config.xTowerInput;

}

FastL1Region::~FastL1Region

(

)

Definition at line 84 of file FastL1Region.cc.

{
}

Member Function Documentation

double FastL1Region::CalcSumE

(

)

Definition at line 690 of file FastL1Region.cc.

References i, and Towers.

Referenced by SetRegionEnergy(), and SumE().

{
  double sume=0;
  for (int i=0; i<16; i++) {
    sume += Towers[i].emEnergy();
    sume += Towers[i].hadEnergy();

  }
  return sume;
}

double FastL1Region::CalcSumEmE

(

)

Definition at line 702 of file FastL1Region.cc.

References i, and Towers.

Referenced by SumEmE().

{
  double sume=0;
  for (int i=0; i<16; i++) {
    sume += Towers[i].emEnergy();
  }
  return sume;
}

double FastL1Region::CalcSumEmEt

(

)

Definition at line 668 of file FastL1Region.cc.

References i, and Towers.

Referenced by SumEmEt().

{
  double sumet=0;
  for (int i=0; i<16; i++) {
    sumet += Towers[i].emEt();
  }

  return sumet;
}

double FastL1Region::CalcSumEt

(

)

Definition at line 656 of file FastL1Region.cc.

References i, and Towers.

Referenced by SetRegionEnergy(), and SumEt().

{
  double sumet=0;
  for (int i=0; i<16; i++) {
    sumet += Towers[i].emEt();
    sumet += Towers[i].hadEt();
  }

  return sumet;
}

double FastL1Region::CalcSumHadE

(

)

Definition at line 712 of file FastL1Region.cc.

References i, and Towers.

Referenced by SumHadE().

{
  double sume=0;
  for (int i=0; i<16; i++) {
    sume += Towers[i].hadEnergy();
  }
  return sume;
}

double FastL1Region::CalcSumHadEt

(

)

Definition at line 679 of file FastL1Region.cc.

References i, and Towers.

Referenced by SumHadEt().

{
  double sumet=0;
  for (int i=0; i<16; i++) {
    sumet += Towers[i].hadEt();
  }

  return sumet;
}

void FastL1Region::Dump

(

void

)

Definition at line 742 of file FastL1Region.cc.

References gather_cfg::cout, GetEastId(), GetId(), GetNEId(), GetNorthId(), GetNWId(), GetSEId(), GetSouthId(), GetSWId(), and GetWestId().

{

  // test tower filling:
  /*
  CaloTowerCollection::const_iterator t;
  int count = 0;
  for (t=Towers.begin(); t!=Towers.end(); t++) {
    std::cout << count << ") " << t->energy() << " | " << t->eta()  << " | " << t->phi() << std::endl;
    count++;
  }
  std::cout << std::endl;
  */

  // test region neighbours:
  std::cout << this->GetNWId() << " "  << this->GetNorthId() << " "  << this->GetNEId() << std::endl;
  std::cout << this->GetWestId() << " "  << this->GetId() << " "  << this->GetEastId() << std::endl;
  std::cout << this->GetSWId() << " "  << this->GetSouthId() << " "  << this->GetSEId() << std::endl;
  std::cout << std::endl;

}

void FastL1Region::FillEMCrystals	(	const CaloTowerConstituentsMap *	theTowerConstituentsMap,
		const CaloTopology *	calotopo,
		const CaloGeometry *	cGeom,
		const EcalRecHitCollection *	ec0,
		const EcalRecHitCollection *	ec1,
		FastL1RegionMap *	m_RMap
	)

Definition at line 120 of file FastL1Region.cc.

References edm::SortedCollection< T, SORT >::begin(), Config, L1Config::CrystalEBThreshold, L1Config::CrystalEEThreshold, CaloTopology::east(), EMCrystalEnergy, L1Config::EMNoiseLevel, edm::SortedCollection< T, SORT >::end(), relval_parameters_module::energy, fgBit, L1Config::FGEEThreshold, GetiEta(), i, CaloTowerDetId::ieta(), EBDetId::ieta(), ieta, CaloTowerDetId::iphi(), EBDetId::iphi(), iphi, L1Config::noFGThreshold, CaloTopology::north(), SetTowerBits(), CaloTopology::south(), EBDetId::tower_ieta(), EBDetId::tower_iphi(), L1Config::TowerEEThreshold, CaloTowerConstituentsMap::towerOf(), Towers, and CaloTopology::west().

{
  //std::vector< std::pair <std::string,std::string> > la;
  //la.resize(2);
  //la[0].first = "ecalRecHit";
  //la[0].second = "EcalRecHitsEB";
  //la[1].first = "ecalRecHit";
  //la[1].second = "EcalRecHitsEE";


  double ethres = Config.CrystalEBThreshold;

  // EB
  //e.getByLabel(la[0].first,la[0].second,ec);
  //e.getByLabel(Config.EmInputs.at(0),ec);

  ethres = Config.CrystalEBThreshold;
  for(EcalRecHitCollection::const_iterator ecItr = ec0->begin();
      ecItr != ec0->end(); ++ecItr) {
    //CaloRecHit recHit = (CaloRecHit)(*ecItr);
    if (ecItr->energy()<ethres) continue;

    EBDetId detId = ecItr->detid();

    //int hiphi = detId.tower_iphi();
    int hieta = detId.tower_ieta();
    int eieta = detId.ieta();
    int eiphi = detId.iphi();
    int crIeta = 999;
    if (hieta>0)
      crIeta = (eieta-1)%5;
    else
      crIeta = 4 + (eieta+1)%5;
    int crIphi = (eiphi - 1)%5;
    
    //const GlobalPoint gP = cGeom->getPosition(detId);

    //CaloTowerDetId towerDetId = theTowerConstituentsMap->towerOf(detId);
    // loop over towers
    for(int i=0;i<16;i++) {
      //int hiphi = m_RMap->convertFromECal_to_HCal_iphi(detId.tower_iphi());
      //int hiphi = m_RMap->convertFromHCal_to_ECal_iphi(detId.tower_iphi());
      int hiphi = detId.tower_iphi();
      if ( !Towers[i].id().iphi()==hiphi ||  !Towers[i].id().ieta()==hieta ) continue;
      EMCrystalEnergy[i][crIeta + 5*crIphi] = ecItr->energy();
    }  
  }
  
  // After having filled crsystal info set all veto bits
  SetTowerBits();
 
  // EE FG bits are filled here!!!
  //e.getByLabel(la[1].first,la[1].second,ec);

  if (GetiEta()==4 || GetiEta()==5 ||  GetiEta()==6 ||  
      GetiEta()==15 || GetiEta()==16 || GetiEta()==17 ) {
    
    //e.getByLabel(Config.EmInputs.at(1),ec);
    ethres = Config.CrystalEEThreshold;
    double towerEnergy[16];
    // loop over towers
    for(int i=0;i<16;i++) {
      fgBit[i] = false; // re-iniate
            
      //if (Towers[i].hadEt()>Config.HadNoiseLevel && Towers[i].emEt()>Config.EMNoiseLevel ) {
      if (Towers[i].emEt()>=Config.EMNoiseLevel ) {
      //if (Towers[i].emEnergy()>Config.EMNoiseLevel ) {
    //towerEnergy[i]  = Towers[i].hadEt() + Towers[i].emEt(); 
    towerEnergy[i]  = Towers[i].hadEnergy() + Towers[i].emEnergy(); 
      } else {
    fgBit[i] = false;
    continue;
      }

      // EB/EE transition area: unset fg bits
      // if (std::abs(Towers[i].id().ieta())==16 || std::abs(Towers[i].id().ieta())==17) {
      // fgBit[i] = false;
      // continue;
      // }
      if (Towers[i].emEt()>Config.noFGThreshold) {
    fgBit[i] = false;
    continue;
      }
      
      //CaloRecHit maxRecHit;
      //CaloRecHit maxRecHit2;
      double maxRecHit=-1.;
      double maxRecHit2=-1.;
      DetId maxDetId;

      double max2En = 0.;
      
      for(EcalRecHitCollection::const_iterator ecItr = ec1->begin();
      ecItr != ec1->end(); ++ecItr) {
    //CaloRecHit recHit = (CaloRecHit)(*ecItr);
    if (ecItr->energy()<ethres) continue;
    
    EEDetId detId = ecItr->detid();
    
    CaloTowerDetId towerDetId = theTowerConstituentsMap->towerOf(detId);
    //int hiphi = m_RMap->convertFromECal_to_HCal_iphi(towerDetId.iphi());
    int hiphi = towerDetId.iphi();
    if (Towers[i].id().iphi()==hiphi && 
        Towers[i].id().ieta()==towerDetId.ieta() ) {    
      if (maxRecHit<ecItr->energy()) {
        maxRecHit = ecItr->energy();
        maxDetId = detId;
      } 
    }
      } 

      std::vector<DetId> westV = calotopo->west(maxDetId);
      std::vector<DetId> eastV = calotopo->east(maxDetId);
      std::vector<DetId> southV = calotopo->south(maxDetId);
      std::vector<DetId> northV = calotopo->north(maxDetId);
      for(EcalRecHitCollection::const_iterator ecItr = ec1->begin();
      ecItr != ec1->end(); ++ecItr) {
    //CaloRecHit recHit = (CaloRecHit)(*ecItr);
    if (ecItr->energy()<ethres) continue;
    
    EEDetId detId = ecItr->detid();
    
    CaloTowerDetId towerDetId = theTowerConstituentsMap->towerOf(detId);
    //int hiphi = m_RMap->convertFromECal_to_HCal_iphi(towerDetId.iphi());
    int hiphi = towerDetId.iphi();
    if (Towers[i].id().iphi()==hiphi && 
        Towers[i].id().ieta()==towerDetId.ieta() ) {
      if ( 
          (!westV.empty() && detId==westV[0]) || 
          (!eastV.empty() && detId==eastV[0]) || 
          (!northV.empty() && detId==northV[0]) || 
          (!southV.empty() && detId==southV[0]) 
          ) {
        if (maxRecHit2<ecItr->energy()) {
          maxRecHit2 = ecItr->energy();
        }   
        max2En += ecItr->energy();
      }
    }
      }  
      
      double eeThres = Config.FGEEThreshold;
      //double totE = maxRecHit.energy() + max2En;
      double totE = maxRecHit + maxRecHit2;
      if (towerEnergy[i]>(Config.TowerEEThreshold)) {
    //double totE = maxRecHit.energy() + maxRecHit2.energy();
    //if (totE/towerEnergy[i]<Config.FGEBThreshold) fgBit[i] = true;    
    if (totE/towerEnergy[i]<eeThres) fgBit[i] = true;    
      }
            
    }
  }
   
 
}

void FastL1Region::FillTower	(	const CaloTower &	t,
		int &	tid,
		edm::ESHandle< CaloGeometry > &	cGeom
	)

Definition at line 292 of file FastL1Region.cc.

References abs, Config, CaloTower::emEt(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), CaloTower::hadEt(), CaloTower::id(), CaloTower::outerEt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), RCTEnergyTrunc(), L1Config::TowerEBThreshold, L1Config::TowerEEThreshold, L1Config::TowerEMLSB, L1Config::TowerHadLSB, L1Config::TowerHBThreshold, L1Config::TowerHEThreshold, and Towers.

{
  double EThres = 0.;
  double HThres = 0.;
  double EBthres = Config.TowerEBThreshold;
  double HBthres = Config.TowerHBThreshold;
  double EEthres = Config.TowerEEThreshold;
  double HEthres = Config.TowerHEThreshold;
  
  if(std::abs(t.eta())<2.322) {
    EThres = EBthres;
  } else {
    EThres = EEthres;
  }
  if(std::abs(t.eta())<2.322) {
    HThres = HBthres;
  } else {
    HThres = HEthres;
  }

  double upperThres = 1024.;
  double emet = RCTEnergyTrunc(t.emEt(),Config.TowerEMLSB,upperThres);
  double hadet = RCTEnergyTrunc(t.hadEt(),Config.TowerHadLSB,upperThres);
  //double eme = RCTEnergyTrunc(t.emEnergy(),Config.TowerEMLSB,upperThres);
  //double hade = RCTEnergyTrunc(t.hadEnergy(),Config.TowerHadLSB,upperThres);

  if ( emet<EThres) emet = 0.;
  if ( hadet<HThres) hadet = 0.;
  //if ( eme<EThres) emet = 0.;
  //if ( hade<HThres) hadet = 0.;
 
  GlobalPoint gP = cGeom->getPosition(t.id());
  math::XYZTLorentzVector lvec(t.px(),t.py(),t.px(),t.energy());
  //Towers[tid] = CaloTower(t);
  //Towers[tid] = CaloTower(t.id(),t.momentum(),emet,hadet,t.outerEt(),0,0);
  Towers[tid] = CaloTower(t.id(),emet,hadet,t.outerEt(),0,0,lvec,gP,gP);
}

void FastL1Region::FillTower_Scaled	(	const CaloTower &	t,
		int &	tid,
		bool	doRCTTrunc,
		edm::ESHandle< CaloGeometry > &	cGeom
	)

Definition at line 333 of file FastL1Region.cc.

References abs, Config, CaloTower::emEnergy(), CaloTower::emEt(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), CaloTower::hadEnergy(), CaloTower::hadEt(), CaloTower::id(), CaloTower::outerEt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), RCTEnergyTrunc(), L1Config::TowerEBScale, L1Config::TowerEBThreshold, L1Config::TowerEEScale, L1Config::TowerEEThreshold, L1Config::TowerEMLSB, L1Config::TowerHadLSB, L1Config::TowerHBScale, L1Config::TowerHBThreshold, L1Config::TowerHEScale, L1Config::TowerHEThreshold, and Towers.

{

  double EThres = 0.;
  double HThres = 0.;
  double EBthres = Config.TowerEBThreshold;
  double HBthres = Config.TowerHBThreshold;
  double EEthres = Config.TowerEEThreshold;
  double HEthres = Config.TowerHEThreshold;
  
  if(std::abs(t.eta())<2.322) {
    EThres = EBthres;
  } else {
    EThres = EEthres;
  }
  if(std::abs(t.eta())<2.322) {
    HThres = HBthres;
  } else {
    HThres = HEthres;
  }

  double emScale = 1.0;
  double hadScale = 1.0;
  //double outerScale = 1.0;

  if (std::abs(t.eta()>1.3050) && std::abs(t.eta())<3.0) {
    hadScale = Config.TowerHEScale;
    emScale = Config.TowerEEScale;
  }
  if (std::abs(t.eta()<1.3050)) {
    hadScale = Config.TowerHBScale;
    emScale = Config.TowerEBScale;
  }

  double emet = emScale * t.emEt();
  double hadet = hadScale * t.hadEt();
  double eme = emScale * t.emEnergy();
  double hade = hadScale * t.hadEnergy();

  if (doRCTTrunc) {
    double upperThres = 1024.;
    emet = RCTEnergyTrunc(emet,Config.TowerEMLSB,upperThres);
    hadet = RCTEnergyTrunc(hadet,Config.TowerHadLSB,upperThres);
    eme = RCTEnergyTrunc(eme,Config.TowerEMLSB,upperThres);
    hade = RCTEnergyTrunc(hade,Config.TowerHadLSB,upperThres);
  }
  if ( emet<EThres) emet = 0.;
  if ( hadet<HThres) hadet = 0.;
  //if ( eme<EThres) emet = 0.;
  //if ( hade<HThres) hadet = 0.;

  /* 
  if (t.emEt()>0. || t.hadEt()>0.) {
    std::cout<<"+++ "
         <<t.emEt()<<" "<<t.hadEt()<<" "
         <<t.eta()<<" "<<t.phi()<<" "
         <<std::endl;
  }
  */

  //Towers[tid] = CaloTower(t);
  //Towers[tid] = CaloTower(t.id(),t.momentum(),emet,hadet,0.,0,0);
  //edm::ESHandle<CaloGeometry> cGeom; 
  //c.get<CaloGeometryRecord>().get(cGeom);    
  GlobalPoint gP = cGeom->getPosition(t.id());
  math::XYZTLorentzVector lvec(t.px(),t.py(),t.px(),t.energy());
  //Towers[tid] = CaloTower(t);
  //Towers[tid] = CaloTower(t.id(),t.momentum(),emet,hadet,t.outerEt(),0,0);
  Towers[tid] = CaloTower(t.id(),emet,hadet,t.outerEt(),0,0,lvec,gP,gP);
  
  //std::cout<<tid<<"  "<<Towers[tid].emEt()<< " " <<Towers[tid].hadEt()<< std::endl;

}

void FastL1Region::FillTowerZero	(	const CaloTower &	t,
		int &	tid
	)

Definition at line 285 of file FastL1Region.cc.

References Towers.

{
  Towers[tid] = CaloTower(t);
  //std::cout<<"--- "<<Towers[tid].emEt()<<" "<<Towers[tid].hadEt()<<std::endl;
}

FastL1BitInfo FastL1Region::getBitInfo ( )

inline

Definition at line 192 of file FastL1Region.h.

References BitInfo.

{ return BitInfo; }

CaloTowerCollection FastL1Region::GetCaloTowers ( )

inline

Definition at line 156 of file FastL1Region.h.

References Towers.

{ return Towers; };

int FastL1Region::GetEastId

(

)

Definition at line 778 of file FastL1Region.cc.

References ieta, and iphi.

Referenced by Dump().

{ if (ieta != 21) return (iphi*22 + ieta+1); else return 999; }

double FastL1Region::GetEMCrystalEnergy ( int  itwr, int  icell  )

inline

Definition at line 158 of file FastL1Region.h.

References EMCrystalEnergy.

{ return EMCrystalEnergy[itwr][icell]; };

bool FastL1Region::GetFGBit ( int  i )

inline

Definition at line 152 of file FastL1Region.h.

References fgBit, and i.

{ if(i>=0 && i<16) {return fgBit[i];} else { return false; } };

bool FastL1Region::GetHCFGBit ( int  i )

inline

Definition at line 154 of file FastL1Region.h.

References hcfgBit, and i.

{  if(i>=0 && i<16) return hcfgBit[i]; else return false; };

bool FastL1Region::GetHOEBit ( int  i )

inline

Definition at line 153 of file FastL1Region.h.

References hOeBit, and i.

{ if(i>=0 && i<16) { return hOeBit[i]; } else { return false; } };

int FastL1Region::GetId ( )

inline

Definition at line 134 of file FastL1Region.h.

References id.

Referenced by Dump().

{ return id; };

int FastL1Region::GetiEta ( )

inline

Definition at line 132 of file FastL1Region.h.

References ieta.

Referenced by FillEMCrystals(), and SetFGBit().

{ return ieta; };

int FastL1Region::GetiPhi ( )

inline

Definition at line 133 of file FastL1Region.h.

References iphi.

{ return iphi; };

double FastL1Region::GetJetE ( )

inline

Definition at line 172 of file FastL1Region.h.

References jetE.

{ return jetE; };

double FastL1Region::GetJetE3x3 ( )

inline

Definition at line 183 of file FastL1Region.h.

References jetE3x3.

{ return jetE3x3; };

double FastL1Region::GetJetEt ( )

inline

Definition at line 173 of file FastL1Region.h.

References jetEt.

{ return jetEt; };

double FastL1Region::GetJetEt3x3 ( )

inline

Definition at line 184 of file FastL1Region.h.

References jetEt3x3.

{ return jetEt3x3; };

bool FastL1Region::GetMIPBit ( )

inline

Definition at line 151 of file FastL1Region.h.

References mipBit.

{ return mipBit; };

int FastL1Region::GetNEId

(

)

Definition at line 807 of file FastL1Region.cc.

References ieta, and iphi.

Referenced by Dump().

{ 
  if (ieta != 21) {
    if (iphi != 17) 
      return ((iphi+1)*22 + ieta+1); 
    else
      return (ieta+1);  
  } else {
    return 999; 
  }
}

int FastL1Region::GetNorthId

(

)

Definition at line 765 of file FastL1Region.cc.

References ieta, and iphi.

Referenced by Dump().

{ if (iphi != 17) return ((iphi+1)*22 + ieta); else return ieta; }

int FastL1Region::GetNWId

(

)

Definition at line 782 of file FastL1Region.cc.

References ieta, and iphi.

Referenced by Dump().

{ 
  if (ieta != 0) {
    if (iphi != 17) 
      return ((iphi+1)*22 + ieta-1); 
    else
      return (ieta-1);  
  } else {
    return 999; 
  }
}

bool FastL1Region::GetQuietBit ( )

inline

Definition at line 150 of file FastL1Region.h.

References quietBit.

{ return quietBit; };

std::pair< double, double > FastL1Region::getRegionCenterEtaPhi

(

const edm::EventSetup &

c

)

Definition at line 723 of file FastL1Region.cc.

References PV3DBase< T, PVType, FrameType >::eta(), eta(), edm::EventSetup::get(), phi, PV3DBase< T, PVType, FrameType >::phi(), and Towers.

{
  edm::ESHandle<CaloGeometry> cGeom; 
  //c.get<IdealGeometryRecord>().get(cGeom);    
  c.get<CaloGeometryRecord>().get(cGeom);    

  const GlobalPoint gP1 = cGeom->getPosition(Towers[5].id());
  //const GlobalPoint gP2 = cGeom->getPosition(Towers[6].id());
  //const GlobalPoint gP3 = cGeom->getPosition(Towers[10].id());

  double eta = gP1.eta();  
  double phi = gP1.phi();
  
  return std::pair<double, double>(eta, phi);
}

int FastL1Region::GetSEId

(

)

Definition at line 820 of file FastL1Region.cc.

References ieta, and iphi.

Referenced by Dump().

{ 
  if (ieta != 21) {
    if (iphi != 0) 
      return ((iphi-1)*22 + ieta+1); 
    else
      return (17*22 + ieta+1);  
  } else {
    return 999; 
  }
}

int FastL1Region::GetSouthId

(

)

Definition at line 770 of file FastL1Region.cc.

References ieta, and iphi.

Referenced by Dump().

{ if (iphi != 0) return ((iphi-1)*22 + ieta); else return (17*22 + ieta); }

int FastL1Region::GetSWId

(

)

Definition at line 795 of file FastL1Region.cc.

References ieta, and iphi.

Referenced by Dump().

{ 
  if (ieta != 0) {
    if (iphi != 0) 
      return ((iphi-1)*22 + ieta-1); 
    else
      return (17*22 + ieta-1);  
  } else {
    return 999; 
  }
}

bool FastL1Region::GetTauBit ( )

inline

Definition at line 149 of file FastL1Region.h.

References tauBit.

{ return tauBit; };

std::pair<int, int> FastL1Region::GetTowerNorthEtaPhi	(	int	ieta,
		int	iphi
	)

int FastL1Region::GetWestId

(

)

Definition at line 774 of file FastL1Region.cc.

References ieta, and iphi.

Referenced by Dump().

{ if (ieta != 0) return (iphi*22 + ieta-1); else return 999; }

int FastL1Region::HighestEmEtTowerID

(

)

Definition at line 624 of file FastL1Region.cc.

References BitInfo, doBitInfo, i, FastL1BitInfo::setHighestEmEtTowerID(), and Towers.

{
  int hid = -1;
  double tmpet=0.;
  for (int i=0; i<16; i++) {
    if ( (Towers[i].emEt()) > tmpet) {
      tmpet = (Towers[i].emEt());
      hid = i;
    }
  }

  if (doBitInfo) BitInfo.setHighestEmEtTowerID (hid);    
  return hid;
}

int FastL1Region::HighestEtTowerID

(

)

Definition at line 607 of file FastL1Region.cc.

References BitInfo, doBitInfo, i, FastL1BitInfo::setHighestEtTowerID(), and Towers.

{
  int hid = -1;
  double tmpet=0.;
  for (int i=0; i<16; i++) {
    if ( (Towers[i].emEt()+Towers[i].hadEt()) > tmpet) {
      tmpet = (Towers[i].emEt()+Towers[i].hadEt());
      hid = i;
    }
  }

  if (doBitInfo) BitInfo.setHighestEtTowerID (hid);      

  return hid;
}

int FastL1Region::HighestHadEtTowerID

(

)

Definition at line 640 of file FastL1Region.cc.

References BitInfo, doBitInfo, i, FastL1BitInfo::setHighestHadEtTowerID(), and Towers.

{
  int hid = -1;
  double tmpet=0.;
  for (int i=0; i<16; i++) {
    if ( (Towers[i].hadEt()) > tmpet) {
      tmpet = (Towers[i].hadEt());
      hid = i;
    }
  }

  if (doBitInfo) BitInfo.setHighestHadEtTowerID (hid);   
  return hid;
}

void FastL1Region::SetDoBitInfo ( bool  doIt )

inline

Definition at line 177 of file FastL1Region.h.

References doBitInfo.

{doBitInfo = doIt;}

void FastL1Region::SetEMCrystalEnergy ( int  itwr, int  icell, double  en  )

inline

Definition at line 157 of file FastL1Region.h.

References EMCrystalEnergy.

{ EMCrystalEnergy[itwr][icell] = en; };

void FastL1Region::SetEtaPhiIndex ( int  eta, int  phi, int  ind  )

inline

Definition at line 129 of file FastL1Region.h.

References eta(), ieta, iphi, and phi.

  { ieta=eta; iphi=phi; id=ind; };

void FastL1Region::SetFGBit	(	int	twrid,
		bool	FGBIT
	)

Definition at line 451 of file FastL1Region.cc.

References fgBit.

{
  fgBit[twrid] = FGBIT;
}

void FastL1Region::SetFGBit ( )

private

Definition at line 467 of file FastL1Region.cc.

References Config, EMCrystalEnergy, L1Config::EMNoiseLevel, fgBit, L1Config::FGEBThreshold, GetiEta(), i, j, L1Config::noFGThreshold, L1Config::TowerEBThreshold, and Towers.

Referenced by SetTowerBits().

{
  double ratioCut = Config.FGEBThreshold;

  double stripEnergy[16][5];
  double duostripEnergy[16][4];
  double towerEnergy[16];

  if (GetiEta()>=7 && GetiEta()<=14) {
    //Barrel
    for (int i=0; i<16; i++) {
      //if (Towers[i].hadEt()>Config.HadNoiseLevel && Towers[i].emEt()>Config.EMNoiseLevel ) {
      if (Towers[i].emEt()>Config.EMNoiseLevel ) {
    //towerEnergy[i]  = Towers[i].hadEt() + Towers[i].emEt(); 
    towerEnergy[i]  = Towers[i].hadEnergy() + Towers[i].emEnergy(); 
      } else {
    fgBit[i] = false;
    continue;
      }

      // EB/EE transition area: unset fg bits
      //if (std::abs(Towers[i].id().ieta())==16 || std::abs(Towers[i].id().ieta())==17) {
      //fgBit[i] = false;
      //continue;
      //}
      if (Towers[i].emEt()>Config.noFGThreshold) {
    fgBit[i] = false;
    continue;
      }

      bool fgflag = false;
    for (int j=0; j<5; j++) {
    stripEnergy[i][j] = EMCrystalEnergy[i][j] + EMCrystalEnergy[i][j+5] + EMCrystalEnergy[i][j+10] + 
      EMCrystalEnergy[i][j+15] + EMCrystalEnergy[i][j+20];
      }
      for (int j=0; j<4; j++) {
    duostripEnergy[i][j] = stripEnergy[i][j] + stripEnergy[i][j+1];
    if (towerEnergy[i]>(Config.TowerEBThreshold)) {
      //std::cout<<duostripEnergy[i][j]<<" |"<<towerEnergy[i]<<" |"<<duostripEnergy[i][j]/towerEnergy[i]<<std::endl;
      if ( (duostripEnergy[i][j] / towerEnergy[i]) > ratioCut) {
        fgflag = true;
      } 
      //std::cout<<duostripEnergy[i][j]<<" | "<<towerEnergy[i]<<": "<<duostripEnergy[i][j]/towerEnergy[i]<<std::endl;   
    }
      }

      if (fgflag) { 
    fgBit[i] = false;
      } else {
    fgBit[i] = true;
      }
      //std::cout<<GetiEta()<<" | "<<i<<": "<<fgBit[i]<<std::endl;
      //std::cout<<"********************************************"<<std::endl;
    }
  } else {
    // Endcap FG bit is already filled in fillEMCrystals()!!! 
  }
  
}

void FastL1Region::SetHCFGBit	(	int	twrid,
		bool	FGBIT
	)

Definition at line 456 of file FastL1Region.cc.

{
  ;
}

void FastL1Region::SetHCFGBit ( )

private

Definition at line 430 of file FastL1Region.cc.

Referenced by SetTowerBits().

{
  // temporary: check definition 
  // if (Tower->hadEt>100GeV) hcfgBit = true; ????
  //for (int i=0; i<16; i++) {
  //}
}

void FastL1Region::SetHOEBit	(	int	twrid,
		bool	FGBIT
	)

Definition at line 461 of file FastL1Region.cc.

References hOeBit.

{
  hOeBit[twrid] = FGBIT;
}

void FastL1Region::SetHOEBit ( )

private

Definition at line 408 of file FastL1Region.cc.

References Config, L1Config::EMNoiseLevel, L1Config::HadNoiseLevel, hOeBit, L1Config::hOeThreshold, i, and Towers.

Referenced by SetTowerBits().

{
  double fracThres = Config.hOeThreshold;

  for (int i=0; i<16; i++) {
    //if (Towers[i].hadEt()>Config.HadNoiseLevel && Towers[i].emEt()>Config.EMNoiseLevel ) {
    if (Towers[i].hadEnergy()>Config.HadNoiseLevel && Towers[i].emEnergy()>Config.EMNoiseLevel ) {
      if((Towers[i].hadEt()/Towers[i].emEt()) > fracThres) {
    hOeBit[i] = true;
      }
    }
  }
}

void FastL1Region::SetJetE ( double  jE )

inline

Definition at line 174 of file FastL1Region.h.

References jetE.

{ jetE = jE; };

void FastL1Region::SetJetE3x3 ( double  jE )

inline

Definition at line 185 of file FastL1Region.h.

References jetE3x3.

{ jetE3x3 = jE; };

void FastL1Region::SetJetEt ( double  jEt )

inline

Definition at line 175 of file FastL1Region.h.

References jetEt.

{ jetEt = jEt; };

void FastL1Region::SetJetEt3x3 ( double  jEt )

inline

Definition at line 186 of file FastL1Region.h.

References jetEt3x3.

{ jetEt3x3 = jEt; };

void FastL1Region::SetMIPBit ( )

private

Definition at line 439 of file FastL1Region.cc.

References hcfgBit, i, mipBit, and quietBit.

Referenced by SetRegionBits().

{
  if (quietBit)
  for (int i=0; i<16; i++) {
    if (hcfgBit) {
      mipBit = true;
      return;
    }
  }
}

void FastL1Region::SetParameters

(

L1Config

iconfig

)

Definition at line 90 of file FastL1Region.cc.

References Config.

{ 
  Config = iconfig;
}

void FastL1Region::SetQuietBit ( )

private

Definition at line 423 of file FastL1Region.cc.

References Config, quietBit, L1Config::QuietRegionThreshold, and SumEt().

Referenced by SetRegionBits().

{
  if (SumEt()<Config.QuietRegionThreshold)
    quietBit = true;
}

void FastL1Region::SetRegionBits

(

edm::Event const &

e

)

Definition at line 103 of file FastL1Region.cc.

References SetMIPBit(), SetQuietBit(), and SetTauBit().

{
  SetTauBit(e);
  SetQuietBit();
  SetMIPBit();
}

void FastL1Region::SetRegionEnergy

(

)

Definition at line 96 of file FastL1Region.cc.

References CalcSumE(), CalcSumEt(), sumE, and sumEt.

{
  sumE = CalcSumE();
  sumEt = CalcSumEt();
}

void FastL1Region::SetTauBit ( edm::Event const &  e )

private

Definition at line 529 of file FastL1Region.cc.

References BitInfo, Config, doBitInfo, L1Config::EMActiveLevel, L1Config::HadActiveLevel, i, FastL1BitInfo::setEmTauVeto(), FastL1BitInfo::setHadTauVeto(), FastL1BitInfo::setIsolationVeto(), FastL1BitInfo::setTauVeto(), tauBit, and Towers.

Referenced by SetRegionBits().

{
  float emThres = Config.EMActiveLevel;
  float hadThres = Config.HadActiveLevel;

  if (doBitInfo) BitInfo.setIsolationVeto(false);   
  if (doBitInfo) BitInfo.setTauVeto(false); 
  if (doBitInfo) BitInfo.setEmTauVeto(false);   
  if (doBitInfo) BitInfo.setHadTauVeto(false);  

  // init pattern containers
  unsigned emEtaPat = 0;
  unsigned emPhiPat = 0;
  unsigned hadEtaPat = 0;
  unsigned hadPhiPat = 0;
  unsigned one = 1;


  // fill hits as bit pattern
  for (int i=0; i<16; i++) {
    if(Towers[i].emEt() > emThres) {
      emEtaPat |= (one << (unsigned)(i%4));
      emPhiPat |= (one << (unsigned)(i/4));
    }

    if( Towers[i].hadEt() > hadThres) {
      hadEtaPat |= (one << (unsigned)(i%4));
      hadPhiPat |= (one << (unsigned)(i/4));
    }

  }

  // Patterns with two or less contiguous bits set are passed
  // rest are vetoed; 5=0101;7=0111;9=1001;10=1010;11=1011;13=1101;14=1110;15=1111
  //  --- Alternate patterns
  //  --- 9=1001;15=1111
  static std::vector<unsigned> vetoPatterns;
  if(vetoPatterns.size() == 0) {
    vetoPatterns.push_back(5);
    vetoPatterns.push_back(7);
    vetoPatterns.push_back(9);
    vetoPatterns.push_back(10);
    vetoPatterns.push_back(11);
    vetoPatterns.push_back(13);
    vetoPatterns.push_back(14);
    vetoPatterns.push_back(15);
  }


  for(std::vector<unsigned>::iterator i = vetoPatterns.begin();
      i != vetoPatterns.end();  i++) {
    unsigned etaPattern = emEtaPat | hadEtaPat;
    unsigned phiPattern = emPhiPat | hadPhiPat;

    //  em pattern
    if(emEtaPat == *i || emPhiPat == *i) {
      if (doBitInfo) BitInfo.setEmTauVeto(true);
    }
    //  had pattern
    if(hadEtaPat == *i || hadPhiPat == *i) {
      if (doBitInfo) BitInfo.setHadTauVeto(true);
    }

    if(etaPattern == *i || phiPattern == *i) // combined pattern
      //if(emEtaPat == *i || emPhiPat == *i || hadEtaPat == *i || hadPhiPat == *i)
      {
    tauBit = true;
    if (doBitInfo) BitInfo.setTauVeto(true);    
    return;
      }  
  }
  
  tauBit = false;
  
}

void FastL1Region::SetTowerBits

(

)

Definition at line 111 of file FastL1Region.cc.

References SetFGBit(), SetHCFGBit(), and SetHOEBit().

Referenced by FillEMCrystals().

{
  SetFGBit();
  SetHOEBit();
  SetHCFGBit();
}

double FastL1Region::SumE ( )

inline

Definition at line 160 of file FastL1Region.h.

References CalcSumE().

{ return CalcSumE(); };

double FastL1Region::SumEmE ( )

inline

Definition at line 162 of file FastL1Region.h.

References CalcSumEmE().

{ return CalcSumEmE(); };

double FastL1Region::SumEmEt ( )

inline

Definition at line 163 of file FastL1Region.h.

References CalcSumEmEt().

{ return CalcSumEmEt(); };

double FastL1Region::SumEt ( )

inline

Definition at line 161 of file FastL1Region.h.

References CalcSumEt().

Referenced by SetQuietBit().

{ return CalcSumEt(); };

double FastL1Region::SumHadE ( )

inline

Definition at line 164 of file FastL1Region.h.

References CalcSumHadE().

{ return CalcSumHadE(); };

double FastL1Region::SumHadEt ( )

inline

Definition at line 165 of file FastL1Region.h.

References CalcSumHadEt().

{ return CalcSumHadEt(); };

Member Data Documentation

FastL1BitInfo FastL1Region::BitInfo

Definition at line 196 of file FastL1Region.h.

Referenced by getBitInfo(), HighestEmEtTowerID(), HighestEtTowerID(), HighestHadEtTowerID(), and SetTauBit().

L1Config FastL1Region::Config

private

Definition at line 242 of file FastL1Region.h.

Referenced by FastL1Region(), FillEMCrystals(), FillTower(), FillTower_Scaled(), SetFGBit(), SetHOEBit(), SetParameters(), SetQuietBit(), and SetTauBit().

bool FastL1Region::doBitInfo

Definition at line 195 of file FastL1Region.h.

Referenced by HighestEmEtTowerID(), HighestEtTowerID(), HighestHadEtTowerID(), SetDoBitInfo(), and SetTauBit().

double FastL1Region::EMCrystalEnergy[16][25]

private

Definition at line 218 of file FastL1Region.h.

Referenced by FastL1Region(), FillEMCrystals(), GetEMCrystalEnergy(), SetEMCrystalEnergy(), and SetFGBit().

bool FastL1Region::fgBit[16]

private

Definition at line 235 of file FastL1Region.h.

Referenced by FastL1Region(), FillEMCrystals(), GetFGBit(), and SetFGBit().

bool FastL1Region::hcfgBit[16]

private

Definition at line 237 of file FastL1Region.h.

Referenced by FastL1Region(), GetHCFGBit(), and SetMIPBit().

bool FastL1Region::hOeBit[16]

private

Definition at line 236 of file FastL1Region.h.

Referenced by FastL1Region(), GetHOEBit(), and SetHOEBit().

int FastL1Region::id

private

Definition at line 227 of file FastL1Region.h.

Referenced by GetId(), and NodeCut.MatchId::match().

int FastL1Region::ieta

private

Definition at line 228 of file FastL1Region.h.

Referenced by FastL1Region(), FillEMCrystals(), GetEastId(), GetiEta(), GetNEId(), GetNorthId(), GetNWId(), GetSEId(), GetSouthId(), GetSWId(), GetWestId(), and SetEtaPhiIndex().

int FastL1Region::iphi

private

Definition at line 229 of file FastL1Region.h.

Referenced by FastL1Region(), FillEMCrystals(), GetEastId(), GetiPhi(), GetNEId(), GetNorthId(), GetNWId(), GetSEId(), GetSouthId(), GetSWId(), GetWestId(), and SetEtaPhiIndex().

double FastL1Region::jetE

private

Definition at line 222 of file FastL1Region.h.

Referenced by FastL1Region(), GetJetE(), and SetJetE().

double FastL1Region::jetE3x3

private

Definition at line 224 of file FastL1Region.h.

Referenced by GetJetE3x3(), and SetJetE3x3().

double FastL1Region::jetEt

private

Definition at line 223 of file FastL1Region.h.

Referenced by FastL1Region(), GetJetEt(), and SetJetEt().

double FastL1Region::jetEt3x3

private

Definition at line 225 of file FastL1Region.h.

Referenced by GetJetEt3x3(), and SetJetEt3x3().

bool FastL1Region::mipBit

private

Definition at line 233 of file FastL1Region.h.

Referenced by FastL1Region(), GetMIPBit(), and SetMIPBit().

bool FastL1Region::quietBit

private

Definition at line 232 of file FastL1Region.h.

Referenced by FastL1Region(), GetQuietBit(), SetMIPBit(), and SetQuietBit().

double FastL1Region::sumE

private

Definition at line 240 of file FastL1Region.h.

Referenced by SetRegionEnergy().

double FastL1Region::sumEt

private

Definition at line 239 of file FastL1Region.h.

Referenced by SetRegionEnergy().

bool FastL1Region::tauBit

private

Definition at line 231 of file FastL1Region.h.

Referenced by FastL1Region(), GetTauBit(), and SetTauBit().

CaloTowerCollection FastL1Region::Towers

private

Definition at line 216 of file FastL1Region.h.

Referenced by CalcSumE(), CalcSumEmE(), CalcSumEmEt(), CalcSumEt(), CalcSumHadE(), CalcSumHadEt(), FastL1Region(), FillEMCrystals(), FillTower(), FillTower_Scaled(), FillTowerZero(), GetCaloTowers(), getRegionCenterEtaPhi(), HighestEmEtTowerID(), HighestEtTowerID(), HighestHadEtTowerID(), SetFGBit(), SetHOEBit(), and SetTauBit().