#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)

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	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

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 24 of file PreshowerClusterShapeProducer.h.

Member Typedef Documentation

◆ Point

typedef math::XYZPoint PreshowerClusterShapeProducer::Point

Definition at line 26 of file PreshowerClusterShapeProducer.h.

Constructor & Destructor Documentation

◆ PreshowerClusterShapeProducer()

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

explicit

Definition at line 37 of file PreshowerClusterShapeProducer.cc.

                                                                                    {
   // 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"));
   caloGeometryToken_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
  
   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;
 }

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

◆ ~PreshowerClusterShapeProducer()

PreshowerClusterShapeProducer::~PreshowerClusterShapeProducer ( )

override

Definition at line 65 of file PreshowerClusterShapeProducer.cc.

65 { delete presh_pi0_algo; }

Member Function Documentation

◆ produce()

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

override

Definition at line 67 of file PreshowerClusterShapeProducer.cc.

                                                                             {
   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.getHandle(caloGeometryToken_);
   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();
 }