#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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > 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 &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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

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

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

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

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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.

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

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

◆ ~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.

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

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

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

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

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::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