#include <L1TdeGCT.h>

Inheritance diagram for L1TdeGCT:

Public Member Functions
	L1TdeGCT (const edm::ParameterSet &)

	~L1TdeGCT ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob (void)

virtual void	endJob ()

Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Member Functions
int	verbose ()

Private Attributes
int	colCount [nGctColl_]

DQMStore *	dbe

edm::InputTag	DEsource_

MonitorElement *	deword [nGctColl_]

MonitorElement *	dword [nGctColl_]

MonitorElement *	errortype [nGctColl_]

MonitorElement *	eta [nGctColl_]

MonitorElement *	etaData [nGctColl_]

MonitorElement *	etaphi [nGctColl_]

MonitorElement *	eword [nGctColl_]

bool	hasRecord_

std::string	histFile_

std::string	histFolder_

MonitorElement *	masked [nGctColl_]

bool	monitorDaemon_

int	nWithCol [nGctColl_]

MonitorElement *	phi [nGctColl_]

MonitorElement *	phiData [nGctColl_]

MonitorElement *	rnk [nGctColl_]

MonitorElement *	rnkData [nGctColl_]

MonitorElement *	sysncand [2]

MonitorElement *	sysrates

int	verbose_

Static Private Attributes
static const int	nGctColl_ = dedefs::GCThfbit-dedefs::GCTisolaem+1

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 27 of file L1TdeGCT.h.

Constructor & Destructor Documentation

L1TdeGCT::L1TdeGCT ( const edm::ParameterSet & iConfig )

explicit

Definition at line 6 of file L1TdeGCT.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), NULL, cppFunctionSkipper::operator, AlCaHLTBitMon_QueryRunRegistry::string, and validate_alignment_devdb10_cfg::verbose.

                                                  {
 
   verbose_ = iConfig.getUntrackedParameter<int>("VerboseFlag",0);
 
   if(verbose())
     std::cout << "L1TdeGCT::L1TdeGCT()...\n" << std::flush;
   
   DEsource_ = iConfig.getParameter<edm::InputTag>("DataEmulCompareSource");
   histFolder_ = iConfig.getUntrackedParameter<std::string>("HistFolder", "L1TEMU/GCTexpert");
   
   dbe = NULL;
   if (iConfig.getUntrackedParameter<bool>("DQMStore", false)) { 
     dbe = edm::Service<DQMStore>().operator->();
     dbe->setVerbose(0);
   }
   
   histFile_ = iConfig.getUntrackedParameter<std::string>("HistFile", "");
   if(iConfig.getUntrackedParameter<bool> ("disableROOToutput", true))
     histFile_ = "";
 
   if (histFile_.size()!=0) {
     edm::LogInfo("OutputRootFile") 
       << "L1TEmulator GCT specific histograms will be saved to " 
       << histFile_.c_str() 
       << std::endl;
   }
 
   if(dbe!=NULL)
     dbe->setCurrentFolder(histFolder_);
 
   hasRecord_=true;
   
   if(verbose())
     std::cout << "L1TdeGCT::L1TdeGCT()...done.\n" << std::flush;
 }

L1TdeGCT::~L1TdeGCT ( )

Definition at line 42 of file L1TdeGCT.cc.

42 {}

Member Function Documentation

void L1TdeGCT::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

protectedvirtual

get the comparison results

get the de candidates

— Fill histograms(me) —

(note see L1Trigger/HardwareValidation/interface/DEtrait.h)

bitset loop

event based rate

object based rate

Implements edm::EDAnalyzer.

Definition at line 219 of file L1TdeGCT.cc.

