RPCEfficiencyPerRingLayer Class Reference

#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	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override final
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) override
	DQMEDHarvester (void)
virtual void	endJob () override final
virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override final
virtual void	endRun (edm::Run const &, edm::EventSetup const &) override
virtual void	produce (edm::Event &, edm::EventSetup const &) override final
virtual	~DQMEDHarvester ()=default
Public Member Functions inherited from edm::one::EDProducer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources, edm::EndLuminosityBlockProducer >
	EDProducer ()=default
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	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 Attributes
MonitorElement *	EfficiencyPerLayer
MonitorElement *	EfficiencyPerRing
std::string	globalFolder_
int	innermostRings_
int	numberOfDisks_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
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.

{ }

Member Function Documentation

void RPCEfficiencyPerRingLayer::beginJob ( void  )

overrideprotectedvirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 24 of file RPCEfficiencyPerRingLayer.cc.

{}

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_, mps_fire::i, innermostRings_, gen::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());

   }

}

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

overrideprotectedvirtual

Reimplemented from DQMEDHarvester.

Definition at line 26 of file RPCEfficiencyPerRingLayer.cc.

{}

Member Data Documentation

MonitorElement* RPCEfficiencyPerRingLayer::EfficiencyPerLayer

private

Definition at line 30 of file RPCEfficiencyPerRingLayer.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencyPerRingLayer::EfficiencyPerRing

private

Definition at line 29 of file RPCEfficiencyPerRingLayer.h.

Referenced by dqmEndJob().

std::string RPCEfficiencyPerRingLayer::globalFolder_

private

Definition at line 35 of file RPCEfficiencyPerRingLayer.h.

Referenced by dqmEndJob(), and RPCEfficiencyPerRingLayer().

int RPCEfficiencyPerRingLayer::innermostRings_

private

Definition at line 33 of file RPCEfficiencyPerRingLayer.h.

Referenced by dqmEndJob(), and RPCEfficiencyPerRingLayer().

int RPCEfficiencyPerRingLayer::numberOfDisks_

private

Definition at line 32 of file RPCEfficiencyPerRingLayer.h.

Referenced by dqmEndJob(), and RPCEfficiencyPerRingLayer().