CastorDigiMonitor.cc

Go to the documentation of this file.

#include "DQM/CastorMonitor/interface/CastorDigiMonitor.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

//****************************************************//
//********** CastorDigiMonitor: ******************//
//********** Author: Dmytro Volyanskyy   *************//
//********** Date  : 29.08.2008 (first version) ******// 
//****************************************************//

//==================================================================//
//======================= Constructor ==============================//
//==================================================================//
CastorDigiMonitor::CastorDigiMonitor() { doPerChannel_ = false;  }


//==================================================================//
//======================= Destructor ===============================//
//==================================================================//
CastorDigiMonitor::~CastorDigiMonitor() {
}


//==================================================================//
//=========================== reset  ===============================//
//==================================================================//
void CastorDigiMonitor::reset(){
}

//==================================================================//
//=========================== setup  ===============================//
//==================================================================//
void CastorDigiMonitor::setup(const edm::ParameterSet& ps, DQMStore* dbe){
  CastorBaseMonitor::setup(ps,dbe);
  baseFolder_ = rootFolder_+"CastorDigiMonitor";
   
  if(fVerbosity>0) std::cout << "CastorDigiMonitor::setup (start)" << std::endl;

  doPerChannel_ = ps.getUntrackedParameter<bool>("DigiPerChannel", false);
  doFCpeds_ = ps.getUntrackedParameter<bool>("DigiInFC", true);

  ievt_=0;

  if ( m_dbe !=NULL ) {
    m_dbe->setCurrentFolder(baseFolder_);
    meEVT_ = m_dbe->bookInt("Digi Task Event Number");
    meEVT_->Fill(ievt_);
        
    m_dbe->setCurrentFolder(baseFolder_);
  
    std::string type = "Castor Digis ADC counts";
    castHists.ALLPEDS =  m_dbe->book1D(type,type,130,0,130);
    
    //type = "Castor Pedestal Mean Reference Values - from CondDB";
    //castHists.PEDESTAL_REFS = m_dbe->book1D(type,type,50,0,50);
    //type = "Castor Pedestal RMS Reference Values - from CondDB";
    //castHists.WIDTH_REFS = m_dbe->book1D(type,type,20,0,10); 
 
}

  else{ 
   if(fVerbosity>0) std::cout << "CastorDigiMonitor::setup - NO DQMStore service" << std::endl; 
  }

 
  outputFile_ = ps.getUntrackedParameter<std::string>("PedestalFile", "");
  if ( outputFile_.size() != 0 ) { if(fVerbosity>0) std::cout << "Castor Pedestal Calibrations will be saved to " << outputFile_.c_str() << std::endl;}


 if(fVerbosity>0) std::cout << "CastorDigiMonitor::setup (end)" << std::endl;

  return;
}




