46 #include <TLorentzVector.h> 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")),
74 TS4TS5EnergyThreshold_(iConfig.getParameter<double>(
"TS4TS5EnergyThreshold"))
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");
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]));
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]));
93 if(iConfig.
existsAs<
double>(
"minR45HitE"))
96 produces<CaloTowerCollection>();
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);
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);
141 descriptions.
add(
"hltHcalTowerNoiseCleaner",desc);
150 using namespace reco;
156 std::set<unsigned int> noisyTowers;
159 edm::LogError(
"HLTHcalTowerNoiseCleaner") <<
"Input Tower Collection is not Valid";
171 edm::LogWarning(
"HLTHcalTowerNoiseCleaner") <<
"Could not find HcalNoiseRBXCollection product named " 178 for(
auto const & rbx : *rbxs_h) {
187 for(
auto const & it : data) {
190 bool passFilter=
true;
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;
198 else if(it.numZeros()>=
minZeros_) passFilter=
false;
201 else if(!it.PassTS4TS5()) passFilter=
false;
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()
225 for( noiseTowersIt = noiseTowers.
begin(); noiseTowersIt != noiseTowers.
end(); noiseTowersIt++){
228 noisyTowers.insert( caloTowerTopology->
denseIndex(
id) );
238 for(inTowersIt = tower_h->
begin(); inTowersIt != tower_h->
end(); inTowersIt++){
241 if(noisyTowers.find( caloTowerTopology->
denseIndex(
id) ) == noisyTowers.end()){
242 OutputTowers->push_back(*inTowersIt);
T getParameter(std::string const &) const
edm::InputTag TowerCollectionTag_
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
bool existsAs(std::string const ¶meterName, bool trackiness=true) const
checks if a parameter exists as a given type
HLTHcalTowerNoiseCleaner(const edm::ParameterSet &)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
std::set< CommonHcalNoiseRBXData, noisedatacomp > noisedataset_t
std::vector< CaloTower >::const_iterator const_iterator
const_iterator end() const
Termination of iteration.
ProductID id() const
Accessor for product ID.
const_iterator begin() const
Initialize an iterator over the RefVector.
std::vector< std::pair< double, double > > TS4TS5UpperCut_
edm::EDGetTokenT< CaloTowerCollection > m_theCaloTowerCollectionToken
edm::SortedCollection< CaloTower > CaloTowerCollection
ParameterDescriptionBase * add(U const &iLabel, T const &value)
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
edm::InputTag HcalNoiseRBXCollectionTag_
const_iterator end() const
uint32_t denseIndex(const DetId &id) const
virtual void produce(edm::Event &, const edm::EventSetup &)
CaloTowerDetId id() const
edm::EDGetTokenT< reco::HcalNoiseRBXCollection > m_theHcalNoiseToken
void add(std::string const &label, ParameterSetDescription const &psetDescription)
~HLTHcalTowerNoiseCleaner()
double TS4TS5EnergyThreshold_
char data[epos_bytes_allocation]
std::vector< std::pair< double, double > > TS4TS5LowerCut_
const_iterator begin() const