#include <RPCEventSummary.h>

Inheritance diagram for RPCEventSummary:

Public Member Functions
	RPCEventSummary (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~RPCEventSummary ()
	Destructor. More...

Public Member Functions inherited from DQMEDHarvester
virtual void	analyze (edm::Event const &, edm::EventSetup const &) final

virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

virtual void	beginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDHarvester (void)

virtual void	endJob () final

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

virtual void	endRun (edm::Run const &, edm::EventSetup const &)

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	beginJob () override

void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &) override

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 Types
enum	RPCQualityFlags { DEAD = 6, PARTIALLY_DEAD =5 }

Private Member Functions
void	clientOperation (DQMStore::IGetter &igetter)

Private Attributes
bool	doEndcapCertification_

bool	enableReportSummary_

std::string	eventInfoPath_

std::pair< int, int >	FEDRange_

std::string	globalFolder_

bool	init_

bool	isIn_

int	lumiCounter_

int	minimumEvents_

int	numberDisk_

int	NumberOfFeds_

bool	offlineDQM_

std::string	prefixDir_

std::string	prefixFolder_

int	prescaleFactor_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 13 of file RPCEventSummary.h.

Member Enumeration Documentation

enum RPCEventSummary::RPCQualityFlags

private

Enumerator
DEAD
PARTIALLY_DEAD

Definition at line 48 of file RPCEventSummary.h.

48 {DEAD = 6, PARTIALLY_DEAD=5};

RPCEventSummary::PARTIALLY_DEAD

Definition: RPCEventSummary.h:48

RPCEventSummary::DEAD

Definition: RPCEventSummary.h:48

Constructor & Destructor Documentation

RPCEventSummary::RPCEventSummary ( const edm::ParameterSet & ps )

Constructor.

Definition at line 16 of file RPCEventSummary.cc.

References doEndcapCertification_, enableReportSummary_, eventInfoPath_, FEDRange_, edm::ParameterSet::getUntrackedParameter(), globalFolder_, minimumEvents_, numberDisk_, NumberOfFeds_, offlineDQM_, prefixFolder_, prescaleFactor_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                           {
   edm::LogVerbatim ("rpceventsummary") << "[RPCEventSummary]: Constructor";
 
   enableReportSummary_ = ps.getUntrackedParameter<bool>("EnableSummaryReport",true);
   prescaleFactor_ =  ps.getUntrackedParameter<int>("PrescaleFactor", 1);
   eventInfoPath_ = ps.getUntrackedParameter<std::string>("EventInfoPath", "RPC/EventInfo");
 
 
   std::string subsystemFolder = ps.getUntrackedParameter<std::string>("RPCFolder", "RPC");
   std::string recHitTypeFolder = ps.getUntrackedParameter<std::string>("RecHitTypeFolder", "AllHits");
   std::string summaryFolder = ps.getUntrackedParameter<std::string>("SummaryFolder", "SummaryHistograms");
 
   globalFolder_  =  subsystemFolder +"/"+  recHitTypeFolder +"/"+ summaryFolder ;
   prefixFolder_  =  subsystemFolder +"/"+  recHitTypeFolder ;
 
   minimumEvents_= ps.getUntrackedParameter<int>("MinimumRPCEvents", 10000);
   numberDisk_ = ps.getUntrackedParameter<int>("NumberOfEndcapDisks", 4);
   doEndcapCertification_ = ps.getUntrackedParameter<bool>("EnableEndcapSummary", false);
 
   FEDRange_.first  = ps.getUntrackedParameter<unsigned int>("MinimumRPCFEDId", 790);
   FEDRange_.second = ps.getUntrackedParameter<unsigned int>("MaximumRPCFEDId", 792);
   
   NumberOfFeds_ =FEDRange_.second -  FEDRange_.first +1;
 
   offlineDQM_ = ps.getUntrackedParameter<bool> ("OfflineDQM",true); 
 
 
 }

RPCEventSummary::~RPCEventSummary ( )

virtual

Destructor.

Definition at line 45 of file RPCEventSummary.cc.

                                  {
   edm::LogVerbatim ("rpceventsummary") << "[RPCEventSummary]: Destructor ";
 }

Member Function Documentation

void RPCEventSummary::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 49 of file RPCEventSummary.cc.

References init_.

                               {
  edm::LogVerbatim ("rpceventsummary") << "[RPCEventSummary]: Begin job ";
  init_ = false;  
 }

void RPCEventSummary::clientOperation ( DQMStore::IGetter & igetter )

private

Definition at line 185 of file RPCEventSummary.cc.

