#include <DQM/SiStripMonitorHardware/plugins/BuildTrackerMap.cc>

Inheritance diagram for BuildTrackerMapPlugin:

Public Member Functions
	BuildTrackerMapPlugin (const edm::ParameterSet &)

	~BuildTrackerMapPlugin () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

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

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

void	endJob () override

void	read (bool aMechView, std::string aFile, const TkDetMap *tkDetMap, std::vector< std::unique_ptr< TkHistoMap >> &aTkMapVec, std::vector< bool > &aValidVec)

void	subtractMap (TkHistoMap aResult, const TkHistoMap aSubtr)

Private Attributes
bool	doDiff_

std::string	fileName_

std::string	fileNameDiff_

std::string	folderName_

std::vector< bool >	isValidMap_

std::vector< bool >	isValidMapDiff_

std::vector< double >	maxVal_

bool	mechanicalView_

std::vector< double >	minVal_

unsigned int	printDebug_

edm::ParameterSet	pset_

std::vector< std::string >	tkHistoMapNameVec_

std::vector< std::unique_ptr< TkHistoMap > >	tkHistoMapVec_

std::vector< std::unique_ptr< TkHistoMap > >	tkHistoMapVecDiff_

std::vector< TrackerMap * >	tkmap_

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

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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: DQM source application to monitor common mode for SiStrip data

Definition at line 55 of file BuildTrackerMap.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 109 of file BuildTrackerMap.cc.

   : fileName_(iConfig.getUntrackedParameter<std::string>("InputFileName","DQMStore.root")),
     mechanicalView_(iConfig.getUntrackedParameter<bool>("MechanicalView",true)),
     folderName_(iConfig.getUntrackedParameter<std::string>("HistogramFolderName","DQMData/")),
     printDebug_(iConfig.getUntrackedParameter<unsigned int>("PrintDebugMessages",1)),
     doDiff_(iConfig.getUntrackedParameter<bool>("DoDifference",false)),
     fileNameDiff_(iConfig.getUntrackedParameter<std::string>("InputFileNameForDiff","DQMStore.root")),
     tkHistoMapNameVec_(iConfig.getUntrackedParameter<std::vector<std::string> >("TkHistoMapNameVec")),
     minVal_(iConfig.getUntrackedParameter<std::vector<double> >("MinValueVec")),
     maxVal_(iConfig.getUntrackedParameter<std::vector<double> >("MaxValueVec")),
     pset_(iConfig.getParameter<edm::ParameterSet>("TkmapParameters"))
 {
 
 
 
 //   for (unsigned int i(0); i<34; i++){
 //     if (i<4) histName_[i] << "TIB/layer_" << i+1 << "/" << tkDetMapName_ << "_TIB_L" << i+1;
 //     else if (i<7)  histName_[i] << "TID/side_1/wheel_" << i-3 << "/" << tkDetMapName_ << "_TIDM_D" << i-3;
 //     else if (i<10)  histName_[i] << "TID/side_2/wheel_" << i-6 << "/" << tkDetMapName_ << "_TIDP_D" << i-6;
 //     else if (i<16) histName_[i] << "TOB/layer_" << i-9 << "/" << tkDetMapName_ << "_TOB_L" << i-9;
 //     else if (i<25) histName_[i] << "TEC/side_1/wheel_" << i-15 << "/" << tkDetMapName_ << "_TECM_W" << i-15;
 //     else if (i<34) histName_[i] << "TEC/side_2/wheel_" << i-24 << "/" << tkDetMapName_ << "_TECP_W" << i-24;
 
 //     std::cout << "histName[" << i << "] =" << histName_[i] << std::endl;
 
 //   } 
   
 }

BuildTrackerMapPlugin::~BuildTrackerMapPlugin ( )

inlineoverride

Definition at line 60 of file BuildTrackerMap.cc.

References analyze(), beginJob(), endJob(), read(), AlCaHLTBitMon_QueryRunRegistry::string, and subtractMap().

60 {}

Member Function Documentation

void BuildTrackerMapPlugin::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 302 of file BuildTrackerMap.cc.

