L1TGT Class Reference

#include <L1TGT.h>

Inheritance diagram for L1TGT:

Public Member Functions
	L1TGT (const edm::ParameterSet &ps)
virtual	~L1TGT ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
void	beginJob ()
void	endJob (void)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Types
enum	activeDAQ {   FDL =0, PSB9, PSB13, PSB14,   PSB15, PSB19, PSB20, PSB21,   GMT }
enum	activeEVM { TCS, FDLEVM }

Private Member Functions
bool	isActive (int word, int bit)

Private Attributes
MonitorElement *	algo_bits
MonitorElement *	algo_bits_corr
MonitorElement *	algo_bits_lumi
MonitorElement *	algo_tt_bits_corr
MonitorElement *	BST_beamMode
MonitorElement *	BST_beamMomentum
MonitorElement *	BST_intensityBeam1
MonitorElement *	BST_intensityBeam2
MonitorElement *	BST_lhcFillNumber
MonitorElement *	BST_MasterStatus
MonitorElement *	BST_turnCountNumber
DQMStore *	dbe
MonitorElement *	dbx_module
MonitorElement *	event_lumi
MonitorElement *	event_number
MonitorElement *	event_type
MonitorElement *	evnum_trignum_lumi
MonitorElement *	gpsfreq
MonitorElement *	gpsfreqlum
MonitorElement *	gpsfreqwide
edm::InputTag	gtEvmSource_
MonitorElement *	gtfe_bx
edm::InputTag	gtSource_
ofstream	logFile_
MonitorElement *	m_monL1PrescaleFactorSet
bool	monitorDaemon_
int	nev_
MonitorElement *	orbit_lumi
std::string	outputFile_
boost::uint64_t	preGps_
boost::uint64_t	preOrb_
MonitorElement *	setupversion_lumi
MonitorElement *	trigger_lumi
MonitorElement *	trigger_number
MonitorElement *	tt_bits
MonitorElement *	tt_bits_corr
MonitorElement *	tt_bits_lumi
bool	verbose_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 39 of file L1TGT.h.

Member Enumeration Documentation

enum L1TGT::activeDAQ

private

Enumerator
FDL
PSB9
PSB13
PSB14
PSB15
PSB19
PSB20
PSB21
GMT

Definition at line 72 of file L1TGT.h.

{ FDL=0, PSB9, PSB13, PSB14, PSB15, PSB19, PSB20, PSB21, GMT };

enum L1TGT::activeEVM

private

Enumerator
TCS
FDLEVM

Definition at line 76 of file L1TGT.h.

{ TCS, FDLEVM };

Constructor & Destructor Documentation

L1TGT::L1TGT

(

const edm::ParameterSet &

ps

)

Definition at line 20 of file L1TGT.cc.

References gather_cfg::cout, dbe, edm::ParameterSet::getUntrackedParameter(), NULL, cmsCodeRules.cppFunctionSkipper::operator, outputFile_, DQMStore::setCurrentFolder(), DQMStore::setVerbose(), and verbose_.

                                   :
    gtSource_(ps.getParameter<InputTag> ("gtSource")),
    gtEvmSource_(ps.getParameter<InputTag> ("gtEvmSource")) {

    // verbosity switch
    verbose_ = ps.getUntrackedParameter<bool> ("verbose", false);

    if (verbose_)
        cout << "L1TGT: constructor...." << endl;

    dbe = NULL;
    if (ps.getUntrackedParameter<bool> ("DQMStore", false)) {
        dbe = Service<DQMStore> ().operator->();
        dbe->setVerbose(0);
    }

    outputFile_ = ps.getUntrackedParameter<string> ("outputFile", "");
    if (outputFile_.size() != 0) {
        cout << "L1T Monitoring histograms will be saved to "
                << outputFile_.c_str() << endl;
    }

    bool disable = ps.getUntrackedParameter<bool> ("disableROOToutput", false);
    if (disable) {
        outputFile_ = "";
    }

    if (dbe != NULL) {
        dbe->setCurrentFolder("L1T/L1TGT");
    }

}

L1TGT::~L1TGT ( )

virtual

Definition at line 53 of file L1TGT.cc.

              {

    // empty
}

Member Function Documentation

void L1TGT::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

