HLTTopPlotter Class Reference

#include <HLTTopPlotter.h>

Inheritance diagram for HLTTopPlotter:

Public Member Functions
void	analyze (const edm::Event &iEvent)
void	begin ()
void	finish ()
	HLTTopPlotter (const edm::ParameterSet &pset, std::string triggerName, std::vector< std::string > moduleNames, MuonSelectionStruct inputSelection, std::string customName, std::vector< std::string > validTriggers, const edm::Run &currentRun, const edm::EventSetup &currentEventSetup)
	Constructor. More...
Public Member Functions inherited from HLTMuonMatchAndPlot
virtual MonitorElement *	bookIt (TString name, TString title, std::vector< double >)
virtual MonitorElement *	bookIt (TString name, TString title, int nbins, float *xBinLowEdges)
TString	calcHistoSuffix (std::string moduleName)
virtual void	endRun (const edm::Run &r, const edm::EventSetup &c)
virtual void	fillPlots (std::vector< MatchStruct > &myRecMatches, std::vector< std::vector< HltFakeStruct > > &myHltFakeCands)
	HLTMuonMatchAndPlot (const edm::ParameterSet &pset, std::string triggerName, std::vector< std::string > moduleNames, MuonSelectionStruct inputSelection, std::string customName, std::vector< std::string > validTriggers, const edm::Run &currentRun, const edm::EventSetup &currentEventSetup)
	Constructor. More...
virtual bool	selectAndMatchMuons (const edm::Event &iEvent, std::vector< MatchStruct > &myRecMatches, std::vector< std::vector< HltFakeStruct > > &myHltFakeCands)
bool	selectAndMatchMuons (const edm::Event &iEvent, std::vector< MatchStruct > &myRecMatches, std::vector< std::vector< HltFakeStruct > > &myHltFakeCands, MuonSelectionStruct muonSelection)
virtual	~HLTMuonMatchAndPlot ()

Private Member Functions
void	sortJets (reco::CaloJetCollection &theJets)

Private Attributes
edm::InputTag	CaloJetInputTag
double	EtaCut_
std::vector< MonitorElement * >	hDeltaRMaxJetLep
std::vector< MonitorElement * >	hJetMultip
int	NJets_
double	PtCut_
std::vector< double >	theJetMParameters