References edmIntegrityCheck::d, DEAD, doEndcapCertification_, eventInfoPath_, MonitorElement::Fill(), DQMStore::IGetter::get(), MonitorElement::getBinContent(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), globalFolder_, mps_fire::i, minimumEvents_, numberDisk_, PARTIALLY_DEAD, prefixFolder_, alignCSCRings::r, alignCSCRings::s, MonitorElement::setBinContent(), w, and y.

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

                                                                {
 
   float  rpcevents = minimumEvents_;
   MonitorElement *  RPCEvents ;
   RPCEvents = igetter.get( prefixFolder_  +"/RPCEvents");  
 
   if(RPCEvents) {
     rpcevents = RPCEvents ->getBinContent(1);
   }
   
 
   if(rpcevents < minimumEvents_) return;
   std::stringstream meName;
   MonitorElement * myMe;
    
   meName.str("");
   meName<<eventInfoPath_ + "/reportSummaryMap";
   MonitorElement * reportMe = igetter.get(meName.str());
   
   MonitorElement * globalMe;
   
   //BARREL
   float barrelFactor = 0;
   
   for(int w = -2 ; w<3; w++){
     
     meName.str("");
     meName<<globalFolder_<<"/RPCChamberQuality_Roll_vs_Sector_Wheel"<<w;
     myMe = igetter.get(meName.str());
     
        if(myMe){      
      float wheelFactor = 0;
      
      for(int s = 1; s<=myMe->getNbinsX() ; s++){
        float sectorFactor = 0;
        int rollInSector = 0;
      
        
      for(int r = 1;r<=myMe->getNbinsY(); r++){
        if((s!=4 && r > 17 ) || ((s ==9 ||s ==10)  && r >15 ) )  continue;
        rollInSector++;
        
        
        if(myMe->getBinContent(s,r) == PARTIALLY_DEAD) sectorFactor+=0.8;
        else if(myMe->getBinContent(s,r) == DEAD )sectorFactor+=0;
        else sectorFactor+=1;    
        
      }
      if(rollInSector!=0)
        sectorFactor = sectorFactor/rollInSector;
      
      if(reportMe)   reportMe->setBinContent(w+8, s, sectorFactor);
      wheelFactor += sectorFactor;
      
      }//end loop on sectors
      
      wheelFactor = wheelFactor/myMe->getNbinsX();
      
      meName.str("");
      meName<<eventInfoPath_ + "/reportSummaryContents/RPC_Wheel"<<w; 
      globalMe=igetter.get(meName.str());
      if(globalMe) globalMe->Fill(wheelFactor);
      
      barrelFactor += wheelFactor;
        }//
      }//end loop on wheel
   
      barrelFactor = barrelFactor/5;
   
 
      float endcapFactor = 0;
      
      if(doEndcapCertification_){
        
        //Endcap
        for(int d = -numberDisk_ ; d<= numberDisk_; d++){
      if (d==0) continue;
      
      meName.str("");
      meName<<globalFolder_<<"/RPCChamberQuality_Ring_vs_Segment_Disk"<<d;
      myMe = igetter.get(meName.str());
      
      if(myMe){      
        float diskFactor = 0;
        
        float sectorFactor[6]= {0,0,0,0,0,0};
        
        for (int i = 0 ;i <6;i++){
          int firstSeg = (i *6 )+1;
          int lastSeg = firstSeg +6;
          int rollInSector = 0; 
          for(int seg = firstSeg; seg< lastSeg ; seg++){
            
            for(int y = 1;y<=myMe->getNbinsY(); y++){
          rollInSector++;
          if(myMe->getBinContent(seg,y) == PARTIALLY_DEAD) sectorFactor[i]+=0.8;
          else if(myMe->getBinContent(seg,y) == DEAD )sectorFactor[i]+=0;
          else sectorFactor[i]+=1;   
          
            }
          }
          sectorFactor[i] = sectorFactor[i]/rollInSector;
        }//end loop on Sectors
        
        
        for (int sec = 0 ; sec<6; sec++){
          diskFactor += sectorFactor[sec];   
          if(reportMe)   {
            if (d<0) reportMe->setBinContent(d+5, sec+1 , sectorFactor[sec]);
            else  reportMe->setBinContent(d+11, sec+1 , sectorFactor[sec]);
          }   
        }
        
        diskFactor = diskFactor/6;
        
        meName.str("");
        meName<<eventInfoPath_ + "/reportSummaryContents/RPC_Disk"<<d; 
        globalMe=igetter.get(meName.str());
        if(globalMe) globalMe->Fill(diskFactor);
        
        endcapFactor += diskFactor;
      }//end loop on disks
      
        }
        
        endcapFactor=endcapFactor/ (numberDisk_ * 2);
        
      } 
      
      //Fill repor summary
      float rpcFactor = barrelFactor;
      if(doEndcapCertification_){ rpcFactor =  ( barrelFactor + endcapFactor)/2; }
      
      globalMe = igetter.get(eventInfoPath_ +"/reportSummary"); 
      if(globalMe) globalMe->Fill(rpcFactor);
      
      
 }

void RPCEventSummary::dqmEndJob	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 180 of file RPCEventSummary.cc.

References clientOperation(), and isIn_.

                                                                                    { 
   
   if(isIn_) { this->clientOperation(igetter);}
 }