References edm::HandleBase::clear(), gather_cfg::cout, eta(), fastmath::etaphi(), HcalObjRepresent::Fill(), dedefs::GCT, dedefs::GCTisolaem, edm::Event::getByLabel(), i, edm::HandleBase::isValid(), j, LogDebug, phi, RPCpg::rate(), L1DataEmulDigi::reset(), evf::utils::sid, dedefs::SystLabel, and validate_alignment_devdb10_cfg::verbose.

                                                                        {
   
   if(!hasRecord_)
     return;
   
   if(verbose())
     std::cout << "L1TdeGCT::analyze()  start\n" << std::flush;
 
   edm::Handle<L1DataEmulRecord> deRecord;
   iEvent.getByLabel(DEsource_, deRecord);
 
   if (!deRecord.isValid()) {
     edm::LogInfo("DataNotFound") 
       << "Cannot find L1DataEmulRecord with label "
       << DEsource_.label() 
       << " Please verify that comparator was successfully executed."
       << " Emulator DQM for GCT will be skipped!"
       << std::endl;
     hasRecord_=false;
     return;
   }
 
   bool isComp = deRecord->get_isComp(GCT);
   if(!isComp) {
     if(verbose()) 
       std::cout << "[L1TdeGCT] Gct information not generated in de-record."
         << " Skiping event!\n" << std::flush;
     return;
   }
 
   int DEncand[2]={0};
   for(int j=0; j<2; j++) 
     DEncand[j] = deRecord->getNCand(GCT,j);
   
   if(verbose()) 
     std::cout << "[L1TdeGCT] ncands" 
           << " data: " << DEncand[0]
           << " emul: " << DEncand[1]
           << std::endl;
 
 
   L1DEDigiCollection deColl;
   deColl = deRecord->getColl();
 
   // extract the GCT comparison digis
   L1DEDigiCollection gctColl;
 
   gctColl.reserve(20);
   gctColl.clear();
 
 
   for(L1DEDigiCollection::const_iterator it=deColl.begin(); 
       it!=deColl.end(); it++) 
     if(!it->empty()) 
       if(it->sid()==GCT)
     gctColl.push_back(*it);
   
   if(verbose()) {
     std::cout << "[L1TdeGCT] record has " << gctColl.size() 
           << " gct de digis\n" << std::flush;
     for(L1DEDigiCollection::const_iterator it=gctColl.begin(); 
     it!=gctColl.end(); it++)
       std::cout << "\t" << *it << std::endl;
   }
 
   const int nullVal = L1DataEmulDigi().reset();
 
   
   // d|e candidate loop
   for(L1DEDigiCollection::const_iterator it=gctColl.begin(); 
       it!=gctColl.end(); it++) {
     
     // sid should be GCT
     int sid = it->sid();
     // cid: GCTisolaem, GCTnoisoem, GCTcenjets, GCTforjets, GCTtaujets
     int cid = it->cid();
 
     if(verbose()) 
       std::cout << "[L1TdeGCT] processing digi "
         << " sys:"  << sid
         << " type:" << cid
         << " \n\t"
         << *it << "\n" 
         << std::flush;
     
     //assert(cid==GCT);
     if(sid!=GCT || it->empty()) {
       LogDebug("L1TdeGCT") << "consistency check failure, non-gct digis!";
       continue;
     }
     
     int type    = it->type();
     double phiv = it->x1();
     double etav = it->x2();
     float rankarr[2]; 
     it->rank(rankarr);
     float rnkv = rankarr[0];
 
     double wei = 1.;
 
     unsigned int mask = (~0x0);
 
     // shift coll type for starting at zero
     int ccid = cid - dedefs::GCTisolaem; 
     if(ccid<0 || ccid >= nGctColl_) {
       LogDebug("L1TdeGCT") << "consistency check failure, col type outbounds:"
                << ccid << "\n";
       ccid=0;
     }
     
     //type: 0:agree 1:loc.agree, 2:loc.disagree, 3:data.only, 4:emul.only
     if(it->type()<4) 
       sysncand[0]->Fill(ccid); 
     if(it->type()<5&&it->type()!=3) 
       sysncand[1]->Fill(ccid);
     
     errortype[ccid]->Fill(type);
 
     wei=1.; if(!type) wei=0.;
     if(etav!=nullVal && phiv!=nullVal)
       etaphi[ccid]->Fill(etav,phiv,wei);
     if(etav!=nullVal)
       eta   [ccid]->Fill(etav,wei);
     if(phiv!=nullVal)
       phi   [ccid]->Fill(phiv,wei);
     rnk[ccid]->Fill(rnkv,wei);    
 
     //exclude e-only cands (only data)
     wei=1.;if(type==4) wei=0.;
     if(etav!=nullVal)
       etaData[ccid]->Fill(etav,wei);
     if(phiv!=nullVal)
       phiData[ccid]->Fill(phiv,wei);
     rnkData[ccid]->Fill(rnkv,wei);
     wei=1;
 
     // GCT trigger bits
     unsigned int word[2];
     it->data(word);
     std::bitset<32> dbits(word[0]);
     std::bitset<32> ebits(word[1]);
     unsigned int dexor = ( (word[0]) ^ (word[1]) );
     //disagreeing bits
     std::bitset<32> debits(dexor);
     //disagreeing bits after masking
     std::bitset<32> dembits( ( (dexor) & (mask) ) );
     
     if(verbose())
       std::cout << "l1degct" 
         << " sid:" << sid << " cid:" << cid << "\n"
         << " data:0x" << std::hex << word[0] << std::dec
         << " bitset:" << dbits
         << "\n"
         << " emul:0x" << std::hex << word[1] << std::dec
         << " bitset:" << ebits
         << "\n"
         << "  xor:0x" << std::hex << dexor << std::dec
         << " bitset:" << debits
         << " bitset:" << ( (dbits) ^ (ebits) )
         << "\n" << std::flush;
 
     for(int ibit=0; ibit<32; ibit++) {
       wei=1.;
       //comparison gives no info if there's only 1 candidate
       if(type==3 || type==4) wei=0.; 
       if(dbits  [ibit]) dword[sid]->Fill(ibit,wei);
       if(ebits  [ibit]) eword[sid]->Fill(ibit,wei);
       if(debits [ibit])deword[sid]->Fill(ibit,wei);
       //if(dembits[ibit])masked[sid]->Fill(ibit,wei);
     }
     wei=1;
     
   }
 
   //error rates per GCT trigger object type
   int hasCol[nGctColl_]={0};
   int nagree[nGctColl_]={0};
   for(L1DEDigiCollection::const_iterator it=gctColl.begin(); 
       it!=gctColl.end(); it++) {
     int ccid = it->cid()-dedefs::GCTisolaem;
     ccid = (ccid<0 || ccid >= nGctColl_) ? 0:ccid;
     hasCol[ccid]++;
     if(!it->type()) 
       nagree[ccid]++;
   }
   for(int i=0; i<nGctColl_; i++) {
     if(!hasCol[i]) continue;
     //nWithCol[i]++;
     //if(nagree[i]<hasCol[i]) colCount[i]++;
     nWithCol[i]+=hasCol[i];//#of objects
     colCount[i]+=nagree[i];//#of agreements
   }
   for(int i=0; i<nGctColl_; i++) {
     int ibin = i+1;
     double rate = nWithCol[i] ? 1.-1.*colCount[i]/nWithCol[i]: 0.;
     sysrates->setBinContent(ibin,rate);
     if(verbose()) {
       std::cout << "[L1TDEMON] analyze rate computation\t\n"
         << " colid:"   << i
         << "(so far)"
         << " nWithCol: " << nWithCol[i]
         << " colCount: " << colCount[i]
         << "(this event)"
         << "hasCol: "    << hasCol[i] 
         << " nagree: "   << nagree[i]
         << " rate:"    << sysrates->getBinContent(ibin) 
         << "\n" << std::flush;
       if(rate>1. || rate<0.)
     std::cout << "problem, error rate for " << SystLabel[i] 
           <<" is "<<sysrates->getBinContent(ibin)
           << "\n" << std::flush;
     }
   }
 
     
   if(verbose())
     std::cout << "L1TdeGCT::analyze() end.\n" << std::flush;
   
 }

