#include <RPCDeadChannelTest.h>

Inheritance diagram for RPCDeadChannelTest:

Public Member Functions
void	beginJob (std::string &)

void	clientOperation ()

void	getMonitorElements (std::vector< MonitorElement * > &, std::vector< RPCDetId > &, std::string &)

void	myBooker (DQMStore::IBooker &)

	RPCDeadChannelTest (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~RPCDeadChannelTest ()
	Destructor. More...

Public Member Functions inherited from RPCClient
virtual	~RPCClient (void)

Private Attributes
MonitorElement *	DEADDisk [10]

MonitorElement *	DEADWheel [5]

std::string	globalFolder_

std::vector< RPCDetId >	myDetIds_

std::vector< MonitorElement * >	myOccupancyMe_

int	numberOfDisks_

int	numberOfRings_

int	prescaleFactor_

bool	useRollInfo_

Detailed Description

Definition at line 10 of file RPCDeadChannelTest.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 13 of file RPCDeadChannelTest.cc.

References edm::ParameterSet::getUntrackedParameter(), numberOfDisks_, numberOfRings_, prescaleFactor_, and useRollInfo_.

                                                                 {
  
   edm::LogVerbatim ("rpcdeadchanneltest") << "[RPCDeadChannelTest]: Constructor";
 
   useRollInfo_ = ps.getUntrackedParameter<bool>("UseRollInfo", false);
 
   prescaleFactor_ = ps.getUntrackedParameter<int>("DiagnosticPrescale", 1);
   numberOfDisks_ = ps.getUntrackedParameter<int>("NumberOfEndcapDisks", 4);
   numberOfRings_ = ps.getUntrackedParameter<int>("NumberOfEndcapRings", 2);
 }

RPCDeadChannelTest::~RPCDeadChannelTest ( )

virtual

Destructor.

Definition at line 24 of file RPCDeadChannelTest.cc.

24 {}

Member Function Documentation

void RPCDeadChannelTest::beginJob ( std::string & workingFolder )

virtual

Implements RPCClient.

Definition at line 26 of file RPCDeadChannelTest.cc.

References globalFolder_.

                                                            {
  edm::LogVerbatim ("rpcdeadchanneltest") << "[RPCDeadChannelTest]: Begin Job";
   globalFolder_ =  workingFolder;
 
 }

void RPCDeadChannelTest::clientOperation ( )

virtual

Implements RPCClient.

Definition at line 47 of file RPCDeadChannelTest.cc.

References DEADDisk, DEADWheel, rpcdqm::utils::detId2RollNr(), MonitorElement::getBinContent(), MonitorElement::getNbinsX(), MonitorElement::getQReport(), QReport::getQTresult(), mps_fire::i, myDetIds_, myOccupancyMe_, NULL, numberOfDisks_, numberOfRings_, RPCDetId::region(), RPCDetId::ring(), RPCDetId::roll(), RPCDetId::sector(), RPCGeomServ::segment(), MonitorElement::setBinContent(), RPCDetId::station(), x, and cuy::yBin.

                                         {
  
   edm::LogVerbatim ("rpcdeadchanneltest") <<"[RPCDeadChannelTest]:Client Operation";
 
 
   MonitorElement * DEAD = NULL;
  
   //Loop on chambers
     for (unsigned int  i = 0 ; i<myOccupancyMe_.size();i++){
       
       RPCDetId & detId = myDetIds_[i];
       MonitorElement * myMe = myOccupancyMe_[i];
 
       if (! myMe ) continue;
 
       const QReport * theOccupancyQReport = myMe->getQReport("DeadChannel_0"); 
  
       float deadFraction = 0.0 ;
 
       if(theOccupancyQReport) {
     
     float qtresult = theOccupancyQReport->getQTresult();
     // std::vector<dqm::me_util::Channel> badChannels = theOccupancyQReport->getBadChannels();
     deadFraction = 1.0 - qtresult;
 
       }else{
     int xBins = myMe->getNbinsX();
     float emptyBins = 0.0;
     for(int x = 1 ; x<= xBins ; x++){if(myMe->getBinContent(x) == 0 ) {emptyBins++;}}
     if (xBins != 0){    deadFraction = emptyBins/xBins;}
       }
      
      if (detId.region()==0)   DEAD = DEADWheel[detId.ring() + 2] ;
      else{
        if(-detId.station()+ numberOfDisks_ >= 0 ){
      
      if(detId.region()<0){
        DEAD  = DEADDisk[-detId.station() + numberOfDisks_];
      }else{
        DEAD = DEADDisk[detId.station() + numberOfDisks_-1];
      }
        }
      }
 
      if (DEAD){
        int xBin,yBin;
        if(detId.region()==0){//Barrel
      xBin= detId.sector();
      rpcdqm::utils rollNumber;
      yBin = rollNumber.detId2RollNr(detId);
        }else{//Endcap
      //get segment number
      RPCGeomServ RPCServ(detId);
      xBin = RPCServ.segment();
      (numberOfRings_ == 3 ? yBin= detId.ring()*3-detId.roll()+1 : yBin= (detId.ring()-1)*3-detId.roll()+1);
      }
        DEAD->setBinContent(xBin,yBin, deadFraction );
 
      }
 
     }//End loop on rolls in given chambers
 
 }

void RPCDeadChannelTest::getMonitorElements	(	std::vector< MonitorElement * > &	meVector,
		std::vector< RPCDetId > &	detIdVector,
		std::string &	clientHistoName
	)

virtual

Implements RPCClient.

Definition at line 33 of file RPCDeadChannelTest.cc.

References mps_fire::i, myDetIds_, myOccupancyMe_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                                                                  {
 
   for (unsigned int i = 0 ; i<meVector.size(); i++){
     
     std::string meName =  meVector[i]->getName();
 
     if(meName.find(clientHistoName) != std::string::npos){
       myOccupancyMe_.push_back(meVector[i]);
       myDetIds_.push_back(detIdVector[i]);
     }  
   }
 }

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

virtual

Implements RPCClient.

Definition at line 111 of file RPCDeadChannelTest.cc.

References DQMStore::IBooker::book2D(), DEADDisk, DEADWheel, globalFolder_, mps_fire::i, rpcdqm::utils::labelXAxisSector(), rpcdqm::utils::labelXAxisSegment(), rpcdqm::utils::labelYAxisRing(), rpcdqm::utils::labelYAxisRoll(), MessageLogger_cfi::limit, numberOfDisks_, numberOfRings_, PFRecoTauDiscriminationByIsolation_cfi::offset, MonitorElement::setBinContent(), DQMStore::IBooker::setCurrentFolder(), useRollInfo_, x, and y.

                                                           {
 
   ibooker.setCurrentFolder( globalFolder_);
   
   std::stringstream histoName;
   
   rpcdqm::utils rpcUtils;
   
   int limit = numberOfDisks_;
   if(numberOfDisks_ < 2) limit = 2;
   
   for (int i = -1 * limit; i<= limit;i++ ){//loop on wheels and disks
     if (i>-3 && i<3){//wheels
       histoName.str("");
       histoName<<"DeadChannelFraction_Roll_vs_Sector_Wheel"<<i;
       DEADWheel[i+2] = ibooker.book2D(histoName.str().c_str(), histoName.str().c_str(), 12, 0.5, 12.5, 21, 0.5, 21.5);
       
       for (int x = 1; x<=12; x++){
     for(int y=1; y<=21; y++){
       DEADWheel[i+2]->setBinContent(x,y,-1);
     }
       }
 
       rpcUtils.labelXAxisSector( DEADWheel[i+2]);
       rpcUtils.labelYAxisRoll( DEADWheel[i+2], 0, i, useRollInfo_);
     }//end wheels
     
     if (i == 0  || i > numberOfDisks_ || i< (-1 * numberOfDisks_)){continue;}
     
     int offset = numberOfDisks_;
     if (i>0) {offset --;} //used to skip case equale to zero
     
     histoName.str("");
     histoName<<"DeadChannelFraction_Ring_vs_Segment_Disk"<<i;
     DEADDisk[i+offset] = ibooker.book2D(histoName.str().c_str(), histoName.str().c_str(),36, 0.5, 36.5, 3*numberOfRings_, 0.5,3*numberOfRings_+ 0.5);
     
     rpcUtils.labelXAxisSegment(DEADDisk[i+offset]);
     rpcUtils.labelYAxisRing(DEADDisk[i+offset], numberOfRings_ ,useRollInfo_);
         
   }//end loop on wheels and disks
   
   
 }