#include <PreshowerClusterShapeProducer.h>

Inheritance diagram for PreshowerClusterShapeProducer:

Public Types
typedef math::XYZPoint	Point

Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Member Functions
	PreshowerClusterShapeProducer (const edm::ParameterSet &ps)
	More...

void	produce (edm::Event &evt, const edm::EventSetup &es) override

	~PreshowerClusterShapeProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
edm::EDGetTokenT< reco::SuperClusterCollection >	endcapSClusterToken_

int	nEvt_

EndcapPiZeroDiscriminatorAlgo *	presh_pi0_algo

edm::EDGetTokenT< EcalRecHitCollection >	preshHitToken_

std::string	PreshowerClusterShapeCollectionX_

std::string	PreshowerClusterShapeCollectionY_

Detailed Description

Definition at line 21 of file PreshowerClusterShapeProducer.h.

Member Typedef Documentation

typedef math::XYZPoint PreshowerClusterShapeProducer::Point

Definition at line 23 of file PreshowerClusterShapeProducer.h.

Constructor & Destructor Documentation

PreshowerClusterShapeProducer::PreshowerClusterShapeProducer ( const edm::ParameterSet & ps )

explicit

Definition at line 40 of file PreshowerClusterShapeProducer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and LogTrace.

                                                                                    {
   // use configuration file to setup input/output collection names
   // Parameters to identify the hit collections
   preshHitToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("preshRecHitProducer"));
   endcapSClusterToken_ =
       consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapSClusterProducer"));
 
   PreshowerClusterShapeCollectionX_ = ps.getParameter<string>("PreshowerClusterShapeCollectionX");
   PreshowerClusterShapeCollectionY_ = ps.getParameter<string>("PreshowerClusterShapeCollectionY");
 
   produces<reco::PreshowerClusterShapeCollection>(PreshowerClusterShapeCollectionX_);
   produces<reco::PreshowerClusterShapeCollection>(PreshowerClusterShapeCollectionY_);
 
   float preshStripECut = ps.getParameter<double>("preshStripEnergyCut");
   int preshNst = ps.getParameter<int>("preshPi0Nstrip");
 
   string debugString = ps.getParameter<string>("debugLevel");
 
   string tmpPath = ps.getUntrackedParameter<string>("pathToWeightFiles", "RecoEcal/EgammaClusterProducers/data/");
 
   presh_pi0_algo = new EndcapPiZeroDiscriminatorAlgo(preshStripECut, preshNst, tmpPath);
 
   LogTrace("EcalClusters") << "PreshowerClusterShapeProducer:presh_pi0_algo class instantiated ";
 
   nEvt_ = 0;
 }

PreshowerClusterShapeProducer::~PreshowerClusterShapeProducer ( )

override

Definition at line 67 of file PreshowerClusterShapeProducer.cc.

67 { delete presh_pi0_algo; }

PreshowerClusterShapeProducer::presh_pi0_algo

EndcapPiZeroDiscriminatorAlgo * presh_pi0_algo

Definition: PreshowerClusterShapeProducer.h:44

Member Function Documentation

void PreshowerClusterShapeProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 69 of file PreshowerClusterShapeProducer.cc.

