SiStripCorrelateNoise Class Reference

#include <DQM/SiStripMonitorSummary/plugins/SiStripCorrelateNoise.cc>

Inheritance diagram for SiStripCorrelateNoise:

Classes
struct	Data

Public Member Functions
	SiStripCorrelateNoise (const edm::ParameterSet &)
	~SiStripCorrelateNoise ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
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
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginRun (const edm::Run &run, const edm::EventSetup &es)
void	checkGainCache (const edm::EventSetup &es)
void	DoAnalysis (const edm::EventSetup &, const SiStripNoises &, SiStripNoises &)
void	DoPlots ()
virtual void	endJob ()
unsigned long long	getGainCache (const edm::EventSetup &eSetup)
float	getGainRatio (const uint32_t &detid, const uint16_t &apv)
TH1F *	getHisto (const long unsigned int &index)
void	getHistos (const uint32_t &detid, const TrackerTopology *tTopo, std::vector< TH1F * > &histos)
float	getMeanNoise (const SiStripNoises::Range &noiseRange, const uint32_t &first, const uint32_t &range)
unsigned long long	getNoiseCache (const edm::EventSetup &eSetup)

Private Attributes
unsigned long long	cacheID_gain
unsigned long long	cacheID_noise
bool	equalGain
TFile *	file
SiStripDetInfoFileReader *	fr
edm::ESHandle< SiStripApvGain >	gainHandle_
SiStripApvGain *	newGain
edm::ESHandle< SiStripNoises >	noiseHandle_
SiStripApvGain *	oldGain
SiStripNoises *	refNoise
uint32_t	theRun
TrackerMap *	tkmap
std::vector< TH1F * >	vTH1

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 50 of file SiStripCorrelateNoise.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file SiStripCorrelateNoise.cc.

References file, fr, Utilities::operator, and tkmap.

                                                                          :
  refNoise(0),oldGain(0),newGain(0),cacheID_noise(0xFFFFFFFF),cacheID_gain(0xFFFFFFFF)
{
  //now do what ever initialization is needed
  if(!edm::Service<SiStripDetInfoFileReader>().isAvailable()){
    edm::LogError("TkLayerMap") << 
      "\n------------------------------------------"
      "\nUnAvailable Service SiStripDetInfoFileReader: please insert in the configuration file an instance like"
      "\n\tprocess.SiStripDetInfoFileReader = cms.Service(\"SiStripDetInfoFileReader\")"
      "\n------------------------------------------";
  }
 
  fr=edm::Service<SiStripDetInfoFileReader>().operator->();
  file = new TFile("correlTest.root","RECREATE");

  file->cd();
  tkmap = new TrackerMap();
}

SiStripCorrelateNoise::~SiStripCorrelateNoise

(

)

Definition at line 29 of file SiStripCorrelateNoise.cc.

{}

Member Function Documentation

virtual void SiStripCorrelateNoise::analyze ( const edm::Event &  , const edm::EventSetup &    )

inlineprivatevirtual

Definition at line 58 of file SiStripCorrelateNoise.h.

References DoAnalysis(), DoPlots(), endJob(), getHisto(), getHistos(), plotFactory::histos, and diffTreeTool::index.

{};

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

privatevirtual

Definition at line 35 of file SiStripCorrelateNoise.cc.

