PatTriggerTagAndProbe Class Reference

Inheritance diagram for PatTriggerTagAndProbe:

Public Member Functions
	PatTriggerTagAndProbe (const edm::ParameterSet &iConfig)
	default constructor More...
	~PatTriggerTagAndProbe ()
	default destructor More...
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup)
	everythin that needs to be done during the event loop More...
virtual void	beginJob ()
	everythin that needs to be done before the event loop More...
virtual void	endJob ()
	everythin that needs to be done after the event loop More...
void	setErrors (TH1D &h, const TH1D &ref)
	helper function to set proper bin errors More...

Private Attributes
std::map< std::string, TH1D * >	histos1D_
	management of 1d histograms More...
std::string	muonMatch_
	input for trigger match objects More...
edm::EDGetTokenT < pat::MuonCollection >	muonsToken_
	input for muons More...
edm::EDGetTokenT < pat::TriggerEvent >	triggerEventToken_
	input for patTriggerEvent More...

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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::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

Definition at line 16 of file PatTriggerTagAndProbe.cc.

Constructor & Destructor Documentation

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

explicit

default constructor

Definition at line 53 of file PatTriggerTagAndProbe.cc.

                                                                              :
  // pat::TriggerEvent
  triggerEventToken_( consumes< pat::TriggerEvent >( iConfig.getParameter< edm::InputTag >( "triggerEvent" ) ) ),
  // muon input collection
  muonsToken_( consumes< pat::MuonCollection >( iConfig.getParameter< edm::InputTag >( "muons" ) ) ),
  // muon match objects
  muonMatch_( iConfig.getParameter< std::string >( "muonMatch" ) ),
  // histogram management
  histos1D_()
{
}

PatTriggerTagAndProbe::~PatTriggerTagAndProbe

(

)

default destructor

Definition at line 65 of file PatTriggerTagAndProbe.cc.

{
}

Member Function Documentation

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

privatevirtual

everythin that needs to be done during the event loop

Implements edm::EDAnalyzer.

Definition at line 85 of file PatTriggerTagAndProbe.cc.

References edm::Event::getByToken(), histos1D_, muonMatch_, patZpeak::muons, muonsToken_, triggerEventToken_, and pat::helper::TriggerMatchHelper::triggerMatchObject().

{
  // trigger event
  edm::Handle< pat::TriggerEvent > triggerEvent;
  iEvent.getByToken( triggerEventToken_, triggerEvent );
  // pat candidate collection
  edm::Handle< pat::MuonCollection > muons;
  iEvent.getByToken( muonsToken_, muons );

  // pat trigger helper to recieve for trigger
  // matching information
  const pat::helper::TriggerMatchHelper matchHelper;

  // ask for trigger accept of HLT_Mu9; otherwise we don't even start
  if(!(triggerEvent->path("HLT_IsoMu17_v5")->wasRun() && triggerEvent->path("HLT_IsoMu17_v5")->wasAccept())){
    return;
  }

  // loop over muon references for the tag muon
  for( size_t idxTag=0; idxTag<muons->size(); ++idxTag){
    const pat::TriggerObjectRef trigRefTag( matchHelper.triggerMatchObject( muons, idxTag, muonMatch_, iEvent, *triggerEvent ) );
    if( trigRefTag.isAvailable() ){
      // loop over muon references for the probe/test muon
      for( size_t idxProbe=0; idxProbe<muons->size() && idxProbe!=idxTag; ++idxProbe){
    histos1D_[ "mass" ]->Fill( (muons->at(idxTag).p4()+muons->at(idxProbe).p4()).mass() );
    if(fabs((muons->at(idxTag).p4()+muons->at(idxProbe).p4()).mass()-90)<5){
      const pat::TriggerObjectRef trigRefProbe( matchHelper.triggerMatchObject( muons, idxProbe, muonMatch_, iEvent, *triggerEvent ) );
      histos1D_[ "probePt"  ]->Fill( muons->at(idxProbe).pt () );
      histos1D_[ "probeEta" ]->Fill( muons->at(idxProbe).eta() );
      if( trigRefProbe.isAvailable() ){
        histos1D_[ "testPt" ]->Fill( muons->at(idxProbe).pt () );
        histos1D_[ "testEta"]->Fill( muons->at(idxProbe).eta() );
      }
    }
      }
    }
  }
}

void PatTriggerTagAndProbe::beginJob ( void  )

privatevirtual

everythin that needs to be done before the event loop

Reimplemented from edm::EDAnalyzer.

Definition at line 69 of file PatTriggerTagAndProbe.cc.

References fileService, histos1D_, and TFileService::make().

{
  edm::Service< TFileService > fileService;

  // mass plot around Z peak
  histos1D_[ "mass"    ] = fileService->make< TH1D >( "mass"  , "Mass_{Z} (GeV)",  90,   30., 120.);
  // pt for test candidate
  histos1D_[ "testPt"  ] = fileService->make< TH1D >( "testPt"  , "p_{T} (GeV)" , 100,    0., 100.);
  // pt for probe candidate
  histos1D_[ "probePt" ] = fileService->make< TH1D >( "probePt" , "p_{T} (GeV)" , 100,    0., 100.);
  // eta for test candidate
  histos1D_[ "testEta" ] = fileService->make< TH1D >( "testEta" , "#eta"        ,  48,  -2.4,  2.4);
  // eta for probe candidate
  histos1D_[ "probeEta"] = fileService->make< TH1D >( "probeEta", "#eta"        ,  48,  -2.4,  2.4);
}

void PatTriggerTagAndProbe::endJob ( void  )

privatevirtual

everythin that needs to be done after the event loop

Reimplemented from edm::EDAnalyzer.

Definition at line 124 of file PatTriggerTagAndProbe.cc.

References histos1D_, and setErrors().

{
  // normalize the entries of the histograms
  histos1D_[ "testPt"  ]->Divide(histos1D_  [ "probePt"  ]);
  setErrors(*histos1D_["testPt" ],*histos1D_[ "probePt"  ]);
  histos1D_[ "testEta" ]->Divide(histos1D_  [ "probeEta" ]);
  setErrors(*histos1D_["testEta"],*histos1D_[ "probeEta" ]);
}

void PatTriggerTagAndProbe::setErrors ( TH1D &  h, const TH1D &  ref  )

private

helper function to set proper bin errors

Definition at line 133 of file PatTriggerTagAndProbe.cc.

References newFWLiteAna::bin, and mathSSE::sqrt().

Referenced by endJob().

{
  for(int bin=0; bin<h.GetNbinsX(); ++bin){
    if(ref.GetBinContent(bin+1)>0){
      h.SetBinError(bin+1, sqrt((h.GetBinContent(bin+1)*(1.-h.GetBinContent(bin+1)))/ref.GetBinContent(bin+1)));
    } else{ h.SetBinError(bin+1, 0.); }
  }
}

Member Data Documentation

std::map< std::string, TH1D* > PatTriggerTagAndProbe::histos1D_

private

management of 1d histograms

Definition at line 42 of file PatTriggerTagAndProbe.cc.

Referenced by analyze(), beginJob(), and endJob().

std::string PatTriggerTagAndProbe::muonMatch_

private

input for trigger match objects

Definition at line 40 of file PatTriggerTagAndProbe.cc.

Referenced by analyze().

edm::EDGetTokenT< pat::MuonCollection > PatTriggerTagAndProbe::muonsToken_

private

input for muons

Definition at line 38 of file PatTriggerTagAndProbe.cc.

Referenced by analyze().

edm::EDGetTokenT< pat::TriggerEvent > PatTriggerTagAndProbe::triggerEventToken_

private

input for patTriggerEvent

Definition at line 36 of file PatTriggerTagAndProbe.cc.

Referenced by analyze().