101 auto clusters_p1 = std::make_unique<reco::PreshowerClusterCollection>();
102 auto clusters_p2 = std::make_unique<reco::PreshowerClusterCollection>();
104 auto superclusters_p = std::make_unique<reco::SuperClusterCollection>();
106 std::unique_ptr<CaloSubdetectorTopology> topology_p;
108 topology_p = std::make_unique<EcalPreshowerTopology>();
119 LogTrace(
"EcalClusters") <<
"PreshowerClusterProducerInfo: ### Total # of preshower RecHits: "<< rechits->
size();
122 std::map<DetId, EcalRecHit> rechits_map;
124 for (it = rechits->
begin(); it != rechits->
end(); it++) {
126 if (it->recoFlag()==1 || it->recoFlag()==14 || (it->recoFlag()<=10 && it->recoFlag()>=5))
continue;
128 rechits_map.insert(std::make_pair(it->id(), *it));
131 std::set<DetId> used_strips;
133 LogTrace(
"EcalClusters") <<
"PreshowerClusterProducerInfo: ### rechits_map of size " << rechits_map.size() <<
" was created!";
140 reco::SuperClusterCollection::const_iterator it_super;
142 for (it_super=SClusts->begin(); it_super!=SClusts->end(); ++it_super) {
150 LogTrace(
"EcalClusters")<<
" superE = " << it_super->energy() <<
" superETA = " << it_super->eta() <<
" superPHI = " << it_super->phi() ;
156 for ( ; bc_iter !=it_super->clustersEnd(); ++bc_iter ) {
160 double X = (*bc_iter)->x();
161 double Y = (*bc_iter)->y();
162 double Z = (*bc_iter)->z();
174 if (status_p1->getStatusCode() == 1) condP1 = 0;
175 if (status_p2->getStatusCode() == 1) condP2 = 0;
176 }
else if (strip1 ==
ESDetId(0))
187 clusters1.push_back(cl1);
194 clusters2.push_back(cl2);
204 LogTrace(
"EcalClusters") <<
" For SC #" << isc-1 <<
", containing " 205 << it_super->clustersSize()
206 <<
" basic clusters, PreshowerClusterAlgo made " 207 << clusters1.size() <<
" in X plane and " 209 <<
" in Y plane " <<
" preshower clusters " ;
212 if(e1+e2 > 1.0
e-10) {
217 if (condP1 == 1 && condP2 == 1) {
219 }
else if (condP1 == 1 && condP2 == 0) {
221 }
else if (condP1 == 0 && condP2 == 1) {
223 }
else if (condP1 == 0 && condP2 == 0) {
229 float E = it_super->energy() + deltaE;
231 LogTrace(
"EcalClusters") <<
" Creating corrected SC ";
234 if (condP1 == 1 && condP2 == 1)
sc.setPreshowerPlanesStatus(0);
235 else if (condP1 == 1 && condP2 == 0)
sc.setPreshowerPlanesStatus(1);
236 else if (condP1 == 0 && condP2 == 1)
sc.setPreshowerPlanesStatus(2);
237 else if (condP1 == 0 && condP2 == 0)
sc.setPreshowerPlanesStatus(3);
240 new_SC.push_back(
sc);
241 LogTrace(
"EcalClusters") <<
" SuperClusters energies: new E = " <<
sc.
energy() <<
" vs. old E =" << it_super->energy();
246 clusters_p1->assign(clusters1.begin(), clusters1.end());
247 clusters_p2->assign(clusters2.begin(), clusters2.end());
251 LogTrace(
"EcalClusters") <<
"Preshower clusters added to the event" ;
254 superclusters_p->assign(new_SC.begin(), new_SC.end());
256 LogTrace(
"EcalClusters") <<
"Corrected SClusters added to the event" ;
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.
bool getByToken(EDGetToken token, Handle< PROD > &result) const
void push_back(Ptr< T > const &iPtr)
const self & getMap() const
std::string preshClusterCollectionY_
Sin< T >::type sin(const T &t)
std::vector< EcalRecHit >::const_iterator const_iterator
PreshowerClusterAlgo * presh_algo
reco::PreshowerCluster makeOneCluster(ESDetId strip, HitsID *used_strips, RecHitsMap *rechits_map, const CaloSubdetectorGeometry *geometry_p, const CaloSubdetectorTopology *topology_p)
std::vector< SuperCluster > SuperClusterCollection
collection of SuperCluser objectr
edm::ESHandle< ESChannelStatus > esChannelStatus_
void setBCRef(const CaloClusterPtr &r)
DetIds of component RecHits – now inherited from CaloCluster.
const_iterator find(uint32_t rawId) const
float energy() const
Energy. Note this is taken from the first SimTrack only.
std::vector< PreshowerCluster > PreshowerClusterCollection
collection of PreshowerCluster objects
double energy() const
cluster energy
float theta() const
Momentum polar angle. Note this is taken from the first SimTrack only.
const_iterator end() const
std::string preshClusterCollectionX_
edm::EDGetTokenT< EcalRecHitCollection > preshHitsToken_
T const * product() const
std::vector< Item >::const_iterator const_iterator
std::string assocSClusterCollection_
edm::EDGetTokenT< reco::SuperClusterCollection > endcapSClusterToken_
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
const_iterator begin() const