1 #ifndef RecoEgamma_EgammaHLTProducers_HLTCaloObjInRegionsProducer_h_ 2 #define RecoEgamma_EgammaHLTProducers_HLTCaloObjInRegionsProducer_h_ 59 EtaPhiRegion(
float iEta,
float iPhi,
float iDR,
float iDEta,
float iDPhi):
60 centreEta_(iEta),centrePhi_(iPhi),maxDeltaR2_(iDR*iDR),
61 maxDEta_(iDEta),maxDPhi_(iDPhi){}
64 return reco::deltaR2(eta,phi,centreEta_,centrePhi_)<maxDeltaR2_ ||
73 virtual void getEtaPhiRegions(
const edm::Event&,std::vector<EtaPhiRegion>&)
const=0;
89 minEt_(para.getParameter<double>(
"minEt")),
90 maxEt_(para.getParameter<double>(
"maxEt")),
91 maxDeltaR_(para.getParameter<double>(
"maxDeltaR")),
92 maxDEta_(para.getParameter<double>(
"maxDEta")),
93 maxDPhi_(para.getParameter<double>(
"maxDPhi")),
94 token_(consumesColl.consumes<T1>(para.getParameter<
edm::InputTag>(
"inputColl"))){}
96 void getEtaPhiRegions(
const edm::Event&,std::vector<EtaPhiRegion>&)
const override;
97 template<
typename T2>
static typename T2::const_iterator
beginIt(
const T2&
coll){
return coll.begin();}
98 template<
typename T2>
static typename T2::const_iterator
endIt(
const T2&
coll){
return coll.end();}
105 template<
typename CaloObjType,
122 static std::unique_ptr<CaloObjCollType>
125 const std::vector<EtaPhiRegion>& regions);
126 static bool validIDForGeom(
const DetId&
id);
134 template<
typename CaloObjType,
typename CaloObjCollType>
137 const std::vector<edm::ParameterSet> etaPhiRegions = para.
getParameter<std::vector<edm::ParameterSet>>(
"etaPhiRegions");
138 for(
auto&
pset : etaPhiRegions){
140 etaPhiRegionData_.emplace_back(createEtaPhiRegionData(type,
pset,consumesCollector()));
143 outputProductNames_=para.
getParameter<std::vector<std::string>>(
"outputProductNames");
144 inputCollTags_=para.
getParameter<std::vector<edm::InputTag>>(
"inputCollTags");
145 if(outputProductNames_.size()!=inputCollTags_.size()){
146 throw cms::Exception(
"InvalidConfiguration") <<
" error outputProductNames and inputCollTags must be the same size, they are "<<outputProductNames_.size()<<
" vs "<<inputCollTags_.size();
148 for (
unsigned int collNr=0; collNr<inputCollTags_.size(); collNr++) {
149 inputTokens_.push_back(consumes<CaloObjCollType>(inputCollTags_[collNr]));
150 produces<CaloObjCollType> (outputProductNames_[collNr]);
154 template<
typename CaloObjType,
typename CaloObjCollType>
158 std::vector<std::string> outputProductNames;
159 outputProductNames.push_back(
"EcalRegionalRecHitsEB");
160 desc.
add<std::vector<std::string>>(
"outputProductNames", outputProductNames);
161 std::vector<edm::InputTag> inputColls;
163 desc.
add<std::vector<edm::InputTag>>(
"inputCollTags", inputColls);
164 std::vector<edm::ParameterSet> etaPhiRegions;
174 etaPhiRegions.push_back(ecalCandPSet);
179 etaPhiRegionDesc.
add<
double>(
"minEt");
180 etaPhiRegionDesc.
add<
double>(
"maxEt");
181 etaPhiRegionDesc.
add<
double>(
"maxDeltaR");
182 etaPhiRegionDesc.
add<
double>(
"maxDEta");
183 etaPhiRegionDesc.
add<
double>(
"maxDPhi");
185 desc.
addVPSet(
"etaPhiRegions",etaPhiRegionDesc,etaPhiRegions);
191 template<
typename CaloObjType,
typename CaloObjCollType>
198 std::vector<EtaPhiRegion> regions;
199 std::for_each(etaPhiRegionData_.begin(),etaPhiRegionData_.end(),
200 [&
event,®ions](
const std::unique_ptr<EtaPhiRegionDataBase>&
input)
201 {
input->getEtaPhiRegions(event,regions);}
204 for(
size_t inputCollNr=0;inputCollNr<inputTokens_.size();inputCollNr++){
206 event.getByToken(inputTokens_[inputCollNr],inputColl);
209 edm::LogError(
"ProductNotFound")<<
"could not get a handle on the "<<
typeid(CaloObjCollType).
name() <<
" named "<< inputCollTags_[inputCollNr].encode() << std::endl;
212 auto outputColl = makeFilteredColl(inputColl,caloGeomHandle,regions);
217 template<
typename CaloObjType,
typename CaloObjCollType>
218 std::unique_ptr<CaloObjCollType>
222 const std::vector<EtaPhiRegion>& regions)
225 auto outputColl = std::make_unique<CaloObjCollType>();
226 if(!inputColl->empty()){
228 if(!regions.empty()){
229 for(
const CaloObjType&
obj : *inputColl){
231 if(objGeom==
nullptr){
237 if(validIDForGeom(
obj.id())) {
238 edm::LogError(
"HLTCaloObjInRegionsProducer") <<
"for an object of type "<<
typeid(CaloObjType).
name()<<
" the geometry returned null for id "<<
DetId(
obj.id()).rawId()<<
" with initial ID "<<
DetId(inputColl->begin()->id()).rawId()<<
" in HLTCaloObjsInRegion, this shouldnt be possible and something has gone wrong, auto accepting hit";
243 float eta = objGeom->getPosition().eta();
244 float phi = objGeom->getPosition().phi();
246 for(
const auto& region : regions){
247 if(region(eta,phi)) {
263 template<
typename CaloObjType,
typename CaloObjCollType>
279 template<
typename CaloObjType,
typename CaloObjCollType>
285 if(type==
"L1EGamma"){
287 }
else if(type==
"L1Jet"){
289 }
else if(type==
"L1Muon"){
291 }
else if(type==
"L1Tau"){
293 }
else if(type==
"RecoEcalCandidate"){
295 }
else if(type==
"RecoChargedCandidate"){
297 }
else if(type==
"Electron"){
301 throw cms::Exception(
"InvalidConfig") <<
" type "<<type<<
" is not recognised, this means the rec-hit you think you are keeping may not be and you should fix this error as it can lead to hard to find efficiency loses"<<std::endl;
306 template<
typename CandCollType>
310 event.getByToken(token_,cands);
312 for(
auto candIt = beginIt(*cands);candIt!=endIt(*cands);++candIt){
313 if(candIt->et() >= minEt_ && (maxEt_<0 || candIt->et() < maxEt_)){
314 regions.push_back(
EtaPhiRegion(candIt->eta(),candIt->phi(),
const_iterator end(int bx) const
T getParameter(std::string const &) const
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
ParameterDescriptionBase * addVPSet(U const &iLabel, ParameterSetDescription const &validator, std::vector< ParameterSet > const &defaults)
std::vector< edm::EDGetTokenT< CaloObjCollType > > inputTokens_
std::vector< edm::InputTag > inputCollTags_
EtaPhiRegionDataBase * createEtaPhiRegionData(const std::string &, const edm::ParameterSet &, edm::ConsumesCollector &&)
def setup(process, global_tag, zero_tesla=False)
~HLTCaloObjInRegionsProducer() override
static BXVector< T2 >::const_iterator beginIt(const BXVector< T2 > &coll)
std::string defaultModuleLabel()
static T2::const_iterator beginIt(const T2 &coll)
static bool validIDForGeom(const DetId &id)
static std::string const input
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
std::vector< std::string > outputProductNames_
static std::unique_ptr< CaloObjCollType > makeFilteredColl(const edm::Handle< CaloObjCollType > &inputColl, const edm::ESHandle< CaloGeometry > &caloGeomHandle, const std::vector< EtaPhiRegion > ®ions)
static T2::const_iterator endIt(const T2 &coll)
void addParameter(std::string const &name, T const &value)
bool operator()(float eta, float phi) const
Abs< T >::type abs(const T &t)
EtaPhiRegionData(const edm::ParameterSet ¶, edm::ConsumesCollector &consumesColl)
SubDetector subDetGeom[21]
ParameterDescriptionBase * add(U const &iLabel, T const &value)
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
HLTCaloObjInRegionsProducer(const edm::ParameterSet &ps)
double deltaPhi(double phi1, double phi2)
edm::EDGetTokenT< T1 > token_
virtual std::shared_ptr< const CaloCellGeometry > getGeometry(const DetId &id) const
Get the cell geometry of a given detector id. Should return false if not found.
void add(std::string const &label, ParameterSetDescription const &psetDescription)
T1 deltaR2(T1 eta1, T2 phi1, T3 eta2, T4 phi2)
static BXVector< T2 >::const_iterator endIt(const BXVector< T2 > &coll)
void getEtaPhiRegions(const edm::Event &, std::vector< EtaPhiRegion > &) const override
const_iterator begin(int bx) const
EtaPhiRegion(float iEta, float iPhi, float iDR, float iDEta, float iDPhi)
std::vector< std::unique_ptr< EtaPhiRegionDataBase > > etaPhiRegionData_
void produce(edm::Event &, const edm::EventSetup &) override