L1GctValidation Class Reference

#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::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
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::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	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
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)
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 () noexcept
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...
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 ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 16 of file L1GctValidation.cc.

    : 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>()) {}

L1GctValidation::~L1GctValidation ( )

override

Definition at line 23 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

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 33 of file L1GctValidation.cc.

References funct::cos(), srCondWrite_cfg::deltaPhi, edm::Event::getByLabel(), edm::EventSetup::getHandle(), metsig::jet, m_energy_tag, m_gctinp_tag, m_hfRingEtScaleToken, m_htMissScaleToken, m_jfParsToken, M_PI, 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));
    }
  }
}

void L1GctValidation::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 210 of file L1GctValidation.cc.

References M_PI, TFileDirectory::make(), TFileService::mkdir(), 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.);
}

void L1GctValidation::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 286 of file L1GctValidation.cc.

{}

Member Data Documentation

edm::InputTag L1GctValidation::m_energy_tag

private

Definition at line 56 of file L1GctValidation.h.

Referenced by analyze().

edm::InputTag L1GctValidation::m_gctinp_tag

private

Definition at line 55 of file L1GctValidation.h.

Referenced by analyze().

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

private

Definition at line 97 of file L1GctValidation.h.

Referenced by analyze().

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

private

Definition at line 96 of file L1GctValidation.h.

Referenced by analyze().

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

private

Definition at line 95 of file L1GctValidation.h.

Referenced by analyze().

TH2F* L1GctValidation::theDPhiVsMissEt

private

Definition at line 78 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theDPhiVsMissHt

private

Definition at line 79 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theHfRing0CountNegativeEta

private

Definition at line 92 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theHfRing0CountPositiveEta

private

Definition at line 91 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theHfRing0EtSumNegativeEta

private

Definition at line 88 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theHfRing0EtSumPositiveEta

private

Definition at line 87 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theHfRing1CountNegativeEta

private

Definition at line 94 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theHfRing1CountPositiveEta

private

Definition at line 93 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theHfRing1EtSumNegativeEta

private

Definition at line 90 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theHfRing1EtSumPositiveEta

private

Definition at line 89 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theHtVsInternalJetsSum

private

Definition at line 81 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theMissEtAngle

private

Definition at line 65 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissEtAngleVsInputRegions

private

Definition at line 73 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theMissEtInGeV

private

Definition at line 64 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theMissEtInLsb

private

Definition at line 60 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissEtMagVsInputRegions

private

Definition at line 72 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissEtVector

private

Definition at line 66 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissEtVsMissHt

private

Definition at line 76 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissEtVsMissHtAngle

private

Definition at line 77 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theMissHtAngle

private

Definition at line 68 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theMissHtInGeV

private

Definition at line 67 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theMissHtInLsb

private

Definition at line 61 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissHtMagVsInputRegions

private

Definition at line 74 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissHtPhiVsInternalJetsSum

private

Definition at line 83 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissHtVector

private

Definition at line 69 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissHtVsInternalJetsSum

private

Definition at line 82 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissHxVsInternalJetsSum

private

Definition at line 84 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theMissHyVsInternalJetsSum

private

Definition at line 85 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theSumEtInGeV

private

Definition at line 62 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theSumEtInLsb

private

Definition at line 58 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH2F* L1GctValidation::theSumEtVsInputRegions

private

Definition at line 71 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theSumHtInGeV

private

Definition at line 63 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().

TH1F* L1GctValidation::theSumHtInLsb

private

Definition at line 59 of file L1GctValidation.h.

Referenced by analyze(), and beginJob().