|
|
|
#include <L1TdeRCT.h>
Definition at line 104 of file L1TdeRCT.h.
| L1TdeRCT::L1TdeRCT | ( | const edm::ParameterSet & | ps | ) |
Definition at line 105 of file L1TdeRCT.cc.
References gather_cfg::cout, dbe, edm::ParameterSet::getUntrackedParameter(), histFolder_, NULL, cmsCodeRules::cppFunctionSkipper::operator, outputFile_, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), singlechannelhistos_, and verbose_.
: rctSourceEmul_( ps.getParameter< InputTag >("rctSourceEmul") ), rctSourceData_( ps.getParameter< InputTag >("rctSourceData") ), ecalTPGData_( ps.getParameter< InputTag >("ecalTPGData") ), hcalTPGData_( ps.getParameter< InputTag >("hcalTPGData") ), gtDigisLabel_( ps.getParameter< InputTag >("gtDigisLabel") ), gtEGAlgoName_ ( ps.getParameter< std::string >("gtEGAlgoName") ), doubleThreshold_ ( ps.getParameter< int >("doubleThreshold") ) { singlechannelhistos_ = ps.getUntrackedParameter < bool > ("singlechannelhistos", false); if (singlechannelhistos_) if(verbose_) std::cout << "L1TdeRCT: single channels histos ON" << std::endl; // verbosity switch verbose_ = ps.getUntrackedParameter < bool > ("verbose", false); if (verbose_) std::cout << "L1TdeRCT: constructor...." << std::endl; dbe = NULL; if (ps.getUntrackedParameter < bool > ("DQMStore", false)) { dbe = Service < DQMStore > ().operator->(); dbe->setVerbose(0); } outputFile_ = ps.getUntrackedParameter < std::string > ("outputFile", ""); if (outputFile_.size() != 0) { if(verbose_) std:: cout << "L1T Monitoring histograms will be saved to " << outputFile_.c_str() << std::endl; } bool disable = ps.getUntrackedParameter < bool > ("disableROOToutput", false); if (disable) { outputFile_ = ""; } histFolder_ = ps.getUntrackedParameter<std::string>("HistFolder", "L1TEMU/L1TdeRCT/"); if (dbe != NULL) { dbe->setCurrentFolder(histFolder_); } }
| L1TdeRCT::~L1TdeRCT | ( | ) | [virtual] |
Definition at line 160 of file L1TdeRCT.cc.
{
}
| void L1TdeRCT::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
Implements edm::EDAnalyzer.
Definition at line 810 of file L1TdeRCT.cc.
References gather_cfg::cout, DivideME1D(), DivideME2D(), doubleThreshold_, ecalTPGData_, funct::false, MonitorElement::Fill(), HcalObjRepresent::Fill(), first, newFWLiteAna::found, edm::Event::getByLabel(), gtDigisLabel_, gtEGAlgoName_, hcalTPGData_, i, edm::HandleBase::isValid(), j, gen::k, edm::InputTag::label(), nev_, notrigCount, PHIBINS, PhiEtaMax, rctBitDataHfPlusTau2D_, rctBitDataMip2D_, rctBitDataOverFlow2D_, rctBitDataQuiet2D_, rctBitDataTauVeto2D_, rctBitEmulHfPlusTau2D_, rctBitEmulMip2D_, rctBitEmulOverFlow2D_, rctBitEmulQuiet2D_, rctBitEmulTauVeto2D_, rctBitHfPlusTauEff2D_, rctBitHfPlusTauIneff2D_, rctBitHfPlusTauOvereff2D_, rctBitMatchedHfPlusTau2D_, rctBitMatchedMip2D_, rctBitMatchedOverFlow2D_, rctBitMatchedQuiet2D_, rctBitMatchedTauVeto2D_, rctBitMipEff2D_, rctBitMipIneff2D_, rctBitMipOvereff2D_, rctBitOverFlowEff2D_, rctBitOverFlowIneff2D_, rctBitOverFlowOvereff2D_, rctBitTauVetoEff2D_, rctBitTauVetoIneff2D_, rctBitTauVetoOvereff2D_, rctBitUnmatchedDataHfPlusTau2D_, rctBitUnmatchedDataMip2D_, rctBitUnmatchedDataOverFlow2D_, rctBitUnmatchedDataQuiet2D_, rctBitUnmatchedDataTauVeto2D_, rctBitUnmatchedEmulHfPlusTau2D_, rctBitUnmatchedEmulMip2D_, rctBitUnmatchedEmulOverFlow2D_, rctBitUnmatchedEmulQuiet2D_, rctBitUnmatchedEmulTauVeto2D_, rctInputTPGEcalOcc_, rctInputTPGEcalRank_, rctInputTPGHcalOcc_, rctInputTPGHcalRank_, rctInputTPGHcalSample_, rctIsoEffChannel_, rctIsoEmDataOcc1D_, rctIsoEmDataOcc_, rctIsoEmEff1_, rctIsoEmEff1Occ1D_, rctIsoEmEff1Occ_, rctIsoEmEff1oneD_, rctIsoEmEff2_, rctIsoEmEff2Occ1D_, rctIsoEmEff2Occ_, rctIsoEmEff2oneD_, rctIsoEmEmulOcc1D_, rctIsoEmEmulOcc_, rctIsoEmIneff1D_, rctIsoEmIneff2_, rctIsoEmIneff2Occ1D_, rctIsoEmIneff2Occ_, rctIsoEmIneff2oneD_, rctIsoEmIneff_, rctIsoEmIneffOcc1D_, rctIsoEmIneffOcc_, rctIsoEmOvereff1D_, rctIsoEmOvereff_, rctIsoEmOvereffOcc1D_, rctIsoEmOvereffOcc_, rctIsoIneffChannel_, rctIsoOvereffChannel_, rctNisoEffChannel_, rctNisoEmDataOcc1D_, rctNisoEmDataOcc_, rctNisoEmEff1_, rctNisoEmEff1Occ1D_, rctNisoEmEff1Occ_, rctNisoEmEff1oneD_, rctNisoEmEff2_, rctNisoEmEff2Occ1D_, rctNisoEmEff2Occ_, rctNisoEmEff2oneD_, rctNisoEmEmulOcc1D_, rctNisoEmEmulOcc_, rctNisoEmIneff1D_, rctNisoEmIneff2_, rctNisoEmIneff2Occ1D_, rctNisoEmIneff2Occ_, rctNisoEmIneff2oneD_, rctNisoEmIneff_, rctNisoEmIneffOcc1D_, rctNisoEmIneffOcc_, rctNisoEmOvereff1D_, rctNisoEmOvereff_, rctNisoEmOvereffOcc1D_, rctNisoEmOvereffOcc_, rctNisoIneffChannel_, rctNisoOvereffChannel_, rctRegDataOcc1D_, rctRegDataOcc2D_, rctRegEff1D_, rctRegEff2D_, rctRegEffChannel_, rctRegEmulOcc1D_, rctRegEmulOcc2D_, rctRegIneff1D_, rctRegIneff2D_, rctRegIneffChannel_, rctRegMatchedOcc1D_, rctRegMatchedOcc2D_, rctRegOvereff1D_, rctRegOvereff2D_, rctRegOvereffChannel_, rctRegSpEff1D_, rctRegSpEff2D_, rctRegSpEffOcc1D_, rctRegSpEffOcc2D_, rctRegSpIneff1D_, rctRegSpIneff2D_, rctRegSpIneffOcc1D_, rctRegSpIneffOcc2D_, rctRegUnmatchedDataOcc1D_, rctRegUnmatchedDataOcc2D_, rctRegUnmatchedEmulOcc1D_, rctRegUnmatchedEmulOcc2D_, rctSourceData_, rctSourceEmul_, singlechannelhistos_, trigCount, trigEff_, trigEffOcc_, trigEffThresh_, trigEffThreshOcc_, trigEffTriggOcc_, trigEffTriggThreshOcc_, triggerAlgoNumbers_, and verbose_.
{
nev_++;
if (verbose_) {
std::cout << "L1TdeRCT: analyze...." << std::endl;
}
// for GT decision word
// edm::ESHandle<L1GtTriggerMenu> menuRcd;
edm::Handle< L1GlobalTriggerReadoutRecord > gtRecord;
// get GT trigger menu, maps algorithms to the bits read out
// c.get<L1GtTriggerMenuRcd>().get(menuRcd) ;
// const L1GtTriggerMenu* menu = menuRcd.product();
// get GT decision word
//e.getByLabel( edm::InputTag("gtDigis"), gtRecord);
e.getByLabel( gtDigisLabel_ , gtRecord );
const DecisionWord dWord = gtRecord->decisionWord(); // this will get the decision word *before* masking disabled bits
int effEGThresholdBitNumber = 999;
if (gtEGAlgoName_ == "L1_SingleEG1")
{
effEGThresholdBitNumber = 46;
}
if (gtEGAlgoName_ == "L1_SingleEG5_0001")
{
effEGThresholdBitNumber = 47;
}
if (gtEGAlgoName_ == "L1_SingleEG8_0001")
{
effEGThresholdBitNumber = 48;
}
if (gtEGAlgoName_ == "L1_SingleEG10_0001")
{
effEGThresholdBitNumber = 49;
}
if (gtEGAlgoName_ == "L1_SingleEG12_0001")
{
effEGThresholdBitNumber = 50;
}
if (gtEGAlgoName_ == "L1_SingleEG15_0001")
{
effEGThresholdBitNumber = 51;
}
if (gtEGAlgoName_ == "L1_SingleEG20_0001")
{
effEGThresholdBitNumber = 52;
}
int algoBitNumber = 0;
bool triggered = false;
bool independent_triggered = false;
DecisionWord::const_iterator algoItr;
for (algoItr = dWord.begin(); algoItr != dWord.end(); algoItr++)
{
if (*algoItr)
{
triggerAlgoNumbers_->Fill(algoBitNumber);
if (algoBitNumber == effEGThresholdBitNumber)
{
triggered = true;// Fill triggered events (numerator) here!
}
if (algoBitNumber <= 45 || algoBitNumber >= 53)
{
independent_triggered = true;// use the muon path only !
}
}
algoBitNumber++;
}
if(triggered)
trigCount++;
else
notrigCount++;
// get TPGs
edm::Handle<EcalTrigPrimDigiCollection> ecalTpData;
edm::Handle<HcalTrigPrimDigiCollection> hcalTpData;
// Get the RCT digis
edm::Handle < L1CaloEmCollection > emData;
edm::Handle < L1CaloRegionCollection > rgnData;
// Get the RCT digis
edm::Handle < L1CaloEmCollection > emEmul;
edm::Handle < L1CaloRegionCollection > rgnEmul;
// need to change to getByLabel
bool doEm = true;
bool doHd = true;
bool doEcal = true;
bool doHcal = true;
// TPG, first try:
e.getByLabel(ecalTPGData_,ecalTpData);
e.getByLabel(hcalTPGData_,hcalTpData);
if (!ecalTpData.isValid()) {
edm::LogInfo("TPG DataNotFound") << "can't find EcalTrigPrimDigiCollection with label "
<< ecalTPGData_.label() ;
if (verbose_)std::cout << "Can not find ecalTpData!" << std::endl ;
doEcal = false ;
}
if(doEcal)
{
for(EcalTrigPrimDigiCollection::const_iterator iEcalTp = ecalTpData->begin(); iEcalTp != ecalTpData->end(); iEcalTp++)
if(iEcalTp->compressedEt() > 0)
{
rctInputTPGEcalRank_ -> Fill(1.*(iEcalTp->compressedEt())) ;
if(iEcalTp->id().ieta() > 0)
rctInputTPGEcalOcc_ -> Fill(1.*(iEcalTp->id().ieta())-0.5,iEcalTp->id().iphi()) ;
else
rctInputTPGEcalOcc_ -> Fill(1.*(iEcalTp->id().ieta())+0.5,iEcalTp->id().iphi()) ;
if(verbose_) std::cout << " ECAL data: Energy: " << iEcalTp->compressedEt() << " eta " << iEcalTp->id().ieta() << " phi " << iEcalTp->id().iphi() << std::endl ;
}
}
if (!hcalTpData.isValid()) {
edm::LogInfo("TPG DataNotFound") << "can't find HcalTrigPrimDigiCollection with label "
<< hcalTPGData_.label() ;
if (verbose_)std::cout << "Can not find hcalTpData!" << std::endl ;
doHcal = false ;
}
if(doHcal)
{
for(HcalTrigPrimDigiCollection::const_iterator iHcalTp = hcalTpData->begin(); iHcalTp != hcalTpData->end(); iHcalTp++)
{
int highSample=0;
int highEt=0;
for (int nSample = 0; nSample < 10; nSample++)
{
if (iHcalTp->sample(nSample).compressedEt() != 0)
{
if(verbose_) std::cout << "HCAL data: Et "
<< iHcalTp->sample(nSample).compressedEt()
<< " fg "
<< iHcalTp->sample(nSample).fineGrain()
<< " ieta " << iHcalTp->id().ieta()
<< " iphi " << iHcalTp->id().iphi()
<< " sample " << nSample
<< std::endl ;
if (iHcalTp->sample(nSample).compressedEt() > highEt)
{
highSample = nSample;
highEt = iHcalTp->sample(nSample).compressedEt() ;
}
}
}
if(highEt != 0)
{
if(iHcalTp->id().ieta() > 0)
rctInputTPGHcalOcc_ -> Fill(1.*(iHcalTp->id().ieta())-0.5,iHcalTp->id().iphi()) ;
else
rctInputTPGHcalOcc_ -> Fill(1.*(iHcalTp->id().ieta())+0.5,iHcalTp->id().iphi()) ;
rctInputTPGHcalSample_ -> Fill(highSample,highEt) ;
rctInputTPGHcalRank_ -> Fill(highEt) ;
}
}
}
e.getByLabel(rctSourceData_,rgnData);
e.getByLabel(rctSourceEmul_,rgnEmul);
if (!rgnData.isValid()) {
edm::LogInfo("DataNotFound") << "can't find L1CaloRegionCollection with label "
<< rctSourceData_.label() ;
if (verbose_)std::cout << "Can not find rgnData!" << std::endl ;
doHd = false;
}
// if ( doHd ) {
if (!rgnEmul.isValid()) {
edm::LogInfo("DataNotFound") << "can't find L1CaloRegionCollection with label "
<< rctSourceEmul_.label() ;
doHd = false;
if (verbose_)std::cout << "Can not find rgnEmul!" << std::endl ;
}
// }
e.getByLabel(rctSourceData_,emData);
e.getByLabel(rctSourceEmul_,emEmul);
if (!emData.isValid()) {
edm::LogInfo("DataNotFound") << "can't find L1CaloEmCollection with label "
<< rctSourceData_.label() ;
if (verbose_)std::cout << "Can not find emData!" << std::endl ;
doEm = false;
}
// if ( doEm ) {
if (!emEmul.isValid()) {
edm::LogInfo("DataNotFound") << "can't find L1CaloEmCollection with label "
<< rctSourceEmul_.label() ;
if (verbose_)std::cout << "Can not find emEmul!" << std::endl ;
doEm = false; return ;
}
// }
// Isolated and non-isolated EM
// StepI: Reset
int nelectrIsoData = 0;
int nelectrNisoData = 0;
int nelectrIsoEmul = 0;
int nelectrNisoEmul = 0;
int electronDataRank[2][PhiEtaMax]={{0}};
int electronDataEta[2][PhiEtaMax]={{0}};
int electronDataPhi[2][PhiEtaMax]={{0}};
int electronEmulRank[2][PhiEtaMax]={{0}};
int electronEmulEta[2][PhiEtaMax]={{0}};
int electronEmulPhi[2][PhiEtaMax]={{0}};
// region/bit arrays
int nRegionData = 0;
int nRegionEmul = 0;
int regionDataRank[PhiEtaMax] = {0};
int regionDataEta [PhiEtaMax] = {0};
int regionDataPhi [PhiEtaMax] = {0};
bool regionDataOverFlow [PhiEtaMax] = {false};
bool regionDataTauVeto [PhiEtaMax] = {false};
bool regionDataMip [PhiEtaMax] = {false};
bool regionDataQuiet [PhiEtaMax] = {false};
bool regionDataHfPlusTau[PhiEtaMax] = {false};
int regionEmulRank[PhiEtaMax] = {0};
int regionEmulEta [PhiEtaMax] = {0};
int regionEmulPhi [PhiEtaMax] = {0};
bool regionEmulOverFlow [PhiEtaMax] = {false};
bool regionEmulTauVeto [PhiEtaMax] = {false};
bool regionEmulMip [PhiEtaMax] = {false};
bool regionEmulQuiet [PhiEtaMax] = {false};
bool regionEmulHfPlusTau[PhiEtaMax] = {false};
if(first)
{
first = false ;
}
// StepII: fill variables
for (L1CaloEmCollection::const_iterator iem = emEmul->begin();
iem != emEmul->end();
iem++)
{
if(iem->rank() >= 1)
{
if(iem->isolated())
{
rctIsoEmEmulOcc_->Fill(iem->regionId().ieta(), iem->regionId().iphi());
// to show bad channles in the 2D efficiency plots
rctIsoEmIneffOcc_->Fill (iem->regionId().ieta(), iem->regionId().iphi(), 0.01);
rctIsoEmEff1Occ_->Fill (iem->regionId().ieta(), iem->regionId().iphi(), 0.01);
int channel;
channel=PHIBINS*iem->regionId().ieta()+iem->regionId().iphi();
rctIsoEmEmulOcc1D_->Fill(channel);
electronEmulRank[0][nelectrIsoEmul]=iem->rank();
electronEmulEta[0][nelectrIsoEmul]=iem->regionId().ieta();
electronEmulPhi[0][nelectrIsoEmul]=iem->regionId().iphi();
nelectrIsoEmul++ ;
}
else
{
rctNisoEmEmulOcc_->Fill(iem->regionId().ieta(), iem->regionId().iphi());
// to show bad channles in the 2D efficiency plots
rctNisoEmIneffOcc_->Fill (iem->regionId().ieta(), iem->regionId().iphi(), 0.01);
rctNisoEmEff1Occ_->Fill (iem->regionId().ieta(), iem->regionId().iphi(), 0.01);
int channel;
//
channel=PHIBINS*iem->regionId().ieta()+iem->regionId().iphi();
rctNisoEmEmulOcc1D_->Fill(channel);
electronEmulRank[1][nelectrNisoEmul]=iem->rank();
electronEmulEta[1][nelectrNisoEmul]=iem->regionId().ieta();
electronEmulPhi[1][nelectrNisoEmul]=iem->regionId().iphi();
nelectrNisoEmul++ ;
}
}
}
for (L1CaloEmCollection::const_iterator iem = emData->begin();
iem != emData->end();
iem++)
{
if(iem->rank() >= 1)
{
if (iem->isolated())
{
rctIsoEmDataOcc_->Fill(iem->regionId().ieta(), iem->regionId().iphi());
// new stuff to avoid 0's in emulator 2D //
// rctIsoEmEmulOcc_->Fill(iem->regionId().ieta(), iem->regionId().iphi(),0.01);
rctIsoEmOvereffOcc_->Fill (iem->regionId().ieta(),
iem->regionId().iphi(), 0.01);
int channel;
channel=PHIBINS*iem->regionId().ieta()+iem->regionId().iphi();
rctIsoEmDataOcc1D_->Fill(channel);
// new stuff to avoid 0's
// rctIsoEmEmulOcc1D_->Fill(channel);
electronDataRank[0][nelectrIsoData]=iem->rank();
electronDataEta[0][nelectrIsoData]=iem->regionId().ieta();
electronDataPhi[0][nelectrIsoData]=iem->regionId().iphi();
nelectrIsoData++ ;
}
else
{
rctNisoEmDataOcc_->Fill(iem->regionId().ieta(), iem->regionId().iphi());
// new stuff to avoid 0's in emulator 2D //
// rctNisoEmEmulOcc_->Fill(iem->regionId().ieta(), iem->regionId().iphi(),0.01);
rctNisoEmOvereffOcc_->Fill (iem->regionId().ieta(),
iem->regionId().iphi(), 0.01);
int channel;
channel=PHIBINS*iem->regionId().ieta()+iem->regionId().iphi();
rctNisoEmDataOcc1D_->Fill(channel);
// new stuff to avoid 0's
// rctNisoEmEmulOcc1D_->Fill(channel);
electronDataRank[1][nelectrNisoData]=iem->rank();
electronDataEta[1][nelectrNisoData]=iem->regionId().ieta();
electronDataPhi[1][nelectrNisoData]=iem->regionId().iphi();
nelectrNisoData++ ;
}
}
}
// fill region/bit arrays for emulator
for(L1CaloRegionCollection::const_iterator ireg = rgnEmul->begin();
ireg != rgnEmul->end();
ireg++)
{
// std::cout << "Emul: " << nRegionEmul << " " << ireg->gctEta() << " " << ireg->gctPhi() << std::endl;
if(ireg->overFlow()) rctBitEmulOverFlow2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->tauVeto()) rctBitEmulTauVeto2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->mip()) rctBitEmulMip2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->quiet()) rctBitEmulQuiet2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->fineGrain()) rctBitEmulHfPlusTau2D_->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->et() > 0)
{
rctRegEmulOcc1D_->Fill(PHIBINS*ireg->gctEta() + ireg->gctPhi());
rctRegEmulOcc2D_->Fill(ireg->gctEta(), ireg->gctPhi());
}
// to show bad channels in 2D efficiency plots:
if(ireg->overFlow()) {
rctBitUnmatchedEmulOverFlow2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
rctBitMatchedOverFlow2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
}
if(ireg->tauVeto()) {
rctBitUnmatchedEmulTauVeto2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
rctBitMatchedTauVeto2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
}
if(ireg->mip()) {
rctBitUnmatchedEmulMip2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
rctBitMatchedMip2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
}
if(ireg->quiet()) {
rctBitUnmatchedEmulQuiet2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
rctBitMatchedQuiet2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
}
if(ireg->fineGrain()) {
rctBitUnmatchedEmulHfPlusTau2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
rctBitMatchedHfPlusTau2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
}
if(ireg->et() > 0) {
rctRegUnmatchedEmulOcc2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
rctRegMatchedOcc2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01);
/* rctRegDeltaEtOcc2D_->Fill (ireg->gctEta(), ireg->gctPhi(), 0.01); */
}
nRegionEmul = PHIBINS * ireg->gctEta() + ireg->gctPhi();
regionEmulRank [nRegionEmul] = ireg->et();
regionEmulEta [nRegionEmul] = ireg->gctEta();
regionEmulPhi [nRegionEmul] = ireg->gctPhi();
regionEmulOverFlow [nRegionEmul] = ireg->overFlow();
regionEmulTauVeto [nRegionEmul] = ireg->tauVeto();
regionEmulMip [nRegionEmul] = ireg->mip();
regionEmulQuiet [nRegionEmul] = ireg->quiet();
regionEmulHfPlusTau[nRegionEmul] = ireg->fineGrain();
}
// fill region/bit arrays for hardware
for(L1CaloRegionCollection::const_iterator ireg = rgnData->begin();
ireg != rgnData->end();
ireg++)
{
// std::cout << "Data: " << nRegionData << " " << ireg->gctEta() << " " << ireg->gctPhi() << std::endl;
if(ireg->overFlow()) rctBitDataOverFlow2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->tauVeto()) rctBitDataTauVeto2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->mip()) rctBitDataMip2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->quiet()) rctBitDataQuiet2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->fineGrain()) rctBitDataHfPlusTau2D_->Fill(ireg->gctEta(), ireg->gctPhi());
if(ireg->et() > 0)
{
rctRegDataOcc1D_ ->Fill(PHIBINS*ireg->gctEta() + ireg->gctPhi());
rctRegDataOcc2D_ ->Fill(ireg->gctEta(), ireg->gctPhi());
}
// to show bad channels in 2D inefficiency:
// if(ireg->overFlow()) rctBitEmulOverFlow2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
// if(ireg->tauVeto()) rctBitEmulTauVeto2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
// if(ireg->mip()) rctBitEmulMip2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
// if(ireg->quiet()) rctBitEmulQuiet2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
// if(ireg->fineGrain()) rctBitEmulHfPlusTau2D_->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
// if(ireg->et() > 0) rctRegEmulOcc2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
if(ireg->overFlow()) rctBitUnmatchedDataOverFlow2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
if(ireg->tauVeto()) rctBitUnmatchedDataTauVeto2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
if(ireg->mip()) rctBitUnmatchedDataMip2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
if(ireg->quiet()) rctBitUnmatchedDataQuiet2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
if(ireg->fineGrain()) rctBitUnmatchedDataHfPlusTau2D_->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
if(ireg->et() > 0) rctRegUnmatchedDataOcc2D_ ->Fill(ireg->gctEta(), ireg->gctPhi(), 0.01);
nRegionData = PHIBINS * ireg->gctEta() + ireg->gctPhi();
regionDataRank [nRegionData] = ireg->et();
regionDataEta [nRegionData] = ireg->gctEta();
regionDataPhi [nRegionData] = ireg->gctPhi();
regionDataOverFlow [nRegionData] = ireg->overFlow();
regionDataTauVeto [nRegionData] = ireg->tauVeto();
regionDataMip [nRegionData] = ireg->mip();
regionDataQuiet [nRegionData] = ireg->quiet();
regionDataHfPlusTau[nRegionData] = ireg->fineGrain();
}
if(verbose_)
{
std::cout << "I found Data! Iso: " << nelectrIsoData << " Niso: " << nelectrNisoData << std::endl ;
for(int i=0; i<nelectrIsoData; i++)
std::cout << " Iso Energy " << electronDataRank[0][i] << " eta " << electronDataEta[0][i] << " phi " << electronDataPhi[0][i] << std::endl ;
for(int i=0; i<nelectrNisoData; i++)
std::cout << " Niso Energy " << electronDataRank[1][i] << " eta " << electronDataEta[1][i] << " phi " << electronDataPhi[1][i] << std::endl ;
std::cout << "I found Emul! Iso: " << nelectrIsoEmul << " Niso: " << nelectrNisoEmul << std::endl ;
for(int i=0; i<nelectrIsoEmul; i++)
std::cout << " Iso Energy " << electronEmulRank[0][i] << " eta " << electronEmulEta[0][i] << " phi " << electronEmulPhi[0][i] << std::endl ;
for(int i=0; i<nelectrNisoEmul; i++)
std::cout << " Niso Energy " << electronEmulRank[1][i] << " eta " << electronEmulEta[1][i] << " phi " << electronEmulPhi[1][i] << std::endl ;
std::cout << "I found Data! Regions: " << PhiEtaMax << std::endl ;
for(int i=0; i<(int)PhiEtaMax; i++)
if(regionDataRank[i] !=0 ) std::cout << " Energy " << regionDataRank[i] << " eta " << regionDataEta[i] << " phi " << regionDataPhi[i] << std::endl ;
std::cout << "I found Emul! Regions: " << PhiEtaMax << std::endl ;
for(int i=0; i<(int)PhiEtaMax; i++)
if(regionEmulRank[i] !=0 ) std::cout << " Energy " << regionEmulRank[i] << " eta " << regionEmulEta[i] << " phi " << regionEmulPhi[i] << std::endl ;
}
// StepIII: calculate and fill
for(int k=0; k<2; k++)
{
int nelectrE, nelectrD;
if(k==0)
{
nelectrE=nelectrIsoEmul;
nelectrD=nelectrIsoData;
}
else
{
nelectrE=nelectrNisoEmul;
nelectrD=nelectrNisoData;
}
for(int i = 0; i < nelectrE; i++)
{
//bool triggered = l1SingleEG2; //false; //HACK until true trigger implimented
double trigThresh = doubleThreshold_; //ditto
if(singlechannelhistos_) {
int chnl=PHIBINS*electronEmulEta[k][i]+electronEmulPhi[k][i];
if(k==1 && independent_triggered) { //non-iso
//std::cout << "eta " << electronEmulEta[k][i] << " phi " << electronEmulPhi[k][i] << " with rank " << electronEmulRank[k][i] << std::endl;
trigEffOcc_[chnl]->Fill(electronEmulRank[k][i]);
// }
if(triggered)
trigEffTriggOcc_[chnl]->Fill(electronEmulRank[k][i]); }
}
//find number of objects with rank above 2x trigger threshold
//and number after requiring a trigger too
if(electronEmulRank[k][i]>=trigThresh){
if(k==1 && independent_triggered) { //non-iso
trigEffThreshOcc_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i]);
trigEffTriggThreshOcc_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.01);
// }
if(triggered)
trigEffTriggThreshOcc_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.98001); }
}
Bool_t found = kFALSE;
for(int j = 0; j < nelectrD; j++)
{
if(electronEmulEta[k][i]==electronDataEta[k][j] &&
electronEmulPhi[k][i]==electronDataPhi[k][j])
{
if(k==0)
{
rctIsoEmEff1Occ_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.98001);
// Weight is for ROOT; when added to initial weight of 0.01, should just exceed 0.99
int chnl;
chnl=PHIBINS*electronEmulEta[k][i]+electronEmulPhi[k][i];
rctIsoEmEff1Occ1D_->Fill(chnl);
if(singlechannelhistos_)
{
int energy_difference;
energy_difference=(electronEmulRank[k][i] - electronDataRank[k][j]);
rctIsoEffChannel_[chnl]->Fill(energy_difference);
}
if(electronEmulRank[k][i]==electronDataRank[k][j])
{
rctIsoEmEff2Occ1D_->Fill(chnl);
rctIsoEmEff2Occ_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.98012);
// Weight is for ROOT; should just exceed 0.99
// NOTE: Weight is different for eff 2 because this isn't filled initially
// for current definition of Eff2 and Ineff2 we need to add additional
// factor 0.99 since we divide over eff1 which is 0.99001 e.g. we use 0.99001**2 !
rctIsoEmIneff2Occ_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.0099);
}
else
{
rctIsoEmIneff2Occ1D_->Fill(chnl);
rctIsoEmIneff2Occ_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.9801);
}
}
else
{
rctNisoEmEff1Occ_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.98001);
// Weight is for ROOT; when added to initial weight of 0.01, should just exceed 0.99
int chnl;
chnl=PHIBINS*electronEmulEta[k][i]+electronEmulPhi[k][i];
rctNisoEmEff1Occ1D_->Fill(chnl);
if(singlechannelhistos_)
{
int energy_difference;
energy_difference=(electronEmulRank[k][i] - electronDataRank[k][j]) ;
rctNisoEffChannel_[chnl]->Fill(energy_difference) ;
}
if(electronEmulRank[k][i]==electronDataRank[k][j])
{
rctNisoEmEff2Occ1D_->Fill(chnl);
rctNisoEmEff2Occ_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.98012);
// Weight is for ROOT; should just exceed 0.99
// NOTE: Weight is different for eff 2 because this isn't filled initially
// see comments fo Iso
rctNisoEmIneff2Occ_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.0099);
}
else
{
rctNisoEmIneff2Occ1D_->Fill(chnl);
rctNisoEmIneff2Occ_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.9801);
}
}
found = kTRUE;
}
}
if(found == kFALSE)
{
if(k==0)
{
rctIsoEmIneffOcc_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.98);
// Weight is for ROOT; when added to initial weight of 0.01, should equal 0.99
int chnl;
chnl=PHIBINS*electronEmulEta[k][i]+electronEmulPhi[k][i];
rctIsoEmIneffOcc1D_->Fill(chnl);
if(singlechannelhistos_)
{
rctIsoIneffChannel_[chnl]->Fill(electronEmulRank[k][i]);
}
}
else
{
rctNisoEmIneffOcc_->Fill(electronEmulEta[k][i], electronEmulPhi[k][i], 0.98);
// Weight is for ROOT; when added to initial weight of 0.01, should equal 0.99
int chnl;
chnl=PHIBINS*electronEmulEta[k][i]+electronEmulPhi[k][i];
rctNisoEmIneffOcc1D_->Fill(chnl);
if(singlechannelhistos_)
{
rctNisoIneffChannel_[chnl]->Fill(electronEmulRank[k][i]);
}
}
}
}
DivideME1D(rctIsoEmEff1Occ1D_, rctIsoEmEmulOcc1D_, rctIsoEmEff1oneD_);
DivideME2D(rctIsoEmEff1Occ_, rctIsoEmEmulOcc_, rctIsoEmEff1_) ;
// DivideME1D(rctIsoEmEff2Occ1D_, rctIsoEmEmulOcc1D_, rctIsoEmEff2oneD_);
// DivideME2D(rctIsoEmEff2Occ_, rctIsoEmEmulOcc_, rctIsoEmEff2_) ;
DivideME1D(rctIsoEmEff2Occ1D_, rctIsoEmEff1Occ1D_, rctIsoEmEff2oneD_);
DivideME2D(rctIsoEmEff2Occ_, rctIsoEmEff1Occ_, rctIsoEmEff2_) ;
// DivideME1D(rctIsoEmIneff2Occ1D_, rctIsoEmEmulOcc1D_, rctIsoEmIneff2oneD_);
// DivideME2D(rctIsoEmIneff2Occ_, rctIsoEmEmulOcc_, rctIsoEmIneff2_) ;
DivideME1D(rctIsoEmIneff2Occ1D_, rctIsoEmEff1Occ1D_, rctIsoEmIneff2oneD_);
DivideME2D(rctIsoEmIneff2Occ_, rctIsoEmEff1Occ_, rctIsoEmIneff2_) ;
DivideME1D(rctNisoEmEff1Occ1D_, rctNisoEmEmulOcc1D_, rctNisoEmEff1oneD_);
DivideME2D(rctNisoEmEff1Occ_, rctNisoEmEmulOcc_, rctNisoEmEff1_);
// DivideME1D(rctNisoEmEff2Occ1D_, rctNisoEmEmulOcc1D_, rctNisoEmEff2oneD_);
// DivideME2D(rctNisoEmEff2Occ_, rctNisoEmEmulOcc_, rctNisoEmEff2_);
DivideME1D(rctNisoEmEff2Occ1D_, rctNisoEmEff1Occ1D_, rctNisoEmEff2oneD_);
DivideME2D(rctNisoEmEff2Occ_, rctNisoEmEff1Occ_, rctNisoEmEff2_);
// DivideME1D(rctNisoEmIneff2Occ1D_, rctNisoEmEmulOcc1D_, rctNisoEmIneff2oneD_);
// DivideME2D(rctNisoEmIneff2Occ_, rctNisoEmEmulOcc_, rctNisoEmIneff2_);
DivideME1D(rctNisoEmIneff2Occ1D_, rctNisoEmEff1Occ1D_, rctNisoEmIneff2oneD_);
DivideME2D(rctNisoEmIneff2Occ_, rctNisoEmEff1Occ_, rctNisoEmIneff2_);
DivideME1D(rctIsoEmIneffOcc1D_, rctIsoEmEmulOcc1D_, rctIsoEmIneff1D_);
DivideME2D(rctIsoEmIneffOcc_, rctIsoEmEmulOcc_, rctIsoEmIneff_);
DivideME1D(rctNisoEmIneffOcc1D_, rctNisoEmEmulOcc1D_, rctNisoEmIneff1D_);
DivideME2D(rctNisoEmIneffOcc_, rctNisoEmEmulOcc_, rctNisoEmIneff_);
DivideME2D(trigEffTriggThreshOcc_, trigEffThreshOcc_, trigEffThresh_);
if(singlechannelhistos_) {
for(int i = 0; i < nelectrE; i++)
{
int chnl=PHIBINS*electronEmulEta[k][i]+electronEmulPhi[k][i];
DivideME1D(trigEffTriggOcc_[chnl], trigEffOcc_[chnl], trigEff_[chnl]);
}
}
for(int i = 0; i < nelectrD; i++)
{
Bool_t found = kFALSE;
for(int j = 0; j < nelectrE; j++)
{
if(electronEmulEta[k][j]==electronDataEta[k][i] &&
electronEmulPhi[k][j]==electronDataPhi[k][i])
{
found = kTRUE;
}
}
if(found == kFALSE)
{
if(k==0)
{
rctIsoEmOvereffOcc_->Fill(electronDataEta[k][i], electronDataPhi[k][i], 0.98);
// Weight is for ROOT; when added to initial weight of 0.01, should equal 0.99
int chnl;
chnl=PHIBINS*electronDataEta[k][i]+electronDataPhi[k][i];
rctIsoEmOvereffOcc1D_->Fill(chnl);
if(singlechannelhistos_)
{
rctIsoOvereffChannel_[chnl]->Fill(electronDataRank[k][i]);
}
}
else
{
rctNisoEmOvereffOcc_->Fill(electronDataEta[k][i], electronDataPhi[k][i], 0.98);
// Weight is for ROOT; when added to initial weight of 0.01, should equal 0.99
int chnl;
chnl=PHIBINS*electronDataEta[k][i]+electronDataPhi[k][i];
rctNisoEmOvereffOcc1D_->Fill(chnl) ;
if(singlechannelhistos_)
{
rctNisoOvereffChannel_[chnl]->Fill(electronDataRank[k][i]);
}
}
}
}
}
// we try new definition of overefficiency:
DivideME1D(rctIsoEmOvereffOcc1D_, rctIsoEmDataOcc1D_, rctIsoEmOvereff1D_);
DivideME2D(rctIsoEmOvereffOcc_, rctIsoEmDataOcc_, rctIsoEmOvereff_);
DivideME1D(rctNisoEmOvereffOcc1D_, rctNisoEmDataOcc1D_, rctNisoEmOvereff1D_);
DivideME2D(rctNisoEmOvereffOcc_, rctNisoEmDataOcc_, rctNisoEmOvereff_);
// calculate region/bit information
for(unsigned int i = 0; i < (int)PhiEtaMax; i++)
{
Bool_t regFound = kFALSE;
Bool_t overFlowFound = kFALSE;
Bool_t tauVetoFound = kFALSE;
Bool_t mipFound = kFALSE;
Bool_t quietFound = kFALSE;
Bool_t hfPlusTauFound = kFALSE;
// for(int j = 0; j < nRegionData; j++)
// {
// if(regionEmulEta[i] == regionDataEta[j] &&
// regionEmulPhi[i] == regionDataPhi[j])
// {
if(regionDataRank[i] >= 1 && regionEmulRank[i] >= 1)
{
int chnl;
chnl = PHIBINS*regionEmulEta[i] + regionEmulPhi[i];
rctRegMatchedOcc1D_->Fill(chnl);
rctRegMatchedOcc2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.98001);
// Weight is for ROOT; when added to initial weight of 0.01, should just exceed 0.99
if(singlechannelhistos_) rctRegEffChannel_[chnl]->Fill(regionEmulRank[i] - regionDataRank[i]);
// see comments for Iso Eff2
if(regionEmulRank[i] == regionDataRank[i])
{
rctRegSpEffOcc1D_->Fill(chnl);
// rctRegSpEffOcc2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.99001);
rctRegSpEffOcc2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.98012);
rctRegSpIneffOcc2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.0099);
}
else
{
rctRegSpIneffOcc1D_->Fill(chnl);
rctRegSpIneffOcc2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.9801);
}
// Weight is for ROOT; should just exceed 0.99
// NOTE: Weight is different for eff 2 because this isn't filled initially
regFound = kTRUE;
}
if(regionEmulOverFlow[i] == true &&
regionDataOverFlow[i] == true)
{
rctBitMatchedOverFlow2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.98001);
overFlowFound = kTRUE;
}
if(regionEmulTauVeto[i] == true &&
regionDataTauVeto[i] == true)
{
rctBitMatchedTauVeto2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.98001);
tauVetoFound = kTRUE;
}
if (regionEmulMip[i] == true && regionDataMip[i] == true) {
rctBitMatchedMip2D_->Fill (regionEmulEta[i], regionEmulPhi[i], 0.98001);
mipFound = kTRUE;
}
if (regionEmulQuiet[i] == true && regionDataQuiet[i] == true) {
rctBitMatchedQuiet2D_->Fill (regionEmulEta[i], regionEmulPhi[i], 0.98001);
quietFound = kTRUE;
}
if (regionEmulHfPlusTau[i] == true && regionDataHfPlusTau[i] == true) {
rctBitMatchedHfPlusTau2D_->Fill (regionEmulEta[i], regionEmulPhi[i], 0.98001);
hfPlusTauFound = kTRUE;
}
// }
// }
if(regFound == kFALSE && regionEmulRank[i] >= 1 )
{
int chnl;
chnl = PHIBINS*regionEmulEta[i] + regionEmulPhi[i];
rctRegUnmatchedEmulOcc1D_->Fill(chnl);
rctRegUnmatchedEmulOcc2D_->Fill(regionEmulEta[i], regionEmulPhi[i],0.98);
// Weight is for ROOT; when added to initial weight of 0.01, should equal 0.99
if(singlechannelhistos_) rctRegIneffChannel_[chnl]->Fill(regionEmulRank[i]);
}
if(overFlowFound == kFALSE && regionEmulOverFlow[i] == true)
{
rctBitUnmatchedEmulOverFlow2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.98);
}
if(tauVetoFound == kFALSE && regionEmulTauVeto[i] == true)
{
rctBitUnmatchedEmulTauVeto2D_->Fill(regionEmulEta[i], regionEmulPhi[i], 0.98);
}
if (mipFound == kFALSE && regionEmulMip[i] == true) {
rctBitUnmatchedEmulMip2D_->Fill (regionEmulEta[i], regionEmulPhi[i], 0.98);
}
if (quietFound == kFALSE && regionEmulQuiet[i] == true) {
rctBitUnmatchedEmulQuiet2D_->Fill (regionEmulEta[i], regionEmulPhi[i], 0.98);
}
if (hfPlusTauFound == kFALSE && regionEmulHfPlusTau[i] == true ) {
rctBitUnmatchedEmulHfPlusTau2D_->Fill (regionEmulEta[i], regionEmulPhi[i], 0.98);
}
}
DivideME1D(rctRegMatchedOcc1D_, rctRegEmulOcc1D_, rctRegEff1D_);
DivideME2D(rctRegMatchedOcc2D_, rctRegEmulOcc2D_, rctRegEff2D_);
// DivideME1D(rctRegSpEffOcc1D_, rctRegEmulOcc1D_, rctRegSpEff1D_);
// DivideME2D(rctRegSpEffOcc2D_, rctRegEmulOcc2D_, rctRegSpEff2D_);
DivideME1D(rctRegSpEffOcc1D_, rctRegMatchedOcc1D_, rctRegSpEff1D_);
DivideME2D(rctRegSpEffOcc2D_, rctRegMatchedOcc2D_, rctRegSpEff2D_);
// DivideME1D(rctRegSpIneffOcc1D_, rctRegEmulOcc1D_, rctRegSpIneff1D_);
// DivideME2D(rctRegSpIneffOcc2D_, rctRegEmulOcc2D_, rctRegSpIneff2D_);
DivideME1D(rctRegSpIneffOcc1D_, rctRegMatchedOcc1D_, rctRegSpIneff1D_);
DivideME2D(rctRegSpIneffOcc2D_, rctRegMatchedOcc2D_, rctRegSpIneff2D_);
DivideME2D(rctBitMatchedOverFlow2D_, rctBitEmulOverFlow2D_, rctBitOverFlowEff2D_);
DivideME2D(rctBitMatchedTauVeto2D_, rctBitEmulTauVeto2D_, rctBitTauVetoEff2D_);
DivideME2D (rctBitMatchedMip2D_, rctBitEmulMip2D_, rctBitMipEff2D_);
// QUIETBIT: To add quiet bit information, uncomment following line:
// DivideME2D (rctBitMatchedQuiet2D_, rctBitEmulQuiet2D_, rctBitQuietEff2D_);
DivideME2D (rctBitMatchedHfPlusTau2D_, rctBitEmulHfPlusTau2D_, rctBitHfPlusTauEff2D_);
DivideME1D (rctRegUnmatchedEmulOcc1D_, rctRegEmulOcc1D_, rctRegIneff1D_);
DivideME2D (rctRegUnmatchedEmulOcc2D_, rctRegEmulOcc2D_, rctRegIneff2D_);
DivideME2D (rctBitUnmatchedEmulOverFlow2D_, rctBitEmulOverFlow2D_, rctBitOverFlowIneff2D_);
DivideME2D (rctBitUnmatchedEmulTauVeto2D_, rctBitEmulTauVeto2D_, rctBitTauVetoIneff2D_);
DivideME2D (rctBitUnmatchedEmulMip2D_, rctBitEmulMip2D_, rctBitMipIneff2D_);
// QUIETBIT: To add quiet bit information, uncomment the following line:
// DivideME2D (rctBitUnmatchedEmulQuiet2D_, rctBitEmulQuiet2D_, rctBitQuietIneff2D_);
DivideME2D (rctBitUnmatchedEmulHfPlusTau2D_, rctBitEmulHfPlusTau2D_, rctBitHfPlusTauIneff2D_);
// for(int i = 0; i < nRegionData; i++)
for (int i = 0; i < (int)PhiEtaMax; i++)
{
Bool_t regFound = kFALSE;
Bool_t overFlowFound = kFALSE;
Bool_t tauVetoFound = kFALSE;
Bool_t mipFound = kFALSE;
Bool_t quietFound = kFALSE;
Bool_t hfPlusTauFound = kFALSE;
// for(int j = 0; j < nRegionEmul; j++)
// {
// if(regionEmulEta[j] == regionDataEta[i] &&
// regionEmulPhi[j] == regionDataPhi[i])
// {
if(regionEmulRank[i] >= 1 && regionDataRank[i] >= 1)
regFound = kTRUE;
if(regionDataOverFlow[i] == true &&
regionEmulOverFlow[i] == true)
overFlowFound = kTRUE;
if(regionDataTauVeto[i] == true &&
regionEmulTauVeto[i] == true)
tauVetoFound = kTRUE;
if (regionDataMip[i] == true && regionEmulMip[i] == true)
mipFound = kTRUE;
if (regionDataQuiet[i] == true && regionEmulQuiet[i] == true)
quietFound = kTRUE;
if (regionDataHfPlusTau[i] == true && regionEmulHfPlusTau[i] == true)
hfPlusTauFound = kTRUE;
// }
// }
if(regFound == kFALSE && regionDataRank[i] >= 1)
{
int chnl;
chnl = PHIBINS*regionDataEta[i] + regionDataPhi[i];
rctRegUnmatchedDataOcc1D_->Fill(chnl);
rctRegUnmatchedDataOcc2D_->Fill(regionDataEta[i], regionDataPhi[i], 0.98);
// Weight is for ROOT; when added to initial weight of 0.01, should equal 0.99
// we try a new definition of overefficiency:
// DivideME1D(rctRegUnmatchedDataOcc1D_, rctRegDataOcc1D_, rctRegOvereff1D_);
// DivideME2D(rctRegUnmatchedDataOcc2D_, rctRegDataOcc2D_, rctRegOvereff2D_);
if(singlechannelhistos_) rctRegOvereffChannel_[chnl]->Fill(regionDataRank[i]);
}
if(overFlowFound == kFALSE && regionDataOverFlow[i] == true )
{
rctBitUnmatchedDataOverFlow2D_->Fill(regionDataEta[i], regionDataPhi[i], 0.98);
}
if(tauVetoFound == kFALSE && regionDataTauVeto[i] == true )
{
rctBitUnmatchedDataTauVeto2D_->Fill(regionDataEta[i], regionDataPhi[i], 0.98);
}
if (mipFound == kFALSE && regionDataMip[i] == true ) {
rctBitUnmatchedDataMip2D_->Fill (regionDataEta[i], regionDataPhi[i], 0.98);
}
if (quietFound == kFALSE && regionDataQuiet[i] == true ) {
rctBitUnmatchedDataQuiet2D_->Fill (regionDataEta[i], regionDataPhi[i], 0.98);
}
if (hfPlusTauFound == kFALSE && regionDataHfPlusTau[i] == true ) {
rctBitUnmatchedDataHfPlusTau2D_->Fill (regionDataEta[i], regionDataPhi[i], 0.98);
}
}
// we try a new definition of overefficiency:
DivideME1D(rctRegUnmatchedDataOcc1D_, rctRegDataOcc1D_, rctRegOvereff1D_);
DivideME2D(rctRegUnmatchedDataOcc2D_, rctRegDataOcc2D_, rctRegOvereff2D_);
DivideME2D(rctBitUnmatchedDataOverFlow2D_, rctBitDataOverFlow2D_, rctBitOverFlowOvereff2D_);
DivideME2D(rctBitUnmatchedDataTauVeto2D_, rctBitDataTauVeto2D_, rctBitTauVetoOvereff2D_);
DivideME2D (rctBitUnmatchedDataMip2D_, rctBitDataMip2D_,
rctBitMipOvereff2D_);
// QUIETBIT: To add quiet bit information, uncomment following 2 lines:
// DivideME2D (rctBitUnmatchedDataQuiet2D_, rctBitDataQuiet2D_,
// rctBitQuietOvereff2D_);
DivideME2D (rctBitUnmatchedDataHfPlusTau2D_, rctBitDataHfPlusTau2D_,
rctBitHfPlusTauOvereff2D_);
}
| void L1TdeRCT::beginJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 164 of file L1TdeRCT.cc.
References DQMStore::book1D(), DQMStore::book2D(), CHNLBINS, CHNLMAX, CHNLMIN, crateFED, dbe, DEBINS, DEMAX, DEMIN, ELBINS, ELMAX, ELMIN, ETABINS, ETAMAX, ETAMIN, fedVectorMonitorLS_, fedVectorMonitorRUN_, MonitorElement::getTH2F(), histFolder_, i, j, m, mergeVDriftHistosByStation::name, nev_, notrigCount, cmsCodeRules::cppFunctionSkipper::operator, PHIBINS, PHIMAX, PHIMIN, rctBitDataHfPlusTau2D_, rctBitDataMip2D_, rctBitDataOverFlow2D_, rctBitDataQuiet2D_, rctBitDataTauVeto2D_, rctBitEmulHfPlusTau2D_, rctBitEmulMip2D_, rctBitEmulOverFlow2D_, rctBitEmulQuiet2D_, rctBitEmulTauVeto2D_, rctBitHfPlusTauEff2D_, rctBitHfPlusTauIneff2D_, rctBitHfPlusTauOvereff2D_, rctBitMatchedHfPlusTau2D_, rctBitMatchedMip2D_, rctBitMatchedOverFlow2D_, rctBitMatchedQuiet2D_, rctBitMatchedTauVeto2D_, rctBitMipEff2D_, rctBitMipIneff2D_, rctBitMipOvereff2D_, rctBitOverFlowEff2D_, rctBitOverFlowIneff2D_, rctBitOverFlowOvereff2D_, rctBitTauVetoEff2D_, rctBitTauVetoIneff2D_, rctBitTauVetoOvereff2D_, rctBitUnmatchedDataHfPlusTau2D_, rctBitUnmatchedDataMip2D_, rctBitUnmatchedDataOverFlow2D_, rctBitUnmatchedDataQuiet2D_, rctBitUnmatchedDataTauVeto2D_, rctBitUnmatchedEmulHfPlusTau2D_, rctBitUnmatchedEmulMip2D_, rctBitUnmatchedEmulOverFlow2D_, rctBitUnmatchedEmulQuiet2D_, rctBitUnmatchedEmulTauVeto2D_, rctInputTPGEcalOcc_, rctInputTPGEcalRank_, rctInputTPGHcalOcc_, rctInputTPGHcalRank_, rctInputTPGHcalSample_, rctIsoEffChannel_, rctIsoEmDataOcc1D_, rctIsoEmDataOcc_, rctIsoEmEff1_, rctIsoEmEff1Occ1D_, rctIsoEmEff1Occ_, rctIsoEmEff1oneD_, rctIsoEmEff2_, rctIsoEmEff2Occ1D_, rctIsoEmEff2Occ_, rctIsoEmEff2oneD_, rctIsoEmEmulOcc1D_, rctIsoEmEmulOcc_, rctIsoEmIneff1D_, rctIsoEmIneff2_, rctIsoEmIneff2Occ1D_, rctIsoEmIneff2Occ_, rctIsoEmIneff2oneD_, rctIsoEmIneff_, rctIsoEmIneffOcc1D_, rctIsoEmIneffOcc_, rctIsoEmOvereff1D_, rctIsoEmOvereff_, rctIsoEmOvereffOcc1D_, rctIsoEmOvereffOcc_, rctIsoIneffChannel_, rctIsoOvereffChannel_, rctNisoEffChannel_, rctNisoEmDataOcc1D_, rctNisoEmDataOcc_, rctNisoEmEff1_, rctNisoEmEff1Occ1D_, rctNisoEmEff1Occ_, rctNisoEmEff1oneD_, rctNisoEmEff2_, rctNisoEmEff2Occ1D_, rctNisoEmEff2Occ_, rctNisoEmEff2oneD_, rctNisoEmEmulOcc1D_, rctNisoEmEmulOcc_, rctNisoEmIneff1D_, rctNisoEmIneff2_, rctNisoEmIneff2Occ1D_, rctNisoEmIneff2Occ_, rctNisoEmIneff2oneD_, rctNisoEmIneff_, rctNisoEmIneffOcc1D_, rctNisoEmIneffOcc_, rctNisoEmOvereff1D_, rctNisoEmOvereff_, rctNisoEmOvereffOcc1D_, rctNisoEmOvereffOcc_, rctNisoIneffChannel_, rctNisoOvereffChannel_, rctRegDataOcc1D_, rctRegDataOcc2D_, rctRegEff1D_, rctRegEff2D_, rctRegEffChannel_, rctRegEmulOcc1D_, rctRegEmulOcc2D_, rctRegIneff1D_, rctRegIneff2D_, rctRegIneffChannel_, rctRegMatchedOcc1D_, rctRegMatchedOcc2D_, rctRegOvereff1D_, rctRegOvereff2D_, rctRegOvereffChannel_, rctRegSpEff1D_, rctRegSpEff2D_, rctRegSpEffOcc1D_, rctRegSpEffOcc2D_, rctRegSpIneff1D_, rctRegSpIneff2D_, rctRegSpIneffOcc1D_, rctRegSpIneffOcc2D_, rctRegUnmatchedDataOcc1D_, rctRegUnmatchedDataOcc2D_, rctRegUnmatchedEmulOcc1D_, rctRegUnmatchedEmulOcc2D_, DQMStore::rmdir(), DQMStore::setCurrentFolder(), singlechannelhistos_, TPGETABINS, TPGETAMAX, TPGETAMIN, TPGPHIBINS, TPGPHIMAX, TPGPHIMIN, TPGRANK, TPGRANKMAX, TPGRANKMIN, trigCount, trigEff_, trigEffOcc_, trigEffThresh_, trigEffThreshOcc_, trigEffTriggOcc_, trigEffTriggThreshOcc_, and triggerAlgoNumbers_.
{
nev_ = 0;
// get hold of back-end interface
DQMStore *dbe = 0;
dbe = Service < DQMStore > ().operator->();
if (dbe) {
dbe->setCurrentFolder(histFolder_);
dbe->rmdir(histFolder_);
}
if (dbe) {
dbe->setCurrentFolder(histFolder_);
triggerAlgoNumbers_ =
dbe->book1D("gtTriggerAlgoNumbers", "gtTriggerAlgoNumbers", 128, -0.5, 127.5);
rctInputTPGEcalOcc_ =
dbe->book2D("rctInputTPGEcalOcc", "rctInputTPGEcalOcc", TPGETABINS, TPGETAMIN,
TPGETAMAX, TPGPHIBINS, TPGPHIMIN, TPGPHIMAX);
rctInputTPGEcalRank_ =
dbe->book1D("rctInputTPGEcalRank", "rctInputTPGEcalRank", TPGRANK, TPGRANKMIN, TPGRANKMAX) ;
rctInputTPGHcalOcc_ =
dbe->book2D("rctInputTPGHcalOcc", "rctInputTPGHcalOcc", TPGETABINS, TPGETAMIN,
TPGETAMAX, TPGPHIBINS, TPGPHIMIN, TPGPHIMAX);
rctInputTPGHcalSample_ =
dbe->book1D("rctInputTPGHcalSample", "rctInputTPGHcalSample", 10, -0.5, 9.5) ;
rctInputTPGHcalRank_ =
dbe->book1D("rctInputTPGHcalRank", "rctInputTPGHcalRank", TPGRANK, TPGRANKMIN, TPGRANKMAX) ;
dbe->setCurrentFolder(histFolder_+"EffCurves/NisoEm/");
trigEffThresh_ =
dbe->book2D("trigEffThresh", "Rank occupancy >= 2x trig thresh",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
dbe->setCurrentFolder(histFolder_+"EffCurves/NisoEm/ServiceData");
trigEffThreshOcc_ =
dbe->book2D("trigEffThreshOcc", "Rank occupancy >= 2x trig thresh",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
trigEffTriggThreshOcc_ =
dbe->book2D("trigEffTriggThreshOcc", "Rank occupancy >= 2x trig thresh, triggered",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
dbe->setCurrentFolder(histFolder_+"IsoEm");
rctIsoEmEff1_ =
dbe->book2D("rctIsoEmEff1", "rctIsoEmEff1", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmEff1oneD_ =
dbe->book1D("rctIsoEmEff1oneD", "rctIsoEmEff1oneD",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmEff2_ =
dbe->book2D("rctIsoEmEff2", "rctIsoEmEff2, energy matching required", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmEff2oneD_ =
dbe->book1D("rctIsoEmEff2oneD", "rctIsoEmEff2oneD, energy matching required",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmIneff2_ =
dbe->book2D("rctIsoEmIneff2", "rctIsoEmIneff2, energy matching required", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmIneff2oneD_ =
dbe->book1D("rctIsoEmIneff2oneD", "rctIsoEmIneff2oneD, energy matching required",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmIneff_ =
dbe->book2D("rctIsoEmIneff", "rctIsoEmIneff", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmIneff1D_ =
dbe->book1D("rctIsoEmIneff1D", "rctIsoEmIneff1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmOvereff_ =
dbe->book2D("rctIsoEmOvereff", "rctIsoEmOvereff", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmOvereff1D_ =
dbe->book1D("rctIsoEmOvereff1D", "rctIsoEmOvereff1D",
CHNLBINS, CHNLMIN, CHNLMAX);
dbe->setCurrentFolder(histFolder_+"IsoEm/ServiceData");
rctIsoEmDataOcc_ =
dbe->book2D("rctIsoEmDataOcc", "rctIsoEmDataOcc", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmDataOcc1D_ =
dbe->book1D("rctIsoEmDataOcc1D", "rctIsoEmDataOcc1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmEmulOcc_ =
dbe->book2D("rctIsoEmEmulOcc", "rctIsoEmEmulOcc", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmEmulOcc1D_ =
dbe->book1D("rctIsoEmEmulOcc1D", "rctIsoEmEmulOcc1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmEff1Occ_ =
dbe->book2D("rctIsoEmEff1Occ", "rctIsoEmEff1Occ", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmEff1Occ1D_ =
dbe->book1D("rctIsoEmEff1Occ1D", "rctIsoEmEff1Occ1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmEff2Occ_ =
dbe->book2D("rctIsoEmEff2Occ", "rctIsoEmEff2Occ", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmEff2Occ1D_ =
dbe->book1D("rctIsoEmEff2Occ1D", "rctIsoEmEff2Occ1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmIneff2Occ_ =
dbe->book2D("rctIsoEmIneff2Occ", "rctIsoEmIneff2Occ", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmIneff2Occ1D_ =
dbe->book1D("rctIsoEmIneff2Occ1D", "rctIsoEmIneff2Occ1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmIneffOcc_ =
dbe->book2D("rctIsoEmIneffOcc", "rctIsoEmIneffOcc", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmIneffOcc1D_ =
dbe->book1D("rctIsoEmIneffOcc1D", "rctIsoEmIneffOcc1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctIsoEmOvereffOcc_ =
dbe->book2D("rctIsoEmOvereffOcc", "rctIsoEmOvereffOcc", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctIsoEmOvereffOcc1D_ =
dbe->book1D("rctIsoEmOvereffOcc1D", "rctIsoEmOvereffOcc1D",
CHNLBINS, CHNLMIN, CHNLMAX);
dbe->setCurrentFolder(histFolder_+"NisoEm");
rctNisoEmEff1_ =
dbe->book2D("rctNisoEmEff1", "rctNisoEmEff1", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmEff1oneD_ =
dbe->book1D("rctNisoEmEff1oneD", "rctNisoEmEff1oneD",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmEff2_ =
dbe->book2D("rctNisoEmEff2", "rctNisoEmEff2, energy matching required", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmEff2oneD_ =
dbe->book1D("rctNisoEmEff2oneD", "rctNisoEmEff2oneD, energy matching required",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmIneff2_ =
dbe->book2D("rctNisoEmIneff2", "rctNisoEmIneff2, energy matching required", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmIneff2oneD_ =
dbe->book1D("rctNisoEmIneff2oneD", "rctNisoEmIneff2oneD, energy matching required",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmIneff_ =
dbe->book2D("rctNisoEmIneff", "rctNisoEmIneff", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmIneff1D_ =
dbe->book1D("rctNisoEmIneff1D", "rctNisoEmIneff1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmOvereff_ =
dbe->book2D("rctNisoEmOvereff", "rctNisoEmOvereff", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmOvereff1D_ =
dbe->book1D("rctNisoEmOvereff1D", "rctNisoEmOvereff1D",
CHNLBINS, CHNLMIN, CHNLMAX);
dbe->setCurrentFolder(histFolder_+"NisoEm/ServiceData");
rctNisoEmDataOcc_ =
dbe->book2D("rctNisoEmDataOcc", "rctNisoEmDataOcc", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmDataOcc1D_ =
dbe->book1D("rctNisoEmDataOcc1D", "rctNisoEmDataOcc1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmEmulOcc_ =
dbe->book2D("rctNisoEmEmulOcc", "rctNisoEmEmulOcc", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmEmulOcc1D_ =
dbe->book1D("rctNisoEmEmulOcc1D", "rctNisoEmEmulOcc1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmEff1Occ_ =
dbe->book2D("rctNisoEmEff1Occ", "rctNisoEmEff1Occ", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmEff1Occ1D_ =
dbe->book1D("rctNisoEmEff1Occ1D", "rctNisoEmEff1Occ1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmEff2Occ_ =
dbe->book2D("rctNisoEmEff2Occ", "rctNisoEmEff2Occ", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmEff2Occ1D_ =
dbe->book1D("rctNisoEmEff2Occ1D", "rctNisoEmEff2Occ1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmIneff2Occ_ =
dbe->book2D("rctNisoEmIneff2Occ", "rctNisoEmIneff2Occ", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmIneff2Occ1D_ =
dbe->book1D("rctNisoEmIneff2Occ1D", "rctNisoEmIneff2Occ1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmIneffOcc_ =
dbe->book2D("rctNisoEmIneffOcc", "rctNisoEmIneffOcc", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmIneffOcc1D_ =
dbe->book1D("rctNisoEmIneffOcc1D", "rctNisoEmIneffOcc1D",
CHNLBINS, CHNLMIN, CHNLMAX);
rctNisoEmOvereffOcc_ =
dbe->book2D("rctNisoEmOvereffOcc", "rctNisoEmOvereffOcc", ETABINS, ETAMIN,
ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctNisoEmOvereffOcc1D_ =
dbe->book1D("rctNisoEmOvereffOcc1D", "rctNisoEmOvereffOcc1D",
CHNLBINS, CHNLMIN, CHNLMAX);
// region information
dbe->setCurrentFolder(histFolder_+"RegionData");
rctRegEff1D_ =
dbe->book1D("rctRegEff1D", "1D region efficiency",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegIneff1D_ =
dbe->book1D("rctRegIneff1D", "1D region inefficiency",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegOvereff1D_ =
dbe->book1D("rctRegOvereff1D", "1D region overefficiency",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegSpEff1D_ =
dbe->book1D("rctRegSpEff1D", "1D region efficiency, energy matching required",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegSpIneff1D_ =
dbe->book1D("rctRegSpIneff1D", "1D region inefficiency, energy matching required",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegEff2D_ =
dbe->book2D("rctRegEff2D", "2D region efficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegIneff2D_ =
dbe->book2D("rctRegIneff2D", "2D region inefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegOvereff2D_ =
dbe->book2D("rctRegOvereff2D", "2D region overefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegSpEff2D_ =
dbe->book2D("rctRegSpEff2D", "2D region efficiency, energy matching required",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegSpIneff2D_ =
dbe->book2D("rctRegSpIneff2D", "2D region inefficiency, energy matching required",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
dbe->setCurrentFolder(histFolder_+"RegionData/ServiceData");
rctRegDataOcc1D_ =
dbe->book1D("rctRegDataOcc1D", "1D region occupancy from data",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegEmulOcc1D_ =
dbe->book1D("rctRegEmulOcc1D", "1D region occupancy from emulator",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegMatchedOcc1D_ =
dbe->book1D("rctRegMatchedOcc1D", "1D region occupancy for matched hits",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegUnmatchedDataOcc1D_ =
dbe->book1D("rctRegUnmatchedDataOcc1D", "1D region occupancy for unmatched hardware hits",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegUnmatchedEmulOcc1D_ =
dbe->book1D("rctRegUnmatchedEmulOcc1D", "1D region occupancy for unmatched emulator hits",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegSpEffOcc1D_ =
dbe->book1D("rctRegSpEffOcc1D", "1D region occupancy for \\Delta E_{T} efficiency",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegSpIneffOcc1D_ =
dbe->book1D("rctRegSpIneffOcc1D", "1D region occupancy for \\Delta E_{T} efficiency ",
CHNLBINS, CHNLMIN, CHNLMAX);
rctRegDataOcc2D_ =
dbe->book2D("rctRegDataOcc2D", "2D region occupancy from hardware",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegEmulOcc2D_ =
dbe->book2D("rctRegEmulOcc2D", "2D region occupancy from emulator",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegMatchedOcc2D_ =
dbe->book2D("rctRegMatchedOcc2D", "2D region occupancy for matched hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegUnmatchedDataOcc2D_ =
dbe->book2D("rctRegUnmatchedDataOcc2D", "2D region occupancy for unmatched hardware hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegUnmatchedEmulOcc2D_ =
dbe->book2D("rctRegUnmatchedEmulOcc2D", "2D region occupancy for unmatched emulator hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
// rctRegDeltaEt2D_ =
// dbe->book2D("rctRegDeltaEt2D", " \\Delta E_{T} for each channel",
// CHNLBINS, CHNLMIN, CHNLMAX, 100, -50., 50.);
rctRegSpEffOcc2D_ =
dbe->book2D("rctRegSpEffOcc2D", "2D region occupancy for \\Delta E_{T} efficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctRegSpIneffOcc2D_ =
dbe->book2D("rctRegSpIneffOcc2D", "2D region occupancy for \\Delta E_{T} inefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
// bit information
dbe->setCurrentFolder(histFolder_+"BitData");
rctBitOverFlowEff2D_ =
dbe->book2D("rctBitOverFlowEff2D", "2D overflow bit efficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitOverFlowIneff2D_ =
dbe->book2D("rctBitOverFlowIneff2D", "2D overflow bit inefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitOverFlowOvereff2D_ =
dbe->book2D("rctBitOverFlowOvereff2D", "2D overflow bit overefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitTauVetoEff2D_ =
dbe->book2D("rctBitTauVetoEff2D", "2D tau veto bit efficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitTauVetoIneff2D_ =
dbe->book2D("rctBitTauVetoIneff2D", "2D tau veto bit inefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitTauVetoOvereff2D_ =
dbe->book2D("rctBitTauVetoOvereff2D", "2D tau veto bit overefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitMipEff2D_ =
dbe->book2D("rctBitMipEff2D", "2D mip bit efficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitMipIneff2D_ =
dbe->book2D("rctBitMipIneff2D", "2D mip bit inefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitMipOvereff2D_ =
dbe->book2D("rctBitMipOvereff2D", "2D mip bit overefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
// QUIETBIT: To add quiet bit information, uncomment following 11 lines:
// rctBitQuietEff2D_ =
// dbe->book2D("rctBitQuietEff2D", "2D quiet bit efficiency",
// ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
// rctBitQuietIneff2D_ =
// dbe->book2D("rctBitQuietIneff2D", "2D quiet bit inefficiency",
// ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
// rctBitQuietOvereff2D_ =
// dbe->book2D("rctBitQuietOvereff2D", "2D quiet bit overefficiency",
// ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitHfPlusTauEff2D_ =
dbe->book2D("rctBitHfPlusTauEff2D", "2D HfPlusTau bit efficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitHfPlusTauIneff2D_ =
dbe->book2D("rctBitHfPlusTauIneff2D", "2D HfPlusTau bit inefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitHfPlusTauOvereff2D_ =
dbe->book2D("rctBitHfPlusTauOvereff2D", "2D HfPlusTau bit overefficiency",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
dbe->setCurrentFolder(histFolder_+"BitData/ServiceData");
rctBitEmulOverFlow2D_ =
dbe->book2D("rctBitEmulOverFlow2D", "2D overflow bit from emulator",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitDataOverFlow2D_ =
dbe->book2D("rctBitDataOverFlow2D", "2D overflow bit from hardware",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitMatchedOverFlow2D_ =
dbe->book2D("rctBitMatchedOverFlow2D", "2D overflow bit for matched hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedEmulOverFlow2D_ =
dbe->book2D("rctBitUnmatchedEmulOverFlow2D", "2D overflow bit for unmatched emulator hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedDataOverFlow2D_ =
dbe->book2D("rctBitUnmatchedDataOverFlow2D", "2D overflow bit for unmatched hardware hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitEmulTauVeto2D_ =
dbe->book2D("rctBitEmulTauVeto2D", "2D tau veto bit from emulator",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitDataTauVeto2D_ =
dbe->book2D("rctBitDataTauVeto2D", "2D tau veto bit from hardware",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitMatchedTauVeto2D_ =
dbe->book2D("rctBitMatchedTauVeto2D", "2D tau veto bit for matched hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedEmulTauVeto2D_ =
dbe->book2D("rctBitUnmatchedEmulTauVeto2D", "2D tau veto bit for unmatched emulator hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedDataTauVeto2D_ =
dbe->book2D("rctBitUnmatchedDataTauVeto2D", "2D tau veto bit for unmatched hardware hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitEmulMip2D_ =
dbe->book2D("rctBitEmulMip2D", "2D mip bit from emulator",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitDataMip2D_ =
dbe->book2D("rctBitDataMip2D", "2D mip bit from hardware",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitMatchedMip2D_ =
dbe->book2D("rctBitMatchedMip2D", "2D mip bit for matched hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedEmulMip2D_ =
dbe->book2D("rctBitUnmatchedEmulMip2D", "2D mip bit for unmatched emulator hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedDataMip2D_ =
dbe->book2D("rctBitUnmatchedDataMip2D", "2D mip bit for unmatched hardware hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitEmulQuiet2D_ =
dbe->book2D("rctBitEmulQuiet2D", "2D quiet bit from emulator",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitDataQuiet2D_ =
dbe->book2D("rctBitDataQuiet2D", "2D quiet bit from hardware",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitMatchedQuiet2D_ =
dbe->book2D("rctBitMatchedQuiet2D", "2D quiet bit for matched hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedEmulQuiet2D_ =
dbe->book2D("rctBitUnmatchedEmulQuiet2D", "2D quiet bit for unmatched emulator hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedDataQuiet2D_ =
dbe->book2D("rctBitUnmatchedDataQuiet2D", "2D quiet bit for unmatched hardware hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitEmulHfPlusTau2D_ =
dbe->book2D("rctBitEmulHfPlusTau2D", "2D HfPlusTau bit from emulator",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitDataHfPlusTau2D_ =
dbe->book2D("rctBitDataHfPlusTau2D", "2D HfPlusTau bit from hardware",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitMatchedHfPlusTau2D_ =
dbe->book2D("rctBitMatchedHfPlusTau2D", "2D HfPlusTau bit for matched hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedEmulHfPlusTau2D_ =
dbe->book2D("rctBitUnmatchedEmulHfPlusTau2D", "2D HfPlusTau bit for unmatched emulator hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
rctBitUnmatchedDataHfPlusTau2D_ =
dbe->book2D("rctBitUnmatchedDataHfPlusTau2D", "2D HfPlusTau bit for unmatched hardware hits",
ETABINS, ETAMIN, ETAMAX, PHIBINS, PHIMIN, PHIMAX);
dbe->setCurrentFolder(histFolder_+"/DBData");
fedVectorMonitorRUN_ = dbe->book2D("rctFedVectorMonitorRUN", "FED Vector Monitor Per Run",90,0,90,2,0,2);
fedVectorMonitorLS_ = dbe->book2D("rctFedVectorMonitorLS", "FED Vector Monitor Per LS",90,0,90,2,0,2);
for(unsigned int i=0;i<90;++i) {
char fed[10];
sprintf(fed,"%d",crateFED[i]);
fedVectorMonitorRUN_->getTH2F()->GetXaxis()->SetBinLabel(i+1,fed);
fedVectorMonitorLS_->getTH2F()->GetXaxis()->SetBinLabel(i+1,fed);
}
fedVectorMonitorRUN_->getTH2F()->GetYaxis()->SetBinLabel(1,"OUT");
fedVectorMonitorRUN_->getTH2F()->GetYaxis()->SetBinLabel(2,"IN");
fedVectorMonitorLS_->getTH2F()->GetYaxis()->SetBinLabel(1,"OUT");
fedVectorMonitorLS_->getTH2F()->GetYaxis()->SetBinLabel(2,"IN");
// for single channels
if(singlechannelhistos_)
{
for(int m=0; m<12; m++)
{
if(m==0) dbe->setCurrentFolder(histFolder_+"IsoEm/ServiceData/Eff1SnglChnls");
if(m==1) dbe->setCurrentFolder(histFolder_+"NisoEm/ServiceData/Eff1SnglChnls");
if(m==2) dbe->setCurrentFolder(histFolder_+"RegionData/ServiceData/EffSnglChnls");
if(m==3) dbe->setCurrentFolder(histFolder_+"IsoEm/ServiceData/IneffSnglChnls");
if(m==4) dbe->setCurrentFolder(histFolder_+"NisoEm/ServiceData/IneffSnglChnls");
if(m==5) dbe->setCurrentFolder(histFolder_+"RegionData/ServiceData/IneffSnglChnls");
if(m==6) dbe->setCurrentFolder(histFolder_+"IsoEm/ServiceData/OvereffSnglChnls");
if(m==7) dbe->setCurrentFolder(histFolder_+"NisoEm/ServiceData/OvereffSnglChnls");
if(m==8) dbe->setCurrentFolder(histFolder_+"RegionData/ServiceData/OvereffSnglChnls");
if(m==9) dbe->setCurrentFolder(histFolder_+"EffCurves/NisoEm/ServiceData/SingleChannels");
if(m==10) dbe->setCurrentFolder(histFolder_+"EffCurves/NisoEm/ServiceData/SingleChannels");
if(m==11) dbe->setCurrentFolder(histFolder_+"EffCurves/NisoEm/ServiceData/SingleChannels");
for(int i=0; i<ETAMAX; i++)
{
for(int j=0; j<PHIMAX; j++)
{
char name[80], channel[80]={""} ;
if(m==0) strcpy(name,"(Eemul-Edata)Chnl") ;
if(m==1) strcpy(name,"(Eemul-Edata)Chnl") ;
if(m==2) strcpy(name,"(Eemul-Edata)Chnl") ;
if(m==3) strcpy(name,"EemulChnl") ;
if(m==4) strcpy(name,"EemulChnl") ;
if(m==5) strcpy(name,"EemulChnl") ;
if(m==6) strcpy(name,"EdataChnl") ;
if(m==7) strcpy(name,"EdataChnl") ;
if(m==8) strcpy(name,"EdataChnl") ;
if(m==9) strcpy(name,"EemulChnlEff") ;
if(m==10) strcpy(name,"EemulChnlTrig") ;
if(m==11) strcpy(name,"EemulChnl") ;
if(i<10 && j<10) sprintf(channel,"_0%d0%d",i,j);
else if(i<10) sprintf(channel,"_0%d%d",i,j);
else if(j<10) sprintf(channel,"_%d0%d",i,j);
else sprintf(channel,"_%d%d",i,j);
strcat(name,channel);
int chnl=PHIBINS*i+j;
if(m==0) rctIsoEffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==1) rctNisoEffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==2) rctRegEffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==3) rctIsoIneffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==4) rctNisoIneffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==5) rctRegIneffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==6) rctIsoOvereffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==7) rctNisoOvereffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==8) rctRegOvereffChannel_[chnl] =
dbe->book1D(name, name, DEBINS, DEMIN, DEMAX);
if(m==9) trigEff_[chnl] =
dbe->book1D(name, name, ELBINS, ELMIN, ELMAX);
if(m==10) trigEffOcc_[chnl] =
dbe->book1D(name, name, ELBINS, ELMIN, ELMAX);
if(m==11) trigEffTriggOcc_[chnl] =
dbe->book1D(name, name, ELBINS, ELMIN, ELMAX);
}
}
}
}
//end of single channels
}
notrigCount=0;
trigCount=0;
}
| void L1TdeRCT::beginLuminosityBlock | ( | const edm::LuminosityBlock & | ls, |
| const edm::EventSetup & | es | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 1816 of file L1TdeRCT.cc.
References fedVectorMonitorLS_, and readFEDVector().
{
readFEDVector(fedVectorMonitorLS_,es);
}
| void L1TdeRCT::beginRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | es | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 1810 of file L1TdeRCT.cc.
References fedVectorMonitorRUN_, and readFEDVector().
{
readFEDVector(fedVectorMonitorRUN_,es);
}
| void L1TdeRCT::DivideME1D | ( | MonitorElement * | numerator, |
| MonitorElement * | denominator, | ||
| MonitorElement * | result | ||
| ) | [protected] |
Definition at line 1799 of file L1TdeRCT.cc.
References MonitorElement::getTH1F().
Referenced by analyze().
| void L1TdeRCT::DivideME2D | ( | MonitorElement * | numerator, |
| MonitorElement * | denominator, | ||
| MonitorElement * | result | ||
| ) | [protected] |
Definition at line 1789 of file L1TdeRCT.cc.
References MonitorElement::getTH2F().
Referenced by analyze().
| void L1TdeRCT::endJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 794 of file L1TdeRCT.cc.
References gather_cfg::cout, dbe, nev_, outputFile_, DQMStore::save(), and verbose_.
{
if (verbose_)
std::cout << "L1TdeRCT: end job...." << std::endl;
LogInfo("EndJob") << "analyzed " << nev_ << " events";
if (outputFile_.size() != 0 && dbe)
dbe->save(outputFile_);
//std::cout << "trig count is " << trigCount << std::endl;
//std::cout << "no trig count is " << notrigCount << std::endl;
return;
}
| void L1TdeRCT::readFEDVector | ( | MonitorElement * | histogram, |
| const edm::EventSetup & | es | ||
| ) | [protected] |
Definition at line 1821 of file L1TdeRCT.cc.
References crateFED, spr::find(), edm::EventSetup::get(), i, RunInfo::m_fed_in, edm::ESHandle< T >::product(), MonitorElement::setBinContent(), and summarizeEdmComparisonLogfiles::summary.
Referenced by beginLuminosityBlock(), and beginRun().
{
// adding fed mask into channel mask
edm::ESHandle<RunInfo> sum;
es.get<RunInfoRcd>().get(sum);
const RunInfo* summary=sum.product();
std::vector<int> caloFeds; // pare down the feds to the intresting ones
const std::vector<int> Feds = summary->m_fed_in;
for(std::vector<int>::const_iterator cf = Feds.begin(); cf != Feds.end(); ++cf){
int fedNum = *cf;
if(fedNum > 600 && fedNum <724)
caloFeds.push_back(fedNum);
}
for(unsigned int i=0;i<90;++i) {
std::vector<int>::iterator fv = std::find(caloFeds.begin(),caloFeds.end(),crateFED[i]);
if(fv!=caloFeds.end()) {
histogram->setBinContent(i+1,2,1);
histogram->setBinContent(i+1,1,0);
}
else
{
histogram->setBinContent(i+1,2,0);
histogram->setBinContent(i+1,1,1);
}
}
}
const int L1TdeRCT::crateFED [static, private] |
{613, 614, 603, 702, 718,
611, 612, 602, 700, 718,
627, 610, 601,716, 722,
625, 626, 609, 714, 722,
623, 624, 608, 712, 722,
621, 622, 607, 710, 720,
619, 620, 606, 708, 720,
617, 618, 605, 706, 720,
615, 616, 604, 704, 718,
631, 632, 648, 703, 719,
629, 630, 647, 701, 719,
645, 628, 646, 717, 723,
643, 644, 654, 715, 723,
641, 642, 653, 713, 723,
639, 640, 652, 711, 721,
637, 638, 651, 709, 721,
635, 636, 650, 707, 721,
633, 634, 649, 705, 719
}
Definition at line 299 of file L1TdeRCT.h.
Referenced by beginJob(), and readFEDVector().
DQMStore* L1TdeRCT::dbe [private] |
Definition at line 134 of file L1TdeRCT.h.
Referenced by beginJob(), endJob(), and L1TdeRCT().
int L1TdeRCT::doubleThreshold_ [private] |
Definition at line 320 of file L1TdeRCT.h.
Referenced by analyze().
edm::InputTag L1TdeRCT::ecalTPGData_ [private] |
Definition at line 316 of file L1TdeRCT.h.
Referenced by analyze().
MonitorElement* L1TdeRCT::fedVectorMonitorLS_ [private] |
Definition at line 301 of file L1TdeRCT.h.
Referenced by beginJob(), and beginLuminosityBlock().
MonitorElement* L1TdeRCT::fedVectorMonitorRUN_ [private] |
Definition at line 300 of file L1TdeRCT.h.
Referenced by beginJob(), and beginRun().
edm::InputTag L1TdeRCT::gtDigisLabel_ [private] |
Definition at line 318 of file L1TdeRCT.h.
Referenced by analyze().
std::string L1TdeRCT::gtEGAlgoName_ [private] |
Definition at line 319 of file L1TdeRCT.h.
Referenced by analyze().
edm::InputTag L1TdeRCT::hcalTPGData_ [private] |
Definition at line 317 of file L1TdeRCT.h.
Referenced by analyze().
std::string L1TdeRCT::histFolder_ [private] |
Definition at line 308 of file L1TdeRCT.h.
Referenced by beginJob(), and L1TdeRCT().
ofstream L1TdeRCT::logFile_ [private] |
Definition at line 312 of file L1TdeRCT.h.
bool L1TdeRCT::monitorDaemon_ [private] |
Definition at line 311 of file L1TdeRCT.h.
int L1TdeRCT::nev_ [private] |
Definition at line 306 of file L1TdeRCT.h.
Referenced by analyze(), beginJob(), and endJob().
int L1TdeRCT::notrigCount [private] |
Definition at line 323 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
std::string L1TdeRCT::outputFile_ [private] |
Definition at line 307 of file L1TdeRCT.h.
Referenced by endJob(), and L1TdeRCT().
MonitorElement* L1TdeRCT::rctBitDataHfPlusTau2D_ [private] |
Definition at line 206 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitDataMip2D_ [private] |
Definition at line 189 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitDataOverFlow2D_ [private] |
Definition at line 173 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitDataQuiet2D_ [private] |
Definition at line 197 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitDataTauVeto2D_ [private] |
Definition at line 181 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitEmulHfPlusTau2D_ [private] |
Definition at line 205 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitEmulMip2D_ [private] |
Definition at line 188 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitEmulOverFlow2D_ [private] |
Definition at line 172 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitEmulQuiet2D_ [private] |
Definition at line 196 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitEmulTauVeto2D_ [private] |
Definition at line 180 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitHfPlusTauEff2D_ [private] |
Definition at line 210 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitHfPlusTauIneff2D_ [private] |
Definition at line 211 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 212 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 207 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitMatchedMip2D_ [private] |
Definition at line 190 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitMatchedOverFlow2D_ [private] |
Definition at line 174 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitMatchedQuiet2D_ [private] |
Definition at line 198 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitMatchedTauVeto2D_ [private] |
Definition at line 182 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitMipEff2D_ [private] |
Definition at line 193 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitMipIneff2D_ [private] |
Definition at line 194 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitMipOvereff2D_ [private] |
Definition at line 195 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitOverFlowEff2D_ [private] |
Definition at line 177 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitOverFlowIneff2D_ [private] |
Definition at line 178 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitOverFlowOvereff2D_ [private] |
Definition at line 179 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitTauVetoEff2D_ [private] |
Definition at line 185 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitTauVetoIneff2D_ [private] |
Definition at line 186 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctBitTauVetoOvereff2D_ [private] |
Definition at line 187 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 209 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 192 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 176 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 200 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 184 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 208 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 191 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 175 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 199 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 183 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctInputTPGEcalOcc_ [private] |
Definition at line 216 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctInputTPGEcalRank_ [private] |
Definition at line 217 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctInputTPGHcalOcc_ [private] |
Definition at line 218 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctInputTPGHcalRank_ [private] |
Definition at line 219 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctInputTPGHcalSample_ [private] |
Definition at line 220 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEffChannel_[396] [private] |
Definition at line 274 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmDataOcc1D_ [private] |
Definition at line 235 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmDataOcc_ [private] |
Definition at line 222 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEff1_ [private] |
Definition at line 229 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEff1Occ1D_ [private] |
Definition at line 237 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEff1Occ_ [private] |
Definition at line 224 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEff1oneD_ [private] |
Definition at line 242 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEff2_ [private] |
Definition at line 230 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEff2Occ1D_ [private] |
Definition at line 238 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEff2Occ_ [private] |
Definition at line 225 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEff2oneD_ [private] |
Definition at line 243 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEmulOcc1D_ [private] |
Definition at line 236 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmEmulOcc_ [private] |
Definition at line 223 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmIneff1D_ [private] |
Definition at line 245 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmIneff2_ [private] |
Definition at line 231 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmIneff2Occ1D_ [private] |
Definition at line 239 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmIneff2Occ_ [private] |
Definition at line 226 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmIneff2oneD_ [private] |
Definition at line 244 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmIneff_ [private] |
Definition at line 232 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmIneffOcc1D_ [private] |
Definition at line 240 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmIneffOcc_ [private] |
Definition at line 227 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmOvereff1D_ [private] |
Definition at line 246 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmOvereff_ [private] |
Definition at line 233 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmOvereffOcc1D_ [private] |
Definition at line 241 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoEmOvereffOcc_ [private] |
Definition at line 228 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoIneffChannel_[396] [private] |
Definition at line 275 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctIsoOvereffChannel_[396] [private] |
Definition at line 276 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEffChannel_[396] [private] |
Definition at line 278 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmDataOcc1D_ [private] |
Definition at line 261 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmDataOcc_ [private] |
Definition at line 248 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEff1_ [private] |
Definition at line 255 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEff1Occ1D_ [private] |
Definition at line 263 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEff1Occ_ [private] |
Definition at line 250 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEff1oneD_ [private] |
Definition at line 268 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEff2_ [private] |
Definition at line 256 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEff2Occ1D_ [private] |
Definition at line 264 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEff2Occ_ [private] |
Definition at line 251 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEff2oneD_ [private] |
Definition at line 269 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEmulOcc1D_ [private] |
Definition at line 262 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmEmulOcc_ [private] |
Definition at line 249 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmIneff1D_ [private] |
Definition at line 271 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmIneff2_ [private] |
Definition at line 257 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmIneff2Occ1D_ [private] |
Definition at line 265 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmIneff2Occ_ [private] |
Definition at line 252 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmIneff2oneD_ [private] |
Definition at line 270 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmIneff_ [private] |
Definition at line 258 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmIneffOcc1D_ [private] |
Definition at line 266 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmIneffOcc_ [private] |
Definition at line 253 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmOvereff1D_ [private] |
Definition at line 272 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmOvereff_ [private] |
Definition at line 259 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmOvereffOcc1D_ [private] |
Definition at line 267 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoEmOvereffOcc_ [private] |
Definition at line 254 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoIneffChannel_[396] [private] |
Definition at line 279 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctNisoOvereffChannel_[396] [private] |
Definition at line 280 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegDataOcc1D_ [private] |
Definition at line 140 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegDataOcc2D_ [private] |
Definition at line 154 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegEff1D_ [private] |
Definition at line 148 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegEff2D_ [private] |
Definition at line 163 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegEffChannel_[396] [private] |
Definition at line 283 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegEmulOcc1D_ [private] |
Definition at line 141 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegEmulOcc2D_ [private] |
Definition at line 155 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegIneff1D_ [private] |
Definition at line 149 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegIneff2D_ [private] |
Definition at line 164 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegIneffChannel_[396] [private] |
Definition at line 284 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegMatchedOcc1D_ [private] |
Definition at line 142 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegMatchedOcc2D_ [private] |
Definition at line 156 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegOvereff1D_ [private] |
Definition at line 150 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegOvereff2D_ [private] |
Definition at line 165 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegOvereffChannel_[396] [private] |
Definition at line 285 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegSpEff1D_ [private] |
Definition at line 151 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegSpEff2D_ [private] |
Definition at line 166 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegSpEffOcc1D_ [private] |
Definition at line 145 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegSpEffOcc2D_ [private] |
Definition at line 160 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegSpIneff1D_ [private] |
Definition at line 152 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegSpIneff2D_ [private] |
Definition at line 167 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegSpIneffOcc1D_ [private] |
Definition at line 146 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::rctRegSpIneffOcc2D_ [private] |
Definition at line 161 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 143 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 157 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 144 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
Definition at line 158 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
edm::InputTag L1TdeRCT::rctSourceData_ [private] |
Definition at line 315 of file L1TdeRCT.h.
Referenced by analyze().
edm::InputTag L1TdeRCT::rctSourceEmul_ [private] |
Definition at line 314 of file L1TdeRCT.h.
Referenced by analyze().
bool L1TdeRCT::singlechannelhistos_ [private] |
Definition at line 310 of file L1TdeRCT.h.
Referenced by analyze(), beginJob(), and L1TdeRCT().
int L1TdeRCT::trigCount [private] |
Definition at line 323 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::trigEff_[396] [private] |
Definition at line 291 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::trigEffOcc_[396] [private] |
Definition at line 292 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::trigEffThresh_ [private] |
Definition at line 288 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::trigEffThreshOcc_ [private] |
Definition at line 289 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::trigEffTriggOcc_[396] [private] |
Definition at line 293 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::trigEffTriggThreshOcc_ [private] |
Definition at line 290 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
MonitorElement* L1TdeRCT::triggerAlgoNumbers_ [private] |
Definition at line 137 of file L1TdeRCT.h.
Referenced by analyze(), and beginJob().
bool L1TdeRCT::verbose_ [private] |
Definition at line 309 of file L1TdeRCT.h.
Referenced by analyze(), endJob(), and L1TdeRCT().
1.7.3