#include <EwkDQM.h>

Inheritance diagram for EwkDQM:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
	Get the analysis. More...

void	beginJob ()
	Inizialize parameters for histo binning. More...

void	beginRun (const edm::Run &, const edm::EventSetup &)

double	calcDeltaPhi (double phi1, double phi2)

void	endJob (void)
	Save the histos. More...

	EwkDQM (const edm::ParameterSet &)
	Constructor. More...

virtual	~EwkDQM ()
	Destructor. More...

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 Attributes
double	eJetMin_

MonitorElement *	h_e1_et

MonitorElement *	h_e1_eta

MonitorElement *	h_e1_phi

MonitorElement *	h_e2_et

MonitorElement *	h_e2_eta

MonitorElement *	h_e2_phi

MonitorElement *	h_e_invWMass

MonitorElement *	h_ee_invMass

MonitorElement *	h_jet2_et

MonitorElement *	h_jet2_eta

MonitorElement *	h_jet2_phi

MonitorElement *	h_jet_count

MonitorElement *	h_jet_et

MonitorElement *	h_jet_eta

MonitorElement *	h_jet_phi

MonitorElement *	h_jet_pt

MonitorElement *	h_m1_eta

MonitorElement *	h_m1_phi

MonitorElement *	h_m1_pt

MonitorElement *	h_m2_eta

MonitorElement *	h_m2_phi

MonitorElement *	h_m2_pt

MonitorElement *	h_m_invWMass

MonitorElement *	h_met

MonitorElement *	h_met_phi

MonitorElement *	h_mumu_invMass

MonitorElement *	h_vertex_chi2

MonitorElement *	h_vertex_d0

MonitorElement *	h_vertex_number

MonitorElement *	h_vertex_numTrks

MonitorElement *	h_vertex_sumTrks

HLTConfigProvider	hltConfigProvider_

bool	isValidHltConfig_

std::string	logTraceName

edm::InputTag	theCaloMETCollectionLabel_

DQMStore *	theDbe

std::vector< std::string >	theElecTriggerPathToPass_

edm::InputTag	theElectronCollectionLabel_

edm::InputTag	theMuonCollectionLabel_

std::vector< std::string >	theMuonTriggerPathToPass_

edm::InputTag	thePFJetCollectionLabel_

edm::InputTag	theTriggerResultsCollection_

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

DQM offline for SMP V+Jets

Date:: 2012/06/28 10:29:38

Revision:: 1.12

Author: Valentina Gori, University of Firenze

Definition at line 27 of file EwkDQM.h.

Constructor & Destructor Documentation

EwkDQM::EwkDQM ( const edm::ParameterSet & parameters )

Constructor.

Definition at line 51 of file EwkDQM.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and cppFunctionSkipper::operator.

                                              {
   eJetMin_     = parameters.getUntrackedParameter<double>("EJetMin", 999999.);
 
   // riguardare questa sintassi
   // Get parameters from configuration file
   theElecTriggerPathToPass_ =
     parameters.getParameter<std::vector<string> >("elecTriggerPathToPass");
   theMuonTriggerPathToPass_ =
     parameters.getParameter<std::vector<string> >("muonTriggerPathToPass");
   //   eleTrigPathNames_ =
   //     parameters.getUntrackedParameter<std::vector<std::string> >("eleTrigPathNames");
   //   muTrigPathNames_ =
   //     parameters.getUntrackedParameter< std::vector<std::string> >("muTrigPathNames");
   theTriggerResultsCollection_ =
     parameters.getParameter<InputTag>("triggerResultsCollection");
   theMuonCollectionLabel_ = parameters.getParameter<InputTag>("muonCollection");
   theElectronCollectionLabel_ =
     parameters.getParameter<InputTag>("electronCollection");
   //  theCaloJetCollectionLabel_   = parameters.getParameter<InputTag>("caloJetCollection");
   thePFJetCollectionLabel_ =
     parameters.getParameter<InputTag>("PFJetCollection");
   theCaloMETCollectionLabel_ = parameters.getParameter<InputTag>("caloMETCollection");
 
   // just to initialize
   isValidHltConfig_ = false;
 
   h_vertex_number = 0;
   h_vertex_chi2 = 0;
   h_vertex_numTrks = 0;
   h_vertex_sumTrks = 0;
   h_vertex_d0 = 0;
 
   h_jet_count = 0;
   h_jet_et = 0;
   h_jet_pt = 0;
   h_jet_eta = 0;
   h_jet_phi = 0;
   h_jet2_et = 0;
   // h_jet2_pt = 0;
   h_jet2_eta = 0;
   h_jet2_phi = 0;
 
   h_e1_et = 0;
   h_e2_et = 0;
   h_e1_eta = 0;
   h_e2_eta = 0;
   h_e1_phi = 0;
   h_e2_phi = 0;
 
   h_m1_pt = 0;
   h_m2_pt = 0;
   h_m1_eta = 0;
   h_m2_eta = 0;
   h_m1_phi = 0;
   h_m2_phi = 0;
 
   // h_t1_et = 0;
   // h_t1_eta = 0;
   // h_t1_phi = 0;
 
   h_met = 0;
   h_met_phi = 0;
 
   h_e_invWMass = 0;
   h_m_invWMass = 0;
   h_mumu_invMass = 0;
   h_ee_invMass = 0;
 
   theDbe = Service<DQMStore>().operator->();
 }

