CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
L1RCTSaveInput Class Reference

#include <L1RCTSaveInput.h>

Inheritance diagram for L1RCTSaveInput:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
 L1RCTSaveInput (const edm::ParameterSet &)
 
 ~L1RCTSaveInput () 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
 
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::vector< ModuleDescription const * > &modules, 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
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Attributes

edm::InputTag ecalDigisLabel
 
std::string fileName
 
edm::InputTag hcalDigisLabel
 
std::ofstream ofs
 
L1RCT * rct
 
L1RCTLookupTables * rctLookupTables
 
bool useDebugTpgScales
 
bool useEcal
 
bool useHcal
 

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 wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
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<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)
 

Detailed Description

src/L1RCTSaveInput/src/L1RCTSaveInput.cc

Description: Saves the input event from TPGs for loading simulated events in hardware

Implementation: Kind of kludgy – should think of a better way in future

Definition at line 39 of file L1RCTSaveInput.h.

Constructor & Destructor Documentation

L1RCTSaveInput::L1RCTSaveInput ( const edm::ParameterSet &  conf)
explicit

Definition at line 45 of file L1RCTSaveInput.cc.

References MessageLogger_cfi::cerr, cmsRelvalreport::exit, fileName, and ofs.

46  : fileName(conf.getUntrackedParameter<std::string>("rctTestInputFile")),
47  rctLookupTables(new L1RCTLookupTables), rct(new L1RCT(rctLookupTables)),
48  useEcal(conf.getParameter<bool>("useEcal")),
49  useHcal(conf.getParameter<bool>("useHcal")),
50  ecalDigisLabel(conf.getParameter<edm::InputTag>("ecalDigisLabel")),
51  hcalDigisLabel(conf.getParameter<edm::InputTag>("hcalDigisLabel")),
52  useDebugTpgScales(conf.getParameter<bool>("useDebugTpgScales")) {
53  ofs.open(fileName.c_str(), std::ios::app);
54  if (!ofs) {
55  std::cerr << "Could not create " << fileName << std::endl;
56  exit(1);
57  }
58 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
L1RCTSaveInput::ecalDigisLabel
edm::InputTag ecalDigisLabel
Definition: L1RCTSaveInput.h:51
MessageLogger_cfi.cerr
cerr
Definition: MessageLogger_cfi.py:547
L1RCTSaveInput::rctLookupTables
L1RCTLookupTables * rctLookupTables
Definition: L1RCTSaveInput.h:47
L1RCTSaveInput::fileName
std::string fileName
Definition: L1RCTSaveInput.h:46
L1RCTSaveInput::ofs
std::ofstream ofs
Definition: L1RCTSaveInput.h:54
L1RCT
Definition: L1RCT.h:20
L1RCTSaveInput::hcalDigisLabel
edm::InputTag hcalDigisLabel
Definition: L1RCTSaveInput.h:52
L1RCTSaveInput::~L1RCTSaveInput ( )
override

Definition at line 60 of file L1RCTSaveInput.cc.

References rct, and rctLookupTables.

60  {
61  if (rct != nullptr)
62  delete rct;
63  if (rctLookupTables != nullptr)
64  delete rctLookupTables;
65 }
L1RCTSaveInput::rctLookupTables
L1RCTLookupTables * rctLookupTables
Definition: L1RCTSaveInput.h:47

Member Function Documentation

void L1RCTSaveInput::analyze ( const edm::Event &  event,
const edm::EventSetup &  eventSetup 
)
override

Definition at line 67 of file L1RCTSaveInput.cc.

References EnergyCorrector::c, gather_cfg::cout, TauDecayModes::dec, L1RCT::digiInput(), MillePedeFileConverter_cfg::e, digitizers_cfi::ecal, EcalBarrel, L1RCT::ecalCompressedET(), ecalDigisLabel, EcalEndcap, L1RCT::ecalFineGrainBit(), edm::EventSetup::get(), EcalTPGScale::getTPGInGeV(), h, digitizers_cfi::hcal, L1RCT::hcalCompressedET(), hcalDigisLabel, CaloTPGTranscoder::hcaletValue(), L1RCT::hcalFineGrainBit(), edm::HandleBase::isValid(), L1RCTLookupTables::lookup(), L1CaloEcalScale::nBinEta, L1CaloHcalScale::nBinEta, L1CaloEcalScale::nBinRank, L1CaloHcalScale::nBinRank, nEvents, ofs, edm::ESHandle< T >::product(), alignCSCRings::r, rct, rctLookupTables, alignCSCRings::s, L1CaloHcalScale::setBin(), L1CaloEcalScale::setBin(), L1RCTLookupTables::setChannelMask(), L1RCTLookupTables::setEcalScale(), EcalTPGScale::setEventSetup(), L1RCTLookupTables::setHcalScale(), L1RCTLookupTables::setL1CaloEtScale(), L1RCTLookupTables::setRCTParameters(), and useDebugTpgScales.

68  {
69  edm::ESHandle<L1RCTParameters> rctParameters;
70  eventSetup.get<L1RCTParametersRcd>().get(rctParameters);
71  const L1RCTParameters *r = rctParameters.product();
72  edm::ESHandle<L1RCTChannelMask> channelMask;
73  eventSetup.get<L1RCTChannelMaskRcd>().get(channelMask);
74  const L1RCTChannelMask *c = channelMask.product();
75  edm::ESHandle<L1CaloEtScale> emScale;
76  eventSetup.get<L1EmEtScaleRcd>().get(emScale);
77  const L1CaloEtScale *s = emScale.product();
78 
79  rctLookupTables->setRCTParameters(r);
80  rctLookupTables->setChannelMask(c);
81  rctLookupTables->setL1CaloEtScale(s);
82 
83  if (useDebugTpgScales) {
84  // use old-style scales
85  edm::ESHandle<CaloTPGTranscoder> transcoder;
86  eventSetup.get<CaloTPGRecord>().get(transcoder);
87  const CaloTPGTranscoder *h_tpg = transcoder.product();
88 
89  EcalTPGScale *e_tpg = new EcalTPGScale();
90  e_tpg->setEventSetup(eventSetup);
91 
92  L1CaloEcalScale *ecalScale = new L1CaloEcalScale();
93  L1CaloHcalScale *hcalScale = new L1CaloHcalScale();
94 
95  // create input scales, werner's code
96  // ECAL
97  std::cout << "ECAL Pos " << L1CaloEcalScale::nBinRank << std::endl;
98  for (unsigned short ieta = 1; ieta <= L1CaloEcalScale::nBinEta; ++ieta) {
99  for (unsigned short irank = 0; irank < L1CaloEcalScale::nBinRank;
100  ++irank) {
101  std::cout << ieta << " " << irank;
102  EcalSubdetector subdet = (ieta <= 17) ? EcalBarrel : EcalEndcap;
103  double etGeVPos =
104  e_tpg->getTPGInGeV(irank, EcalTrigTowerDetId(1, // +ve eta
105  subdet, ieta,
106  1)); // dummy phi value
107  ecalScale->setBin(irank, ieta, 1, etGeVPos);
108  std::cout << etGeVPos << ", ";
109  }
110  std::cout << std::endl;
111  }
112  std::cout << std::endl;
113 
114  std::cout << "ECAL Neg" << std::endl;
115  for (unsigned short ieta = 1; ieta <= L1CaloEcalScale::nBinEta; ++ieta) {
116  for (unsigned short irank = 0; irank < L1CaloEcalScale::nBinRank;
117  ++irank) {
118  EcalSubdetector subdet = (ieta <= 17) ? EcalBarrel : EcalEndcap;
119 
120  std::cout << ieta << " " << irank;
121  double etGeVNeg =
122  e_tpg->getTPGInGeV(irank,
123  EcalTrigTowerDetId(-1, // -ve eta
124  subdet, ieta,
125  2)); // dummy phi value
126  ecalScale->setBin(irank, ieta, -1, etGeVNeg);
127  std::cout << etGeVNeg << ", ";
128  }
129  std::cout << std::endl;
130  }
131  std::cout << std::endl;
132 
133  // HCAL
134  std::cout << "HCAL" << std::endl;
135  for (unsigned short ieta = 1; ieta <= L1CaloHcalScale::nBinEta; ++ieta) {
136  for (unsigned short irank = 0; irank < L1CaloHcalScale::nBinRank;
137  ++irank) {
138  double etGeV = h_tpg->hcaletValue(ieta, irank);
139 
140  hcalScale->setBin(irank, ieta, 1, etGeV);
141  hcalScale->setBin(irank, ieta, -1, etGeV);
142  std::cout << etGeV << ", ";
143  std::cout << std::endl;
144  }
145  std::cout << std::endl;
146  }
147 
148  // set the input scales
149  rctLookupTables->setEcalScale(ecalScale);
150  rctLookupTables->setHcalScale(hcalScale);
151 
152  delete e_tpg;
153 
154  } else {
155  edm::ESHandle<L1CaloHcalScale> hcalScale;
156  eventSetup.get<L1CaloHcalScaleRcd>().get(hcalScale);
157  const L1CaloHcalScale *h = hcalScale.product();
158  edm::ESHandle<L1CaloEcalScale> ecalScale;
159  eventSetup.get<L1CaloEcalScaleRcd>().get(ecalScale);
160  const L1CaloEcalScale *e = ecalScale.product();
161 
162  rctLookupTables->setHcalScale(h);
163  rctLookupTables->setEcalScale(e);
164  }
165 
166  edm::Handle<EcalTrigPrimDigiCollection> ecal;
167  edm::Handle<HcalTrigPrimDigiCollection> hcal;
168  event.getByLabel(ecalDigisLabel, ecal);
169  event.getByLabel(hcalDigisLabel, hcal);
170  EcalTrigPrimDigiCollection ecalColl;
171  HcalTrigPrimDigiCollection hcalColl;
172  if (ecal.isValid()) {
173  ecalColl = *ecal;
174  }
175  if (hcal.isValid()) {
176  hcalColl = *hcal;
177  }
178  rct->digiInput(ecalColl, hcalColl);
179  static int nEvents = 0;
180  if (nEvents == 0) {
181  ofs << "Crate = 0-17" << std::endl
182  << "Card = 0-7 within the crate" << std::endl
183  << "Tower = 0-31 covers 4 x 8 covered by the card" << std::endl
184  << "EMAddr(0:8) = EMFGBit(0:0)+CompressedEMET(1:8)" << std::endl
185  << "HDAddr(0:8) = HDFGBit(0:0)+CompressedHDET(1:8) - note: "
186  "HDFGBit(0:0) is not part of the hardware LUT address"
187  << std::endl
188  << "LutOut(0:17)= "
189  "LinearEMET(0:6)+HoEFGVetoBit(7:7)+LinearJetET(8:16)+ActivityBit(17:"
190  "17)"
191  << std::endl
192  << "Event"
193  << "\t"
194  << "Crate"
195  << "\t"
196  << "Card"
197  << "\t"
198  << "Tower"
199  << "\t"
200  << "EMAddr"
201  << "\t"
202  << "HDAddr"
203  << "\t"
204  << "LUTOut" << std::endl;
205  }
206  if (nEvents < 64) {
207  for (unsigned short iCrate = 0; iCrate < 18; iCrate++) {
208  for (unsigned short iCard = 0; iCard < 7; iCard++) {
209  // tower numbered from 0-31
210  for (unsigned short iTower = 0; iTower < 32; iTower++) {
211  unsigned short ecal = rct->ecalCompressedET(iCrate, iCard, iTower);
212  unsigned short hcal = rct->hcalCompressedET(iCrate, iCard, iTower);
213  unsigned short fgbit = rct->ecalFineGrainBit(iCrate, iCard, iTower);
214  unsigned short mubit = rct->hcalFineGrainBit(iCrate, iCard, iTower);
215  unsigned long lutOutput =
216  rctLookupTables->lookup(ecal, hcal, fgbit, iCrate, iCard, iTower);
217  ofs << std::hex << nEvents << "\t" << iCrate << "\t" << iCard << "\t"
218  << iTower << "\t" << ecal * 2 + fgbit << "\t" << hcal * 2 + mubit
219  << "\t" << lutOutput << std::dec << std::endl;
220  }
221  }
222  }
223  }
224  nEvents++;
225 }
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
EcalTPGScale::setEventSetup
void setEventSetup(const edm::EventSetup &evtSetup)
Definition: EcalTPGScale.cc:19
L1RCT::hcalCompressedET
unsigned short hcalCompressedET(int crate, int card, int tower)
Definition: L1RCT.h:74
L1RCTSaveInput::ecalDigisLabel
edm::InputTag ecalDigisLabel
Definition: L1RCTSaveInput.h:51
L1CaloHcalScale::setBin
void setBin(unsigned short rank, unsigned short eta, short etaSign, double et)
set scale element; use this to create non-linear scales
Definition: L1CaloHcalScale.cc:36
L1RCTSaveInput::rctLookupTables
L1RCTLookupTables * rctLookupTables
Definition: L1RCTSaveInput.h:47
L1CaloEcalScale::setBin
void setBin(unsigned short rank, unsigned short eta, short etaSign, double et)
set scale element; use this to create non-linear scales
Definition: L1CaloEcalScale.cc:36
L1RCT::hcalFineGrainBit
unsigned short hcalFineGrainBit(int crate, int card, int tower)
Definition: L1RCT.h:77
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
EcalTPGScale::getTPGInGeV
double getTPGInGeV(const EcalTriggerPrimitiveDigi &tpDigi)
Definition: EcalTPGScale.cc:24
L1RCTLookupTables::setHcalScale
void setHcalScale(const L1CaloHcalScale *hcalScale)
Definition: L1RCTLookupTables.h:35
L1RCT::ecalFineGrainBit
unsigned short ecalFineGrainBit(int crate, int card, int tower)
Definition: L1RCT.h:71
L1RCT::ecalCompressedET
unsigned short ecalCompressedET(int crate, int card, int tower)
Definition: L1RCT.h:68
digitizers_cfi.hcal
hcal
Definition: digitizers_cfi.py:25
L1CaloHcalScale::nBinRank
static const unsigned short nBinRank
Definition: L1CaloHcalScale.h:33
digitizers_cfi.ecal
ecal
Definition: digitizers_cfi.py:22
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
L1RCTLookupTables::setL1CaloEtScale
void setL1CaloEtScale(const L1CaloEtScale *etScale)
Definition: L1RCTLookupTables.h:39
L1RCTSaveInput::ofs
std::ofstream ofs
Definition: L1RCTSaveInput.h:54
L1RCTLookupTables::setEcalScale
void setEcalScale(const L1CaloEcalScale *ecalScale)
Definition: L1RCTLookupTables.h:41
L1RCTLookupTables::setRCTParameters
void setRCTParameters(const L1RCTParameters *rctParameters)
Definition: L1RCTLookupTables.h:23
L1CaloEcalScale::nBinRank
static const unsigned short nBinRank
Definition: L1CaloEcalScale.h:33
L1RCT::digiInput
void digiInput(const EcalTrigPrimDigiCollection &ecalCollection, const HcalTrigPrimDigiCollection &hcalCollection)
Definition: L1RCT.cc:109
edm::EventSetup::get
T get() const
Definition: EventSetup.h:71
L1CaloHcalScale::nBinEta
static const unsigned short nBinEta
Definition: L1CaloHcalScale.h:34
L1CaloEcalScale::nBinEta
static const unsigned short nBinEta
Definition: L1CaloEcalScale.h:34
CaloTPGTranscoder::hcaletValue
virtual double hcaletValue(const int &ieta, const int &iphi, const int &version, const int &compressedValue) const =0
gather_cfg.cout
cout
Definition: gather_cfg.py:144
nEvents
UInt_t nEvents
Definition: hcalCalib.cc:41
EcalSubdetector
EcalSubdetector
Definition: EcalSubdetector.h:10
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
L1RCTLookupTables::setChannelMask
void setChannelMask(const L1RCTChannelMask *channelMask)
Definition: L1RCTLookupTables.h:27
L1RCTLookupTables::lookup
unsigned int lookup(unsigned short ecalInput, unsigned short hcalInput, unsigned short fgbit, unsigned short crtNo, unsigned short crdNo, unsigned short twrNo) const
Definition: L1RCTLookupTables.cc:23
L1RCTSaveInput::hcalDigisLabel
edm::InputTag hcalDigisLabel
Definition: L1RCTSaveInput.h:52

Member Data Documentation

edm::InputTag L1RCTSaveInput::ecalDigisLabel
private

Definition at line 51 of file L1RCTSaveInput.h.

Referenced by analyze().

std::string L1RCTSaveInput::fileName
private

Definition at line 46 of file L1RCTSaveInput.h.

Referenced by L1RCTSaveInput().

edm::InputTag L1RCTSaveInput::hcalDigisLabel
private

Definition at line 52 of file L1RCTSaveInput.h.

Referenced by analyze().

std::ofstream L1RCTSaveInput::ofs
private

Definition at line 54 of file L1RCTSaveInput.h.

Referenced by analyze(), and L1RCTSaveInput().

L1RCT* L1RCTSaveInput::rct
private

Definition at line 48 of file L1RCTSaveInput.h.

Referenced by analyze(), and ~L1RCTSaveInput().

L1RCTLookupTables* L1RCTSaveInput::rctLookupTables
private

Definition at line 47 of file L1RCTSaveInput.h.

Referenced by analyze(), and ~L1RCTSaveInput().

bool L1RCTSaveInput::useDebugTpgScales
private

Definition at line 53 of file L1RCTSaveInput.h.

Referenced by analyze().

bool L1RCTSaveInput::useEcal
private

Definition at line 49 of file L1RCTSaveInput.h.

bool L1RCTSaveInput::useHcal
private

Definition at line 50 of file L1RCTSaveInput.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11