protectedvirtual

get Global Trigger algo and technical trigger bit statistics

Implements edm::EDAnalyzer.

Definition at line 291 of file L1TGT.cc.

References L1GtfeWord::activeBoards(), algo_bits, algo_bits_corr, algo_bits_lumi, L1GtfeExtWord::beamMode(), L1GtfeExtWord::beamMomentum(), BST_beamMode, BST_beamMomentum, BST_intensityBeam1, BST_intensityBeam2, BST_lhcFillNumber, BST_MasterStatus, BST_turnCountNumber, L1GtfeExtWord::bstMasterStatus(), L1TcsWord::bxNr(), L1GtPsbWord::bxNr(), L1GtFdlWord::bxNr(), L1GtfeWord::bxNr(), gather_cfg::cout, dbx_module, event_lumi, event_number, event_type, L1TcsWord::eventNr(), evnum_trignum_lumi, FDL, MonitorElement::Fill(), edm::Event::getByLabel(), GMT, gpsfreq, gpsfreqlum, gpsfreqwide, L1GtfeExtWord::gpsTime(), gtEvmSource_, gtfe_bx, L1GtFdlWord::gtPrescaleFactorIndexAlgo(), gtSource_, i, isActive(), edm::HandleBase::isValid(), j, edm::InputTag::label(), L1GtfeExtWord::lhcFillNumber(), L1GtPsbWord::localBxNr(), L1GtFdlWord::localBxNr(), edm::EventBase::luminosityBlock(), m_monL1PrescaleFactorSet, funct::master(), nev_, orbit_lumi, L1TcsWord::orbitNr(), L1TcsWord::partTrigNr(), preGps_, preOrb_, PSB9, L1GtfeWord::setupVersion(), setupversion_lumi, TCS, L1GtfeExtWord::totalIntensityBeam1(), L1GtfeExtWord::totalIntensityBeam2(), trigger_lumi, trigger_number, L1TcsWord::triggerType(), tt_bits, tt_bits_corr, tt_bits_lumi, L1GtfeExtWord::turnCountNumber(), and verbose_.

                                                       {

    nev_++;

    if (verbose_) {
        cout << "L1TGT: analyze...." << endl;
    }

    // initialize bx's to invalid value
    int tcsBx = -1;
    int gtfeEvmBx = -1;

    // open EVM readout record if available
    Handle<L1GlobalTriggerEvmReadoutRecord> gtEvmReadoutRecord;
    e.getByLabel(gtEvmSource_, gtEvmReadoutRecord);

    if (!gtEvmReadoutRecord.isValid()) {
        edm::LogInfo("DataNotFound")
                << "can't find L1GlobalTriggerEvmReadoutRecord with label "
                << gtSource_.label();
    } else {

        // get all info from the EVM record if available and fill the histograms

        const L1GtfeWord& gtfeEvmWord = gtEvmReadoutRecord->gtfeWord();
        const L1GtfeExtWord& gtfeEvmExtWord = gtEvmReadoutRecord->gtfeWord();

        gtfeEvmBx = gtfeEvmWord.bxNr();
        int gtfeEvmActiveBoards = gtfeEvmWord.activeBoards();

        if (isActive(gtfeEvmActiveBoards, TCS)) { // if TCS present in the record

            const L1TcsWord& tcsWord = gtEvmReadoutRecord->tcsWord();
            tcsBx = tcsWord.bxNr();

            event_type->Fill(tcsWord.triggerType());
            orbit_lumi->Fill(e.luminosityBlock(), tcsWord.orbitNr());

            trigger_number->Fill(tcsWord.partTrigNr());
            event_number->Fill(tcsWord.eventNr());

            trigger_lumi->Fill(e.luminosityBlock(), tcsWord.partTrigNr());
            event_lumi->Fill(e.luminosityBlock(), tcsWord.eventNr());
            evnum_trignum_lumi->Fill(e.luminosityBlock(),
                    double(tcsWord.eventNr()) / double(tcsWord.partTrigNr()));

            boost::uint16_t master = gtfeEvmExtWord.bstMasterStatus();
            boost::uint32_t turnCount = gtfeEvmExtWord.turnCountNumber();
            boost::uint32_t lhcFill = gtfeEvmExtWord.lhcFillNumber();
            boost::uint16_t beam = gtfeEvmExtWord.beamMode();
            boost::uint16_t momentum = gtfeEvmExtWord.beamMomentum();
            boost::uint32_t intensity1 = gtfeEvmExtWord.totalIntensityBeam1();
            boost::uint32_t intensity2 = gtfeEvmExtWord.totalIntensityBeam2();

            BST_MasterStatus->Fill(e.luminosityBlock(), (float) (master));
            BST_turnCountNumber->Fill(e.luminosityBlock(), (float) (turnCount));
            BST_lhcFillNumber->Fill((float) (lhcFill % 1000));
            BST_beamMode->Fill(e.luminosityBlock(), (float) (beam));

            BST_beamMomentum->Fill(e.luminosityBlock(), (float) (momentum));
            BST_intensityBeam1->Fill(e.luminosityBlock(), (float) (intensity1));
            BST_intensityBeam2->Fill(e.luminosityBlock(), (float) (intensity2));

            if (verbose_) {
                cout << " check mode = " << beam << "    momentum " << momentum
                        << " int2 " << intensity2 << endl;
            }

            boost::uint64_t orb = tcsWord.orbitNr();
            boost::uint64_t gpsr = gtfeEvmExtWord.gpsTime();
            boost::uint64_t gpshi = (gpsr >> 32) & 0xffffffff;
            boost::uint64_t gpslo = gpsr & 0xffffffff;
            boost::uint64_t gps = gpshi * 1000000 + gpslo;
            //  cout << "  gpsr = " << hex << gpsr << " hi=" << gpshi << " lo=" << gpslo << " gps=" << gps << endl;

            Long64_t delorb = orb - preOrb_;
            Long64_t delgps = gps - preGps_;
            Double_t freq = -1.;

            if (delgps > 0) {
                freq = ((Double_t)(delorb)) * 3564. / ((Double_t)(delgps));
            }

            if (delorb > 0) {
                gpsfreq->Fill(freq);
                gpsfreqwide->Fill(freq);
                gpsfreqlum->Fill(e.luminosityBlock(), freq);
                if (verbose_) {
                    if (freq > 200.) {
                        cout << " preOrb_ = " << preOrb_ << " orb=" << orb
                                << " delorb=" << delorb << hex << " preGps_="
                                << preGps_ << " gps=" << gps << dec
                                << " delgps=" << delgps << " freq=" << freq
                                << endl;

                    }
                }
            }

            preGps_ = gps;
            preOrb_ = orb;

        }
    }

    // open GT DAQ readout record - exit if failed
    Handle<L1GlobalTriggerReadoutRecord> gtReadoutRecord;
    e.getByLabel(gtSource_, gtReadoutRecord);

    if (!gtReadoutRecord.isValid()) {
        edm::LogInfo("DataNotFound")
                << "can't find L1GlobalTriggerReadoutRecord with label "
                << gtSource_.label();
        return;
    }

    // initialize bx's to invalid value
    int gtfeBx = -1;
    int fdlBx[2] = { -1, -1 };
    int psbBx[2][7] = { { -1, -1, -1, -1, -1, -1, -1 }, { -1, -1, -1, -1, -1,
            -1, -1 } };
    int gmtBx = -1;

    // get info from GTFE DAQ record
    const L1GtfeWord& gtfeWord = gtReadoutRecord->gtfeWord();
    gtfeBx = gtfeWord.bxNr();
    gtfe_bx->Fill(gtfeBx);
    setupversion_lumi->Fill(e.luminosityBlock(), gtfeWord.setupVersion());
    int gtfeActiveBoards = gtfeWord.activeBoards();

    // look for GMT readout collection from the same source if GMT active
    if (isActive(gtfeActiveBoards, GMT)) {
        edm::Handle<L1MuGMTReadoutCollection> gmtReadoutCollection;
        e.getByLabel(gtSource_, gmtReadoutCollection);

        if (gmtReadoutCollection.isValid()) {
            gmtBx = gmtReadoutCollection->getRecord().getBxNr();
        }
    }

    // get info from FDL if active (including decision word)
    if (isActive(gtfeActiveBoards, FDL)) {
        const L1GtFdlWord& fdlWord = gtReadoutRecord->gtFdlWord();
        fdlBx[0] = fdlWord.bxNr();
        fdlBx[1] = fdlWord.localBxNr();

        const DecisionWord& gtDecisionWord = gtReadoutRecord->decisionWord();
        const TechnicalTriggerWord& gtTTWord =
                gtReadoutRecord->technicalTriggerWord();

        int dbitNumber = 0;
        DecisionWord::const_iterator GTdbitItr;
        algo_bits->Fill(-1.); // fill underflow to normalize
        for (GTdbitItr = gtDecisionWord.begin(); GTdbitItr
                != gtDecisionWord.end(); GTdbitItr++) {
            if (*GTdbitItr) {
                algo_bits->Fill(dbitNumber);
                algo_bits_lumi->Fill(e.luminosityBlock(), dbitNumber);
                int dbitNumber1 = 0;
                DecisionWord::const_iterator GTdbitItr1;
                for (GTdbitItr1 = gtDecisionWord.begin(); GTdbitItr1
                        != gtDecisionWord.end(); GTdbitItr1++) {
                    if (*GTdbitItr1)
                        algo_bits_corr->Fill(dbitNumber, dbitNumber1);
                    dbitNumber1++;
                }
                int tbitNumber1 = 0;
                TechnicalTriggerWord::const_iterator GTtbitItr1;
                for (GTtbitItr1 = gtTTWord.begin(); GTtbitItr1
                        != gtTTWord.end(); GTtbitItr1++) {
                    if (*GTtbitItr1)
                        tt_bits_corr->Fill(dbitNumber, tbitNumber1);
                    tbitNumber1++;
                }
            }
            dbitNumber++;
        }

        int tbitNumber = 0;
        TechnicalTriggerWord::const_iterator GTtbitItr;
        tt_bits->Fill(-1.); // fill underflow to normalize
        for (GTtbitItr = gtTTWord.begin(); GTtbitItr != gtTTWord.end(); GTtbitItr++) {
            if (*GTtbitItr) {
                tt_bits->Fill(tbitNumber);
                tt_bits_lumi->Fill(e.luminosityBlock(), tbitNumber);
                int tbitNumber1 = 0;
                TechnicalTriggerWord::const_iterator GTtbitItr1;
                for (GTtbitItr1 = gtTTWord.begin(); GTtbitItr1
                        != gtTTWord.end(); GTtbitItr1++) {
                    if (*GTtbitItr1)
                        tt_bits_corr->Fill(tbitNumber, tbitNumber1);
                    tbitNumber1++;
                }
            }
            tbitNumber++;
        }


        // fill the index of actual prescale factor set
        // the index for technical triggers and algorithm trigger is the same (constraint in L1 GT TS)
        // so we read only pfIndexAlgoTrig (boost::uint16_t)

        const float pfIndexAlgoTrig = static_cast<float>(fdlWord.gtPrescaleFactorIndexAlgo());
        m_monL1PrescaleFactorSet->Fill(e.luminosityBlock(), pfIndexAlgoTrig);


    }

    // get info from active PSB's
    int ibit = PSB9; // first psb
    // for now hardcode psb id's - TODO - get them from Vasile's board maps...
    int psbID[7] = { 0xbb09, 0xbb0d, 0xbb0e, 0xbb0f, 0xbb13, 0xbb14, 0xbb15 };
    for (int i = 0; i < 7; i++) {
        if (isActive(gtfeActiveBoards, ibit)) {
            L1GtPsbWord psbWord = gtReadoutRecord->gtPsbWord(psbID[i]);
            psbBx[0][i] = psbWord.bxNr();
            psbBx[1][i] = psbWord.localBxNr();
        }
        ibit++;
    }

    //fill the dbx histo
    if (gtfeEvmBx > -1)
        dbx_module->Fill(0., gtfeEvmBx - gtfeBx);
    if (tcsBx > -1)
        dbx_module->Fill(1., tcsBx - gtfeBx);
    for (int i = 0; i < 2; i++) {
        if (fdlBx[i] > -1)
            dbx_module->Fill(2. + i, fdlBx[i] - gtfeBx);
    }
    for (int j = 0; j < 7; j++) {
        for (int i = 0; i < 2; i++) {
            if (psbBx[i][j] > -1)
                dbx_module->Fill(4. + i + 2 * j, psbBx[i][j] - gtfeBx);
        }
    }
    if (gmtBx > -1)
        dbx_module->Fill(18., gmtBx - gtfeBx);

}

