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

Inheritance diagram for TestResolution:

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

	~TestResolution ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

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

virtual void	endJob ()

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

Private Attributes
int	eventCounter_

TFile *	outputFile_

std::auto_ptr< ResolutionFunction >	resolutionFunction_

TProfile *	sigmaPt_

edm::InputTag	theMuonLabel_

int	theMuonType_

std::string	theRootFileName_

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

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

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

Definition at line 54 of file TestResolution.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 20 of file TestResolution.cc.

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

                                                              :
   theMuonLabel_( iConfig.getParameter<edm::InputTag>( "MuonLabel" ) ),
   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;
 }

TestResolution::~TestResolution ( )

Definition at line 36 of file TestResolution.cc.

References MultipleCompare::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

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 63 of file TestResolution.cc.

References gather_cfg::cout, eventCounter_, fillMuonCollection(), edm::Event::getByLabel(), patZpeak::muons, resolutionFunction_, sigmaPt_, theMuonLabel_, theMuonType_, and testEve_cfg::tracks.

                                                                                 {
   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.getByLabel (theMuonLabel_, glbMuons);
     muons = fillMuonCollection(*glbMuons);
   }
   else if (theMuonType_==2) { // StandaloneMuons
     Handle<reco::TrackCollection> saMuons;
     iEvent.getByLabel (theMuonLabel_, saMuons);
     muons = fillMuonCollection(*saMuons);
   }
   else if (theMuonType_==3) { // Tracker tracks
     Handle<reco::TrackCollection> tracks;
     iEvent.getByLabel (theMuonLabel_, 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), recMuon->pt());
 
   }
 }

virtual void TestResolution::endJob ( void )

inlineprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file TestResolution.h.

61 {};

template<typename T >

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

inlineprivate

Definition at line 63 of file TestResolution.h.

References metsig::muon, patZpeak::muons, and mathSSE::sqrt().

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