RPCChamberQuality Class Reference

#include <RPCChamberQuality.h>

Inheritance diagram for RPCChamberQuality:

Public Member Functions
	RPCChamberQuality (const edm::ParameterSet &ps)
virtual	~RPCChamberQuality ()
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
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Protected Member Functions
void	beginJob ()
void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)
Protected Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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	chamberQualityState {   GoodState = 1, OffState =2, NoisyStripState = 3, NoisyRollState = 4,   PartiallyDeadState =5, DeadState =6, BadShapeState =7 }

Private Member Functions
void	fillMonitorElements (DQMStore::IGetter &)
void	myBooker (DQMStore::IBooker &)
void	performeClientOperation (std::string, int, MonitorElement *, DQMStore::IGetter &)

Private Attributes
bool	enableDQMClients_
bool	init_
int	lumiCounter_
int	minEvents
int	numberOfDisks_
int	numLumBlock_
bool	offlineDQM_
std::string	prefixDir_
int	prescaleFactor_
MonitorElement *	RpcEvents
std::string	summaryDir_
bool	useRollInfo_

Static Private Attributes
static const std::string	regions_ [3] = {"EndcapNegative","Barrel","EndcapPositive"}
static const std::string	xLabels_ [7] = {"Good", "OFF", "Nois.St","Nois.Ch","Part.Dead","Dead","Bad.Shape"}

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
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 11 of file RPCChamberQuality.h.

Member Enumeration Documentation

enum RPCChamberQuality::chamberQualityState

private

Enumerator
GoodState
OffState
NoisyStripState
NoisyRollState
PartiallyDeadState
DeadState
BadShapeState

Definition at line 30 of file RPCChamberQuality.h.

{ GoodState= 1 , OffState =2, NoisyStripState= 3, NoisyRollState= 4 , PartiallyDeadState=5 , DeadState=6,BadShapeState=7 };

Constructor & Destructor Documentation

RPCChamberQuality::RPCChamberQuality

(

const edm::ParameterSet &

ps

)

Definition at line 13 of file RPCChamberQuality.cc.

References enableDQMClients_, edm::ParameterSet::getUntrackedParameter(), init_, lumiCounter_, minEvents, numberOfDisks_, offlineDQM_, prefixDir_, prescaleFactor_, AlCaHLTBitMon_QueryRunRegistry::string, summaryDir_, and useRollInfo_.

                                                              {
  edm::LogVerbatim ("rpcchamberquality") << "[RPCChamberQuality]: Constructor";
  
  prescaleFactor_ =  ps.getUntrackedParameter<int>("PrescaleFactor", 5);

  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");

  summaryDir_ =  subsystemFolder +"/"+  recHitTypeFolder +"/"+ summaryFolder ;
  prefixDir_ =  subsystemFolder +"/"+  recHitTypeFolder  ;

  enableDQMClients_ = ps.getUntrackedParameter<bool> ("EnableRPCDqmClient",true); 

  minEvents = ps.getUntrackedParameter<int>("MinimumRPCEvents", 10000);
  numberOfDisks_ = ps.getUntrackedParameter<int>("NumberOfEndcapDisks", 4);
  useRollInfo_ = ps.getUntrackedParameter<bool> ("UseRollInfo",false); 
  offlineDQM_ = ps.getUntrackedParameter<bool> ("OfflineDQM",true); 

  init_ = false; 
  lumiCounter_ = 0;
}

RPCChamberQuality::~RPCChamberQuality ( )

virtual

Definition at line 36 of file RPCChamberQuality.cc.

                                     {
  edm::LogVerbatim ("rpcchamberquality") << "[RPCChamberQuality]: Destructor ";
}

Member Function Documentation

void RPCChamberQuality::beginJob ( void  )

protectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 40 of file RPCChamberQuality.cc.

                                {
  edm::LogVerbatim ("rpcchamberquality") << "[RPCChamberQuality]: Begin job ";
}

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

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 68 of file RPCChamberQuality.cc.

