#include <QcdPhotonsDQM.h>

Inheritance diagram for QcdPhotonsDQM:

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 &)

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

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

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

virtual	~QcdPhotonsDQM ()
	Destructor. More...

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

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Attributes
MonitorElement *	h_deltaEt_photon_jet

MonitorElement *	h_deltaPhi_jet_jet2

MonitorElement *	h_deltaPhi_photon_jet

MonitorElement *	h_deltaPhi_photon_jet2

MonitorElement *	h_deltaR_jet_jet2

MonitorElement *	h_deltaR_photon_jet2

MonitorElement *	h_jet2_eta

MonitorElement *	h_jet2_pt

MonitorElement *	h_jet2_ptOverPhotonEt

MonitorElement *	h_jet_count

MonitorElement *	h_jet_eta

MonitorElement *	h_jet_pt

MonitorElement *	h_photon_count_bar

MonitorElement *	h_photon_count_end

MonitorElement *	h_photon_et

MonitorElement *	h_photon_et_jetco

MonitorElement *	h_photon_et_jetcs

MonitorElement *	h_photon_et_jetfo

MonitorElement *	h_photon_et_jetfs

MonitorElement *	h_photon_et_ratio_co_cs

MonitorElement *	h_photon_et_ratio_co_fo

MonitorElement *	h_photon_et_ratio_co_fs

MonitorElement *	h_photon_et_ratio_cs_fo

MonitorElement *	h_photon_et_ratio_cs_fs

MonitorElement *	h_photon_et_ratio_fo_fs

MonitorElement *	h_photon_eta

MonitorElement *	h_triggers_passed

HLTConfigProvider	hltConfigProvider_

bool	isValidHltConfig_

std::string	logTraceName

int	num_events_in_run

edm::InputTag	theCaloJetCollectionLabel

DQMStore *	theDbe

std::string	theHltMenu

double	theMinCaloJetPt

double	theMinPhotonEt

edm::InputTag	thePhotonCollectionLabel

double	thePlotJetMaxEta

double	thePlotPhotonMaxEt

double	thePlotPhotonMaxEta

std::vector< std::string >	thePlotTheseTriggersToo

bool	theRequirePhotonFound

std::string	theTriggerPathToPass

std::string	theTriggerResultsCollection

edm::InputTag	theVertexCollectionLabel

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)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

DQM offline for QCD-Photons

Date:: 2010/06/16 15:53:52

Revision:: 1.15

Author: Michael B. Anderson, University of Wisconsin Madison

Definition at line 27 of file QcdPhotonsDQM.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 65 of file QcdPhotonsDQM.cc.

References edm::ParameterSet::getParameter().

                                                            {
   // Get parameters from configuration file
   theTriggerPathToPass        = parameters.getParameter<string>("triggerPathToPass");
   thePlotTheseTriggersToo     = parameters.getParameter<vector<string> >("plotTheseTriggersToo");
   theHltMenu                  = parameters.getParameter<string>("hltMenu");
   theTriggerResultsCollection = parameters.getParameter<string>("triggerResultsCollection");
   thePhotonCollectionLabel    = parameters.getParameter<InputTag>("photonCollection");
   theCaloJetCollectionLabel   = parameters.getParameter<InputTag>("caloJetCollection");
   theVertexCollectionLabel    = parameters.getParameter<InputTag>("vertexCollection");
   theMinCaloJetPt             = parameters.getParameter<double>("minCaloJetPt");
   theMinPhotonEt              = parameters.getParameter<double>("minPhotonEt");
   theRequirePhotonFound       = parameters.getParameter<bool>("requirePhotonFound");
   thePlotPhotonMaxEt          = parameters.getParameter<double>("plotPhotonMaxEt");
   thePlotPhotonMaxEta         = parameters.getParameter<double>("plotPhotonMaxEta");
   thePlotJetMaxEta            = parameters.getParameter<double>("plotJetMaxEta");
   // just to initialize
   isValidHltConfig_ = false;
 
 }

QcdPhotonsDQM::~QcdPhotonsDQM ( )

virtual

Destructor.

Definition at line 85 of file QcdPhotonsDQM.cc.

85 {

86 }

