#include <DQMOffline/DQMHOAlCaRecoStream/src/DQMHOAlCaRecoStream.cc>

Inheritance diagram for DQMHOAlCaRecoStream:

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

	~DQMHOAlCaRecoStream () override

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

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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 edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

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

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

Definition at line 10 of file DQMHOAlCaRecoStream.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 64 of file DQMHOAlCaRecoStream.cc.

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

                                                                        :
   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);
 }

DQMHOAlCaRecoStream::~DQMHOAlCaRecoStream ( )

override

Definition at line 81 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

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

overrideprivate

Definition at line 96 of file DQMHOAlCaRecoStream.cc.

References funct::abs(), runTauDisplay::dr, MonitorElement::Fill(), edm::Event::getByToken(), hDirCosine, hHOTime, hMuonEta, hMuonMom, hMuonMultipl, hMuonPhi, hoCalibVariableCollectionTag, hPedRing, hSignal3x3, hSigRing, createfilelist::int, edm::HandleBase::isValid(), gen::k, cmsBatch::log, LogDebug, m_highRadPosInMuch, Nevents, Nmuons, funct::pow(), mathSSE::return(), and funct::tan().

 {
   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) { 
 }

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 161 of file DQMHOAlCaRecoStream.cc.

References DQMStore::IBooker::book1D(), folderName_, hDirCosine, hHOTime, hMuonEta, hMuonMom, hMuonMultipl, hMuonPhi, hPedRing, hSignal3x3, hSigRing, mps_fire::i, gen::k, m_highEdge, m_lowEdge, m_nbins, dataset::name, Nevents, Nmuons, MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), and fftjetcommon_cfi::title.

 {
   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;
 
 }