#include <L1Trigger/L1GlobalCaloTrigger/plugins/L1GctValidation.cc>

Inheritance diagram for L1GctValidation:

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

	~L1GctValidation () override

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

	EDConsumerBase (EDConsumerBase &&)=default

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

void	beginJob () override

void	endJob () override

Private Attributes
edm::InputTag	m_energy_tag

edm::InputTag	m_gctinp_tag

edm::ESGetToken< L1CaloEtScale, L1HfRingEtScaleRcd >	m_hfRingEtScaleToken

edm::ESGetToken< L1CaloEtScale, L1HtMissScaleRcd >	m_htMissScaleToken

edm::ESGetToken< L1GctJetFinderParams, L1GctJetFinderParamsRcd >	m_jfParsToken

TH2F *	theDPhiVsMissEt

TH2F *	theDPhiVsMissHt

TH1F *	theHfRing0CountNegativeEta

TH1F *	theHfRing0CountPositiveEta

TH1F *	theHfRing0EtSumNegativeEta

TH1F *	theHfRing0EtSumPositiveEta

TH1F *	theHfRing1CountNegativeEta

TH1F *	theHfRing1CountPositiveEta

TH1F *	theHfRing1EtSumNegativeEta

TH1F *	theHfRing1EtSumPositiveEta

TH2F *	theHtVsInternalJetsSum

TH1F *	theMissEtAngle

TH2F *	theMissEtAngleVsInputRegions

TH1F *	theMissEtInGeV

TH1F *	theMissEtInLsb

TH2F *	theMissEtMagVsInputRegions

TH2F *	theMissEtVector

TH2F *	theMissEtVsMissHt

TH2F *	theMissEtVsMissHtAngle

TH1F *	theMissHtAngle

TH1F *	theMissHtInGeV

TH1F *	theMissHtInLsb

TH2F *	theMissHtMagVsInputRegions

TH2F *	theMissHtPhiVsInternalJetsSum

TH2F *	theMissHtVector

TH2F *	theMissHtVsInternalJetsSum

TH2F *	theMissHxVsInternalJetsSum

TH2F *	theMissHyVsInternalJetsSum

TH1F *	theSumEtInGeV

TH1F *	theSumEtInLsb

TH2F *	theSumEtVsInputRegions

TH1F *	theSumHtInGeV

TH1F *	theSumHtInLsb

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 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<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: produces standard plots of Gct output quantities to enable validation of global event quantities in particular

Definition at line 43 of file L1GctValidation.h.

Constructor & Destructor Documentation

◆ L1GctValidation()

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

explicit

Definition at line 16 of file L1GctValidation.cc.

References TFileService::kSharedResource.

     : m_gctinp_tag(iConfig.getUntrackedParameter<edm::InputTag>("rctInputTag", edm::InputTag("rctDigis"))),
       m_energy_tag(iConfig.getUntrackedParameter<edm::InputTag>("gctInputTag", edm::InputTag("gctDigis"))),
       m_jfParsToken(esConsumes<L1GctJetFinderParams, L1GctJetFinderParamsRcd>()),
       m_htMissScaleToken(esConsumes<L1CaloEtScale, L1HtMissScaleRcd>()),
       m_hfRingEtScaleToken(esConsumes<L1CaloEtScale, L1HfRingEtScaleRcd>()) {
   usesResource(TFileService::kSharedResource);
 }

◆ ~L1GctValidation()

L1GctValidation::~L1GctValidation ( )

override