References enableDQMClients_, fillMonitorElements(), init_, myBooker(), and offlineDQM_.

                                                                                   {

  edm::LogVerbatim ("rpcchamberquality") <<"[RPCChamberQuality]: End Job";
  
  if( enableDQMClients_ ) {
    edm::LogVerbatim ("rpcchamberquality") <<"[RPCChamberQuality]: Performing DQM client operation";
    if( offlineDQM_ && !init_ ) {this->myBooker(ibooker);}
    if(init_){this->fillMonitorElements(igetter);}
  }
}

void RPCChamberQuality::dqmEndLuminosityBlock ( DQMStore::IBooker &  ibooker, DQMStore::IGetter &  igetter, edm::LuminosityBlock const &  iLumi, edm::EventSetup const &  iSetup  )

protectedvirtual

Reimplemented from DQMEDHarvester.

Definition at line 46 of file RPCChamberQuality.cc.

References enableDQMClients_, fillMonitorElements(), lumiCounter_, myBooker(), offlineDQM_, and prescaleFactor_.

                                                                                                                                                                {

  edm::LogVerbatim ("rpcchamberquality") <<"[RPCChamberQuality]: End Lumi Block";
 
  if( enableDQMClients_  &&  !offlineDQM_ ) {
 
    edm::LogVerbatim ("rpcchamberquality") <<"[RPCChamberQuality]: Performing DQM client operation";  

    if (lumiCounter_==0){ //book histograms during at firt LB end
      this->myBooker(ibooker); 
      lumiCounter_ = prescaleFactor_;
    }
    
    if (lumiCounter_%prescaleFactor_ == 0){//Fill histograms every N=prescaleFactor_ LBs
      this->fillMonitorElements(igetter);
    }
  }
  lumiCounter_++;
}

void RPCChamberQuality::fillMonitorElements ( DQMStore::IGetter &  igetter )

private

Definition at line 145 of file RPCChamberQuality.cc.

References python.tagInventory::entries, DQMStore::IGetter::get(), MonitorElement::getBinContent(), MonitorElement::getEntries(), i, init_, minEvents, NULL, numberOfDisks_, performeClientOperation(), prefixDir_, alignCSCRings::r, regions_, MonitorElement::Reset(), RpcEvents, MonitorElement::setBinContent(), edmLumisInFiles::summary, summaryDir_, and x.

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

                                                                    {

  std::stringstream meName;
   
  meName.str("");
  meName<<prefixDir_<<"/RPCEvents"; 
  int rpcEvents=minEvents;
  RpcEvents = igetter.get(meName.str());
  
  if(RpcEvents) rpcEvents= (int)RpcEvents->getBinContent(1);
  
  if(rpcEvents >= minEvents){

    init_ = true;
    
    MonitorElement * summary[3];
    
    for(int r = 0 ; r < 3 ; r++) {    
      meName.str("");
      meName<<summaryDir_<<"/RPCChamberQuality_"<<RPCChamberQuality::regions_[r]; 
      summary[r] = igetter.get(meName.str());
      
      if( summary[r] != 0 ) summary[r]->Reset();
    }
    
    //Barrel
    for (int wheel=-2; wheel<3; wheel++) { // loop by Wheels
      meName.str("");
      meName<<"Roll_vs_Sector_Wheel"<<wheel;
      
      this->performeClientOperation(meName.str(), 0 , summary[1],igetter);
    } // loop by Wheels
    
    
    // Endcap
    for(int i=-numberOfDisks_; i<=numberOfDisks_; i++) {//loop on Disks
      if(i==0) continue;
      
      meName.str("");
      meName<<"Ring_vs_Segment_Disk"<<i;

      if(i<0) this->performeClientOperation(meName.str(), -1 , summary[0],igetter);
      else this->performeClientOperation(meName.str(), 1 , summary[2],igetter);
    }//loop on Disks
      
    MonitorElement * RpcOverview = NULL;
    meName.str("");
    meName<<summaryDir_<<"/RPC_System_Quality_Overview"; 
    RpcOverview = igetter.get(meName.str());
    RpcOverview->Reset();

    if(RpcOverview) {//Fill Overview ME
      for(int r = 0 ; r< 3; r++) {
    if (summary[r] == 0 ) continue;
    double entries = summary[r]->getEntries();
    if(entries == 0) continue;
    for (int x = 1; x <= 7; x++) {
      RpcOverview->setBinContent(x,r+1,(summary[r]->getBinContent(x)/entries));
       } 
      } 
    } //loop by LimiBloks
  }
} 

