#include <DQMHOAlCaRecoStream.h>

Inheritance diagram for DQMHOAlCaRecoStream:

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

	~DQMHOAlCaRecoStream () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Private Attributes
std::string	folderName_

MonitorElement *	hDirCosine

MonitorElement *	hHOTime

MonitorElement *	hMuonEta

MonitorElement *	hMuonMom

MonitorElement *	hMuonMultipl

MonitorElement *	hMuonPhi

edm::EDGetTokenT< HOCalibVariableCollection >	hoCalibVariableCollectionTag

MonitorElement *	hPedRing [5]

MonitorElement *	hSignal3x3 [9]

MonitorElement *	hSigRing [5]

double	m_highEdge

double	m_highRadPosInMuch

double	m_lowEdge

double	m_lowRadPosInMuch

int	m_nbins

double	m_sigmaValue

int	Nevents

int	Nmuons

bool	saveToFile_

std::string	theRootFileName

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

DQMOffline/DQMHOAlCaRecoStream/src/DQMHOAlCaRecoStream.cc

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 9 of file DQMHOAlCaRecoStream.h.

Constructor & Destructor Documentation

◆ DQMHOAlCaRecoStream()

DQMHOAlCaRecoStream::DQMHOAlCaRecoStream ( const edm::ParameterSet & iConfig )

explicit

Definition at line 57 of file DQMHOAlCaRecoStream.cc.

     : hoCalibVariableCollectionTag(
           consumes<HOCalibVariableCollection>(iConfig.getParameter<edm::InputTag>("hoCalibVariableCollectionTag"))) {
   // now do what ever initialization is needed
  
   theRootFileName = iConfig.getUntrackedParameter<string>("RootFileName", "tmp.root");
   folderName_ = iConfig.getUntrackedParameter<string>("folderName");
   m_sigmaValue = iConfig.getUntrackedParameter<double>("sigmaval", 0.2);
   m_lowRadPosInMuch = iConfig.getUntrackedParameter<double>("lowradposinmuch", 400.0);
   m_highRadPosInMuch = iConfig.getUntrackedParameter<double>("highradposinmuch", 480.0);
   m_lowEdge = iConfig.getUntrackedParameter<double>("lowedge", -2.0);
   m_highEdge = iConfig.getUntrackedParameter<double>("highedge", 6.0);
   m_nbins = iConfig.getUntrackedParameter<int>("nbins", 40);
   saveToFile_ = iConfig.getUntrackedParameter<bool>("saveToFile", false);
 }

References folderName_, edm::ParameterSet::getUntrackedParameter(), m_highEdge, m_highRadPosInMuch, m_lowEdge, m_lowRadPosInMuch, m_nbins, m_sigmaValue, saveToFile_, and theRootFileName.

◆ ~DQMHOAlCaRecoStream()

DQMHOAlCaRecoStream::~DQMHOAlCaRecoStream ( )

override

Definition at line 73 of file DQMHOAlCaRecoStream.cc.

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

Member Function Documentation

◆ analyze()

void DQMHOAlCaRecoStream::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 83 of file DQMHOAlCaRecoStream.cc.

                                                                                      {
   using namespace edm;
  
   Nevents++;
  
   edm::Handle<HOCalibVariableCollection> HOCalib;
   bool isCosMu = true;
  
   iEvent.getByToken(hoCalibVariableCollectionTag, HOCalib);
  
   if (!HOCalib.isValid()) {
     LogDebug("") << "DQMHOAlCaRecoStream:: Error! can't get HOCalib product!" << std::endl;
     return;
   }
  
   if (isCosMu) {
     hMuonMultipl->Fill((*HOCalib).size(), 1.);
     if (!(*HOCalib).empty()) {
       for (HOCalibVariableCollection::const_iterator hoC = (*HOCalib).begin(); hoC != (*HOCalib).end(); hoC++) {
         // OK!!!!
         float okt = 2.;
         double okx = std::pow((*hoC).trkvx, okt) + std::pow((*hoC).trkvy, okt);
         double dr = std::pow(okx, 0.5);
         if (dr < m_lowRadPosInMuch || dr > m_highRadPosInMuch)
           continue;
  
         if ((*hoC).isect < 0)
           continue;
         if (fabs((*hoC).trkth - acos(-1.) / 2) < 0.000001)
           continue;
         int ieta = int((std::abs((*hoC).isect) % 10000) / 100.) - 30;
  
         if (std::abs(ieta) >= 16)
           continue;
  
         Nmuons++;
  
         hMuonMom->Fill((*hoC).trkmm, 1.0);
         hMuonEta->Fill(-log(tan((*hoC).trkth / 2.0)), 1.0);
         hMuonPhi->Fill((*hoC).trkph, 1.0);
         hDirCosine->Fill((*hoC).hoang, 1.0);
         hHOTime->Fill((*hoC).htime, 1.0);
  
         double energy = (*hoC).hosig[4];
         double pedval = (*hoC).hocro;
         int iring = 0;
         if (ieta >= -15 && ieta <= -11) {
           iring = -2;
         } else if (ieta >= -10 && ieta <= -5) {
           iring = -1;
         } else if (ieta >= 5 && ieta <= 10) {
           iring = 1;
         } else if (ieta >= 11 && ieta <= 15) {
           iring = 2;
         }
  
         hSigRing[iring + 2]->Fill(energy, 1.0);
         hPedRing[iring + 2]->Fill(pedval, 1.0);
  
         for (int k = 0; k < 9; k++) {
           hSignal3x3[k]->Fill((*hoC).hosig[k]);
         }
       }  // for (HOCalibVariableCollection::const_iterator hoC=(*HOCalib).begin()
     }    // if ((*HOCalib).size() >0 ) {
   }      // if (isCosMu) {
 }