References gather_cfg::cout, doDiff_, fileName_, fileNameDiff_, single_neutrino_cfg::firstEvent, edm::EventSetup::get(), mps_fire::i, isValidMap_, isValidMapDiff_, mechanicalView_, edm::ESHandle< T >::product(), pset_, read(), tkHistoMapNameVec_, tkHistoMapVec_, tkHistoMapVecDiff_, and tkmap_.

Referenced by ~BuildTrackerMapPlugin().

 {
 
   static bool firstEvent = true;
 
   edm::ESHandle<SiStripFedCabling> fedcabling;
   iSetup.get<SiStripFedCablingRcd>().get(fedcabling );
   
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   if ( tkHistoMapVec_.empty() && ( ! tkHistoMapNameVec_.empty() ) ) {
     edm::ESHandle<TkDetMap> tkDetMapHandle;
     iSetup.get<TrackerTopologyRcd>().get(tkDetMapHandle);
     const TkDetMap* tkDetMap = tkDetMapHandle.product();
     read(mechanicalView_, fileName_, tkDetMap, tkHistoMapVec_, isValidMap_);
     if (doDiff_) read(mechanicalView_, fileNameDiff_, tkDetMap, tkHistoMapVecDiff_, isValidMapDiff_);
   }
 
   if (firstEvent) {
     for (unsigned int i(0); i<tkHistoMapNameVec_.size(); i++){
       tkmap_.push_back(new TrackerMap(pset_,&(*fedcabling),tTopo));
     }
 
   }
 
   firstEvent = false;
 
   std::cout << "End of analyze method: tkmap_ size = " << tkmap_.size() << std::endl;
 
 }//analyze method

void BuildTrackerMapPlugin::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 338 of file BuildTrackerMap.cc.

Referenced by ~BuildTrackerMapPlugin().

 {
 
 }

void BuildTrackerMapPlugin::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 345 of file BuildTrackerMap.cc.

References gather_cfg::cout, mps_fire::i, isValidMap_, maxVal_, minVal_, TrackerMap::save(), TrackerMap::save_as_fedtrackermap(), TrackerMap::setPalette(), TrackerMap::showPalette(), AlCaHLTBitMon_QueryRunRegistry::string, subtractMap(), tkHistoMapNameVec_, tkHistoMapVec_, tkHistoMapVecDiff_, and tkmap_.

Referenced by ~BuildTrackerMapPlugin().

 {
   //edm::ESHandle<SiStripFedCabling> pDD1;
   //iSetup.get<SiStripFedCablingRcd>().get(pDD1);
   std::cout << "Processing endjob with " << tkHistoMapNameVec_.size()<< " elements." << std::endl;
 
   assert (minVal_.size() == tkHistoMapNameVec_.size());
   assert (maxVal_.size() == tkHistoMapNameVec_.size());
 
   for (unsigned int i(0); i<tkHistoMapNameVec_.size(); i++){
 
     std::cout << "Processing element " << i << ": " << tkHistoMapNameVec_.at(i) << std::endl;
     std::cout << "Min, max = " << minVal_.at(i) << " " << maxVal_.at(i) << std::endl;
 
     TrackerMap* lTkMap = tkmap_.at(i);
 
     if (!lTkMap) {
       std::cout << "tkmap_ is NULL for element " << i << "... continuing ..." << std::endl;
       continue;
     }
 
     subtractMap(tkHistoMapVec_.at(i).get(),tkHistoMapVecDiff_.at(i).get());
 
 
     //(pset_,pDD1); 
     lTkMap->setPalette(1);
     lTkMap->showPalette(true);
     if (!tkHistoMapVec_.at(i) || !isValidMap_.at(i)) {
       std::cout << "Warning, tkHistoMap is invalid for element " << i << "... continuing ..." << std::endl;
       continue;
     }
     tkHistoMapVec_.at(i)->dumpInTkMap(lTkMap);
 
     //to print all figures to create fancy view
     //lTkMap->printall(true,0,255,tkHistoMapNameVec_.at(i));
     lTkMap->save(true,
          minVal_.at(i),
          maxVal_.at(i),
          tkHistoMapNameVec_.at(i)+std::string(".png"));
     lTkMap->save_as_fedtrackermap(true,
                   minVal_.at(i),
                   maxVal_.at(i),
                   tkHistoMapNameVec_.at(i)+std::string("_FED.png"));
 
   }
 
 }