References cacheID_noise, checkGainCache(), dir, DoAnalysis(), DoPlots(), file, edm::EventSetup::get(), getNoiseCache(), noiseHandle_, edm::ESHandle< T >::product(), refNoise, edm::RunBase::run(), and theRun.

                                                                         {

  if(getNoiseCache(es)==cacheID_noise )
    return;

  edm::LogInfo("") << "[SiStripCorrelateNoise::beginRun]  cacheID_noise " << cacheID_noise << std::endl; 
  
  es.get<SiStripNoisesRcd>().get(noiseHandle_);
  SiStripNoises * aNoise= new SiStripNoises(*noiseHandle_.product());
  
  //Check if gain is the same from one noise iov to the other, otherwise cache the new gain (and the old one) to rescale

  checkGainCache(es);

  if(cacheID_noise!=0xFFFFFFFF){
    char dir[128];
    theRun=run.run();
    sprintf(dir,"Run_%d",theRun);
    file->cd("");
    file->mkdir(dir);
    file->cd(dir);
    DoAnalysis(es,*noiseHandle_.product(),*refNoise);
    DoPlots(); 
  }

  cacheID_noise=getNoiseCache(es);  
  if(refNoise!=0)
    delete refNoise;
  refNoise=aNoise;
}

void SiStripCorrelateNoise::checkGainCache ( const edm::EventSetup &  es )

private

Definition at line 67 of file SiStripCorrelateNoise.cc.

References cacheID_gain, equalGain, gainHandle_, edm::EventSetup::get(), getGainCache(), newGain, oldGain, and edm::ESHandle< T >::product().

Referenced by beginRun(), and getGainCache().

                                                            {
  equalGain=true;
  if(getGainCache(es)!=cacheID_gain ){
    es.get<SiStripApvGainRcd>().get(gainHandle_);
    if(oldGain!=0)
      delete oldGain;
    
    oldGain = newGain;
    newGain = new SiStripApvGain(*gainHandle_.product());

    if(cacheID_gain!=0xFFFFFFFF)
      equalGain=false;
   cacheID_gain=getGainCache(es);
   edm::LogInfo("") << "[SiStripCorrelateNoise::checkGainCache]  cacheID_gain " << cacheID_gain << std::endl; 
  }
}

void SiStripCorrelateNoise::DoAnalysis ( const edm::EventSetup &  es, const SiStripNoises &  _Noise, SiStripNoises &  refNoise  )

private

Definition at line 116 of file SiStripCorrelateNoise.cc.

References mps_alisetup::collection, equalGain, TrackerMap::fill(), HcalObjRepresent::Fill(), edm::EventSetup::get(), getGainRatio(), getHistos(), plotFactory::histos, mps_fire::i, digi_MixPreMix_cfi::Noise, edm::ESHandle< T >::product(), refNoise, digi_MixPreMix_cfi::strip, tkmap, and relativeConstraints::value.

Referenced by analyze(), and beginRun().

                                                                                                              {
  SiStripNoises Noise = _Noise;
  typedef std::vector<SiStripNoises::ratioData> collection; 
  collection divNoise=Noise/refNoise;

  edm::LogInfo("") << "[Doanalysis]";

  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  es.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology* const tTopo = tTopoHandle.product();

  std::vector<TH1F *>histos;

  collection::const_iterator iter=divNoise.begin();
  collection::const_iterator iterE=divNoise.end();

  float value;
  float gainRatio=1.;
  //Divide result by d
  for(;iter!=iterE;++iter){
    getHistos(iter->detid, tTopo, histos);
    
    size_t strip=0, stripE= iter->values.size();
    size_t apvNb=7;

    for (;strip<stripE;++strip){       
      if(!equalGain && strip/128!=apvNb){
    apvNb=strip/128;
    if(apvNb<6)
      gainRatio=getGainRatio(iter->detid,apvNb);
    else
      edm::LogInfo("") << "[Doanalysis] detid " << iter->detid << " strip " << strip << " apvNb " << apvNb;
      }
      //edm::LogInfo("") << "[Doanalysis] detid " << iter->detid << " strip " << strip << " value " << iter->values[strip];
      value=iter->values[strip]*gainRatio;
      tkmap->fill(iter->detid,value);
      for(size_t i=0;i<histos.size();++i)
    histos[i]->Fill(value);
    }
  }
}

void SiStripCorrelateNoise::DoPlots ( )

private

Definition at line 85 of file SiStripCorrelateNoise.cc.

