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

Inheritance diagram for DQMHOAlCaRecoStream:

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

	~DQMHOAlCaRecoStream ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

virtual void	beginJob ()

virtual void	endJob ()

Private Attributes
DQMStore *	dbe_

std::string	folderName_

MonitorElement *	hDirCosine

MonitorElement *	hHOTime

MonitorElement *	hMuonEta

MonitorElement *	hMuonMom

MonitorElement *	hMuonMultipl

MonitorElement *	hMuonPhi

edm::InputTag	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::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)

static void	prevalidate (ConfigurationDescriptions &)

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

CurrentProcessingContext const *	currentContext () const

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)

Detailed Description

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

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

Definition at line 5 of file DQMHOAlCaRecoStream.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 66 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(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 ( )

Definition at line 83 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
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 98 of file DQMHOAlCaRecoStream.cc.

References abs, relval_parameters_module::energy, MonitorElement::Fill(), edm::Event::getByLabel(), hDirCosine, hHOTime, hMuonEta, hMuonMom, hMuonMultipl, hMuonPhi, hoCalibVariableCollectionTag, hPedRing, hSignal3x3, hSigRing, gen::k, create_public_lumi_plots::log, LogDebug, m_highRadPosInMuch, Nevents, Nmuons, funct::pow(), hitfit::return, and funct::tan().

 {
   using namespace edm;
   
   Nevents++;
   
   edm::Handle<HOCalibVariableCollection>HOCalib;
   bool isCosMu = true;
   
   iEvent.getByLabel(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).size() >0 ) {
       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::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 163 of file DQMHOAlCaRecoStream.cc.

References DQMStore::book1D(), dbe_, folderName_, hDirCosine, hHOTime, hMuonEta, hMuonMom, hMuonMultipl, hMuonPhi, hPedRing, hSignal3x3, hSigRing, i, j, gen::k, m_highEdge, m_lowEdge, m_nbins, mergeVDriftHistosByStation::name, Nevents, Nmuons, cppFunctionSkipper::operator, MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), and indexGen::title.

 {
   dbe_ = edm::Service<DQMStore>().operator->();
   dbe_->setCurrentFolder(folderName_);
 
   char title[200];
   char name[200];
 
   hMuonMom = dbe_->book1D("hMuonMom", "Muon momentum (GeV)", 50, -100, 100);
   hMuonMom ->setAxisTitle("Muon momentum (GeV)",1);
   
   hMuonEta = dbe_->book1D("hMuonEta", "Pseudo-rapidity of muon", 50, -1.5, 1.5);
   hMuonEta ->setAxisTitle("Pseudo-rapidity of muon",1);
   
   hMuonPhi = dbe_->book1D("hMuonPhi", "Azimuthal angle of muon", 24, -acos(-1), acos(-1));
   hMuonPhi ->setAxisTitle("Azimuthal angle of muon",1);
     
   hMuonMultipl = dbe_->book1D("hMuonMultipl", "Muon Multiplicity", 10, 0.5, 10.5); 
   hMuonMultipl ->setAxisTitle("Muon Multiplicity",1);
 
   hDirCosine = dbe_->book1D("hDirCosine", "Direction Cosine of muon at HO tower", 50, -1., 1.);
   hDirCosine ->setAxisTitle("Direction Cosine of muon at HO tower",1);
 
   hHOTime = dbe_->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] = dbe_->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] = dbe_->book1D(name, title, m_nbins, m_lowEdge, m_highEdge);
     hPedRing[i]->setAxisTitle(title,1);
   }
 
   //  hSigRingm1 = dbe_->book1D("hSigRingm1", "HO signal in Ring-1", m_nbins, m_lowEdge, m_highEdge);
   //  hSigRingm1->setAxisTitle("HO signal in Ring-1",1);
 
   //  hSigRing00 = dbe_->book1D("hSigRing00", "HO signal in Ring_0", m_nbins, m_lowEdge, m_highEdge);
   //  hSigRing00->setAxisTitle("HO signal in Ring_0",1);
 
   //  hSigRingp1 = dbe_->book1D("hSigRingp1", "HO signal in Ring-1", m_nbins, m_lowEdge, m_highEdge);
   //  hSigRingp1->setAxisTitle("HO signal in Ring+1",1);
 
   //  hSigRingp2 = dbe_->book1D("hSigRingp2", "HO signal in Ring-2", m_nbins, m_lowEdge, m_highEdge);
   //  hSigRingp2->setAxisTitle("HO signal in Ring+2",1);
   
   //  hPedRingm2 = dbe_->book1D("hPedRingm2", "HO pedestal in Ring-2", m_nbins, m_lowEdge, m_highEdge);
   //  hPedRingm1 = dbe_->book1D("hPedRingm1", "HO pedestal in Ring-1", m_nbins, m_lowEdge, m_highEdge);
   //  hPedRing00 = dbe_->book1D("hPedRing00", "HO pedestal in Ring_0", m_nbins, m_lowEdge, m_highEdge);
   //  hPedRingp1 = dbe_->book1D("hPedRingp1", "HO pedestal in Ring-1", m_nbins, m_lowEdge, m_highEdge);
   //  hPedRingp2 = dbe_->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] = dbe_->book1D(title, title, m_nbins, m_lowEdge, m_highEdge);
       hSignal3x3[k]->setAxisTitle(title,1);
     }
   }
 
   Nevents = 0;
   Nmuons = 0;
 
 }

void DQMHOAlCaRecoStream::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 236 of file DQMHOAlCaRecoStream.cc.

References dbe_, MonitorElement::getTH1F(), hDirCosine, hHOTime, hMuonEta, hMuonMom, hMuonPhi, hPedRing, hSignal3x3, hSigRing, gen::k, max(), Nmuons, DQMStore::save(), saveToFile_, pileupReCalc_HLTpaths::scale, and theRootFileName.

                             {
   if (saveToFile_) {
 
     //    double scale = 1./max(1,Nevents);
     double scale = 1./max(1,Nmuons);
     hMuonMom->getTH1F()->Scale(scale);
     hMuonEta->getTH1F()->Scale(scale);
     hMuonPhi->getTH1F()->Scale(scale);
     hDirCosine->getTH1F()->Scale(scale);
     hHOTime->getTH1F()->Scale(scale);
     
     //    scale = 1./max(1,Nmuons);
     for (int k=0; k<5; k++) {
       hSigRing[k]->getTH1F()->Scale(scale);
       hPedRing[k]->getTH1F()->Scale(scale);
     }
     
     for (int k=0; k<9; k++) {
       hSignal3x3[k]->getTH1F()->Scale(scale);
     }
 
     dbe_->save(theRootFileName); 
   }
 
 }