26 using namespace egHLT;
156 if(errCode==0) errCode =
getHandles(edmEvent,setup);
159 if(errCode==0) errCode =
setTrigInfo(edmEvent, offEvent, c);
204 egHLTOffEles.clear();
206 for(reco::GsfElectronCollection::const_iterator gsfIter=
recoEles_->begin(); gsfIter!=
recoEles_->end();++gsfIter){
207 if(!gsfIter->ecalDrivenSeed())
continue;
213 && std::fabs( nVit->z()<24.0)
214 &&
sqrt(nVit->x()*nVit->x() + nVit->y()*nVit->y())<2.0){nVertex++;}
229 egHLTOffEles.push_back(
OffEle(*gsfIter,clusShapeData,isolData,hltData,eventData));
232 OffEle& ele = egHLTOffEles.back();
236 std::vector<std::pair<TrigCodes::TrigBitSet,int> >trigCutsCutCodes;
259 else isolData.
hltHad = 0.;
268 ecalIsolAlgoEE.getEtSum(&ele);
269 else isolData.
hltEm = 0.;
278 double e5x5 = ele.
e5x5();
289 const DetId seedDetId = seedClus.hitsAndFractions()[0].first;
297 }
else clusShapeData.
r9 = -1.;
306 }
else clusShapeData.
r9 = -1.;
320 math::XYZPoint scPosWRTVtx(scClus->x()-bsPos.x(), scClus->y()-bsPos.y() , scClus->z()-ctfTrack->vz());
321 hltData.
dEtaIn = fabs(scPosWRTVtx.eta()-ctfTrack->eta());
324 float deltaPhi=fabs(ctfTrack->outerPosition().phi()-scClus->phi());
325 if(deltaPhi>6.283185308) deltaPhi -= 6.283185308;
326 if(deltaPhi>3.141592654) deltaPhi = 6.283185308-
deltaPhi;
358 egHLTOffPhos.clear();
360 for(reco::PhotonCollection::const_iterator phoIter=
recoPhos_->begin(); phoIter!=
recoPhos_->end();++phoIter){
371 egHLTOffPhos.push_back(
OffPho(*phoIter,clusShapeData,isolData,hltData));
372 OffPho& pho = egHLTOffPhos.back();
376 std::vector<std::pair<TrigCodes::TrigBitSet,int> >trigCutsCutCodes;
401 else isolData.
hltHad = 0.;
408 ecalIsolAlgoEE.getEtSum(&pho);
409 else isolData.
hltEm = 0.;
417 double e5x5 = pho.
e5x5();
425 clusShapeData.
r9 = pho.
r9();
431 const DetId seedDetId = seedClus.hitsAndFractions()[0].first;
void setTrigCutsCutCodes(const std::vector< std::pair< TrigCodes::TrigBitSet, int > > &trigCutsCutCodes)
T getParameter(std::string const &) const
double hltEMIsolEtaSliceEB_
edm::Handle< reco::TrackCollection > isolTrks_
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
bool isNonnull() const
Checks for non-null.
The single EDProduct to be saved for each event (AOD case)
reco::SuperClusterRef superCluster() const
Ref to SuperCluster.
int fillOffPhoVec(std::vector< OffPho > &offPhos)
double hltPhoTrkIsolRSpan_
TrigBitSet getCode(const char *descript) const
edm::EDGetTokenT< CaloTowerCollection > caloTowersToken
float hcalTowerSumEtConeDR03() const
Hcal isolation sum.
void setLooseCutCode(int code)
bool hltPhoTrkIsolCountTrks_
void fillClusShapeData(const reco::GsfElectron &ele, OffEle::ClusShapeData &clusShapeData)
edm::Handle< reco::PhotonCollection > recoPhos_
std::vector< std::pair< TrigCodes::TrigBitSet, OffEgSel > > trigCuts_
int getCutCode(const OffEle &ele, int cutMask=~0x0) const
int setTrigInfo(const edm::Event &edmEvent, egHLT::OffEvt &offEvent, const TrigCodes &trigCodes)
double hltEMIsolOuterCone_
std::vector< Track > TrackCollection
collection of Tracks
float e1x5() const
Shower shape variables.
edm::EDGetTokenT< HBHERecHitCollection > hbheHitsToken
std::vector< std::string > l1PreScaledPaths_
EgammaHLTTrackIsolation * hltPhoTrkIsolAlgo_
std::vector< Vertex > VertexCollection
collection of Vertex objects
std::vector< std::pair< std::string, std::string > > l1PreAndSeedFilters_
double hltEleTrkIsolOuterCone_
math::XYZPoint Point
point in the space
void fillIsolData(const reco::GsfElectron &ele, OffEle::IsolData &isolData)
int photonTrackCount(const reco::RecoCandidate *const recocand, const reco::TrackCollection *isoTracks, bool useVertex)
float dr03HcalDepth2TowerSumEt() const
void fillHLTData(const reco::GsfElectron &ele, OffEle::HLTData &hltData)
float trkSumPtHollowConeDR03() const
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
edm::EDGetTokenT< reco::TrackCollection > isolTrkToken
double hltEMIsolInnerConeEE_
std::vector< GsfElectron > GsfElectronCollection
collection of GsfElectron objects
void setEvtTrigBits(TrigCodes::TrigBitSet bits)
edm::EDGetTokenT< reco::GsfElectronCollection > electronsToken
std::vector< edm::ParameterSet > trigCutParams_
double hltEleTrkIsolInnerCone_
TrackRef closestCtfTrackRef() const
edm::Handle< HBHERecHitCollection > hbheHits_
double hltPhoTrkIsolPtMin_
double hltEleTrkIsolZSpan_
edm::EDGetTokenT< trigger::TriggerEvent > triggerSummaryToken
void fillHLTDataPho(const reco::Photon &pho, OffPho::HLTData &hltData)
float sigmaIetaIeta() const
unsigned int triggerIndex(std::string const &name) const
double hltHadIsolInnerCone_
float electronPtSum(const reco::Track *const tr, const reco::TrackCollection *isoTracks)
Get Pt sum of tracks inside an isolation cone for electrons.
void setup(const edm::ParameterSet &)
virtual SuperClusterRef superCluster() const
reference to a SuperCluster
double hltPhoTrkIsolZSpan_
std::vector< std::pair< std::string, int > > hltFiltersUsedWithNrCandsCut_
edm::EDGetTokenT< EcalRecHitCollection > ecalRecHitsEEToken
double hltEleTrkIsolRSpan_
float sigmaIetaIeta() const
double hltHadIsolOuterCone_
edm::EDGetTokenT< reco::PhotonCollection > photonsToken
float dr03TkSumPt() const
static bool getHandle(const edm::Event &event, const edm::EDGetTokenT< T > &token, edm::Handle< T > &handle)
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken
edm::Handle< HFRecHitCollection > hfHits_
float dr03HcalDepth1TowerSumEt() const
void setLooseCutCode(int code)
edm::Handle< std::vector< reco::CaloJet > > recoJets_
float sigmaEtaEta() const
double hltEleTrkIsolPtMin_
int nTrkHollowConeDR03() const
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
void setCutCode(int code)
edm::ESHandle< MagneticField > magField_
edm::Handle< trigger::TriggerEvent > trigEvt_
edm::ESHandle< CaloGeometry > caloGeom_
const std::vector< OffPho > & phos() const
double hltPhoTrkIsolInnerCone_
void setCutCode(int code)
double hltEMIsolEtaSliceEE_
double hltEMIsolInnerConeEB_
edm::Handle< reco::GsfElectronCollection > recoEles_
edm::EDGetTokenT< edm::TriggerResults > trigResultsToken
T const * product() const
edm::EDGetTokenT< HFRecHitCollection > hfHitsToken
edm::Handle< edm::TriggerResults > trigResults_
edm::ESHandle< CaloTopology > caloTopology_
XYZPointD XYZPoint
point in space with cartesian internal representation
int makeOffEvt(const edm::Event &edmEvent, const edm::EventSetup &setup, egHLT::OffEvt &offEvent, const TrigCodes &trigCodes)
std::vector< Photon > PhotonCollection
collectin of Photon objects
T const * product() const
edm::EDGetTokenT< reco::CaloJetCollection > caloJetsToken
std::vector< std::string > l1PreScaledFilters_
void setJets(edm::Handle< std::vector< reco::CaloJet > > jets)
float dr03EcalRecHitSumEt() const
edm::EDGetTokenT< EcalRecHitCollection > ecalRecHitsEBToken
edm::Handle< EcalRecHitCollection > eeRecHits_
const std::vector< OffEle > & eles() const
void setTrigCutsCutCodes(const std::vector< std::pair< TrigCodes::TrigBitSet, int > > &trigCutsCutCodes)
int getHandles(const edm::Event &event, const edm::EventSetup &setup)
void setup(const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)
double hltPhoTrkIsolOuterCone_
float ecalRecHitSumEtConeDR03() const
Isolation variables in cone dR=0.3.
edm::Handle< CaloTowerCollection > caloTowers_
edm::Handle< EcalRecHitCollection > ebRecHits_
std::vector< std::string > hltFiltersUsed_
edm::Handle< reco::VertexCollection > recoVertices_
edm::ESHandle< EcalSeverityLevelAlgo > ecalSeverityLevel_
void setupTriggers(const HLTConfigProvider &config, const std::vector< std::string > &hltFiltersUsed, const TrigCodes &trigCodes)
EgammaHLTTrackIsolation * hltEleTrkIsolAlgo_
edm::Handle< reco::BeamSpot > beamSpot_
void setup(std::vector< TH2F > &depth, std::string name, std::string units="")
tuple size
Write out results.
float photonPtSum(const reco::RecoCandidate *const recocand, const reco::TrackCollection *isoTracks, bool useVertex)
std::vector< CaloJet > CaloJetCollection
collection of CaloJet objects
float sigmaEtaEta() const
edm::EDGetTokenT< reco::VertexCollection > vertexToken
int fillOffEleVec(std::vector< OffEle > &offEles)
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
std::bitset< maxNrBits_ > TrigBitSet