L1TdeGCT Class Reference

#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 ()

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
CurrentProcessingContext const *	currentContext () const

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, 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.

{}

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, scaleCards::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(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), 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;
}

int L1TdeGCT::verbose ( )

inlineprivate

Definition at line 48 of file L1TdeGCT.h.

References verbose_.

{return verbose_;}

Member Data Documentation

int L1TdeGCT::colCount[nGctColl_]

private

Definition at line 64 of file L1TdeGCT.h.

DQMStore* L1TdeGCT::dbe

private

Definition at line 57 of file L1TdeGCT.h.

edm::InputTag L1TdeGCT::DEsource_

private

Definition at line 43 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::deword[nGctColl_]

private

Definition at line 83 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::dword[nGctColl_]

private

Definition at line 81 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::errortype[nGctColl_]

private

Definition at line 70 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::eta[nGctColl_]

private

Definition at line 73 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::etaData[nGctColl_]

private

Definition at line 76 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::etaphi[nGctColl_]

private

Definition at line 72 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::eword[nGctColl_]

private

Definition at line 82 of file L1TdeGCT.h.

bool L1TdeGCT::hasRecord_

private

Definition at line 44 of file L1TdeGCT.h.

std::string L1TdeGCT::histFile_

private

Definition at line 51 of file L1TdeGCT.h.

std::string L1TdeGCT::histFolder_

private

Definition at line 54 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::masked[nGctColl_]

private

Definition at line 84 of file L1TdeGCT.h.

bool L1TdeGCT::monitorDaemon_

private

Definition at line 58 of file L1TdeGCT.h.

const int L1TdeGCT::nGctColl_ = dedefs::GCThfbit-dedefs::GCTisolaem+1

staticprivate

Definition at line 61 of file L1TdeGCT.h.

int L1TdeGCT::nWithCol[nGctColl_]

private

Definition at line 65 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::phi[nGctColl_]

private

Definition at line 74 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::phiData[nGctColl_]

private

Definition at line 77 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::rnk[nGctColl_]

private

Definition at line 75 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::rnkData[nGctColl_]

private

Definition at line 78 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::sysncand[2]

private

Definition at line 69 of file L1TdeGCT.h.

MonitorElement* L1TdeGCT::sysrates

private

Definition at line 68 of file L1TdeGCT.h.

int L1TdeGCT::verbose_

private

Definition at line 47 of file L1TdeGCT.h.

Referenced by verbose().