#include <MuonAnalysis/MomentumScaleCalibration/plugins/TestResolution.cc>

Inheritance diagram for TestResolution:

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

	~TestResolution () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	endJob () override

template<typename T >
std::vector< reco::LeafCandidate >	fillMuonCollection (const std::vector< T > &tracks)

Private Attributes
int	eventCounter_

edm::EDGetTokenT< reco::MuonCollection >	glbMuonsToken_

TFile *	outputFile_

std::unique_ptr< ResolutionFunction >	resolutionFunction_

edm::EDGetTokenT< reco::TrackCollection >	saMuonsToken_

TProfile *	sigmaPt_

edm::InputTag	theMuonLabel_

int	theMuonType_

std::string	theRootFileName_

edm::EDGetTokenT< reco::TrackCollection >	tracksToken_

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

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 53 of file TestResolution.h.

Constructor & Destructor Documentation

◆ TestResolution()

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

explicit

Definition at line 20 of file TestResolution.cc.

     : theMuonLabel_(iConfig.getParameter<edm::InputTag>("MuonLabel")),
       glbMuonsToken_(mayConsume<reco::MuonCollection>(theMuonLabel_)),
       saMuonsToken_(mayConsume<reco::TrackCollection>(theMuonLabel_)),
       tracksToken_(mayConsume<reco::TrackCollection>(theMuonLabel_)),
       theMuonType_(iConfig.getParameter<int>("MuonType")),
       theRootFileName_(iConfig.getUntrackedParameter<std::string>("OutputFileName")) {
   //now do what ever initialization is needed
   outputFile_ = new TFile(theRootFileName_.c_str(), "RECREATE");
   outputFile_->cd();
   sigmaPt_ = new TProfile("sigmaPtOverPt", "sigmaPt/Pt vs muon Pt", 1000, 0, 100);
   eventCounter_ = 0;
   // Create the corrector and set the parameters
   resolutionFunction_.reset(
       new ResolutionFunction(iConfig.getUntrackedParameter<std::string>("ResolutionsIdentifier")));
   std::cout << "resolutionFunction_ = " << &*resolutionFunction_ << std::endl;
 }

References gather_cfg::cout, eventCounter_, edm::ParameterSet::getUntrackedParameter(), outputFile_, resolutionFunction_, sigmaPt_, AlCaHLTBitMon_QueryRunRegistry::string, and theRootFileName_.

◆ ~TestResolution()

TestResolution::~TestResolution ( )

override

Definition at line 38 of file TestResolution.cc.

                                 {
   outputFile_->cd();
   TCanvas canvas("sigmaPtOverPt", "sigmaPt/Pt vs muon Pt", 1000, 800);
   canvas.cd();
   sigmaPt_->GetXaxis()->SetTitle("Pt(GeV)");
   //   TLegend * legend = new TLegend(0.7,0.71,0.98,1.);
   //   legend->SetTextSize(0.02);
   //   legend->SetFillColor(0); // Have a white background
   //   legend->AddEntry(uncorrectedPt_, "original pt");
   //   legend->AddEntry(correctedPt_, "corrected pt");
   sigmaPt_->Draw();
   //   legend->Draw("same");
  
   canvas.Write();
   sigmaPt_->Write();
   outputFile_->Close();
  
   std::cout << "Total analyzed events = " << eventCounter_ << std::endl;
 }

References svgfig::canvas(), gather_cfg::cout, eventCounter_, outputFile_, and sigmaPt_.

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 63 of file TestResolution.cc.

                                                                                 {
   using namespace edm;
  
   ++eventCounter_;
   if (eventCounter_ % 100 == 0) {
     std::cout << "Event number " << eventCounter_ << std::endl;
   }
  
   // Take the reco-muons, depending on the type selected in the cfg
   // --------------------------------------------------------------
  
   std::vector<reco::LeafCandidate> muons;
  
   if (theMuonType_ == 1) {  // GlobalMuons
     Handle<reco::MuonCollection> glbMuons;
     iEvent.getByToken(glbMuonsToken_, glbMuons);
     muons = fillMuonCollection(*glbMuons);
   } else if (theMuonType_ == 2) {  // StandaloneMuons
     Handle<reco::TrackCollection> saMuons;
     iEvent.getByToken(saMuonsToken_, saMuons);
     muons = fillMuonCollection(*saMuons);
   } else if (theMuonType_ == 3) {  // Tracker tracks
     Handle<reco::TrackCollection> tracks;
     iEvent.getByToken(tracksToken_, tracks);
     muons = fillMuonCollection(*tracks);
   }
  
   // Loop on the recMuons
   std::vector<reco::LeafCandidate>::const_iterator recMuon = muons.begin();
   for (; recMuon != muons.end(); ++recMuon) {
     // Fill the histogram with uncorrected pt values
     sigmaPt_->Fill(resolutionFunction_->sigmaPt(*recMuon, 0), recMuon->pt());
   }
 }