void RPCEventSummary::dqmEndLuminosityBlock	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		edm::LuminosityBlock const &	lb,
		edm::EventSetup const &	setup
	)

overrideprotectedvirtual

Reimplemented from DQMEDHarvester.

Definition at line 54 of file RPCEventSummary.cc.

References CentralityFilter_cfi::BinLabel, DQMStore::IBooker::book2D(), DQMStore::IBooker::bookFloat(), clientOperation(), eventInfoPath_, f, FEDRange_, MonitorElement::Fill(), edm::EventSetup::find(), edm::eventsetup::heterocontainer::HCTypeTag::findType(), edm::EventSetup::get(), mps_fire::i, init_, isIn_, lumiCounter_, RunInfo::m_fed_in, numberDisk_, offlineDQM_, prescaleFactor_, MonitorElement::setBinContent(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                                                                           {
  
   edm::LogVerbatim ("rpceventsummary") << "[RPCEventSummary]: Begin run";
   
   
   if(!init_){
     lumiCounter_ = prescaleFactor_ ;
     
     edm::eventsetup::EventSetupRecordKey recordKey(edm::eventsetup::EventSetupRecordKey::TypeTag::findType("RunInfoRcd"));
     
     int defaultValue = 1;
     isIn_ = true;
 
     if(0 != setup.find( recordKey ) ) {
       defaultValue = -1;
       //get fed summary information
       edm::ESHandle<RunInfo> sumFED;
       setup.get<RunInfoRcd>().get(sumFED);    
       std::vector<int> FedsInIds= sumFED->m_fed_in;   
       unsigned int f = 0;
       isIn_ = false;
       while(!isIn_ && f < FedsInIds.size()) {
     int fedID=FedsInIds[f];
     //make sure fed id is in allowed range  
     if(fedID>=FEDRange_.first && fedID<=FEDRange_.second) {
       defaultValue = 1;
       isIn_ = true;
     } 
       f++;
       }   
     }   
     
     MonitorElement* me;
     ibooker.setCurrentFolder(eventInfoPath_);
     
     //a global summary float [0,1] providing a global summary of the status 
     //and showing the goodness of the data taken by the the sub-system 
     std::string histoName="reportSummary";
     me=0;
     me = ibooker.bookFloat(histoName);
     me->Fill(defaultValue);
     
     //TH2F ME providing a mapof values[0-1] to show if problems are localized or distributed
     me=0;    
     me = ibooker.book2D("reportSummaryMap", "RPC Report Summary Map", 15, -7.5, 7.5, 12, 0.5 ,12.5);
     
     //customize the 2d histo
     std::stringstream BinLabel;
     for (int i= 1 ; i<=15; i++){
       BinLabel.str("");
       if(i<13){
     BinLabel<<"Sec"<<i;
     me->setBinLabel(i,BinLabel.str(),2);
       } 
       
       BinLabel.str("");
       if(i<5)
     BinLabel<<"Disk"<<i-5;
       else if(i>11)
     BinLabel<<"Disk"<<i-11;
       else if(i==11 || i==5)
     BinLabel.str("");
       else
     BinLabel<<"Wheel"<<i-8;
       
       me->setBinLabel(i,BinLabel.str(),1);
     }
     
     //fill the histo with "1" --- just for the moment
     for(int i=1; i<=15; i++){
       for (int j=1; j<=12; j++ ){
     if(i==5 || i==11 || (j>6 && (i<6 || i>10)))    
       me->setBinContent(i,j,-1);//bins that not correspond to subdetector parts
     else
       me->setBinContent(i,j,defaultValue);
       }
     }
     
     if(numberDisk_ < 4)
       for (int j=1; j<=12; j++ ){
     me->setBinContent(1,j,-1);//bins that not correspond to subdetector parts
     me->setBinContent(15,j,-1);
       }
     
     //the reportSummaryContents folder containins a collection of ME floats [0-1] (order of 5-10)
     // which describe the behavior of the respective subsystem sub-components.
     ibooker.setCurrentFolder(eventInfoPath_+ "/reportSummaryContents");
     
     std::stringstream segName;
     std::vector<std::string> segmentNames;
     for(int i=-2; i<=2; i++){
       segName.str("");
       segName<<"RPC_Wheel"<<i;
       segmentNames.push_back(segName.str());
     }
     
     for(int i=-numberDisk_; i<=numberDisk_; i++){
       if(i==0) continue;
       segName.str("");
       segName<<"RPC_Disk"<<i;
       segmentNames.push_back(segName.str());
     }
     
     
     for(unsigned int i=0; i<segmentNames.size(); i++){
       me =0;
       me = ibooker.bookFloat(segmentNames[i]);
       me->Fill(defaultValue);
     }
     
     lumiCounter_ = prescaleFactor_;
     init_ = true;
   }
 
   
   
   if(isIn_ && !offlineDQM_  && lumiCounter_%prescaleFactor_ == 0 ){
     this->clientOperation(igetter);
   }
 
   lumiCounter_++;
 
 }