#include <PFTester.h>

Inheritance diagram for PFTester:

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

virtual void	beginJob ()

virtual void	endJob ()

	PFTester (const edm::ParameterSet &)

virtual	~PFTester ()

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
DQMStore *	dbe_

edm::EDGetTokenT < reco::PFCandidateCollection >	inputPFlowLabel_tok_

std::map< std::string, MonitorElement * >	me

std::string	outputFile_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

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 19 of file PFTester.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 40 of file PFTester.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
 
   inputPFlowLabel_tok_         = consumes<reco::PFCandidateCollection> (iConfig.getParameter<std::string>("InputPFlowLabel") );
   outputFile_                  = iConfig.getUntrackedParameter<std::string>("OutputFile");
 
   if (outputFile_.size() > 0)
     edm::LogInfo("OutputInfo") << " ParticleFLow Task histograms will be saved to '" << outputFile_.c_str() << "'";
   else edm::LogInfo("OutputInfo") << " ParticleFlow Task histograms will NOT be saved";
 
 }

PFTester::~PFTester ( )

virtual

Definition at line 52 of file PFTester.cc.

52 { }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 95 of file PFTester.cc.

References reco::LeafCandidate::charge(), reco::LeafCandidate::et(), reco::LeafCandidate::eta(), edm::Event::getByToken(), reco::PFCandidate::particleId(), reco::LeafCandidate::pdgId(), reco::LeafCandidate::phi(), and edm::Handle< T >::product().

 {
   
   // Data to Retrieve from the Event
   const PFCandidateCollection *pflow_candidates;
 
   // ==========================================================
   // Retrieve!
   // ==========================================================
 
   { 
 
     // Get Particle Flow Candidates
     Handle<PFCandidateCollection> pflow_hnd;
     iEvent.getByToken(inputPFlowLabel_tok_, pflow_hnd);
     pflow_candidates = pflow_hnd.product();
 
   }
 
   if (!pflow_candidates) {
 
     edm::LogInfo("OutputInfo") << " failed to retrieve data required by ParticleFlow Task";
     edm::LogInfo("OutputInfo") << " ParticleFlow Task cannot continue...!";
     return;
 
   }
 
   // ==========================================================
   // Analyze!
   // ==========================================================
 
   // Loop Over Particle Flow Candidates
   PFCandidateCollection::const_iterator pf;
   for (pf = pflow_candidates->begin(); pf != pflow_candidates->end(); pf++) {
 
     const PFCandidate *particle = &(*pf);
 
     // Fill Histograms for Candidate Methods
     me["CandidateEt"]->Fill(particle->et());
     me["CandidateEta"]->Fill(particle->eta());
     me["CandidatePhi"]->Fill(particle->phi());
     me["CandidateCharge"]->Fill(particle->charge());
     me["CandidatePdgId"]->Fill(particle->pdgId());
 
     // Fill Histograms for PFCandidate Specific Methods
     me["PFCandidateType"]->Fill(particle->particleId());
     // particle->elementsInBlocks(); 
     
    // Get the PFBlock and Elements
    // JW: Returns vector of blocks now ,TO BE FIXED ----
     /*PFBlock block = *(particle->block());
     OwnVector<PFBlockElement> elements = block.elements();
     int numElements = elements.size();
     int numTrackElements = 0;
     int numPS1Elements = 0;
     int numPS2Elements = 0;
     int numECALElements = 0;
     int numHCALElements = 0;
     int numMuonElements = 0;
 
     // Loop over Elements in Block
     OwnVector<PFBlockElement>::const_iterator element;
     for (element = elements.begin(); element != elements.end(); element++) {
 
       int element_type = element->type();
       // Element is a Tracker Track
       if (element_type == PFBlockElement::TRACK) {
 
         // Get General Information about the Track
         PFRecTrack track = *(element->trackRefPF());
         me["TrackCharge"]->Fill(track.charge());
         me["TrackNumPoints"]->Fill(track.nTrajectoryPoints());
         me["TrackNumMeasurements"]->Fill(track.nTrajectoryMeasurements());
 
         // Loop Over Points in the Track
         vector<PFTrajectoryPoint> points = track.trajectoryPoints();
         vector<PFTrajectoryPoint>::iterator point;
         for (point = points.begin(); point != points.end(); point++) {
           int point_layer = point->layer();
           double x = point->positionXYZ().x();
           double y = point->positionXYZ().y();
           double z = point->positionXYZ().z();
           //switch (point_layer) {
           //case PFTrajectoryPoint::ClosestApproach:
           // Fill the Track's D0
           if (point_layer == PFTrajectoryPoint::ClosestApproach) {
             me["TrackImpactParameter"]->Fill(sqrt(x*x + y*y + z*z));
           }
         }
         numTrackElements++;
       }
 
       // Element is an ECAL Cluster
       else if (element_type == PFBlockElement::ECAL) {
         numECALElements++;
       }
       // Element is a HCAL Cluster
       else if (element_type == PFBlockElement::HCAL) {
         numHCALElements++;
       }
       // Element is a Muon Track
       else if (element_type == PFBlockElement::MUON) {
         numMuonElements++;
       } 
       // Fill the Respective Elements Sizes
       me["NumElements"]->Fill(numElements);
       me["NumTrackElements"]->Fill(numTrackElements);
       me["NumPS1Elements"]->Fill(numPS1Elements);
       me["NumPS2Elements"]->Fill(numPS2Elements);
       me["NumECALElements"]->Fill(numECALElements);
       me["NumHCALElements"]->Fill(numHCALElements);
       me["NumMuonElements"]->Fill(numMuonElements);
     } ----------------------------------------------  */
 
   }
 
 }