References gather_cfg::cout, eventCounter_, fillMuonCollection(), glbMuonsToken_, iEvent, PDWG_BPHSkim_cff::muons, resolutionFunction_, saMuonsToken_, sigmaPt_, theMuonType_, PDWG_EXOHSCP_cff::tracks, and tracksToken_.

◆ endJob()

void TestResolution::endJob ( void )

inlineoverrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 60 of file TestResolution.h.

60 {};

◆ fillMuonCollection()

template<typename T >

std::vector<reco::LeafCandidate> TestResolution::fillMuonCollection ( const std::vector< T > & tracks )

inlineprivate

Definition at line 62 of file TestResolution.h.

                                                                               {
     std::vector<reco::LeafCandidate> muons;
     typename std::vector<T>::const_iterator track;
     for (track = tracks.begin(); track != tracks.end(); ++track) {
       // Where 0.011163612 is the squared muon mass.
       reco::Particle::LorentzVector mu(
           track->px(), track->py(), track->pz(), sqrt(track->p() * track->p() + 0.011163612));
       reco::LeafCandidate muon(track->charge(), mu);
       // Store muon
       // ----------
       muons.push_back(muon);
     }
     return muons;
   }

References amptDefaultParameters_cff::mu, HLT_2018_cff::muon, PDWG_BPHSkim_cff::muons, mathSSE::sqrt(), HLT_2018_cff::track, and PDWG_EXOHSCP_cff::tracks.

Referenced by analyze().

Member Data Documentation

◆ eventCounter_

int TestResolution::eventCounter_

private

Definition at line 92 of file TestResolution.h.

Referenced by analyze(), TestResolution(), and ~TestResolution().

◆ glbMuonsToken_

edm::EDGetTokenT<reco::MuonCollection> TestResolution::glbMuonsToken_

private

Definition at line 82 of file TestResolution.h.

Referenced by analyze().

◆ outputFile_

TFile* TestResolution::outputFile_

private

Definition at line 88 of file TestResolution.h.

Referenced by TestResolution(), and ~TestResolution().

◆ resolutionFunction_

std::unique_ptr<ResolutionFunction> TestResolution::resolutionFunction_

private

Definition at line 94 of file TestResolution.h.

Referenced by analyze(), and TestResolution().

◆ saMuonsToken_

edm::EDGetTokenT<reco::TrackCollection> TestResolution::saMuonsToken_

private

Definition at line 83 of file TestResolution.h.

Referenced by analyze().

◆ sigmaPt_

TProfile* TestResolution::sigmaPt_

private

Definition at line 90 of file TestResolution.h.

Referenced by analyze(), TestResolution(), and ~TestResolution().

◆ theMuonLabel_

edm::InputTag TestResolution::theMuonLabel_

private

Definition at line 81 of file TestResolution.h.

◆ theMuonType_

int TestResolution::theMuonType_

private

Definition at line 86 of file TestResolution.h.

Referenced by analyze().

◆ theRootFileName_

std::string TestResolution::theRootFileName_

private

Definition at line 87 of file TestResolution.h.

Referenced by TestResolution().

◆ tracksToken_

edm::EDGetTokenT<reco::TrackCollection> TestResolution::tracksToken_

private

Definition at line 84 of file TestResolution.h.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TestResolution()

◆ ~TestResolution()

Member Function Documentation

◆ analyze()

◆ endJob()

◆ fillMuonCollection()

Member Data Documentation

◆ eventCounter_

◆ glbMuonsToken_

◆ outputFile_

◆ resolutionFunction_

◆ saMuonsToken_

◆ sigmaPt_

◆ theMuonLabel_

◆ theMuonType_

◆ theRootFileName_

◆ tracksToken_