ECALpedestalPCLworker Class Reference

#include <ECALpedestalPCLworker.h>

Inheritance diagram for ECALpedestalPCLworker:

Public Member Functions
	ECALpedestalPCLworker (const edm::ParameterSet &)
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	endJob () override

Static Private Member Functions
static bool	adc_compare (uint16_t a, uint16_t b)

Private Attributes
bool	checkSignal_
edm::EDGetTokenT< EBDigiCollection >	digiTokenEB_
edm::EDGetTokenT< EEDigiCollection >	digiTokenEE_
std::string	dqmDir_
bool	dynamicBooking_
int	fixedBookingCenterBin_
std::vector< MonitorElement * >	meEB_
std::vector< MonitorElement * >	meEE_
int	nBins_
uint32_t	pedestalSamples_
bool	requireStableBeam_
uint32_t	sThresholdEB_
uint32_t	sThresholdEE_
edm::EDGetTokenT< TCDSRecord >	tcdsToken_

Detailed Description

Description: Fill DQM histograms with pedestals. Intended to be used on laser data from the TestEnablesEcalHcal dataset

Definition at line 35 of file ECALpedestalPCLworker.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file ECALpedestalPCLworker.cc.

References checkSignal_, digiTokenEB_, digiTokenEE_, dqmDir_, dynamicBooking_, fixedBookingCenterBin_, edm::ParameterSet::getParameter(), nBins_, pedestalSamples_, requireStableBeam_, sThresholdEB_, sThresholdEE_, AlCaHLTBitMon_QueryRunRegistry::string, and tcdsToken_.

{
    edm::InputTag digiTagEB= iConfig.getParameter<edm::InputTag>("BarrelDigis");
    edm::InputTag digiTagEE= iConfig.getParameter<edm::InputTag>("EndcapDigis");

    digiTokenEB_ = consumes<EBDigiCollection>(digiTagEB);
    digiTokenEE_ = consumes<EEDigiCollection>(digiTagEE);

    pedestalSamples_ = iConfig.getParameter<uint32_t>("pedestalSamples");
    checkSignal_     = iConfig.getParameter<bool>("checkSignal");
    sThresholdEB_      = iConfig.getParameter<uint32_t>("sThresholdEB");
    sThresholdEE_      = iConfig.getParameter<uint32_t>("sThresholdEE");

    dynamicBooking_ = iConfig.getParameter<bool>("dynamicBooking");
    fixedBookingCenterBin_ = iConfig.getParameter<int>("fixedBookingCenterBin");
    nBins_          = iConfig.getParameter<int>("nBins");
    dqmDir_         = iConfig.getParameter<std::string>("dqmDir");

    edm::InputTag tcdsRecord= iConfig.getParameter<edm::InputTag>("tcdsRecord");
    tcdsToken_       = consumes<TCDSRecord>(tcdsRecord);
    requireStableBeam_ = iConfig.getParameter<bool>("requireStableBeam");
}

Member Function Documentation

static bool ECALpedestalPCLworker::adc_compare ( uint16_t  a, uint16_t  b  )

inlinestaticprivate

Definition at line 72 of file ECALpedestalPCLworker.h.

Referenced by analyze().

{ return ( a&0xFFF )  < (b&0xFFF);}

void ECALpedestalPCLworker::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 42 of file ECALpedestalPCLworker.cc.

References adc_compare(), edm::DataFrame::begin(), edm::DataFrameContainer::begin(), checkSignal_, digiTokenEB_, digiTokenEE_, edm::DataFrame::end(), edm::DataFrameContainer::end(), HcalObjRepresent::Fill(), EcalDataFrame::frame(), BSTRecord::getBeamMode(), TCDSRecord::getBST(), edm::Event::getByToken(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), meEB_, meEE_, pedestalSamples_, requireStableBeam_, sThresholdEB_, sThresholdEE_, and tcdsToken_.

{
    using namespace edm;

    Handle<EBDigiCollection> pDigiEB;
    iEvent.getByToken(digiTokenEB_,pDigiEB);

    Handle<EEDigiCollection> pDigiEE;
    iEvent.getByToken(digiTokenEE_,pDigiEE);


    // Only Events with stable beam
 
    if (requireStableBeam_){
        edm::Handle<TCDSRecord> tcdsData;
        iEvent.getByToken(tcdsToken_,tcdsData);
        int beamMode = tcdsData->getBST().getBeamMode();
        if (beamMode != BSTRecord::BeamMode::STABLE) return;
    }

    for (EBDigiCollection::const_iterator pDigi=pDigiEB->begin(); pDigi!=pDigiEB->end(); ++pDigi){

        EBDetId id = pDigi->id();
        uint32_t hashedId = id.hashedIndex();

        EBDataFrame digi( *pDigi );
        
        if (checkSignal_){
            uint16_t maxdiff = *std::max_element(digi.frame().begin(), digi.frame().end(), adc_compare )  -
                               *std::min_element(digi.frame().begin(), digi.frame().end(), adc_compare );
            if ( maxdiff> sThresholdEB_ ) continue; // assume there is signal in this frame
        }

        //for (auto& mgpasample : digi.frame()) meEB_[hashedId]->Fill(mgpasample&0xFFF);
        for (edm::DataFrame::iterator mgpasample = digi.frame().begin();
             mgpasample!=digi.frame().begin()+pedestalSamples_;
             ++mgpasample )
            meEB_[hashedId]->Fill(*mgpasample&0xFFF);

    } // eb digis



    for (EEDigiCollection::const_iterator pDigi=pDigiEE->begin(); pDigi!=pDigiEE->end(); ++pDigi){

        EEDetId id = pDigi->id();
        uint32_t hashedId = id.hashedIndex();

        EEDataFrame digi( *pDigi );

        if (checkSignal_){
            uint16_t maxdiff = *std::max_element(digi.frame().begin(), digi.frame().end(), adc_compare )  -
                               *std::min_element(digi.frame().begin(), digi.frame().end(), adc_compare );
            if ( maxdiff> sThresholdEE_ ) continue; // assume there is signal in this frame
        }

        //for (auto& mgpasample : digi.frame()) meEE_[hashedId]->Fill(mgpasample&0xFFF);
        for (edm::DataFrame::iterator mgpasample = digi.frame().begin();
             mgpasample!=digi.frame().begin()+pedestalSamples_;
             ++mgpasample )
            meEE_[hashedId]->Fill(*mgpasample&0xFFF);

    } // ee digis




}

