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

Inheritance diagram for DQMHOAlCaRecoStream:

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

	~DQMHOAlCaRecoStream () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

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
	)

overrideprivate

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