DPGAnalysis/SiStripTools/plugins/APVCyclePhaseMonitor.cc - master - githubmirror/cmssw

DPGAnalysis/SiStripTools/plugins/APVCyclePhaseMonitor.cc
286 строк11 KB
Christopher Jones
Use thread safe modules in DPGAnalysis/SiStripTools
01 дек 2021, 17:52
01 дек 2021, 17:5271fbca3
// -*- C++ -*-
//
// Package:    SiStripTools
// Class:      APVCyclePhaseMonitor
//
/**\class APVCyclePhaseMonitor APVCyclePhaseMonitor.cc DPGAnalysis/SiStripTools/plugins/APVCyclePhaseMonitor.cc

 Description: <one line class summary>

 Implementation:
     <Notes on implementation>
*/
//
// Original Author:  Andrea Venturi
//         Created:  Tue Jul 19 11:56:00 CEST 2009
//
//

// system include files
#include <memory>

// user include files
#include "TH1F.h"
#include "TProfile.h"

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/Run.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"

#include "FWCore/Utilities/interface/InputTag.h"

#include "DPGAnalysis/SiStripTools/interface/APVCyclePhaseCollection.h"

#include "DPGAnalysis/SiStripTools/interface/RunHistogramManager.h"

//
// class decleration
//

class APVCyclePhaseMonitor : public edm::one::EDAnalyzer<edm::one::WatchRuns, edm::one::SharedResources> {
public:
  explicit APVCyclePhaseMonitor(const edm::ParameterSet&);
  ~APVCyclePhaseMonitor() override;

private:
  void beginJob() override;
  void beginRun(const edm::Run&, const edm::EventSetup&) override;
  void endRun(const edm::Run&, const edm::EventSetup&) override;
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void endJob() override;

  // ----------member data ---------------------------