void ECALpedestalPCLworker::beginJob ( void  )

overrideprivate

Definition at line 113 of file ECALpedestalPCLworker.cc.

{
}

void ECALpedestalPCLworker::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  run, edm::EventSetup const &  es  )

overrideprivate

Definition at line 133 of file ECALpedestalPCLworker.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), EBDetId::detIdFromDenseIndex(), EEDetId::detIdFromDenseIndex(), dqmDir_, dynamicBooking_, EcalCondObjectContainer< T >::find(), fixedBookingCenterBin_, edm::EventSetup::get(), mps_fire::i, createfilelist::int, EBDetId::kSizeForDenseIndexing, EEDetId::kSizeForDenseIndexing, SiStripPI::max, meEB_, meEE_, min(), nBins_, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                            {

    ibooker.cd();
    ibooker.setCurrentFolder(dqmDir_);

    edm::ESHandle<EcalPedestals> peds;
    es.get<EcalPedestalsRcd>().get(peds);
    

    for ( uint32_t i = 0 ; i< EBDetId::kSizeForDenseIndexing; ++i){

        
        ibooker.setCurrentFolder(dqmDir_+"/EB/"+std::to_string(int(i/100)));

        std::string hname = "eb_" + std::to_string(i);
        DetId id = EBDetId::detIdFromDenseIndex(i);
        int centralBin = fixedBookingCenterBin_;
         
        if (dynamicBooking_){
            centralBin =  int ((peds->find(id))->mean_x12) ;
        }

        int min = centralBin - nBins_/2;
        int max = centralBin + nBins_/2;

        meEB_.push_back(ibooker.book1D(hname,hname,nBins_,min,max));
    }

    for ( uint32_t i = 0 ; i< EEDetId::kSizeForDenseIndexing; ++i){

        ibooker.setCurrentFolder(dqmDir_+"/EE/"+std::to_string(int(i/100)));

        std::string hname = "ee_" + std::to_string(i);

        DetId id = EEDetId::detIdFromDenseIndex(i);
        int centralBin = fixedBookingCenterBin_;
         
        if (dynamicBooking_){
            centralBin =  int ((peds->find(id))->mean_x12) ;
        }

        int min = centralBin - nBins_/2;
        int max = centralBin + nBins_/2;

        meEE_.push_back(ibooker.book1D(hname,hname,nBins_,min,max));

    }

}

void ECALpedestalPCLworker::endJob ( void  )

overrideprivate

Definition at line 119 of file ECALpedestalPCLworker.cc.

Referenced by o2olib.O2ORunMgr::executeJob().

{
}

void ECALpedestalPCLworker::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 125 of file ECALpedestalPCLworker.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                  {
    edm::ParameterSetDescription desc;
    desc.setUnknown();
    descriptions.addDefault(desc);
}

Member Data Documentation

bool ECALpedestalPCLworker::checkSignal_

private

Definition at line 58 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and ECALpedestalPCLworker().

edm::EDGetTokenT<EBDigiCollection> ECALpedestalPCLworker::digiTokenEB_

private

Definition at line 50 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and ECALpedestalPCLworker().

edm::EDGetTokenT<EEDigiCollection> ECALpedestalPCLworker::digiTokenEE_

private

Definition at line 51 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and ECALpedestalPCLworker().

std::string ECALpedestalPCLworker::dqmDir_

private

Definition at line 67 of file ECALpedestalPCLworker.h.

Referenced by bookHistograms(), and ECALpedestalPCLworker().

bool ECALpedestalPCLworker::dynamicBooking_

private

Definition at line 64 of file ECALpedestalPCLworker.h.

Referenced by bookHistograms(), and ECALpedestalPCLworker().

int ECALpedestalPCLworker::fixedBookingCenterBin_

private

Definition at line 65 of file ECALpedestalPCLworker.h.

Referenced by bookHistograms(), and ECALpedestalPCLworker().

std::vector<MonitorElement *> ECALpedestalPCLworker::meEB_

private

Definition at line 54 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement *> ECALpedestalPCLworker::meEE_

private

Definition at line 55 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and bookHistograms().

int ECALpedestalPCLworker::nBins_

private

Definition at line 66 of file ECALpedestalPCLworker.h.

Referenced by bookHistograms(), and ECALpedestalPCLworker().

uint32_t ECALpedestalPCLworker::pedestalSamples_

private

Definition at line 57 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and ECALpedestalPCLworker().

bool ECALpedestalPCLworker::requireStableBeam_

private

Definition at line 69 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and ECALpedestalPCLworker().

uint32_t ECALpedestalPCLworker::sThresholdEB_

private

Definition at line 59 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and ECALpedestalPCLworker().

uint32_t ECALpedestalPCLworker::sThresholdEE_

private

Definition at line 62 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and ECALpedestalPCLworker().

edm::EDGetTokenT<TCDSRecord> ECALpedestalPCLworker::tcdsToken_

private

Definition at line 52 of file ECALpedestalPCLworker.h.

Referenced by analyze(), and ECALpedestalPCLworker().