EwkDQM::~EwkDQM ( )

virtual

Destructor.

Definition at line 122 of file EwkDQM.cc.

122 {

123 }

Member Function Documentation

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

virtual

Get the analysis.

Implements edm::EDAnalyzer.

Definition at line 247 of file EwkDQM.cc.

References calcDeltaPhi(), funct::cos(), reco::e1, reco::e2, reco::LeafCandidate::et(), reco::LeafCandidate::eta(), edm::Event::getByLabel(), i, getHLTprescales::index, edm::HandleBase::isValid(), j, metsig::jet, LogTrace, RPCpg::mu, reco::LeafCandidate::phi(), Geom::pi(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), trigNames, findQualityFiles::v, and GoodVertex_cfg::vertexCollection.

                                                                   {
   // short-circuit if hlt problems
   if (!isValidHltConfig_)
     return;
 
   LogTrace(logTraceName) << "Analysis of event # ";
   // Did it pass certain HLT path?
   Handle<TriggerResults> HLTresults;
   iEvent.getByLabel(theTriggerResultsCollection_, HLTresults);
   if (!HLTresults.isValid())
     return;
 
   const edm::TriggerNames & trigNames = iEvent.triggerNames(*HLTresults);
 
   // a temporary, until we have a list of triggers of interest
   std::vector<std::string> eleTrigPathNames;
   std::vector<std::string> muTrigPathNames;
 
   // eleTrigPathNames.push_back(theElecTriggerPathToPass_);
   // muTrigPathNames.push_back(theMuonTriggerPathToPass_);
   // end of temporary
 
   bool passed_electron_HLT = false;
   bool passed_muon_HLT     = false;
   for (unsigned int i = 0; i < HLTresults->size(); i++) {
     const std::string trigName = trigNames.triggerName(i);
     // check if triggerName matches electronPath
     for (unsigned int index = 0;
          index < theElecTriggerPathToPass_.size() && !passed_electron_HLT;
          index++) {
       // 0 if found, pos if not
       size_t trigPath = trigName.find(theElecTriggerPathToPass_[index]);
       if (trigPath == 0) {
         //      cout << "MuonTrigger passed (=trigName): " << trigName <<endl;
         passed_electron_HLT = HLTresults->accept(i);
       }
     }
     // check if triggerName matches muonPath
     for (unsigned int index = 0;
          index < theMuonTriggerPathToPass_.size() && !passed_muon_HLT;
          index++) {
       // 0 if found, pos if not
       size_t trigPath = trigName.find(theMuonTriggerPathToPass_[index]);
       if (trigPath == 0) {
         //      cout << "MuonTrigger passed (=trigName): " << trigName <<endl;
         passed_muon_HLT = HLTresults->accept(i);
       }
     }
   }
 
   // we are interested in events with a valid electron or muon
   if (!(passed_electron_HLT || passed_muon_HLT))
     return;
 
   // Vertex information
   Handle<VertexCollection> vertexHandle;
   iEvent.getByLabel("offlinePrimaryVertices", vertexHandle);
   if (!vertexHandle.isValid())
     return;
   VertexCollection vertexCollection = *(vertexHandle.product());
   VertexCollection::const_iterator v = vertexCollection.begin();
   int vertex_number     = vertexCollection.size();
   double vertex_chi2    = v->normalizedChi2();  // v->chi2();
   double vertex_d0      = sqrt(v->x()*v->x()+v->y()*v->y());
   double vertex_numTrks = v->tracksSize();
   double vertex_sumTrks = 0.0;
   // std::cout << "vertex_d0=" << vertex_d0 << "\n";
   // double vertex_ndof    = v->ndof();cout << "ndof="<<vertex_ndof<<endl;
   for (Vertex::trackRef_iterator vertex_curTrack = v->tracks_begin();
        vertex_curTrack != v->tracks_end(); vertex_curTrack++)
     vertex_sumTrks += (*vertex_curTrack)->pt();
 
   // Missing ET
   Handle< View<MET> > caloMETCollection;
   iEvent.getByLabel(theCaloMETCollectionLabel_, caloMETCollection);
   if (!caloMETCollection.isValid())
     return;
   float missing_et = caloMETCollection->begin()->et();
   float met_phi    = caloMETCollection->begin()->phi();
 
 
   // grab "gaussian sum fitting" electrons
   Handle<GsfElectronCollection> electronCollection;
   iEvent.getByLabel(theElectronCollectionLabel_, electronCollection);
   if (!electronCollection.isValid())
     return;
 
   // Find the highest and 2nd highest electron
   float electron_et   = -8.0;
   float electron_eta  = -8.0;
   float electron_phi  = -8.0;
   float electron2_et  = -9.0;
   float electron2_eta = -9.0;
   float electron2_phi = -9.0;
   float ee_invMass    = -9.0;
   TLorentzVector e1, e2;
 
   // If it passed electron HLT and the collection was found, find electrons near Z mass
   if (passed_electron_HLT) {
     for (reco::GsfElectronCollection::const_iterator recoElectron = electronCollection->begin();
          recoElectron != electronCollection->end(); recoElectron++) {
       // Require electron to pass some basic cuts
       if (recoElectron->et() < 20 || fabs(recoElectron->eta()) > 2.5)
         continue;
 
       // Tighter electron cuts
       if (recoElectron->deltaPhiSuperClusterTrackAtVtx() > 0.58 ||
           recoElectron->deltaEtaSuperClusterTrackAtVtx() > 0.01 ||
           recoElectron->sigmaIetaIeta() > 0.027)
         continue;
 
       if (recoElectron->et() > electron_et) {
         electron2_et  = electron_et;  // 2nd highest gets values from current highest
         electron2_eta = electron_eta;
         electron2_phi = electron_phi;
         electron_et   = recoElectron->et();  // 1st highest gets values from new highest
         electron_eta  = recoElectron->eta();
         electron_phi  = recoElectron->phi();
         e1 = TLorentzVector(recoElectron->momentum().x(),
                             recoElectron->momentum().y(),
                             recoElectron->momentum().z(),
                             recoElectron->p());
       } else if (recoElectron->et() > electron2_et) {
         electron2_et  = recoElectron->et();
         electron2_eta = recoElectron->eta();
         electron2_phi = recoElectron->phi();
         e2 = TLorentzVector(recoElectron->momentum().x(),
                             recoElectron->momentum().y(),
                             recoElectron->momentum().z(),
                             recoElectron->p());
       }
     }  // end of loop over electrons
     if (electron2_et > 0.0) {
       TLorentzVector pair = e1+e2;
       ee_invMass = pair.M();
     }
   }  // end of "are electrons valid"
 
 
 
   // Take the STA muon container
   Handle<MuonCollection> muonCollection;
   iEvent.getByLabel(theMuonCollectionLabel_, muonCollection);
   if (!muonCollection.isValid())
     return;
 
   // Find the highest pt muons
   float mm_invMass = -9.0;
   float muon_pt    = -9.0;
   float muon_eta   = -9.0;
   float muon_phi   = -9.0;
   float muon2_pt   = -9.0;
   float muon2_eta  = -9.0;
   float muon2_phi  = -9.0;
   TLorentzVector m1, m2;
 
   if (passed_muon_HLT) {
     for (reco::MuonCollection::const_iterator recoMuon = muonCollection->begin();
          recoMuon != muonCollection->end(); recoMuon++) {
       // Require muon to pass some basic cuts
       if (recoMuon->pt() < 20 || !recoMuon->isGlobalMuon())
         continue;
       // Some tighter muon cuts
       if (recoMuon->globalTrack()->normalizedChi2() > 10)
         continue;
 
       if (recoMuon->pt() > muon_pt) {
         muon2_pt  = muon_pt;  // 2nd highest gets values from current highest
         muon2_eta = muon_eta;
         muon2_phi = muon_phi;
         muon_pt   = recoMuon->pt();  // 1st highest gets values from new highest
         muon_eta  = recoMuon->eta();
         muon_phi  = recoMuon->phi();
         m1 = TLorentzVector(recoMuon->momentum().x(),
                             recoMuon->momentum().y(),
                             recoMuon->momentum().z(),
                             recoMuon->p());
       } else if (recoMuon->pt() > muon2_pt) {
         muon2_pt  = recoMuon->pt();
         muon2_eta = recoMuon->eta();
         muon2_phi = recoMuon->phi();
         m2 = TLorentzVector(recoMuon->momentum().x(),
                             recoMuon->momentum().y(),
                             recoMuon->momentum().z(),
                             recoMuon->p());
       }
     }
   }
   if (muon2_pt > 0.0) {
     TLorentzVector pair = m1+m2;
     mm_invMass = pair.M();
   }
 
 
   // Find the highest et jet
 
   //  Handle<CaloJetCollection> caloJetCollection;
   Handle<View<Jet> > PFJetCollection;
   //  iEvent.getByLabel (theCaloJetCollectionLabel,caloJetCollection);
   iEvent.getByLabel(thePFJetCollectionLabel_, PFJetCollection);
   //  if ( !caloJetCollection.isValid() ) return;
   if (!PFJetCollection.isValid())
     return;
 
   unsigned int muonCollectionSize = muonCollection->size();
   // unsigned int jetCollectionSize = jetCollection->size();
   unsigned int PFJetCollectionSize = PFJetCollection->size();
   int jet_count = 0;
   // int LEADJET=-1;  double max_pt=0;
 
 
   float jet_et    = -80.0;
   float jet_pt    = -80.0;  // prova
   float jet_eta   = -80.0;  // now USED
   float jet_phi   = -80.0;  // now USED
   float jet2_et   = -90.0;
   float jet2_eta  = -90.0;  // now USED
   float jet2_phi  = -90.0;  // now USED
   //  for (CaloJetCollection::const_iterator i_calojet = caloJetCollection->begin();
   //       i_calojet != caloJetCollection->end(); i_calojet++) {
   //  for (PFJetCollection::const_iterator i_pfjet = PFJetCollection->begin();
   //       i_pfjet != PFJetCollection->end(); i_pfjet++) {
   //  float jet_current_et = i_calojet->et();
   //  float jet_current_et = i_pfjet->et();            // e` identico a jet.et()
   //    jet_count++;
 
   // cleaning: va messo prima del riempimento dell'istogramma // This is in order to use PFJets
   for (unsigned int i = 0; i < PFJetCollectionSize; i++) {
     const Jet& jet = PFJetCollection->at(i);
     // la classe "jet" viene definita qui!!!
     double minDistance = 99999;
     for (unsigned int j = 0; j < muonCollectionSize; j++) {
       const Muon& mu = muonCollection->at(j);
       double distance = sqrt((mu.eta() - jet.eta()) * (mu.eta() - jet.eta()) +
                              (mu.phi() - jet.phi()) * (mu.phi() - jet.phi()));
       if (minDistance > distance)
         minDistance = distance;
     }
     if (minDistance < 0.3)
       continue;  // 0.3 is the isolation cone around the muon
     // se la distanza muone-cono del jet e` minore di 0.3, passo avanti e non conteggio il mio jet
 
     // If it overlaps with ELECTRON, it is not a jet
     if (electron_et > 0.0 &&
         fabs(jet.eta() - electron_eta) < 0.2 &&
         calcDeltaPhi(jet.phi(), electron_phi) < 0.2)
       continue;
     if (electron2_et > 0.0 &&
         fabs(jet.eta() - electron2_eta) < 0.2 &&
         calcDeltaPhi(jet.phi(), electron2_phi) < 0.2)
       continue;
 
     // provo a cambiare la parte degli elettroni in modo simmetrico alla parte per i muoni
 
     // ...
     // ...
 
 
     // if it has too low Et, throw away
     if (jet.et() < eJetMin_)
       continue;
     jet_count++;
 
     // ovvero: incrementa jet_count se:
     //   - non c'e un muone entro 0.3 di distanza dal cono del jet;
     //   - se il jet non si sovrappone ad un elettrone;
     //   - se l'energia trasversa e` maggiore della soglia impostata (15?)
 
     // if(jet.et()>max_pt) { LEADJET=i; max_pt=jet.et();}
     // se l'energia del jet e` maggiore di max_pt, diventa "i"
     // l'indice del jet piu` energetico e max_pt la sua energia
 
     // riguardare questo!!!
     // fino ad ora, jet_et era inizializzato a -8.0
     if (jet.et() > jet_et) {
       jet2_et  = jet_et;  // 2nd highest jet gets et from current highest
       // perche` prende l'energia del primo jet??
       jet2_eta = jet_eta;  // now USED
       jet2_phi = jet_phi;  // now USED
       // jet_et = i_calojet->et(); // current highest jet gets
       // et from the new highest
       jet_et   = jet.et();  // current highest jet gets et from the new highest
       // ah, ok! lo riaggiorna solo dopo!
       jet_pt   = jet.pt();   // e` il pT del leading jet
       jet_eta  = jet.eta();  // now USED
       jet_phi  = jet.phi() * (Geom::pi() / 180.);  // now USED
     } else if (jet.et() > jet2_et) {
       //      jet2_et  = i_calojet->et();
       jet2_et  = jet.et();
       //      jet2_eta = i_calojet->eta();  // UNUSED
       //      jet2_phi = i_calojet->phi();  // UNUSED
       jet2_eta = jet.eta();  // now USED
       jet2_phi = jet.phi();  // now USED
     }
     // questo elseif funziona
   }
 
 
 
   //                 Fill Histograms                                            //
 
   bool fill_e1  = false;
   bool fill_e2  = false;
   bool fill_m1  = false;
   bool fill_m2  = false;
   bool fill_met = false;
 
   // Was Z->ee found?
   if (ee_invMass > 0.0) {
     h_ee_invMass->Fill(ee_invMass);
     fill_e1 = true;
     fill_e2 = true;
   }
 
   // Was Z->mu mu found?
   if (mm_invMass > 0.0) {
     h_mumu_invMass->Fill(mm_invMass);
     fill_m1 = true;
     fill_m2 = true;
     h_jet2_et ->Fill(jet2_et);
   }
 
   // Was W->e nu found?
   if (electron_et > 0.0 && missing_et > 20.0) {
     float dphiW  = fabs(met_phi-electron_phi);
     float W_mt_e = sqrt(2 * missing_et * electron_et * (1 - cos(dphiW)));
     h_e_invWMass->Fill(W_mt_e);
     fill_e1  = true;
     fill_met = true;
   }
 
   // Was W->mu nu found?
   if (muon_pt > 0.0 && missing_et > 20.0) {
     float dphiW = fabs(met_phi - muon_phi);
     float W_mt_m = sqrt(2 * missing_et * muon_pt * (1 - cos(dphiW)));
     h_m_invWMass->Fill(W_mt_m);
     fill_m1  = true;
     fill_met = true;
   }
 
   if (jet_et > -10.0) {
     h_jet_et->Fill(jet_et);
     h_jet_count->Fill(jet_count);
   }
 
   if (jet_pt > 0.) {
     h_jet_pt->Fill(jet_pt);
   }
 
   if (jet_eta > -50.) {
     h_jet_eta->Fill(jet_eta);
   }
 
   if (jet_phi > -10.) {
     h_jet_phi->Fill(jet_phi);
   }
 
   if (jet2_et > -10.0) {
     h_jet2_et->Fill(jet2_et);
   }
 
   // if (jet2_pt>0.) {
   //  h_jet2_pt   ->Fill(jet2_pt);
   // }
 
   if (jet2_eta > -50.) {
     h_jet2_eta->Fill(jet2_eta);
   }
 
   if (jet2_phi > -10.) {
     h_jet2_phi->Fill(jet2_phi);
   }
 
 
 
   if (fill_e1 || fill_m1) {
     h_vertex_number->Fill(vertex_number);
     h_vertex_chi2->Fill(vertex_chi2);
     h_vertex_d0->Fill(vertex_d0);
     h_vertex_numTrks->Fill(vertex_numTrks);
     h_vertex_sumTrks->Fill(vertex_sumTrks);
   }
 
   if (fill_e1) {
     h_e1_et->Fill(electron_et);
     h_e1_eta->Fill(electron_eta);
     h_e1_phi->Fill(electron_phi);
   }
   if (fill_e2) {
     h_e2_et->Fill(electron2_et);
     h_e2_eta->Fill(electron2_eta);
     h_e2_phi->Fill(electron2_phi);
   }
   if (fill_m1) {
     h_m1_pt->Fill(muon_pt);
     h_m1_eta->Fill(muon_eta);
     h_m1_phi->Fill(muon_phi);
   }
   if (fill_m2) {
     h_m2_pt->Fill(muon2_pt);
     h_m2_eta->Fill(muon2_eta);
     h_m2_phi->Fill(muon2_phi);
   }
   if (fill_met) {
     h_met->Fill(missing_et);
     h_met_phi->Fill(met_phi);
   }
 }