void RPCChamberQuality::myBooker ( DQMStore::IBooker &  ibooker )

private

Definition at line 81 of file RPCChamberQuality.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), ztail::d, init_, rpcdqm::utils::labelXAxisSector(), rpcdqm::utils::labelXAxisSegment(), rpcdqm::utils::labelYAxisRing(), rpcdqm::utils::labelYAxisRoll(), numberOfDisks_, alignCSCRings::r, regions_, MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), summaryDir_, useRollInfo_, w, x, and xLabels_.

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

                                                          {

  init_ = true;
 
  ibooker.setCurrentFolder(summaryDir_);

  MonitorElement* me;
  std::stringstream histoName;
  
  rpcdqm::utils rpcUtils;

  for (int r = 0 ; r < 3; r++){

    histoName.str("");
    histoName<<"RPCChamberQuality_"<<regions_[r]; 
    me = ibooker.book1D(histoName.str().c_str(), histoName.str().c_str(),  7, 0.5, 7.5);
    
    for (int x = 1; x <8 ; x++) {me->setBinLabel(x, xLabels_[x-1]);}
  }


  histoName.str("");
  histoName<<"RPC_System_Quality_Overview"; 
  me = ibooker.book2D(histoName.str().c_str(), histoName.str().c_str(),  7, 0.5, 7.5, 3, 0.5, 3.5);
  me->setBinLabel(1, "E+", 2);
  me->setBinLabel(2, "B", 2);
  me->setBinLabel(3, "E-", 2);
    
  for (int x = 1; x <8 ; x++) {me->setBinLabel(x, xLabels_[x-1]);}
  
  for(int w=-2; w<3;w++){//Loop on wheels
    
    histoName.str("");
    histoName<<"RPCChamberQuality_Roll_vs_Sector_Wheel"<<w;    
    me = ibooker.book2D(histoName.str().c_str(), histoName.str().c_str(),  12, 0.5, 12.5, 21, 0.5, 21.5);

    rpcUtils.labelXAxisSector( me);
    rpcUtils.labelYAxisRoll(me, 0, w, useRollInfo_ );

    histoName.str("");
    histoName<<"RPCChamberQuality_Distribution_Wheel"<<w;    
    me = ibooker.book1D(histoName.str().c_str(), histoName.str().c_str(),  7, 0.5, 7.5);

    for (int x = 1; x <8; x++) {me->setBinLabel(x, xLabels_[x-1]);}        
  }//end loop on wheels

  for(int d= -numberOfDisks_; d<= numberOfDisks_ ; d++) { // Loop on disk
    if(d==0) continue; 
      histoName.str("");
      histoName<<"RPCChamberQuality_Ring_vs_Segment_Disk"<<d;       //  2D histo for RPC Qtest
      me = ibooker.book2D(histoName.str().c_str(), histoName.str().c_str(),  36, 0.5, 36.5, 6, 0.5, 6.5);
      rpcUtils.labelXAxisSegment(me);
      rpcUtils.labelYAxisRing(me, 2, useRollInfo_ );
      
      histoName.str("");
      histoName<<"RPCChamberQuality_Distribution_Disk"<<d;    
      me = ibooker.book1D(histoName.str().c_str(), histoName.str().c_str(),  7, 0.5, 7.5);
      
      for (int x = 1; x <8 ; x++) {me->setBinLabel(x, xLabels_[x-1]);}
  } 
}

