#include <PFTester.h>

Inheritance diagram for PFTester:

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

void	beginJob () override

void	endJob () override

	PFTester (const edm::ParameterSet &)

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

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

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

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
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_.empty())
     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 ( )

override

Definition at line 52 of file PFTester.cc.

52 { }

Member Function Documentation

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

override

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(), packedPFCandidateRefMixer_cfi::pf, 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 )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 55 of file PFTester.cc.

References dbe_, M_PI, and Utilities::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 )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 214 of file PFTester.cc.

References dbe_.

 {
 
   // Store the DAQ Histograms
   if (!outputFile_.empty() && 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