CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Attributes
L1RCTSaveInput Class Reference

#include <src/L1RCTSaveInput/src/L1RCTSaveInput.cc>

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

Public Member Functions

virtual void analyze (const edm::Event &, const edm::EventSetup &)
 
 L1RCTSaveInput (const edm::ParameterSet &)
 
 ~L1RCTSaveInput ()
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

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
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
- 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 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

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 38 of file L1RCTSaveInput.h.

Constructor & Destructor Documentation

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

Definition at line 45 of file L1RCTSaveInput.cc.

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

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

Definition at line 63 of file L1RCTSaveInput.cc.

References rct, and rctLookupTables.

64 {
65  if(rct != 0) delete rct;
66  if(rctLookupTables != 0) delete rctLookupTables;
67 }
L1RCTSaveInput::rctLookupTables
L1RCTLookupTables * rctLookupTables
Definition: L1RCTSaveInput.h:45

Member Function Documentation

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

Implements edm::EDAnalyzer.

Definition at line 70 of file L1RCTSaveInput.cc.

References EnergyCorrector::c, gather_cfg::cout, TauDecayModes::dec, L1RCT::digiInput(), alignCSCRings::e, patCandidatesForDimuonsSequences_cff::ecal, EcalBarrel, L1RCT::ecalCompressedET(), ecalDigisLabel, EcalEndcap, L1RCT::ecalFineGrainBit(), edm::EventSetup::get(), EcalTPGScale::getTPGInGeV(), h, patCandidatesForDimuonsSequences_cff::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< class >::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.

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

Member Data Documentation

edm::InputTag L1RCTSaveInput::ecalDigisLabel
private

Definition at line 49 of file L1RCTSaveInput.h.

Referenced by analyze().

std::string L1RCTSaveInput::fileName
private

Definition at line 44 of file L1RCTSaveInput.h.

Referenced by L1RCTSaveInput().

edm::InputTag L1RCTSaveInput::hcalDigisLabel
private

Definition at line 50 of file L1RCTSaveInput.h.

Referenced by analyze().

std::ofstream L1RCTSaveInput::ofs
private

Definition at line 52 of file L1RCTSaveInput.h.

Referenced by analyze(), and L1RCTSaveInput().

L1RCT* L1RCTSaveInput::rct
private

Definition at line 46 of file L1RCTSaveInput.h.

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

L1RCTLookupTables* L1RCTSaveInput::rctLookupTables
private

Definition at line 45 of file L1RCTSaveInput.h.

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

bool L1RCTSaveInput::useDebugTpgScales
private

Definition at line 51 of file L1RCTSaveInput.h.

Referenced by analyze().

bool L1RCTSaveInput::useEcal
private

Definition at line 47 of file L1RCTSaveInput.h.

bool L1RCTSaveInput::useHcal
private

Definition at line 48 of file L1RCTSaveInput.h.

Generated for CMSSW Reference Manual by doxygen 1.8.5