void PFTester::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 55 of file PFTester.cc.

References dbe_, M_PI, and cppFunctionSkipper::operator.

 {
 
   // get ahold of back-end interface
   dbe_ = edm::Service<DQMStore>().operator->();
   
   if (dbe_) {
 
     dbe_->setCurrentFolder("PFTask/PFCandidates");
 
     me["CandidateEt"] = dbe_->book1D("CandidateEt","CandidateEt",1000,0,1000);
     me["CandidateEta"] = dbe_->book1D("CandidateEta","CandidateEta",200,-5,5);
     me["CandidatePhi"] = dbe_->book1D("CandidatePhi","CandidatePhi",200,-M_PI,M_PI);
     me["CandidateCharge"] = dbe_->book1D("CandidateCharge","CandidateCharge",5,-2,2);
     me["PFCandidateType"] = dbe_->book1D("PFCandidateType","PFCandidateType",10,0,10);
 
     dbe_->setCurrentFolder("PFTask/PFBlocks");
 
     me["NumElements"] = dbe_->book1D("NumElements","NumElements",25,0,25);
     me["NumTrackElements"] = dbe_->book1D("NumTrackElements","NumTrackElements",5,0,5);
     me["NumPS1Elements"] = dbe_->book1D("NumPS1Elements","NumPS1Elements",5,0,5);
     me["NumPS2Elements"] = dbe_->book1D("NumPS2Elements","NumPS2Elements",5,0,5);
     me["NumECALElements"] = dbe_->book1D("NumECALElements","NumECALElements",5,0,5);
     me["NumHCALElements"] = dbe_->book1D("NumHCALElements","NumHCALElements",5,0,5);
     me["NumMuonElements"] = dbe_->book1D("NumMuonElements","NumMuonElements",5,0,5);
 
     dbe_->setCurrentFolder("PFTask/PFTracks");
 
     me["TrackCharge"] = dbe_->book1D("TrackCharge","TrackCharge",5,-2,2);
     me["TrackNumPoints"] = dbe_->book1D("TrackNumPoints","TrackNumPoints",100,0,100);
     me["TrackNumMeasurements"] = dbe_->book1D("TrackNumMeasurements","TrackNumMeasurements",100,0,100);
     me["TrackImpactParameter"] = dbe_->book1D("TrackImpactParameter","TrackImpactParameter",1000,0,1);
 
     dbe_->setCurrentFolder("PFTask/PFClusters");
 
   }
 
 }

void PFTester::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 214 of file PFTester.cc.

References dbe_.

 {
 
   // Store the DAQ Histograms
   if (outputFile_.size() > 0 && dbe_)
     dbe_->save(outputFile_);
 
 }

Member Data Documentation

DQMStore* PFTester::dbe_

private

Definition at line 32 of file PFTester.h.

edm::EDGetTokenT<reco::PFCandidateCollection> PFTester::inputPFlowLabel_tok_

private

Definition at line 37 of file PFTester.h.

std::map<std::string, MonitorElement*> PFTester::me

private

Definition at line 33 of file PFTester.h.

std::string PFTester::outputFile_

private

Definition at line 36 of file PFTester.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation