Inheritance diagram for PreshowerClusterProducer:

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
	PreshowerClusterProducer (const edm::ParameterSet &ps)

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

void	set (const edm::EventSetup &es)

	~PreshowerClusterProducer () 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
double	aEta_ [4]

double	alpha0_

double	alpha1_

double	alpha2_

double	alpha3_

std::string	assocSClusterCollection_

double	bEta_ [4]

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

edm::EDGetTokenT< reco::SuperClusterCollection >	endcapSClusterToken_

edm::ESHandle< ESChannelStatus >	esChannelStatus_

edm::ESGetToken< ESChannelStatus, ESChannelStatusRcd >	esChannelStatusToken_

edm::ESHandle< ESEEIntercalibConstants >	esEEInterCalib_

edm::ESGetToken< ESEEIntercalibConstants, ESEEIntercalibConstantsRcd >	esEEInterCalibToken_

edm::ESHandle< ESGain >	esgain_

edm::ESGetToken< ESGain, ESGainRcd >	esGainToken_

edm::ESHandle< ESMIPToGeVConstant >	esMIPToGeV_

edm::ESGetToken< ESMIPToGeVConstant, ESMIPToGeVConstantRcd >	esMIPToGeVToken_

edm::ESHandle< ESMissingEnergyCalibration >	esMissingECalib_

edm::ESGetToken< ESMissingEnergyCalibration, ESMissingEnergyCalibrationRcd >	esMissingECalibToken_

double	etThresh_

double	gamma0_

double	gamma1_

double	gamma2_

double	gamma3_

double	mip_

int	nEvt_

PreshowerClusterAlgo *	presh_algo

float	preshClustECut

std::string	preshClusterCollectionX_

std::string	preshClusterCollectionY_

edm::EDGetTokenT< EcalRecHitCollection >	preshHitsToken_

int	preshNclust_

Detailed Description

Definition at line 47 of file PreshowerClusterProducer.cc.

Member Typedef Documentation

◆ Point

typedef math::XYZPoint PreshowerClusterProducer::Point

Definition at line 49 of file PreshowerClusterProducer.cc.

Constructor & Destructor Documentation

◆ PreshowerClusterProducer()

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

explicit

Definition at line 111 of file PreshowerClusterProducer.cc.

                                                                             {
   // use configuration file to setup input/output collection names
   preshHitsToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("preshRecHitProducer"));
  
   // Name of a SuperClusterCollection to make associations:
   endcapSClusterToken_ =
       consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapSClusterProducer"));
  
   esGainToken_ = esConsumes<ESGain, ESGainRcd>();
   esMIPToGeVToken_ = esConsumes<ESMIPToGeVConstant, ESMIPToGeVConstantRcd>();
   esEEInterCalibToken_ = esConsumes<ESEEIntercalibConstants, ESEEIntercalibConstantsRcd>();
   esMissingECalibToken_ = esConsumes<ESMissingEnergyCalibration, ESMissingEnergyCalibrationRcd>();
   esChannelStatusToken_ = esConsumes<ESChannelStatus, ESChannelStatusRcd>();
   caloGeometryToken_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
  
   // Output collections:
   preshClusterCollectionX_ = ps.getParameter<std::string>("preshClusterCollectionX");
   preshClusterCollectionY_ = ps.getParameter<std::string>("preshClusterCollectionY");
   preshNclust_ = ps.getParameter<int>("preshNclust");
  
   etThresh_ = ps.getParameter<double>("etThresh");
  
   assocSClusterCollection_ = ps.getParameter<std::string>("assocSClusterCollection");
  
   produces<reco::PreshowerClusterCollection>(preshClusterCollectionX_);
   produces<reco::PreshowerClusterCollection>(preshClusterCollectionY_);
   produces<reco::SuperClusterCollection>(assocSClusterCollection_);
  
   float preshStripECut = ps.getParameter<double>("preshStripEnergyCut");
   int preshSeededNst = ps.getParameter<int>("preshSeededNstrip");
   preshClustECut = ps.getParameter<double>("preshClusterEnergyCut");
  
   presh_algo = new PreshowerClusterAlgo(preshStripECut, preshClustECut, preshSeededNst);
  
   nEvt_ = 0;
 }

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~PreshowerClusterProducer()