void RPCChamberQuality::performeClientOperation ( std::string  MESufix, int  region, MonitorElement *  quality, DQMStore::IGetter &  igetter  )

private

Definition at line 214 of file RPCChamberQuality.cc.

References BadShapeState, DeadState, MonitorElement::Fill(), HcalObjRepresent::Fill(), DQMStore::IGetter::get(), MonitorElement::getBinContent(), GoodState, AnalysisDataFormats_SUSYBSMObjects::hv, VarParsing::mult, NoisyRollState, NoisyStripState, NULL, OffState, PartiallyDeadState, MonitorElement::Reset(), summaryDir_, x, and detailsBasic3DVector::y.

Referenced by fillMonitorElements().

                                                                                                                                 {


  MonitorElement * RCQ=NULL;  
  MonitorElement * RCQD=NULL; 
  
  MonitorElement * DEAD=NULL;  
  MonitorElement * CLS=NULL;
  MonitorElement * MULT=NULL;
  MonitorElement * NoisySt=NULL;
  MonitorElement * Chip=NULL;
  MonitorElement * HV=NULL;
  MonitorElement * LV=NULL;
  std::stringstream meName; 

  meName.str("");
  meName<<summaryDir_<<"/RPCChamberQuality_"<<MESufix; 
  RCQ = igetter.get(meName.str());
  //  if (RCQ)  RCQ->Reset();


  int pos = MESufix.find_last_of("_");
  meName.str("");
  meName<<summaryDir_<<"/RPCChamberQuality_Distribution"<<MESufix.substr(pos); 
  RCQD = igetter.get(meName.str());
 if (RCQD) RCQD->Reset();
  
  //get HV Histo
  meName.str("");                        
  meName<<summaryDir_<<"/HVStatus_"<<MESufix;
  HV = igetter.get(meName.str());   
  //get LV Histo
  meName.str("");                     
  meName<<summaryDir_<<"/LVStatus_"<<MESufix; 
  LV = igetter.get(meName.str());
  //Dead 
  meName.str("");
  meName << summaryDir_<<"/DeadChannelFraction_"<<MESufix;
  DEAD = igetter.get(meName.str());
  //ClusterSize
  meName.str("");
  meName<<summaryDir_<<"/ClusterSizeIn1Bin_"<<MESufix;
  CLS = igetter.get(meName.str());
  //NoisyStrips
  meName.str("");
  meName<<summaryDir_<<"/RPCNoisyStrips_"<<MESufix;
  NoisySt = igetter.get(meName.str());
  //Multiplicity
  meName.str("");
  meName<<summaryDir_<<"/NumberOfDigi_Mean_"<<MESufix;
  MULT = igetter.get(meName.str());
  //Asymetry
  meName.str("");
  meName<<summaryDir_<<"/AsymmetryLeftRight_"<<MESufix;
  Chip = igetter.get(meName.str());        
  
  int xBinMax, yBinMax;    

  if (region != 0) xBinMax = 37;
  else xBinMax = 13;

  for(int x=1; x<xBinMax; x++) {
    if (region != 0 )  {
      yBinMax = 7;
    }else {  
      if(x==4) yBinMax=22;
      else if(x==9 || x==11) yBinMax=16;
      else yBinMax=18;
    }  
    for(int y=1; y<yBinMax; y++) {
      int hv=1;
      int lv=1;
      float dead =0;
      float firstbin= 0;
      float noisystrips = 0;
      float mult = 0;
      float asy = 0;
      chamberQualityState chamberState = GoodState;
       
      if(HV) hv = (int)HV ->getBinContent(x,y);
      if(LV) lv = (int)LV ->getBinContent(x,y);
        
      if( hv!=1 || lv!=1) { 
    chamberState = OffState;
      }else {
    if(DEAD) dead= DEAD -> getBinContent(x,y);
    if(dead>=0.80 ) {  
      chamberState = DeadState;
    }else if (0.33<=dead && dead<0.80 ){
      chamberState = PartiallyDeadState;
    }else {        
      if(CLS ) firstbin = CLS -> getBinContent(x,y);
      if(firstbin >= 0.88) {
        chamberState = NoisyStripState;
      } else {   
        if(NoisySt)  noisystrips = NoisySt -> getBinContent(x,y);
        if (noisystrips > 0){ 
          chamberState = NoisyStripState;
        }else {  
          if(MULT) mult = MULT -> getBinContent(x,y);
          if(mult>=6) {  
        chamberState = NoisyRollState;
          }else {  
        if (Chip) asy = Chip->getBinContent(x,y);
        if(asy>0.35) {  
          chamberState  = BadShapeState;
        }else {  
          chamberState  = GoodState;
        }
          }
        } 
      }
    }
      }
      if (RCQ)  RCQ -> setBinContent(x,y, chamberState);
      if (RCQD)   RCQD -> Fill(chamberState); 
      if (quality)   quality->Fill(chamberState); 
    }
  }
  return;
} 