//==================================================================//
//=========================== processEvent  ========================//
//==================================================================//
void CastorDigiMonitor::processEvent(const CastorDigiCollection& castorDigis, const CastorDbService& cond)
{
  
 if(fVerbosity>0) std::cout << "==>CastorDigiMonitor::processEvent !!!"<< std::endl;

 
  meEVT_->Fill(ievt_);
  
  //if(!shape_) shape_ = cond.getCastorShape(); // this one is generic


  if(!m_dbe) { 
    if(fVerbosity>0) std::cout<<"CastorDigiMonitor::processEvent DQMStore is not instantiated!!!"<<std::endl;  
    return; 
  }


  CaloSamples tool;  
 
  if(castorDigis.size()>0) {

   for (CastorDigiCollection::const_iterator j=castorDigis.begin(); j!=castorDigis.end(); j++){
      const CastorDataFrame digi = (const CastorDataFrame)(*j); 
 

       detID_.clear(); capID_.clear(); pedVals_.clear();



      // const CastorCalibrations& calibrations = cond.getCastorCalibrations(digi.id().rawId());
      // const CastorPedestal* ped              = cond.getPedestal(digi.id()); 
      // const CastorPedestalWidth* pedw        = cond.getPedestalWidth(digi.id());
       // const CastorPedestal* ped = cond.getPedestal(digi.id()); 
       // const CastorPedestalWidth* pedw = cond.getPedestalWidth(digi.id());
      
       /*
      if(doFCpeds_){
    channelCoder_ = cond.getCastorCoder(digi.id());
    CastorCoderDb coderDB(*channelCoder_, *shape_);
    coderDB.adc2fC(digi,tool);
      }
     
     
       for(int capID=0; capID<4; capID++){
        float pedvalue=0;    
        if(ped) pedvalue=ped->getValue(capID);
            castHists.PEDESTAL_REFS->Fill(pedvalue);
            PEDESTAL_REFS->Fill(pedvalue);
           float width=0;
       if(pedw) width = pedw->getWidth(capID);
           castHists.WIDTH_REFS->Fill(width);
           WIDTH_REFS->Fill(width);
    }
     */
     
 
       
       //      if(ievt_ %1000 == 0 )           // PK: skip limited number of events
      //   { 
      for (int i=0; i<digi.size(); i++) {
    if(doFCpeds_) pedVals_.push_back(tool[i]); // default is FALSE
    else pedVals_.push_back(digi.sample(i).adc());
    detID_.push_back(digi.id());
    capID_.push_back(digi.sample(i).capid());
    castHists.ALLPEDS->Fill(pedVals_[i]);
      }
      
      //      }      

      //      if( ievt_%100 == 0 && doPerChannel_) perChanHists(detID_,capID_,pedVals_,castHists.PEDVALS, baseFolder_);
      if( doPerChannel_) perChanHists(detID_,capID_,pedVals_,castHists.PEDVALS, baseFolder_); // PK: no special event selection done

    }
  } 
   else {
    if(fVerbosity>0) std::cout << "CastorPSMonitor::processEvent NO Castor Digis !!!" << std::endl;
  }

 if (showTiming) { 
      cpu_timer.stop(); std::cout << " TIMER::CastorDigi -> " << cpu_timer.cpuTime() << std::endl; 
      cpu_timer.reset(); cpu_timer.start();  
    }

  ievt_++;
  return;
}

void CastorDigiMonitor::done(){
  return;
}


//==================================================================//
//======================= perChanHists  ============================//
//==================================================================//
void CastorDigiMonitor::perChanHists( std::vector<HcalCastorDetId> detID, std::vector<int> capID, std::vector<float> peds,
                          std::map<HcalCastorDetId, std::map<int, MonitorElement*> > &toolP,  
                                          std::string baseFolder) 
 {
  
  if(m_dbe) m_dbe->setCurrentFolder(baseFolder);

  for(unsigned int d=0; d<detID.size(); d++){
    HcalCastorDetId detid = detID[d];
    int capid = capID[d];
    float pedVal = peds[d];
    bool gotit=false;
    if(REG[detid]) gotit=true;
    
    if(gotit){
      std::map<int, MonitorElement*> _mei = toolP[detid];
      if(_mei[capid]==NULL){
    if(fVerbosity>0) std::cout<<"CastorDigiMonitor::perChanHists  This histo is NULL!!??"<< std::endl;
      }
      else _mei[capid]->Fill(pedVal);
      
    }
    else{
      if(m_dbe){
    std::map<int,MonitorElement*> insertP; //-- Pedestal values in ADC
    
    for(int i=0; i<4; i++){
      char name[1024];
      sprintf(name,"Castor Digi Value (ADC) zside=%d module=%d sector=%d CAPID=%d",
          detid.zside(),detid.module(),detid.sector(),i);      
      insertP[i] =  m_dbe->book1D(name,name,10,-0.5,9.5);
      
    }
    
    insertP[capid]->Fill(pedVal);
    toolP[detid] = insertP;
      }
      REG[detid] = true;
    }
  }
}