void BuildTrackerMapPlugin::read	(	bool	aMechView,
		std::string	aFile,
		const TkDetMap *	tkDetMap,
		std::vector< std::unique_ptr< TkHistoMap >> &	aTkMapVec,
		std::vector< bool > &	aValidVec
	)

private

Definition at line 143 of file BuildTrackerMap.cc.

References gather_cfg::cout, TrackerOfflineValidation_Dqm_cff::dirName, cmsRelvalreport::exit, folderName_, mps_fire::i, TkHistoMap::loadTkHistoMap(), eostools::move(), DQMStore::open(), Utilities::operator, AlCaHLTBitMon_QueryRunRegistry::string, and tkHistoMapNameVec_.

Referenced by analyze(), edmIntegrityCheck.PublishToFileSystem::get(), Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::goto(), Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::setFilterBranches(), and ~BuildTrackerMapPlugin().

 {
   
   DQMStore * lDqmStore = edm::Service<DQMStore>().operator->();
   lDqmStore->open(aFile);  
 
   unsigned int nHists = tkHistoMapNameVec_.size();
   aTkMapVec.reserve(nHists);
   aValidVec.reserve(nHists);
 
   std::string dirName = folderName_;
   if (dirName == "") {
     dirName += "Run ";
     dirName += aFile.substr(aFile.find_last_of("_")+5,6);
     dirName += "/SiStrip/Run summary";
     std::cout << " -- DirName = " << dirName << std::endl;
    }
 
   //lDqmStore->setCurrentFolder(dirName);
   //lDqmStore->showDirStructure();
   
   unsigned int nFailTot=0;
   unsigned int nTotTot = 0;
   for (unsigned int i(0); i<nHists; i++){
 
     std::unique_ptr<TkHistoMap> tkHistoMap{new TkHistoMap(tkDetMap)};
 
     tkHistoMap->loadTkHistoMap(dirName,tkHistoMapNameVec_.at(i),aMechView);
     
     std::vector<MonitorElement*>& lMaps = tkHistoMap->getAllMaps();
 
     std::cout << " -- map " << i << ", nHistos = " << lMaps.size() << std::endl;
     unsigned int nFail=0;
     unsigned int nTot=0;
  
     for (unsigned int im(0); im<lMaps.size(); im++){
       if (!lMaps[im]) {
     std::cout << " -- Failed to get element " << im << " for map " << i << std::endl;
     nFail++;
     nFailTot++;
       }
       nTot++;
       nTotTot++;
     }
 
     if (nFail == nTot) aValidVec[i] = false;
     aTkMapVec.emplace_back(std::move(tkHistoMap));
   }
 
   if (nFailTot < nTotTot) std::cout << " - " << nTotTot-nFailTot << "/" << nTotTot 
                   << " histomaps read with success for file ." << aFile << std::endl;
   else {
     std::cout << " - Failed to read any map for file " << aFile << ". Exiting line ... " << __LINE__ << std::endl;
      exit(1);
   }
 
 //   //get list of detid for which |deltaRMS(APV0-APV1)|>1
 //   unsigned int lHistoNumber = 35;
 //   TkDetMap lTkDetMap;
 //   std::ofstream list,listRms0,listRms1;
 //   list.open("./cmBadModuleList.dat",std::ios::out);
 //   listRms0.open("./cmBadModuleList_rms0.dat",std::ios::out);
 //   listRms1.open("./cmBadModuleList_rms1.dat",std::ios::out);
 //   if (!list || !listRms0 || !listRms1) {
 //     std::cout << "Warning, can't open output file to write bad module list !" << std::endl;
 //     exit(1);
 //   }
 
 //   TCanvas *lCan = new TCanvas("lCan","",1);
 //   TH1F *p_deltaMean = new TH1F("p_deltaMean",";CM_{mean}(APV0)-CM_{mean}(APV1)",500,-2,2);
 //   TH1F *p_deltaRMS = new TH1F("p_deltaRMS",";CM_{RMS}(APV0)-CM_{RMS}(APV1)",500,0,3);
 //   TH1F *p_MeanAPV0 = new TH1F("p_MeanAPV0",";CM_{mean}(APV0)",500,100,140);
 //   //TH1F *p_MeanAPV1 = new TH1F("p_MeanAPV1",";CM_{mean}(APV1)",500,100,140);
 //   TH1F *p_RMSAPV0 = new TH1F("p_RMSAPV0",";CM_{RMS}(APV0)",500,0,10);
 //   //TH1F *p_RMSAPV1 = new TH1F("p_RMSAPV1",";CM_{RMS}(APV1)",500,0,10);
 
 
 
 //   gStyle->SetOptStat(1111111);
 
 //   for(unsigned int layer=1;layer<lHistoNumber;++layer){
 //     std::vector<uint32_t> dets;
 //     lTkDetMap.getDetsForLayer(layer,dets);
 //     for(size_t i=0;i<dets.size();++i){
 //       if(dets[i]>0){
 //  //if(tkHistoMap[5]->getEntries(dets[i])>0 && tkHistoMap[5]->getValue(dets[i])) {
 //  if(nHists > 3){
 //    if (tkHistoMap[3]->getValue(dets[i]) > 1) {
 //      list << dets[i] << " " << tkHistoMap[3]->getValue(dets[i]) << std::endl;
 //    }
 //  p_deltaRMS->Fill(tkHistoMap[3]->getValue(dets[i]));
 //  }
 //  p_MeanAPV0->Fill(tkHistoMap[0]->getValue(dets[i]));
 //  //p_MeanAPV1->Fill(tkHistoMap[1]->getValue(dets[i]));
 //  p_RMSAPV0->Fill(tkHistoMap[1]->getValue(dets[i]));
 //  if (tkHistoMap[1]->getValue(dets[i]) > 2)
 //    listRms0 << dets[i] << " " << tkHistoMap[1]->getValue(dets[i]) << std::endl;
 //  //p_RMSAPV1->Fill(tkHistoMap[3]->getValue(dets[i]));
 //  //if (tkHistoMap[3]->getValue(dets[i]) > 2)
 //  //listRms1 << dets[i] << " " << tkHistoMap[3]->getValue(dets[i]) << std::endl;
 
 //  if(nHists > 2) p_deltaMean->Fill(tkHistoMap[2]->getValue(dets[i]));
 //       }
 //     }
 //   }
 //   list.close();
 //   listRms0.close();
 //   listRms1.close();
 
 //   lCan->cd();
 //   p_deltaRMS->Draw();
 //   //lCan->Print("./deltaRMStotal.png");
 //   lCan->Print("./deltaRMStotal.C");
 
 //   p_deltaMean->Draw();
 //   lCan->Update();
 //   lCan->Print("./deltaMeantotal.C");
 
 //   TPaveStats *statBox[2] = {0,0};  
 //   statBox[0] = (TPaveStats*)p_MeanAPV0->FindObject("stats");
 //   //statBox[1] = (TPaveStats*)p_MeanAPV1->FindObject("stats");
 
 //   p_MeanAPV0->Draw();
 //   //p_MeanAPV1->SetLineColor(2);
 //   //p_MeanAPV1->Draw("same");
 //   if (statBox[0]) statBox[0]->Draw("same");
 //   if (statBox[1]) { 
 //     statBox[1]->SetLineColor(2);
 //     statBox[1]->SetTextColor(2);
 //     statBox[1]->Draw("same");
 //   }
 //   lCan->Update();
 //   lCan->Print("./meanAPVstotal.C");
 
 //   statBox[0] = (TPaveStats*)p_RMSAPV0->FindObject("stats");
 //   //statBox[1] = (TPaveStats*)p_RMSAPV1->FindObject("stats");
 
 //   p_RMSAPV0->Draw();
 //   //p_RMSAPV1->SetLineColor(2);
 //   //p_RMSAPV1->Draw("same");
 //   if (statBox[0]) statBox[0]->Draw("same");
 //   if (statBox[1]) { 
 //     statBox[1]->SetLineColor(2);
 //     statBox[1]->SetTextColor(2);
 //     statBox[1]->Draw("same");
 //   }
 //   lCan->Update();
 //   lCan->Print("./rmsAPVstotal.C");
 
 
 }

void BuildTrackerMapPlugin::subtractMap	(	TkHistoMap *	aResult,
		const TkHistoMap *	aSubtr
	)

private

Definition at line 394 of file BuildTrackerMap.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, and TkHistoMap::getAllMaps().

Referenced by endJob(), and ~BuildTrackerMapPlugin().

 {
   
   std::vector<MonitorElement*>& lMaps = aResult->getAllMaps();
   const std::vector<MonitorElement*>& lMapsDiff = aSubtr->getAllMaps();
 
   assert(lMaps.size() == lMapsDiff.size());  
 
     for (unsigned int im(0); im<lMaps.size(); im++){
       if (!lMaps[im] || !lMapsDiff[im]) {
     std::cout << " -- Failed to get element " << im << " for maps." << std::endl;
       }
       else {
     (lMaps[im]->getTProfile2D())->Add(lMapsDiff[im]->getTProfile2D(),-1);
       }
     }
 
 }

Member Data Documentation

bool BuildTrackerMapPlugin::doDiff_

private

Definition at line 85 of file BuildTrackerMap.cc.

Referenced by analyze().

std::string BuildTrackerMapPlugin::fileName_

private

Definition at line 77 of file BuildTrackerMap.cc.

Referenced by analyze().

std::string BuildTrackerMapPlugin::fileNameDiff_

