/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/DataFormats/Luminosity/src/LumiSummary.cc

Go to the documentation of this file.

// $Id: LumiSummary.cc,v 1.21 2010/10/12 10:54:59 xiezhen Exp $

#include "DataFormats/Luminosity/interface/LumiSummary.h"

#include <iomanip>
#include <ostream>
#include <iostream>
float
LumiSummary::avgInsDelLumi()const{ 
  return avginsdellumi_;
}
float
LumiSummary::intgDelLumi()const{
  return avginsdellumi_*float(this->lumiSectionLength());
}
float
LumiSummary::avgInsDelLumiErr()const{ 
  return  avginsdellumierr_;
}
float 
LumiSummary::intgRecLumi()const{
  return this->avgInsRecLumi() *float(this->lumiSectionLength());
}
short
LumiSummary::lumiSecQual()const {
  return lumisecqual_; 
}
unsigned long long
LumiSummary::deadcount() const{
  return deadcount_;
}
float 
LumiSummary::deadFrac() const {
  //definition: deadcount/bizerocount
  //if no trigger data, return deadfraction 1.0,mask out this LS
  //if bitzerocount=0, return -1.0 meaning no beam
  if (l1data_.size()==0) return 1.0;
  if (l1data_.begin()->ratecount==0) return -1.0;
  return float(deadcount_)/float(l1data_.begin()->ratecount*l1data_.begin()->prescale);
}
float 
LumiSummary::liveFrac() const { 
  //1-deadfraction
  //else if deadfraction<0 meaning no beam, live fraction=0
  //
  if (deadFrac()<0) return 0;
  return 1-deadFrac();
}
float
LumiSummary::lumiSectionLength() const {
  //numorbits*3564*25e-09
  return numorbit_*3564.0*25.0*10e-9;
}
unsigned int 
LumiSummary::lsNumber() const{
  return lsnumber_; 
}
unsigned int
LumiSummary::startOrbit() const{
  return startorbit_; 
}
unsigned int
LumiSummary::numOrbit() const{
  return numorbit_;
}
bool 
LumiSummary::isValid() const {
  return (lumiversion_!="-1"); 
}
LumiSummary::L1   
LumiSummary::l1info(const std::string& name) const{
  for(std::vector<L1>::const_iterator it=l1data_.begin();it!=l1data_.end();++it){
    if(it->triggername==name) return *it;
  }
  return LumiSummary::L1();
}
LumiSummary::L1  
LumiSummary::l1info(unsigned int idx)const{
  return l1data_.at(idx);
}
LumiSummary::HLT  
LumiSummary::hltinfo(unsigned int idx) const {
  return hltdata_.at(idx);
}
LumiSummary::HLT  
LumiSummary::hltinfo(const std::string& pathname) const {
  for(std::vector<HLT>::const_iterator it=hltdata_.begin();it!=hltdata_.end();++it){
    if(it->pathname==pathname) return *it;
  }
  return LumiSummary::HLT();
}
size_t
LumiSummary::nTriggerLine()const{
  return l1data_.size();
}
size_t
LumiSummary::nHLTPath()const{
  return hltdata_.size();
}
std::vector<std::string>
LumiSummary::HLTPaths()const{
  std::vector<std::string> result;
  for(std::vector<HLT>::const_iterator it=hltdata_.begin();it!=hltdata_.end();++it){
    result.push_back(it->pathname);
  }
  return result;
}
float
LumiSummary::avgInsRecLumi() const {
  return avginsdellumi_ * liveFrac(); 
}
float
LumiSummary::avgInsRecLumiErr() const {
  return avginsdellumierr_ * liveFrac(); 
}
bool
LumiSummary::isProductEqual(LumiSummary const& next) const {
  return (avginsdellumi_ == next.avginsdellumi_ &&
          avginsdellumierr_ == next.avginsdellumierr_ &&
          lumisecqual_ == next.lumisecqual_ &&
          deadcount_ == next.deadcount_ &&
          lsnumber_ == next.lsnumber_ &&
          startorbit_== next.startorbit_ &&
          numorbit_==next.numorbit_&&
          l1data_.size() == next.l1data_.size() &&
          hltdata_.size() == next.hltdata_.size() &&
          lumiversion_ == next.lumiversion_ );
}
std::string 
LumiSummary::lumiVersion()const{
  return lumiversion_;
}
void
LumiSummary::setLumiVersion(const std::string& lumiversion){
  lumiversion_=lumiversion;
}
void 
LumiSummary::setLumiData(float instlumi,float instlumierr,short lumiquality){
  avginsdellumi_=instlumi;
  avginsdellumierr_=instlumierr;
  lumisecqual_=lumiquality;
}
void 
LumiSummary::setDeadtime(unsigned long long deadcount){
  deadcount_=deadcount;
}
void 
LumiSummary::setlsnumber(unsigned int lsnumber){
  lsnumber_=lsnumber;
}
void 
LumiSummary::setOrbitData(unsigned int startorbit,unsigned int numorbit){
  startorbit_=startorbit;
  numorbit_=numorbit;
}
void 
LumiSummary::swapL1Data(std::vector<L1>& l1data){
  l1data_.swap(l1data);
}
void 
LumiSummary::swapHLTData(std::vector<HLT>& hltdata){
  hltdata_.swap(hltdata);
}
void 
LumiSummary::copyL1Data(const std::vector<L1>& l1data){
  l1data_.assign(l1data.begin(),l1data.end());
}
void 
LumiSummary::copyHLTData(const std::vector<HLT>& hltdata){
  hltdata_.assign(hltdata.begin(),hltdata.end());
}
std::ostream& operator<<(std::ostream& s, const LumiSummary& lumiSummary) {
  s << "\nDumping LumiSummary\n\n";
  if(!lumiSummary.isValid()){
    s << " === Invalid Lumi values === \n";
  }
  s << "  lumiVersion = " << lumiSummary.lumiVersion()  << "\n";
  s << "  avgInsDelLumi = " << lumiSummary.avgInsDelLumi() << "\n";
  s << "  avgInsDelLumiErr = " << lumiSummary.avgInsDelLumiErr() << "\n";
  s << "  lumiSecQual = " << lumiSummary.lumiSecQual() << "\n";
  s << "  deadCount = " << lumiSummary.deadcount() << "\n";
  s << "  deadFrac = " << (float)lumiSummary.deadFrac() << "\n";
  s << "  liveFrac = " << (float)lumiSummary.liveFrac() << "\n";
  s << "  lsNumber = " << lumiSummary.lsNumber() << "\n";
  s << "  startOrbit = " << lumiSummary.startOrbit() <<"\n";
  s << "  numOrbit = " << lumiSummary.numOrbit() <<"\n";
  s << "  avgInsRecLumi = " << lumiSummary.avgInsRecLumi() << "\n";
  s << "  avgInsRecLumiErr = "  << lumiSummary.avgInsRecLumiErr() << "\n\n";
  s << std::setw(15) << "l1name";
  s << std::setw(15) << "l1count";
  s << std::setw(15) << "l1prescale";
  s << "\n";
  size_t nTriggers=lumiSummary.nTriggerLine();
  size_t nHLTPath=lumiSummary.nHLTPath();
  for(unsigned int i = 0; i < nTriggers; ++i) {
    s << std::setw(15);
    s << lumiSummary.l1info(i).triggername;
    
    s << std::setw(15);
    s << lumiSummary.l1info(i).ratecount;
    
    s << std::setw(15);
    s << lumiSummary.l1info(i).prescale;
    s<<"\n";
  }
  s << std::setw(15) << "hltpath";
  s << std::setw(15) << "hltcount";
  s << std::setw(15) << "hltprescale";
  s << std::setw(15) << "hltinput";
  s << "\n";
  for(unsigned int i = 0; i < nHLTPath; ++i) {
    s << std::setw(15);
    s << lumiSummary.hltinfo(i).pathname;
    s << std::setw(15);
    s << lumiSummary.hltinfo(i).ratecount;
    s << std::setw(15);
    s << lumiSummary.hltinfo(i).prescale;
    s << std::setw(15);
    s << lumiSummary.hltinfo(i).inputcount;
    s << "\n";
  }
  return s << "\n";
}