References edm::SortedCollection< T, SORT >::begin(), DetId::Ecal, EcalPreshower, edm::SortedCollection< T, SORT >::end(), JetChargeProducer_cfi::exp, geometry, edm::EventSetup::get(), edm::Event::getByToken(), CaloGeometry::getSubdetectorGeometry(), edm::EventBase::id(), LogDebug, LogTrace, eostools::move(), edm::Handle< T >::product(), edm::Event::put(), TrackInfoProducer_cfi::rechits, reco::PreshowerClusterShape::setSCRef(), funct::sin(), and edm::SortedCollection< T, SORT >::size().

                                                                             {
   ostringstream ostr;  // use this stream for all messages in produce
 
   LogTrace("EcalClusters") << "\n .......  Event " << evt.id() << " with Number = " << nEvt_ + 1
                            << " is analyzing ....... ";
 
   Handle<EcalRecHitCollection> pRecHits;
   Handle<SuperClusterCollection> pSuperClusters;
 
   // get the ECAL -> Preshower geometry and topology:
   ESHandle<CaloGeometry> geoHandle;
   es.get<CaloGeometryRecord>().get(geoHandle);
   const CaloSubdetectorGeometry* geometry = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
   const CaloSubdetectorGeometry*& geometry_p = geometry;
 
   // create a unique_ptr to a PreshowerClusterShapeCollection
   auto ps_cl_for_pi0_disc_x = std::make_unique<reco::PreshowerClusterShapeCollection>();
   auto ps_cl_for_pi0_disc_y = std::make_unique<reco::PreshowerClusterShapeCollection>();
 
   std::unique_ptr<CaloSubdetectorTopology> topology_p;
   if (geometry)
     topology_p = std::make_unique<EcalPreshowerTopology>();
 
   // fetch the Preshower product (RecHits)
   evt.getByToken(preshHitToken_, pRecHits);
   // pointer to the object in the product
   const EcalRecHitCollection* rechits = pRecHits.product();
 
   LogTrace("EcalClusters") << "PreshowerClusterShapeProducer: ### Total # of preshower RecHits: " << rechits->size();
 
   //  if ( rechits->size() <= 0 ) return;
 
   // make the map of Preshower rechits:
   map<DetId, EcalRecHit> rechits_map;
   EcalRecHitCollection::const_iterator it;
   for (it = rechits->begin(); it != rechits->end(); it++) {
     rechits_map.insert(make_pair(it->id(), *it));
   }
 
   LogTrace("EcalClusters") << "PreshowerClusterShapeProducer: ### Preshower RecHits_map of size " << rechits_map.size()
                            << " was created!";
 
   reco::PreshowerClusterShapeCollection ps_cl_x, ps_cl_y;
 
   //make cycle over Photon Collection
   int SC_index = 0;
   //  Handle<PhotonCollection> correctedPhotonHandle;
   //  evt.getByLabel(photonCorrCollectionProducer_, correctedPhotonCollection_ , correctedPhotonHandle);
   //  const PhotonCollection corrPhoCollection = *(correctedPhotonHandle.product());
   //  cout << " Photon Collection size : " << corrPhoCollection.size() << endl;
 
   evt.getByToken(endcapSClusterToken_, pSuperClusters);
   const reco::SuperClusterCollection* SClusts = pSuperClusters.product();
   LogTrace("EcalClusters") << "### Total # Endcap Superclusters: " << SClusts->size();
 
   SuperClusterCollection::const_iterator it_s;
   for (it_s = SClusts->begin(); it_s != SClusts->end(); it_s++) {
     SuperClusterRef it_super(reco::SuperClusterRef(pSuperClusters, SC_index));
 
     float SC_eta = it_super->eta();
 
     LogTrace("EcalClusters") << "PreshowerClusterShapeProducer: superCl_E = " << it_super->energy()
                              << " superCl_Et = " << it_super->energy() * sin(2 * atan(exp(-it_super->eta())))
                              << " superCl_Eta = " << SC_eta << " superCl_Phi = " << it_super->phi();
 
     if (fabs(SC_eta) >= 1.65 && fabs(SC_eta) <= 2.5) {  //  Use Preshower region only
       if (geometry) {
         const GlobalPoint pointSC(it_super->x(), it_super->y(), it_super->z());  // get the centroid of the SC
         LogTrace("EcalClusters") << "SC centroind = " << pointSC;
 
         // Get the Preshower 2-planes RecHit vectors associated with the given SC
 
         DetId tmp_stripX = (dynamic_cast<const EcalPreshowerGeometry*>(geometry_p))->getClosestCellInPlane(pointSC, 1);
         DetId tmp_stripY = (dynamic_cast<const EcalPreshowerGeometry*>(geometry_p))->getClosestCellInPlane(pointSC, 2);
         ESDetId stripX = (tmp_stripX == DetId(0)) ? ESDetId(0) : ESDetId(tmp_stripX);
         ESDetId stripY = (tmp_stripY == DetId(0)) ? ESDetId(0) : ESDetId(tmp_stripY);
 
         vector<float> vout_stripE1 = presh_pi0_algo->findPreshVector(stripX, &rechits_map, topology_p.get());
         vector<float> vout_stripE2 = presh_pi0_algo->findPreshVector(stripY, &rechits_map, topology_p.get());
 
         LogTrace("EcalClusters") << "PreshowerClusterShapeProducer : ES Energy vector associated to the given SC = ";
         for (int k1 = 0; k1 < 11; k1++) {
           LogTrace("EcalClusters") << vout_stripE1[k1] << " ";
         }
 
         for (int k1 = 0; k1 < 11; k1++) {
           LogTrace("EcalClusters") << vout_stripE2[k1] << " ";
         }
 
         reco::PreshowerClusterShape ps1 = reco::PreshowerClusterShape(vout_stripE1, 1);
         ps1.setSCRef(it_super);
         ps_cl_x.push_back(ps1);
 
         reco::PreshowerClusterShape ps2 = reco::PreshowerClusterShape(vout_stripE2, 2);
         ps2.setSCRef(it_super);
         ps_cl_y.push_back(ps2);
       }
       SC_index++;
     }  // end of cycle over Endcap SC
   }
   // put collection of PreshowerClusterShape in the Event:
   ps_cl_for_pi0_disc_x->assign(ps_cl_x.begin(), ps_cl_x.end());
   ps_cl_for_pi0_disc_y->assign(ps_cl_y.begin(), ps_cl_y.end());
 
   evt.put(std::move(ps_cl_for_pi0_disc_x), PreshowerClusterShapeCollectionX_);
   evt.put(std::move(ps_cl_for_pi0_disc_y), PreshowerClusterShapeCollectionY_);
   LogTrace("EcalClusters") << "PreshowerClusterShapeCollection added to the event";
 
   nEvt_++;
 
   LogDebug("PiZeroDiscriminatorDebug") << ostr.str();
 }