104 std::unique_ptr<const CaloSubdetectorTopology>
topology;
120 std::vector<RectangularEtaPhiRegion>
regions;
130 auto uhits = std::make_unique<EcalUncalibratedRecHitCollection>();
133 if (!(urhcH[i].isValid())) {
134 edm::LogError(
"ProductNotFound")<<
"could not get a handle on the EcalRecHitCollection! (" <<
hitLabels[
i].encode() <<
")" << std::endl;
139 if (!uncalibRecHits->
empty()) {
140 if ((*uncalibRecHits)[0].id().subdetId() ==
EcalBarrel) {
142 topology = std::make_unique<EcalBarrelTopology>(*geoHandle);
143 }
else if ((*uncalibRecHits)[0].
id().subdetId() ==
EcalEndcap) {
145 topology = std::make_unique<EcalEndcapTopology>(*geoHandle);
146 }
else if ((*uncalibRecHits)[0].
id().subdetId() ==
EcalPreshower) {
148 topology = std::make_unique<EcalPreshowerTopology>();
149 }
else throw(std::runtime_error(
"\n\nProducer encountered invalied ecalhitcollection type.\n\n"));
151 if(!regions.empty()) {
154 for (it = uncalibRecHits->
begin(); it != uncalibRecHits->
end(); it++){
155 auto this_cell = geometry_p->
getGeometry(it->id());
157 std::vector<RectangularEtaPhiRegion>::const_iterator region;
158 for (region=regions.begin(); region!=regions.end(); region++) {
159 if (region->inRegion(this_cell->etaPos(),this_cell->phiPos())) {
160 uhits->push_back(*it);
174 for (
unsigned int i=0; i<
hitLabels.size(); i++) {
175 auto hits = std::make_unique<EcalRecHitCollection>();
178 if (!(rhcH[i].isValid())) {
179 edm::LogError(
"ProductNotFound")<<
"could not get a handle on the EcalRecHitCollection! (" <<
hitLabels[
i].encode() <<
")" << std::endl;
184 if (!recHits->
empty()) {
185 if ((*recHits)[0].id().subdetId() ==
EcalBarrel) {
187 topology = std::make_unique<EcalBarrelTopology>(*geoHandle);
188 }
else if ((*recHits)[0].
id().subdetId() ==
EcalEndcap) {
190 topology = std::make_unique<EcalEndcapTopology>(*geoHandle);
193 topology = std::make_unique<EcalPreshowerTopology>();
194 }
else throw(std::runtime_error(
"\n\nProducer encountered invalied ecalhitcollection type.\n\n"));
196 if(!regions.empty()) {
198 for (it = recHits->
begin(); it != recHits->
end(); it++){
199 auto this_cell = geometry_p->
getGeometry(it->id());
201 std::vector<RectangularEtaPhiRegion>::const_iterator region;
202 for (region=regions.begin(); region!=regions.end(); region++) {
203 if (region->inRegion(this_cell->etaPos(),this_cell->phiPos())) {
204 hits->push_back(*it);
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
const std::vector< edm::InputTag > hitLabels
CaloTopology const * topology(0)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
const std::vector< std::string > productLabels
std::vector< T >::const_iterator const_iterator
const edm::EDGetTokenT< T1Collection > l1TokenIsolated_
std::vector< edm::EDGetTokenT< EcalRecHitCollection > > hitTokens
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
const_iterator end() const
T const * product() const
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 getEtaPhiRegions(std::vector< RectangularEtaPhiRegion > *, T1Collection, const L1CaloGeometry &, bool)
const edm::EDGetTokenT< T1Collection > l1TokenNonIsolated_
const double l1LowerThrIgnoreIsolation_
const_iterator begin() const
std::vector< edm::EDGetTokenT< EcalUncalibratedRecHitCollection > > uncalibHitTokens