#include <RPCEfficiencyPerRingLayer.h>

Inheritance diagram for RPCEfficiencyPerRingLayer:

Public Member Functions
	RPCEfficiencyPerRingLayer (const edm::ParameterSet &iConfig)
	Constructor. More...

virtual	~RPCEfficiencyPerRingLayer ()
	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 ()

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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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::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 Attributes
MonitorElement *	EfficiencyPerLayer

MonitorElement *	EfficiencyPerRing

std::string	globalFolder_

int	innermostRings_

int	numberOfDisks_

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 10 of file RPCEfficiencyPerRingLayer.h.

Constructor & Destructor Documentation

RPCEfficiencyPerRingLayer::RPCEfficiencyPerRingLayer ( const edm::ParameterSet & iConfig )

Constructor.

Definition at line 13 of file RPCEfficiencyPerRingLayer.cc.

References edm::ParameterSet::getUntrackedParameter(), globalFolder_, innermostRings_, numberOfDisks_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                               {
 
   globalFolder_ = ps.getUntrackedParameter<std::string>("GlobalFolder", "RPC/RPCEfficiency/");
   numberOfDisks_ = ps.getUntrackedParameter<int>("NumberOfEndcapDisks", 4);
   innermostRings_ = ps.getUntrackedParameter<int>("NumberOfInnermostEndcapRings", 2);
 
 }

RPCEfficiencyPerRingLayer::~RPCEfficiencyPerRingLayer ( )

virtual

Destructor.

Definition at line 22 of file RPCEfficiencyPerRingLayer.cc.

22 { }

Member Function Documentation

void RPCEfficiencyPerRingLayer::beginJob ( void )

protectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 24 of file RPCEfficiencyPerRingLayer.cc.

24 {}

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

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 28 of file RPCEfficiencyPerRingLayer.cc.

References funct::abs(), DQMStore::IBooker::book1D(), EfficiencyPerLayer, EfficiencyPerRing, DQMStore::IGetter::get(), MonitorElement::getBinContent(), MonitorElement::getTH1(), globalFolder_, i, innermostRings_, j, relval_2017::k, NULL, numberOfDisks_, MonitorElement::setBinContent(), MonitorElement::setBinError(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), and mathSSE::sqrt().

                                                                                               {
 
 
   ibooker.setCurrentFolder(globalFolder_);
   EfficiencyPerRing = ibooker.book1D("EfficiencyPerRing","Efficiency per Ring",12,0.5,12.5);
   EfficiencyPerLayer = ibooker.book1D("EfficiencyPerLayer","Efficiency per Layer",6,0.5,6.5);
 
 
 
    std::stringstream meName1;
    MonitorElement * myMe1;
 
    std::stringstream meName2;
    MonitorElement * myMe2;
 
    std::stringstream binName;
    binName.str("");
 
    meName1.str("");
    meName1<<globalFolder_;
 
    meName2.str("");
    meName2<<globalFolder_;
 
    int k = 0;
 
    //EfficiencyPerRing
 
    for(int i=(-1*numberOfDisks_ ); i<=numberOfDisks_ ; i++){
     if(i == 0) continue;
     for(int j=innermostRings_; j<=3; j++){
         meName1.str("");
         meName2.str("");
         if(i < 0){
                 meName1<<globalFolder_<<"Azimutal/ExGregDistroDm"<<abs(i)<<"R"<<j;
                 meName2<<globalFolder_<<"Azimutal/OcGregDistroDm"<<abs(i)<<"R"<<j;
             //std::cout<<meName1.str()<<std::endl;
             //std::cout<<meName2.str()<<std::endl;
             binName.str("");
             binName<<"RE-"<<i<<"/"<<j;
             k++;}
         else if(i >0){
                 meName1<<globalFolder_<<"Azimutal/ExGregDistroD"<<abs(i)<<"R"<<j;
                 meName2<<globalFolder_<<"Azimutal/OcGregDistroD"<<abs(i)<<"R"<<j;
             //std::cout<<meName1.str()<<std::endl;
             //std::cout<<meName2.str()<<std::endl;
             binName.str("");
             binName<<"RE+"<<abs(i)<<"/"<<j;
             k++;}
 
         myMe1 = NULL; myMe1 = NULL;
             myMe1 = igetter.get(meName1.str());
             myMe2 = igetter.get(meName2.str());
 
         if(myMe1 && myMe2){
           TH1 * histo1 = myMe1->getTH1();
           TH1 * histo2 = myMe2->getTH1();
           
           int exg = histo1->Integral(); 
           int obg = histo2->Integral();
           double eff = 0;
           double err = 0;
           if(exg != 0){ 
             eff = (double)obg/(double)exg; 
             err = sqrt(eff*(1-eff)/(double)exg);}
           
           if(k%2 != 0 && k < 7){ 
             EfficiencyPerRing->setBinContent(k+1, eff);
             EfficiencyPerRing->setBinError(k+1, err);
             EfficiencyPerRing->setBinLabel(k+1, binName.str());}
           else if(k%2 == 0 && k < 7){
             EfficiencyPerRing->setBinContent(k-1, eff);
             EfficiencyPerRing->setBinError(k-1, err);
             EfficiencyPerRing->setBinLabel(k-1, binName.str());}
           else if(k >= 7){
             EfficiencyPerRing->setBinContent(k, eff);
             EfficiencyPerRing->setBinError(k, err);
             EfficiencyPerRing->setBinLabel(k, binName.str());}
         }
     }
    }
 
    //EfficiencyPerLayer
    
    for(int i=1; i<=6; i++){
      
     int Exp = 0;
     int Obs = 0;
     double eff = 0;
     double err = 0;
 
     for(int j=-2; j<=2; j++){
         meName1.str("");
         meName2.str("");
             if(j < 0){ 
             meName1<<globalFolder_<<"BarrelPerLayer/ExpLayerWm"<<abs(j);
                 meName2<<globalFolder_<<"BarrelPerLayer/ObsLayerWm"<<abs(j);}
         else{
             meName1<<globalFolder_<<"BarrelPerLayer/ExpLayerW"<<j;
                 meName2<<globalFolder_<<"BarrelPerLayer/ObsLayerW"<<j;}
         myMe1 = NULL; myMe1 = NULL;
         myMe1 = igetter.get(meName1.str());
             myMe2 = igetter.get(meName2.str());
         if(myMe1 && myMe2){
             Exp += myMe1->getBinContent(i);
             Obs += myMe2->getBinContent(i);
         }
     }
 
     if(Exp != 0){
         eff = (double)Obs/(double)Exp;
         err = sqrt(eff*(1-eff)/(double)Exp);}
     EfficiencyPerLayer->setBinContent(i, eff);
     EfficiencyPerLayer->setBinError(i, err);
     binName.str("");
     binName<<"Layer "<<i;
     EfficiencyPerLayer->setBinLabel(i, binName.str());
 
    }
 
 }