L1TRPCTF Class Reference

#include <L1TRPCTF.h>

Inheritance diagram for L1TRPCTF:

Classes
struct	BxDelays

Public Member Functions
	L1TRPCTF (const edm::ParameterSet &ps)
virtual	~L1TRPCTF ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
void	beginJob (void)
void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context)
void	endJob (void)
void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
void	endRun (const edm::Run &r, const edm::EventSetup &c)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Member Functions
void	fillRateHistos (int orbit, bool flush=false)

Private Attributes
MonitorElement *	m_bxDiff
DQMStore *	m_dbe
std::set< unsigned long long int >	m_globBX
long long int	m_lastUsedBxInBxdiff
MonitorElement *	m_muonsEtaPhi [3]
MonitorElement *	m_qualVsEta [3]
MonitorElement *	m_rateAvg
int	m_rateBinSize
L1TRateHelper::L1TRateHelper	m_rateHelper
MonitorElement *	m_rateMax
MonitorElement *	m_rateMin
int	m_rateNoOfBins
int	m_rateUpdateTime
bool	monitorDaemon_
int	nev_
int	nevRPC_
std::string	output_dir_
std::string	outputFile_
MonitorElement *	rpctfbx
MonitorElement *	rpctfchargevalue [3]
MonitorElement *	rpctfcratesynchro [12]
MonitorElement *	rpctfetavalue [3]
MonitorElement *	rpctfntrack_b [3]
MonitorElement *	rpctfntrack_e [3]
MonitorElement *	rpctfphivalue [3]
MonitorElement *	rpctfptvalue [3]
MonitorElement *	rpctfquality [3]
edm::InputTag	rpctfSource_
bool	verbose_

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)

Detailed Description

Definition at line 45 of file L1TRPCTF.h.

Constructor & Destructor Documentation

L1TRPCTF::L1TRPCTF

(

const edm::ParameterSet &

ps

)

Definition at line 21 of file L1TRPCTF.cc.

References gather_cfg::cout, edm::ParameterSet::getUntrackedParameter(), m_dbe, NULL, cmsCodeRules.cppFunctionSkipper::operator, output_dir_, outputFile_, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), and verbose_.

  : rpctfSource_( ps.getParameter< InputTag >("rpctfSource") ),
//    digiSource_( ps.getParameter< InputTag >("rpctfRPCDigiSource") ),
//   m_rpcDigiFine(false),
//    m_useRpcDigi(true),
   m_rateUpdateTime( ps.getParameter< int >("rateUpdateTime") ),
   m_rateBinSize( ps.getParameter< int >("rateBinSize") ),
   m_rateNoOfBins( ps.getParameter< int >("rateNoOfBins") ),
   m_lastUsedBxInBxdiff(0),
   output_dir_ (ps.getUntrackedParameter<string>("output_dir") )
//    m_rpcDigiWithBX0(0),
//    m_rpcDigiWithBXnon0(0)

 {

  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);

  if(verbose_) cout << "L1TRPCTF: constructor...." << endl;


  m_dbe = NULL;
  if ( ps.getUntrackedParameter<bool>("DQMStore", false) ) 
  {
    m_dbe = Service<DQMStore>().operator->();
    m_dbe->setVerbose(0);
  }

  outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
  if ( outputFile_.size() != 0 ) {
    cout << "L1T Monitoring histograms will be saved to " << outputFile_.c_str() << endl;
  }

  bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false);
  if(disable){
    outputFile_="";
  }


  if ( m_dbe !=NULL ) {
    m_dbe->setCurrentFolder(output_dir_);
  }



}

L1TRPCTF::~L1TRPCTF ( )

virtual

Definition at line 68 of file L1TRPCTF.cc.

{
}

Member Function Documentation

void L1TRPCTF::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 291 of file L1TRPCTF.cc.

