EcalIntercalib method to read DBv2 xml files and create sqlite

cms-sw · Jul 3, 2017 · 78dbf41 · 78dbf41
1 parent fc51b5b
commit 78dbf41
Show file tree

Hide file tree

Showing 2 changed files with 77 additions and 18 deletions.
diff --git a/CondTools/Ecal/interface/EcalIntercalibHandler.h b/CondTools/Ecal/interface/EcalIntercalibHandler.h
@@ -61,6 +61,7 @@ namespace popcon
 			~EcalIntercalibHandler(); 
 
 			void getNewObjects();
+			void readXML(const std::string& filename, EcalFloatCondObjectContainer& record);
 
 			std::string id() const { return m_name;}
 			EcalCondDBInterface* econn;

diff --git a/CondTools/Ecal/src/EcalIntercalibHandler.cc b/CondTools/Ecal/src/EcalIntercalibHandler.cc
@@ -3,9 +3,10 @@
 #include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
 #include "CondTools/Ecal/interface/EcalFloatCondObjectContainerXMLTranslator.h"
 
-
 #include<iostream>
 
+const Int_t kEBChannels = 61200, kEEChannels = 14648;
+
 popcon::EcalIntercalibHandler::EcalIntercalibHandler(const edm::ParameterSet & ps)
   :    m_name(ps.getUntrackedParameter<std::string>("name","EcalIntercalibHandler")) {
 
@@ -175,28 +176,26 @@ void popcon::EcalIntercalibHandler::getNewObjects()
 
     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);
-
+      std::cout << "Generating popcon record for run " << irun << "..." << std::flush;
+      std::cout << "going to open file "<<file_ << std::flush;
 
-	    Time_t snc= (Time_t) irun ;
 
-	    popcon::PopConSourceHandler<EcalIntercalibConstants>::m_to_transfer.push_back(
-											  std::make_pair(payload,snc));
+      //      EcalCondHeader   header;
+      EcalIntercalibConstants * payload = new EcalIntercalibConstants;
+      //	    EcalFloatCondObjectContainerXMLTranslator::readXML(file_,header,*payload);   // DBv1
+      readXML(file_,*payload);   //DBv2
 
-	    ss << "Run=" << irun << "_Magnet_changed_"<<std::endl; 
-	    m_userTextLog = ss.str()+";";
+      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()+";";
+      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;
@@ -208,3 +207,62 @@ void popcon::EcalIntercalibHandler::getNewObjects()
   }
 	std::cout << "Ecal - > end of getNewObjects -----------\n";
 }
+
+void popcon::EcalIntercalibHandler::readXML(const std::string& file_,
+					    EcalFloatCondObjectContainer& record){
+  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;
+  }
+}