Member Function Documentation

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

virtual

Get the analysis.

Implements edm::EDAnalyzer.

Definition at line 171 of file QcdPhotonsDQM.cc.

References abs, Geom::deltaPhi(), deltaR(), edm::SortedCollection< T, SORT >::find(), flags, edm::EventSetup::get(), edm::Event::getByLabel(), EcalClusterLazyTools::getMaximum(), i, edm::HandleBase::isValid(), LogTrace, edm::Handle< T >::product(), findQualityFiles::size, cond::rpcobgas::time, and GoodVertex_cfg::vertexCollection.

                                                                          {
   num_events_in_run++;
 
   // short-circuit if hlt problems
   if( ! isValidHltConfig_ ) return;
   
   LogTrace(logTraceName)<<"Analysis of event # ";
 
   // Did event pass HLT paths?
   Handle<TriggerResults> HLTresults;
   iEvent.getByLabel(InputTag(theTriggerResultsCollection, "", theHltMenu), HLTresults); 
   if (!HLTresults.isValid()) return;
 
   unsigned int triggerIndex; // index of trigger path
   bool passed_HLT;
 
   // See if event passed trigger paths
   //  increment that bin in the trigger plot
   for (unsigned int i=0; i<thePlotTheseTriggersToo.size(); i++) {
     passed_HLT = false;
     triggerIndex = hltConfigProvider_.triggerIndex(thePlotTheseTriggersToo[i]);
     if (triggerIndex < HLTresults->size()) passed_HLT = HLTresults->accept(triggerIndex);
     if (passed_HLT) h_triggers_passed->Fill(i);
   }
 
   // Quit if the event did not pass the HLT path we care about
   passed_HLT = false;
   triggerIndex = hltConfigProvider_.triggerIndex(theTriggerPathToPass); // index of trigger path
   if (triggerIndex < HLTresults->size()) passed_HLT = HLTresults->accept(triggerIndex);
   if (!passed_HLT) return;
 
 
   // Does event have valid vertex?
   // Get the primary event vertex
   Handle<VertexCollection> vertexHandle;
   iEvent.getByLabel(theVertexCollectionLabel, vertexHandle);
   VertexCollection vertexCollection = *(vertexHandle.product());
   double vtx_ndof = -1.0;
   double vtx_z    = 0.0;
   bool   vtx_isFake = true;
   if (vertexCollection.size()>0) {
     vtx_ndof = vertexCollection.begin()->ndof();
     vtx_z    = vertexCollection.begin()->z();
     vtx_isFake = false;
   }
   if (vtx_isFake || fabs(vtx_z)>15 || vtx_ndof<4) return;
 
 
   // Did the event pass certain L1 Technical Trigger bits?
   // It's probably beam halo
   //  TODO: ADD code
 
   // grab photons
   Handle<PhotonCollection> photonCollection;
   iEvent.getByLabel(thePhotonCollectionLabel, photonCollection);
 
   // If photon collection is empty, exit
   if (!photonCollection.isValid()) return;
 
   // For finding spikes
   Handle<EcalRecHitCollection> EBReducedRecHits;
   iEvent.getByLabel("reducedEcalRecHitsEB", EBReducedRecHits);
   Handle<EcalRecHitCollection> EEReducedRecHits;
   iEvent.getByLabel("reducedEcalRecHitsEE", EEReducedRecHits); 
   EcalClusterLazyTools lazyTool(iEvent, iSetup, InputTag("reducedEcalRecHitsEB"), InputTag("reducedEcalRecHitsEE") );
   
 
   // Find the highest et "decent" photon
   float photon_et  = -9.0;
   float photon_eta = -9.0;
   float photon_phi = -9.0;
   bool  photon_passPhotonID = false;
   bool  found_lead_pho = false;
   int   photon_count_bar = 0;
   int   photon_count_end = 0;
   // Assumption: reco photons are ordered by Et
   for (PhotonCollection::const_iterator recoPhoton = photonCollection->begin(); recoPhoton!=photonCollection->end(); recoPhoton++){
 
     // stop looping over photons once we get to too low Et
     if ( recoPhoton->et() < theMinPhotonEt ) break;
 
     //  Ignore ECAL Spikes
     const reco::CaloClusterPtr  seed = recoPhoton->superCluster()->seed();
     DetId id = lazyTool.getMaximum(*seed).first; // Cluster shape variables
     float time  = -999., outOfTimeChi2 = -999., chi2 = -999.;
     int   flags=-1, severity = -1; 
     const EcalRecHitCollection & rechits = ( recoPhoton->isEB() ? *EBReducedRecHits : *EEReducedRecHits); 
     EcalRecHitCollection::const_iterator it = rechits.find( id );
     if( it != rechits.end() ) {
       time = it->time(); 
       outOfTimeChi2 = it->outOfTimeChi2();
       chi2 = it->chi2();
       flags = it->recoFlag();
 
       edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
       iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
       severity = sevlv->severityLevel( id, rechits);
     }
     bool isNotSpike = ((recoPhoton->isEB() && (severity!=3 && severity!=4 ) && (flags != 2) ) || recoPhoton->isEE());
     if (!isNotSpike) continue;  // move on to next photon
     // END of determining ECAL Spikes
 
     bool pho_current_passPhotonID = false;
     bool pho_current_isEB = recoPhoton->isEB();
     bool pho_current_isEE = recoPhoton->isEE();
 
     if ( pho_current_isEB && (recoPhoton->sigmaIetaIeta() < 0.01 || recoPhoton->hadronicOverEm() < 0.05) ) {
       // Photon object in barrel passes photon ID
       pho_current_passPhotonID = true;
       photon_count_bar++;
     } else if ( pho_current_isEE && (recoPhoton->hadronicOverEm() < 0.05) ) {
       // Photon object in endcap passes photon ID
       pho_current_passPhotonID = true;
       photon_count_end++;
     }
 
     if (!found_lead_pho) {
       found_lead_pho = true;
       photon_passPhotonID = pho_current_passPhotonID;
       photon_et  = recoPhoton->et();
       photon_eta = recoPhoton->eta();
       photon_phi = recoPhoton->phi();
     }
   }
   
   // If user requires a photon to be found, but none is, return.
   //   theRequirePhotonFound should pretty much always be set to 'True'
   //    except when running on qcd monte carlo just to see the jets.
   if ( theRequirePhotonFound && (!photon_passPhotonID || photon_et<theMinPhotonEt) ) return;
 
 
   // Find the highest et jet
   Handle<CaloJetCollection> caloJetCollection;
   iEvent.getByLabel (theCaloJetCollectionLabel,caloJetCollection);
   if (!caloJetCollection.isValid()) return;
 
   float jet_pt    = -8.0;
   float jet_eta   = -8.0;
   float jet_phi   = -8.0;
   int   jet_count = 0;
   float jet2_pt   = -9.0;
   float jet2_eta  = -9.0;
   float jet2_phi  = -9.0;
   // Assumption: jets are ordered by Et
   for (CaloJetCollection::const_iterator i_calojet = caloJetCollection->begin(); i_calojet != caloJetCollection->end(); i_calojet++) {
 
     float jet_current_pt = i_calojet->pt();
 
     // don't care about jets that overlap with the lead photon
     if ( deltaR(i_calojet->eta(), i_calojet->phi(), photon_eta, photon_phi) < 0.5 ) continue;
     // stop looping over jets once we get to too low Et
     if (jet_current_pt < theMinCaloJetPt) break;
 
     jet_count++;
     if (jet_current_pt > jet_pt) {
       jet2_pt  = jet_pt;  // 2nd highest jet get's et from current highest
       jet2_eta = jet_eta;
       jet2_phi = jet_phi;
       jet_pt   = jet_current_pt; // current highest jet gets et from the new highest
       jet_eta  = i_calojet->eta();
       jet_phi  = i_calojet->phi();
     } else if (jet_current_pt > jet2_pt) {
       jet2_pt  = jet_current_pt;
       jet2_eta = i_calojet->eta();
       jet2_phi = i_calojet->phi();
     }
   }
 
 
   // Fill histograms if a jet found
   // NOTE: if a photon was required to be found, but wasn't
   //        we wouldn't have made it to this point in the code
   if ( jet_pt > 0.0 ) {
 
     // Photon Plots
     h_photon_et       ->Fill( photon_et  );
     h_photon_eta      ->Fill( photon_eta );
     h_photon_count_bar->Fill( photon_count_bar );
     h_photon_count_end->Fill( photon_count_end );
 
     // Photon Et hists for different orientations to the jet
     if ( fabs(photon_eta)<1.45 && photon_passPhotonID ) {  // Lead photon is in barrel
       if (fabs(jet_eta)<1.45){                          //   jet is in barrel
     if (photon_eta*jet_eta>0) {
       h_photon_et_jetcs->Fill(photon_et);
     } else {
       h_photon_et_jetco->Fill(photon_et);
     }
       } else if (jet_eta>1.55 && jet_eta<2.5) {         // jet is in endcap
     if (photon_eta*jet_eta>0) {
       h_photon_et_jetfs->Fill(photon_et);
     } else {
       h_photon_et_jetfo->Fill(photon_et);
     }
       }
     } // END of Lead Photon is in Barrel
 
     // Jet Plots
     h_jet_pt       ->Fill( jet_pt     );
     h_jet_eta      ->Fill( jet_eta    );
     h_jet_count    ->Fill( jet_count  );
     h_deltaPhi_photon_jet   ->Fill( abs(deltaPhi(photon_phi, jet_phi)) );
     if ( abs(deltaPhi(photon_phi,jet_phi))>2.8 ) h_deltaEt_photon_jet->Fill( (photon_et-jet_pt)/photon_et );
 
     // 2nd Highest Jet Plots
     if ( jet2_pt  > 0.0 ) {
       h_jet2_pt             ->Fill( jet2_pt  );
       h_jet2_eta            ->Fill( jet2_eta );
       h_jet2_ptOverPhotonEt ->Fill( jet2_pt/photon_et );
       h_deltaPhi_photon_jet2->Fill( abs(deltaPhi(photon_phi, jet2_phi)) );
       h_deltaPhi_jet_jet2   ->Fill( abs(deltaPhi(   jet_phi, jet2_phi)) );
       h_deltaR_jet_jet2     ->Fill( deltaR(   jet_eta,    jet_phi, jet2_eta, jet2_phi) );
       h_deltaR_photon_jet2  ->Fill( deltaR(photon_eta, photon_phi, jet2_eta, jet2_phi) );
     }
   } 
   // End of Filling histograms
 }

