#include <EcalIntercalibHandler.h>

Inheritance diagram for popcon::EcalIntercalibHandler:

Public Member Functions
	EcalIntercalibHandler (edm::ParameterSet const &)

void	getNewObjects ()

std::string	id () const

void	readXML (const std::string &filename, EcalFloatCondObjectContainer &record)

	~EcalIntercalibHandler ()

Public Member Functions inherited from popcon::PopConSourceHandler< EcalIntercalibConstants >
void	convertFromOld ()

Summary *	dummySummary (typename OldContainer::value_type const &) const

void	initialize (const cond::persistency::Session &dbSession, cond::TagInfo_t const &tagInfo, cond::LogDBEntry_t const &logDBEntry)

Ref	lastPayload () const

cond::LogDBEntry_t const &	logDBEntry () const

std::pair< Container const *, std::string const >	operator() (const cond::persistency::Session &session, cond::TagInfo_t const &tagInfo, cond::LogDBEntry_t const &logDBEntry) const

	PopConSourceHandler ()

Container const &	returnData ()

void	sort ()

cond::TagInfo_t const &	tagInfo () const

std::string const &	userTextLog () const

virtual	~PopConSourceHandler ()

Public Attributes
EcalCondDBInterface *	econn

Private Attributes
std::string	m_file_highfield

std::string	m_file_lowfield

unsigned int	m_firstRun

std::string	m_gentag

unsigned int	m_lastRun

std::string	m_location

std::string	m_locationsource

std::string	m_name

std::string	m_pass

std::string	m_sid

std::string	m_user

double	m_value_highfield

const EcalIntercalibConstants *	myintercalib

Additional Inherited Members
Public Types inherited from popcon::PopConSourceHandler< EcalIntercalibConstants >
typedef std::vector< Triplet >	Container

typedef std::vector< std::pair< EcalIntercalibConstants *, cond::Time_t > >	OldContainer

typedef PopConSourceHandler< EcalIntercalibConstants >	self

typedef cond::Summary	Summary

typedef cond::Time_t	Time_t

typedef EcalIntercalibConstants	value_type

Protected Member Functions inherited from popcon::PopConSourceHandler< EcalIntercalibConstants >
int	add (value_type payload, Summary summary, Time_t time)

cond::persistency::Session &	dbSession () const

Protected Attributes inherited from popcon::PopConSourceHandler< EcalIntercalibConstants >
OldContainer	m_to_transfer

std::string	m_userTextLog

Detailed Description

Definition at line 56 of file EcalIntercalibHandler.h.

Constructor & Destructor Documentation

popcon::EcalIntercalibHandler::EcalIntercalibHandler ( edm::ParameterSet const & )

Definition at line 10 of file EcalIntercalibHandler.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_file_highfield, m_file_lowfield, m_firstRun, m_gentag, m_lastRun, m_location, m_locationsource, m_pass, m_sid, m_user, m_value_highfield, and AlCaHLTBitMon_QueryRunRegistry::string.

   :    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");
     m_value_highfield= ps.getUntrackedParameter< double >("Value_Bon");
     //  m_value_highfield= 0.75585;
 
 
 
         std::cout << m_sid<<"/"<<m_user<<"/"<<m_location<<"/"<<m_gentag   << std::endl;
 
 }

popcon::EcalIntercalibHandler::~EcalIntercalibHandler ( )

Definition at line 33 of file EcalIntercalibHandler.cc.

34 {

35 }

Member Function Documentation

void popcon::EcalIntercalibHandler::getNewObjects ( )

virtual

Implements popcon::PopConSourceHandler< EcalIntercalibConstants >.

Definition at line 38 of file EcalIntercalibHandler.cc.

References gather_cfg::cout, createfilelist::dataset, econn, Exception, EcalCondDBInterface::fetchDataSet(), EcalCondDBInterface::fetchValidDataSet(), objects.autophobj::float, RunDCSMagnetDat::getMagnetCurrent(), RunIOV::getRunNumber(), cond::TagInfo_t::lastInterval, popcon::PopConSourceHandler< EcalIntercalibConstants >::lastPayload(), m_file_highfield, m_file_lowfield, m_pass, m_sid, m_user, popcon::PopConSourceHandler< EcalIntercalibConstants >::m_userTextLog, m_value_highfield, jetCorrFactors_cfi::payload, DetId::rawId(), readXML(), lumiCalc2::runmax, cond::TagInfo_t::size, findQualityFiles::size, AlCaHLTBitMon_QueryRunRegistry::string, popcon::PopConSourceHandler< EcalIntercalibConstants >::tagInfo(), and EEDetId::validDetId().

 {
 
   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;
   bool something_to_transfer = false;
   bool magnet_high = true; 
   if(tagInfo().size) {
     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. ; // relies on it being a float.
                                           // instead should perhaps
                                           // protect the next if when
                                               // the EEDetId isn't valid?
 
     int iX=50;
     int iY=5;
     int iZ=-1;
 
 
     float the_value_high=(float)m_value_highfield;
     std::cout << "The value for high field at EE x/y/z= 50/5/-1 is " << the_value_high << std::endl;
 
     if (EEDetId::validDetId(iX,iY,iZ))
       {
         EEDetId eedetidpos(iX,iY,iZ);
         
         EcalIntercalibConstants::const_iterator it =ped_db->find(eedetidpos.rawId());
         
 
         the_cal = (*it);
 
       }
     
     if(the_cal!= the_value_high) magnet_high=false;
   }  // check if there is already a payload
   else something_to_transfer = true;
 
   // 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;
     if(tagInfo().size) {
 
       if(mag_cur>7000. && magnet_high ) {
         
         std::cout << " the magnet is ON and the constants are for magnet ON " << std::endl; 
         
       } else if(mag_cur>7000. && !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<6000. && 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<6000. && !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;
         
       }
     }
     else {
       if(mag_cur>7000.)
     std::cout <<" first payload, the magnet is ON " << std::endl;
       else if( mag_cur<6000.) {
     std::cout <<" first payload, the magnet is 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);   // DBv1
       readXML(file_,*payload);   //DBv2
         
       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;
   }  // irun > max_since
   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";
 }

std::string popcon::EcalIntercalibHandler::id ( ) const

inlinevirtual

Implements popcon::PopConSourceHandler< EcalIntercalibConstants >.

Definition at line 66 of file EcalIntercalibHandler.h.

66 { return m_name;}

popcon::EcalIntercalibHandler::m_name

std::string m_name

Definition: EcalIntercalibHandler.h:83

void popcon::EcalIntercalibHandler::readXML	(	const std::string &	filename,
		EcalFloatCondObjectContainer &	record
	)

Definition at line 211 of file EcalIntercalibHandler.cc.

References gather_cfg::cout, cmsRelvalreport::exit, runEdmFileComparison::found, mps_fire::i, kEBChannels, kEEChannels, AlCaHLTBitMon_QueryRunRegistry::string, EBDetId::unhashIndex(), EEDetId::unhashIndex(), and heppy_batch::val.

Referenced by getNewObjects().

                                                              {
   std::string dummyLine, bid;
   std::ifstream fxml;
   fxml.open(file_);
   if(!fxml.is_open()) {
     std::cout << "ERROR : cannot open file " << file_ << std::endl;
     exit (1);
   }
   // header
   for( int i=0; i< 6; i++) {
     getline(fxml, dummyLine);   // skip first lines
     //  std::cout << dummyLine << std::endl;
   }
   fxml >> bid;
   //  std::cout << bid << std::endl;
   std::string stt = bid.substr(7,5);
   std::istringstream iEB(stt);
   int nEB;
   iEB >> nEB;
   if(nEB != kEBChannels) {
     std::cout << " strange number of EB channels " << nEB << std::endl;
     exit(-1);
   }
   fxml >> bid;   // <item_version>0</item_version>
   for (int iChannel = 0; iChannel < kEBChannels; iChannel++) {
     EBDetId myEBDetId = EBDetId::unhashIndex(iChannel);
     fxml >> bid;
     std::size_t found = bid.find("</");
     stt = bid.substr(6, found - 6);
     float val = std::stof(stt);
     record[myEBDetId] = val;
   }
   for( int i=0; i< 5; i++) {
     getline(fxml, dummyLine);   // skip first lines
     //  std::cout << dummyLine << std::endl;
   }
   fxml >> bid;
   //    cout << bid << endl;
   stt = bid.substr(7,5);
   std::istringstream iEE(stt);
   int nEE;
   iEE >> nEE;
   if(nEE != kEEChannels) {
     std::cout << " strange number of EE channels " << nEE << std::endl;
     exit(-1);
   }
   fxml >> bid;   // <item_version>0</item_version>
   // now endcaps
   for (int iChannel = 0; iChannel < kEEChannels; iChannel++) {
     EEDetId myEEDetId = EEDetId::unhashIndex(iChannel);
     fxml >> bid;
     std::size_t found = bid.find("</");
     stt = bid.substr(6, found - 6);
     float val = std::stof(stt);
     record[myEEDetId] = val;
   }
 }