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

	DQMEDHarvester (edm::ParameterSet const &iConfig)

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

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

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &es) 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::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

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

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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	checkAndRetrieveTopology (const edm::EventSetup &setup)

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	IsGoodLandauFit (double *FitResults)

bool	produceTagFilter (const MonitorElement *Charge_Vs_Index)

int	statCollectionFromMode (const char *tag) 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_

int	CalibrationLevel

bool	doChargeMonitorPerPlane

bool	doStoreOnDB

std::vector< std::string >	dqm_tag_

unsigned int	GOOD

std::string	m_calibrationMode

std::string	m_DQMdir

std::string	m_Record

unsigned int	MASKED

double	MinNrEntries

int	NPixelDets

int	NStripAPVs

double	tagCondition_GoodFrac

double	tagCondition_NClusters

edm::ESHandle< TrackerGeometry >	tkGeom_

const TrackerTopology *	tTopo_

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
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<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

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

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

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)

template<Transition B>
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 (std::string instanceName)

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

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Protected Attributes inherited from DQMEDHarvester
DQMStore *	dqmstore_

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 45 of file SiStripGainsPCLHarvester.h.

Constructor & Destructor Documentation

◆ SiStripGainsPCLHarvester()

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

explicit

Definition at line 37 of file SiStripGainsPCLHarvester.cc.

     : doStoreOnDB(false),
       GOOD(0),
       BAD(0),
       MASKED(0),
       NStripAPVs(0),
       NPixelDets(0),
       bareTkGeomPtr_(nullptr),
       tTopo_(nullptr) {
   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");
  
   //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
 }

References CalibrationLevel, doChargeMonitorPerPlane, dqm_tag_, edm::ParameterSet::getUntrackedParameter(), m_calibrationMode, m_DQMdir, m_Record, MinNrEntries, AlCaHLTBitMon_QueryRunRegistry::string, tagCondition_GoodFrac, tagCondition_NClusters, and VChargeHisto.

Member Function Documentation

◆ algoComputeMPVandGain()

void SiStripGainsPCLHarvester::algoComputeMPVandGain ( const MonitorElement * Charge_Vs_Index )

private

Definition at line 378 of file SiStripGainsPCLHarvester.cc.

                                                                                           {
   unsigned int I = 0;
   TH1F* Proj = nullptr;
   double FitResults[6];
   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 = (TH1F*)(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);
       TH1F* Proj2 = (TH1F*)(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);
         TH1F* Proj2 = (TH1F*)(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);
     }
  
     getPeakOfLandau(Proj, FitResults);
     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();
  
     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");
 }

References sistrip::APV, APVsColl, BAD, CalibrationLevel, getPeakOfLandau(), GOOD, Exhume::I, mps_fire::i, IsGoodLandauFit(), m_calibrationMode, MASKED, and mps_setup::stdout.

Referenced by dqmEndJob().

◆ beginRun()

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

override

Definition at line 70 of file SiStripGainsPCLHarvester.cc.

                                                                                       {
   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
   this->checkAndRetrieveTopology(iSetup);
  
   edm::ESHandle<SiStripGain> gainHandle;
   iSetup.get<SiStripGainRcd>().get(gainHandle);
   if (!gainHandle.isValid()) {
     edm::LogError("SiStripGainPCLHarvester") << "gainHandle is not valid\n";
     exit(0);
   }
  
   edm::ESHandle<SiStripQuality> SiStripQuality_;
   iSetup.get<SiStripQualityRcd>().get(SiStripQuality_);
  
   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 = SiStripQuality_->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;
   }
 }

References a, sistrip::APV, APVsCollOrdered, checkAndRetrieveTopology(), checkBookAPVColls(), beamvalidation::exit(), edm::EventSetup::get(), get, SiStripGain::getApvGain(), SiStripGain::getNumberOfTags(), SiStripGain::getRange(), SiStripQuality::IsApvBad(), edm::ESHandleBase::isValid(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, and mps_setup::stdout.

◆ checkAndRetrieveTopology()

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

privatevirtual

Definition at line 628 of file SiStripGainsPCLHarvester.cc.

                                                                                   {
   if (!tTopo_) {
     edm::ESHandle<TrackerTopology> TopoHandle;
     setup.get<TrackerTopologyRcd>().get(TopoHandle);
     tTopo_ = TopoHandle.product();
   }
 }

References get, edm::ESHandle< T >::product(), singleTopDQM_cfi::setup, and tTopo_.

Referenced by beginRun().

◆ checkBookAPVColls()

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

privatevirtual

Definition at line 518 of file SiStripGainsPCLHarvester.cc.

                                                                         {
   es.get<TrackerDigiGeometryRecord>().get(tkGeom_);
   const TrackerGeometry* newBareTkGeomPtr = &(*tkGeom_);
   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->APVId = j;
           APV->SubDet = SubDet;
           APV->FitMPV = -1;
           APV->FitMPVErr = -1;
           APV->FitWidth = -1;
           APV->FitWidthErr = -1;
           APV->FitChi2 = -1;
           APV->FitNorm = -1;
           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->APVId = (j << 3 | i);
             APV->SubDet = SubDet;
             APV->FitMPV = -1;
             APV->FitMPVErr = -1;
             APV->FitWidth = -1;
             APV->FitWidthErr = -1;
             APV->FitChi2 = -1;
             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;
 }

References sistrip::APV, APVsColl, APVsCollOrdered, bareTkGeomPtr_, TrackerGeometry::dets(), edm::EventSetup::get(), get, mps_fire::i, dqmiolumiharvest::j, PixelTopology::ncolumns(), NPixelDets, PixelTopology::nrows(), NStripAPVs, StripTopology::nstrips(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), ntupleEnum::SubDet, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, tkGeom_, and StripSubdetector::TOB.

Referenced by beginRun().

◆ dqmEndJob()

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

overrideprivatevirtual

Implements DQMEDHarvester.

Definition at line 120 of file SiStripGainsPCLHarvester.cc.

                                                                                              {
   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->writeOne(theAPVGains.get(), 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);
 }

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

◆ endRun()

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

override

Definition at line 719 of file SiStripGainsPCLHarvester.cc.

719 {}

◆ fillDescriptions()

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

static

Definition at line 712 of file SiStripGainsPCLHarvester.cc.

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

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

◆ gainQualityMonitor()

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

private

Definition at line 161 of file SiStripGainsPCLHarvester.cc.

                                                                                                {
   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 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(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_);
       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);
   }
 }

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(), ntupleEnum::SubDet, APVGain::subdetectorId(), APVGain::subdetectorPlane(), APVGain::subdetectorSide(), tTopo_, VChargeHisto, and z.

Referenced by dqmEndJob().

◆ getNewObject()

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

private

Definition at line 666 of file SiStripGainsPCLHarvester.cc.

                                                                                                           {
   std::unique_ptr<SiStripApvGain> obj = std::unique_ptr<SiStripApvGain>(new 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;
 }

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

Referenced by dqmEndJob().

◆ getPeakOfLandau()

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

private

Definition at line 482 of file SiStripGainsPCLHarvester.cc.

                                                                                                                      {
   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;
  
   // 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);
 }

References MinNrEntries, and tools::TF1.

Referenced by algoComputeMPVandGain().

◆ IsGoodLandauFit()

bool SiStripGainsPCLHarvester::IsGoodLandauFit ( double * FitResults )

private

Definition at line 508 of file SiStripGainsPCLHarvester.cc.

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

Referenced by algoComputeMPVandGain().

◆ produceTagFilter()

bool SiStripGainsPCLHarvester::produceTagFilter ( const MonitorElement * Charge_Vs_Index )

private

Definition at line 637 of file SiStripGainsPCLHarvester.cc.

                                                                                      {
   // 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;
 }

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

Referenced by getNewObject().

◆ statCollectionFromMode()

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

private

Member Data Documentation

◆ APVsColl

std::unordered_map<unsigned int, std::shared_ptr<stAPVGain> > SiStripGainsPCLHarvester::APVsColl

private

Definition at line 96 of file SiStripGainsPCLHarvester.h.

Referenced by algoComputeMPVandGain(), and checkBookAPVColls().

◆ APVsCollOrdered

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

private

Definition at line 95 of file SiStripGainsPCLHarvester.h.

Referenced by beginRun(), checkBookAPVColls(), gainQualityMonitor(), and getNewObject().

◆ BAD

unsigned int SiStripGainsPCLHarvester::BAD

private

Definition at line 71 of file SiStripGainsPCLHarvester.h.

Referenced by algoComputeMPVandGain(), and produceTagFilter().

◆ bareTkGeomPtr_

const TrackerGeometry* SiStripGainsPCLHarvester::bareTkGeomPtr_

private

Definition at line 92 of file SiStripGainsPCLHarvester.h.

Referenced by checkBookAPVColls().

◆ CalibrationLevel

int SiStripGainsPCLHarvester::CalibrationLevel

private

Definition at line 89 of file SiStripGainsPCLHarvester.h.

Referenced by algoComputeMPVandGain(), and SiStripGainsPCLHarvester().

◆ doChargeMonitorPerPlane

bool SiStripGainsPCLHarvester::doChargeMonitorPerPlane

private

Charge monitor per detector plane

Definition at line 69 of file SiStripGainsPCLHarvester.h.

Referenced by gainQualityMonitor(), and SiStripGainsPCLHarvester().

◆ doStoreOnDB

bool SiStripGainsPCLHarvester::doStoreOnDB

private

Definition at line 68 of file SiStripGainsPCLHarvester.h.

Referenced by dqmEndJob(), and getNewObject().

◆ dqm_tag_

std::vector<std::string> SiStripGainsPCLHarvester::dqm_tag_

private

Definition at line 87 of file SiStripGainsPCLHarvester.h.

Referenced by dqmEndJob(), and SiStripGainsPCLHarvester().

◆ GOOD

unsigned int SiStripGainsPCLHarvester::GOOD

private

Definition at line 70 of file SiStripGainsPCLHarvester.h.

Referenced by algoComputeMPVandGain(), and produceTagFilter().

◆ m_calibrationMode

std::string SiStripGainsPCLHarvester::m_calibrationMode

private

Type of statistics for the calibration

Definition at line 84 of file SiStripGainsPCLHarvester.h.

Referenced by algoComputeMPVandGain(), dqmEndJob(), produceTagFilter(), and SiStripGainsPCLHarvester().

◆ m_DQMdir

std::string SiStripGainsPCLHarvester::m_DQMdir

private

DQM folder hosting the charge statistics and the monitor plots

Definition at line 83 of file SiStripGainsPCLHarvester.h.

Referenced by dqmEndJob(), and SiStripGainsPCLHarvester().

◆ m_Record

std::string SiStripGainsPCLHarvester::m_Record

private

Definition at line 81 of file SiStripGainsPCLHarvester.h.

Referenced by dqmEndJob(), and SiStripGainsPCLHarvester().

◆ MASKED

unsigned int SiStripGainsPCLHarvester::MASKED

private

Definition at line 72 of file SiStripGainsPCLHarvester.h.

Referenced by algoComputeMPVandGain().

◆ MinNrEntries

double SiStripGainsPCLHarvester::MinNrEntries

private

Definition at line 79 of file SiStripGainsPCLHarvester.h.

Referenced by getPeakOfLandau(), and SiStripGainsPCLHarvester().

◆ NPixelDets

int SiStripGainsPCLHarvester::NPixelDets

private

Definition at line 78 of file SiStripGainsPCLHarvester.h.

Referenced by checkBookAPVColls().

◆ NStripAPVs

int SiStripGainsPCLHarvester::NStripAPVs

private

Definition at line 77 of file SiStripGainsPCLHarvester.h.

Referenced by checkBookAPVColls(), and produceTagFilter().

◆ tagCondition_GoodFrac

double SiStripGainsPCLHarvester::tagCondition_GoodFrac

private

Definition at line 75 of file SiStripGainsPCLHarvester.h.

Referenced by produceTagFilter(), and SiStripGainsPCLHarvester().

◆ tagCondition_NClusters

double SiStripGainsPCLHarvester::tagCondition_NClusters

private

Definition at line 74 of file SiStripGainsPCLHarvester.h.

Referenced by produceTagFilter(), and SiStripGainsPCLHarvester().

◆ tkGeom_

edm::ESHandle<TrackerGeometry> SiStripGainsPCLHarvester::tkGeom_

private

Definition at line 91 of file SiStripGainsPCLHarvester.h.

Referenced by checkBookAPVColls().

◆ tTopo_

const TrackerTopology* SiStripGainsPCLHarvester::tTopo_

private

Definition at line 93 of file SiStripGainsPCLHarvester.h.

Referenced by checkAndRetrieveTopology(), and gainQualityMonitor().

◆ VChargeHisto

std::vector<std::string> SiStripGainsPCLHarvester::VChargeHisto

private

Charge monitor plots to be output

Definition at line 85 of file SiStripGainsPCLHarvester.h.

Referenced by gainQualityMonitor(), and SiStripGainsPCLHarvester().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SiStripGainsPCLHarvester()

Member Function Documentation

◆ algoComputeMPVandGain()

◆ beginRun()

◆ checkAndRetrieveTopology()

◆ checkBookAPVColls()

◆ dqmEndJob()

◆ endRun()

◆ fillDescriptions()

◆ gainQualityMonitor()

◆ getNewObject()

◆ getPeakOfLandau()

◆ IsGoodLandauFit()

◆ produceTagFilter()

◆ statCollectionFromMode()

Member Data Documentation

◆ APVsColl

◆ APVsCollOrdered

◆ BAD

◆ bareTkGeomPtr_

◆ CalibrationLevel

◆ doChargeMonitorPerPlane

◆ doStoreOnDB

◆ dqm_tag_

◆ GOOD

◆ m_calibrationMode

◆ m_DQMdir

◆ m_Record

◆ MASKED

◆ MinNrEntries

◆ NPixelDets

◆ NStripAPVs

◆ tagCondition_GoodFrac

◆ tagCondition_NClusters

◆ tkGeom_

◆ tTopo_

◆ VChargeHisto