EcalIntercalibHandler.cc

Go to the documentation of this file.

#include "CondTools/Ecal/interface/EcalIntercalibHandler.h"
#include "FWCore/ParameterSet/interface/ParameterSetfwd.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "CondTools/Ecal/interface/EcalFloatCondObjectContainerXMLTranslator.h"


#include<iostream>

popcon::EcalIntercalibHandler::EcalIntercalibHandler(const edm::ParameterSet & ps)
  :    m_name(ps.getUntrackedParameter<std::string>("name","EcalIntercalibHandler")) {

    std::cout << "EcalIntercalib Source handler constructor\n" << std::endl;
        m_firstRun=static_cast<unsigned int>(atoi( ps.getParameter<std::string>("firstRun").c_str()));
        m_lastRun=static_cast<unsigned int>(atoi( ps.getParameter<std::string>("lastRun").c_str()));
        m_sid= ps.getParameter<std::string>("OnlineDBSID");
        m_user= ps.getParameter<std::string>("OnlineDBUser");
        m_pass= ps.getParameter<std::string>("OnlineDBPassword");
        m_locationsource= ps.getParameter<std::string>("LocationSource");
        m_location=ps.getParameter<std::string>("Location");
        m_gentag=ps.getParameter<std::string>("GenTag");
        m_file_lowfield= ps.getParameter<std::string>("FileLowField");
        m_file_highfield= ps.getParameter<std::string>("FileHighField");



        std::cout << m_sid<<"/"<<m_user<<"/"<<m_location<<"/"<<m_gentag   << std::endl;


}

popcon::EcalIntercalibHandler::~EcalIntercalibHandler()
{
}


void popcon::EcalIntercalibHandler::getNewObjects()
{
  
  std::cout << "------- Ecal - > getNewObjects\n";
  
  std::ostringstream ss; 
  ss<<"ECAL ";
  
  unsigned int max_since=0;
  max_since=static_cast<unsigned int>(tagInfo().lastInterval.first);
  std::cout << "max_since : "  << max_since << std::endl;
  Ref ped_db = lastPayload();
  
  // we parse the last record in the DB and check if it is low or high field 
  
  std::cout << "retrieved last payload "  << std::endl;
  
  
  EcalIntercalibConstant the_cal = 0.;
  
  int iX=50;
  int iY=5;
  int iZ=-1;
  float the_value_high=0.75585;
  
  if (EEDetId::validDetId(iX,iY,iZ))
    {
      EEDetId eedetidpos(iX,iY,iZ);
      
      EcalIntercalibConstants::const_iterator it =ped_db->find(eedetidpos.rawId());
      
      
      the_cal = (*it);
      
    }
  
  bool magnet_high=true; 
  if(the_cal!= the_value_high) magnet_high=false; 
  
  
  // here we connect to the online DB to check the value of the magnetic field 
  
  std::cout << "Connecting to ONLINE DB ... " << std::endl;
  econn = new EcalCondDBInterface( m_sid, m_user, m_pass );
  std::cout << "Connection done" << std::endl;
  
  if (!econn)
    {
      std::cout << " Problem with OMDS: connection parameters " <<m_sid <<"/"<<m_user<<"/"<<m_pass<<std::endl;
      throw cms::Exception("OMDS not available");
    } 
  
  
  std::cout << "Retrieving last run from ONLINE DB ... " << std::endl;
  std::map<EcalLogicID, RunDat> rundat;
  RunIOV rp ;
  run_t runmax=10000000;
  std::string location_p5="P5_Co";
  econn->fetchValidDataSet(&rundat , &rp, location_p5 ,runmax);
  
  unsigned long long  irun=(unsigned long long) rp.getRunNumber();
  
  std::cout<< "retrieved run number "<< irun <<std::endl;  
  
  if(irun>max_since) {
    
    
    // retrieve from last value data record 
    // always call this method at first run
    
    std::map<EcalLogicID, RunDCSMagnetDat> dataset;
    
    econn->fetchDataSet(&dataset, &rp);
    
    if (!dataset.size()) {
      throw(std::runtime_error("Zero rows read back"));
    } else {
      std::cout<< "retrieved magnet current"<<std::endl;  
    }
    
    
    float mag_cur=0;
    
    std::map<  EcalLogicID, RunDCSMagnetDat >::iterator it;
    for (it=dataset.begin(); it!=dataset.end(); ++it){
      
      RunDCSMagnetDat  a_mag = (*it).second;
      mag_cur= a_mag.getMagnetCurrent();
      
    }
    
    
    std::string file_=m_file_highfield;
    bool something_to_transfer=false;
    if(mag_cur>5000. && magnet_high ) {
      
      std::cout << " the magnet is ON and the constants are for magnet ON " << std::endl; 
      
    } else if(mag_cur>5000. && !magnet_high ) {
      something_to_transfer=true;
      std::cout << " the magnet is ON and the constants are for magnet OFF " << std::endl; 
      std::cout << " I transfer the ON constants "<< std::endl; 
      file_=m_file_highfield;
      
    } else if(mag_cur<3000. && magnet_high ) {
      something_to_transfer=true;
      std::cout << " the magnet is OFF and the constants are for magnet ON "<< std::endl;
      std::cout << " I transfer the OFF constants "<< std::endl;
      file_=m_file_lowfield;
      
    } else if( mag_cur<3000. && !magnet_high ){
      
      std::cout << " the magnet is OFF and the constants are for magnet OFF "<< std::endl;
      file_=m_file_lowfield;
      
    } else {
      
      std::cout << " the magnet is in a strange situation I do nothing ... just be patient "<< std::endl;
      
    }
    
    
    if(something_to_transfer){
      
      std::cout << "Generating popcon record for run " << irun << "..." << std::flush;
      std::cout << "going to open file "<<file_ << std::flush;
      
      
      EcalCondHeader   header;
      EcalIntercalibConstants * payload = new EcalIntercalibConstants;
      EcalFloatCondObjectContainerXMLTranslator::readXML(file_,header,*payload);
      
      
      Time_t snc= (Time_t) irun ;
      
      popcon::PopConSourceHandler<EcalIntercalibConstants>::m_to_transfer.push_back(
                                            std::make_pair(payload,snc));
      
      ss << "Run=" << irun << "_Magnet_changed_"<<std::endl; 
      m_userTextLog = ss.str()+";";
      
      
    } else {
      std::cout << "Run " << irun << " nothing sent to the DB"<< std::endl;
      
      ss<< "Run=" << irun << "_Magnet_NOT_changed_"<<std::endl; 
      m_userTextLog = ss.str()+";";
    }
    
    
    delete econn;
  } else {
    std::cout << "Run " << irun << " nothing sent to the DB"<< std::endl;
    ss<< "Run=" << irun << "_no_new_runs_"<<std::endl; 
    m_userTextLog = ss.str()+";";
    
    
  }
  
  
  
  std::cout << "Ecal - > end of getNewObjects -----------\n";
  
}