References patCaloMETCorrections_cff::C, dir, file, mps_fire::i, dataset::outName, TrackerMap::save(), theRun, tkmap, and vTH1.

Referenced by analyze(), and beginRun().

                              {
  TCanvas *C=new TCanvas();
  C->Divide(2,2);
  
  char outName[128];
  sprintf(outName,"Run_%d.png",theRun);
  for(size_t i=0;i<vTH1.size();i++)
    if(vTH1[i]!=0){
      if(i%100==0){
    C->cd(i/100);
    vTH1[i]->SetLineColor(i/100);
    vTH1[i]->Draw();
    C->cd(i/100)->SetLogy();
      }
      vTH1[i]->Write();
    }
  
  C->Print(outName);
  delete C;
  
  vTH1.clear();
  file->cd("");

  char dir[128];
  sprintf(dir,"Run_%d_TkMap.png",theRun);
  tkmap->save(false,0,5,dir);
  delete tkmap;
  tkmap = new TrackerMap();
}

void SiStripCorrelateNoise::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 237 of file SiStripCorrelateNoise.cc.

References file.

Referenced by analyze().

                              {
  file->Write();
  file->Close();


}

unsigned long long SiStripCorrelateNoise::getGainCache ( const edm::EventSetup &  eSetup )

inlineprivate

Definition at line 67 of file SiStripCorrelateNoise.h.

References checkGainCache(), plotBeamSpotDB::first, edm::EventSetup::get(), getGainRatio(), and getMeanNoise().

Referenced by checkGainCache().

{ return eSetup.get<SiStripApvGainRcd>().cacheIdentifier();}

float SiStripCorrelateNoise::getGainRatio ( const uint32_t &  detid, const uint16_t &  apv  )

private

Definition at line 160 of file SiStripCorrelateNoise.cc.

References SiStripApvGain::getApvGain(), SiStripApvGain::getRange(), newGain, and oldGain.

Referenced by DoAnalysis(), and getGainCache().

                                                                             {

  SiStripApvGain::Range oldRange=oldGain->getRange(detid); 
  SiStripApvGain::Range newRange=newGain->getRange(detid); 

  if(oldRange.first==oldRange.second ||
     newRange.first==newRange.second)
    return 1.;

  return oldGain->getApvGain(apv,oldRange)/newGain->getApvGain(apv,newRange);

}

TH1F * SiStripCorrelateNoise::getHisto ( const long unsigned int &  index )

private

Definition at line 213 of file SiStripCorrelateNoise.cc.

References diffTreeTool::index, dataset::name, AlCaHLTBitMon_QueryRunRegistry::string, theRun, and vTH1.

Referenced by analyze(), and getHistos().

                                                             {
  if(vTH1.size()<index+1)
    vTH1.resize(index+1,0);
  
  if(vTH1[index]==0){
    char name[128];
    std::string SubD;
    if(index<200)
      SubD="TIB";
    else if(index<300)
      SubD="TID";
    else if(index<400)
      SubD="TOB";
    else 
      SubD="TEC";
    sprintf(name,"%d_%lu__%s",theRun,index,SubD.c_str());
    edm::LogInfo("")<<"[getHisto] creating index " << index << std::endl;
    vTH1[index]=new TH1F(name,name,200,-0.5,10.5);
  }
  
  return vTH1[index];
}

void SiStripCorrelateNoise::getHistos ( const uint32_t &  detid, const TrackerTopology *  tTopo, std::vector< TH1F * > &  histos  )

private

Definition at line 184 of file SiStripCorrelateNoise.cc.

References a, getHisto(), diffTreeTool::index, DetId::subdetId(), TrackerTopology::tecSide(), TrackerTopology::tecWheel(), TrackerTopology::tibLayer(), TrackerTopology::tidSide(), TrackerTopology::tidWheel(), and TrackerTopology::tobLayer().