  edm::EDGetTokenT<APVCyclePhaseCollection> _apvphasecollectionToken;
  std::vector<std::string> _selectedparts;
  std::vector<std::string> _selectedvectorparts;
  const unsigned int m_maxLS;
  const unsigned int m_LSfrac;
  RunHistogramManager m_rhm;
  std::map<std::string, TH1F*> _hphases;
  std::map<std::string, TH1F**> _hselectedphases;
  std::map<std::string, TH1F**> _hselectedphasesvector;
  std::map<std::string, TH1F**> _hselectedphasessize;
  std::map<std::string, TProfile*> _hphasevsorbit;
  std::map<std::string, TProfile**> _hselectedphasevsorbit;
  std::map<std::string, TProfile**> _hselectedphasevectorvsorbit;
  unsigned int _nevents;
};

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
APVCyclePhaseMonitor::APVCyclePhaseMonitor(const edm::ParameterSet& iConfig)
    : _apvphasecollectionToken(
          consumes<APVCyclePhaseCollection>(iConfig.getParameter<edm::InputTag>("apvCyclePhaseCollection"))),
      _selectedparts(
          iConfig.getUntrackedParameter<std::vector<std::string> >("selectedPartitions", std::vector<std::string>())),
      _selectedvectorparts(iConfig.getUntrackedParameter<std::vector<std::string> >("selectedVectorPartitions",
                                                                                    std::vector<std::string>())),
      m_maxLS(iConfig.getUntrackedParameter<unsigned int>("maxLSBeforeRebin", 125)),
      m_LSfrac(iConfig.getUntrackedParameter<unsigned int>("startingLSFraction", 16)),
      m_rhm(consumesCollector()),
      _hphases(),
      _hselectedphases(),
      _hselectedphasesvector(),
      _hselectedphasessize(),
      _hphasevsorbit(),
      _hselectedphasevsorbit(),
      _hselectedphasevectorvsorbit(),
      _nevents(0) {
  //now do what ever initialization is needed

  usesResource(TFileService::kSharedResource);

  edm::LogInfo("UsedAPVCyclePhaseCollection")
      << " APVCyclePhaseCollection " << iConfig.getParameter<edm::InputTag>("apvCyclePhaseCollection") << " used";

  for (std::vector<std::string>::const_iterator part = _selectedparts.begin(); part != _selectedparts.end(); ++part) {
    char hname[300];

    sprintf(hname, "selected_phase_%s", part->c_str());
    edm::LogInfo("SelectedTH1FBeingBooked") << "TH1F " << hname << " being booked";
    _hselectedphases[*part] = m_rhm.makeTH1F(hname, hname, 70, -0.5, 69.5);

    sprintf(hname, "selected_phasevsorbit_%s", part->c_str());
    edm::LogInfo("SelectedTProfileBeingBooked") << "TProfile " << hname << " being booked";
    _hselectedphasevsorbit[*part] = m_rhm.makeTProfile(hname, hname, m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
  }

  for (std::vector<std::string>::const_iterator part = _selectedvectorparts.begin(); part != _selectedvectorparts.end();
       ++part) {
    char hname[300];

    sprintf(hname, "selected_phase_vector_%s", part->c_str());
    edm::LogInfo("SelectedVectTH1FBeingBooked") << "TH1F " << hname << " being booked";
    _hselectedphasesvector[*part] = m_rhm.makeTH1F(hname, hname, 70, -0.5, 69.5);

    sprintf(hname, "selected_phase_size_%s", part->c_str());
    edm::LogInfo("SelectedVectSizeTH1FBeingBooked") << "TH1F " << hname << " being booked";
    _hselectedphasessize[*part] = m_rhm.makeTH1F(hname, hname, 10, -0.5, 9.5);

    sprintf(hname, "selected_phasevectorvsorbit_%s", part->c_str());
    edm::LogInfo("SelectedVectTProfileBeingBooked") << "TProfile " << hname << " being booked";
    _hselectedphasevectorvsorbit[*part] = m_rhm.makeTProfile(hname, hname, m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
  }
}

APVCyclePhaseMonitor::~APVCyclePhaseMonitor() {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

//
// member functions
//

// ------------ method called to for each event  ------------
void APVCyclePhaseMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  using namespace edm;

  _nevents++;

  edm::Handle<APVCyclePhaseCollection> apvphases;
  iEvent.getByToken(_apvphasecollectionToken, apvphases);

  // improve the matchin between default and actual partitions

  edm::Service<TFileService> tfserv;

  for (std::map<std::string, int>::const_iterator phase = apvphases->get().begin(); phase != apvphases->get().end();
       ++phase) {
    if (_hphases.find(phase->first) == _hphases.end()) {
      char dirname[300];
      sprintf(dirname, "run_%d", iEvent.run());
      TFileDirectory subrun = tfserv->mkdir(dirname);

      char hname[300];

      sprintf(hname, "phase_%s", phase->first.c_str());
      edm::LogInfo("TH1FBeingBooked") << "TH1F " << hname << " being booked";
      _hphases[phase->first] = subrun.make<TH1F>(hname, hname, 70, -0.5, 69.5);
      _hphases[phase->first]->GetXaxis()->SetTitle("BX mod 70");
      _hphases[phase->first]->GetYaxis()->SetTitle("Events");

      sprintf(hname, "phasevsorbit_%s", phase->first.c_str());
      edm::LogInfo("TProfileBeingBooked") << "TProfile " << hname << " being booked";
      _hphasevsorbit[phase->first] = subrun.make<TProfile>(hname, hname, m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
      _hphasevsorbit[phase->first]->SetCanExtend(TH1::kXaxis);
      _hphasevsorbit[phase->first]->GetXaxis()->SetTitle("time [orbit#]");
      _hphasevsorbit[phase->first]->GetYaxis()->SetTitle("Phase");
    }
    _hphases[phase->first]->Fill(phase->second);
    _hphasevsorbit[phase->first]->Fill(iEvent.orbitNumber(), phase->second);
  }

  // selected partitions

  for (std::vector<std::string>::const_iterator part = _selectedparts.begin(); part != _selectedparts.end(); ++part) {
    if (_hselectedphases.find(*part) != _hselectedphases.end() && _hselectedphases[*part] && *_hselectedphases[*part]) {
      (*_hselectedphases[*part])->Fill(apvphases->getPhase(*part));
    }
    if (_hselectedphasevsorbit.find(*part) != _hselectedphasevsorbit.end() && _hselectedphasevsorbit[*part] &&
        *_hselectedphasevsorbit[*part]) {
      (*_hselectedphasevsorbit[*part])->Fill(iEvent.orbitNumber(), apvphases->getPhase(*part));
    }
  }

  for (std::vector<std::string>::const_iterator part = _selectedvectorparts.begin(); part != _selectedvectorparts.end();
       ++part) {
    const std::vector<int> phases = apvphases->getPhases(*part);

    if (_hselectedphasessize.find(*part) != _hselectedphasessize.end() && _hselectedphasessize[*part] &&
        *_hselectedphasessize[*part]) {
      (*_hselectedphasessize[*part])->Fill(phases.size());
    }

    for (std::vector<int>::const_iterator phase = phases.begin(); phase != phases.end(); ++phase) {
      if (_hselectedphasesvector.find(*part) != _hselectedphasesvector.end() && _hselectedphasesvector[*part] &&
          *_hselectedphasesvector[*part]) {
        (*_hselectedphasesvector[*part])->Fill(*phase);
      }
      if (_hselectedphasevectorvsorbit.find(*part) != _hselectedphasevectorvsorbit.end() &&
          _hselectedphasevectorvsorbit[*part] && *_hselectedphasevectorvsorbit[*part]) {
        (*_hselectedphasevectorvsorbit[*part])->Fill(iEvent.orbitNumber(), *phase);
      }
    }
  }
}

void APVCyclePhaseMonitor::beginRun(const edm::Run& iRun, const edm::EventSetup&) {
  _hphases.clear();
  _hphasevsorbit.clear();

  m_rhm.beginRun(iRun);

  for (std::map<std::string, TH1F**>::const_iterator hist = _hselectedphases.begin(); hist != _hselectedphases.end();
       ++hist) {
    if (*(hist->second)) {
      (*(hist->second))->GetXaxis()->SetTitle("BX mod 70");
      (*(hist->second))->GetYaxis()->SetTitle("Events");
    }
  }
  for (std::map<std::string, TProfile**>::const_iterator prof = _hselectedphasevsorbit.begin();
       prof != _hselectedphasevsorbit.end();
       ++prof) {
    if (*(prof->second)) {
      (*(prof->second))->SetCanExtend(TH1::kXaxis);
      (*(prof->second))->GetXaxis()->SetTitle("time [orbit#]");
      (*(prof->second))->GetYaxis()->SetTitle("Phase");
    }
  }
  for (std::map<std::string, TH1F**>::const_iterator hist = _hselectedphasesvector.begin();
       hist != _hselectedphasesvector.end();
       ++hist) {
    if (*(hist->second)) {
      (*(hist->second))->GetXaxis()->SetTitle("BX mod 70");
      (*(hist->second))->GetYaxis()->SetTitle("Events");
    }
  }
  for (std::map<std::string, TH1F**>::const_iterator hist = _hselectedphasessize.begin();
       hist != _hselectedphasessize.end();
       ++hist) {
    if (*(hist->second)) {
      (*(hist->second))->GetXaxis()->SetTitle("Number of Phases");
      (*(hist->second))->GetYaxis()->SetTitle("Events");
    }
  }
  for (std::map<std::string, TProfile**>::const_iterator prof = _hselectedphasevectorvsorbit.begin();
       prof != _hselectedphasevectorvsorbit.end();
       ++prof) {
    if (*(prof->second)) {
      (*(prof->second))->SetCanExtend(TH1::kXaxis);
      (*(prof->second))->GetXaxis()->SetTitle("time [orbit#]");
      (*(prof->second))->GetYaxis()->SetTitle("Phase");
    }
  }
}

void APVCyclePhaseMonitor::endRun(const edm::Run& iRun, const edm::EventSetup&) {}

// ------------ method called once each job just before starting event loop  ------------
void APVCyclePhaseMonitor::beginJob() {}

// ------------ method called once each job just after ending the event loop  ------------
void APVCyclePhaseMonitor::endJob() { edm::LogInfo("EndOfJob") << _nevents << " analyzed events"; }

//define this as a plug-in
DEFINE_FWK_MODULE(APVCyclePhaseMonitor);