CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalTrigPrimDigiProducer.cc
Go to the documentation of this file.
22 
23 #include <algorithm>
24 
25 
26 
27 
29 :
30  theAlgo_(ps.getParameter<bool>("peakFilter"),
31  ps.getParameter<std::vector<double> >("weights"),
32  ps.getParameter<int>("latency"),
33  ps.getParameter<uint32_t>("FG_threshold"),
34  ps.getParameter<uint32_t>("ZS_threshold"),
35  ps.getParameter<int>("numberOfSamples"),
36  ps.getParameter<int>("numberOfPresamples"),
37  ps.getParameter<int>("numberOfSamplesHF"),
38  ps.getParameter<int>("numberOfPresamplesHF"),
39  ps.getParameter<uint32_t>("MinSignalThreshold"),
40  ps.getParameter<uint32_t>("PMTNoiseThreshold")
41  ),
42  inputLabel_(ps.getParameter<std::vector<edm::InputTag> >("inputLabel")),
43  inputTagFEDRaw_(ps.getParameter<edm::InputTag> ("InputTagFEDRaw")),
44  runZS_(ps.getParameter<bool>("RunZS")),
45  runFrontEndFormatError_(ps.getParameter<bool>("FrontEndFormatError"))
46 {
47  HFEMB_ = false;
48  if(ps.exists("LSConfig"))
49  {
51  HFEMB_ = LongShortCut_.getParameter<bool>("HcalFeatureHFEMBit");
52  MinLongEnergy_ = LongShortCut_.getParameter<double>("Min_Long_Energy"); //minimum long energy
53  MinShortEnergy_ = LongShortCut_.getParameter<double>("Min_Short_Energy"); //minimum short energy
54  LongShortSlope_ = LongShortCut_.getParameter<double>("Long_vrs_Short_Slope"); //slope of the line that cuts are based on
55  LongShortOffset_ = LongShortCut_.getParameter<double>("Long_Short_Offset"); //offset of line
56  }
57  // register for data access
59  tok_raw_ = consumes<FEDRawDataCollection>(inputTagFEDRaw_);
60  }
61  tok_hbhe_ = consumes<HBHEDigiCollection>(inputLabel_[0]);
62  tok_hf_ = consumes<HFDigiCollection>(inputLabel_[1]);
63 
64  produces<HcalTrigPrimDigiCollection>();
65  theAlgo_.setPeakFinderAlgorithm(ps.getParameter<int>("PeakFinderAlgorithm"));
66 }
67 
68 
70 
71  // Step A: get the conditions, for the decoding
72  edm::ESHandle<HcalTPGCoder> inputCoder;
73  eventSetup.get<HcalTPGRecord>().get(inputCoder);
74 
76  eventSetup.get<CaloTPGRecord>().get(outTranscoder);
77 
79  eventSetup.get<HcalLutMetadataRcd>().get(lutMetadata);
80  float rctlsb = lutMetadata->getRctLsb();
81 
83  eventSetup.get<CaloGeometryRecord>().get(pG);
84 
85  // Step B: Create empty output
86  std::auto_ptr<HcalTrigPrimDigiCollection> result(new HcalTrigPrimDigiCollection());
87 
90 
91  iEvent.getByToken(tok_hbhe_,hbheDigis);
92  iEvent.getByToken(tok_hf_,hfDigis);
93 
94  // protect here against missing input collections
95  // there is no protection in HcalTriggerPrimitiveAlgo
96 
97  if (!hbheDigis.isValid()) {
98  edm::LogInfo("HcalTrigPrimDigiProducer")
99  << "\nWarning: HBHEDigiCollection with input tag "
100  << inputLabel_[0]
101  << "\nrequested in configuration, but not found in the event."
102  << "\nQuit returning empty product." << std::endl;
103 
104  // put empty HcalTrigPrimDigiCollection in the event
105  iEvent.put(result);
106 
107  return;
108  }
109 
110  if (!hfDigis.isValid()) {
111  edm::LogInfo("HcalTrigPrimDigiProducer")
112  << "\nWarning: HFDigiCollection with input tag "
113  << inputLabel_[1]
114  << "\nrequested in configuration, but not found in the event."
115  << "\nQuit returning empty product." << std::endl;
116 
117  // put empty HcalTrigPrimDigiCollection in the event
118  iEvent.put(result);
119 
120  return;
121  }
122 
123 
125  eventSetup.get<HcalDbRecord> ().get(pSetup);
126 
127  HcalFeatureBit* hfembit = 0;
128 
129  if(HFEMB_)
130  {
131  hfembit = new HcalFeatureHFEMBit(MinShortEnergy_, MinLongEnergy_, LongShortSlope_, LongShortOffset_, *pSetup); //inputs values that cut will be based on
132  theAlgo_.run(inputCoder.product(), outTranscoder->getHcalCompressor().get(),
133  *hbheDigis, *hfDigis, *result, &(*pG), rctlsb, hfembit);
134 
135  }// Step C: Invoke the algorithm, passing in inputs and getting back outputs.
136  else
137  {
138  theAlgo_.run(inputCoder.product(), outTranscoder->getHcalCompressor().get(),
139  *hbheDigis, *hfDigis, *result, &(*pG), rctlsb);
140  }
141 
142  // Step C: Invoke the algorithm, passing in inputs and getting back outputs.
143 
144 
145  // Step C.1: Run FE Format Error / ZS for real data.
147 
148 
149  const HcalElectronicsMap *emap = pSetup->getHcalMapping();
150 
152  iEvent.getByToken(tok_raw_, fedHandle);
153 
154  if (fedHandle.isValid() && emap != 0) {
155  theAlgo_.runFEFormatError(fedHandle.product(), emap, *result);
156  } else {
157  edm::LogInfo("HcalTrigPrimDigiProducer")
158  << "\nWarning: FEDRawDataCollection with input tag "
159  << inputTagFEDRaw_
160  << "\nrequested in configuration, but not found in the event."
161  << "\nQuit returning empty product." << std::endl;
162 
163  // produce empty HcalTrigPrimDigiCollection and put it in the event
164  std::auto_ptr < HcalTrigPrimDigiCollection > emptyResult(
166 
167  iEvent.put(emptyResult);
168 
169  return;
170  }
171 
172  }
173 
174  if (runZS_) theAlgo_.runZS(*result);
175 
176  // edm::LogInfo("HcalTrigPrimDigiProducer") << "HcalTrigPrims: " << result->size();
177 
178  // Step D: Put outputs into event
179  iEvent.put(result);
180 }
181 
182 
T getParameter(std::string const &) const
edm::InputTag inputTagFEDRaw_
input tag for FEDRawDataCollection
T getUntrackedParameter(std::string const &, T const &) const
void runFEFormatError(const FEDRawDataCollection *rawraw, const HcalElectronicsMap *emap, HcalTrigPrimDigiCollection &result)
HcalTrigPrimDigiProducer(const edm::ParameterSet &ps)
edm::SortedCollection< HcalTriggerPrimitiveDigi > HcalTrigPrimDigiCollection
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
bool exists(std::string const &parameterName) const
checks if a parameter exists
HcalTriggerPrimitiveAlgo theAlgo_
edm::EDGetTokenT< HFDigiCollection > tok_hf_
int iEvent
Definition: GenABIO.cc:230
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:121
tuple result
Definition: query.py:137
bool isValid() const
Definition: HandleBase.h:75
void runZS(HcalTrigPrimDigiCollection &tp)
T const * product() const
Definition: Handle.h:81
const T & get() const
Definition: EventSetup.h:56
T const * product() const
Definition: ESHandle.h:86
virtual void produce(edm::Event &e, const edm::EventSetup &c)
edm::EDGetTokenT< FEDRawDataCollection > tok_raw_
std::vector< edm::InputTag > inputLabel_
input tags for HCAL digis
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
void run(const HcalTPGCoder *incoder, const HcalTPGCompressor *outcoder, const HBHEDigiCollection &hbheDigis, const HFDigiCollection &hfDigis, HcalTrigPrimDigiCollection &result, const HcalTrigTowerGeometry *trigTowerGeometry, float rctlsb, const HcalFeatureBit *LongvrsShortCut=0)