102 auto clusters_p1 = std::make_unique<reco::PreshowerClusterCollection>();
103 auto clusters_p2 = std::make_unique<reco::PreshowerClusterCollection>();
105 auto superclusters_p = std::make_unique<reco::SuperClusterCollection>();
120 LogTrace(
"EcalClusters") <<
"PreshowerClusterProducerInfo: ### Total # of preshower RecHits: "<< rechits->
size();
123 std::map<DetId, EcalRecHit> rechits_map;
125 for (it = rechits->
begin(); it != rechits->
end(); it++) {
127 if (it->recoFlag()==1 || it->recoFlag()==14 || (it->recoFlag()<=10 && it->recoFlag()>=5))
continue;
129 rechits_map.insert(std::make_pair(it->id(), *it));
132 std::set<DetId> used_strips;
134 LogTrace(
"EcalClusters") <<
"PreshowerClusterProducerInfo: ### rechits_map of size " << rechits_map.size() <<
" was created!";
141 reco::SuperClusterCollection::const_iterator it_super;
143 for (it_super=SClusts->begin(); it_super!=SClusts->end(); ++it_super) {
151 LogTrace(
"EcalClusters")<<
" superE = " << it_super->energy() <<
" superETA = " << it_super->eta() <<
" superPHI = " << it_super->phi() ;
157 for ( ; bc_iter !=it_super->clustersEnd(); ++bc_iter ) {
161 double X = (*bc_iter)->x();
162 double Y = (*bc_iter)->y();
163 double Z = (*bc_iter)->z();
175 if (status_p1->getStatusCode() == 1) condP1 = 0;
176 if (status_p2->getStatusCode() == 1) condP2 = 0;
177 }
else if (strip1 ==
ESDetId(0))
188 clusters1.push_back(cl1);
195 clusters2.push_back(cl2);
205 LogTrace(
"EcalClusters") <<
" For SC #" << isc-1 <<
", containing " 206 << it_super->clustersSize()
207 <<
" basic clusters, PreshowerClusterAlgo made " 208 << clusters1.size() <<
" in X plane and " 210 <<
" in Y plane " <<
" preshower clusters " ;
213 if(e1+e2 > 1.0
e-10) {
218 if (condP1 == 1 && condP2 == 1) {
220 }
else if (condP1 == 1 && condP2 == 0) {
222 }
else if (condP1 == 0 && condP2 == 1) {
224 }
else if (condP1 == 0 && condP2 == 0) {
230 float E = it_super->energy() + deltaE;
232 LogTrace(
"EcalClusters") <<
" Creating corrected SC ";
235 if (condP1 == 1 && condP2 == 1)
sc.setPreshowerPlanesStatus(0);
236 else if (condP1 == 1 && condP2 == 0)
sc.setPreshowerPlanesStatus(1);
237 else if (condP1 == 0 && condP2 == 1)
sc.setPreshowerPlanesStatus(2);
238 else if (condP1 == 0 && condP2 == 0)
sc.setPreshowerPlanesStatus(3);
241 new_SC.push_back(
sc);
242 LogTrace(
"EcalClusters") <<
" SuperClusters energies: new E = " <<
sc.
energy() <<
" vs. old E =" << it_super->energy();
247 clusters_p1->assign(clusters1.begin(), clusters1.end());
248 clusters_p2->assign(clusters2.begin(), clusters2.end());
252 LogTrace(
"EcalClusters") <<
"Preshower clusters added to the event" ;
255 superclusters_p->assign(new_SC.begin(), new_SC.end());
257 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, 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
ESHandle< TrackerGeometry > geometry
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