void QcdPhotonsDQM::beginJob ( void )

virtual

Inizialize parameters for histo binning.

Reimplemented from edm::EDAnalyzer.

Definition at line 89 of file QcdPhotonsDQM.cc.

References i, LogTrace, and cmsCodeRules.cppFunctionSkipper::operator.

                              {
  
   logTraceName = "QcdPhotonAnalyzer";
 
   LogTrace(logTraceName)<<"Parameters initialization";
   theDbe = Service<DQMStore>().operator->();
  
   theDbe->setCurrentFolder("Physics/QcdPhotons");  // Use folder with name of PAG
 
   std::stringstream aStringStream;
   std::string aString;
   aStringStream << theMinCaloJetPt;
   aString = aStringStream.str();
 
   // Monitor of triggers passed
   int numOfTriggersToMonitor = thePlotTheseTriggersToo.size();
   h_triggers_passed = theDbe->book1D("h_triggers_passed", "Events passing these trigger paths", numOfTriggersToMonitor, 0, numOfTriggersToMonitor);
   for (int i=0; i<numOfTriggersToMonitor; i++) {
     h_triggers_passed->setBinLabel(i+1,thePlotTheseTriggersToo[i]);
   }
 
   // Keep the number of plots and number of bins to a minimum!
   h_photon_et           = theDbe->book1D("h_photon_et",       "#gamma with highest E_{T};E_{T}(#gamma) (GeV)", 20, 0., thePlotPhotonMaxEt);
   h_photon_eta          = theDbe->book1D("h_photon_eta",      "#gamma with highest E_{T};#eta(#gamma)", 40, -thePlotPhotonMaxEta, thePlotPhotonMaxEta);
   h_photon_count_bar    = theDbe->book1D("h_photon_count_bar","Number of #gamma's passing selection (Barrel);Number of #gamma's", 8, -0.5, 7.5);
   h_photon_count_end    = theDbe->book1D("h_photon_count_end","Number of #gamma's passing selection (Endcap);Number of #gamma's", 8, -0.5, 7.5);
 
   h_jet_pt             = theDbe->book1D("h_jet_pt",   "Jet with highest p_{T} (from "+theCaloJetCollectionLabel.label()+");p_{T}(1^{st} jet) (GeV)",    20, 0., thePlotPhotonMaxEt);
   h_jet_eta             = theDbe->book1D("h_jet_eta", "Jet with highest p_{T} (from "+theCaloJetCollectionLabel.label()+");#eta(1^{st} jet)", 20, -thePlotJetMaxEta, thePlotJetMaxEta);
   h_deltaPhi_photon_jet = theDbe->book1D("h_deltaPhi_photon_jet", "#Delta#phi between Highest E_{T} #gamma and jet;#Delta#phi(#gamma,1^{st} jet)", 20, 0, 3.1415926);
   h_deltaPhi_jet_jet2   = theDbe->book1D("h_deltaPhi_jet_jet2", "#Delta#phi between Highest E_{T} jet and 2^{nd} jet;#Delta#phi(1^{st} jet,2^{nd} jet)", 20, 0, 3.1415926);
   h_deltaEt_photon_jet  = theDbe->book1D("h_deltaEt_photon_jet",  "(E_{T}(#gamma)-p_{T}(jet))/E_{T}(#gamma) when #Delta#phi(#gamma,1^{st} jet) > 2.8;#DeltaE_{T}(#gamma,1^{st} jet)/E_{T}(#gamma)", 20, -1.0, 1.0);
   h_jet_count           = theDbe->book1D("h_jet_count", "Number of "+theCaloJetCollectionLabel.label()+" (p_{T} > "+aString+" GeV);Number of Jets", 8, -0.5, 7.5);
   h_jet2_pt             = theDbe->book1D("h_jet2_pt",   "Jet with 2^{nd} highest p_{T} (from "+theCaloJetCollectionLabel.label()+");p_{T}(2^{nd} jet) (GeV)",    20, 0., thePlotPhotonMaxEt);
   h_jet2_eta            = theDbe->book1D("h_jet2_eta", "Jet with 2^{nd} highest p_{T} (from "+theCaloJetCollectionLabel.label()+");#eta(2^{nd} jet)", 20, -thePlotJetMaxEta, thePlotJetMaxEta);
   h_jet2_ptOverPhotonEt = theDbe->book1D("h_jet2_ptOverPhotonEt", "p_{T}(2^{nd} highest jet) / E_{T}(#gamma);p_{T}(2^{nd} Jet)/E_{T}(#gamma)", 20, 0.0, 4.0);
   h_deltaPhi_photon_jet2 = theDbe->book1D("h_deltaPhi_photon_jet2","#Delta#phi between Highest E_{T} #gamma and 2^{nd} highest jet;#Delta#phi(#gamma,2^{nd} jet)", 20, 0, 3.1415926);
   h_deltaR_jet_jet2      = theDbe->book1D("h_deltaR_jet_jet2", "#DeltaR between Highest Jet and 2^{nd} Highest;#DeltaR(1^{st} jet,2^{nd} jet)", 30, 0, 6.0);
   h_deltaR_photon_jet2   = theDbe->book1D("h_deltaR_photon_jet2", "#DeltaR between Highest E_{T} #gamma and 2^{nd} jet;#DeltaR(#gamma, 2^{nd} jet)", 30, 0, 6.0);
 
   // Photon Et for different jet configurations
   Float_t bins_et[] = {15,20,30,50,80};
   int num_bins_et = 4;
   h_photon_et_jetcs = theDbe->book1D("h_photon_et_jetcs", "#gamma with highest E_{T} (#eta(jet)<1.45, #eta(#gamma)#eta(jet)>0);E_{T}(#gamma) (GeV)", num_bins_et, bins_et);
   h_photon_et_jetco = theDbe->book1D("h_photon_et_jetco", "#gamma with highest E_{T} (#eta(jet)<1.45, #eta(#gamma)#eta(jet)<0);E_{T}(#gamma) (GeV)", num_bins_et, bins_et);
   h_photon_et_jetfs = theDbe->book1D("h_photon_et_jetfs", "#gamma with highest E_{T} (1.55<#eta(jet)<2.5, #eta(#gamma)#eta(jet)>0);E_{T}(#gamma) (GeV)", num_bins_et, bins_et);
   h_photon_et_jetfo = theDbe->book1D("h_photon_et_jetfo", "#gamma with highest E_{T} (1.55<#eta(jet)<2.5, #eta(#gamma)#eta(jet)<0);E_{T}(#gamma) (GeV)", num_bins_et, bins_et);
   h_photon_et_jetcs->getTH1F()->Sumw2();
   h_photon_et_jetco->getTH1F()->Sumw2();
   h_photon_et_jetfs->getTH1F()->Sumw2();
   h_photon_et_jetfo->getTH1F()->Sumw2();
   // Ratio of the above Photon Et distributions
   h_photon_et_ratio_co_cs = theDbe->book1D("h_photon_et_ratio_00_co_cs", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)<0) / D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)>0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
   h_photon_et_ratio_fo_fs = theDbe->book1D("h_photon_et_ratio_01_fo_fs", "D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)<0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)>0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
   h_photon_et_ratio_cs_fs = theDbe->book1D("h_photon_et_ratio_02_cs_fs", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)>0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)>0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
   h_photon_et_ratio_co_fs = theDbe->book1D("h_photon_et_ratio_03_co_fs", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)<0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)>0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
   h_photon_et_ratio_cs_fo = theDbe->book1D("h_photon_et_ratio_04_cs_fo", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)>0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)<0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
   h_photon_et_ratio_co_fo = theDbe->book1D("h_photon_et_ratio_05_co_fo", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)<0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)<0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
   h_photon_et_ratio_co_cs->getTH1F()->Sumw2();
   h_photon_et_ratio_fo_fs->getTH1F()->Sumw2();
   h_photon_et_ratio_cs_fs->getTH1F()->Sumw2();
   h_photon_et_ratio_co_fs->getTH1F()->Sumw2();
   h_photon_et_ratio_cs_fo->getTH1F()->Sumw2();
   h_photon_et_ratio_co_fo->getTH1F()->Sumw2();
 }