void EwkDQM::beginJob ( void )

virtual

Inizialize parameters for histo binning.

Reimplemented from edm::EDAnalyzer.

Definition at line 126 of file EwkDQM.cc.

References LogTrace, and pi.

                       {
   char chtitle[256] = "";
   const size_t title_s = sizeof(chtitle);
 
   logTraceName = "EwkAnalyzer";
 
   LogTrace(logTraceName) << "Parameters initialization";
   theDbe->setCurrentFolder("Physics/EwkDQM");  // Use folder with name of PAG
 
   const float pi = 4*atan(1);
 
   // Keep the number of plots and number of bins to a minimum!
   h_vertex_number = theDbe->book1D("vertex_number",
                                    "Number of event vertices in collection",
                                    10, -0.5, 9.5);
   h_vertex_chi2  = theDbe->book1D("vertex_chi2",
                                   "Event Vertex #chi^{2}/n.d.o.f.",
                                   20, 0.0, 2.0);
   h_vertex_numTrks = theDbe->book1D("vertex_numTrks",
                                     "Event Vertex, number of tracks",
                                     20, -0.5, 59.5);
   h_vertex_sumTrks = theDbe->book1D("vertex_sumTrks",
                                     "Event Vertex, sum of track pt",
                                     20, 0.0, 100.0);
   h_vertex_d0 = theDbe->book1D("vertex_d0",
                                "Event Vertex d0", 20, 0.0, 0.05);
 
   snprintf(chtitle, title_s, "Number of %s (E_{T} > 15 GeV);Number of Jets",
            thePFJetCollectionLabel_.label().data());
   h_jet_count = theDbe->book1D("jet_count", chtitle, 8, -0.5, 7.5);
 
   snprintf(chtitle, title_s, "Leading jet E_{T} (from %s);E_{T}(1^{st} jet) (GeV)",
            thePFJetCollectionLabel_.label().data());
   h_jet_et = theDbe->book1D("jet_et", chtitle,    20, 0., 200.0);
 
   snprintf(chtitle, title_s, "Leading jet p_{T} (from %s);p_{T}(1^{st} jet) (GeV/c)",
            thePFJetCollectionLabel_.label().data());
   h_jet_pt = theDbe->book1D("jet_pt", chtitle,  20, 0., 200.0);
 
   snprintf(chtitle, title_s,  "Leading jet #eta (from %s); #eta (1^{st} jet)",
            thePFJetCollectionLabel_.label().data());
   h_jet_eta = theDbe->book1D("jet_eta", chtitle, 20, -10., 10.0);
   snprintf(chtitle, title_s, "Leading jet #phi (from %s); #phi(1^{st} jet)",
            thePFJetCollectionLabel_.label().data());
   h_jet_phi = theDbe->book1D("jet_phi", chtitle, 22, -1.1*pi, 1.1*pi);
 
   snprintf(chtitle, title_s, "2^{nd} leading jet E_{T} (from %s);E_{T}(2^{nd} jet) (GeV)",
            thePFJetCollectionLabel_.label().data());
   h_jet2_et = theDbe->book1D("jet2_et", chtitle,  20, 0., 200.0);
   // snprintf(chtitle, title_s, "2^{nd} leading jet p_{T} (from %s);p_{T}(2^{nd} jet) (GeV/c)",
   // thePFJetCollectionLabel_.label().data());
   // h_jet2_pt = theDbe->book1D("jet2_pt", chtitle,  20, 0., 200.0);
 
   snprintf(chtitle, title_s,  "2^{nd} leading jet #eta (from %s); #eta (2^{nd} jet)",
            thePFJetCollectionLabel_.label().data());
   h_jet2_eta = theDbe->book1D("jet2_eta", chtitle,  20, -10., 10.0);
 
   snprintf(chtitle, title_s, "2^{nd} leading jet #phi (from %s); #phi(2^{nd} jet)",
            thePFJetCollectionLabel_.label().data());
   h_jet2_phi = theDbe->book1D("jet2_phi", chtitle,  22, -1.1*pi, 1.1*pi);
 
   h_e1_et  = theDbe->book1D("e1_et",  "E_{T} of Leading Electron;E_{T} (GeV)",
                             20, 0.0, 100.0);
   h_e2_et  = theDbe->book1D("e2_et",  "E_{T} of Second Electron;E_{T} (GeV)",
                             20, 0.0, 100.0);
   h_e1_eta = theDbe->book1D("e1_eta", "#eta of Leading Electron;#eta",
                             20, -4.0, 4.0);
   h_e2_eta = theDbe->book1D("e2_eta", "#eta of Second Electron;#eta",
                             20, -4.0, 4.0);
   h_e1_phi = theDbe->book1D("e1_phi", "#phi of Leading Electron;#phi",
                             22, -1.1*pi, 1.1*pi);
   h_e2_phi = theDbe->book1D("e2_phi", "#phi of Second Electron;#phi",
                             22, -1.1*pi, 1.1*pi);
   h_m1_pt  = theDbe->book1D("m1_pt",  "p_{T} of Leading Muon;p_{T}(1^{st} #mu) (GeV)",
                             20, 0.0, 100.0);
   h_m2_pt  = theDbe->book1D("m2_pt",  "p_{T} of Second Muon;p_{T}(2^{nd} #mu) (GeV)",
                             20, 0.0, 100.0);
   h_m1_eta = theDbe->book1D("m1_eta", "#eta of Leading Muon;#eta(1^{st} #mu)",
                             20, -4.0, 4.0);
   h_m2_eta = theDbe->book1D("m2_eta", "#eta of Second Muon;#eta(2^{nd} #mu)",
                             20, -4.0, 4.0);
   h_m1_phi = theDbe->book1D("m1_phi", "#phi of Leading Muon;#phi(1^{st} #mu)",
                             20, (-1. - 1./10.)*pi, (1. + 1./10.)*pi);
   h_m2_phi = theDbe->book1D("m2_phi", "#phi of Second Muon;#phi(2^{nd} #mu)",
                             20, (-1. - 1./10.)*pi, (1. + 1./10.)*pi);
   //  h_t1_et = theDbe->book1D("t1_et", "E_{T} of Leading Tau;E_{T} (GeV)",
   //  20, 0.0 , 100.0);
   //  h_t1_eta = theDbe->book1D("t1_eta", "#eta of Leading Tau;#eta",
   //  20, -4.0, 4.0);
   //  h_t1_phi = theDbe->book1D("t1_phi", "#phi of Leading Tau;#phi",
   //  20, -4.0, 4.0);
   snprintf(chtitle, title_s, "Missing E_{T} (%s); GeV",
            theCaloMETCollectionLabel_.label().data());
   h_met          = theDbe->book1D("met", chtitle, 20, 0.0, 100);
   h_met_phi      = theDbe->book1D("met_phi", "Missing E_{T} #phi;#phi(MET)",
                                   22, (-1. - 1./10.)*pi, (1. + 1./10.)*pi);
 
   h_e_invWMass   = theDbe->book1D("we_invWMass", "W-> e #nu Transverse Mass;M_{T} (GeV)",
                                   20, 0.0, 140.0);
   h_m_invWMass   = theDbe->book1D("wm_invWMass", "W-> #mu #nu Transverse Mass;M_{T} (GeV)",
                                   20, 0.0, 140.0);
   h_mumu_invMass = theDbe->book1D("z_mm_invMass", "#mu#mu Invariant Mass;InvMass (GeV)",
                                   20, 40.0, 140.0);
   h_ee_invMass   = theDbe->book1D("z_ee_invMass", "ee Invariant Mass;InvMass (Gev)",
                                   20, 40.0, 140.0);
 }