References funct::abs(), flavorHistoryFilter_cfi::dr, HCALHighEnergyHPDFilter_cfi::energy, dqm::impl::MonitorElement::Fill(), hDirCosine, hHOTime, hMuonEta, hMuonMom, hMuonMultipl, hMuonPhi, hoCalibVariableCollectionTag, hPedRing, hSignal3x3, hSigRing, LEDCalibrationChannels::ieta, iEvent, createfilelist::int, edm::HandleBase::isValid(), dqmdumpme::k, dqm-mbProfile::log, LogDebug, m_highRadPosInMuch, Nevents, Nmuons, funct::pow(), and funct::tan().

◆ bookHistograms()

void DQMHOAlCaRecoStream::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	irun,
		edm::EventSetup const &	isetup
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 153 of file DQMHOAlCaRecoStream.cc.

                                                                       {
   ibooker.setCurrentFolder(folderName_);
  
   char title[200];
   char name[200];
  
   hMuonMom = ibooker.book1D("hMuonMom", "Muon momentum (GeV)", 50, -100, 100);
   hMuonMom->setAxisTitle("Muon momentum (GeV)", 1);
  
   hMuonEta = ibooker.book1D("hMuonEta", "Pseudo-rapidity of muon", 50, -1.5, 1.5);
   hMuonEta->setAxisTitle("Pseudo-rapidity of muon", 1);
  
   hMuonPhi = ibooker.book1D("hMuonPhi", "Azimuthal angle of muon", 24, -acos(-1), acos(-1));
   hMuonPhi->setAxisTitle("Azimuthal angle of muon", 1);
  
   hMuonMultipl = ibooker.book1D("hMuonMultipl", "Muon Multiplicity", 10, 0.5, 10.5);
   hMuonMultipl->setAxisTitle("Muon Multiplicity", 1);
  
   hDirCosine = ibooker.book1D("hDirCosine", "Direction Cosine of muon at HO tower", 50, -1., 1.);
   hDirCosine->setAxisTitle("Direction Cosine of muon at HO tower", 1);
  
   hHOTime = ibooker.book1D("hHOTime", "HO time distribution", 60, -20, 100.);
   hHOTime->setAxisTitle("HO time distribution", 1);
  
   for (int i = 0; i < 5; i++) {
     sprintf(name, "hSigRing_%i", i - 2);
     sprintf(title, "HO signal in Ring_%i", i - 2);
     hSigRing[i] = ibooker.book1D(name, title, m_nbins, m_lowEdge, m_highEdge);
     hSigRing[i]->setAxisTitle(title, 1);
  
     sprintf(name, "hPedRing_%i", i - 2);
     sprintf(title, "HO Pedestal in Ring_%i", i - 2);
     hPedRing[i] = ibooker.book1D(name, title, m_nbins, m_lowEdge, m_highEdge);
     hPedRing[i]->setAxisTitle(title, 1);
   }
  
   //  hSigRingm1 = ibooker.book1D("hSigRingm1", "HO signal in Ring-1", m_nbins,
   //  m_lowEdge, m_highEdge); hSigRingm1->setAxisTitle("HO signal in Ring-1",1);
  
   //  hSigRing00 = ibooker.book1D("hSigRing00", "HO signal in Ring_0", m_nbins,
   //  m_lowEdge, m_highEdge); hSigRing00->setAxisTitle("HO signal in Ring_0",1);
  
   //  hSigRingp1 = ibooker.book1D("hSigRingp1", "HO signal in Ring-1", m_nbins,
   //  m_lowEdge, m_highEdge); hSigRingp1->setAxisTitle("HO signal in Ring+1",1);
  
   //  hSigRingp2 = ibooker.book1D("hSigRingp2", "HO signal in Ring-2", m_nbins,
   //  m_lowEdge, m_highEdge); hSigRingp2->setAxisTitle("HO signal in Ring+2",1);
  
   //  hPedRingm2 = ibooker.book1D("hPedRingm2", "HO pedestal in Ring-2",
   //  m_nbins, m_lowEdge, m_highEdge); hPedRingm1 = ibooker.book1D("hPedRingm1",
   //  "HO pedestal in Ring-1", m_nbins, m_lowEdge, m_highEdge); hPedRing00 =
   //  ibooker.book1D("hPedRing00", "HO pedestal in Ring_0", m_nbins, m_lowEdge,
   //  m_highEdge); hPedRingp1 = ibooker.book1D("hPedRingp1", "HO pedestal in
   //  Ring-1", m_nbins, m_lowEdge, m_highEdge); hPedRingp2 =
   //  ibooker.book1D("hPedRingp2", "HO pedestal in Ring-2", m_nbins, m_lowEdge,
   //  m_highEdge);
  
   for (int i = -1; i <= 1; i++) {
     for (int j = -1; j <= 1; j++) {
       int k = 3 * (i + 1) + j + 1;
  
       sprintf(title, "hSignal3x3_deta%i_dphi%i", i, j);
       hSignal3x3[k] = ibooker.book1D(title, title, m_nbins, m_lowEdge, m_highEdge);
       hSignal3x3[k]->setAxisTitle(title, 1);
     }
   }
  
   Nevents = 0;
   Nmuons = 0;
 }