Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 12:54:39

0001 // -*- C++ -*-
0002 //
0003 // Package:    ScalersRecover
0004 // Class:      ScalersRecover
0005 // 
0006 /**\class ScalersRecover ScalersRecover.cc DataFormats/ScalersRecover/src/ScalersRecover.cc
0007 
0008  Description: EDAnalyzer to fetch trigger scalers and convert it to SQL
0009 
0010  Implementation:
0011    This module may be used to recover L1 trigger scalers data from 
0012    CMS raw data files, and convert it to a series of SQL INSERT statements 
0013    that can be used to back-populate the corresponding L1 database 
0014    tables.   Should be performed on a run-by-run basis as necessary.
0015    First data file processed MUST contain data from the first lumi 
0016    sections.  In general, the files should be in lumi section order.
0017 
0018    We recommend running the job on MinBias RECO files.   If you 
0019    run with RAW files, you will have to include the ScalerRawToDigi 
0020    conversion module.
0021 
0022    The resulting SQL commands will be contained in a file 
0023    named scalers.sql
0024 */
0025 //
0026 // Original Author:  William Badgett
0027 //         Created:  Mon May 24 14:45:17 CEST 2010
0028 //
0029 //
0030 
0031 
0032 // system include files
0033 #include <memory>
0034 
0035 // user include files
0036 #include "FWCore/Framework/interface/Frameworkfwd.h"
0037 #include "FWCore/Framework/interface/EDAnalyzer.h"
0038 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0039 
0040 #include "FWCore/Framework/interface/Event.h"
0041 #include "FWCore/Framework/interface/MakerMacros.h"
0042 
0043 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0044 #include "DataFormats/Scalers/interface/Level1TriggerScalers.h"
0045 #include <ctime>
0046 #include "DataFormats/Scalers/interface/TimeSpec.h"
0047 //
0048 // class declaration
0049 //
0050 
0051 class ScalersRecover : public edm::EDAnalyzer 
0052 {
0053   public:
0054     explicit ScalersRecover(const edm::ParameterSet&);
0055     ~ScalersRecover();
0056 
0057 
0058   private:
0059     virtual void beginJob() ;
0060     virtual void analyze(const edm::Event&, const edm::EventSetup&);
0061     virtual void endJob() ;
0062 
0063     int lastLumiSection;
0064     FILE * sql;
0065 };
0066 
0067 // Constructor
0068 ScalersRecover::ScalersRecover(const edm::ParameterSet& iConfig)
0069 
0070 {
0071   sql = NULL;
0072   sql = fopen("scalers.sql","w");
0073 }
0074 
0075 // Destructor
0076 ScalersRecover::~ScalersRecover()
0077 {
0078   if ( sql != NULL ) { fclose(sql);}
0079 }
0080 
0081 // ------------ method called to for each event  ------------
0082 void ScalersRecover::analyze(const edm::Event& iEvent, 
0083                  const edm::EventSetup& iSetup)
0084 {
0085   using namespace edm;
0086   char heure [32];
0087   char sNanos [16];
0088   struct tm *hora;
0089 
0090   edm::Handle<Level1TriggerScalersCollection> data;
0091   bool ok = iEvent.getByLabel("scalersRawToDigi",data);
0092 
0093   if ( !ok ) 
0094   {
0095     LogError("ScalersRecover") << 
0096       "Could not find Level1TriggerScalersCollection";
0097     return;
0098   }
0099 
0100   if ( data->size() < 1 )
0101   {
0102     LogError("ScalersRecover") << 
0103       "Could not find Level1TriggerScalers element from Collection";
0104     return;
0105   }
0106 
0107   Level1TriggerScalersCollection::const_iterator triggerScalers = 
0108     data->begin();
0109 
0110 
0111   int lumiSection = triggerScalers->lumiSegmentNrLumiSeg();
0112   if ( ( ( lastLumiSection==-1 ) && ( lumiSection == 1 )) || 
0113        ( ( lastLumiSection>0 )   && ( lastLumiSection != lumiSection )))
0114   {
0115     timespec zeit = triggerScalers->collectionTimeLumiSeg();
0116     time_t seconds = zeit.tv_sec;
0117     long int nanos = zeit.tv_nsec;
0118 
0119     hora = gmtime(&seconds);
0120     strftime(heure,sizeof(heure),"%Y.%m.%d %H:%M:%S", hora);
0121     sprintf(sNanos,"%9.9d", (int)nanos);
0122 
0123     std::ostringstream insert;
0124     insert <<  
0125       "INSERT INTO LEVEL1_TRIGGER_CONDITIONS (RUNNUMBER,LUMISEGMENTNR,TIME,TIME_NS,TIME_STAMP" <<
0126       ",TRIGGERSPHYSICSGENERATEDFDL" << 
0127       ",TRIGGERSPHYSICSLOST" << 
0128       ",TRIGGERSPHYSICSLOSTBEAMACTIVE" << 
0129       ",TRIGGERSPHYSICSLOSTBEAMINACTI" << 
0130       ",L1ASPHYSICS" << 
0131       ",L1ASRANDOM" << 
0132       ",L1ASTEST" << 
0133       ",L1ASCALIBRATION" << 
0134       ",DEADTIME" << 
0135       ",DEADTIMEBEAMACTIVE" << 
0136       ",DEADTIMEBEAMACTIVETRIGGERRULE" << 
0137       ",DEADTIMEBEAMACTIVECALIBRATION" << 
0138       ",DEADTIMEBEAMACTIVEPRIVATEORBI" << 
0139       ",DEADTIMEBEAMACTIVEPARTITIONCO" << 
0140       ",DEADTIMEBEAMACTIVETIMESLOT" << 
0141       ") VALUES (" << iEvent.run() << 
0142       "," << lumiSection << 
0143       "," << zeit.tv_sec << 
0144       "," << nanos << 
0145       ",TO_TIMESTAMP('" << heure << "." << sNanos <<
0146       "','YYYY.MM.DD HH24:MI:SS.FF')" << 
0147       "," << triggerScalers->triggersPhysicsGeneratedFDL() <<
0148       "," << triggerScalers->triggersPhysicsLost() << 
0149       "," << triggerScalers->triggersPhysicsLostBeamActive() << 
0150       "," << triggerScalers->triggersPhysicsLostBeamInactive() << 
0151       "," << triggerScalers->l1AsPhysics() << 
0152       "," << triggerScalers->l1AsRandom() << 
0153       "," << triggerScalers->l1AsTest() << 
0154       "," << triggerScalers->l1AsCalibration() << 
0155       "," << triggerScalers->deadtime() << 
0156       "," << triggerScalers->deadtimeBeamActive() << 
0157       "," << triggerScalers->deadtimeBeamActiveTriggerRules() << 
0158       "," << triggerScalers->deadtimeBeamActiveCalibration() << 
0159       "," << triggerScalers->deadtimeBeamActivePrivateOrbit() << 
0160       "," << triggerScalers->deadtimeBeamActivePartitionController() << 
0161       "," << triggerScalers->deadtimeBeamActiveTimeSlot() << 
0162       ");" ;
0163     
0164     if ( sql != NULL ) { fprintf(sql,"%s\n", insert.str().c_str());}
0165     
0166     std::vector<unsigned int> algo = triggerScalers->gtAlgoCounts();
0167     int length = algo.size();
0168     for ( int i=0; i<length; i++)
0169     {
0170       std::ostringstream ainsert;
0171       ainsert << "INSERT INTO LEVEL1_TRIGGER_ALGO_CONDITIONS (RUNNUMBER,BIT,LUMISEGMENTNR,TIME,TIME_NS,TIME_STAMP,GTALGOCOUNTS) VALUES (" <<
0172     iEvent.run() << 
0173     "," << i << 
0174     "," << lumiSection << 
0175     "," << zeit.tv_sec << 
0176     "," << nanos << 
0177     ",TO_TIMESTAMP('" << heure << "." <<  sNanos << 
0178     "','YYYY.MM.DD HH24:MI:SS.FF')," 
0179           << algo[i] << ");";
0180       
0181       if ( sql != NULL ) { fprintf(sql,"%s\n", ainsert.str().c_str());}
0182     }
0183 
0184     std::vector<unsigned int> tech = triggerScalers->gtAlgoCounts();
0185     length = tech.size();
0186     for ( int i=0; i<length; i++)
0187     {
0188       std::ostringstream tinsert;
0189       tinsert << "INSERT INTO LEVEL1_TRIGGER_TECH_CONDITIONS (RUNNUMBER,BIT,LUMISEGMENTNR,TIME,TIME_NS,TIME_STAMP,GTTECHCOUNTS) VALUES (" <<
0190         iEvent.run() << 
0191         "," << i << 
0192         "," << lumiSection << 
0193         "," << zeit.tv_sec << 
0194         "," << nanos << 
0195         ",TO_TIMESTAMP('" << heure << "." <<  sNanos << 
0196         "','YYYY.MM.DD HH24:MI:SS.FF')," 
0197         << tech[i] << ");";
0198 
0199       if ( sql != NULL ) { fprintf(sql,"%s\n", tinsert.str().c_str());}
0200     }
0201     lastLumiSection = lumiSection;
0202   }
0203 }
0204 
0205 
0206 void ScalersRecover::beginJob()
0207 {
0208 }
0209 
0210 void ScalersRecover::endJob() 
0211 {
0212 }
0213 
0214 //define this as a plug-in
0215 DEFINE_FWK_MODULE(ScalersRecover);