AlcaBeamMonitorClient Class Reference

#include <AlcaBeamMonitorClient.h>

Inheritance diagram for AlcaBeamMonitorClient:

Public Member Functions
	AlcaBeamMonitorClient (const edm::ParameterSet &)
	~AlcaBeamMonitorClient ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

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)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

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_

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

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

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 112 of file AlcaBeamMonitorClient.cc.

                                                                                 {
}

void AlcaBeamMonitorClient::beginJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file AlcaBeamMonitorClient.cc.

                                     {
}

void AlcaBeamMonitorClient::beginLuminosityBlock ( const edm::LuminosityBlock &  iLumi, const edm::EventSetup &  iSetup  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 108 of file AlcaBeamMonitorClient.cc.

                                                                                                       {
}

void AlcaBeamMonitorClient::beginRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 82 of file AlcaBeamMonitorClient.cc.

References dbe_, histosMap_, Association::map, 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  )

protectedvirtual

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  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 130 of file AlcaBeamMonitorClient.cc.

References newFWLiteAna::bin, DQMStore::book1D(), dbe_, first, histosMap_, edm::RunBase::id(), Association::map, 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.

Referenced by Mixins._TypedParameterizable::clone(), Types.PSet::clone(), and Mixins._TypedParameterizable::copy().

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().