void L1TdeGCT::beginJob ( void )

protectedvirtual

labeling

Reimplemented from edm::EDAnalyzer.

Definition at line 45 of file L1TdeGCT.cc.

References DQMStore::book1D(), DQMStore::book2D(), gather_cfg::cout, eta(), fastmath::etaphi(), i, j, cppFunctionSkipper::operator, phi, DQMStore::rmdir(), DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, indexGen::title, and validate_alignment_devdb10_cfg::verbose.

                        {
 
   if(verbose())
     std::cout << "L1TdeGCT::beginJob()  start\n" << std::flush;
 
   DQMStore* dbe = 0;
   dbe = edm::Service<DQMStore>().operator->();
   if(dbe) {
     dbe->setCurrentFolder(histFolder_);
     dbe->rmdir(histFolder_);
   }
 
   // (em) iso, no-iso, (jets) cen, for, tau
   std::string cLabel[nGctColl_]= 
     {"IsoEM", "NoisoEM", "CenJet", "ForJet", "TauJet", "HT", "MET", "ET", "MHT", "HFSums", "HFCnts"};
   const int nerr  = 5; 
   const int nbit = 32;
   
   if(dbe) {
     dbe->setCurrentFolder(histFolder_);
 
     // book histograms here 
 
     const int    phiNBins = 18  ;
     const double phiMinim = -0.5;
     const double phiMaxim = 17.5;
     const int    etaNBins = 22  ;
     const double etaMinim = -0.5;
     const double etaMaxim = 21.5;
     const int    rnkNBins = 63;
     const double rnkMinim = 0.5;
     const double rnkMaxim = 63.5;
 
     sysrates = dbe->book1D("sysrates","RATE OF COMPARISON FAILURES",nGctColl_, 0, nGctColl_ );
 
     for(int j=0; j<2; j++) {
       std::string lbl("sysncand"); 
       lbl += (j==0?"Data":"Emul");
       std::string title("GCT OBJECT MULTIPLICITY ");
       title += (j==0?"(DATA)":"(EMULATOR)");
       sysncand[j] = dbe->book1D(lbl.data(),title.data(),nGctColl_, 0, nGctColl_ );
     }
     
     for(int j=0; j<nGctColl_; j++) {
       
       dbe->setCurrentFolder(std::string(histFolder_+"/"+cLabel[j]));
       
       std::string lbl("");
       std::string title("");
       lbl.clear();
       title.clear();
       lbl+=cLabel[j];lbl+="ErrorFlag"; 
       title+=cLabel[j];title+=" ErrorFlag"; 
       errortype[j] = dbe->book1D(lbl.data(),title.data(), nerr, 0, nerr);
       
       lbl.clear();
       title.clear();
       lbl+=cLabel[j];lbl+="Eta"; 
       title+=cLabel[j];title+=" ETA OF COMPARISON FAILURES"; 
       eta[j] = dbe->book1D(lbl.data(),title.data(),
                etaNBins, etaMinim, etaMaxim);
       lbl.clear();
       title.clear();
       lbl+=cLabel[j];lbl+="Phi"; 
       title+=cLabel[j];title+=" PHI OF COMPARISON FAILURES"; 
       phi[j] = dbe->book1D(lbl.data(),title.data(),
                phiNBins, phiMinim, phiMaxim);
 
       lbl.clear();
       title.clear();
       lbl+=cLabel[j];lbl+="Etaphi"; 
       title+=cLabel[j];title+=" ETA PHI OF COMPARISON FAILURES"; 
       etaphi[j] = dbe->book2D(lbl.data(),title.data(), 
                   etaNBins, etaMinim, etaMaxim,
                   phiNBins, phiMinim, phiMaxim
                   );
       lbl.clear();
       title.clear();
       lbl+=cLabel[j];lbl+="Rank";
       title+=cLabel[j];title+=" RANK OF COMPARISON FAILURES"; 
       rnk[j] = dbe->book1D(lbl.data(),title.data(),
                    rnkNBins, rnkMinim, rnkMaxim);
       //
       lbl.clear();
       title.clear();
       lbl+=cLabel[j];lbl+="Eta"; lbl+="Data";
       title+=cLabel[j];title+=" ETA (DATA)"; 
       etaData[j] = dbe->book1D(lbl.data(),title.data(),
                    etaNBins, etaMinim, etaMaxim);
       lbl.clear();
       title.clear();
       lbl+=cLabel[j];lbl+="Phi";  lbl+="Data";
       title+=cLabel[j];title+=" PHI (DATA)"; 
       phiData[j] = dbe->book1D(lbl.data(),title.data(),
                    phiNBins, phiMinim, phiMaxim);
 
       lbl.clear();
       title.clear();
       lbl+=cLabel[j];lbl+="Rank";  lbl+="Data";
       title+=cLabel[j];title+=" RANK (DATA)"; 
       rnkData[j] = dbe->book1D(lbl.data(),title.data(),
                    rnkNBins, rnkMinim, rnkMaxim);
       lbl.clear();
       lbl+=cLabel[j];lbl+="Dword"; 
       dword[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
       lbl.clear();
       lbl+=cLabel[j];lbl+="Eword"; 
       eword[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
       lbl.clear();
       lbl+=cLabel[j];lbl+="DEword"; 
       deword[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
       //lbl.clear();
       //lbl+=cLabel[j];lbl+="Masked"; 
       //masked[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
     }
 
   }
   
   std::string errLabel[nerr]= {
     "Agree", "Loc. Agree", "L.Disagree", "Data only", "Emul only"
   };
   
   for(int i=0; i<nGctColl_; i++) {
     sysrates   ->setBinLabel(i+1,cLabel[i]);
     sysncand[0]->setBinLabel(i+1,cLabel[i]);
     sysncand[1]->setBinLabel(i+1,cLabel[i]);
   }
 
   for(int i=0; i<nGctColl_; i++) {
     for(int j=0; j<nerr; j++) {
       errortype[i]->setBinLabel(j+1,errLabel[j]);
     }
   }
   
   for(int i=0; i<nGctColl_; i++) {
     etaphi [i]->setAxisTitle("GCT #eta",1);
     etaphi [i]->setAxisTitle("GCT #phi",2);
     eta    [i]->setAxisTitle("GCT #eta");
     phi    [i]->setAxisTitle("GCT #phi");
     rnk    [i]->setAxisTitle("Rank");
     etaData[i]->setAxisTitle("GCT #eta");
     phiData[i]->setAxisTitle("GCT #phi");
     rnkData[i]->setAxisTitle("Rank");
     dword  [i]->setAxisTitle("trigger data word bit");
     eword  [i]->setAxisTitle("trigger data word bit");
     deword [i]->setAxisTitle("trigger data word bit");
     //masked [i]->setAxisTitle("trigger data word bit");
   }
 
   for(int i=0; i<nGctColl_; i++) {
     colCount[i]=0;
     nWithCol[i]=0;
   }
   
   if(verbose())
     std::cout << "L1TdeGCT::beginJob()  end.\n" << std::flush;
 }

void L1TdeGCT::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 205 of file L1TdeGCT.cc.

References gather_cfg::cout, and validate_alignment_devdb10_cfg::verbose.

                  {
   if(verbose())
     std::cout << "L1TdeGCT::endJob()...\n" << std::flush;
 
   if(histFile_.size()!=0  && dbe) 
     dbe->save(histFile_);
 
   if(verbose())
     std::cout << "L1TdeGCT::endJob()  end.\n" << std::flush;
 }