void QcdPhotonsDQM::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	iSetup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 159 of file QcdPhotonsDQM.cc.

                                                                            {
 
   // passed as parameter to HLTConfigProvider::init(), not yet used
   bool isConfigChanged = false;
 
   // isValidHltConfig_ used to short-circuit analyze() in case of problems
   isValidHltConfig_ = hltConfigProvider_.init( r, iSetup, theHltMenu, isConfigChanged );
 
   num_events_in_run = 0;
 }

void QcdPhotonsDQM::endJob ( void )

virtual

Save the histos.

Reimplemented from edm::EDAnalyzer.

Definition at line 400 of file QcdPhotonsDQM.cc.

400 {}

void QcdPhotonsDQM::endRun	(	const edm::Run &	run,
		const edm::EventSetup &	es
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 402 of file QcdPhotonsDQM.cc.

                                                                      {
   if (num_events_in_run>0) { 
     h_triggers_passed->getTH1F()->Scale(1.0/num_events_in_run);
   }
   h_photon_et_ratio_co_cs->getTH1F()->Divide( h_photon_et_jetco->getTH1F(), h_photon_et_jetcs->getTH1F() );
   h_photon_et_ratio_fo_fs->getTH1F()->Divide( h_photon_et_jetfo->getTH1F(), h_photon_et_jetfs->getTH1F() );
   h_photon_et_ratio_cs_fs->getTH1F()->Divide( h_photon_et_jetcs->getTH1F(), h_photon_et_jetfs->getTH1F() );
   h_photon_et_ratio_co_fs->getTH1F()->Divide( h_photon_et_jetco->getTH1F(), h_photon_et_jetfs->getTH1F() );
   h_photon_et_ratio_cs_fo->getTH1F()->Divide( h_photon_et_jetcs->getTH1F(), h_photon_et_jetfo->getTH1F() );
   h_photon_et_ratio_co_fo->getTH1F()->Divide( h_photon_et_jetco->getTH1F(), h_photon_et_jetfo->getTH1F() );
 }

Member Data Documentation

MonitorElement* QcdPhotonsDQM::h_deltaEt_photon_jet

private

Definition at line 91 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_deltaPhi_jet_jet2

private

Definition at line 90 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_deltaPhi_photon_jet

private

Definition at line 89 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_deltaPhi_photon_jet2

private

Definition at line 95 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_deltaR_jet_jet2

private

Definition at line 96 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_deltaR_photon_jet2

private

Definition at line 97 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_jet2_eta

private

Definition at line 94 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_jet2_pt

private

Definition at line 93 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_jet2_ptOverPhotonEt

private

Definition at line 92 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_jet_count

private

Definition at line 88 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_jet_eta

private

Definition at line 87 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_jet_pt

private

Definition at line 86 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_count_bar

private

Definition at line 84 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_count_end

private

Definition at line 85 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et

private

Definition at line 82 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_jetco

private

Definition at line 100 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_jetcs

private

Definition at line 99 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_jetfo

private

Definition at line 102 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_jetfs

private

Definition at line 101 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_ratio_co_cs

private

Definition at line 104 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_ratio_co_fo

private

Definition at line 109 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_ratio_co_fs

private

Definition at line 107 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_ratio_cs_fo

private

Definition at line 108 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_ratio_cs_fs

private

Definition at line 106 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_et_ratio_fo_fs

private

Definition at line 105 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_photon_eta

private

Definition at line 83 of file QcdPhotonsDQM.h.

MonitorElement* QcdPhotonsDQM::h_triggers_passed

private

Definition at line 81 of file QcdPhotonsDQM.h.

HLTConfigProvider QcdPhotonsDQM::hltConfigProvider_

private

Definition at line 57 of file QcdPhotonsDQM.h.

bool QcdPhotonsDQM::isValidHltConfig_

private

Definition at line 58 of file QcdPhotonsDQM.h.

std::string QcdPhotonsDQM::logTraceName

private

Definition at line 61 of file QcdPhotonsDQM.h.

int QcdPhotonsDQM::num_events_in_run

private

Definition at line 78 of file QcdPhotonsDQM.h.

edm::InputTag QcdPhotonsDQM::theCaloJetCollectionLabel

private

Definition at line 69 of file QcdPhotonsDQM.h.

DQMStore* QcdPhotonsDQM::theDbe

private

Definition at line 55 of file QcdPhotonsDQM.h.

std::string QcdPhotonsDQM::theHltMenu

private

Definition at line 66 of file QcdPhotonsDQM.h.

double QcdPhotonsDQM::theMinCaloJetPt

private

Definition at line 71 of file QcdPhotonsDQM.h.

double QcdPhotonsDQM::theMinPhotonEt

private

Definition at line 72 of file QcdPhotonsDQM.h.

edm::InputTag QcdPhotonsDQM::thePhotonCollectionLabel

private

Definition at line 68 of file QcdPhotonsDQM.h.

double QcdPhotonsDQM::thePlotJetMaxEta

private

Definition at line 76 of file QcdPhotonsDQM.h.

double QcdPhotonsDQM::thePlotPhotonMaxEt

private

Definition at line 74 of file QcdPhotonsDQM.h.

double QcdPhotonsDQM::thePlotPhotonMaxEta

private

Definition at line 75 of file QcdPhotonsDQM.h.

std::vector<std::string> QcdPhotonsDQM::thePlotTheseTriggersToo

private

Definition at line 65 of file QcdPhotonsDQM.h.

bool QcdPhotonsDQM::theRequirePhotonFound

private

Definition at line 73 of file QcdPhotonsDQM.h.

std::string QcdPhotonsDQM::theTriggerPathToPass

private

Definition at line 64 of file QcdPhotonsDQM.h.

std::string QcdPhotonsDQM::theTriggerResultsCollection

private

Definition at line 67 of file QcdPhotonsDQM.h.

edm::InputTag QcdPhotonsDQM::theVertexCollectionLabel

private

Definition at line 70 of file QcdPhotonsDQM.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation