#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	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

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

	DQMEDHarvester ()

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

void	endJob () final

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

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

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

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

	~DQMEDHarvester () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

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

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)

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

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

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 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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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

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

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

Protected Attributes inherited from DQMEDHarvester
edm::EDPutTokenT< DQMToken >	lumiToken_

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 ( const edm::ParameterSet & ps )

explicit

Definition at line 37 of file SiStripGainsPCLHarvester.cc.

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

                                                                            :
   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
 
 }

Member Function Documentation

void SiStripGainsPCLHarvester::algoComputeMPVandGain ( const MonitorElement * Charge_Vs_Index )

private

Definition at line 338 of file SiStripGainsPCLHarvester.cc.

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

Referenced by dqmEndJob().

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

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

override

Definition at line 74 of file SiStripGainsPCLHarvester.cc.

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

 {
   using namespace edm;
 
   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 = 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;      
   }
 }

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

privatevirtual

Definition at line 566 of file SiStripGainsPCLHarvester.cc.

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

Referenced by beginRun().

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

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

privatevirtual

Definition at line 456 of file SiStripGainsPCLHarvester.cc.

References sistrip::APV, APVsColl, APVsCollOrdered, bareTkGeomPtr_, TrackerGeometry::dets(), edm::EventSetup::get(), mps_fire::i, 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().

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

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

overrideprivatevirtual

Implements DQMEDHarvester.

Definition at line 112 of file SiStripGainsPCLHarvester.cc.

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

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

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

override

Definition at line 661 of file SiStripGainsPCLHarvester.cc.

661 {

662 }

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

static

Definition at line 652 of file SiStripGainsPCLHarvester.cc.

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

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

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

private

Definition at line 156 of file SiStripGainsPCLHarvester.cc.

References a, sistrip::APV, APVsCollOrdered, stringResolutionProvider_cfi::bin, DQMStore::IBooker::book1DD(), DQMStore::IBooker::book2DD(), DetId::DetId(), doChargeMonitorPerPlane, MillePedeFileConverter_cfg::e, reco::tau::disc::Eta(), APVGain::FetchMonitor(), MonitorElement::Fill(), plotBeamSpotDB::first, g4SimHits_cfi::Gain, h, mps_fire::i, APVGain::monHnames(), colinearityKinematic::Phi, dttmaxenums::R, edm::second(), DQMStore::IBooker::setCurrentFolder(), ntupleEnum::SubDet, APVGain::subdetectorId(), APVGain::subdetectorPlane(), APVGain::subdetectorSide(), 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 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);
 
   }
 }

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

private

Definition at line 607 of file SiStripGainsPCLHarvester.cc.

References a, sistrip::APV, APVsCollOrdered, doStoreOnDB, LogDebug, GetRecoTauVFromDQM_MC_cff::obj, and produceTagFilter().

Referenced by dqmEndJob().

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

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

private

Definition at line 417 of file SiStripGainsPCLHarvester.cc.

References MinNrEntries.

Referenced by algoComputeMPVandGain().

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

bool SiStripGainsPCLHarvester::IsGoodLandauFit ( double * FitResults )

private

Definition at line 445 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;   
 }

bool SiStripGainsPCLHarvester::produceTagFilter ( const MonitorElement * Charge_Vs_Index )

private

Definition at line 575 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; 
 }