void EwkDQM::beginRun	(	const edm::Run &	theRun,
		const edm::EventSetup &	theSetup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 236 of file EwkDQM.cc.

References AlCaHLTBitMon_QueryRunRegistry::string.

                                                                          {
   // passed as parameter to HLTConfigProvider::init(), not yet used
   bool isConfigChanged = false;
 
   // isValidHltConfig_ used to short-circuit analyze() in case of problems
   const std::string hltProcessName(theTriggerResultsCollection_.process());
   isValidHltConfig_ = hltConfigProvider_.init(theRun, theSetup,
                                               hltProcessName, isConfigChanged);
 }

double EwkDQM::calcDeltaPhi	(	double	phi1,
		double	phi2
	)

Definition at line 672 of file EwkDQM.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi.

                                                     {
   double deltaPhi = phi1 - phi2;
 
   if (deltaPhi < 0)
     deltaPhi = -deltaPhi;
 
   if (deltaPhi > 3.1415926)
     deltaPhi = 2 * 3.1415926 - deltaPhi;
 
   return deltaPhi;
 }

void EwkDQM::endJob ( void )

virtual

Save the histos.

Reimplemented from edm::EDAnalyzer.

Definition at line 668 of file EwkDQM.cc.

668 {}

Member Data Documentation

double EwkDQM::eJetMin_

private

Definition at line 74 of file EwkDQM.h.

MonitorElement* EwkDQM::h_e1_et

private

Definition at line 94 of file EwkDQM.h.

MonitorElement* EwkDQM::h_e1_eta

private

Definition at line 96 of file EwkDQM.h.

MonitorElement* EwkDQM::h_e1_phi

private

Definition at line 98 of file EwkDQM.h.

MonitorElement* EwkDQM::h_e2_et

private

Definition at line 95 of file EwkDQM.h.

MonitorElement* EwkDQM::h_e2_eta

private

Definition at line 97 of file EwkDQM.h.

MonitorElement* EwkDQM::h_e2_phi

private

Definition at line 99 of file EwkDQM.h.

MonitorElement* EwkDQM::h_e_invWMass

private

Definition at line 115 of file EwkDQM.h.

MonitorElement* EwkDQM::h_ee_invMass

private

Definition at line 118 of file EwkDQM.h.

MonitorElement* EwkDQM::h_jet2_et

private

Definition at line 89 of file EwkDQM.h.

MonitorElement* EwkDQM::h_jet2_eta

private

Definition at line 91 of file EwkDQM.h.

MonitorElement* EwkDQM::h_jet2_phi

private

Definition at line 92 of file EwkDQM.h.

MonitorElement* EwkDQM::h_jet_count

private

Definition at line 83 of file EwkDQM.h.

MonitorElement* EwkDQM::h_jet_et

private

Definition at line 84 of file EwkDQM.h.

MonitorElement* EwkDQM::h_jet_eta

private

Definition at line 86 of file EwkDQM.h.

MonitorElement* EwkDQM::h_jet_phi

private

Definition at line 87 of file EwkDQM.h.

MonitorElement* EwkDQM::h_jet_pt

private

Definition at line 85 of file EwkDQM.h.

MonitorElement* EwkDQM::h_m1_eta

private

Definition at line 103 of file EwkDQM.h.

MonitorElement* EwkDQM::h_m1_phi

private

Definition at line 105 of file EwkDQM.h.

MonitorElement* EwkDQM::h_m1_pt

private

Definition at line 101 of file EwkDQM.h.

MonitorElement* EwkDQM::h_m2_eta

private

Definition at line 104 of file EwkDQM.h.

MonitorElement* EwkDQM::h_m2_phi

private

Definition at line 106 of file EwkDQM.h.

MonitorElement* EwkDQM::h_m2_pt

private

Definition at line 102 of file EwkDQM.h.

MonitorElement* EwkDQM::h_m_invWMass

private

Definition at line 116 of file EwkDQM.h.

MonitorElement* EwkDQM::h_met

private

Definition at line 112 of file EwkDQM.h.

MonitorElement* EwkDQM::h_met_phi

private

Definition at line 113 of file EwkDQM.h.

MonitorElement* EwkDQM::h_mumu_invMass

private

Definition at line 117 of file EwkDQM.h.

MonitorElement* EwkDQM::h_vertex_chi2

private

Definition at line 78 of file EwkDQM.h.

MonitorElement* EwkDQM::h_vertex_d0

private

Definition at line 81 of file EwkDQM.h.

MonitorElement* EwkDQM::h_vertex_number

private

Definition at line 77 of file EwkDQM.h.

MonitorElement* EwkDQM::h_vertex_numTrks

private

Definition at line 79 of file EwkDQM.h.

MonitorElement* EwkDQM::h_vertex_sumTrks

private

Definition at line 80 of file EwkDQM.h.

HLTConfigProvider EwkDQM::hltConfigProvider_

private

Definition at line 58 of file EwkDQM.h.

bool EwkDQM::isValidHltConfig_

private

Definition at line 59 of file EwkDQM.h.

std::string EwkDQM::logTraceName

private

Definition at line 56 of file EwkDQM.h.

edm::InputTag EwkDQM::theCaloMETCollectionLabel_

private

Definition at line 72 of file EwkDQM.h.

DQMStore* EwkDQM::theDbe

private

Definition at line 54 of file EwkDQM.h.

std::vector<std::string> EwkDQM::theElecTriggerPathToPass_

private

Definition at line 63 of file EwkDQM.h.

edm::InputTag EwkDQM::theElectronCollectionLabel_

private

Definition at line 69 of file EwkDQM.h.

edm::InputTag EwkDQM::theMuonCollectionLabel_

private

Definition at line 68 of file EwkDQM.h.

std::vector<std::string> EwkDQM::theMuonTriggerPathToPass_

private

Definition at line 64 of file EwkDQM.h.

edm::InputTag EwkDQM::thePFJetCollectionLabel_

private

Definition at line 71 of file EwkDQM.h.

edm::InputTag EwkDQM::theTriggerResultsCollection_

private

Definition at line 67 of file EwkDQM.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation