#include <SiStripGainsPCLHarvester.h>

Inheritance diagram for SiStripGainsPCLHarvester:

Public Member Functions
void	beginRun (edm::Run const &run, edm::EventSetup const &isetup) override

void	endRun (edm::Run const &run, edm::EventSetup const &isetup) override

	SiStripGainsPCLHarvester (const edm::ParameterSet &ps)

Public Member Functions inherited from DQMEDHarvester
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginJob () override

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &) final

void	beginRun (edm::Run const &run, edm::EventSetup const &) override

virtual void	dqmAnalyze (DQMStore::IBooker &, DQMStore::IGetter &, edm::Event const &, edm::EventSetup const &)

	DQMEDHarvester (edm::ParameterSet const &iConfig)

	DQMEDHarvester ()

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

virtual void	dqmEndRun (DQMStore::IBooker &, DQMStore::IGetter &, edm::Run const &, edm::EventSetup const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &es) final

void	endProcessBlockProduce (edm::ProcessBlock &) final

void	endRun (edm::Run const &, edm::EventSetup const &) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &es) final

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::EndProcessBlockProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDProducerBase () override

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 ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	algoComputeMPVandGain (const MonitorElement *Charge_Vs_Index)

virtual void	checkBookAPVColls (const edm::EventSetup &setup)

void	dqmEndJob (DQMStore::IBooker &ibooker_, DQMStore::IGetter &igetter_) override

void	gainQualityMonitor (DQMStore::IBooker &ibooker_, const MonitorElement *Charge_Vs_Index) const

std::unique_ptr< SiStripApvGain >	getNewObject (const MonitorElement *Charge_Vs_Index)

void	getPeakOfLandau (TH1 InputHisto, double FitResults, double LowRange=50, double HighRange=5400, bool gaussianConvolution=false)

bool	IsGoodLandauFit (double *FitResults)

bool	produceTagFilter (const MonitorElement *Charge_Vs_Index)

int	statCollectionFromMode (const char *tag) const

void	storeGainsTree (const TAxis *chVsIdxXaxis) const

Private Attributes
std::unordered_map< unsigned int, std::shared_ptr< stAPVGain > >	APVsColl

std::vector< std::shared_ptr< stAPVGain > >	APVsCollOrdered

unsigned int	BAD

const TrackerGeometry *	bareTkGeomPtr_ = nullptr

int	CalibrationLevel

bool	doChargeMonitorPerPlane

bool	doStoreOnDB

std::vector< std::string >	dqm_tag_

bool	fit_dataDrivenRange_ = false

bool	fit_gaussianConvolution_ = false

bool	fit_gaussianConvolutionTOBL56_ = false

edm::ESGetToken< SiStripGain, SiStripGainRcd >	gainToken_

unsigned int	GOOD

std::string	m_calibrationMode

std::string	m_DQMdir

std::string	m_Record

unsigned int	MASKED

double	MinNrEntries

int	NPixelDets

int	NStripAPVs

edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	qualityToken_

bool	storeGainsTree_

double	tagCondition_GoodFrac

double	tagCondition_NClusters

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

std::unique_ptr< TrackerTopology >	tTopo_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoTokenBR_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoTokenER_

std::vector< std::string >	VChargeHisto

Additional Inherited Members
Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore	DQMStore

typedef dqm::harvesting::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Protected Attributes inherited from DQMEDHarvester
DQMStore *	dqmstore_

edm::GetterOfProducts< DQMToken >	jobmegetter_

edm::EDPutTokenT< DQMToken >	jobToken_

edm::GetterOfProducts< DQMToken >	lumimegetter_

edm::EDPutTokenT< DQMToken >	lumiToken_

edm::GetterOfProducts< DQMToken >	runmegetter_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Description: Harvests output of SiStripGainsPCLWorker and creates histograms and Gains Payload

Definition at line 50 of file SiStripGainsPCLHarvester.h.

Constructor & Destructor Documentation

◆ SiStripGainsPCLHarvester()

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

explicit

Definition at line 34 of file SiStripGainsPCLHarvester.cc.

References CalibrationLevel, doChargeMonitorPerPlane, dqm_tag_, fit_dataDrivenRange_, fit_gaussianConvolution_, fit_gaussianConvolutionTOBL56_, gainToken_, edm::ParameterSet::getUntrackedParameter(), m_calibrationMode, m_DQMdir, m_Record, MinNrEntries, qualityToken_, storeGainsTree_, AlCaHLTBitMon_QueryRunRegistry::string, tagCondition_GoodFrac, tagCondition_NClusters, tkGeomToken_, tTopoTokenBR_, tTopoTokenER_, and VChargeHisto.

     : doStoreOnDB(false), GOOD(0), BAD(0), MASKED(0), NStripAPVs(0), NPixelDets(0) {
   m_Record = ps.getUntrackedParameter<std::string>("Record", "SiStripApvGainRcd");
   CalibrationLevel = ps.getUntrackedParameter<int>("CalibrationLevel", 0);
   MinNrEntries = ps.getUntrackedParameter<double>("minNrEntries", 20);
   m_DQMdir = ps.getUntrackedParameter<std::string>("DQMdir", "AlCaReco/SiStripGains");
   m_calibrationMode = ps.getUntrackedParameter<std::string>("calibrationMode", "StdBunch");
   tagCondition_NClusters = ps.getUntrackedParameter<double>("NClustersForTagProd", 2E8);
   tagCondition_GoodFrac = ps.getUntrackedParameter<double>("GoodFracForTagProd", 0.95);
   doChargeMonitorPerPlane = ps.getUntrackedParameter<bool>("doChargeMonitorPerPlane", false);
   VChargeHisto = ps.getUntrackedParameter<std::vector<std::string>>("ChargeHisto");
   fit_gaussianConvolution_ = ps.getUntrackedParameter<bool>("FitGaussianConvolution", false);
   fit_gaussianConvolutionTOBL56_ = ps.getUntrackedParameter<bool>("FitGaussianConvolutionTOBL5L6", false);
   fit_dataDrivenRange_ = ps.getUntrackedParameter<bool>("FitDataDrivenRange", false);
   storeGainsTree_ = ps.getUntrackedParameter<bool>("StoreGainsTree", false);
 
   //Set the monitoring element tag and store
   dqm_tag_.reserve(7);
   dqm_tag_.clear();
   dqm_tag_.push_back("StdBunch");    // statistic collection from Standard Collision Bunch @ 3.8 T
   dqm_tag_.push_back("StdBunch0T");  // statistic collection from Standard Collision Bunch @ 0 T
   dqm_tag_.push_back("AagBunch");    // statistic collection from First Collision After Abort Gap @ 3.8 T
   dqm_tag_.push_back("AagBunch0T");  // statistic collection from First Collision After Abort Gap @ 0 T
   dqm_tag_.push_back("IsoMuon");     // statistic collection from Isolated Muon @ 3.8 T
   dqm_tag_.push_back("IsoMuon0T");   // statistic collection from Isolated Muon @ 0 T
   dqm_tag_.push_back("Harvest");     // statistic collection: Harvest
 
   tTopoTokenBR_ = esConsumes<edm::Transition::BeginRun>();
   tTopoTokenER_ = esConsumes<edm::Transition::EndRun>();
   tkGeomToken_ = esConsumes<edm::Transition::BeginRun>();
   gainToken_ = esConsumes<edm::Transition::BeginRun>();
   qualityToken_ = esConsumes<edm::Transition::BeginRun>();
 }

Member Function Documentation

◆ algoComputeMPVandGain()

void SiStripGainsPCLHarvester::algoComputeMPVandGain ( const MonitorElement * Charge_Vs_Index )

private

Definition at line 490 of file SiStripGainsPCLHarvester.cc.

References A, sistrip::APV, APVsColl, cms::cuda::assert(), B, BAD, CalibrationLevel, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), EcalPhiSymFlatTableProducers_cfi::fill, fit_dataDrivenRange_, fit_gaussianConvolution_, fit_gaussianConvolutionTOBL56_, getPeakOfLandau(), GOOD, Exhume::I, mps_fire::i, IsGoodLandauFit(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, m_calibrationMode, MASKED, mps_setup::stdout, StripSubdetector::TOB, and tTopo_.

Referenced by dqmEndJob().

                                                                                           {
   unsigned int I = 0;
   TH1D* Proj = nullptr;
   static constexpr double DEF_F = -9999.;
   double FitResults[6] = {DEF_F, DEF_F, DEF_F, DEF_F, DEF_F, DEF_F};
   double MPVmean = 300;
 
   if (Charge_Vs_Index == nullptr) {
     edm::LogError("SiStripGainsPCLHarvester")
         << "Harvesting: could not execute algoComputeMPVandGain method because " << m_calibrationMode
         << " statistics cannot be retrieved.\n"
         << "Please check if input contains " << m_calibrationMode << " data." << std::endl;
     return;
   }
 
   TH2S* chvsidx = (Charge_Vs_Index)->getTH2S();
 
   printf("Progressing Bar              :0%%       20%%       40%%       60%%       80%%       100%%\n");
   printf("Fitting Charge Distribution  :");
   int TreeStep = APVsColl.size() / 50;
 
   for (auto it = APVsColl.begin(); it != APVsColl.end(); it++, I++) {
     if (I % TreeStep == 0) {
       printf(".");
       fflush(stdout);
     }
     std::shared_ptr<stAPVGain> APV = it->second;
     if (APV->Bin < 0)
       APV->Bin = chvsidx->GetXaxis()->FindBin(APV->Index);
 
     if (APV->isMasked) {
       APV->Gain = APV->PreviousGain;
       MASKED++;
       continue;
     }
 
     Proj = chvsidx->ProjectionY(
         "", chvsidx->GetXaxis()->FindBin(APV->Index), chvsidx->GetXaxis()->FindBin(APV->Index), "e");
     if (!Proj)
       continue;
 
     if (CalibrationLevel == 0) {
     } else if (CalibrationLevel == 1) {
       int SecondAPVId = APV->APVId;
       if (SecondAPVId % 2 == 0) {
         SecondAPVId = SecondAPVId + 1;
       } else {
         SecondAPVId = SecondAPVId - 1;
       }
       std::shared_ptr<stAPVGain> APV2 = APVsColl[(APV->DetId << 4) | SecondAPVId];
       if (APV2->Bin < 0)
         APV2->Bin = chvsidx->GetXaxis()->FindBin(APV2->Index);
       TH1D* Proj2 = chvsidx->ProjectionY("", APV2->Bin, APV2->Bin, "e");
       if (Proj2) {
         Proj->Add(Proj2, 1);
         delete Proj2;
       }
     } else if (CalibrationLevel == 2) {
       for (unsigned int i = 0; i < 16; i++) {  //loop up to 6APV for Strip and up to 16 for Pixels
         auto tmpit = APVsColl.find((APV->DetId << 4) | i);
         if (tmpit == APVsColl.end())
           continue;
         std::shared_ptr<stAPVGain> APV2 = tmpit->second;
         if (APV2->DetId != APV->DetId || APV2->APVId == APV->APVId)
           continue;
         if (APV2->Bin < 0)
           APV2->Bin = chvsidx->GetXaxis()->FindBin(APV2->Index);
         TH1D* Proj2 = chvsidx->ProjectionY("", APV2->Bin, APV2->Bin, "e");
         if (Proj2) {
           Proj->Add(Proj2, 1);
           delete Proj2;
         }
       }
     } else {
       CalibrationLevel = 0;
       printf("Unknown Calibration Level, will assume %i\n", CalibrationLevel);
     }
 
     std::pair<double, double> fitRange{50., 5400.};
     if (fit_dataDrivenRange_) {
       fitRange = findFitRange(Proj);
     }
 
     const bool isTOBL5L6 =
         (DetId{APV->DetId}.subdetId() == StripSubdetector::TOB) && (tTopo_->tobLayer(APV->DetId) > 4);
     getPeakOfLandau(Proj,
                     FitResults,
                     fitRange.first,
                     fitRange.second,
                     (isTOBL5L6 ? fit_gaussianConvolutionTOBL56_ : fit_gaussianConvolution_));
 
     // throw if the fit results are not set
     assert(FitResults[0] != DEF_F);
 
     APV->FitMPV = FitResults[0];
     APV->FitMPVErr = FitResults[1];
     APV->FitWidth = FitResults[2];
     APV->FitWidthErr = FitResults[3];
     APV->FitChi2 = FitResults[4];
     APV->FitNorm = FitResults[5];
     APV->NEntries = Proj->GetEntries();
 
     // fall back to legacy fit in case of  very low chi2 probability
     if (APV->FitChi2 <= 0.1) {
       edm::LogInfo("SiStripGainsPCLHarvester")
           << "fit failed on this APV:" << APV->DetId << ":" << APV->APVId << " !" << std::endl;
 
       std::fill(FitResults, FitResults + 6, 0);
       fitRange = std::make_pair(50., 5400.);
 
       APV->FitGrade = fitgrade::B;
 
       getPeakOfLandau(Proj, FitResults, fitRange.first, fitRange.second, false);
 
       APV->FitMPV = FitResults[0];
       APV->FitMPVErr = FitResults[1];
       APV->FitWidth = FitResults[2];
       APV->FitWidthErr = FitResults[3];
       APV->FitChi2 = FitResults[4];
       APV->FitNorm = FitResults[5];
       APV->NEntries = Proj->GetEntries();
     } else {
       APV->FitGrade = fitgrade::A;
     }
 
     if (IsGoodLandauFit(FitResults)) {
       APV->Gain = APV->FitMPV / MPVmean;
       if (APV->SubDet > 2)
         GOOD++;
     } else {
       APV->Gain = APV->PreviousGain;
       if (APV->SubDet > 2)
         BAD++;
     }
     if (APV->Gain <= 0)
       APV->Gain = 1;
 
     delete Proj;
   }
   printf("\n");
 }

◆ beginRun()

void SiStripGainsPCLHarvester::beginRun	(	edm::Run const &	run,
		edm::EventSetup const &	isetup
	)

override

Definition at line 70 of file SiStripGainsPCLHarvester.cc.

References a, sistrip::APV, APVsCollOrdered, checkBookAPVColls(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), beamvalidation::exit(), gainToken_, edm::EventSetup::getData(), edm::EventSetup::getHandle(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, qualityToken_, and mps_setup::stdout.

                                                                                       {
   using namespace edm;
   static constexpr float defaultGainTick = 690. / 640.;
 
   this->checkBookAPVColls(iSetup);  // check whether APV colls are booked and do so if not yet done
 
   const auto gainHandle = iSetup.getHandle(gainToken_);
   if (!gainHandle.isValid()) {
     edm::LogError("SiStripGainPCLHarvester") << "gainHandle is not valid\n";
     exit(0);
   }
 
   const auto& stripQuality = iSetup.getData(qualityToken_);
 
   for (unsigned int a = 0; a < APVsCollOrdered.size(); a++) {
     std::shared_ptr<stAPVGain> APV = APVsCollOrdered[a];
 
     if (APV->SubDet == PixelSubdetector::PixelBarrel || APV->SubDet == PixelSubdetector::PixelEndcap)
       continue;
 
     APV->isMasked = stripQuality.IsApvBad(APV->DetId, APV->APVId);
 
     if (gainHandle->getNumberOfTags() != 2) {
       edm::LogError("SiStripGainPCLHarvester") << "NUMBER OF GAIN TAG IS EXPECTED TO BE 2\n";
       fflush(stdout);
       exit(0);
     };
     float newPreviousGain = gainHandle->getApvGain(APV->APVId, gainHandle->getRange(APV->DetId, 1), 1);
     if (APV->PreviousGain != 1 and newPreviousGain != APV->PreviousGain)
       edm::LogWarning("SiStripGainPCLHarvester") << "WARNING: ParticleGain in the global tag changed\n";
     APV->PreviousGain = newPreviousGain;
 
     float newPreviousGainTick =
         APV->isMasked ? defaultGainTick : gainHandle->getApvGain(APV->APVId, gainHandle->getRange(APV->DetId, 0), 0);
     if (APV->PreviousGainTick != 1 and newPreviousGainTick != APV->PreviousGainTick) {
       edm::LogWarning("SiStripGainPCLHarvester")
           << "WARNING: TickMarkGain in the global tag changed\n"
           << std::endl
           << " APV->SubDet: " << APV->SubDet << " APV->APVId:" << APV->APVId << std::endl
           << " APV->PreviousGainTick: " << APV->PreviousGainTick << " newPreviousGainTick: " << newPreviousGainTick
           << std::endl;
     }
     APV->PreviousGainTick = newPreviousGainTick;
   }
 }

◆ checkBookAPVColls()

void SiStripGainsPCLHarvester::checkBookAPVColls ( const edm::EventSetup & setup )

privatevirtual

Definition at line 688 of file SiStripGainsPCLHarvester.cc.

References sistrip::APV, APVsColl, APVsCollOrdered, bareTkGeomPtr_, edm::EventSetup::getData(), mps_fire::i, dqmiolumiharvest::j, PixelTopology::ncolumns(), NONE, NPixelDets, PixelTopology::nrows(), NStripAPVs, StripTopology::nstrips(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, tkGeomToken_, StripSubdetector::TOB, and tTopoTokenBR_.

Referenced by beginRun().

                                                                         {
   auto newBareTkGeomPtr = &es.getData(tkGeomToken_);
   auto bareTkTopoPtr = &es.getData(tTopoTokenBR_);
   if (newBareTkGeomPtr == bareTkGeomPtr_)
     return;  // already filled APVColls, nothing changed
 
   if (!bareTkGeomPtr_) {  // pointer not yet set: called the first time => fill the APVColls
     auto const& Det = newBareTkGeomPtr->dets();
 
     unsigned int Index = 0;
 
     for (unsigned int i = 0; i < Det.size(); i++) {
       DetId Detid = Det[i]->geographicalId();
       int SubDet = Detid.subdetId();
 
       if (SubDet == StripSubdetector::TIB || SubDet == StripSubdetector::TID || SubDet == StripSubdetector::TOB ||
           SubDet == StripSubdetector::TEC) {
         auto DetUnit = dynamic_cast<const StripGeomDetUnit*>(Det[i]);
         if (!DetUnit)
           continue;
 
         const StripTopology& Topo = DetUnit->specificTopology();
         unsigned int NAPV = Topo.nstrips() / 128;
 
         for (unsigned int j = 0; j < NAPV; j++) {
           auto APV = std::make_shared<stAPVGain>();
           APV->Index = Index;
           APV->Bin = -1;
           APV->DetId = Detid.rawId();
           APV->Side = 0;
 
           if (SubDet == StripSubdetector::TID) {
             APV->Side = bareTkTopoPtr->tidSide(Detid);
           } else if (SubDet == StripSubdetector::TEC) {
             APV->Side = bareTkTopoPtr->tecSide(Detid);
           }
 
           APV->APVId = j;
           APV->SubDet = SubDet;
           APV->FitMPV = -1;
           APV->FitMPVErr = -1;
           APV->FitWidth = -1;
           APV->FitWidthErr = -1;
           APV->FitChi2 = -1;
           APV->FitNorm = -1;
           APV->FitGrade = fitgrade::NONE;
           APV->Gain = -1;
           APV->PreviousGain = 1;
           APV->PreviousGainTick = 1;
           APV->x = DetUnit->position().basicVector().x();
           APV->y = DetUnit->position().basicVector().y();
           APV->z = DetUnit->position().basicVector().z();
           APV->Eta = DetUnit->position().basicVector().eta();
           APV->Phi = DetUnit->position().basicVector().phi();
           APV->R = DetUnit->position().basicVector().transverse();
           APV->Thickness = DetUnit->surface().bounds().thickness();
           APV->NEntries = 0;
           APV->isMasked = false;
 
           APVsCollOrdered.push_back(APV);
           APVsColl[(APV->DetId << 4) | APV->APVId] = APV;
           Index++;
           NStripAPVs++;
         }  // loop on APVs
       }  // if is Strips
     }  // loop on dets
 
     for (unsigned int i = 0; i < Det.size();
          i++) {  //Make two loop such that the Pixel information is added at the end --> make transition simpler
       DetId Detid = Det[i]->geographicalId();
       int SubDet = Detid.subdetId();
       if (SubDet == PixelSubdetector::PixelBarrel || SubDet == PixelSubdetector::PixelEndcap) {
         auto DetUnit = dynamic_cast<const PixelGeomDetUnit*>(Det[i]);
         if (!DetUnit)
           continue;
 
         const PixelTopology& Topo = DetUnit->specificTopology();
         unsigned int NROCRow = Topo.nrows() / (80.);
         unsigned int NROCCol = Topo.ncolumns() / (52.);
 
         for (unsigned int j = 0; j < NROCRow; j++) {
           for (unsigned int i = 0; i < NROCCol; i++) {
             auto APV = std::make_shared<stAPVGain>();
             APV->Index = Index;
             APV->Bin = -1;
             APV->DetId = Detid.rawId();
             APV->Side = 0;
             APV->APVId = (j << 3 | i);
             APV->SubDet = SubDet;
             APV->FitMPV = -1;
             APV->FitMPVErr = -1;
             APV->FitWidth = -1;
             APV->FitWidthErr = -1;
             APV->FitChi2 = -1;
             APV->FitGrade = fitgrade::NONE;
             APV->Gain = -1;
             APV->PreviousGain = 1;
             APV->PreviousGainTick = 1;
             APV->x = DetUnit->position().basicVector().x();
             APV->y = DetUnit->position().basicVector().y();
             APV->z = DetUnit->position().basicVector().z();
             APV->Eta = DetUnit->position().basicVector().eta();
             APV->Phi = DetUnit->position().basicVector().phi();
             APV->R = DetUnit->position().basicVector().transverse();
             APV->Thickness = DetUnit->surface().bounds().thickness();
             APV->isMasked = false;  //SiPixelQuality_->IsModuleBad(Detid.rawId());
             APV->NEntries = 0;
 
             APVsCollOrdered.push_back(APV);
             APVsColl[(APV->DetId << 4) | APV->APVId] = APV;
             Index++;
             NPixelDets++;
 
           }  // loop on ROC cols
         }  // loop on ROC rows
       }  // if Pixel
     }  // loop on Dets
   }  //if (!bareTkGeomPtr_) ...
   bareTkGeomPtr_ = newBareTkGeomPtr;
 }

◆ dqmEndJob()

void SiStripGainsPCLHarvester::dqmEndJob	(	DQMStore::IBooker &	ibooker_,
		DQMStore::IGetter &	igetter_
	)

overrideprivatevirtual

Implements DQMEDHarvester.

Definition at line 117 of file SiStripGainsPCLHarvester.cc.

References algoComputeMPVandGain(), cond::service::PoolDBOutputService::currentTime(), doStoreOnDB, dqm_tag_, spr::find(), gainQualityMonitor(), dqm::implementation::IGetter::get(), getNewObject(), dqm::legacy::MonitorElement::getTH2S(), edm::Service< T >::isAvailable(), m_calibrationMode, m_DQMdir, m_Record, storeGainsTree(), storeGainsTree_, AlCaHLTBitMon_QueryRunRegistry::string, and cond::service::PoolDBOutputService::writeOneIOV().

                                                                                              {
   edm::LogInfo("SiStripGainsPCLHarvester") << "Starting harvesting statistics" << std::endl;
 
   std::string DQM_dir = m_DQMdir;
 
   std::string stag = *(std::find(dqm_tag_.begin(), dqm_tag_.end(), m_calibrationMode));
   if (!stag.empty() && stag[0] != '_')
     stag.insert(0, 1, '_');
 
   std::string cvi = DQM_dir + std::string("/Charge_Vs_Index") + stag;
 
   MonitorElement* Charge_Vs_Index = igetter_.get(cvi);
 
   if (Charge_Vs_Index == nullptr) {
     edm::LogError("SiStripGainsPCLHarvester")
         << "Harvesting: could not retrieve " << cvi.c_str() << ", statistics will not be summed!" << std::endl;
   } else {
     edm::LogInfo("SiStripGainsPCLHarvester")
         << "Harvesting " << (Charge_Vs_Index)->getTH2S()->GetEntries() << " more clusters" << std::endl;
   }
 
   algoComputeMPVandGain(Charge_Vs_Index);
   std::unique_ptr<SiStripApvGain> theAPVGains = this->getNewObject(Charge_Vs_Index);
 
   // write out the APVGains record
   edm::Service<cond::service::PoolDBOutputService> poolDbService;
 
   if (doStoreOnDB) {
     if (poolDbService.isAvailable())
       poolDbService->writeOneIOV(*theAPVGains, poolDbService->currentTime(), m_Record);
     else
       throw std::runtime_error("PoolDBService required.");
   } else {
     edm::LogInfo("SiStripGainsPCLHarvester") << "Will not produce payload!" << std::endl;
   }
 
   //Collect the statistics for monitoring and validation
   gainQualityMonitor(ibooker_, Charge_Vs_Index);
 
   if (storeGainsTree_) {
     if (Charge_Vs_Index != nullptr) {
       storeGainsTree(Charge_Vs_Index->getTH2S()->GetXaxis());
     } else {
       edm::LogError("SiStripGainsPCLHarvester")
           << "Harvesting: could not retrieve " << cvi.c_str() << "Tree won't be stored" << std::endl;
     }
   }
 }

◆ endRun()

void SiStripGainsPCLHarvester::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	isetup
	)

override

Definition at line 957 of file SiStripGainsPCLHarvester.cc.

References edm::EventSetup::getData(), tTopo_, and tTopoTokenER_.

                                                                                     {
   if (!tTopo_) {
     tTopo_ = std::make_unique<TrackerTopology>(isetup.getData(tTopoTokenER_));
   }
 }

◆ fillDescriptions()

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

static

Definition at line 950 of file SiStripGainsPCLHarvester.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

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

◆ gainQualityMonitor()

void SiStripGainsPCLHarvester::gainQualityMonitor	(	DQMStore::IBooker &	ibooker_,
		const MonitorElement *	Charge_Vs_Index
	)		const

private

Definition at line 167 of file SiStripGainsPCLHarvester.cc.

References a, sistrip::APV, APVsCollOrdered, newFWLiteAna::bin, dqm::implementation::IBooker::book1DD(), dqm::implementation::IBooker::book2DD(), doChargeMonitorPerPlane, MillePedeFileConverter_cfg::e, APVGain::FetchMonitor(), dqm::impl::MonitorElement::Fill(), dqmdumpme::first, IntegrityClient_cfi::Gain, h, mps_fire::i, APVGain::monHnames(), beam_dqm_sourceclient-live_cfg::monitor, VtxSmearedParameters_cfi::Phi, dttmaxenums::R, edm::second(), dqm::implementation::NavigatorBase::setCurrentFolder(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), APVGain::subdetectorId(), APVGain::subdetectorPlane(), APVGain::subdetectorSide(), APVGain::thickness(), tTopo_, VChargeHisto, and z.

Referenced by dqmEndJob().

                                                                                                {
   ibooker_.setCurrentFolder("AlCaReco/SiStripGainsHarvesting/");
 
   std::vector<APVGain::APVmon> new_charge_histos;
   std::vector<std::pair<std::string, std::string>> cnames =
       APVGain::monHnames(VChargeHisto, doChargeMonitorPerPlane, "newG2");
   for (unsigned int i = 0; i < cnames.size(); i++) {
     MonitorElement* monitor = ibooker_.book1DD((cnames[i]).first, (cnames[i]).second.c_str(), 100, 0., 1000.);
     int thick = APVGain::thickness((cnames[i]).first);
     int id = APVGain::subdetectorId((cnames[i]).first);
     int side = APVGain::subdetectorSide((cnames[i]).first);
     int plane = APVGain::subdetectorPlane((cnames[i]).first);
     new_charge_histos.push_back(APVGain::APVmon(thick, id, side, plane, monitor));
   }
 
   int MPVbin = 300;
   float MPVmin = 0.;
   float MPVmax = 600.;
 
   MonitorElement* MPV_Vs_EtaTIB =
       ibooker_.book2DD("MPVvsEtaTIB", "MPV vs Eta TIB", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_EtaTID =
       ibooker_.book2DD("MPVvsEtaTID", "MPV vs Eta TID", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_EtaTOB =
       ibooker_.book2DD("MPVvsEtaTOB", "MPV vs Eta TOB", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_EtaTEC =
       ibooker_.book2DD("MPVvsEtaTEC", "MPV vs Eta TEC", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_EtaTECthin =
       ibooker_.book2DD("MPVvsEtaTEC1", "MPV vs Eta TEC-thin", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_EtaTECthick =
       ibooker_.book2DD("MPVvsEtaTEC2", "MPV vs Eta TEC-thick", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
 
   MonitorElement* MPV_Vs_PhiTIB =
       ibooker_.book2DD("MPVvsPhiTIB", "MPV vs Phi TIB", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_PhiTID =
       ibooker_.book2DD("MPVvsPhiTID", "MPV vs Phi TID", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_PhiTOB =
       ibooker_.book2DD("MPVvsPhiTOB", "MPV vs Phi TOB", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_PhiTEC =
       ibooker_.book2DD("MPVvsPhiTEC", "MPV vs Phi TEC", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_PhiTECthin =
       ibooker_.book2DD("MPVvsPhiTEC1", "MPV vs Phi TEC-thin ", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
   MonitorElement* MPV_Vs_PhiTECthick =
       ibooker_.book2DD("MPVvsPhiTEC2", "MPV vs Phi TEC-thick", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
 
   MonitorElement* NoMPVfit = ibooker_.book2DD("NoMPVfit", "Modules with bad Landau Fit", 350, -350, 350, 240, 0, 120);
   MonitorElement* NoMPVmasked = ibooker_.book2DD("NoMPVmasked", "Masked Modules", 350, -350, 350, 240, 0, 120);
 
   MonitorElement* Gains = ibooker_.book1DD("Gains", "Gains", 300, 0, 2);
   MonitorElement* MPVs = ibooker_.book1DD("MPVs", "MPVs", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVs320 = ibooker_.book1DD("MPV_320", "MPV 320 thickness", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVs500 = ibooker_.book1DD("MPV_500", "MPV 500 thickness", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTIB = ibooker_.book1DD("MPV_TIB", "MPV TIB", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTID = ibooker_.book1DD("MPV_TID", "MPV TID", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTIDP = ibooker_.book1DD("MPV_TIDP", "MPV TIDP", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTIDM = ibooker_.book1DD("MPV_TIDM", "MPV TIDM", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTOB = ibooker_.book1DD("MPV_TOB", "MPV TOB", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTEC = ibooker_.book1DD("MPV_TEC", "MPV TEC", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTECP = ibooker_.book1DD("MPV_TECP", "MPV TECP", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTECM = ibooker_.book1DD("MPV_TECM", "MPV TECM", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTECthin = ibooker_.book1DD("MPV_TEC1", "MPV TEC thin", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTECthick = ibooker_.book1DD("MPV_TEC2", "MPV TEC thick", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTECP1 = ibooker_.book1DD("MPV_TECP1", "MPV TECP thin ", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTECP2 = ibooker_.book1DD("MPV_TECP2", "MPV TECP thick", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTECM1 = ibooker_.book1DD("MPV_TECM1", "MPV TECM thin", MPVbin, MPVmin, MPVmax);
   MonitorElement* MPVsTECM2 = ibooker_.book1DD("MPV_TECM2", "MPV TECM thick", MPVbin, MPVmin, MPVmax);
 
   MonitorElement* MPVError = ibooker_.book1DD("MPVError", "MPV Error", 150, 0, 150);
   MonitorElement* MPVErrorVsMPV = ibooker_.book2DD("MPVErrorVsMPV", "MPV Error vs MPV", 300, 0, 600, 150, 0, 150);
   MonitorElement* MPVErrorVsEta = ibooker_.book2DD("MPVErrorVsEta", "MPV Error vs Eta", 50, -3.0, 3.0, 150, 0, 150);
   MonitorElement* MPVErrorVsPhi = ibooker_.book2DD("MPVErrorVsPhi", "MPV Error vs Phi", 50, -3.4, 3.4, 150, 0, 150);
   MonitorElement* MPVErrorVsN = ibooker_.book2DD("MPVErrorVsN", "MPV Error vs N", 500, 0, 1000, 150, 0, 150);
 
   MonitorElement* DiffWRTPrevGainTIB = ibooker_.book1DD("DiffWRTPrevGainTIB", "Diff w.r.t. PrevGain TIB", 250, 0, 2);
   MonitorElement* DiffWRTPrevGainTID = ibooker_.book1DD("DiffWRTPrevGainTID", "Diff w.r.t. PrevGain TID", 250, 0, 2);
   MonitorElement* DiffWRTPrevGainTOB = ibooker_.book1DD("DiffWRTPrevGainTOB", "Diff w.r.t. PrevGain TOB", 250, 0, 2);
   MonitorElement* DiffWRTPrevGainTEC = ibooker_.book1DD("DiffWRTPrevGainTEC", "Diff w.r.t. PrevGain TEC", 250, 0, 2);
 
   MonitorElement* GainVsPrevGainTIB =
       ibooker_.book2DD("GainVsPrevGainTIB", "Gain vs PrevGain TIB", 100, 0, 2, 100, 0, 2);
   MonitorElement* GainVsPrevGainTID =
       ibooker_.book2DD("GainVsPrevGainTID", "Gain vs PrevGain TID", 100, 0, 2, 100, 0, 2);
   MonitorElement* GainVsPrevGainTOB =
       ibooker_.book2DD("GainVsPrevGainTOB", "Gain vs PrevGain TOB", 100, 0, 2, 100, 0, 2);
   MonitorElement* GainVsPrevGainTEC =
       ibooker_.book2DD("GainVsPrevGainTEC", "Gain vs PrevGain TEC", 100, 0, 2, 100, 0, 2);
 
   for (unsigned int a = 0; a < APVsCollOrdered.size(); a++) {
     std::shared_ptr<stAPVGain> APV = APVsCollOrdered[a];
     if (APV == nullptr)
       continue;
 
     unsigned int Index = APV->Index;
     unsigned int SubDet = APV->SubDet;
     unsigned int DetId = APV->DetId;
     float z = APV->z;
     float Eta = APV->Eta;
     float R = APV->R;
     float Phi = APV->Phi;
     float Thickness = APV->Thickness;
     double FitMPV = APV->FitMPV;
     double FitMPVErr = APV->FitMPVErr;
     double Gain = APV->Gain;
     double NEntries = APV->NEntries;
     double PreviousGain = APV->PreviousGain;
 
     if (SubDet < 3)
       continue;  // avoid to loop over Pixel det id
 
     if (Gain != 1.) {
       std::vector<MonitorElement*> charge_histos = APVGain::FetchMonitor(new_charge_histos, DetId, tTopo_.get());
 
       if (!Charge_Vs_Index)
         continue;
       TH2S* chvsidx = (Charge_Vs_Index)->getTH2S();
       int bin = chvsidx->GetXaxis()->FindBin(Index);
       TH1D* Proj = chvsidx->ProjectionY("proj", bin, bin);
       for (int binId = 0; binId < Proj->GetXaxis()->GetNbins(); binId++) {
         double new_charge = Proj->GetXaxis()->GetBinCenter(binId) / Gain;
         if (Proj->GetBinContent(binId) != 0.) {
           for (unsigned int h = 0; h < charge_histos.size(); h++) {
             TH1D* chisto = (charge_histos[h])->getTH1D();
             for (int e = 0; e < Proj->GetBinContent(binId); e++)
               chisto->Fill(new_charge);
           }
         }
       }
     }
 
     if (FitMPV <= 0.) {  // No fit of MPV
       if (APV->isMasked)
         NoMPVmasked->Fill(z, R);
       else
         NoMPVfit->Fill(z, R);
 
     } else {  // Fit of MPV
       if (FitMPV > 0.)
         Gains->Fill(Gain);
 
       MPVs->Fill(FitMPV);
       if (Thickness < 0.04)
         MPVs320->Fill(FitMPV);
       if (Thickness > 0.04)
         MPVs500->Fill(FitMPV);
 
       MPVError->Fill(FitMPVErr);
       MPVErrorVsMPV->Fill(FitMPV, FitMPVErr);
       MPVErrorVsEta->Fill(Eta, FitMPVErr);
       MPVErrorVsPhi->Fill(Phi, FitMPVErr);
       MPVErrorVsN->Fill(NEntries, FitMPVErr);
 
       if (SubDet == 3) {
         MPV_Vs_EtaTIB->Fill(Eta, FitMPV);
         MPV_Vs_PhiTIB->Fill(Phi, FitMPV);
         MPVsTIB->Fill(FitMPV);
 
       } else if (SubDet == 4) {
         MPV_Vs_EtaTID->Fill(Eta, FitMPV);
         MPV_Vs_PhiTID->Fill(Phi, FitMPV);
         MPVsTID->Fill(FitMPV);
         if (Eta < 0.)
           MPVsTIDM->Fill(FitMPV);
         if (Eta > 0.)
           MPVsTIDP->Fill(FitMPV);
 
       } else if (SubDet == 5) {
         MPV_Vs_EtaTOB->Fill(Eta, FitMPV);
         MPV_Vs_PhiTOB->Fill(Phi, FitMPV);
         MPVsTOB->Fill(FitMPV);
 
       } else if (SubDet == 6) {
         MPV_Vs_EtaTEC->Fill(Eta, FitMPV);
         MPV_Vs_PhiTEC->Fill(Phi, FitMPV);
         MPVsTEC->Fill(FitMPV);
         if (Eta < 0.)
           MPVsTECM->Fill(FitMPV);
         if (Eta > 0.)
           MPVsTECP->Fill(FitMPV);
         if (Thickness < 0.04) {
           MPV_Vs_EtaTECthin->Fill(Eta, FitMPV);
           MPV_Vs_PhiTECthin->Fill(Phi, FitMPV);
           MPVsTECthin->Fill(FitMPV);
           if (Eta > 0.)
             MPVsTECP1->Fill(FitMPV);
           if (Eta < 0.)
             MPVsTECM1->Fill(FitMPV);
         }
         if (Thickness > 0.04) {
           MPV_Vs_EtaTECthick->Fill(Eta, FitMPV);
           MPV_Vs_PhiTECthick->Fill(Phi, FitMPV);
           MPVsTECthick->Fill(FitMPV);
           if (Eta > 0.)
             MPVsTECP2->Fill(FitMPV);
           if (Eta < 0.)
             MPVsTECM2->Fill(FitMPV);
         }
       }
     }
 
     if (SubDet == 3 && PreviousGain != 0.)
       DiffWRTPrevGainTIB->Fill(Gain / PreviousGain);
     else if (SubDet == 4 && PreviousGain != 0.)
       DiffWRTPrevGainTID->Fill(Gain / PreviousGain);
     else if (SubDet == 5 && PreviousGain != 0.)
       DiffWRTPrevGainTOB->Fill(Gain / PreviousGain);
     else if (SubDet == 6 && PreviousGain != 0.)
       DiffWRTPrevGainTEC->Fill(Gain / PreviousGain);
 
     if (SubDet == 3)
       GainVsPrevGainTIB->Fill(PreviousGain, Gain);
     else if (SubDet == 4)
       GainVsPrevGainTID->Fill(PreviousGain, Gain);
     else if (SubDet == 5)
       GainVsPrevGainTOB->Fill(PreviousGain, Gain);
     else if (SubDet == 6)
       GainVsPrevGainTEC->Fill(PreviousGain, Gain);
   }
 }

◆ getNewObject()

std::unique_ptr< SiStripApvGain > SiStripGainsPCLHarvester::getNewObject ( const MonitorElement * Charge_Vs_Index )

private

Definition at line 839 of file SiStripGainsPCLHarvester.cc.

References a, sistrip::APV, APVsCollOrdered, doStoreOnDB, LogDebug, getGTfromDQMFile::obj, produceTagFilter(), and isotrackApplyRegressor::range.

Referenced by dqmEndJob().

                                                                                                           {
   std::unique_ptr<SiStripApvGain> obj = std::make_unique<SiStripApvGain>();
 
   if (!produceTagFilter(Charge_Vs_Index)) {
     edm::LogWarning("SiStripGainsPCLHarvester")
         << "getNewObject -> will not produce a paylaod because produceTagFilter returned false " << std::endl;
     return obj;
   } else {
     doStoreOnDB = true;
   }
 
   std::vector<float> theSiStripVector;
   unsigned int PreviousDetId = 0;
   for (unsigned int a = 0; a < APVsCollOrdered.size(); a++) {
     std::shared_ptr<stAPVGain> APV = APVsCollOrdered[a];
     if (APV == nullptr) {
       printf("Bug\n");
       continue;
     }
     if (APV->SubDet <= 2)
       continue;
     if (APV->DetId != PreviousDetId) {
       if (!theSiStripVector.empty()) {
         SiStripApvGain::Range range(theSiStripVector.begin(), theSiStripVector.end());
         if (!obj->put(PreviousDetId, range))
           printf("Bug to put detId = %i\n", PreviousDetId);
       }
       theSiStripVector.clear();
       PreviousDetId = APV->DetId;
     }
     theSiStripVector.push_back(APV->Gain);
 
     LogDebug("SiStripGainsPCLHarvester") << " DetId: " << APV->DetId << " APV:   " << APV->APVId
                                          << " Gain:  " << APV->Gain << std::endl;
   }
   if (!theSiStripVector.empty()) {
     SiStripApvGain::Range range(theSiStripVector.begin(), theSiStripVector.end());
     if (!obj->put(PreviousDetId, range))
       printf("Bug to put detId = %i\n", PreviousDetId);
   }
 
   return obj;
 }

◆ getPeakOfLandau()

void SiStripGainsPCLHarvester::getPeakOfLandau	(	TH1 *	InputHisto,
		double *	FitResults,
		double	LowRange = `50`,
		double	HighRange = `5400`,
		bool	gaussianConvolution = `false`
	)

private

Definition at line 633 of file SiStripGainsPCLHarvester.cc.

References personalPlayback::fp, MinNrEntries, and pfDeepBoostedJetPreprocessParams_cfi::sv.

Referenced by algoComputeMPVandGain().

                                                                                                       {
   // undo defaults (checked for assertion)
   FitResults[0] = -0.5;  //MPV
   FitResults[1] = 0;     //MPV error
   FitResults[2] = -0.5;  //Width
   FitResults[3] = 0;     //Width error
   FitResults[4] = -0.5;  //Fit Chi2/NDF
   FitResults[5] = 0;     //Normalization
 
   if (InputHisto->GetEntries() < MinNrEntries)
     return;
 
   if (gaussianConvolution) {
     // perform fit with landau convoluted with a gaussian
     Double_t fr[2] = {LowRange, HighRange};
     Double_t sv[4] = {25., 300., InputHisto->Integral(), 40.};
     Double_t pllo[4] = {0.5, 100., 1.0, 0.4};
     Double_t plhi[4] = {100., 500., 1.e7, 100.};
     Double_t fp[4], fpe[4];
     Double_t chisqr;
     Int_t ndf;
     auto fitsnr = langaufit(dynamic_cast<TH1D*>(InputHisto), fr, sv, pllo, plhi, fp, fpe, &chisqr, &ndf);
     FitResults[0] = fitsnr->GetMaximumX();  //MPV
     FitResults[1] = fpe[1];                 //MPV error // FIXME add error propagation
     FitResults[2] = fp[0];                  //Width
     FitResults[3] = fpe[0];                 //Width error
     FitResults[4] = chisqr / ndf;           //Fit Chi2/NDF
     FitResults[5] = fp[2];
   } else {
     // perform fit with standard landau
     TF1 MyLandau("MyLandau", "landau", LowRange, HighRange);
     MyLandau.SetParameter(1, 300);
     InputHisto->Fit(&MyLandau, "0QR WW");
     // MPV is parameter 1 (0=constant, 1=MPV, 2=Sigma)
     FitResults[0] = MyLandau.GetParameter(1);                     //MPV
     FitResults[1] = MyLandau.GetParError(1);                      //MPV error
     FitResults[2] = MyLandau.GetParameter(2);                     //Width
     FitResults[3] = MyLandau.GetParError(2);                      //Width error
     FitResults[4] = MyLandau.GetChisquare() / MyLandau.GetNDF();  //Fit Chi2/NDF
     FitResults[5] = MyLandau.GetParameter(0);
   }
 }

◆ IsGoodLandauFit()

bool SiStripGainsPCLHarvester::IsGoodLandauFit ( double * FitResults )

private

Definition at line 678 of file SiStripGainsPCLHarvester.cc.

Referenced by algoComputeMPVandGain().

                                                                  {
   if (FitResults[0] <= 0)
     return false;
   //   if(FitResults[1] > MaxMPVError   )return false;
   //   if(FitResults[4] > MaxChi2OverNDF)return false;
   return true;
 }

◆ produceTagFilter()

bool SiStripGainsPCLHarvester::produceTagFilter ( const MonitorElement * Charge_Vs_Index )

private

Definition at line 810 of file SiStripGainsPCLHarvester.cc.

References BAD, GOOD, funct::integral(), m_calibrationMode, NStripAPVs, tagCondition_GoodFrac, and tagCondition_NClusters.

Referenced by getNewObject().

                                                                                      {
   // The goal of this function is to check wether or not there is enough statistics
   // to produce a meaningful tag for the DB
 
   if (Charge_Vs_Index == nullptr) {
     edm::LogError("SiStripGainsPCLHarvester")
         << "produceTagFilter -> Return false: could not retrieve the " << m_calibrationMode << " statistics.\n"
         << "Please check if input contains " << m_calibrationMode << " data." << std::endl;
     return false;
   }
 
   float integral = (Charge_Vs_Index)->getTH2S()->Integral();
   if ((Charge_Vs_Index)->getTH2S()->Integral(0, NStripAPVs + 1, 0, 99999) < tagCondition_NClusters) {
     edm::LogWarning("SiStripGainsPCLHarvester")
         << "calibrationMode  -> " << m_calibrationMode << "\n"
         << "produceTagFilter -> Return false: Statistics is too low : " << integral << std::endl;
     return false;
   }
   if ((1.0 * GOOD) / (GOOD + BAD) < tagCondition_GoodFrac) {
     edm::LogWarning("SiStripGainsPCLHarvester")
         << "calibrationMode  -> " << m_calibrationMode << "\n"
         << "produceTagFilter ->  Return false: ratio of GOOD/TOTAL is too low: " << (1.0 * GOOD) / (GOOD + BAD)
         << std::endl;
     return false;
   }
   return true;
 }

◆ statCollectionFromMode()

int SiStripGainsPCLHarvester::statCollectionFromMode ( const char * tag ) const

private

◆ storeGainsTree()

void SiStripGainsPCLHarvester::storeGainsTree ( const TAxis * chVsIdxXaxis ) const

private

Definition at line 883 of file SiStripGainsPCLHarvester.cc.

References APVsCollOrdered, AlignmentPI::t_x, AlignmentPI::t_y, and AlignmentPI::t_z.

Referenced by dqmEndJob().

                                                                              {
   unsigned int t_Index, t_Bin, t_DetId;
   unsigned char t_APVId, t_SubDet;
   float t_x, t_y, t_z, t_Eta, t_R, t_Phi, t_Thickness;
   float t_FitMPV, t_FitMPVErr, t_FitWidth, t_FitWidthErr, t_FitChi2NDF, t_FitNorm, t_FitGrade;
   double t_Gain, t_PrevGain, t_PrevGainTick, t_NEntries;
   bool t_isMasked;
   auto tree = edm::Service<TFileService>()->make<TTree>("APVGain", "APVGain");
   tree->Branch("Index", &t_Index, "Index/i");
   tree->Branch("Bin", &t_Bin, "Bin/i");
   tree->Branch("DetId", &t_DetId, "DetId/i");
   tree->Branch("APVId", &t_APVId, "APVId/b");
   tree->Branch("SubDet", &t_SubDet, "SubDet/b");
   tree->Branch("x", &t_x, "x/F");
   tree->Branch("y", &t_y, "y/F");
   tree->Branch("z", &t_z, "z/F");
   tree->Branch("Eta", &t_Eta, "Eta/F");
   tree->Branch("R", &t_R, "R/F");
   tree->Branch("Phi", &t_Phi, "Phi/F");
   tree->Branch("Thickness", &t_Thickness, "Thickness/F");
   tree->Branch("FitMPV", &t_FitMPV, "FitMPV/F");
   tree->Branch("FitMPVErr", &t_FitMPVErr, "FitMPVErr/F");
   tree->Branch("FitWidth", &t_FitWidth, "FitWidth/F");
   tree->Branch("FitWidthErr", &t_FitWidthErr, "FitWidthErr/F");
   tree->Branch("FitChi2NDF", &t_FitChi2NDF, "FitChi2NDF/F");
   tree->Branch("FitNorm", &t_FitNorm, "FitNorm/F");
   tree->Branch("FitGrade", &t_FitGrade, "FitGrade/F");
   tree->Branch("Gain", &t_Gain, "Gain/D");
   tree->Branch("PrevGain", &t_PrevGain, "PrevGain/D");
   tree->Branch("PrevGainTick", &t_PrevGainTick, "PrevGainTick/D");
   tree->Branch("NEntries", &t_NEntries, "NEntries/D");
   tree->Branch("isMasked", &t_isMasked, "isMasked/O");
 
   for (const auto& iAPV : APVsCollOrdered) {
     if (iAPV) {
       t_Index = iAPV->Index;
       t_Bin = chVsIdxXaxis->FindBin(iAPV->Index);
       t_DetId = iAPV->DetId;
       t_APVId = iAPV->APVId;
       t_SubDet = iAPV->SubDet;
       t_x = iAPV->x;
       t_y = iAPV->y;
       t_z = iAPV->z;
       t_Eta = iAPV->Eta;
       t_R = iAPV->R;
       t_Phi = iAPV->Phi;
       t_Thickness = iAPV->Thickness;
       t_FitMPV = iAPV->FitMPV;
       t_FitMPVErr = iAPV->FitMPVErr;
       t_FitWidth = iAPV->FitWidth;
       t_FitWidthErr = iAPV->FitWidthErr;
       t_FitChi2NDF = iAPV->FitChi2;
       t_FitNorm = iAPV->FitNorm;
       t_FitGrade = iAPV->FitGrade;
       t_Gain = iAPV->Gain;
       t_PrevGain = iAPV->PreviousGain;
       t_PrevGainTick = iAPV->PreviousGainTick;
       t_NEntries = iAPV->NEntries;
       t_isMasked = iAPV->isMasked;
 
       tree->Fill();
     }
   }
 }