private

Definition at line 86 of file BuildTrackerMap.cc.

Referenced by analyze().

std::string BuildTrackerMapPlugin::folderName_

private

Definition at line 81 of file BuildTrackerMap.cc.

Referenced by read().

std::vector<bool> BuildTrackerMapPlugin::isValidMap_

private

Definition at line 96 of file BuildTrackerMap.cc.

Referenced by analyze(), and endJob().

std::vector<bool> BuildTrackerMapPlugin::isValidMapDiff_

private

Definition at line 97 of file BuildTrackerMap.cc.

Referenced by analyze().

std::vector<double> BuildTrackerMapPlugin::maxVal_

private

Definition at line 94 of file BuildTrackerMap.cc.

Referenced by endJob().

bool BuildTrackerMapPlugin::mechanicalView_

private

Definition at line 79 of file BuildTrackerMap.cc.

Referenced by analyze().

std::vector<double> BuildTrackerMapPlugin::minVal_

private

Definition at line 93 of file BuildTrackerMap.cc.

Referenced by endJob().

unsigned int BuildTrackerMapPlugin::printDebug_

private

Definition at line 83 of file BuildTrackerMap.cc.

edm::ParameterSet BuildTrackerMapPlugin::pset_

private

Definition at line 99 of file BuildTrackerMap.cc.

Referenced by analyze().

std::vector<std::string> BuildTrackerMapPlugin::tkHistoMapNameVec_

private

Definition at line 92 of file BuildTrackerMap.cc.

Referenced by analyze(), endJob(), and read().

std::vector<std::unique_ptr<TkHistoMap> > BuildTrackerMapPlugin::tkHistoMapVec_

private

Definition at line 88 of file BuildTrackerMap.cc.

Referenced by analyze(), and endJob().

std::vector<std::unique_ptr<TkHistoMap> > BuildTrackerMapPlugin::tkHistoMapVecDiff_

private

Definition at line 89 of file BuildTrackerMap.cc.

Referenced by analyze(), and endJob().

std::vector<TrackerMap*> BuildTrackerMapPlugin::tkmap_

private

Definition at line 100 of file BuildTrackerMap.cc.

Referenced by analyze(), and endJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation