CMS 3D CMS Logo

HLTHcalTowerNoiseCleaner.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Class: HLTHcalTowerNoiseCleaner
4 //
12 //
13 // Original Author: Alexander Mott
14 // Created: Mon Nov 21 11:32:00 CEST 2011
15 //
16 //
17 //
18 
19 
25 
27 
37 
41 
42 #include <iostream>
43 #include <string>
44 #include <fstream>
45 #include <TVector3.h>
46 #include <TLorentzVector.h>
47 #include <set>
48 
50 
51 //#include <Point.h>
52 
54  : HcalNoiseRBXCollectionTag_(iConfig.getParameter<edm::InputTag>("HcalNoiseRBXCollection")),
55  TowerCollectionTag_(iConfig.getParameter<edm::InputTag>("CaloTowerCollection")),
56  severity_(iConfig.getParameter<int> ("severity")),
57  maxNumRBXs_(iConfig.getParameter<int>("maxNumRBXs")),
58  numRBXsToConsider_(iConfig.getParameter<int>("numRBXsToConsider")),
59  needEMFCoincidence_(iConfig.getParameter<bool>("needEMFCoincidence")),
60  minRBXEnergy_(iConfig.getParameter<double>("minRBXEnergy")),
61  minRatio_(iConfig.getParameter<double>("minRatio")),
62  maxRatio_(iConfig.getParameter<double>("maxRatio")),
63  minHPDHits_(iConfig.getParameter<int>("minHPDHits")),
64  minRBXHits_(iConfig.getParameter<int>("minRBXHits")),
65  minHPDNoOtherHits_(iConfig.getParameter<int>("minHPDNoOtherHits")),
66  minZeros_(iConfig.getParameter<int>("minZeros")),
67  minHighEHitTime_(iConfig.getParameter<double>("minHighEHitTime")),
68  maxHighEHitTime_(iConfig.getParameter<double>("maxHighEHitTime")),
69  maxRBXEMF_(iConfig.getParameter<double>("maxRBXEMF")),
70  minRecHitE_(iConfig.getParameter<double>("minRecHitE")),
71  minLowHitE_(iConfig.getParameter<double>("minLowHitE")),
72  minHighHitE_(iConfig.getParameter<double>("minHighHitE")),
73  minR45HitE_(5.0),
74  TS4TS5EnergyThreshold_(iConfig.getParameter<double>("TS4TS5EnergyThreshold"))
75 {
76 
77  std::vector<double> TS4TS5UpperThresholdTemp = iConfig.getParameter<std::vector<double> >("TS4TS5UpperThreshold");
78  std::vector<double> TS4TS5UpperCutTemp = iConfig.getParameter<std::vector<double> >("TS4TS5UpperCut");
79  std::vector<double> TS4TS5LowerThresholdTemp = iConfig.getParameter<std::vector<double> >("TS4TS5LowerThreshold");
80  std::vector<double> TS4TS5LowerCutTemp = iConfig.getParameter<std::vector<double> >("TS4TS5LowerCut");
81 
82  for(int i = 0; i < (int)TS4TS5UpperThresholdTemp.size() && i < (int)TS4TS5UpperCutTemp.size(); i++)
83  TS4TS5UpperCut_.push_back(std::pair<double, double>(TS4TS5UpperThresholdTemp[i], TS4TS5UpperCutTemp[i]));
84  sort(TS4TS5UpperCut_.begin(), TS4TS5UpperCut_.end());
85 
86  for(int i = 0; i < (int)TS4TS5LowerThresholdTemp.size() && i < (int)TS4TS5LowerCutTemp.size(); i++)
87  TS4TS5LowerCut_.push_back(std::pair<double, double>(TS4TS5LowerThresholdTemp[i], TS4TS5LowerCutTemp[i]));
88  sort(TS4TS5LowerCut_.begin(), TS4TS5LowerCut_.end());
89 
90  m_theHcalNoiseToken = consumes<reco::HcalNoiseRBXCollection>(HcalNoiseRBXCollectionTag_);
91  m_theCaloTowerCollectionToken = consumes<CaloTowerCollection>(TowerCollectionTag_);
92 
93  if(iConfig.existsAs<double>("minR45HitE"))
94  minR45HitE_ = iConfig.getParameter<double>("minR45HitE");
95 
96  produces<CaloTowerCollection>();
97 }
98 
99 
101 
102 void
105  desc.add<edm::InputTag>("HcalNoiseRBXCollection",edm::InputTag("hltHcalNoiseInfoProducer"));
106  desc.add<edm::InputTag>("CaloTowerCollection",edm::InputTag("hltTowerMakerForAll"));
107  desc.add<double>("maxTowerNoiseEnergyFraction",0.5);
108  desc.add<int>("severity",1);
109  desc.add<int>("maxNumRBXs",2);
110  desc.add<int>("numRBXsToConsider",2);
111  desc.add<bool>("needEMFCoincidence",true);
112  desc.add<double>("minRBXEnergy",50.0);
113  desc.add<double>("minRatio",-999.);
114  desc.add<double>("maxRatio",999.);
115  desc.add<int>("minHPDHits",17);
116  desc.add<int>("minRBXHits",999);
117  desc.add<int>("minHPDNoOtherHits",10);
118  desc.add<int>("minZeros",10);
119  desc.add<double>("minHighEHitTime",-9999.0);
120  desc.add<double>("maxHighEHitTime",9999.0);
121  desc.add<double>("maxRBXEMF",0.02);
122  desc.add<double>("minRecHitE",1.5);
123  desc.add<double>("minLowHitE",10.0);
124  desc.add<double>("minHighHitE",25.0);
125  desc.add<double>("minR45HitE",5.0);
126  desc.add<double>("TS4TS5EnergyThreshold",50.0);
127 
128  double TS4TS5UpperThresholdArray[5] = {70, 90, 100, 400, 4000 };
129  double TS4TS5UpperCutArray[5] = {1, 0.8, 0.75, 0.72, 0.72};
130  double TS4TS5LowerThresholdArray[7] = {100, 120, 150, 200, 300, 400, 500};
131  double TS4TS5LowerCutArray[7] = {-1, -0.7, -0.4, -0.2, -0.08, 0, 0.1};
132  std::vector<double> TS4TS5UpperThreshold(TS4TS5UpperThresholdArray, TS4TS5UpperThresholdArray+5);
133  std::vector<double> TS4TS5UpperCut(TS4TS5UpperCutArray, TS4TS5UpperCutArray+5);
134  std::vector<double> TS4TS5LowerThreshold(TS4TS5LowerThresholdArray, TS4TS5LowerThresholdArray+7);
135  std::vector<double> TS4TS5LowerCut(TS4TS5LowerCutArray, TS4TS5LowerCutArray+7);
136 
137  desc.add<std::vector<double> >("TS4TS5UpperThreshold", TS4TS5UpperThreshold);
138  desc.add<std::vector<double> >("TS4TS5UpperCut", TS4TS5UpperCut);
139  desc.add<std::vector<double> >("TS4TS5LowerThreshold", TS4TS5LowerThreshold);
140  desc.add<std::vector<double> >("TS4TS5LowerCut", TS4TS5LowerCut);
141  descriptions.add("hltHcalTowerNoiseCleaner",desc);
142 }
143 
144 //
145 // member functions
146 //
147 
149 {
150  using namespace reco;
151 
152  //get the calo MET / MHT
155 
156  std::set<unsigned int> noisyTowers;
157 
158  if(not tower_h.isValid()){ //No towers MET, don't do anything and accept the event
159  edm::LogError("HLTHcalTowerNoiseCleaner") << "Input Tower Collection is not Valid";
160  return;
161  }
162 
163  //get the calotower topology
164  edm::ESHandle<CaloTowerTopology> caloTowerTopology;
165  iSetup.get<HcalRecNumberingRecord>().get(caloTowerTopology);
166 
167  // get the RBXs produced by RecoMET/METProducers/HcalNoiseInfoProducer
169  iEvent.getByToken(m_theHcalNoiseToken,rbxs_h);
170  if(!rbxs_h.isValid()) {
171  edm::LogWarning("HLTHcalTowerNoiseCleaner") << "Could not find HcalNoiseRBXCollection product named "
172  << HcalNoiseRBXCollectionTag_ << "." << std::endl;
173  severity_=0;
174  }
175 
176  // create a sorted set of the RBXs, ordered by energy
178  for(auto const & rbx : *rbxs_h) {
181  data.insert(d);
182  }
183 
184  // data is now sorted by RBX energy
185  // only consider top N=numRBXsToConsider_ energy RBXs
186  if(severity_>0){
187  for(auto const & it : data) {
188 
189 
190  bool passFilter=true;
191  bool passEMF=true;
192  if(it.energy()>minRBXEnergy_) {
193  if(it.validRatio() && it.ratio()<minRatio_) passFilter=false;
194  else if(it.validRatio() && it.ratio()>maxRatio_) passFilter=false;
195  else if(it.numHPDHits()>=minHPDHits_) passFilter=false;
196  else if(it.numRBXHits()>=minRBXHits_) passFilter=false;
197  else if(it.numHPDNoOtherHits()>=minHPDNoOtherHits_) passFilter=false;
198  else if(it.numZeros()>=minZeros_) passFilter=false;
199  else if(it.minHighEHitTime()<minHighEHitTime_) passFilter=false;
200  else if(it.maxHighEHitTime()>maxHighEHitTime_) passFilter=false;
201  else if(!it.PassTS4TS5()) passFilter=false;
202 
203  if(it.RBXEMF()<maxRBXEMF_){
204  passEMF=false;
205  }
206  }
207 
208  if((needEMFCoincidence_ && !passEMF && !passFilter) ||
209  (!needEMFCoincidence_ && !passFilter)) { // check for noise
210  LogDebug("") << "HLTHcalTowerNoiseCleaner debug: Found a noisy RBX: "
211  << "energy=" << it.energy() << "; "
212  << "ratio=" << it.ratio() << "; "
213  << "# RBX hits=" << it.numRBXHits() << "; "
214  << "# HPD hits=" << it.numHPDHits() << "; "
215  << "# Zeros=" << it.numZeros() << "; "
216  << "min time=" << it.minHighEHitTime() << "; "
217  << "max time=" << it.maxHighEHitTime() << "; "
218  << "passTS4TS5=" << it.PassTS4TS5() << "; "
219  << "RBX EMF=" << it.RBXEMF()
220  << std::endl;
221  // add calotowers associated with this RBX to the noise list
222  edm::RefVector<CaloTowerCollection> noiseTowers = it.rbxTowers();
224  //add these calotowers to the noisy list
225  for( noiseTowersIt = noiseTowers.begin(); noiseTowersIt != noiseTowers.end(); noiseTowersIt++){
226  edm::Ref<edm::SortedCollection<CaloTower> > tower_ref = *noiseTowersIt;
227  CaloTowerDetId id = tower_ref->id();
228  noisyTowers.insert( caloTowerTopology->denseIndex(id) );
229  }}
230  } // done with noise loop
231  }//if(severity_>0)
232 
233  //output collection
234  std::unique_ptr<CaloTowerCollection> OutputTowers(new CaloTowerCollection() );
235 
237 
238  for(inTowersIt = tower_h->begin(); inTowersIt != tower_h->end(); inTowersIt++){
239  const CaloTower & tower = (*inTowersIt);
240  CaloTowerDetId id = tower.id();
241  if(noisyTowers.find( caloTowerTopology->denseIndex(id) ) == noisyTowers.end()){ // the tower is not noisy
242  OutputTowers->push_back(*inTowersIt);
243  }
244  }
245  iEvent.put(std::move(OutputTowers));
246 
247 }
248 
249 
#define LogDebug(id)
T getParameter(std::string const &) const
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:186
HLTHcalTowerNoiseCleaner(const edm::ParameterSet &)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
std::set< CommonHcalNoiseRBXData, noisedatacomp > noisedataset_t
std::vector< CaloTower >::const_iterator const_iterator
const_iterator end() const
Termination of iteration.
Definition: RefVector.h:253
ProductID id() const
Accessor for product ID.
Definition: Ref.h:258
const_iterator begin() const
Initialize an iterator over the RefVector.
Definition: RefVector.h:248
std::vector< std::pair< double, double > > TS4TS5UpperCut_
int iEvent
Definition: GenABIO.cc:230
edm::EDGetTokenT< CaloTowerCollection > m_theCaloTowerCollectionToken
edm::SortedCollection< CaloTower > CaloTowerCollection
Definition: CaloTowerDefs.h:16
ParameterDescriptionBase * add(U const &iLabel, T const &value)
bool isValid() const
Definition: HandleBase.h:74
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
const_iterator end() const
uint32_t denseIndex(const DetId &id) const
virtual void produce(edm::Event &, const edm::EventSetup &)
CaloTowerDetId id() const
Definition: CaloTower.h:103
const T & get() const
Definition: EventSetup.h:56
edm::EDGetTokenT< reco::HcalNoiseRBXCollection > m_theHcalNoiseToken
void add(std::string const &label, ParameterSetDescription const &psetDescription)
fixed size matrix
HLT enums.
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
std::vector< std::pair< double, double > > TS4TS5LowerCut_
def move(src, dest)
Definition: eostools.py:510
const_iterator begin() const