void L1TGT::beginJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 58 of file L1TGT.cc.

References algo_bits, algo_bits_corr, algo_bits_lumi, algo_tt_bits_corr, DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), BST_beamMode, BST_beamMomentum, BST_intensityBeam1, BST_intensityBeam2, BST_lhcFillNumber, BST_MasterStatus, BST_turnCountNumber, dbe, dbx_module, event_lumi, event_number, event_type, evnum_trignum_lumi, gpsfreq, gpsfreqlum, gpsfreqwide, gtfe_bx, m_monL1PrescaleFactorSet, nev_, cmsCodeRules.cppFunctionSkipper::operator, orbit_lumi, preGps_, preOrb_, DQMStore::rmdir(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), setupversion_lumi, trigger_lumi, trigger_number, tt_bits, tt_bits_corr, and tt_bits_lumi.

                     {

    nev_ = 0;
    preGps_ = 0;
    preOrb_ = 0;

    // get hold of back-end interface
    DQMStore* dbe = 0;
    dbe = Service<DQMStore> ().operator->();

    if (dbe) {
        dbe->setCurrentFolder("L1T/L1TGT");
        dbe->rmdir("L1T/L1TGT");
    }

    if (dbe) {
        dbe->setCurrentFolder("L1T/L1TGT");

        algo_bits = dbe->book1D("algo_bits", "GT algo bits", 128, -0.5, 127.5);
        algo_bits->setAxisTitle("algorithm bits", 1);

        algo_bits_corr = dbe->book2D("algo_bits_corr",
                "GT algo bit correlation",
                128, -0.5, 127.5, 128, -0.5, 127.5);
        algo_bits_corr->setAxisTitle("algorithm bits", 1);
        algo_bits_corr->setAxisTitle("algorithm bits", 2);

        tt_bits = dbe->book1D("tt_bits",
                "GT technical trigger bits",
                64, -0.5, 63.5);
        tt_bits->setAxisTitle("technical trigger bits", 1);

        tt_bits_corr = dbe->book2D("tt_bits_corr",
                "GT tech. trig. bit correlation",
                64, -0.5, 63.5, 64, -0.5, 63.5);
        tt_bits_corr->setAxisTitle("technical trigger bits", 1);
        tt_bits_corr->setAxisTitle("technical trigger bits", 2);

        algo_tt_bits_corr = dbe->book2D("algo_tt_bits_corr",
                "GT algo tech. trig. bit correlation",
                128, -0.5, 127.5, 64, -0.5, 63.5);
        algo_tt_bits_corr->setAxisTitle("algorithm bits", 1);
        algo_tt_bits_corr->setAxisTitle("technical trigger bits", 2);

        algo_bits_lumi = dbe->book2D("algo_bits_lumi",
                "GT algo bit rate per lumi segment",
                250, 0., 250., 128, -0.5, 127.5);
        algo_bits_lumi->setAxisTitle("luminosity segment", 1);
        algo_bits_lumi->setAxisTitle("algorithm bits", 2);

        tt_bits_lumi = dbe->book2D("tt_bits_lumi",
                "GT tech. trig. bit rate per lumi segment",
                250, 0., 250., 64, -0.5, 63.5);
        tt_bits_lumi->setAxisTitle("luminosity segment", 1);
        tt_bits_lumi->setAxisTitle("technical trigger bits", 2);

        event_type = dbe->book1D("event_type", "GT event type", 10, -0.5, 9.5);
        event_type->setAxisTitle("event type", 1);
        event_type->setBinLabel(2, "Physics", 1);
        event_type->setBinLabel(3, "Calibration", 1);
        event_type->setBinLabel(4, "Random", 1);
        event_type->setBinLabel(6, "Traced", 1);
        event_type->setBinLabel(7, "Test", 1);
        event_type->setBinLabel(8, "Error", 1);

        event_number = dbe->book1D("event_number",
                "GT Event number (from last resync)",
                100, 0., 50000.);
        event_number->setAxisTitle("event number", 1);

        event_lumi = dbe->bookProfile("event_lumi",
                "GT Event number (from last resync) vs lumi section",
                250, 0., 250., 100, -0.1, 1.e15, "s");
        event_lumi->setAxisTitle("luminosity segment", 1);
        event_lumi->setAxisTitle("event number", 2);

        trigger_number = dbe->book1D("trigger_number",
                "GT Trigger number (from start run)",
                100, 0., 50000.);
        trigger_number->setAxisTitle("trigger number", 1);

        trigger_lumi = dbe->bookProfile("trigger_lumi",
                "GT Trigger number (from start run) vs lumi section",
                250, 0., 250., 100, -0.1, 1.e15, "s");
        trigger_lumi->setAxisTitle("luminosity segment", 1);
        trigger_lumi->setAxisTitle("trigger number", 2);

        evnum_trignum_lumi = dbe->bookProfile("evnum_trignum_lumi",
                "GT Event/Trigger number ratio vs lumi section",
                250, 0., 250., 100, -0.1, 2., "s");
        evnum_trignum_lumi->setAxisTitle("luminosity segment", 1);
        evnum_trignum_lumi->setAxisTitle("event/trigger number ratio", 2);

        orbit_lumi = dbe->bookProfile("orbit_lumi",
                "GT orbit number vs lumi section",
                250, 0., 250., 100, -0.1, 1.e15, "s");
        orbit_lumi->setAxisTitle("luminosity segment", 1);
        orbit_lumi->setAxisTitle("orbit number", 2);

        setupversion_lumi = dbe->bookProfile("setupversion_lumi",
                "GT setup version vs lumi section",
                250, 0., 250., 100, -0.1, 1.e10, "i");
        setupversion_lumi->setAxisTitle("luminosity segment", 1);
        setupversion_lumi->setAxisTitle("prescale stup version", 2);

        gtfe_bx = dbe->book1D("gtfe_bx", "GTFE Bx number", 3600, 0., 3600.);
        gtfe_bx->setAxisTitle("GTFE BX number", 1);

        dbx_module = dbe->bookProfile("dbx_module",
                "delta Bx of GT modules wrt. GTFE",
                20, 0., 20., 100, -4000., 4000., "i");
        dbx_module->setAxisTitle("GT crate module", 1);
        dbx_module->setAxisTitle("Module Bx - GTFE Bx", 2);
        dbx_module->setBinLabel(1, "GTFEevm", 1);
        dbx_module->setBinLabel(2, "TCS", 1);
        dbx_module->setBinLabel(3, "FDL", 1);
        dbx_module->setBinLabel(4, "FDLloc", 1);
        dbx_module->setBinLabel(5, "PSB9", 1);
        dbx_module->setBinLabel(6, "PSB9loc", 1);
        dbx_module->setBinLabel(7, "PSB13", 1);
        dbx_module->setBinLabel(8, "PSB13loc", 1);
        dbx_module->setBinLabel(9, "PSB14", 1);
        dbx_module->setBinLabel(10, "PSB14loc", 1);
        dbx_module->setBinLabel(11, "PSB15", 1);
        dbx_module->setBinLabel(12, "PSB15loc", 1);
        dbx_module->setBinLabel(13, "PSB19", 1);
        dbx_module->setBinLabel(14, "PSB19loc", 1);
        dbx_module->setBinLabel(15, "PSB20", 1);
        dbx_module->setBinLabel(16, "PSB20loc", 1);
        dbx_module->setBinLabel(17, "PSB21", 1);
        dbx_module->setBinLabel(18, "PSB21loc", 1);
        dbx_module->setBinLabel(19, "GMT", 1);

        BST_MasterStatus = dbe->book2D("BST_MasterStatus",
                "BST master status over lumi segment",
                250, 0., 250., 6, -1., 5.);
        BST_MasterStatus->setAxisTitle("luminosity segment", 1);
        BST_MasterStatus->setAxisTitle("BST master status", 2);
        BST_MasterStatus->setBinLabel(2, "Master Beam 1", 2);
        BST_MasterStatus->setBinLabel(3, "Master Beam 2", 2);

        BST_turnCountNumber = dbe->book2D("BST_turnCountNumber",
                "BST turn count over lumi segment",
                250, 0., 250., 250, 0., 4.3e9);
        BST_turnCountNumber->setAxisTitle("luminosity segment", 1);
        BST_turnCountNumber->setAxisTitle("BST turn count number", 2);

        BST_lhcFillNumber = dbe->book1D("BST_lhcFillNumber",
                "BST LHC fill number % 1000", 1000, 0., 1000.);
        BST_lhcFillNumber->setAxisTitle("BST LHC fill number modulo 1000");

        BST_beamMode = dbe->book2D("BST_beamMode",
                "BST beam mode over lumi segment",
                250, 0., 250., 25, 1., 26.);
        BST_beamMode->setAxisTitle("luminosity segment", 1);
        BST_beamMode->setAxisTitle("mode", 2);
        BST_beamMode->setBinLabel(1, "no mode", 2);
        BST_beamMode->setBinLabel(2, "setup", 2);
        BST_beamMode->setBinLabel(3, "inj pilot", 2);
        BST_beamMode->setBinLabel(4, "inj intr", 2);
        BST_beamMode->setBinLabel(5, "inj nomn", 2);
        BST_beamMode->setBinLabel(6, "pre ramp", 2);
        BST_beamMode->setBinLabel(7, "ramp", 2);
        BST_beamMode->setBinLabel(8, "flat top", 2);
        BST_beamMode->setBinLabel(9, "squeeze", 2);
        BST_beamMode->setBinLabel(10, "adjust", 2);
        BST_beamMode->setBinLabel(11, "stable", 2);
        BST_beamMode->setBinLabel(12, "unstable", 2);
        BST_beamMode->setBinLabel(13, "beam dump", 2);
        BST_beamMode->setBinLabel(14, "ramp down", 2);
        BST_beamMode->setBinLabel(15, "recovery", 2);
        BST_beamMode->setBinLabel(16, "inj dump", 2);
        BST_beamMode->setBinLabel(17, "circ dump", 2);
        BST_beamMode->setBinLabel(18, "abort", 2);
        BST_beamMode->setBinLabel(19, "cycling", 2);
        BST_beamMode->setBinLabel(20, "warn beam dump", 2);
        BST_beamMode->setBinLabel(21, "no beam", 2);

        BST_beamMomentum = dbe->book2D("BST_beamMomentum",
                "BST beam momentum",
                250, 0., 250., 100, 0., 7200.);
        BST_beamMomentum->setAxisTitle("luminosity segment", 1);
        BST_beamMomentum->setAxisTitle("beam momentum", 2);

        gpsfreq = dbe->book1D("gpsfreq", "clock frequency measured by GPS",
                1000, 39.95, 40.2);
        gpsfreq->setAxisTitle("CMS clock frequency (MHz)");

        gpsfreqwide = dbe->book1D("gpsfreqwide",
                "clock frequency measured by GPS", 1000, -2., 200.);
        gpsfreqwide->setAxisTitle("CMS clock frequency (MHz)");

        gpsfreqlum = dbe->book2D("gpsfreqlum",
                "clock frequency measured by GPS",
                250, 0., 250., 100, 39.95, 40.2);
        gpsfreqlum->setAxisTitle("luminosity segment", 1);
        gpsfreqlum->setAxisTitle("CMS clock frequency (MHz)", 2);

        BST_intensityBeam1 = dbe->book2D("BST_intensityBeam1",
                "intensity beam 1",
                250, 0., 250., 100, 0., 100.);
        BST_intensityBeam1->setAxisTitle("luminosity segment", 1);
        BST_intensityBeam1->setAxisTitle("beam intensity", 2);

        BST_intensityBeam2 = dbe->book2D("BST_intensityBeam2",
                "intensity beam 2",
                250, 0., 250., 100, 0., 100.);
        BST_intensityBeam2->setAxisTitle("luminosity segment", 1);
        BST_intensityBeam2->setAxisTitle("beam intensity", 2);

        m_monL1PrescaleFactorSet = dbe->book2D("L1PrescaleFactorSet",
                "Index of L1 prescale factor set",
                250, 0., 250., 25, 0., 25.);
        m_monL1PrescaleFactorSet->setAxisTitle("Luminosity segment", 1);
        m_monL1PrescaleFactorSet->setAxisTitle("L1 PF set index", 2);

    }
}