Additional Inherited Members
Public Attributes inherited from HLTMuonMatchAndPlot
std::vector< std::vector < HltFakeStruct > >	hltFakeCands
std::vector< MatchStruct >	recMatches
Protected Member Functions inherited from HLTMuonMatchAndPlot
virtual bool	applyTrackSelection (MuonSelectionStruct mySelection, reco::Muon candMuon)
virtual bool	applyTriggerSelection (MuonSelectionStruct mySelection, const edm::Event &event)
int	findGenMatch (double eta, double phi, double maxDeltaR, std::vector< MatchStruct > matches)
const reco::Candidate *	findMother (const reco::Candidate *)
int	findRecMatch (double eta, double phi, double maxdeltaR, std::vector< MatchStruct > matches)
void	getAodTriggerObjectsForModule (edm::InputTag collectionTag, edm::Handle< trigger::TriggerEvent > aodTriggerEvent, trigger::TriggerObjectCollection trigObjs, std::vector< trigger::TriggerObject > &foundObjects, MuonSelectionStruct muonSelection)
reco::TrackRef	getCandTrackRef (MuonSelectionStruct mySelection, reco::Muon candMuon)
int	getCharge (int pdgId)
void	moveOverflow (MonitorElement *myElement)
Protected Attributes inherited from HLTMuonMatchAndPlot
std::vector< MonitorElement * >	allHltCandEta
std::vector< MonitorElement * >	allHltCandPhi
std::vector< MonitorElement * >	allHltCandPt
reco::BeamSpot	beamSpot
edm::InputTag	BeamSpotInputTag
std::vector< MonitorElement * >	booked1DMonitorElements
bool	createStandAloneHistos
DQMStore *	dbe_
int	eventNumber
std::vector< MonitorElement * >	fakeHltCandEta
std::vector< MonitorElement * >	fakeHltCandPhi
std::vector< MonitorElement * >	fakeHltCandPt
bool	foundBeamSpot
std::vector< MonitorElement * >	hBeamSpotZ0Rec
std::vector< MonitorElement * >	hChargeFlipMatched
std::vector< MonitorElement * >	hDeltaRMatched
std::vector< MonitorElement * >	hIsolationRec
std::string	histoFileName
int	HLT_PLOT_OFFSET
edm::InputTag	HltAodInputTag
edm::InputTag	HltRawInputTag
std::vector< MonitorElement * >	hMatchedDeltaPhi
MonitorElement *	hNumObjects
MonitorElement *	hNumOrphansRec
std::vector< MonitorElement * >	hPassCharge
std::vector< MonitorElement * >	hPassD0BeamRec
std::vector< MonitorElement * >	hPassD0Rec
std::vector< MonitorElement * >	hPassEtaRec
std::vector< MonitorElement * >	hPassExaclyOneMuonMaxPtRec
std::vector< MonitorElement * >	hPassMatchPtRec
std::vector< MonitorElement * >	hPassMaxPtRec
std::vector< MonitorElement * >	hPassPhiRec
std::vector< MonitorElement * >	hPassPtRec
std::vector< MonitorElement * >	hPassPtRecExactlyOne
std::vector< MonitorElement * >	hPassZ0BeamRec
std::vector< MonitorElement * >	hPassZ0Rec
std::vector< MonitorElement * >	hPhiVsEtaRec
std::vector< MonitorElement * >	hResoEtaAodRec
std::vector< MonitorElement * >	hResoPhiAodRec
std::vector< MonitorElement * >	hResoPtAodRec
bool	includeOverflow
bool	isIsolatedPath
bool	isL1Path
bool	isL2Path
bool	isL3Path
bool	makeNtuple
std::string	matchType
MonitorElement *	meNumberOfEvents
MuonSelectionStruct	mySelection
int	numBinsInPtHisto
unsigned int	numHltLabels
unsigned int	numL1Cands
float	ptBins [100]
std::vector< MonitorElement * >	rawMatchHltCandEta
std::vector< MonitorElement * >	rawMatchHltCandPhi
std::vector< MonitorElement * >	rawMatchHltCandPt
edm::InputTag	RecoMuonInputTag
bool	requireL1SeedForHLTPaths
std::vector< std::string >	selectedValidTriggers
std::string	theAodL1Label
std::string	theAodL2Label
std::vector< double >	theChargeFlipParameters
std::vector< double >	theChargeParameters
std::vector< double >	theD0Parameters
std::vector< double >	theDRParameters
std::vector< double >	theEtaParameters
TFile *	theFile
std::vector< std::string >	theHltCollectionLabels
std::string	theHltProcessName
std::vector< double >	theIsolationParameters
std::string	theL1CollectionLabel
double	theL1DrCut
std::string	theL1SeedModuleForHLTPath
double	theL2DrCut
double	theL3DrCut
double	theMaxEtaCut
std::vector< double >	theMaxPtParameters
double	theMinPtCut
int	theMotherParticleId
std::vector< double >	theNSigmas
TNtuple *	theNtuple
std::string	theNtupleFileName
float	theNtuplePars [100]
std::string	theNtuplePath
unsigned int	theNumberOfObjects
std::vector< double >	thePhiEtaParameters2d
std::vector< double >	thePhiParameters
std::vector< double >	thePhiParameters0Pi
std::vector< double >	thePtParameters
std::vector< double >	theResParameters
std::string	theTriggerName
std::vector< double >	theZ0Parameters
edm::InputTag	TriggerResultLabel
bool	useAod
bool	useFullDebugInformation
bool	useMuonFromReco
bool	useOldLabels
Static Protected Attributes inherited from HLTMuonMatchAndPlot
static const int	NEG_CHARGE = -1
static const int	POS_CHARGE = 1

Detailed Description

Get L1/HLT efficiency/rate plots Documentation available on the CMS TWiki: https://twiki.cern.ch/twiki/bin/view/CMS/MuonHLTOfflinePerformance

Date:

2010/03/16 14:36:39

Revision:

1.5

Author

