CMS 3D CMS Logo
Public Member Functions | Protected Member Functions | Private Types | Private Attributes

AlcaBeamMonitorClient Class Reference

#include <AlcaBeamMonitorClient.h>

Inheritance diagram for AlcaBeamMonitorClient:
edm::EDAnalyzer

List of all members.

Public Member Functions

 AlcaBeamMonitorClient (const edm::ParameterSet &)
 ~AlcaBeamMonitorClient ()

Protected Member Functions

void analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup)
void beginJob (void)
void beginLuminosityBlock (const edm::LuminosityBlock &iLumi, const edm::EventSetup &iSetup)
void beginRun (const edm::Run &iRun, const edm::EventSetup &iSetup)
void endJob (const edm::LuminosityBlock &iLumi, const edm::EventSetup &iSetup)
void endLuminosityBlock (const edm::LuminosityBlock &iLumi, const edm::EventSetup &iSetup)
void endRun (const edm::Run &iRun, const edm::EventSetup &iSetup)

Private Types

typedef std::map< std::string,
std::map< std::string,
std::map< std::string,
MonitorElement * > > > 		HistosContainer
typedef std::map< std::string,
std::map< std::string,
std::map< std::string, int > > > 		PositionContainer

Private Attributes

DQMStore * dbe_
std::multimap< std::string,
std::string > 		histoByCategoryNames_
HistosContainer histosMap_
std::string monitorName_
int numberOfValuesToSave_
edm::ParameterSet parameters_
PositionContainer positionsMap_
std::map
< edm::LuminosityBlockNumber_t,
std::vector< double > > 		valuesMap_
std::vector< std::string > varNamesV_

Detailed Description

*

Date:
2010/09/24 06:36:05
Revision:
1.2
Author:
Lorenzo Uplegger/FNAL

Definition at line 26 of file AlcaBeamMonitorClient.h.

Member Typedef Documentation

typedef std::map<std::string,std::map<std::string,std::map<std::string,MonitorElement*> > > AlcaBeamMonitorClient::HistosContainer [private]

Definition at line 43 of file AlcaBeamMonitorClient.h.

typedef std::map<std::string,std::map<std::string,std::map<std::string,int> > > AlcaBeamMonitorClient::PositionContainer [private]

Definition at line 46 of file AlcaBeamMonitorClient.h.

Constructor & Destructor Documentation

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

Definition at line 26 of file AlcaBeamMonitorClient.cc.

References dbe_, histoByCategoryNames_, histosMap_, monitorName_, numberOfValuesToSave_, cmsCodeRules::cppFunctionSkipper::operator, positionsMap_, and varNamesV_.

                                                                     :
  parameters_         (ps),
  monitorName_        (parameters_.getUntrackedParameter<string>("MonitorName","YourSubsystemName")),
  numberOfValuesToSave_ (0)
{
  dbe_            = Service<DQMStore>().operator->();
  
  if (monitorName_ != "" ) monitorName_ = monitorName_+"/" ;
  

  varNamesV_.push_back("x");
  varNamesV_.push_back("y");
  varNamesV_.push_back("z");
  varNamesV_.push_back("sigmaX");
  varNamesV_.push_back("sigmaY");
  varNamesV_.push_back("sigmaZ");

/*
  histoByCategoryNames_.insert( pair<string,string>("run",        "Coordinate"));
  histoByCategoryNames_.insert( pair<string,string>("run",        "PrimaryVertex fit-DataBase"));
  histoByCategoryNames_.insert( pair<string,string>("run",        "PrimaryVertex fit-BeamFit"));
  histoByCategoryNames_.insert( pair<string,string>("run",        "PrimaryVertex fit-Scalers"));
  histoByCategoryNames_.insert( pair<string,string>("run",        "PrimaryVertex-DataBase"));
  histoByCategoryNames_.insert( pair<string,string>("run",        "PrimaryVertex-BeamFit"));
  histoByCategoryNames_.insert( pair<string,string>("run",        "PrimaryVertex-Scalers"));
*/
  histoByCategoryNames_.insert( pair<string,string>("lumi",       "Lumibased BeamSpotFit"));  
  histoByCategoryNames_.insert( pair<string,string>("lumi",       "Lumibased PrimaryVertex"));
  histoByCategoryNames_.insert( pair<string,string>("lumi",       "Lumibased DataBase"));     
  histoByCategoryNames_.insert( pair<string,string>("lumi",       "Lumibased Scalers"));      
  histoByCategoryNames_.insert( pair<string,string>("lumi",       "Lumibased PrimaryVertex-DataBase fit"));
  histoByCategoryNames_.insert( pair<string,string>("lumi",       "Lumibased PrimaryVertex-Scalers fit"));
  histoByCategoryNames_.insert( pair<string,string>("validation", "Lumibased Scalers-DataBase fit"));
  histoByCategoryNames_.insert( pair<string,string>("validation", "Lumibased PrimaryVertex-DataBase"));
  histoByCategoryNames_.insert( pair<string,string>("validation", "Lumibased PrimaryVertex-Scalers"));


  for(vector<string>::iterator itV=varNamesV_.begin(); itV!=varNamesV_.end(); itV++){
    for(multimap<string,string>::iterator itM=histoByCategoryNames_.begin(); itM!=histoByCategoryNames_.end(); itM++){
      histosMap_[*itV][itM->first][itM->second] = 0;
      positionsMap_[*itV][itM->first][itM->second] = 3*numberOfValuesToSave_;//value, error, ok 
      ++numberOfValuesToSave_;
    }
  }
}
AlcaBeamMonitorClient::~AlcaBeamMonitorClient ( )

Definition at line 73 of file AlcaBeamMonitorClient.cc.

                                              {
}

Member Function Documentation

void AlcaBeamMonitorClient::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
) [protected, virtual]

Implements edm::EDAnalyzer.

Definition at line 112 of file AlcaBeamMonitorClient.cc.

                                                                                 {
}
void AlcaBeamMonitorClient::beginJob ( void  ) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file AlcaBeamMonitorClient.cc.

                                     {
}
void AlcaBeamMonitorClient::beginLuminosityBlock ( const edm::LuminosityBlock &  iLumi,
const edm::EventSetup &  iSetup 
) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 108 of file AlcaBeamMonitorClient.cc.

                                                                                                       {
}
void AlcaBeamMonitorClient::beginRun ( const edm::Run &  iRun,
const edm::EventSetup &  iSetup 
) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 82 of file AlcaBeamMonitorClient.cc.

References dbe_, histosMap_, monitorName_, and DQMStore::removeElement().

                                                                               {
  for(HistosContainer::iterator itM=histosMap_.begin(); itM!=histosMap_.end(); itM++){
    for(map<string,map<string,MonitorElement*> >::iterator itMM=itM->second.begin(); itMM!=itM->second.end(); itMM++){
      if(itMM->first != "run"){
        for(map<string,MonitorElement*>::iterator itMMM=itMM->second.begin(); itMMM!=itMM->second.end(); itMMM++){
          if( itMMM->second != 0){
            if(itMM->first == "lumi"){
              dbe_->removeElement(monitorName_+"Debug",itMMM->second->getName());
            }
            else if(itMM->first == "validation"){
              dbe_->removeElement(monitorName_+"Validation",itMMM->second->getName());
            }
            else{
              LogInfo("AlcaBeamMonitorClient") 
                << "Unrecognized category " << itMM->first;
                //assert(0);        
            }
            itMMM->second = 0;
          }
        }
      }
    }
  }
}
void AlcaBeamMonitorClient::endJob ( const edm::LuminosityBlock &  iLumi,
const edm::EventSetup &  iSetup 
) [protected]

Definition at line 346 of file AlcaBeamMonitorClient.cc.

                                                                                        {
}
void AlcaBeamMonitorClient::endLuminosityBlock ( const edm::LuminosityBlock &  iLumi,
const edm::EventSetup &  iSetup 
) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 117 of file AlcaBeamMonitorClient.cc.