PreshowerClusterProducer::~PreshowerClusterProducer ( )

override

Definition at line 148 of file PreshowerClusterProducer.cc.

148 { delete presh_algo; }

Member Function Documentation

◆ produce()

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

override

Definition at line 150 of file PreshowerClusterProducer.cc.

                                                                              {
   edm::Handle<EcalRecHitCollection> pRecHits;
   edm::Handle<reco::SuperClusterCollection> pSuperClusters;
  
   // get the ECAL geometry:
   edm::ESHandle<CaloGeometry> geoHandle = es.getHandle(caloGeometryToken_);
  
   // retrieve ES-EE intercalibration constants and channel status
   set(es);
   const ESChannelStatus* channelStatus = esChannelStatus_.product();
  
   const CaloSubdetectorGeometry* geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
  
   // create unique_ptr to a PreshowerClusterCollection
   auto clusters_p1 = std::make_unique<reco::PreshowerClusterCollection>();
   auto clusters_p2 = std::make_unique<reco::PreshowerClusterCollection>();
   // create new collection of corrected super clusters
   auto superclusters_p = std::make_unique<reco::SuperClusterCollection>();
  
   std::unique_ptr<CaloSubdetectorTopology> topology_p;
   if (geometry_p)
     topology_p = std::make_unique<EcalPreshowerTopology>();
  
   // fetch the product (pSuperClusters)
   evt.getByToken(endcapSClusterToken_, pSuperClusters);
   const reco::SuperClusterCollection* SClusts = pSuperClusters.product();
  
   // fetch the product (RecHits)
   evt.getByToken(preshHitsToken_, pRecHits);
   // pointer to the object in the product
   const EcalRecHitCollection* rechits = pRecHits.product();  // EcalRecHitCollection hit_collection = *rhcHandle;
  
   LogTrace("EcalClusters") << "PreshowerClusterProducerInfo: ### Total # of preshower RecHits: " << rechits->size();
  
   // make the map of rechits:
   std::map<DetId, EcalRecHit> rechits_map;
   EcalRecHitCollection::const_iterator it;
   for (it = rechits->begin(); it != rechits->end(); it++) {
     // remove bad ES rechits
     if (it->recoFlag() == 1 || it->recoFlag() == 14 || (it->recoFlag() <= 10 && it->recoFlag() >= 5))
       continue;
     //Make the map of DetID, EcalRecHit pairs
     rechits_map.insert(std::make_pair(it->id(), *it));
   }
   // The set of used DetID's for a given event:
   std::set<DetId> used_strips;
   used_strips.clear();
   LogTrace("EcalClusters") << "PreshowerClusterProducerInfo: ### rechits_map of size " << rechits_map.size()
                            << " was created!";
  
   reco::PreshowerClusterCollection clusters1, clusters2;  // output collection of corrected PCs
   reco::SuperClusterCollection new_SC;                    // output collection of corrected SCs
  
   //make cycle over super clusters
   reco::SuperClusterCollection::const_iterator it_super;
   int isc = 0;
   for (it_super = SClusts->begin(); it_super != SClusts->end(); ++it_super) {
     float e1 = 0;
     float e2 = 0;
     float deltaE = 0;
  
     reco::CaloClusterPtrVector new_BC;
     ++isc;
     LogTrace("EcalClusters") << " superE = " << it_super->energy() << " superETA = " << it_super->eta()
                              << " superPHI = " << it_super->phi();
  
     int nBC = 0;
     int condP1 = 1;  // 0: dead channel; 1: active channel
     int condP2 = 1;
     reco::CaloCluster_iterator bc_iter = it_super->clustersBegin();
     for (; bc_iter != it_super->clustersEnd(); ++bc_iter) {
       if (geometry_p) {
         // Get strip position at intersection point of the line EE - Vertex:
         double X = (*bc_iter)->x();
         double Y = (*bc_iter)->y();
         double Z = (*bc_iter)->z();
         const GlobalPoint point(X, Y, Z);
  
         DetId tmp1 = (dynamic_cast<const EcalPreshowerGeometry*>(geometry_p))->getClosestCellInPlane(point, 1);
         DetId tmp2 = (dynamic_cast<const EcalPreshowerGeometry*>(geometry_p))->getClosestCellInPlane(point, 2);
         ESDetId strip1 = (tmp1 == DetId(0)) ? ESDetId(0) : ESDetId(tmp1);
         ESDetId strip2 = (tmp2 == DetId(0)) ? ESDetId(0) : ESDetId(tmp2);
  
         if (nBC == 0) {
           if (strip1 != ESDetId(0) && strip2 != ESDetId(0)) {
             ESChannelStatusMap::const_iterator status_p1 = channelStatus->getMap().find(strip1);
             ESChannelStatusMap::const_iterator status_p2 = channelStatus->getMap().find(strip2);
             if (status_p1->getStatusCode() == 1)
               condP1 = 0;
             if (status_p2->getStatusCode() == 1)
               condP2 = 0;
           } else if (strip1 == ESDetId(0))
             condP1 = 0;
           else if (strip2 == ESDetId(0))
             condP2 = 0;
         }
  
         // Get a vector of ES clusters (found by the PreshSeeded algorithm) associated with a given EE basic cluster.
         for (int i = 0; i < preshNclust_; i++) {
           reco::PreshowerCluster cl1 =
               presh_algo->makeOneCluster(strip1, &used_strips, &rechits_map, geometry_p, topology_p.get());
           cl1.setBCRef(*bc_iter);
           if (cl1.energy() > preshClustECut) {
             clusters1.push_back(cl1);
             e1 += cl1.energy();
           }
           reco::PreshowerCluster cl2 =
               presh_algo->makeOneCluster(strip2, &used_strips, &rechits_map, geometry_p, topology_p.get());
           cl2.setBCRef(*bc_iter);
  
           if (cl2.energy() > preshClustECut) {
             clusters2.push_back(cl2);
             e2 += cl2.energy();
           }
         }  // end of cycle over ES clusters
       }
  
       new_BC.push_back(*bc_iter);
       nBC++;
     }  // end of cycle over BCs
  
     LogTrace("EcalClusters") << " For SC #" << isc - 1 << ", containing " << it_super->clustersSize()
                              << " basic clusters, PreshowerClusterAlgo made " << clusters1.size() << " in X plane and "
                              << clusters2.size() << " in Y plane "
                              << " preshower clusters ";
  
     // update energy of the SuperCluster
     if (e1 + e2 > 1.0e-10) {
       e1 = e1 / mip_;  // GeV to #MIPs
       e2 = e2 / mip_;
  
       if (condP1 == 1 && condP2 == 1) {
         deltaE = gamma0_ * (e1 + alpha0_ * e2);
       } else if (condP1 == 1 && condP2 == 0) {
         deltaE = gamma1_ * (e1 + alpha1_ * e2);
       } else if (condP1 == 0 && condP2 == 1) {
         deltaE = gamma2_ * (e1 + alpha2_ * e2);
       } else if (condP1 == 0 && condP2 == 0) {
         deltaE = gamma3_ * (e1 + alpha3_ * e2);
       }
     }
  
     //corrected Energy
     float E = it_super->energy() + deltaE;
  
     LogTrace("EcalClusters") << " Creating corrected SC ";
  
     reco::SuperCluster sc(E, it_super->position(), it_super->seed(), new_BC, deltaE);
     if (condP1 == 1 && condP2 == 1)
       sc.setPreshowerPlanesStatus(0);
     else if (condP1 == 1 && condP2 == 0)
       sc.setPreshowerPlanesStatus(1);
     else if (condP1 == 0 && condP2 == 1)
       sc.setPreshowerPlanesStatus(2);
     else if (condP1 == 0 && condP2 == 0)
       sc.setPreshowerPlanesStatus(3);
  
     if (sc.energy() * sin(sc.position().theta()) > etThresh_)
       new_SC.push_back(sc);
     LogTrace("EcalClusters") << " SuperClusters energies: new E = " << sc.energy()
                              << " vs. old E =" << it_super->energy();
  
   }  // end of cycle over SCs
  
   // copy the preshower clusters into collections and put in the Event:
   clusters_p1->assign(clusters1.begin(), clusters1.end());
   clusters_p2->assign(clusters2.begin(), clusters2.end());
   // put collection of preshower clusters to the event
   evt.put(std::move(clusters_p1), preshClusterCollectionX_);
   evt.put(std::move(clusters_p2), preshClusterCollectionY_);
   LogTrace("EcalClusters") << "Preshower clusters added to the event";
  
   // put collection of corrected super clusters to the event
   superclusters_p->assign(new_SC.begin(), new_SC.end());
   evt.put(std::move(superclusters_p), assocSClusterCollection_);
   LogTrace("EcalClusters") << "Corrected SClusters added to the event";
  
   nEvt_++;
 }

References GetRecoTauVFromDQM_MC_cff::cl2, EgHLTOffHistBins_cfi::deltaE, MillePedeFileConverter_cfg::e, StorageManager_cfg::e1, DetId::Ecal, EcalPreshower, reco::CaloCluster::energy(), ESCondObjectContainer< T >::find(), edm::Event::getByToken(), edm::EventSetup::getHandle(), ESCondObjectContainer< T >::getMap(), CaloGeometry::getSubdetectorGeometry(), mps_fire::i, LogTrace, eostools::move(), point, edm::Handle< T >::product(), edm::PtrVector< T >::push_back(), edm::Event::put(), HI_PhotonSkim_cff::rechits, reco::PreshowerCluster::setBCRef(), reco::SuperCluster::setPreshowerPlanesStatus(), funct::sin(), X, BeamSpotPI::Y, and BeamSpotPI::Z.

◆ set()

void PreshowerClusterProducer::set ( const edm::EventSetup & es )

Definition at line 330 of file PreshowerClusterProducer.cc.

                                                           {
   esgain_ = es.getHandle(esGainToken_);
   const ESGain* gain = esgain_.product();
  
   double ESGain = gain->getESGain();
  
   esMIPToGeV_ = es.getHandle(esMIPToGeVToken_);
   const ESMIPToGeVConstant* mipToGeV = esMIPToGeV_.product();
  
   mip_ = (ESGain == 1) ? mipToGeV->getESValueLow() : mipToGeV->getESValueHigh();
  
   esChannelStatus_ = es.getHandle(esChannelStatusToken_);
  
   esEEInterCalib_ = es.getHandle(esEEInterCalibToken_);
   const ESEEIntercalibConstants* esEEInterCalib = esEEInterCalib_.product();
  
   // both planes work
   gamma0_ = (ESGain == 1) ? 0.02 : esEEInterCalib->getGammaHigh0();
   alpha0_ = (ESGain == 1) ? esEEInterCalib->getAlphaLow0() : esEEInterCalib->getAlphaHigh0();
  
   // only first plane works
   gamma1_ = (ESGain == 1) ? (0.02 * esEEInterCalib->getGammaLow1()) : esEEInterCalib->getGammaHigh1();
   alpha1_ = (ESGain == 1) ? esEEInterCalib->getAlphaLow1() : esEEInterCalib->getAlphaHigh1();
  
   // only second plane works
   gamma2_ = (ESGain == 1) ? (0.02 * esEEInterCalib->getGammaLow2()) : esEEInterCalib->getGammaHigh2();
   alpha2_ = (ESGain == 1) ? esEEInterCalib->getAlphaLow2() : esEEInterCalib->getAlphaHigh2();
  
   // both planes do not work
   gamma3_ = (ESGain == 1) ? 0.02 : esEEInterCalib->getGammaHigh3();
   alpha3_ = (ESGain == 1) ? esEEInterCalib->getAlphaLow3() : esEEInterCalib->getAlphaHigh3();
  
   esMissingECalib_ = es.getHandle(esMissingECalibToken_);
   const ESMissingEnergyCalibration* esMissingECalib = esMissingECalib_.product();
  
   // |eta| < 1.9
   aEta_[0] = esMissingECalib->getConstAEta0();
   bEta_[0] = esMissingECalib->getConstBEta0();
  
   // 1.9 < |eta| < 2.1
   aEta_[1] = esMissingECalib->getConstAEta1();
   bEta_[1] = esMissingECalib->getConstBEta1();
  
   // 2.1 < |eta| < 2.3
   aEta_[2] = esMissingECalib->getConstAEta2();
   bEta_[2] = esMissingECalib->getConstBEta2();
  
   // 2.3 < |eta| < 2.5
   aEta_[3] = esMissingECalib->getConstAEta3();
   bEta_[3] = esMissingECalib->getConstBEta3();
 }