BuildTrackerMapPlugin Class Reference

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

Inheritance diagram for BuildTrackerMapPlugin:

Public Member Functions
	BuildTrackerMapPlugin (const edm::ParameterSet &)
	~BuildTrackerMapPlugin ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()
void	read (bool aMechView, std::string aFile, std::vector< TkHistoMap * > &aTkMapVec, std::vector< bool > &aValidVec)
void	subtractMap (TkHistoMap *&aResult, 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< TkHistoMap * >	tkHistoMapVec_
std::vector< TkHistoMap * >	tkHistoMapVecDiff_
std::vector< TrackerMap * >	tkmap_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
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 54 of file BuildTrackerMap.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 107 of file BuildTrackerMap.cc.

References doDiff_, fileName_, fileNameDiff_, isValidMap_, isValidMapDiff_, mechanicalView_, read(), tkHistoMapVec_, and tkHistoMapVecDiff_.

  : 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"))
{

  read(mechanicalView_,fileName_,tkHistoMapVec_,isValidMap_);
  if (doDiff_) read(mechanicalView_,fileNameDiff_,tkHistoMapVecDiff_,isValidMapDiff_);


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

(

)

Definition at line 138 of file BuildTrackerMap.cc.

References doDiff_, tkHistoMapVec_, and tkHistoMapVecDiff_.

{

  tkHistoMapVec_.clear();
  if (doDiff_) tkHistoMapVecDiff_.clear();
}

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 310 of file BuildTrackerMap.cc.

References gather_cfg::cout, create_complex_test_file0_cfg::firstEvent, edm::EventSetup::get(), i, pset_, tkHistoMapNameVec_, and tkmap_.

{

  static bool firstEvent = true;

  edm::ESHandle<SiStripFedCabling> fedcabling;
  iSetup.get<SiStripFedCablingRcd>().get(fedcabling );
  
  if (firstEvent) {
    for (unsigned int i(0); i<tkHistoMapNameVec_.size(); i++){
      tkmap_.push_back(new TrackerMap(pset_,&(*fedcabling)));
    }

  }

  firstEvent = false;

  std::cout << "End of analyze method: tkmap_ size = " << tkmap_.size() << std::endl;

}//analyze method

void BuildTrackerMapPlugin::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 334 of file BuildTrackerMap.cc.

{

}

void BuildTrackerMapPlugin::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 341 of file BuildTrackerMap.cc.

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

{
  //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),tkHistoMapVecDiff_.at(i));


    //(pset_,pDD1); 
    lTkMap->setPalette(1);
    lTkMap->showPalette(1);
    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, std::vector< TkHistoMap * > &  aTkMapVec, std::vector< bool > &  aValidVec  )

private

Definition at line 151 of file BuildTrackerMap.cc.

References gather_cfg::cout, TrackerOfflineValidation_Dqm_cff::dirName, cmsRelvalreport::exit, folderName_, i, DQMStore::open(), cppFunctionSkipper::operator, AlCaHLTBitMon_QueryRunRegistry::string, and tkHistoMapNameVec_.

Referenced by BuildTrackerMapPlugin(), Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::goto(), and Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::setFilterBranches().

{
  
  DQMStore * lDqmStore = edm::Service<DQMStore>().operator->();
  lDqmStore->open(aFile);  
  std::vector<TkHistoMap *> tkHistoMap;

  unsigned int nHists = tkHistoMapNameVec_.size();
  tkHistoMap.resize(nHists,0);
  aValidVec.resize(nHists,true);

  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++){

    tkHistoMap[i] = new TkHistoMap();

    tkHistoMap[i]->loadTkHistoMap(dirName,tkHistoMapNameVec_.at(i),aMechView);
    
    std::vector<MonitorElement*>& lMaps = tkHistoMap[i]->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.push_back(tkHistoMap[i]);
  }

  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, TkHistoMap *&  aSubtr  )

private

Definition at line 390 of file BuildTrackerMap.cc.

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

Referenced by endJob().

{
  
  std::vector<MonitorElement*>& lMaps = aResult->getAllMaps();
  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 83 of file BuildTrackerMap.cc.

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

std::string BuildTrackerMapPlugin::fileName_

private

Definition at line 75 of file BuildTrackerMap.cc.

Referenced by BuildTrackerMapPlugin().

std::string BuildTrackerMapPlugin::fileNameDiff_

private

Definition at line 84 of file BuildTrackerMap.cc.

Referenced by BuildTrackerMapPlugin().

std::string BuildTrackerMapPlugin::folderName_

private

Definition at line 79 of file BuildTrackerMap.cc.

Referenced by read().

std::vector<bool> BuildTrackerMapPlugin::isValidMap_

private

Definition at line 94 of file BuildTrackerMap.cc.

Referenced by BuildTrackerMapPlugin(), and endJob().

std::vector<bool> BuildTrackerMapPlugin::isValidMapDiff_

private

Definition at line 95 of file BuildTrackerMap.cc.

Referenced by BuildTrackerMapPlugin().

std::vector<double> BuildTrackerMapPlugin::maxVal_

private

Definition at line 92 of file BuildTrackerMap.cc.

Referenced by endJob().

bool BuildTrackerMapPlugin::mechanicalView_

private

Definition at line 77 of file BuildTrackerMap.cc.

Referenced by BuildTrackerMapPlugin().

std::vector<double> BuildTrackerMapPlugin::minVal_

private

Definition at line 91 of file BuildTrackerMap.cc.

Referenced by endJob().

unsigned int BuildTrackerMapPlugin::printDebug_

private

Definition at line 81 of file BuildTrackerMap.cc.

edm::ParameterSet BuildTrackerMapPlugin::pset_

private

Definition at line 97 of file BuildTrackerMap.cc.

Referenced by analyze().

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

private

Definition at line 90 of file BuildTrackerMap.cc.

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

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

private

Definition at line 86 of file BuildTrackerMap.cc.

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

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

private

Definition at line 87 of file BuildTrackerMap.cc.

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

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

private

Definition at line 98 of file BuildTrackerMap.cc.

Referenced by analyze(), and endJob().