Inheritance diagram for PreshowerClusterShapeProducer:

Public Types
typedef math::XYZPoint	Point

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

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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 33 of file PreshowerClusterShapeProducer.cc.

Member Typedef Documentation

◆ Point

typedef math::XYZPoint PreshowerClusterShapeProducer::Point

Definition at line 35 of file PreshowerClusterShapeProducer.cc.

Constructor & Destructor Documentation

◆ PreshowerClusterShapeProducer()

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

explicit

Definition at line 68 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 96 of file PreshowerClusterShapeProducer.cc.

96 { delete presh_pi0_algo; }

Member Function Documentation

◆ produce()

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

override

Definition at line 98 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();
 }