#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::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) 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::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

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

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

Definition at line 51 of file TestResolution.cc.

Constructor & Destructor Documentation

◆ TestResolution()

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

explicit

Definition at line 97 of file TestResolution.cc.

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

     : 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_ =
       std::make_unique<ResolutionFunction>(iConfig.getUntrackedParameter<std::string>("ResolutionsIdentifier"));
   std::cout << "resolutionFunction_ = " << &*resolutionFunction_ << std::endl;
 }

◆ ~TestResolution()

TestResolution::~TestResolution ( )

override

Definition at line 115 of file TestResolution.cc.

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

                                 {
   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;
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 140 of file TestResolution.cc.

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

                                                                                 {
   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());
   }
 }

◆ fillMuonCollection()

template<typename T >

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

inlineprivate

Definition at line 59 of file TestResolution.cc.

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

Referenced by analyze().

                                                                               {
     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;
   }

Member Data Documentation

◆ eventCounter_

int TestResolution::eventCounter_

private

Definition at line 89 of file TestResolution.cc.

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

◆ glbMuonsToken_

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

private

Definition at line 79 of file TestResolution.cc.

Referenced by analyze().

◆ outputFile_

TFile* TestResolution::outputFile_

private

Definition at line 85 of file TestResolution.cc.

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

◆ resolutionFunction_

std::unique_ptr<ResolutionFunction> TestResolution::resolutionFunction_

private

Definition at line 91 of file TestResolution.cc.

Referenced by analyze(), and TestResolution().

◆ saMuonsToken_

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

private

Definition at line 80 of file TestResolution.cc.

Referenced by analyze().

◆ sigmaPt_

TProfile* TestResolution::sigmaPt_

private

Definition at line 87 of file TestResolution.cc.

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

◆ theMuonLabel_

edm::InputTag TestResolution::theMuonLabel_

private

Definition at line 78 of file TestResolution.cc.

◆ theMuonType_

int TestResolution::theMuonType_

private

Definition at line 83 of file TestResolution.cc.

Referenced by analyze().

◆ theRootFileName_

std::string TestResolution::theRootFileName_

private

Definition at line 84 of file TestResolution.cc.

Referenced by TestResolution().

◆ tracksToken_

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

private

Definition at line 81 of file TestResolution.cc.

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

◆ fillMuonCollection()

Member Data Documentation

◆ eventCounter_

◆ glbMuonsToken_

◆ outputFile_

◆ resolutionFunction_

◆ saMuonsToken_

◆ sigmaPt_

◆ theMuonLabel_

◆ theMuonType_

◆ theRootFileName_

◆ tracksToken_