References L1TRateHelper::L1TRateHelper::addOrbit(), edm::EventBase::bunchCrossing(), L1TRPCTF::BxDelays::bx, gather_cfg::cout, diffTreeTool::diff, L1TRPCTF::BxDelays::eta_t, MonitorElement::Fill(), fillRateHistos(), edm::Event::getByLabel(), L1MuGMTReadoutCollection::getRecords(), i, edm::HandleBase::isValid(), j, edm::InputTag::label(), m_bxDiff, m_globBX, m_lastUsedBxInBxdiff, m_muonsEtaPhi, m_qualVsEta, m_rateHelper, min, nev_, edm::EventBase::orbitNumber(), L1TRPCTF::BxDelays::phi_p, edm::Handle< T >::product(), rpctfbx, rpctfchargevalue, rpctfcratesynchro, rpctfetavalue, rpctfntrack_b, rpctfntrack_e, rpctfphivalue, rpctfptvalue, rpctfquality, rpctfSource_, and verbose_.

{
  nev_++; 
  if(verbose_) cout << "L1TRPCTF: analyze...." << endl;

  edm::Handle<L1MuGMTReadoutCollection> pCollection;
  e.getByLabel(rpctfSource_,pCollection);
  
  if (!pCollection.isValid()) {
    edm::LogInfo("DataNotFound") << "can't find L1MuGMTReadoutCollection with label "
                   << rpctfSource_.label() ;
    return;
  }

  
  L1MuGMTReadoutCollection const* gmtrc = pCollection.product();
  vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
  vector<L1MuGMTReadoutRecord>::const_iterator RRItr;

  std::vector<int> nrpctftrack_b(3,0);
  std::vector<int> nrpctftrack_e(3,0);
  bool rpcCandsPresentInEvent = false;

  vector<L1TRPCTF::BxDelays> all_bxdelays;


  for( RRItr = gmt_records.begin() ;
       RRItr != gmt_records.end() ;
       RRItr++ ) 
  {
    
   if (verbose_) cout << "Readout Record " << RRItr->getBxInEvent() << endl;
   
   vector<vector<L1MuRegionalCand> > brlAndFwdCands;
   brlAndFwdCands.push_back(RRItr->getBrlRPCCands());
   brlAndFwdCands.push_back(RRItr->getFwdRPCCands());
  
   int beIndex = 0;
   vector<vector<L1MuRegionalCand> >::iterator RPCTFCands = brlAndFwdCands.begin();
   for(; RPCTFCands!= brlAndFwdCands.end(); ++RPCTFCands)
   {
      
      for( vector<L1MuRegionalCand>::const_iterator 
          ECItr = RPCTFCands->begin() ;
          ECItr != RPCTFCands->end() ;
          ++ECItr ) 
      {
  
        int bxindex = 1 ; // bx == 0
        if (ECItr->bx() > 0) bxindex = 2;
        if (ECItr->bx() < 0) bxindex = 0;
        
        if (!ECItr->empty()) {
          
          rpcCandsPresentInEvent = true;
          
          if (beIndex == 0) ++nrpctftrack_b[bxindex];
          if (beIndex == 1) ++nrpctftrack_e[bxindex];
    
          if (verbose_) cout << "RPCTFCand bx " << ECItr->bx() << endl;
          
          int tower = ECItr->eta_packed();
          if (tower > 16) {
            tower = - ( (~tower & 63) + 1);
          }

          rpctfbx->Fill(ECItr->bx());
    
          rpctfetavalue[bxindex]->Fill(tower);
          if (verbose_) cout << "\tRPCTFCand eta value " << ECItr->etaValue() << endl;
  
          rpctfphivalue[bxindex]->Fill(ECItr->phi_packed());
          if (verbose_) cout << "\tRPCTFCand phi value " << ECItr->phiValue() << endl;
    
          rpctfptvalue[bxindex]->Fill(ECItr->ptValue());
          if (verbose_) cout << "\tRPCTFCand pt value " << ECItr->ptValue()<< endl;
    
          rpctfchargevalue[bxindex]->Fill(ECItr->chargeValue());
          if (verbose_) cout << "\tRPCTFCand charge value " << ECItr->chargeValue() << endl;
    
          rpctfquality[bxindex]->Fill(ECItr->quality());
          if (verbose_) cout << "\tRPCTFCand quality " << ECItr->quality() << endl;
          

          m_qualVsEta[bxindex]->Fill(tower, ECItr->quality());
          m_muonsEtaPhi[bxindex]->Fill(tower, ECItr->phi_packed());

          BxDelays bx_del;
          bx_del.bx = ECItr->bx();
          bx_del.eta_t = tower;
          bx_del.phi_p = ECItr->phi_packed();
          all_bxdelays.push_back(bx_del);
          
        } // if !empty
      } // end candidates iteration
      ++beIndex;
   } // end brl/endcap iteration
  } // end GMT records iteration

  for (int bxI = 0; bxI < 3; ++bxI){
    rpctfntrack_b[bxI]->Fill(nrpctftrack_b[bxI]);
    rpctfntrack_e[bxI]->Fill(nrpctftrack_e[bxI]);
  }
  
  
  if (rpcCandsPresentInEvent) {
    m_rateHelper.addOrbit(e.orbitNumber());
        
    unsigned int globBx = e.orbitNumber()*3564+e.bunchCrossing();
    if (m_globBX.find(globBx)==m_globBX.end()) m_globBX.insert(globBx);
    if (m_globBX.size()>1020){
      long long int diff = *m_globBX.begin()-m_lastUsedBxInBxdiff; // first entry will go to overflow bin, ignore
      m_bxDiff->Fill(diff);
      m_lastUsedBxInBxdiff = *m_globBX.begin();
      m_globBX.erase(m_globBX.begin());
    
    }
    
    
    
  }


   for(unsigned int i = 0; i < all_bxdelays.size(); i++) {

     int sector= ((all_bxdelays[i].phi_p+ 142)%144)/12;
     if (sector>11 || sector < 0) continue;
     int eta_tower = all_bxdelays[i].eta_t;
     for(unsigned int j = 0; j < all_bxdelays.size(); j++) {
       if(i == j) continue;
       int sector2= ((all_bxdelays[j].phi_p + 142)%144)/12;
 
       int distance_cut = 1;
       int distance = ((sector+12)-sector2)%12;
       distance = min(distance, 11-distance);
       if(distance<distance_cut) continue;
 
       int bxDiff = all_bxdelays[i].bx-all_bxdelays[j].bx;
       rpctfcratesynchro[sector]->Fill(bxDiff,eta_tower);
     }

  }
  fillRateHistos(e.orbitNumber());
  
  
    
}

void L1TRPCTF::beginJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 72 of file L1TRPCTF.cc.

References b, DQMStore::book1D(), DQMStore::book2D(), i, prof2calltree::l, m_bxDiff, m_dbe, m_muonsEtaPhi, m_qualVsEta, m_rateAvg, m_rateMax, m_rateMin, m_rateNoOfBins, nev_, nevRPC_, connectstrParser::o, cmsCodeRules.cppFunctionSkipper::operator, output_dir_, DQMStore::rmdir(), rpctfbx, rpctfchargevalue, rpctfcratesynchro, rpctfetavalue, rpctfntrack_b, rpctfntrack_e, rpctfphivalue, rpctfptvalue, rpctfquality, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().

{

  nev_ = 0;
  nevRPC_ = 0;

  // get hold of back-end interface
  m_dbe = Service<DQMStore>().operator->();

  if ( m_dbe ) {
    m_dbe->setCurrentFolder(output_dir_);
    m_dbe->rmdir(output_dir_);
  }


  if ( m_dbe ) 
  {
    

    ostringstream oDir; oDir<< output_dir_ << "/CrateSynchroHistograms/";
    m_dbe->setCurrentFolder(oDir.str());
    for( unsigned int i = 0; i < 12; i++) {
      
       ostringstream o; o<<"RPCTF_crate_"<<i<<"_synchro";
       rpctfcratesynchro[i] = m_dbe->book2D(o.str(), o.str(), 5, -2.5, 2.5, 33, -16.5, 16.5);
       for (int bx = -2; bx < 3; ++bx){
          ostringstream b; b<<"BX="<<bx;
          rpctfcratesynchro[i]->setBinLabel(bx+3, b.str(),1);
       }
       rpctfcratesynchro[i]->setAxisTitle("Tower",2);
     
    }
    m_dbe->setCurrentFolder(output_dir_);
    
    rpctfetavalue[1] = m_dbe->book1D("RPCTF_eta_value_bx0", 
       "RPCTF eta value bx=0", 33, -16.5, 16.5 ) ;
    rpctfetavalue[2] = m_dbe->book1D("RPCTF_eta_value_bx+", 
       "RPCTF eta value bx>0", 33, -16.5, 16.5 ) ;
    rpctfetavalue[0] = m_dbe->book1D("RPCTF_eta_value_bx-", 
       "RPCTF eta value bx<0", 33, -16.5, 16.5 ) ;
    
    rpctfphivalue[1] = m_dbe->book1D("RPCTF_phi_value_bx0", 
       "RPCTF phi value bx=0", 144, -0.5, 143.5) ;
    rpctfphivalue[2] = m_dbe->book1D("RPCTF_phi_value_bx+", 
       "RPCTF phi value bx>0", 144, -0.5, 143.5 ) ;
    rpctfphivalue[0] = m_dbe->book1D("RPCTF_phi_value_bx-", 
       "RPCTF phi value bx<0", 144, -0.5, 143.5 ) ;
      
       
       
    rpctfptvalue[1] = m_dbe->book1D("RPCTF_pt_value_bx0", 
                                    "RPCTF pt value bx=0", 160, -0.5, 159.5 );
    rpctfptvalue[2] = m_dbe->book1D("RPCTF_pt_value_bx+", 
                                    "RPCTF pt value bx>0", 160, -0.5, 159.5 );
    rpctfptvalue[0] = m_dbe->book1D("RPCTF_pt_value_bx-", 
                                    "RPCTF pt value bx<0", 160, -0.5, 159.5 );
    
    
    rpctfchargevalue[1] = m_dbe->book1D("RPCTF_charge_value_bx0", 
                                        "RPCTF charge value bx=0", 3, -1.5, 1.5 ) ;
    rpctfchargevalue[2] = m_dbe->book1D("RPCTF_charge_value_bx+", 
                                        "RPCTF charge value bx>0", 3, -1.5, 1.5 ) ;
    rpctfchargevalue[0] = m_dbe->book1D("RPCTF_charge_value_bx-", 
                                        "RPCTF charge value bx<01", 3, -1.5, 1.5 ) ;

    rpctfquality[1] = m_dbe->book1D("RPCTF_quality", 
                                    "RPCTF quality bx=0", 6, -0.5, 5.5 ) ;
    rpctfquality[2] = m_dbe->book1D("RPCTF_quality_bx+", 
                                    "RPCTF quality bx>0", 6, -0.5, 5.5 ) ;
    rpctfquality[0] = m_dbe->book1D("RPCTF_quality_bx-", 
                                    "RPCTF quality bx<0", 6, -0.5, 5.5 ) ;

    rpctfntrack_b[1] = m_dbe->book1D("RPCTF_ntrack_brl_bx0", 
                                     "RPCTF number of tracks - barrel, bx=0", 5, -0.5, 4.5 ) ;
    rpctfntrack_b[2] = m_dbe->book1D("RPCTF_ntrack_brl_bx+", 
                                     "RPCTF number of tracks - barrel, bx>0", 5, -0.5, 4.5 ) ;
    rpctfntrack_b[0] = m_dbe->book1D("RPCTF_ntrack_brl_bx-", 
                                     "RPCTF number of tracks - barrel, bx<0", 5, -0.5, 4.5 ) ;
    
    
           
    rpctfntrack_e[1] = m_dbe->book1D("RPCTF_ntrack_fwd_bx0", 
                                     "RPCTF number of tracks - endcap, bx=0", 5, -0.5, 4.5 ) ;
    rpctfntrack_e[2] = m_dbe->book1D("RPCTF_ntrack_fwd_bx+", 
                                     "RPCTF number of tracks - endcap, bx>0", 5, -0.5, 4.5 ) ;
    rpctfntrack_e[0] = m_dbe->book1D("RPCTF_ntrack_fwd_bx-", 
                                     "RPCTF number of tracks - endcap, bx<0", 5, -0.5, 4.5 ) ;

    
                                        
       

    m_qualVsEta[1] = m_dbe->book2D("RPCTF_quality_vs_eta_bx0", 
                              "RPCTF quality vs eta, bx=0", 
                               33, -16.5, 16.5,
                               6, -0.5, 5.5); // Currently only 0...3 quals are possible
    m_qualVsEta[2] = m_dbe->book2D("RPCTF_quality_vs_eta_bx+", 
                                   "RPCTF quality vs eta, bx>0", 
                                   33, -16.5, 16.5,
                                   6, -0.5, 5.5); // Currently only 0...3 quals are possible
    m_qualVsEta[0] = m_dbe->book2D("RPCTF_quality_vs_eta_bx-", 
                                   "RPCTF quality vs eta, bx<0", 
                                   33, -16.5, 16.5,
                                   6, -0.5, 5.5); // Currently only 0...3 quals are possible
    
    
        
    m_muonsEtaPhi[1] = m_dbe->book2D("RPCTF_muons_eta_phi_bx0", 
                                  "RPCTF occupancy(eta,phi), bx=0",  
                                  33, -16.5, 16.5,
                                  144,  -0.5, 143.5);
    m_muonsEtaPhi[2] = m_dbe->book2D("RPCTF_muons_eta_phi_bx+", 
                                     "RPCTF occupancy(eta,phi), bx>0",  
                                     33, -16.5, 16.5,
                                     144,  -0.5, 143.5);
    m_muonsEtaPhi[0] = m_dbe->book2D("RPCTF_muons_eta_phi_bx-", 
                                     "RPCTF occupancy(eta,phi), bx<0",  
                                     33, -16.5, 16.5,
                                     144,  -0.5, 143.5);

    rpctfbx = m_dbe->book1D("RPCTF_bx", 
                            "RPCTF bx distribiution", 7, -3.5, 3.5 );
    
    //axis labels
    for (int l = 0; l<3; ++l){
      m_muonsEtaPhi[l]->setAxisTitle("tower",1);
      m_qualVsEta[l]->setAxisTitle("tower");
      rpctfetavalue[l]->setAxisTitle("tower");
      
      m_muonsEtaPhi[l]->setAxisTitle("phi",2);
      rpctfphivalue[l]->setAxisTitle("phi");
    }
    
    // set phi bin labels
    for (int i = 0; i < 12 ; ++i ){
       //float lPhi  = (30./360)*i*2*3.14;
      int lPhi  = 30*i;
      int lBin = int((30./360)*i*144)+1;
      std::stringstream ss;
      ss << "phi=" <<lPhi;
      for (int l = 0; l<3; ++l){
        rpctfphivalue[l]->setBinLabel(lBin,ss.str());
        m_muonsEtaPhi[l]->setBinLabel(lBin,ss.str(), 2);
      }
    }

    /*
    // set TC numbers on phi axis
    for (int tc = 0; tc < 12 ; ++tc ){
      int lBin  = (tc*12+3+1)%144;
      std::stringstream ss;
      ss << "TC" <<tc;
      for (int l = 0; l<3; ++l){
        rpctfphivalue[l]->setBinLabel(lBin,ss.str());
        m_muonsEtaPhi[l]->setBinLabel(lBin,ss.str(), 2);
      }
  }*/

        
    // set eta bin labels
    for (int i = -16; i < 17 ; ++i ){
      std::stringstream ss;
      ss << i;
      for (int l = 0; l<3; ++l){
        rpctfetavalue[l]->setBinLabel(i+17, ss.str());
        m_muonsEtaPhi[l]->setBinLabel(i+17, ss.str(), 1);
        m_qualVsEta[l]->setBinLabel(i+17, ss.str());
      }
    }

    
    
    m_rateMin = m_dbe->book1D("RPCTF_rate_min",
                              "RPCTrigger - minimal rate", m_rateNoOfBins, 0, m_rateNoOfBins); 
    
    m_rateMax = m_dbe->book1D("RPCTF_rate_max",
                              "RPCTrigger - peak rate", m_rateNoOfBins, 0, m_rateNoOfBins);

    m_rateAvg = m_dbe->book1D("RPCTF_rate_avg",
                              "RPCTrigger - average rate", m_rateNoOfBins, 0, m_rateNoOfBins); 
                              
    m_bxDiff = m_dbe->book1D("RPCTF_bx_diff",
                  "RPCTrigger - bx difference", 12000, -.5, 11999.5); 
   


  }   // if (m_dbe)
}

void L1TRPCTF::beginLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg, const edm::EventSetup &  context  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 506 of file L1TRPCTF.cc.

{
//    m_rpcDigiWithBX0=0;
//    m_rpcDigiWithBXnon0=0;
//    m_bxs.clear();
//    m_useRpcDigi = true;

                          
}

void L1TRPCTF::endJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 279 of file L1TRPCTF.cc.

References gather_cfg::cout, m_dbe, nev_, outputFile_, DQMStore::save(), and verbose_.

{
  
  if(verbose_) cout << "L1TRPCTF: end job...." << endl;
  LogInfo("EndJob") << "analyzed " << nev_ << " events"; 

  if ( outputFile_.size() != 0  && m_dbe ) m_dbe->save(outputFile_);
    
  return;

}

void L1TRPCTF::endLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg, const edm::EventSetup &  c  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 518 of file L1TRPCTF.cc.

{

}

void L1TRPCTF::endRun ( const edm::Run &  r, const edm::EventSetup &  c  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 261 of file L1TRPCTF.cc.

References diffTreeTool::diff, MonitorElement::Fill(), fillRateHistos(), m_bxDiff, m_globBX, and m_lastUsedBxInBxdiff.

                                                              {
  
      fillRateHistos(0,true);


     // fixme, norm iteration would be better
     while (m_globBX.begin() !=  m_globBX.end() ) {
        long long int  diff = *m_globBX.begin()-m_lastUsedBxInBxdiff; // first entry will go to overflow bin, ignore
        m_bxDiff->Fill(diff);
        m_lastUsedBxInBxdiff = *m_globBX.begin();
        m_globBX.erase(m_globBX.begin());

     }

     
}

void L1TRPCTF::fillRateHistos ( int  orbit, bool  flush = false  )

private

Fills rate histos.

Definition at line 442 of file L1TRPCTF.cc.

References newFWLiteAna::bin, L1TRateHelper::L1TRateHelper::getEarliestTime(), L1TRateHelper::L1TRateHelper::getTimeForOrbit(), m_rateAvg, m_rateBinSize, m_rateHelper, m_rateMax, m_rateMin, m_rateUpdateTime, max(), min, L1TEmulatorMonitor_cff::p, L1TRateHelper::L1TRateHelper::removeAndGetRateForEarliestTime(), and MonitorElement::setBinContent().

Referenced by analyze(), and endRun().

{
  
  static bool flushed = false;

  if (flushed) {
    LogWarning("L1TRPCTF") << "Rate histos allready flushed \n";
  }
    
  if (flush) flushed = true;
  

  int nbinsUsed = 0;
  do  {
    nbinsUsed = 0;
    int et = m_rateHelper.getEarliestTime();
    if (et==-1) break;
    
    if ( (( m_rateHelper.getTimeForOrbit(orbit) - et > m_rateUpdateTime+m_rateBinSize)  // 1 minute bins
            && m_rateUpdateTime!=-1) || flush  )
    {
      
      
      int startTimeInMinutes=et/m_rateBinSize; 
      int bin = 0;
      std::pair<int, int> p; 
      int max = 0, min = 0;
      float avg=0;
      int curTimeInMinutes=startTimeInMinutes;
      while (curTimeInMinutes==startTimeInMinutes){
        p = m_rateHelper.removeAndGetRateForEarliestTime(); 
        if (p.first < 0) break; // no more items to analize, go fill histos
        
        if (nbinsUsed==0) {
          bin = p.first/m_rateBinSize+1;
          max = p.second;
          min = p.second;
        } else {
          if (max < p.second) max = p.second;
          if (min > p.second) min = p.second;
        }
        
        ++nbinsUsed;
        avg+=p.second;
        curTimeInMinutes=m_rateHelper.getEarliestTime()/m_rateBinSize;
      }
      
      
      avg/=m_rateBinSize;
      
      if (nbinsUsed > 0){
        m_rateAvg->setBinContent(bin,avg); // smallest possible value in f.first is 1
        m_rateMin->setBinContent(bin,min);
        m_rateMax->setBinContent(bin,max);
      }
    
    }
  } while (flush);



}

Member Data Documentation

MonitorElement* L1TRPCTF::m_bxDiff

private

Definition at line 95 of file L1TRPCTF.h.

Referenced by analyze(), beginJob(), and endRun().

DQMStore* L1TRPCTF::m_dbe

private

Definition at line 78 of file L1TRPCTF.h.

Referenced by beginJob(), endJob(), and L1TRPCTF().

std::set<unsigned long long int> L1TRPCTF::m_globBX

private

Definition at line 97 of file L1TRPCTF.h.

Referenced by analyze(), and endRun().

long long int L1TRPCTF::m_lastUsedBxInBxdiff

private

Definition at line 116 of file L1TRPCTF.h.

Referenced by analyze(), and endRun().

MonitorElement* L1TRPCTF::m_muonsEtaPhi[3]

private

Definition at line 89 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::m_qualVsEta[3]

private

Definition at line 88 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::m_rateAvg

private

Definition at line 92 of file L1TRPCTF.h.

Referenced by beginJob(), and fillRateHistos().

int L1TRPCTF::m_rateBinSize

private

Definition at line 114 of file L1TRPCTF.h.

Referenced by fillRateHistos().

L1TRateHelper::L1TRateHelper L1TRPCTF::m_rateHelper

private

Definition at line 100 of file L1TRPCTF.h.

Referenced by analyze(), and fillRateHistos().

MonitorElement* L1TRPCTF::m_rateMax

private

Definition at line 94 of file L1TRPCTF.h.

Referenced by beginJob(), and fillRateHistos().

MonitorElement* L1TRPCTF::m_rateMin

private

Definition at line 93 of file L1TRPCTF.h.

Referenced by beginJob(), and fillRateHistos().

int L1TRPCTF::m_rateNoOfBins

private

Definition at line 115 of file L1TRPCTF.h.

Referenced by beginJob().

int L1TRPCTF::m_rateUpdateTime

private

Definition at line 113 of file L1TRPCTF.h.

Referenced by fillRateHistos().

bool L1TRPCTF::monitorDaemon_

private

Definition at line 109 of file L1TRPCTF.h.

int L1TRPCTF::nev_

private

Definition at line 105 of file L1TRPCTF.h.

Referenced by analyze(), beginJob(), and endJob().

int L1TRPCTF::nevRPC_

private

Definition at line 106 of file L1TRPCTF.h.

Referenced by beginJob().

std::string L1TRPCTF::output_dir_

private

Definition at line 117 of file L1TRPCTF.h.

Referenced by beginJob(), and L1TRPCTF().

std::string L1TRPCTF::outputFile_

private

Definition at line 107 of file L1TRPCTF.h.

Referenced by endJob(), and L1TRPCTF().

MonitorElement* L1TRPCTF::rpctfbx

private

Definition at line 87 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::rpctfchargevalue[3]

private

Definition at line 83 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::rpctfcratesynchro[12]

private

Definition at line 96 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::rpctfetavalue[3]

private

Definition at line 80 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::rpctfntrack_b[3]

private

Definition at line 85 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::rpctfntrack_e[3]

private

Definition at line 86 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::rpctfphivalue[3]

private

Definition at line 81 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::rpctfptvalue[3]

private

Definition at line 82 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TRPCTF::rpctfquality[3]

private

Definition at line 84 of file L1TRPCTF.h.

Referenced by analyze(), and beginJob().

edm::InputTag L1TRPCTF::rpctfSource_

private

Definition at line 103 of file L1TRPCTF.h.

Referenced by analyze().

bool L1TRPCTF::verbose_

private

Definition at line 108 of file L1TRPCTF.h.

Referenced by analyze(), endJob(), and L1TRPCTF().