Referenced by analyze(), and DoAnalysis().

                                                                                                             {
  
  histos.clear();
  
  int subdet=-999; int component=-999;
  SiStripDetId a(detid);
  if ( a.subdetId() == 3 ){
    subdet=0;
    component=tTopo->tibLayer(detid);
  } else if ( a.subdetId() == 4 ) {
    subdet=1;
    component=tTopo->tidSide(detid)==2?tTopo->tidWheel(detid):tTopo->tidWheel(detid)+3;
  } else if ( a.subdetId() == 5 ) {
    subdet=2;
    component=tTopo->tobLayer(detid);
  } else if ( a.subdetId() == 6 ) {
    subdet=3;
    component=tTopo->tecSide(detid)==2?tTopo->tecWheel(detid):tTopo->tecWheel(detid)+9;
  } 
  
  int index=100+subdet*100+component;
  

  histos.push_back(getHisto(100+100*subdet));
  histos.push_back(getHisto(index));
  
}

float SiStripCorrelateNoise::getMeanNoise ( const SiStripNoises::Range &  noiseRange, const uint32_t &  first, const uint32_t &  range  )

private

Definition at line 174 of file SiStripCorrelateNoise.cc.

References SiStripNoises::getNoise(), RecoTauDiscriminantConfiguration::mean, and noiseHandle_.

Referenced by getGainCache().

                                                                                                                         {
  
  float mean=0;
  for (size_t istrip=firstStrip;istrip<firstStrip+range;istrip++){
    mean+=noiseHandle_->getNoise(istrip,noiseRange);
  }
  return mean/(1.*range);
}

unsigned long long SiStripCorrelateNoise::getNoiseCache ( const edm::EventSetup &  eSetup )

inlineprivate

Definition at line 66 of file SiStripCorrelateNoise.h.

References edm::EventSetup::get().

Referenced by beginRun().

{ return eSetup.get<SiStripNoisesRcd>().cacheIdentifier();}

Member Data Documentation

unsigned long long SiStripCorrelateNoise::cacheID_gain

private

Definition at line 98 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache().

unsigned long long SiStripCorrelateNoise::cacheID_noise

private

Definition at line 97 of file SiStripCorrelateNoise.h.

Referenced by beginRun().

bool SiStripCorrelateNoise::equalGain

private

Definition at line 88 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and DoAnalysis().

TFile* SiStripCorrelateNoise::file

private

Definition at line 91 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), DoPlots(), endJob(), ztee.GZipLog::finish(), and SiStripCorrelateNoise().

SiStripDetInfoFileReader* SiStripCorrelateNoise::fr

private

Definition at line 80 of file SiStripCorrelateNoise.h.

Referenced by SiStripCorrelateNoise().

edm::ESHandle<SiStripApvGain> SiStripCorrelateNoise::gainHandle_

private

Definition at line 81 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache().

SiStripApvGain * SiStripCorrelateNoise::newGain

private

Definition at line 87 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and getGainRatio().

edm::ESHandle<SiStripNoises> SiStripCorrelateNoise::noiseHandle_

private

Definition at line 82 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), and getMeanNoise().

SiStripApvGain* SiStripCorrelateNoise::oldGain

private

Definition at line 87 of file SiStripCorrelateNoise.h.

Referenced by checkGainCache(), and getGainRatio().

SiStripNoises* SiStripCorrelateNoise::refNoise

private

Definition at line 85 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), and DoAnalysis().

uint32_t SiStripCorrelateNoise::theRun

private

Definition at line 84 of file SiStripCorrelateNoise.h.

Referenced by beginRun(), DoPlots(), and getHisto().

TrackerMap* SiStripCorrelateNoise::tkmap

private

Definition at line 94 of file SiStripCorrelateNoise.h.

Referenced by DoAnalysis(), DoPlots(), and SiStripCorrelateNoise().

std::vector<TH1F*> SiStripCorrelateNoise::vTH1

private

Definition at line 92 of file SiStripCorrelateNoise.h.

Referenced by DoPlots(), and getHisto().