References dbe_, DQMStore::get(), MonitorElement::getTProfile(), i, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), monitorName_, numberOfValuesToSave_, tmp, and valuesMap_.

                                                                                                     {
  MonitorElement * tmp = 0;     
  tmp = dbe_->get(monitorName_+"Service/hHistoLumiValues");
  if(!tmp){
    return;
  }
  valuesMap_[iLumi.id().luminosityBlock()] = vector<double>();
  for(int i=0; i<3*numberOfValuesToSave_; i++){
    valuesMap_[iLumi.id().luminosityBlock()].push_back(tmp->getTProfile()->GetBinContent(i+1)); 
  }
}
void AlcaBeamMonitorClient::endRun ( const edm::Run &  iRun,
const edm::EventSetup &  iSetup 
) [protected, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 130 of file AlcaBeamMonitorClient.cc.

References newFWLiteAna::bin, DQMStore::book1D(), dbe_, first, histosMap_, edm::RunBase::id(), max(), min, monitorName_, mergeVDriftHistosByStation::name, positionsMap_, edm::RunID::run(), edm::second(), DQMStore::setCurrentFolder(), indexGen::title, and valuesMap_.

                                                                            {

  // use this in case any LS is missing.
  if(valuesMap_.size() == 0){
    LogInfo("AlcaBeamMonitorClient") 
      << "The histogram " << monitorName_+"Service/hHistoLumiValues which contains all the values has not been found in any lumi!";
    return;
  }
  int lastLumi = (--valuesMap_.end())->first;
  int firstLumi = valuesMap_.begin()->first;

  // create and cd into new folder
  dbe_->setCurrentFolder(monitorName_+"Validation");

  LogInfo("AlcaBeamMonitorClient") 
    << "End of run " << iRun.id().run() << "(" << firstLumi << "-" << lastLumi << ")";

  string name;
  string title;
  //x,y,z,sigmaX,sigmaY,sigmaZ
  for(HistosContainer::iterator itM=histosMap_.begin(); itM!=histosMap_.end(); itM++){
    for(map<string,map<string,MonitorElement*> >::iterator itMM=itM->second.begin(); itMM!=itM->second.end(); itMM++){
      if(itMM->first != "run"){
        for(map<string,MonitorElement*>::iterator itMMM=itMM->second.begin(); itMMM!=itMM->second.end(); itMMM++){
          name = string("h") + itM->first + itMMM->first;
          title = itM->first + "_{0} " + itMMM->first;
          if(itMM->first == "lumi"){
            dbe_->setCurrentFolder(monitorName_+"Debug");
            itMMM->second = dbe_->book1D(name,title,lastLumi-firstLumi+1,firstLumi-0.5,lastLumi+0.5);
          }
          else if(itMM->first == "validation" && itMMM->first == "Lumibased Scalers-DataBase fit"){
            dbe_->setCurrentFolder(monitorName_+"Validation");
            itMMM->second = dbe_->book1D(name,title,lastLumi-firstLumi+1,firstLumi-0.5,lastLumi+0.5);
          }
          else if(itMM->first == "validation" && itMMM->first != "Lumibased Scalers-DataBase fit" && (itM->first == "x" || itM->first == "y" || itM->first == "z")){
//LORE    else if(itMM->first == "validation" && itMMM->first != "Lumibased Scalers-DataBase fit" && (itM->first == "x" || itM->first == "y") ){
            dbe_->setCurrentFolder(monitorName_+"Validation");
            itMMM->second = dbe_->book1D(name,title,lastLumi-firstLumi+1,firstLumi-0.5,lastLumi+0.5);
//LORE      itMMM->second = dbe_->bookProfile(name,title,lastLumi-firstLumi+1,firstLumi-0.5,lastLumi+0.5,-0.1,0.1,"");
          }
//LORE    else if(itMM->first == "validation" && itMMM->first != "Lumibased Scalers-DataBase fit" && itM->first == "z" ){
//LORE      dbe_->setCurrentFolder(monitorName_+"Validation");
//LORE      itMMM->second = dbe_->book1D(name,title,lastLumi-firstLumi+1,firstLumi-0.5,lastLumi+0.5);
//LORE      itMMM->second = dbe_->bookProfile(name,title,lastLumi-firstLumi+1,firstLumi-0.5,lastLumi+0.5,-1,1,"");
//        }
//        else if(itMM->first == "validation" && itMMM->first == "Lumibased Scalers-DataBase" && (itM->first == "sigmaX" || itM->first == "sigmaY") ){
//          dbe_->setCurrentFolder(monitorName_+"Validation");
//          itMMM->second = dbe_->bookProfile(name,title,lastLumi-firstLumi+1,firstLumi-0.5,lastLumi+0.5,-0.1,0.1,"");
//        }
//        else if(itMM->first == "validation" && itMMM->first == "Lumibased Scalers-DataBase" && (itM->first == "sigmaZ") ){
//          dbe_->setCurrentFolder(monitorName_+"Validation");
//          itMMM->second = dbe_->bookProfile(name,title,lastLumi-firstLumi+1,firstLumi-0.5,lastLumi+0.5,-10,10,"");
//        }
//        else if(itMM->first == "validation" && itMMM->first != "Lumibased Scalers-DataBase" && (itM->first == "sigmaX" || itM->first == "sigmaY" || itM->first == "sigmaZ") ){
          else if(itMM->first == "validation" && (itM->first == "sigmaX" || itM->first == "sigmaY" || itM->first == "sigmaZ") ){
            dbe_->setCurrentFolder(monitorName_+"Validation");
            itMMM->second = 0;
          }
          else{
            LogInfo("AlcaBeamMonitorClient") 
              << "Unrecognized category " << itMM->first;
              //assert(0);          
          }
          if(itMMM->second != 0){
            if(itMMM->first.find('-') != string::npos){                                     
              itMMM->second->setAxisTitle(string("#Delta ") + itM->first + "_{0} (cm)",2);  
            }
            else{
              itMMM->second->setAxisTitle(itM->first + "_{0} (cm)",2);  
            }
            itMMM->second->setAxisTitle("Lumisection",1);
          }
        }
      }           
    }
  }

  unsigned int bin=0;
  for(HistosContainer::iterator itH=histosMap_.begin(); itH!=histosMap_.end(); itH++){
    for(map<string, map<string,MonitorElement*> >::iterator itHH=itH->second.begin(); itHH!=itH->second.end(); itHH++){
      for(map<string,MonitorElement*>::iterator itHHH=itHH->second.begin(); itHHH!=itHH->second.end(); itHHH++){
        for(map<LuminosityBlockNumber_t,vector<double> >::iterator itVal = valuesMap_.begin(); itVal != valuesMap_.end(); itVal++){
          if(itHHH->second != 0){
//          cout << positionsMap_[itH->first][itHH->first][itHHH->first] << endl;
            if(itVal->second[positionsMap_[itH->first][itHH->first][itHHH->first]+2] == 1){
              bin = itHHH->second->getTH1()->FindBin(itVal->first);
              itHHH->second->setBinContent(bin,itVal->second[positionsMap_[itH->first][itHH->first][itHHH->first]]);
              itHHH->second->setBinError(bin,itVal->second[positionsMap_[itH->first][itHH->first][itHHH->first]+1]);
            }
          }
        }
      }
    }
  }


  const double bigNumber = 1000000.;
  for(HistosContainer::iterator itH=histosMap_.begin(); itH!=histosMap_.end(); itH++){
    for(map<string, map<string,MonitorElement*> >::iterator itHH=itH->second.begin(); itHH!=itH->second.end(); itHH++){
      double min = bigNumber;
      double max = -bigNumber;
      double minDelta = bigNumber;
      double maxDelta = -bigNumber;
//      double minDeltaProf = bigNumber;
//      double maxDeltaProf = -bigNumber;
      if(itHH->first != "run"){
        for(map<string,MonitorElement*>::iterator itHHH=itHH->second.begin(); itHHH!=itHH->second.end(); itHHH++){
          if(itHHH->second != 0){
            for(int bin=1; bin<=itHHH->second->getTH1()->GetNbinsX(); bin++){
              if(itHHH->second->getTH1()->GetBinError(bin) != 0 || itHHH->second->getTH1()->GetBinContent(bin) != 0){
                if(itHHH->first == "Lumibased BeamSpotFit" 
                || itHHH->first == "Lumibased PrimaryVertex" 
                || itHHH->first == "Lumibased DataBase" 
                || itHHH->first == "Lumibased Scalers"){
                  if(min > itHHH->second->getTH1()->GetBinContent(bin)){                                                   
                    min = itHHH->second->getTH1()->GetBinContent(bin);
                  }                                                                                                        
                  if(max < itHHH->second->getTH1()->GetBinContent(bin)){                                                   
                    max = itHHH->second->getTH1()->GetBinContent(bin);
                  }                                                                                                        
                }
                else if(itHHH->first == "Lumibased PrimaryVertex-DataBase fit" 
                || itHHH->first == "Lumibased PrimaryVertex-Scalers fit"
                || itHHH->first == "Lumibased Scalers-DataBase fit"
                || itHHH->first == "Lumibased PrimaryVertex-DataBase" 
                || itHHH->first == "Lumibased PrimaryVertex-Scalers"
                ){
                  if(minDelta > itHHH->second->getTH1()->GetBinContent(bin)){
                    minDelta = itHHH->second->getTH1()->GetBinContent(bin);
                  }
                  if(maxDelta < itHHH->second->getTH1()->GetBinContent(bin)){
                    maxDelta = itHHH->second->getTH1()->GetBinContent(bin);
                  }
                }
//              else if(itHHH->first == "Lumibased PrimaryVertex-DataBase" 
//              || itHHH->first == "Lumibased PrimaryVertex-Scalers"
//              ){
//                if(minDelta > itHHH->second->getTProfile()->GetBinContent(bin)){
//                  minDelta = itHHH->second->getTProfile()->GetBinContent(bin);
//                }
//                if(maxDelta < itHHH->second->getTProfile()->GetBinContent(bin)){
//                  maxDelta = itHHH->second->getTProfile()->GetBinContent(bin);
//                }
//              }
                else{
                  LogInfo("AlcaBeamMonitorClient")
                    << "The histosMap_ have a histogram named " << itHHH->first << " that I can't recognize in this loop!";
                  //assert(0);

                }
              }
            }
          }
        }
        for(map<string,MonitorElement*>::iterator itHHH=itHH->second.begin(); itHHH!=itHH->second.end(); itHHH++){
//        LogInfo("AlcaBeamMonitorClient")
//          << itH->first << itHHH->first << " max-min=" << max-min << " delta=" << maxDelta-minDelta;
          if(itHHH->second != 0){
            if(itHHH->first == "Lumibased BeamSpotFit" 
            || itHHH->first == "Lumibased PrimaryVertex" 
            || itHHH->first == "Lumibased DataBase" 
            || itHHH->first == "Lumibased Scalers"){
              if((max == -bigNumber && min == bigNumber) || max-min == 0){
                itHHH->second->getTH1()->SetMinimum(itHHH->second->getTH1()->GetMinimum()-0.01);
                itHHH->second->getTH1()->SetMaximum(itHHH->second->getTH1()->GetMaximum()+0.01);
              }
              else{
                itHHH->second->getTH1()->SetMinimum(min-0.1*(max-min));
                itHHH->second->getTH1()->SetMaximum(max+0.1*(max-min));
              }
            }
            else if(itHHH->first == "Lumibased PrimaryVertex-DataBase fit" 
            || itHHH->first == "Lumibased PrimaryVertex-Scalers fit"
            || itHHH->first == "Lumibased Scalers-DataBase fit"
            || itHHH->first == "Lumibased PrimaryVertex-DataBase"
            || itHHH->first == "Lumibased PrimaryVertex-Scalers"
            ){
              if((maxDelta == -bigNumber && minDelta == bigNumber) || maxDelta-minDelta == 0){
                itHHH->second->getTH1()->SetMinimum(itHHH->second->getTH1()->GetMinimum()-0.01);
                itHHH->second->getTH1()->SetMaximum(itHHH->second->getTH1()->GetMaximum()+0.01);
              }
              else{
                itHHH->second->getTH1()->SetMinimum(minDelta-2*(maxDelta-minDelta));
                itHHH->second->getTH1()->SetMaximum(maxDelta+2*(maxDelta-minDelta));
              }
            }
//          else if(itHHH->first == "Lumibased PrimaryVertex-DataBase" 
//          || itHHH->first == "Lumibased PrimaryVertex-Scalers"
//          ){
//            if((maxDelta == -bigNumber && minDelta == bigNumber) || maxDelta-minDelta == 0){
//              itHHH->second->getTProfile()->SetMinimum(itHHH->second->getTProfile()->GetMinimum()-0.01);
//              itHHH->second->getTProfile()->SetMaximum(itHHH->second->getTProfile()->GetMaximum()+0.01);
//            }
//            else{
//              itHHH->second->getTProfile()->SetMinimum(minDelta-5*(maxDelta-minDelta));
//              itHHH->second->getTProfile()->SetMaximum(maxDelta+5*(maxDelta-minDelta));
 //             }
//          }
            else{
              LogInfo("AlcaBeamMonitorClient")
                << "The histosMap_ have a histogram named " << itHHH->first << " that I can't recognize in this loop!";
              //assert(0);

            }
          }
        }
      }
    }
  }


}

Member Data Documentation

DQMStore* AlcaBeamMonitorClient::dbe_ [private]

Definition at line 54 of file AlcaBeamMonitorClient.h.

Referenced by AlcaBeamMonitorClient(), beginRun(), endLuminosityBlock(), and endRun().

std::multimap<std::string,std::string> AlcaBeamMonitorClient::histoByCategoryNames_ [private]

Definition at line 59 of file AlcaBeamMonitorClient.h.

Referenced by AlcaBeamMonitorClient().

HistosContainer AlcaBeamMonitorClient::histosMap_ [private]

Definition at line 57 of file AlcaBeamMonitorClient.h.

Referenced by AlcaBeamMonitorClient(), beginRun(), and endRun().

std::string AlcaBeamMonitorClient::monitorName_ [private]

Definition at line 50 of file AlcaBeamMonitorClient.h.

Referenced by AlcaBeamMonitorClient(), beginRun(), endLuminosityBlock(), and endRun().

int AlcaBeamMonitorClient::numberOfValuesToSave_ [private]

Definition at line 53 of file AlcaBeamMonitorClient.h.

Referenced by AlcaBeamMonitorClient(), and endLuminosityBlock().

edm::ParameterSet AlcaBeamMonitorClient::parameters_ [private]

Definition at line 49 of file AlcaBeamMonitorClient.h.

PositionContainer AlcaBeamMonitorClient::positionsMap_ [private]

Definition at line 61 of file AlcaBeamMonitorClient.h.

Referenced by AlcaBeamMonitorClient(), and endRun().

std::map<edm::LuminosityBlockNumber_t,std::vector<double> > AlcaBeamMonitorClient::valuesMap_ [private]

Definition at line 60 of file AlcaBeamMonitorClient.h.

Referenced by endLuminosityBlock(), and endRun().

std::vector<std::string> AlcaBeamMonitorClient::varNamesV_ [private]

Definition at line 58 of file AlcaBeamMonitorClient.h.

Referenced by AlcaBeamMonitorClient().

Generated for CMSSW Reference Manual by doxygen 1.7.3