#include <SiStripGainsPCLHarvester.h>

Inheritance diagram for SiStripGainsPCLHarvester:

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

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

	SiStripGainsPCLHarvester (const edm::ParameterSet &ps)

Public Member Functions inherited from DQMEDHarvester
virtual void	analyze (edm::Event const &, edm::EventSetup const &) final

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

	DQMEDHarvester (void)

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

virtual void	endJob () final

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

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	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 &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

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

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
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)

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

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	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_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

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)

Detailed Description

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

Definition at line 44 of file SiStripGainsPCLHarvester.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file SiStripGainsPCLHarvester.cc.

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

                                                                            :
   doStoreOnDB(false),
   GOOD(0),
   BAD(0),
   MASKED(0),
   NStripAPVs(0),
   NPixelDets(0),
   bareTkGeomPtr_(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);
   
   //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 146 of file SiStripGainsPCLHarvester.cc.

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

Referenced by dqmEndJob().

                                                                                      {
 
   unsigned int I=0;
   TH1F* Proj = NULL;
   double FitResults[6];
   double MPVmean = 300;
 
   if ( Charge_Vs_Index==0 ) {
     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
	)

virtual

Reimplemented from DQMEDHarvester.

Definition at line 70 of file SiStripGainsPCLHarvester.cc.

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

 {
   using namespace edm;
 
   this->checkBookAPVColls(iSetup); // check whether APV colls are booked and do so if not yet done
 
   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::checkBookAPVColls ( const edm::EventSetup & setup )

privatevirtual

Definition at line 264 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 || 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_
	)

privatevirtual

Implements DQMEDHarvester.

Definition at line 107 of file SiStripGainsPCLHarvester.cc.

References algoComputeMPVandGain(), cond::service::PoolDBOutputService::currentTime(), doStoreOnDB, dqm_tag_, spr::find(), 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.size()!=0 && stag[0]!='_') stag.insert(0,1,'_');
 
   std::string cvi      = DQM_dir + std::string("/Charge_Vs_Index") + stag;
      
   MonitorElement* Charge_Vs_Index           = igetter_.get(cvi.c_str());
   
   if (Charge_Vs_Index==0) {
     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;  
   }
 }

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

virtual

Reimplemented from DQMEDHarvester.

Definition at line 460 of file SiStripGainsPCLHarvester.cc.

460 {

461 }

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

static

Definition at line 451 of file SiStripGainsPCLHarvester.cc.

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

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

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

private

Definition at line 406 of file SiStripGainsPCLHarvester.cc.

References a, sistrip::APV, APVsCollOrdered, doStoreOnDB, LogDebug, NULL, MuonAssociatorByHits_cfi::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==NULL){ 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 225 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 253 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 374 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==0 ) {
     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; 
 }