Member Data Documentation

bool RPCChamberQuality::enableDQMClients_

private

Definition at line 34 of file RPCChamberQuality.h.

Referenced by dqmEndJob(), dqmEndLuminosityBlock(), and RPCChamberQuality().

bool RPCChamberQuality::init_

private

Definition at line 46 of file RPCChamberQuality.h.

Referenced by dqmEndJob(), fillMonitorElements(), myBooker(), and RPCChamberQuality().

int RPCChamberQuality::lumiCounter_

private

Definition at line 32 of file RPCChamberQuality.h.

Referenced by dqmEndLuminosityBlock(), and RPCChamberQuality().

int RPCChamberQuality::minEvents

private

Definition at line 48 of file RPCChamberQuality.h.

Referenced by fillMonitorElements(), and RPCChamberQuality().

int RPCChamberQuality::numberOfDisks_

private

Definition at line 44 of file RPCChamberQuality.h.

Referenced by fillMonitorElements(), myBooker(), and RPCChamberQuality().

int RPCChamberQuality::numLumBlock_

private

Definition at line 49 of file RPCChamberQuality.h.

bool RPCChamberQuality::offlineDQM_

private

Definition at line 35 of file RPCChamberQuality.h.

Referenced by dqmEndJob(), dqmEndLuminosityBlock(), and RPCChamberQuality().

std::string RPCChamberQuality::prefixDir_

private

Definition at line 39 of file RPCChamberQuality.h.

Referenced by fillMonitorElements(), and RPCChamberQuality().

int RPCChamberQuality::prescaleFactor_

private

Definition at line 43 of file RPCChamberQuality.h.

Referenced by dqmEndLuminosityBlock(), and RPCChamberQuality().

const std::string RPCChamberQuality::regions_ = {"EndcapNegative","Barrel","EndcapPositive"}

staticprivate

Definition at line 41 of file RPCChamberQuality.h.

Referenced by fillMonitorElements(), and myBooker().

MonitorElement* RPCChamberQuality::RpcEvents

private

Definition at line 29 of file RPCChamberQuality.h.

Referenced by fillMonitorElements().

std::string RPCChamberQuality::summaryDir_

private

Definition at line 39 of file RPCChamberQuality.h.

Referenced by fillMonitorElements(), myBooker(), performeClientOperation(), and RPCChamberQuality().

bool RPCChamberQuality::useRollInfo_

private

Definition at line 42 of file RPCChamberQuality.h.

Referenced by myBooker(), and RPCChamberQuality().

const std::string RPCChamberQuality::xLabels_ = {"Good", "OFF", "Nois.St","Nois.Ch","Part.Dead","Dead","Bad.Shape"}

staticprivate

Definition at line 40 of file RPCChamberQuality.h.

Referenced by myBooker().