CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
L1GctValidation Class Reference

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

Inheritance diagram for L1GctValidation:
edm::one::EDAnalyzer< edm::one::SharedResources > edm::one::EDAnalyzerBase edm::EDConsumerBase

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
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
const EDAnalyzeroperator= (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< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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, L1HfRingEtScaleRcdm_hfRingEtScaleToken
 
edm::ESGetToken< L1CaloEtScale, L1HtMissScaleRcdm_htMissScaleToken
 
edm::ESGetToken< L1GctJetFinderParams, L1GctJetFinderParamsRcdm_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< Bconsumes (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.

17  : m_gctinp_tag(iConfig.getUntrackedParameter<edm::InputTag>("rctInputTag", edm::InputTag("rctDigis"))),
18  m_energy_tag(iConfig.getUntrackedParameter<edm::InputTag>("gctInputTag", edm::InputTag("gctDigis"))),
19  m_jfParsToken(esConsumes<L1GctJetFinderParams, L1GctJetFinderParamsRcd>()),
20  m_htMissScaleToken(esConsumes<L1CaloEtScale, L1HtMissScaleRcd>()),
21  m_hfRingEtScaleToken(esConsumes<L1CaloEtScale, L1HfRingEtScaleRcd>()) {
22  usesResource(TFileService::kSharedResource);
23 }
static const std::string kSharedResource
Definition: TFileService.h:76
edm::ESGetToken< L1CaloEtScale, L1HtMissScaleRcd > m_htMissScaleToken
T getUntrackedParameter(std::string const &, T const &) const
edm::ESGetToken< L1GctJetFinderParams, L1GctJetFinderParamsRcd > m_jfParsToken
edm::InputTag m_energy_tag
edm::InputTag m_gctinp_tag
edm::ESGetToken< L1CaloEtScale, L1HfRingEtScaleRcd > m_hfRingEtScaleToken

◆ ~L1GctValidation()

L1GctValidation::~L1GctValidation ( )
override

Definition at line 25 of file L1GctValidation.cc.

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

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

35  {
36  using namespace edm;
37 
38  // Get the scales from the event setup
40 
41  double lsbForEt = jfPars.product()->getRgnEtLsbGeV();
42  double lsbForHt = jfPars.product()->getHtLsbGeV();
43 
44  unsigned httJetThreshold = static_cast<int>(jfPars.product()->getHtJetEtThresholdGeV() / lsbForHt);
45  unsigned htmJetThreshold = static_cast<int>(jfPars.product()->getMHtJetEtThresholdGeV() / lsbForHt);
46 
49 
50  // Get the Gct energy sums from the event
52  iEvent.getByLabel(m_energy_tag, sumEtColl);
54  iEvent.getByLabel(m_energy_tag, sumHtColl);
56  iEvent.getByLabel(m_energy_tag, missEtColl);
58  iEvent.getByLabel(m_energy_tag, missHtColl);
59 
60  // Get the input calo regions from the event (for checking MEt)
62  iEvent.getByLabel(m_gctinp_tag, inputColl);
63 
64  // Get the internal jet data from the event (for checking Ht)
65  Handle<L1GctInternJetDataCollection> internalJetsColl;
66  iEvent.getByLabel(m_energy_tag, internalJetsColl);
67 
68  double etTot = 0.0;
69  for (L1GctEtTotalCollection::const_iterator jbx = sumEtColl->begin(); jbx != sumEtColl->end(); jbx++) {
70  if (jbx->bx() == 0) {
71  etTot = static_cast<double>(jbx->et());
72  }
73  }
74 
75  double etHad = 0.0;
76  for (L1GctEtHadCollection::const_iterator jbx = sumHtColl->begin(); jbx != sumHtColl->end(); jbx++) {
77  if (jbx->bx() == 0) {
78  etHad = static_cast<double>(jbx->et());
79  }
80  }
81 
82  double etMiss = 0.0;
83  double etMAng = 0.0;
84  for (L1GctEtMissCollection::const_iterator jbx = missEtColl->begin(); jbx != missEtColl->end(); jbx++) {
85  if (jbx->bx() == 0) {
86  etMiss = static_cast<double>(jbx->et());
87  int phibin = jbx->phi();
88  if (phibin >= 36)
89  phibin -= 72;
90  double etMPhi = static_cast<double>(phibin);
91 
92  etMAng = (etMPhi + 0.5) * M_PI / 36.;
93  }
94  }
95 
96  double etTotFromRegions = 0.0;
97  double exTotFromRegions = 0.0;
98  double eyTotFromRegions = 0.0;
99  for (L1CaloRegionCollection::const_iterator jrg = inputColl->begin(); jrg != inputColl->end(); jrg++) {
100  if (jrg->bx() == 0) {
101  double rgEt = static_cast<double>(jrg->et()) * lsbForEt;
102  double rgPhibin = static_cast<double>(jrg->id().iphi());
103  double rgPh = (rgPhibin + 0.5) * M_PI / 9.;
104 
105  etTotFromRegions += rgEt;
106  exTotFromRegions += rgEt * cos(rgPh);
107  eyTotFromRegions += rgEt * sin(rgPh);
108  }
109  }
110 
111  double htMissGct = 0.0;
112  double htMissAng = 0.0;
113  double htMissGeV = 0.0;
114  for (L1GctHtMissCollection::const_iterator jbx = missHtColl->begin(); jbx != missHtColl->end(); jbx++) {
115  if (jbx->bx() == 0) {
116  htMissGct = static_cast<double>(jbx->et());
117  htMissGeV = htMissScale->et(jbx->et());
118  int phibin = jbx->phi();
119  if (phibin >= 9)
120  phibin -= 18;
121  double htMPhi = static_cast<double>(phibin);
122  htMissAng = (htMPhi + 0.5) * M_PI / 9.;
123  }
124  }
125 
126  double htFromJets = 0.0;
127  double hxFromJets = 0.0;
128  double hyFromJets = 0.0;
129  for (L1GctInternJetDataCollection::const_iterator jet = internalJetsColl->begin(); jet != internalJetsColl->end();
130  jet++) {
131  if (jet->bx() == 0 && !jet->empty()) {
132  unsigned jetEtGct = jet->et();
133  double jetEt = static_cast<double>(jetEtGct);
134  int phibin = jet->regionId().iphi();
135  if (phibin >= 9)
136  phibin -= 18;
137  // The phi bin centres are at 0, 20, 40, ... degrees
138  double jetAng = (static_cast<double>(phibin)) * M_PI / 9.;
139  if (jetEtGct > httJetThreshold) {
140  htFromJets += jetEt;
141  }
142  if (jetEtGct > htmJetThreshold) {
143  hxFromJets += jetEt * cos(jetAng);
144  hyFromJets += jetEt * sin(jetAng);
145  }
146  }
147  }
148 
149  double dPhiMetMht = deltaPhi(etMAng, htMissAng);
150 
151  theSumEtInLsb->Fill(etTot);
152  theSumHtInLsb->Fill(etHad);
153  theMissEtInLsb->Fill(etMiss);
154  theMissHtInLsb->Fill(htMissGct);
155  theSumEtInGeV->Fill(etTot * lsbForEt);
156  theSumHtInGeV->Fill(etHad * lsbForHt);
157  theMissEtInGeV->Fill(etMiss * lsbForEt);
158  theMissEtAngle->Fill(etMAng);
159  theMissEtVector->Fill(etMiss * lsbForEt * cos(etMAng), etMiss * lsbForEt * sin(etMAng));
160  if (htMissGct < 126.5) {
161  theMissHtInGeV->Fill(htMissGeV);
162  theMissHtAngle->Fill(htMissAng);
163  theMissHtVector->Fill(htMissGeV * cos(htMissAng), htMissGeV * sin(htMissAng));
164  }
165 
166  theSumEtVsInputRegions->Fill(etTot * lsbForEt, etTotFromRegions);
167  theMissEtMagVsInputRegions->Fill(etMiss * lsbForEt,
168  sqrt(exTotFromRegions * exTotFromRegions + eyTotFromRegions * eyTotFromRegions));
169  theMissEtAngleVsInputRegions->Fill(etMAng, atan2(-eyTotFromRegions, -exTotFromRegions));
170  theMissHtMagVsInputRegions->Fill(htMissGeV,
171  sqrt(exTotFromRegions * exTotFromRegions + eyTotFromRegions * eyTotFromRegions));
172 
173  theMissEtVsMissHt->Fill(etMiss * lsbForEt, htMissGeV);
174  theMissEtVsMissHtAngle->Fill(etMAng, htMissAng);
175  theDPhiVsMissEt->Fill(dPhiMetMht, etMiss * lsbForEt);
176  theDPhiVsMissHt->Fill(dPhiMetMht, htMissGeV);
177 
178  theHtVsInternalJetsSum->Fill(etHad * lsbForHt, htFromJets * lsbForHt);
179  if (htMissGct < 126.5) {
180  theMissHtVsInternalJetsSum->Fill(htMissGeV, sqrt(hxFromJets * hxFromJets + hyFromJets * hyFromJets) * lsbForHt);
181  theMissHtPhiVsInternalJetsSum->Fill(htMissAng, atan2(-hyFromJets, -hxFromJets));
182  theMissHxVsInternalJetsSum->Fill(htMissGeV * cos(htMissAng), hxFromJets * lsbForHt);
183  theMissHyVsInternalJetsSum->Fill(htMissGeV * sin(htMissAng), hyFromJets * lsbForHt);
184  }
185 
186  // Get minbias trigger quantities from HF
189  iEvent.getByLabel(m_energy_tag, HFEtSumsColl);
190  iEvent.getByLabel(m_energy_tag, HFCountsColl);
191 
192  for (L1GctHFRingEtSumsCollection::const_iterator es = HFEtSumsColl->begin(); es != HFEtSumsColl->end(); es++) {
193  if (es->bx() == 0) {
194  theHfRing0EtSumPositiveEta->Fill(hfRingEtScale->et(es->etSum(0)));
195  theHfRing0EtSumNegativeEta->Fill(hfRingEtScale->et(es->etSum(1)));
196  theHfRing1EtSumPositiveEta->Fill(hfRingEtScale->et(es->etSum(2)));
197  theHfRing1EtSumNegativeEta->Fill(hfRingEtScale->et(es->etSum(3)));
198  }
199  }
200 
201  for (L1GctHFBitCountsCollection::const_iterator bc = HFCountsColl->begin(); bc != HFCountsColl->end(); bc++) {
202  if (bc->bx() == 0) {
203  theHfRing0CountPositiveEta->Fill(bc->bitCount(0));
204  theHfRing0CountNegativeEta->Fill(bc->bitCount(1));
205  theHfRing1CountPositiveEta->Fill(bc->bitCount(2));
206  theHfRing1CountNegativeEta->Fill(bc->bitCount(3));
207  }
208  }
209 }
TH1F * theHfRing1EtSumPositiveEta
TH2F * theMissHxVsInternalJetsSum
TH2F * theMissHtVsInternalJetsSum
double getRgnEtLsbGeV() const
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
TH2F * theHtVsInternalJetsSum
double getHtJetEtThresholdGeV() const
edm::ESGetToken< L1CaloEtScale, L1HtMissScaleRcd > m_htMissScaleToken
TH2F * theMissHyVsInternalJetsSum
TH2F * theMissEtVsMissHt
TH1F * theHfRing1CountPositiveEta
TH1F * theHfRing0EtSumPositiveEta
int iEvent
Definition: GenABIO.cc:224
T const * product() const
Definition: ESHandle.h:86
TH1F * theHfRing0CountNegativeEta
double et(const uint16_t rank) const
convert from rank to physically meaningful quantity
T sqrt(T t)
Definition: SSEVec.h:19
edm::ESGetToken< L1GctJetFinderParams, L1GctJetFinderParamsRcd > m_jfParsToken
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
TH2F * theSumEtVsInputRegions
edm::InputTag m_energy_tag
TH2F * theMissEtAngleVsInputRegions
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:130
TH1F * theHfRing1EtSumNegativeEta
edm::InputTag m_gctinp_tag
#define M_PI
TH1F * theHfRing1CountNegativeEta
HLT enums.
edm::ESGetToken< L1CaloEtScale, L1HfRingEtScaleRcd > m_hfRingEtScaleToken
TH1F * theHfRing0CountPositiveEta
TH2F * theMissHtPhiVsInternalJetsSum
TH1F * theHfRing0EtSumNegativeEta
TH2F * theMissEtVsMissHtAngle
TH2F * theMissEtMagVsInputRegions
double getMHtJetEtThresholdGeV() const
TH2F * theMissHtMagVsInputRegions

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

212  {
214 
215  TFileDirectory dir0 = fs->mkdir("L1GctEtSums");
216 
217  theSumEtInLsb = dir0.make<TH1F>("SumEtInLsb", "Total Et (GCT units)", 128, 0., 2048.);
218  theSumHtInLsb = dir0.make<TH1F>("SumHtInLsb", "Total Ht (GCT units)", 128, 0., 2048.);
219  theMissEtInLsb = dir0.make<TH1F>("MissEtInLsb", "Missing Et magnitude (GCT units)", 128, 0., 1024.);
220  theMissHtInLsb = dir0.make<TH1F>("MissHtInLsb", "Missing Ht magnitude (GCT units)", 128, 0., 127.);
221  theSumEtInGeV = dir0.make<TH1F>("SumEtInGeV", "Total Et (in GeV)", 100, 0., 1000.);
222  theSumHtInGeV = dir0.make<TH1F>("SumHtInGeV", "Total Ht (in GeV)", 100, 0., 1000.);
223  theMissEtInGeV = dir0.make<TH1F>("MissEtInGeV", "Missing Et magnitude (in GeV)", 100, 0., 500.);
224  theMissEtAngle = dir0.make<TH1F>("MissEtAngle", "Missing Et angle", 72, -M_PI, M_PI);
225  theMissEtVector = dir0.make<TH2F>("MissEtVector", "Missing Ex vs Missing Ey", 100, -100., 100., 100, -100., 100.);
226  theMissHtInGeV = dir0.make<TH1F>("MissHtInGeV", "Missing Ht magnitude (in GeV)", 100, 0., 500.);
227  theMissHtAngle = dir0.make<TH1F>("MissHtAngle", "Missing Ht angle", 72, -M_PI, M_PI);
228  theMissHtVector = dir0.make<TH2F>("MissHtVector", "Missing Hx vs Missing Hy", 100, -100., 100., 100, -100., 100.);
230  dir0.make<TH2F>("SumEtVsInputRegions", "Total Et vs sum of input regions", 100, 0., 1000., 100, 0., 1000.);
231  theMissEtMagVsInputRegions = dir0.make<TH2F>(
232  "MissEtMagVsInputRegions", "Missing Et magnitude vs sum of input regions", 100, 0., 500., 100, 0., 500.);
233  theMissEtAngleVsInputRegions = dir0.make<TH2F>(
234  "MissEtAngleVsInputRegions", "Missing Et angle vs sum of input regions", 72, -M_PI, M_PI, 72, -M_PI, M_PI);
235  theMissHtMagVsInputRegions = dir0.make<TH2F>(
236  "MissHtMagVsInputRegions", "Missing Ht magnitude vs sum of input regions", 100, 0., 500., 100, 0., 500.);
237  theMissEtVsMissHt = dir0.make<TH2F>("MissEtVsMissHt", "Missing Et vs Missing Ht", 100, 0., 500., 100, 0., 500.);
238  theMissEtVsMissHtAngle = dir0.make<TH2F>(
239  "MissEtVsMissHtAngle", "Angle correlation Missing Et vs Missing Ht", 72, -M_PI, M_PI, 72, -M_PI, M_PI);
241  dir0.make<TH2F>("theDPhiVsMissEt", "Angle difference MET-MHT vs MET magnitude", 72, -M_PI, M_PI, 100, 0., 500.);
243  dir0.make<TH2F>("theDPhiVsMissHt", "Angle difference MET-MHT vs MHT magnitude", 72, -M_PI, M_PI, 100, 0., 500.);
244 
245  theHtVsInternalJetsSum = dir0.make<TH2F>(
246  "HtVsInternalJetsSum", "Ht vs scalar sum of jet Et values (in GeV)", 128, 0., 2048., 128, 0., 2048.);
247  theMissHtVsInternalJetsSum = dir0.make<TH2F>(
248  "MissHtVsInternalJetsSum", "Missing Ht vs vector sum of jet Et values (in GeV)", 128, 0., 512., 128, 0., 512.);
249  theMissHtPhiVsInternalJetsSum = dir0.make<TH2F>("MissHtPhiVsInternalJetsSum",
250  "Angle correlation Missing Ht vs vector sum of jet Et values",
251  72,
252  -M_PI,
253  M_PI,
254  72,
255  -M_PI,
256  M_PI);
257  theMissHxVsInternalJetsSum = dir0.make<TH2F>("MissHxVsInternalJetsSum",
258  "Missing Ht x component vs sum of jet Et values (in GeV)",
259  128,
260  -256.,
261  256.,
262  128,
263  -256.,
264  256.);
265  theMissHyVsInternalJetsSum = dir0.make<TH2F>("MissHyVsInternalJetsSum",
266  "Missing Ht y component vs sum of jet Et values (in GeV)",
267  128,
268  -256.,
269  256.,
270  128,
271  -256.,
272  256.);
273 
274  TFileDirectory dir1 = fs->mkdir("L1GctHfSumsAndJetCounts");
275 
276  // Minimum bias triggers from Hf inner rings
277  theHfRing0EtSumPositiveEta = dir1.make<TH1F>("HfRing0EtSumPositiveEta", "Hf Inner Ring0 Et eta+", 60, 0., 30.);
278  theHfRing0EtSumNegativeEta = dir1.make<TH1F>("HfRing0EtSumNegativeEta", "Hf Inner Ring0 Et eta-", 60, 0., 30.);
279  theHfRing1EtSumPositiveEta = dir1.make<TH1F>("HfRing1EtSumPositiveEta", "Hf Inner Ring1 Et eta+", 60, 0., 30.);
280  theHfRing1EtSumNegativeEta = dir1.make<TH1F>("HfRing1EtSumNegativeEta", "Hf Inner Ring1 Et eta-", 60, 0., 30.);
281  theHfRing0CountPositiveEta = dir1.make<TH1F>("HfRing0CountPositiveEta", "Hf Threshold bits Ring0 eta+", 20, 0., 20.);
282  theHfRing0CountNegativeEta = dir1.make<TH1F>("HfRing0CountNegativeEta", "Hf Threshold bits Ring0 eta-", 20, 0., 20.);
283  theHfRing1CountPositiveEta = dir1.make<TH1F>("HfRing1CountPositiveEta", "Hf Threshold bits Ring1 eta+", 20, 0., 20.);
284  theHfRing1CountNegativeEta = dir1.make<TH1F>("HfRing1CountNegativeEta", "Hf Threshold bits Ring1 eta-", 20, 0., 20.);
285 }
TH1F * theHfRing1EtSumPositiveEta
TH2F * theMissHxVsInternalJetsSum
TH2F * theMissHtVsInternalJetsSum
TH2F * theHtVsInternalJetsSum
TH2F * theMissHyVsInternalJetsSum
TH2F * theMissEtVsMissHt
TH1F * theHfRing1CountPositiveEta
T * make(const Args &...args) const
make new ROOT object
TH1F * theHfRing0EtSumPositiveEta
TH1F * theHfRing0CountNegativeEta
TH2F * theSumEtVsInputRegions
TH2F * theMissEtAngleVsInputRegions
TH1F * theHfRing1EtSumNegativeEta
#define M_PI
TH1F * theHfRing1CountNegativeEta
TH1F * theHfRing0CountPositiveEta
TH2F * theMissHtPhiVsInternalJetsSum
TH1F * theHfRing0EtSumNegativeEta
TH2F * theMissEtVsMissHtAngle
TH2F * theMissEtMagVsInputRegions
TH2F * theMissHtMagVsInputRegions

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