EgammaHLTPixelMatchVarProducer Class Reference

Inheritance diagram for EgammaHLTPixelMatchVarProducer:

Public Member Functions
std::array< float, 4 >	calS2 (const reco::ElectronSeed &seed, int charge) const
	EgammaHLTPixelMatchVarProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~EgammaHLTPixelMatchVarProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
egPM::Param< reco::ElectronSeed >	dPhi1Para_
egPM::Param< reco::ElectronSeed >	dPhi2Para_
egPM::Param< reco::ElectronSeed >	dRZ2Para_
size_t	nrHits_
const edm::EDGetTokenT< reco::ElectronSeedCollection >	pixelSeedsToken_
int	productsToWrite_
const edm::EDGetTokenT< reco::RecoEcalCandidateCollection >	recoEcalCandidateToken_

Additional Inherited Members
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

Detailed Description

Definition at line 101 of file EgammaHLTPixelMatchVarProducer.cc.

Constructor & Destructor Documentation

EgammaHLTPixelMatchVarProducer::EgammaHLTPixelMatchVarProducer ( const edm::ParameterSet &  config )

explicit

Definition at line 125 of file EgammaHLTPixelMatchVarProducer.cc.

References nrHits_, and productsToWrite_.

                                                                                            : 
  recoEcalCandidateToken_(consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"))),
  pixelSeedsToken_(consumes<reco::ElectronSeedCollection>(config.getParameter<edm::InputTag>("pixelSeedsProducer"))),
  dPhi1Para_(config.getParameter<edm::ParameterSet>("dPhi1SParams")),
  dPhi2Para_(config.getParameter<edm::ParameterSet>("dPhi2SParams")),
  dRZ2Para_(config.getParameter<edm::ParameterSet>("dRZ2SParams")),
  productsToWrite_(config.getParameter<int>("productsToWrite")),
  nrHits_(4)
  
{
  //register your products  
  produces < reco::RecoEcalCandidateIsolationMap >("s2");
  if(productsToWrite_>=1){
    produces < reco::RecoEcalCandidateIsolationMap >("dPhi1BestS2");
    produces < reco::RecoEcalCandidateIsolationMap >("dPhi2BestS2");
    produces < reco::RecoEcalCandidateIsolationMap >("dzBestS2");
  }
  if(productsToWrite_>=2){
    //note for product names we start from index 1
    for(size_t hitNr=1;hitNr<=nrHits_;hitNr++){
      produces < reco::RecoEcalCandidateIsolationMap >("dPhi"+std::to_string(hitNr));
      produces < reco::RecoEcalCandidateIsolationMap >("dPhi"+std::to_string(hitNr)+"Info");
      produces < reco::RecoEcalCandidateIsolationMap >("dRZ"+std::to_string(hitNr));
      produces < reco::RecoEcalCandidateIsolationMap >("dRZ"+std::to_string(hitNr)+"Info");
    } 
    produces < reco::RecoEcalCandidateIsolationMap >("nrClus");
    produces < reco::RecoEcalCandidateIsolationMap >("seedClusEFrac");
    produces < reco::RecoEcalCandidateIsolationMap >("phiWidth");
    produces < reco::RecoEcalCandidateIsolationMap >("etaWidth");
  }

}

EgammaHLTPixelMatchVarProducer::~EgammaHLTPixelMatchVarProducer ( )

override

Definition at line 158 of file EgammaHLTPixelMatchVarProducer.cc.

{}

Member Function Documentation

std::array< float, 4 > EgammaHLTPixelMatchVarProducer::calS2	(	const reco::ElectronSeed &	seed,
		int	charge
	)		const

Definition at line 309 of file EgammaHLTPixelMatchVarProducer.cc.

References DEFINE_FWK_MODULE, dPhi1Para_, dPhi2Para_, reco::ElectronSeed::dPhiNeg(), reco::ElectronSeed::dPhiPos(), dRZ2Para_, reco::ElectronSeed::dRZNeg(), reco::ElectronSeed::dRZPos(), and indexGen::s2.

Referenced by produce().

{
  const float dPhi1Const = dPhi1Para_(seed);    
  const float dPhi2Const = dPhi2Para_(seed);
  const float dRZ2Const = dRZ2Para_(seed);
  
  float dPhi1 = (charge <0 ? seed.dPhiNeg(0) : seed.dPhiPos(0))/dPhi1Const;
  float dPhi2 = (charge <0 ? seed.dPhiNeg(1) : seed.dPhiPos(1))/dPhi2Const;
  float dRz2 = (charge <0 ? seed.dRZNeg(1) : seed.dRZPos(1))/dRZ2Const;
  
  float s2 = dPhi1*dPhi1+dPhi2*dPhi2+dRz2*dRz2;
  return std::array<float,4>{{s2,dPhi1,dPhi2,dRz2}}; 
}

void EgammaHLTPixelMatchVarProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 161 of file EgammaHLTPixelMatchVarProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::ifValue(), eostools::move(), and or.

                                                                                                {

  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>(("recoEcalCandidateProducer"), edm::InputTag("hltL1SeededRecoEcalCandidate"));
  desc.add<edm::InputTag>(("pixelSeedsProducer"), edm::InputTag("electronPixelSeeds"));
  
  edm::ParameterSetDescription varParamDesc;
  edm::ParameterSetDescription binParamDesc;
  
  auto binDescCases = 
    "AbsEtaClus" >> 
    (edm::ParameterDescription<double>("xMin",0.0,true) and
     edm::ParameterDescription<double>("xMax",3.0,true) and
     edm::ParameterDescription<int>("yMin",0,true) and
     edm::ParameterDescription<int>("yMax",99999,true) and
     edm::ParameterDescription<std::string>("funcType","pol0",true) and
     edm::ParameterDescription<std::vector<double>>("funcParams",{0.},true)) or
    "AbsEtaClusPhi" >>
    (edm::ParameterDescription<double>("xMin",0.0,true) and
     edm::ParameterDescription<double>("xMax",3.0,true) and
     edm::ParameterDescription<int>("yMin",0,true) and
     edm::ParameterDescription<int>("yMax",99999,true) and
     edm::ParameterDescription<std::string>("funcType","pol0",true) and
     edm::ParameterDescription<std::vector<double>>("funcParams",{0.},true)) or 
     "AbsEtaClusEt" >>
    (edm::ParameterDescription<double>("xMin",0.0,true) and
     edm::ParameterDescription<double>("xMax",3.0,true) and
     edm::ParameterDescription<int>("yMin",0,true) and
     edm::ParameterDescription<int>("yMax",99999,true) and
     edm::ParameterDescription<std::string>("funcType","pol0",true) and
     edm::ParameterDescription<std::vector<double>>("funcParams",{0.},true));
  
  binParamDesc.ifValue(edm::ParameterDescription<std::string>("binType","AbsEtaClus",true), std::move(binDescCases));
  
  
  varParamDesc.addVPSet("bins",binParamDesc);
  desc.add("dPhi1SParams",varParamDesc);
  desc.add("dPhi2SParams",varParamDesc);
  desc.add("dRZ2SParams",varParamDesc);
  desc.add<int>("productsToWrite",0);
  descriptions.add(("hltEgammaHLTPixelMatchVarProducer"), desc);  
}

void EgammaHLTPixelMatchVarProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 204 of file EgammaHLTPixelMatchVarProducer.cc.

References reco::ElectronSeed::caloCluster(), calS2(), edm::RefToBase< T >::castTo(), reco::ElectronSeed::dPhiBest(), reco::ElectronSeed::dRZBest(), edm::EventSetup::get(), edm::Event::getByToken(), edm::HandleBase::isValid(), SiStripPI::max, eostools::move(), nrHits_, pixelSeedsToken_, productsToWrite_, edm::Event::put(), and recoEcalCandidateToken_.

                                                                                         {
  
  // Get the HLT filtered objects
  edm::Handle<reco::RecoEcalCandidateCollection> recoEcalCandHandle;
  iEvent.getByToken(recoEcalCandidateToken_,recoEcalCandHandle);


  edm::Handle<reco::ElectronSeedCollection> pixelSeedsHandle;
  iEvent.getByToken(pixelSeedsToken_,pixelSeedsHandle);

  if(!recoEcalCandHandle.isValid()) return;
  else if(!pixelSeedsHandle.isValid()){
    auto s2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
    for(unsigned int candNr = 0; candNr<recoEcalCandHandle->size(); candNr++) {
      reco::RecoEcalCandidateRef candRef(recoEcalCandHandle,candNr);
      s2Map->insert(candRef,0);
    }
    iEvent.put(std::move(s2Map),"s2");
    return;
  }

  edm::ESHandle<TrackerTopology> trackerTopoHandle;
  iSetup.get<TrackerTopologyRcd>().get(trackerTopoHandle);
  
  auto dPhi1BestS2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
  auto dPhi2BestS2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
  auto dzBestS2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
  auto s2Map = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
  
  auto nrClusMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle);
  auto seedClusEFracMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle); 
  auto phiWidthMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle); 
  auto etaWidthMap = std::make_unique<reco::RecoEcalCandidateIsolationMap>(recoEcalCandHandle); 
  
  std::vector<PixelData> pixelData;
  for(size_t hitNr=0;hitNr<nrHits_;hitNr++){
    pixelData.emplace_back(PixelData("dPhi",hitNr,&reco::ElectronSeed::dPhiBest,recoEcalCandHandle));
    pixelData.emplace_back(PixelData("dRZ",hitNr,&reco::ElectronSeed::dRZBest,recoEcalCandHandle));
  }

  for(unsigned int candNr = 0; candNr<recoEcalCandHandle->size(); candNr++) {

    reco::RecoEcalCandidateRef candRef(recoEcalCandHandle,candNr);
    reco::SuperClusterRef candSCRef = candRef->superCluster();
    
    std::array<float,4> bestS2{{std::numeric_limits<float>::max(),std::numeric_limits<float>::max(),std::numeric_limits<float>::max(),std::numeric_limits<float>::max()}};
    for(auto & seed : *pixelSeedsHandle){
      const edm::RefToBase<reco::CaloCluster>& pixelClusterRef = seed.caloCluster() ;
      reco::SuperClusterRef pixelSCRef = pixelClusterRef.castTo<reco::SuperClusterRef>() ;
      if(&(*candSCRef) ==  &(*pixelSCRef)){
    
    std::array<float,4> s2Data = calS2(seed,-1);
    std::array<float,4> s2DataPos = calS2(seed,+1);
    if(s2Data[0]<bestS2[0]) bestS2=s2Data;
    if(s2DataPos[0]<bestS2[0]) bestS2=s2DataPos;
    
    if(productsToWrite_>=2){
      for(auto& pixelDatum : pixelData){
        pixelDatum.fill(seed);
      }
    }
      }
    }

   
    s2Map->insert(candRef,bestS2[0]);
    if(productsToWrite_>=1){
      dPhi1BestS2Map->insert(candRef,bestS2[1]);
      dPhi2BestS2Map->insert(candRef,bestS2[2]);
      dzBestS2Map->insert(candRef,bestS2[3]);
    }
    if(productsToWrite_>=2){
       nrClusMap->insert(candRef,candSCRef->clustersSize());
       float seedClusEFrac = candSCRef->rawEnergy()>0 ? candSCRef->seed()->energy() / candSCRef->rawEnergy() : 0.;
       //       std::cout <<"cand "<<candSCRef->energy()<<" E Corr "<<candSCRef->correctedEnergyUncertainty()<<" "<<candSCRef->correctedEnergy()<<" width "<<candSCRef->phiWidth()<<std::endl;
       //  float seedClusEFrac = candSCRef->phiWidth();
       seedClusEFracMap->insert(candRef,seedClusEFrac);
       float phiWidth = candSCRef->phiWidth();
       float etaWidth = candSCRef->etaWidth();
       phiWidthMap->insert(candRef,phiWidth);
       etaWidthMap->insert(candRef,etaWidth);
       
       for(auto& pixelDatum : pixelData){
     pixelDatum.fill(candRef);
       }
    } 
  }

  iEvent.put(std::move(s2Map),"s2");
  if(productsToWrite_>=1){
    iEvent.put(std::move(dPhi1BestS2Map),"dPhi1BestS2");
    iEvent.put(std::move(dPhi2BestS2Map),"dPhi2BestS2");
    iEvent.put(std::move(dzBestS2Map),"dzBestS2");
  }
  if(productsToWrite_>=2){
    for(auto& pixelDatum : pixelData){
      pixelDatum.putInto(iEvent);
    }   
    iEvent.put(std::move(nrClusMap),"nrClus");
    iEvent.put(std::move(seedClusEFracMap),"seedClusEFrac");
    iEvent.put(std::move(phiWidthMap),"phiWidth");
    iEvent.put(std::move(etaWidthMap),"etaWidth");
  }
}

Member Data Documentation

egPM::Param<reco::ElectronSeed> EgammaHLTPixelMatchVarProducer::dPhi1Para_

private

Definition at line 117 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

egPM::Param<reco::ElectronSeed> EgammaHLTPixelMatchVarProducer::dPhi2Para_

private

Definition at line 118 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

egPM::Param<reco::ElectronSeed> EgammaHLTPixelMatchVarProducer::dRZ2Para_

private

Definition at line 119 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by calS2().

size_t EgammaHLTPixelMatchVarProducer::nrHits_

private

Definition at line 122 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by EgammaHLTPixelMatchVarProducer(), and produce().

const edm::EDGetTokenT<reco::ElectronSeedCollection> EgammaHLTPixelMatchVarProducer::pixelSeedsToken_

private

Definition at line 115 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by produce().

int EgammaHLTPixelMatchVarProducer::productsToWrite_

private

Definition at line 121 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by EgammaHLTPixelMatchVarProducer(), and produce().

const edm::EDGetTokenT<reco::RecoEcalCandidateCollection> EgammaHLTPixelMatchVarProducer::recoEcalCandidateToken_

private

Definition at line 114 of file EgammaHLTPixelMatchVarProducer.cc.

Referenced by produce().