J. Slaunwhite (modified from above

Get L1/HLT efficiency/rate plots Documentation available on the CMS TWiki: https://twiki.cern.ch/twiki/bin/view/CMS/MuonHLTOfflinePerformance

Date:

2010/03/16 14:36:40

Revision:

1.5

Author

J. Slaunwhite (modified from above

Definition at line 57 of file HLTTopPlotter.h.

Constructor & Destructor Documentation

HLTTopPlotter::HLTTopPlotter	(	const edm::ParameterSet &	pset,
		std::string	triggerName,
		std::vector< std::string >	moduleNames,
		MuonSelectionStruct	inputSelection,
		std::string	customName,
		std::vector< std::string >	validTriggers,
		const edm::Run &	currentRun,
		const edm::EventSetup &	currentEventSetup
	)

Constructor.

Definition at line 62 of file HLTTopPlotter.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and LogTrace.

  : HLTMuonMatchAndPlot(pset, triggerName, moduleNames, inputSelection, customName, validTriggers, currentRun, currentEventSetup)
    
{

  CaloJetInputTag       = pset.getParameter<edm::InputTag>("CaloJetInputTag");
  EtaCut_               = pset.getUntrackedParameter<double>("EtaCut");
  PtCut_                = pset.getUntrackedParameter<double>("PtCut");   
  NJets_                = pset.getUntrackedParameter<int>("NJets");   
  theJetMParameters    = pset.getUntrackedParameter< vector<double> >("JetMParameters");

  
  LogTrace ("HLTMuonVal") << "\n\n Inside HLTTopPlotter Constructor";
  LogTrace ("HLTMuonVal") << "The trigger name is " << triggerName
                          << "and we've done all the other intitializations";

  LogTrace ("HLTMuonVal") << "exiting constructor\n\n";

}

Member Function Documentation

void HLTTopPlotter::analyze ( const edm::Event &  iEvent )

virtual

Reimplemented from HLTMuonMatchAndPlot.

Definition at line 100 of file HLTTopPlotter.cc.

References HLTMuonMatchAndPlot::analyze(), reco::deltaR(), edm::Event::getByLabel(), i, edm::HandleBase::isValid(), j, metsig::jet, analyzePatCleaning_cfg::jets, LogTrace, edm::Handle< T >::product(), and ExpressReco_HICollisions_FallBack::pt.

{

  LogTrace ("HLTMuonVal") << "Inside of TopPlotter analyze method!"
                          << "calling my match and plot module's analyze..."
                          << endl;

  // Make sure you are valid before proceeding


  // Do some top specific selection, then call the muon matching
  // if the event looks top-like

  LogTrace ("HLTMuonVal") << "Do top-specific selection" << endl;
  

  // get calo jet collection
  Handle<CaloJetCollection> jetsHandle;
  iEvent.getByLabel(CaloJetInputTag, jetsHandle);

  int n_jets_20 = 0;
  CaloJetCollection selectedJets;
  
  if (jetsHandle.isValid()) {
    LogTrace ("HLTMuonVal") << "Looking in jet collection" << endl;
    //Jet Collection to use
    
    // Raw jets
    const CaloJetCollection *jets = jetsHandle.product();
    CaloJetCollection::const_iterator jet;
    
    //int n_jets_20=0;

    // The parameters for the n jets should be 

    
    for (jet = jets->begin(); jet != jets->end(); jet++){        
      // if (fabs(jet->eta()) <2.4 && jet->et() > 20) n_jets_20++; 
      if (fabs(jet->eta()) <EtaCut_ && jet->et() > PtCut_) {
        n_jets_20++;
        selectedJets.push_back((*jet));
      }
      
    } 
    
  }
  
 

  // sort the result
  sortJets(selectedJets);

  
  
  LogTrace ("HLTMuonVal") << "Number of jets in this event = "
                          << n_jets_20
                          << endl;

 // if (n_jets_20 <= 1 ) {
  if (n_jets_20 < NJets_ ) {
    LogTrace ("HLTMuonVal") << "Not enought jets in this event, skipping it"
                            << endl;

    return;
  }
  
 

  //
  //     Call the other analyze method
  //
  


  LogTrace("HLTMuonVal") << "Calling analyze for muon ana" << endl;
  HLTMuonMatchAndPlot::analyze(iEvent);


  LogTrace ("HLTMuonVal") << "TOPPLOT: returned from muon ana, now in top module"
                          << endl
                          << "muon ana stored size rec muons = "
    //<< myMuonAna->recMatches.size()
                          << endl;
  
  vector<HLTMuonMatchAndPlot::MatchStruct>::const_iterator iRecMuon;

  int numCands = 0;
  for ( unsigned int i  = 0;
        i < recMatches.size();
        i++ ) {


    LogTrace ("HLTMuonVal") << "Cand " << numCands
                            << ", Pt = "
                            << recMatches[i].recCand->pt()
                            << ", eta = "
                            << recMatches[i].recCand->eta()
                            << ", phi = "
                            << recMatches[i].recCand->phi()
                            << endl;

    
    double deltaRLeadJetLep = reco::deltaR (recMatches[i].recCand->eta(), recMatches[i].recCand->phi(),
                                       selectedJets[0].eta(), selectedJets[0].phi());

    //
    //   Fill Plots for All
    //
    
    hDeltaRMaxJetLep[0]->Fill(deltaRLeadJetLep);
    hJetMultip[0]->Fill(n_jets_20);


    //
    //   Fill Plots for L1
    //

    
    if ( (recMatches[i].l1Cand.pt() > 0) && ((useFullDebugInformation) || (isL1Path)) ) {
      hDeltaRMaxJetLep[1]->Fill(deltaRLeadJetLep);
      hJetMultip[1]->Fill(n_jets_20);
    }

    //
    //   Fill Plots for HLT
    //
    
    for ( size_t j = 0; j < recMatches[i].hltCands.size(); j++ ) {
      if ( recMatches[i].hltCands[j].pt() > 0 ) {
        // you've found it!
        hDeltaRMaxJetLep[j+HLT_PLOT_OFFSET]->Fill(deltaRLeadJetLep);
    hJetMultip[j+HLT_PLOT_OFFSET]->Fill(n_jets_20);
        
      }
    }
    

    
    numCands++;
  }
  

  LogTrace ("HLTMuonVal") << "-----End of top plotter analyze method-----" << endl;
} // Done filling histograms

void HLTTopPlotter::begin ( void  )

virtual

Reimplemented from HLTMuonMatchAndPlot.

Definition at line 256 of file HLTTopPlotter.cc.

References HLTMuonMatchAndPlot::begin(), DQMStore::cd(), dbe_, h, i, testEve_cfg::level, LogTrace, and DQMStore::setCurrentFolder().

{

  TString myLabel, newFolder;
  vector<TH1F*> h;

  LogTrace ("HLTMuonVal") << "Inside begin for top analyzer" << endl;

  
  LogTrace ("HLTMuonVal") << "Calling begin for muon analyzer" << endl;
  HLTMuonMatchAndPlot::begin();

  LogTrace ("HLTMuonVal") << "Continuing with top analyzer begin" << endl;

  if ( dbe_ ) {
    dbe_->cd();
    dbe_->setCurrentFolder("HLT/Muon");

    
    // JMS I think this is trimming all L1 names to
    // to be L1Filtered
    myLabel = theL1CollectionLabel;
    myLabel = myLabel(myLabel.Index("L1"),myLabel.Length());
    myLabel = myLabel(0,myLabel.Index("Filtered")+8);


    // JMS Old way of doing things
    //newFolder = "HLT/Muon/Distributions/" + theTriggerName;
    newFolder = "HLT/Muon/Distributions/" + theTriggerName + "/" + mySelection.customLabel;

    
    
    dbe_->setCurrentFolder( newFolder.Data() );
    
    vector<string> binLabels;
    binLabels.push_back( theL1CollectionLabel.c_str() );
    for ( size_t i = 0; i < theHltCollectionLabels.size(); i++ )
      binLabels.push_back( theHltCollectionLabels[i].c_str() );


    //------- Define labels for plots -------
    
    if (useOldLabels) { 
      myLabel = theL1CollectionLabel;
      myLabel = myLabel(myLabel.Index("L1"),myLabel.Length());
      myLabel = myLabel(0,myLabel.Index("Filtered")+8);
    } else {
      myLabel = "L1Filtered";
    }

    //------ Definte the plots themselves------------------------


    //
    //         ALL + L1 plots 
    //
    
    hDeltaRMaxJetLep.push_back (bookIt("topDeltaRMaxJetLep_All", "delta R between muon and highest pt jet", theDRParameters));
    if (useFullDebugInformation || isL1Path) hDeltaRMaxJetLep.push_back (bookIt("topDeltaRMaxJetLep_" + myLabel, "delta R between muon and highest pt jet", theDRParameters));

   hJetMultip.push_back (bookIt("topJetMultip_All", "Jet multiplicity", theJetMParameters));
    if (useFullDebugInformation || isL1Path) hJetMultip.push_back(bookIt("topJetMultip_" + myLabel, "Jet multiplicity", theJetMParameters));


    //
    //         ALL + L1 plots 
    //
    
    
    // we won't enter this loop if we don't have an hlt label
    // we won't have an hlt label is this is a l1 path
    for ( unsigned int i = 0; i < theHltCollectionLabels.size(); i++ ) {

      if (useOldLabels) {
        myLabel = theHltCollectionLabels[i];
        TString level = ( myLabel.Contains("L2") ) ? "L2" : "L3";
        myLabel = myLabel(myLabel.Index(level),myLabel.Length());
        myLabel = myLabel(0,myLabel.Index("Filtered")+8);
      } else {
        TString tempString = theHltCollectionLabels[i];
        TString level = ( tempString.Contains("L2") ) ? "L2" : "L3";
        myLabel = level + "Filtered";
      } // end if useOldLabels
    

      // Book for L2, L3
      hDeltaRMaxJetLep.push_back (bookIt("topDeltaRMaxJetLep_" + myLabel, "delta R between muon and highest pt jet", theDRParameters)) ;
      hJetMultip.push_back (bookIt("topJetMultip_" + myLabel, "Jet Multiplicity",   theJetMParameters)) ;
      
      
      
    }// end for each collection label

  }// end if dbe_ exists

}// end begin method

void HLTTopPlotter::finish ( )

virtual

Reimplemented from HLTMuonMatchAndPlot.

Definition at line 88 of file HLTTopPlotter.cc.

References HLTMuonMatchAndPlot::finish().

{

  // you could do something else in here
  // but for now, just do what the base class
  // would have done
  
  HLTMuonMatchAndPlot::finish();
}

void HLTTopPlotter::sortJets ( reco::CaloJetCollection &  theJets )

private

Definition at line 370 of file HLTTopPlotter.cc.

References ExpressReco_HICollisions_FallBack::et, and LogTrace.

                                                         {


  LogTrace ("HLTMuonVal") << "Sorting Jets" << endl;
  
  // bubble sort jets
  
  for ( unsigned int iJet = 0;
        iJet < theJets.size();          
        iJet ++) {

    for ( unsigned int jJet = iJet;
          jJet < theJets.size();
          jJet ++ ){

      if ( theJets[jJet].et() > theJets[iJet].et() ) {
        reco::CaloJet tmpJet = theJets[iJet];
        theJets[iJet] = theJets[jJet];
        theJets[jJet] = tmpJet;        
      }
      
    }// end for each jJet

  }// end for each iJet

  for ( unsigned int iJet = 0;
        iJet != theJets.size();
        iJet ++ ) {

    LogTrace ("HLTMuonVal") << "Jet # " << iJet
                            << "  Et =  " << theJets[iJet].et()
                            << endl;

  }
  
}

Member Data Documentation

edm::InputTag HLTTopPlotter::CaloJetInputTag

private

Definition at line 91 of file HLTTopPlotter.h.

double HLTTopPlotter::EtaCut_

private

Definition at line 89 of file HLTTopPlotter.h.

std::vector<MonitorElement*> HLTTopPlotter::hDeltaRMaxJetLep

private

Definition at line 84 of file HLTTopPlotter.h.

std::vector<MonitorElement*> HLTTopPlotter::hJetMultip

private

Definition at line 85 of file HLTTopPlotter.h.

int HLTTopPlotter::NJets_

private

Definition at line 90 of file HLTTopPlotter.h.

double HLTTopPlotter::PtCut_

private

Definition at line 89 of file HLTTopPlotter.h.

std::vector<double> HLTTopPlotter::theJetMParameters

private

Definition at line 93 of file HLTTopPlotter.h.