void L1TGT::endJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 277 of file L1TGT.cc.

References gather_cfg::cout, dbe, nev_, outputFile_, DQMStore::save(), and verbose_.

                       {

    if (verbose_) {
        cout << "L1TGT: end job...." << endl;
    }

    LogInfo("EndJob") << "analyzed " << nev_ << " events";

    if (outputFile_.size() != 0 && dbe)
        dbe->save(outputFile_);

    return;
}

bool L1TGT::isActive ( int  word, int  bit  )

private

Definition at line 534 of file L1TGT.cc.

Referenced by analyze().

                                      {
    if (word & (1 << bit))
        return true;
    return false;
}

Member Data Documentation

MonitorElement* L1TGT::algo_bits

private

Definition at line 81 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::algo_bits_corr

private

Definition at line 82 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::algo_bits_lumi

private

Definition at line 86 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::algo_tt_bits_corr

private

Definition at line 85 of file L1TGT.h.

Referenced by beginJob().

MonitorElement* L1TGT::BST_beamMode

private

Definition at line 104 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::BST_beamMomentum

private

Definition at line 105 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::BST_intensityBeam1

private

Definition at line 106 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::BST_intensityBeam2

private

Definition at line 107 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::BST_lhcFillNumber

private

Definition at line 103 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::BST_MasterStatus

private

Definition at line 101 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::BST_turnCountNumber

private

Definition at line 102 of file L1TGT.h.

Referenced by analyze(), and beginJob().

DQMStore* L1TGT::dbe

private

Definition at line 79 of file L1TGT.h.

Referenced by beginJob(), endJob(), and L1TGT().

MonitorElement* L1TGT::dbx_module

private

Definition at line 99 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::event_lumi

private

Definition at line 91 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::event_number

private

Definition at line 90 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::event_type

private

Definition at line 88 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::evnum_trignum_lumi

private

Definition at line 94 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::gpsfreq

private

Definition at line 108 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::gpsfreqlum

private

Definition at line 110 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::gpsfreqwide

private

Definition at line 109 of file L1TGT.h.

Referenced by analyze(), and beginJob().

edm::InputTag L1TGT::gtEvmSource_

private

Definition at line 123 of file L1TGT.h.

Referenced by analyze().

MonitorElement* L1TGT::gtfe_bx

private

Definition at line 98 of file L1TGT.h.

Referenced by analyze(), and beginJob().

edm::InputTag L1TGT::gtSource_

private

Definition at line 122 of file L1TGT.h.

Referenced by analyze().

ofstream L1TGT::logFile_

private

Definition at line 121 of file L1TGT.h.

MonitorElement* L1TGT::m_monL1PrescaleFactorSet

private

Definition at line 112 of file L1TGT.h.

Referenced by analyze(), and beginJob().

bool L1TGT::monitorDaemon_

private

Definition at line 120 of file L1TGT.h.

int L1TGT::nev_

private

Definition at line 117 of file L1TGT.h.

Referenced by analyze(), beginJob(), and endJob().

MonitorElement* L1TGT::orbit_lumi

private

Definition at line 95 of file L1TGT.h.

Referenced by analyze(), and beginJob().

std::string L1TGT::outputFile_

private

Definition at line 118 of file L1TGT.h.

Referenced by endJob(), and L1TGT().

boost::uint64_t L1TGT::preGps_

private

Definition at line 125 of file L1TGT.h.

Referenced by analyze(), and beginJob().

boost::uint64_t L1TGT::preOrb_

private

Definition at line 126 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::setupversion_lumi

private

Definition at line 96 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::trigger_lumi

private

Definition at line 93 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::trigger_number

private

Definition at line 92 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::tt_bits

private

Definition at line 83 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::tt_bits_corr

private

Definition at line 84 of file L1TGT.h.

Referenced by analyze(), and beginJob().

MonitorElement* L1TGT::tt_bits_lumi

private

Definition at line 87 of file L1TGT.h.

Referenced by analyze(), and beginJob().

bool L1TGT::verbose_

private

Definition at line 119 of file L1TGT.h.

Referenced by analyze(), endJob(), and L1TGT().