#include <ComparisonPlots/HCALGPUAnalyzer/plugins/HCALGPUAnalyzer.cc>

Inheritance diagram for HCALGPUAnalyzer:

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

	~HCALGPUAnalyzer () override=default

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

void	beginJob () override

void	endJob () override

Private Attributes
int	depth

int	depthGPU

edm::Service< TFileService >	FileService

TH2F *	hEnergy_2dM0

TH2F *	hEnergy_2dMahi

TH1F *	hEnergy_cpu

TH1F *	hEnergy_cpugpu

TH1F *	hEnergy_cpugpu_rel

TH1F *	hEnergy_gpu

TH1F *	hEnergyM0_cpu

TH1F *	hEnergyM0_gpu

edm::EDGetTokenT< HBHERecHitCollection >	hRhToken

edm::EDGetTokenT< HBHERecHitCollection >	hRhTokenGPU

TH2F *	hTime_2dMahi

TH1F *	hTime_cpu

TH1F *	hTime_gpu

double	iEta

double	iEtaGPU

double	iPhi

double	iPhiGPU

TH2F *	Matched

double	Method0Energy

double	Method0EnergyGPU

double	RecHitEnergy

double	RecHitEnergyGPU

double	RecHitTime

double	RecHitTimeGPU

TH2F *	Unmatched

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 54 of file HCALGPUAnalyzer.cc.

Constructor & Destructor Documentation

◆ HCALGPUAnalyzer()

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

explicit

Definition at line 110 of file HCALGPUAnalyzer.cc.

References edm::ParameterSet::getUntrackedParameter(), and reco::Matched.

                                                                {
   usesResource("TFileService");
 
   hRhToken = consumes<HBHERecHitCollection>(iConfig.getUntrackedParameter<string>("HBHERecHits", "hbheprereco"));
   hRhTokenGPU = consumes<HBHERecHitCollection>(
       iConfig.getUntrackedParameter<string>("HBHERecHits", "hcalCPURecHitsProducer:recHitsLegacyHBHE"));
 
   //
 
   hEnergy_2dM0 = FileService->make<TH2F>("hEnergy_2dM0", "hEnergy_2dM0", 1000, 0., 100., 1000, 0., 100.);
   hEnergy_2dM0->GetXaxis()->SetTitle("Cpu M0 Energy");
   hEnergy_2dM0->GetYaxis()->SetTitle("GPU M0 Energy");
 
   hEnergy_2dMahi = FileService->make<TH2F>("hEnergy_2dMahi", "hEnergy_2dMahi", 1000, 0., 100., 1000, 0., 100.);
   hEnergy_2dMahi->GetXaxis()->SetTitle("CPU Energy");
   hEnergy_2dMahi->GetYaxis()->SetTitle("GPU Energy");
 
   hTime_2dMahi = FileService->make<TH2F>("hTime_2dMahi", "hTime_2dMahi", 250, -12.5, 12.5, 250, -12.5, 12.5);
   hTime_2dMahi->GetXaxis()->SetTitle("Mahi Time CPU");
   hTime_2dMahi->GetYaxis()->SetTitle("Mahi Time GPU");
 
   //
 
   hEnergyM0_cpu = FileService->make<TH1F>("hEnergyM0_cpu", "hEnergyM0_cpu", 100, 0., 100.);
   hEnergyM0_cpu->GetXaxis()->SetTitle("CPU Energy");
 
   hEnergy_cpu = FileService->make<TH1F>("hEnergy_cpu", "hEnergy_cpu", 50, 0., 50.);
   hEnergy_cpu->GetXaxis()->SetTitle("CPU Energy");
 
   hEnergy_gpu = FileService->make<TH1F>("hEnergy_gpu", "hEnergy_gpu", 50, 0., 50.);
   hEnergy_gpu->GetXaxis()->SetTitle("GPU Energy");
 
   //
 
   hEnergy_cpugpu = FileService->make<TH1F>("hEnergy_cpugpu", "hEnergy_cpugpu", 500, -2.5, 2.5);
   hEnergy_cpugpu->GetXaxis()->SetTitle("GPU Energy - CPU Energy [GeV]");
   hEnergy_cpugpu->GetYaxis()->SetTitle("# RecHits");
 
   hEnergy_cpugpu_rel =
       FileService->make<TH1F>("hEnergy_cpugpu_rel", "hEnergy_cpugpu_rel ( E > 0.005 GeV)", 500, -2.5, 2.5);
   hEnergy_cpugpu_rel->GetXaxis()->SetTitle("(GPU Energy - CPU Energy) / CPU energy");
   hEnergy_cpugpu_rel->GetYaxis()->SetTitle("# RecHits");
 
   //
 
   hTime_cpu = FileService->make<TH1F>("hTime_cpu", "hTime_cpu", 50, -25., 25.);
   hTime_cpu->GetXaxis()->SetTitle("CPU Time");
 
   hTime_gpu = FileService->make<TH1F>("hTime_gpu", "hTime_gpu", 50, -25., 25.);
   hTime_gpu->GetXaxis()->SetTitle("GPU Time");
 
   Unmatched = FileService->make<TH2F>("Unmatched", "Unmatched (eta,phi)", 100, -50., 50., 85, 0., 85.);
   Matched = FileService->make<TH2F>("Matched", "Matched (eta,phi)", 100, -50., 50., 85, 0., 85.);
 
   //now do what ever initialization is needed
 }

◆ ~HCALGPUAnalyzer()

HCALGPUAnalyzer::~HCALGPUAnalyzer ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 172 of file HCALGPUAnalyzer.cc.

References hcalRecHitTable_cff::depth, hcalRecHitTable_cff::energy, mps_fire::i, l1tPhase2CaloJetEmulator_cfi::iEta, HcalDetId::ieta(), iEvent, createfilelist::int, HcalDetId::iphi(), dqmiolumiharvest::j, reco::Matched, DetId::rawId(), edm::SortedCollection< T, SORT >::size(), and hcalRecHitTable_cff::time.

                                                                                  {
   using namespace edm;
 
   // Read events
   Handle<HBHERecHitCollection> hRecHits;
   iEvent.getByToken(hRhToken, hRecHits);
 
   Handle<HBHERecHitCollection> hRecHitsGPU;
   iEvent.getByToken(hRhTokenGPU, hRecHitsGPU);
 
   // Loop over all rechits in one event
   for (int i = 0; i < (int)hRecHits->size(); i++) {
     // get ID information for the reconstructed hit
     HcalDetId detID_rh = (*hRecHits)[i].id().rawId();
 
     // ID information can get us detector coordinates
     depth = (*hRecHits)[i].id().depth();
     iEta = detID_rh.ieta();
     iPhi = detID_rh.iphi();
 
     // get some variables
     Method0Energy = (*hRecHits)[i].eraw();
     RecHitEnergy = (*hRecHits)[i].energy();
     RecHitTime = (*hRecHits)[i].time();
 
     hEnergy_cpu->Fill(RecHitEnergy);
     hTime_cpu->Fill(RecHitTime);
 
     /*
      cout << "Run " << i << ": ";
      cout << "Method0Energy: " << Method0Energy;
      cout << "RecHitEnergy: " << RecHitEnergy;
      cout << "depth: " << depth;
      cout << "iEta: " << iEta;
      cout << "iPhi: " << iPhi;
      cout << "RecHitTime" << RecHitTime;
      */
   }
 
   for (int i = 0; i < (int)hRecHitsGPU->size(); i++) {
     // get ID information for the reconstructed hit
     HcalDetId detID_rh = (*hRecHitsGPU)[i].id().rawId();
 
     // ID information can get us detector coordinates
     depthGPU = (*hRecHitsGPU)[i].id().depth();
     iEtaGPU = detID_rh.ieta();
     iPhiGPU = detID_rh.iphi();
 
     // get some variables
     Method0EnergyGPU = (*hRecHitsGPU)[i].eraw();
     RecHitEnergyGPU = (*hRecHitsGPU)[i].energy();
     RecHitTimeGPU = (*hRecHitsGPU)[i].time();
 
     hEnergy_gpu->Fill(RecHitEnergyGPU);
     hTime_gpu->Fill(RecHitTimeGPU);
 
     /*
      cout << "Run " << i << ": ";
      cout << "Method0Energy: " << Method0EnergyGPU;
      cout << "RecHitEnergy: " << RecHitEnergyGPU;
      cout << "depth: " << depthGPU;
      cout << "iEta: " << iEtaGPU;
      cout << "iPhi: " << iPhiGPU;
      cout << "RecHitTime" << RecHitTimeGPU;
      */
   }
 
   // Loop over all rechits in one event
   for (int i = 0; i < (int)hRecHits->size(); i++) {
     HcalDetId detID_rh = (*hRecHits)[i].id().rawId();
 
     bool unmatched = true;
     //     cout << "--------------------------------------------------------" << endl;
 
     for (int j = 0; j < (int)hRecHitsGPU->size(); j++) {
       HcalDetId detID_gpu = (*hRecHitsGPU)[j].id().rawId();
 
       if ((detID_rh == detID_gpu)) {
         /*
      cout << "Mtime(cpu)" << (*hRecHits)[i].time() << endl; 
      cout << "     Mtime(gpu)" << (*hRecHitsGPU)[j].time() << endl;
 
      cout << "M0E(cpu)" << (*hRecHits)[i].eraw() << endl; 
      cout << "     M0E(gpu)" << (*hRecHitsGPU)[j].eraw() << endl;
      */
 
         auto relValue = ((*hRecHitsGPU)[j].energy() - (*hRecHits)[i].energy()) / (*hRecHits)[i].energy();
 
         hEnergy_2dM0->Fill((*hRecHits)[i].eraw(), (*hRecHitsGPU)[j].eraw());
         hEnergy_2dMahi->Fill((*hRecHits)[i].energy(), (*hRecHitsGPU)[j].energy());
         hEnergy_cpugpu->Fill((*hRecHitsGPU)[j].energy() - (*hRecHits)[i].energy());
         if ((*hRecHits)[i].energy() > 0.005)
           hEnergy_cpugpu_rel->Fill(relValue);
         hTime_2dMahi->Fill((*hRecHits)[i].time(), (*hRecHitsGPU)[j].time());
 
         /*
      if((relValue < - 0.9) and ((*hRecHits)[i].energy()>0.005)) {
        cout << "----------------------------------"<< endl;
        cout << " detID = " << detID_rh.rawId() << endl;
        cout << "ME(cpu)" << (*hRecHits)[i].energy() << endl; 
        cout << "     ME(gpu)" << (*hRecHitsGPU)[j].energy() << endl;
      }
      */
 
         Matched->Fill(detID_rh.ieta(), detID_rh.iphi());
 
         unmatched = false;
       }
     }
 
 
     if (unmatched) {
       Unmatched->Fill(detID_rh.ieta(), detID_rh.iphi());
       //       cout << "   recHit not matched ="  << detID_rh << "  E(raw)=" << (*hRecHits)[i].eraw() << " E=" << (*hRecHits)[i].energy() << endl;
     }
   }
 }

◆ beginJob()

void HCALGPUAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 292 of file HCALGPUAnalyzer.cc.

292 {}

◆ endJob()

void HCALGPUAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 295 of file HCALGPUAnalyzer.cc.

295 {}

◆ fillDescriptions()

void HCALGPUAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 298 of file HCALGPUAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                  {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

Member Data Documentation

◆ depth

int HCALGPUAnalyzer::depth

private

Definition at line 75 of file HCALGPUAnalyzer.cc.

Referenced by dqm_interfaces.DirID::__eq__(), dqm_interfaces.DirID::__repr__(), and dqm_interfaces.DirWalkerDB::run().

◆ depthGPU

int HCALGPUAnalyzer::depthGPU

private

Definition at line 75 of file HCALGPUAnalyzer.cc.

◆ FileService

edm::Service<TFileService> HCALGPUAnalyzer::FileService

private

Definition at line 93 of file HCALGPUAnalyzer.cc.

◆ hEnergy_2dM0

TH2F* HCALGPUAnalyzer::hEnergy_2dM0

private

Definition at line 78 of file HCALGPUAnalyzer.cc.

◆ hEnergy_2dMahi

TH2F* HCALGPUAnalyzer::hEnergy_2dMahi

private

Definition at line 77 of file HCALGPUAnalyzer.cc.

◆ hEnergy_cpu

TH1F* HCALGPUAnalyzer::hEnergy_cpu

private

Definition at line 83 of file HCALGPUAnalyzer.cc.

◆ hEnergy_cpugpu

TH1F* HCALGPUAnalyzer::hEnergy_cpugpu

private

Definition at line 85 of file HCALGPUAnalyzer.cc.

◆ hEnergy_cpugpu_rel

TH1F* HCALGPUAnalyzer::hEnergy_cpugpu_rel

private

Definition at line 86 of file HCALGPUAnalyzer.cc.

◆ hEnergy_gpu

TH1F* HCALGPUAnalyzer::hEnergy_gpu

private

Definition at line 84 of file HCALGPUAnalyzer.cc.

◆ hEnergyM0_cpu

TH1F* HCALGPUAnalyzer::hEnergyM0_cpu

private

Definition at line 87 of file HCALGPUAnalyzer.cc.

◆ hEnergyM0_gpu

TH1F* HCALGPUAnalyzer::hEnergyM0_gpu

private

Definition at line 88 of file HCALGPUAnalyzer.cc.

◆ hRhToken

edm::EDGetTokenT<HBHERecHitCollection> HCALGPUAnalyzer::hRhToken

private

Definition at line 95 of file HCALGPUAnalyzer.cc.

◆ hRhTokenGPU

edm::EDGetTokenT<HBHERecHitCollection> HCALGPUAnalyzer::hRhTokenGPU

private

Definition at line 96 of file HCALGPUAnalyzer.cc.

◆ hTime_2dMahi

TH2F* HCALGPUAnalyzer::hTime_2dMahi

private

Definition at line 79 of file HCALGPUAnalyzer.cc.

◆ hTime_cpu

TH1F* HCALGPUAnalyzer::hTime_cpu

private

Definition at line 89 of file HCALGPUAnalyzer.cc.

◆ hTime_gpu

TH1F* HCALGPUAnalyzer::hTime_gpu

private

Definition at line 90 of file HCALGPUAnalyzer.cc.

◆ iEta

double HCALGPUAnalyzer::iEta

private

Definition at line 73 of file HCALGPUAnalyzer.cc.

◆ iEtaGPU

double HCALGPUAnalyzer::iEtaGPU

private

Definition at line 73 of file HCALGPUAnalyzer.cc.

◆ iPhi

double HCALGPUAnalyzer::iPhi

private

Definition at line 74 of file HCALGPUAnalyzer.cc.

◆ iPhiGPU

double HCALGPUAnalyzer::iPhiGPU

private

Definition at line 74 of file HCALGPUAnalyzer.cc.

◆ Matched

TH2F* HCALGPUAnalyzer::Matched

private

Definition at line 82 of file HCALGPUAnalyzer.cc.

◆ Method0Energy

double HCALGPUAnalyzer::Method0Energy

private

Definition at line 70 of file HCALGPUAnalyzer.cc.

◆ Method0EnergyGPU

double HCALGPUAnalyzer::Method0EnergyGPU

private

Definition at line 70 of file HCALGPUAnalyzer.cc.

◆ RecHitEnergy

double HCALGPUAnalyzer::RecHitEnergy

private

Definition at line 71 of file HCALGPUAnalyzer.cc.

◆ RecHitEnergyGPU

double HCALGPUAnalyzer::RecHitEnergyGPU

private

Definition at line 71 of file HCALGPUAnalyzer.cc.

◆ RecHitTime

double HCALGPUAnalyzer::RecHitTime

private

Definition at line 72 of file HCALGPUAnalyzer.cc.

◆ RecHitTimeGPU

double HCALGPUAnalyzer::RecHitTimeGPU

private

Definition at line 72 of file HCALGPUAnalyzer.cc.

◆ Unmatched

TH2F* HCALGPUAnalyzer::Unmatched

private

Definition at line 81 of file HCALGPUAnalyzer.cc.

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HCALGPUAnalyzer()

◆ ~HCALGPUAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ fillDescriptions()

Member Data Documentation

◆ depth

◆ depthGPU

◆ FileService

◆ hEnergy_2dM0

◆ hEnergy_2dMahi

◆ hEnergy_cpu

◆ hEnergy_cpugpu

◆ hEnergy_cpugpu_rel

◆ hEnergy_gpu

◆ hEnergyM0_cpu

◆ hEnergyM0_gpu

◆ hRhToken

◆ hRhTokenGPU

◆ hTime_2dMahi

◆ hTime_cpu

◆ hTime_gpu

◆ iEta

◆ iEtaGPU

◆ iPhi

◆ iPhiGPU

◆ Matched

◆ Method0Energy

◆ Method0EnergyGPU

◆ RecHitEnergy

◆ RecHitEnergyGPU

◆ RecHitTime

◆ RecHitTimeGPU

◆ Unmatched