Definition at line 25 of file L1GctValidation.cc.

                                   {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 35 of file L1GctValidation.cc.

References funct::cos(), SiPixelRawToDigiRegional_cfi::deltaPhi, L1CaloEtScale::et(), edm::EventSetup::getHandle(), L1GctJetFinderParams::getHtJetEtThresholdGeV(), L1GctJetFinderParams::getHtLsbGeV(), L1GctJetFinderParams::getMHtJetEtThresholdGeV(), L1GctJetFinderParams::getRgnEtLsbGeV(), iEvent, HLT_2023v12_cff::inputColl, metsig::jet, m_energy_tag, m_gctinp_tag, m_hfRingEtScaleToken, m_htMissScaleToken, m_jfParsToken, M_PI, edm::ESHandle< T >::product(), funct::sin(), mathSSE::sqrt(), theDPhiVsMissEt, theDPhiVsMissHt, theHfRing0CountNegativeEta, theHfRing0CountPositiveEta, theHfRing0EtSumNegativeEta, theHfRing0EtSumPositiveEta, theHfRing1CountNegativeEta, theHfRing1CountPositiveEta, theHfRing1EtSumNegativeEta, theHfRing1EtSumPositiveEta, theHtVsInternalJetsSum, theMissEtAngle, theMissEtAngleVsInputRegions, theMissEtInGeV, theMissEtInLsb, theMissEtMagVsInputRegions, theMissEtVector, theMissEtVsMissHt, theMissEtVsMissHtAngle, theMissHtAngle, theMissHtInGeV, theMissHtInLsb, theMissHtMagVsInputRegions, theMissHtPhiVsInternalJetsSum, theMissHtVector, theMissHtVsInternalJetsSum, theMissHxVsInternalJetsSum, theMissHyVsInternalJetsSum, theSumEtInGeV, theSumEtInLsb, theSumEtVsInputRegions, theSumHtInGeV, and theSumHtInLsb.

                                                                                  {
   using namespace edm;
 
   // Get the scales from the event setup
   ESHandle<L1GctJetFinderParams> jfPars = iSetup.getHandle(m_jfParsToken);
 
   double lsbForEt = jfPars.product()->getRgnEtLsbGeV();
   double lsbForHt = jfPars.product()->getHtLsbGeV();
 
   unsigned httJetThreshold = static_cast<int>(jfPars.product()->getHtJetEtThresholdGeV() / lsbForHt);
   unsigned htmJetThreshold = static_cast<int>(jfPars.product()->getMHtJetEtThresholdGeV() / lsbForHt);
 
   ESHandle<L1CaloEtScale> htMissScale = iSetup.getHandle(m_htMissScaleToken);
   ESHandle<L1CaloEtScale> hfRingEtScale = iSetup.getHandle(m_hfRingEtScaleToken);
 
   // Get the Gct energy sums from the event
   Handle<L1GctEtTotalCollection> sumEtColl;
   iEvent.getByLabel(m_energy_tag, sumEtColl);
   Handle<L1GctEtHadCollection> sumHtColl;
   iEvent.getByLabel(m_energy_tag, sumHtColl);
   Handle<L1GctEtMissCollection> missEtColl;
   iEvent.getByLabel(m_energy_tag, missEtColl);
   Handle<L1GctHtMissCollection> missHtColl;
   iEvent.getByLabel(m_energy_tag, missHtColl);
 
   // Get the input calo regions from the event (for checking MEt)
   Handle<L1CaloRegionCollection> inputColl;
   iEvent.getByLabel(m_gctinp_tag, inputColl);
 
   // Get the internal jet data from the event (for checking Ht)
   Handle<L1GctInternJetDataCollection> internalJetsColl;
   iEvent.getByLabel(m_energy_tag, internalJetsColl);
 
   double etTot = 0.0;
   for (L1GctEtTotalCollection::const_iterator jbx = sumEtColl->begin(); jbx != sumEtColl->end(); jbx++) {
     if (jbx->bx() == 0) {
       etTot = static_cast<double>(jbx->et());
     }
   }
 
   double etHad = 0.0;
   for (L1GctEtHadCollection::const_iterator jbx = sumHtColl->begin(); jbx != sumHtColl->end(); jbx++) {
     if (jbx->bx() == 0) {
       etHad = static_cast<double>(jbx->et());
     }
   }
 
   double etMiss = 0.0;
   double etMAng = 0.0;
   for (L1GctEtMissCollection::const_iterator jbx = missEtColl->begin(); jbx != missEtColl->end(); jbx++) {
     if (jbx->bx() == 0) {
       etMiss = static_cast<double>(jbx->et());
       int phibin = jbx->phi();
       if (phibin >= 36)
         phibin -= 72;
       double etMPhi = static_cast<double>(phibin);
 
       etMAng = (etMPhi + 0.5) * M_PI / 36.;
     }
   }
 
   double etTotFromRegions = 0.0;
   double exTotFromRegions = 0.0;
   double eyTotFromRegions = 0.0;
   for (L1CaloRegionCollection::const_iterator jrg = inputColl->begin(); jrg != inputColl->end(); jrg++) {
     if (jrg->bx() == 0) {
       double rgEt = static_cast<double>(jrg->et()) * lsbForEt;
       double rgPhibin = static_cast<double>(jrg->id().iphi());
       double rgPh = (rgPhibin + 0.5) * M_PI / 9.;
 
       etTotFromRegions += rgEt;
       exTotFromRegions += rgEt * cos(rgPh);
       eyTotFromRegions += rgEt * sin(rgPh);
     }
   }
 
   double htMissGct = 0.0;
   double htMissAng = 0.0;
   double htMissGeV = 0.0;
   for (L1GctHtMissCollection::const_iterator jbx = missHtColl->begin(); jbx != missHtColl->end(); jbx++) {
     if (jbx->bx() == 0) {
       htMissGct = static_cast<double>(jbx->et());
       htMissGeV = htMissScale->et(jbx->et());
       int phibin = jbx->phi();
       if (phibin >= 9)
         phibin -= 18;
       double htMPhi = static_cast<double>(phibin);
       htMissAng = (htMPhi + 0.5) * M_PI / 9.;
     }
   }
 
   double htFromJets = 0.0;
   double hxFromJets = 0.0;
   double hyFromJets = 0.0;
   for (L1GctInternJetDataCollection::const_iterator jet = internalJetsColl->begin(); jet != internalJetsColl->end();
        jet++) {
     if (jet->bx() == 0 && !jet->empty()) {
       unsigned jetEtGct = jet->et();
       double jetEt = static_cast<double>(jetEtGct);
       int phibin = jet->regionId().iphi();
       if (phibin >= 9)
         phibin -= 18;
       // The phi bin centres are at 0, 20, 40, ... degrees
       double jetAng = (static_cast<double>(phibin)) * M_PI / 9.;
       if (jetEtGct > httJetThreshold) {
         htFromJets += jetEt;
       }
       if (jetEtGct > htmJetThreshold) {
         hxFromJets += jetEt * cos(jetAng);
         hyFromJets += jetEt * sin(jetAng);
       }
     }
   }
 
   double dPhiMetMht = deltaPhi(etMAng, htMissAng);
 
   theSumEtInLsb->Fill(etTot);
   theSumHtInLsb->Fill(etHad);
   theMissEtInLsb->Fill(etMiss);
   theMissHtInLsb->Fill(htMissGct);
   theSumEtInGeV->Fill(etTot * lsbForEt);
   theSumHtInGeV->Fill(etHad * lsbForHt);
   theMissEtInGeV->Fill(etMiss * lsbForEt);
   theMissEtAngle->Fill(etMAng);
   theMissEtVector->Fill(etMiss * lsbForEt * cos(etMAng), etMiss * lsbForEt * sin(etMAng));
   if (htMissGct < 126.5) {
     theMissHtInGeV->Fill(htMissGeV);
     theMissHtAngle->Fill(htMissAng);
     theMissHtVector->Fill(htMissGeV * cos(htMissAng), htMissGeV * sin(htMissAng));
   }
 
   theSumEtVsInputRegions->Fill(etTot * lsbForEt, etTotFromRegions);
   theMissEtMagVsInputRegions->Fill(etMiss * lsbForEt,
                                    sqrt(exTotFromRegions * exTotFromRegions + eyTotFromRegions * eyTotFromRegions));
   theMissEtAngleVsInputRegions->Fill(etMAng, atan2(-eyTotFromRegions, -exTotFromRegions));
   theMissHtMagVsInputRegions->Fill(htMissGeV,
                                    sqrt(exTotFromRegions * exTotFromRegions + eyTotFromRegions * eyTotFromRegions));
 
   theMissEtVsMissHt->Fill(etMiss * lsbForEt, htMissGeV);
   theMissEtVsMissHtAngle->Fill(etMAng, htMissAng);
   theDPhiVsMissEt->Fill(dPhiMetMht, etMiss * lsbForEt);
   theDPhiVsMissHt->Fill(dPhiMetMht, htMissGeV);
 
   theHtVsInternalJetsSum->Fill(etHad * lsbForHt, htFromJets * lsbForHt);
   if (htMissGct < 126.5) {
     theMissHtVsInternalJetsSum->Fill(htMissGeV, sqrt(hxFromJets * hxFromJets + hyFromJets * hyFromJets) * lsbForHt);
     theMissHtPhiVsInternalJetsSum->Fill(htMissAng, atan2(-hyFromJets, -hxFromJets));
     theMissHxVsInternalJetsSum->Fill(htMissGeV * cos(htMissAng), hxFromJets * lsbForHt);
     theMissHyVsInternalJetsSum->Fill(htMissGeV * sin(htMissAng), hyFromJets * lsbForHt);
   }
 
   // Get minbias trigger quantities from HF
   Handle<L1GctHFRingEtSumsCollection> HFEtSumsColl;
   Handle<L1GctHFBitCountsCollection> HFCountsColl;
   iEvent.getByLabel(m_energy_tag, HFEtSumsColl);
   iEvent.getByLabel(m_energy_tag, HFCountsColl);
 
   for (L1GctHFRingEtSumsCollection::const_iterator es = HFEtSumsColl->begin(); es != HFEtSumsColl->end(); es++) {
     if (es->bx() == 0) {
       theHfRing0EtSumPositiveEta->Fill(hfRingEtScale->et(es->etSum(0)));
       theHfRing0EtSumNegativeEta->Fill(hfRingEtScale->et(es->etSum(1)));
       theHfRing1EtSumPositiveEta->Fill(hfRingEtScale->et(es->etSum(2)));
       theHfRing1EtSumNegativeEta->Fill(hfRingEtScale->et(es->etSum(3)));
     }
   }
 
   for (L1GctHFBitCountsCollection::const_iterator bc = HFCountsColl->begin(); bc != HFCountsColl->end(); bc++) {
     if (bc->bx() == 0) {
       theHfRing0CountPositiveEta->Fill(bc->bitCount(0));
       theHfRing0CountNegativeEta->Fill(bc->bitCount(1));
       theHfRing1CountPositiveEta->Fill(bc->bitCount(2));
       theHfRing1CountNegativeEta->Fill(bc->bitCount(3));
     }
   }
 }

◆ beginJob()

void L1GctValidation::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 212 of file L1GctValidation.cc.

References L1TDiffHarvesting_cfi::dir1, compareTotals::fs, M_PI, TFileDirectory::make(), theDPhiVsMissEt, theDPhiVsMissHt, theHfRing0CountNegativeEta, theHfRing0CountPositiveEta, theHfRing0EtSumNegativeEta, theHfRing0EtSumPositiveEta, theHfRing1CountNegativeEta, theHfRing1CountPositiveEta, theHfRing1EtSumNegativeEta, theHfRing1EtSumPositiveEta, theHtVsInternalJetsSum, theMissEtAngle, theMissEtAngleVsInputRegions, theMissEtInGeV, theMissEtInLsb, theMissEtMagVsInputRegions, theMissEtVector, theMissEtVsMissHt, theMissEtVsMissHtAngle, theMissHtAngle, theMissHtInGeV, theMissHtInLsb, theMissHtMagVsInputRegions, theMissHtPhiVsInternalJetsSum, theMissHtVector, theMissHtVsInternalJetsSum, theMissHxVsInternalJetsSum, theMissHyVsInternalJetsSum, theSumEtInGeV, theSumEtInLsb, theSumEtVsInputRegions, theSumHtInGeV, and theSumHtInLsb.

                                {
   edm::Service<TFileService> fs;
 
   TFileDirectory dir0 = fs->mkdir("L1GctEtSums");
 
   theSumEtInLsb = dir0.make<TH1F>("SumEtInLsb", "Total Et (GCT units)", 128, 0., 2048.);
   theSumHtInLsb = dir0.make<TH1F>("SumHtInLsb", "Total Ht (GCT units)", 128, 0., 2048.);
   theMissEtInLsb = dir0.make<TH1F>("MissEtInLsb", "Missing Et magnitude (GCT units)", 128, 0., 1024.);
   theMissHtInLsb = dir0.make<TH1F>("MissHtInLsb", "Missing Ht magnitude (GCT units)", 128, 0., 127.);
   theSumEtInGeV = dir0.make<TH1F>("SumEtInGeV", "Total Et (in GeV)", 100, 0., 1000.);
   theSumHtInGeV = dir0.make<TH1F>("SumHtInGeV", "Total Ht (in GeV)", 100, 0., 1000.);
   theMissEtInGeV = dir0.make<TH1F>("MissEtInGeV", "Missing Et magnitude (in GeV)", 100, 0., 500.);
   theMissEtAngle = dir0.make<TH1F>("MissEtAngle", "Missing Et angle", 72, -M_PI, M_PI);
   theMissEtVector = dir0.make<TH2F>("MissEtVector", "Missing Ex vs Missing Ey", 100, -100., 100., 100, -100., 100.);
   theMissHtInGeV = dir0.make<TH1F>("MissHtInGeV", "Missing Ht magnitude (in GeV)", 100, 0., 500.);
   theMissHtAngle = dir0.make<TH1F>("MissHtAngle", "Missing Ht angle", 72, -M_PI, M_PI);
   theMissHtVector = dir0.make<TH2F>("MissHtVector", "Missing Hx vs Missing Hy", 100, -100., 100., 100, -100., 100.);
   theSumEtVsInputRegions =
       dir0.make<TH2F>("SumEtVsInputRegions", "Total Et vs sum of input regions", 100, 0., 1000., 100, 0., 1000.);
   theMissEtMagVsInputRegions = dir0.make<TH2F>(
       "MissEtMagVsInputRegions", "Missing Et magnitude vs sum of input regions", 100, 0., 500., 100, 0., 500.);
   theMissEtAngleVsInputRegions = dir0.make<TH2F>(
       "MissEtAngleVsInputRegions", "Missing Et angle vs sum of input regions", 72, -M_PI, M_PI, 72, -M_PI, M_PI);
   theMissHtMagVsInputRegions = dir0.make<TH2F>(
       "MissHtMagVsInputRegions", "Missing Ht magnitude vs sum of input regions", 100, 0., 500., 100, 0., 500.);
   theMissEtVsMissHt = dir0.make<TH2F>("MissEtVsMissHt", "Missing Et vs Missing Ht", 100, 0., 500., 100, 0., 500.);
   theMissEtVsMissHtAngle = dir0.make<TH2F>(
       "MissEtVsMissHtAngle", "Angle correlation Missing Et vs Missing Ht", 72, -M_PI, M_PI, 72, -M_PI, M_PI);
   theDPhiVsMissEt =
       dir0.make<TH2F>("theDPhiVsMissEt", "Angle difference MET-MHT vs MET magnitude", 72, -M_PI, M_PI, 100, 0., 500.);
   theDPhiVsMissHt =
       dir0.make<TH2F>("theDPhiVsMissHt", "Angle difference MET-MHT vs MHT magnitude", 72, -M_PI, M_PI, 100, 0., 500.);
 
   theHtVsInternalJetsSum = dir0.make<TH2F>(
       "HtVsInternalJetsSum", "Ht vs scalar sum of jet Et values (in GeV)", 128, 0., 2048., 128, 0., 2048.);
   theMissHtVsInternalJetsSum = dir0.make<TH2F>(
       "MissHtVsInternalJetsSum", "Missing Ht vs vector sum of jet Et values (in GeV)", 128, 0., 512., 128, 0., 512.);
   theMissHtPhiVsInternalJetsSum = dir0.make<TH2F>("MissHtPhiVsInternalJetsSum",
                                                   "Angle correlation Missing Ht vs vector sum of jet Et values",
                                                   72,
                                                   -M_PI,
                                                   M_PI,
                                                   72,
                                                   -M_PI,
                                                   M_PI);
   theMissHxVsInternalJetsSum = dir0.make<TH2F>("MissHxVsInternalJetsSum",
                                                "Missing Ht x component vs sum of jet Et values (in GeV)",
                                                128,
                                                -256.,
                                                256.,
                                                128,
                                                -256.,
                                                256.);
   theMissHyVsInternalJetsSum = dir0.make<TH2F>("MissHyVsInternalJetsSum",
                                                "Missing Ht y component vs sum of jet Et values (in GeV)",
                                                128,
                                                -256.,
                                                256.,
                                                128,
                                                -256.,
                                                256.);
 
   TFileDirectory dir1 = fs->mkdir("L1GctHfSumsAndJetCounts");
 
   // Minimum bias triggers from Hf inner rings
   theHfRing0EtSumPositiveEta = dir1.make<TH1F>("HfRing0EtSumPositiveEta", "Hf Inner Ring0 Et eta+", 60, 0., 30.);
   theHfRing0EtSumNegativeEta = dir1.make<TH1F>("HfRing0EtSumNegativeEta", "Hf Inner Ring0 Et eta-", 60, 0., 30.);
   theHfRing1EtSumPositiveEta = dir1.make<TH1F>("HfRing1EtSumPositiveEta", "Hf Inner Ring1 Et eta+", 60, 0., 30.);
   theHfRing1EtSumNegativeEta = dir1.make<TH1F>("HfRing1EtSumNegativeEta", "Hf Inner Ring1 Et eta-", 60, 0., 30.);
   theHfRing0CountPositiveEta = dir1.make<TH1F>("HfRing0CountPositiveEta", "Hf Threshold bits Ring0 eta+", 20, 0., 20.);
   theHfRing0CountNegativeEta = dir1.make<TH1F>("HfRing0CountNegativeEta", "Hf Threshold bits Ring0 eta-", 20, 0., 20.);
   theHfRing1CountPositiveEta = dir1.make<TH1F>("HfRing1CountPositiveEta", "Hf Threshold bits Ring1 eta+", 20, 0., 20.);
   theHfRing1CountNegativeEta = dir1.make<TH1F>("HfRing1CountNegativeEta", "Hf Threshold bits Ring1 eta-", 20, 0., 20.);
 }

◆ endJob()

void L1GctValidation::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 288 of file L1GctValidation.cc.

288 {}

Member Data Documentation

◆ m_energy_tag

edm::InputTag L1GctValidation::m_energy_tag

private

Definition at line 56 of file L1GctValidation.h.

Referenced by analyze().

◆ m_gctinp_tag

edm::InputTag L1GctValidation::m_gctinp_tag

private

Definition at line 55 of file L1GctValidation.h.

Referenced by analyze().

◆ m_hfRingEtScaleToken

edm::ESGetToken<L1CaloEtScale, L1HfRingEtScaleRcd> L1GctValidation::m_hfRingEtScaleToken

private

Definition at line 97 of file L1GctValidation.h.

Referenced by analyze().

◆ m_htMissScaleToken

edm::ESGetToken<L1CaloEtScale, L1HtMissScaleRcd> L1GctValidation::m_htMissScaleToken

private

Definition at line 96 of file L1GctValidation.h.

Referenced by analyze().

◆ m_jfParsToken

edm::ESGetToken<L1GctJetFinderParams, L1GctJetFinderParamsRcd> L1GctValidation::m_jfParsToken

private

Definition at line 95 of file L1GctValidation.h.

Referenced by analyze().

◆ theDPhiVsMissEt

TH2F* L1GctValidation::theDPhiVsMissEt

private

Definition at line 78 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theDPhiVsMissHt

TH2F* L1GctValidation::theDPhiVsMissHt

private

Definition at line 79 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHfRing0CountNegativeEta

TH1F* L1GctValidation::theHfRing0CountNegativeEta

private

Definition at line 92 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHfRing0CountPositiveEta

TH1F* L1GctValidation::theHfRing0CountPositiveEta

private

Definition at line 91 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHfRing0EtSumNegativeEta

TH1F* L1GctValidation::theHfRing0EtSumNegativeEta

private

Definition at line 88 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHfRing0EtSumPositiveEta

TH1F* L1GctValidation::theHfRing0EtSumPositiveEta

private

Definition at line 87 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHfRing1CountNegativeEta

TH1F* L1GctValidation::theHfRing1CountNegativeEta

private

Definition at line 94 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHfRing1CountPositiveEta

TH1F* L1GctValidation::theHfRing1CountPositiveEta

private

Definition at line 93 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHfRing1EtSumNegativeEta

TH1F* L1GctValidation::theHfRing1EtSumNegativeEta

private

Definition at line 90 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHfRing1EtSumPositiveEta

TH1F* L1GctValidation::theHfRing1EtSumPositiveEta

private

Definition at line 89 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theHtVsInternalJetsSum

TH2F* L1GctValidation::theHtVsInternalJetsSum

private

Definition at line 81 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissEtAngle

TH1F* L1GctValidation::theMissEtAngle

private

Definition at line 65 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissEtAngleVsInputRegions

TH2F* L1GctValidation::theMissEtAngleVsInputRegions

private

Definition at line 73 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissEtInGeV

TH1F* L1GctValidation::theMissEtInGeV

private

Definition at line 64 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissEtInLsb

TH1F* L1GctValidation::theMissEtInLsb

private

Definition at line 60 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissEtMagVsInputRegions

TH2F* L1GctValidation::theMissEtMagVsInputRegions

private

Definition at line 72 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissEtVector

TH2F* L1GctValidation::theMissEtVector

private

Definition at line 66 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissEtVsMissHt

TH2F* L1GctValidation::theMissEtVsMissHt

private

Definition at line 76 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissEtVsMissHtAngle

TH2F* L1GctValidation::theMissEtVsMissHtAngle

private

Definition at line 77 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHtAngle

TH1F* L1GctValidation::theMissHtAngle

private

Definition at line 68 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHtInGeV

TH1F* L1GctValidation::theMissHtInGeV

private

Definition at line 67 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHtInLsb

TH1F* L1GctValidation::theMissHtInLsb

private

Definition at line 61 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHtMagVsInputRegions

TH2F* L1GctValidation::theMissHtMagVsInputRegions

private

Definition at line 74 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHtPhiVsInternalJetsSum

TH2F* L1GctValidation::theMissHtPhiVsInternalJetsSum

private

Definition at line 83 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHtVector

TH2F* L1GctValidation::theMissHtVector

private

Definition at line 69 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHtVsInternalJetsSum

TH2F* L1GctValidation::theMissHtVsInternalJetsSum

private

Definition at line 82 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHxVsInternalJetsSum

TH2F* L1GctValidation::theMissHxVsInternalJetsSum

private

Definition at line 84 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theMissHyVsInternalJetsSum

TH2F* L1GctValidation::theMissHyVsInternalJetsSum

private

Definition at line 85 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theSumEtInGeV

TH1F* L1GctValidation::theSumEtInGeV

private

Definition at line 62 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theSumEtInLsb

TH1F* L1GctValidation::theSumEtInLsb

private

Definition at line 58 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theSumEtVsInputRegions

TH2F* L1GctValidation::theSumEtVsInputRegions

private

Definition at line 71 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theSumHtInGeV

TH1F* L1GctValidation::theSumHtInGeV

private

Definition at line 63 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

◆ theSumHtInLsb

TH1F* L1GctValidation::theSumHtInLsb

private

Definition at line 59 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ L1GctValidation()

◆ ~L1GctValidation()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

Member Data Documentation

◆ m_energy_tag

◆ m_gctinp_tag

◆ m_hfRingEtScaleToken

◆ m_htMissScaleToken

◆ m_jfParsToken

◆ theDPhiVsMissEt

◆ theDPhiVsMissHt

◆ theHfRing0CountNegativeEta

◆ theHfRing0CountPositiveEta

◆ theHfRing0EtSumNegativeEta

◆ theHfRing0EtSumPositiveEta

◆ theHfRing1CountNegativeEta

◆ theHfRing1CountPositiveEta

◆ theHfRing1EtSumNegativeEta

◆ theHfRing1EtSumPositiveEta

◆ theHtVsInternalJetsSum

◆ theMissEtAngle

◆ theMissEtAngleVsInputRegions

◆ theMissEtInGeV

◆ theMissEtInLsb

◆ theMissEtMagVsInputRegions

◆ theMissEtVector

◆ theMissEtVsMissHt

◆ theMissEtVsMissHtAngle

◆ theMissHtAngle

◆ theMissHtInGeV

◆ theMissHtInLsb

◆ theMissHtMagVsInputRegions

◆ theMissHtPhiVsInternalJetsSum

◆ theMissHtVector

◆ theMissHtVsInternalJetsSum

◆ theMissHxVsInternalJetsSum

◆ theMissHyVsInternalJetsSum

◆ theSumEtInGeV

◆ theSumEtInLsb

◆ theSumEtVsInputRegions

◆ theSumHtInGeV

◆ theSumHtInLsb