CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
HLTTauRefProducer Class Reference

#include <HLTTauRefProducer.h>

Inheritance diagram for HLTTauRefProducer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 HLTTauRefProducer (const edm::ParameterSet &)
 
virtual void produce (edm::Event &, const edm::EventSetup &)
 
 ~HLTTauRefProducer ()
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducer ()
 
- Public Member Functions inherited from edm::ProducerBase
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription
const &)> 
registrationCallback () const
 used by the fwk to register list of products More...
 
virtual ~ProducerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 
itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Types

typedef math::XYZTLorentzVectorD LorentzVector
 
typedef std::vector
< LorentzVector
LorentzVectorCollection
 

Private Member Functions

void doElectrons (edm::Event &, const edm::EventSetup &)
 
void doElectronsFromZ (edm::Event &, const edm::EventSetup &, std::auto_ptr< LorentzVectorCollection > &)
 
void doJets (edm::Event &, const edm::EventSetup &)
 
void doMuons (edm::Event &, const edm::EventSetup &)
 
void doPFTaus (edm::Event &, const edm::EventSetup &)
 
void doPhotons (edm::Event &, const edm::EventSetup &)
 
void doTowers (edm::Event &, const edm::EventSetup &)
 
double ElectronTrkIsolation (const reco::TrackCollection *, const reco::GsfElectron &)
 

Private Attributes

bool doElecFromZ_
 
bool doElectrons_
 
bool doJets_
 
bool doMuons_
 
bool doPFTaus_
 
bool doPhotons_
 
bool doTowers_
 
edm::EDGetTokenT
< reco::TrackCollection
e_ctfTrackCollection_
 
edm::InputTag e_ctfTrackCollectionSrc_
 
bool e_doID_
 
bool e_doTrackIso_
 
double e_FromZet_
 
edm::EDGetTokenT
< reco::ElectronIDAssociationCollection
e_idAssocProd_
 
double e_isoMaxSumPt_
 
double e_lipCut_
 
double e_maxIsoDR_
 
double e_minIsoDR_
 
double e_trackMinPt_
 
double e_zMmax_
 
double e_zMmin_
 
edm::EDGetTokenT
< reco::GsfElectronCollection
Electrons_
 
double etaMax
 
edm::EDGetTokenT
< reco::CaloJetCollection
Jets_
 
edm::EDGetTokenT
< reco::MuonCollection
Muons_
 
std::vector< edm::EDGetTokenT
< reco::PFTauDiscriminator > > 
PFTauDis_
 
edm::EDGetTokenT
< reco::PFTauCollection
PFTaus_
 
double photonEcalIso_
 
edm::EDGetTokenT
< reco::PhotonCollection
Photons_
 
double ptMinElectron_
 
double ptMinJet_
 
double ptMinMuon_
 
double ptMinPFTau_
 
double ptMinPhoton_
 
double ptMinTower_
 
double towerIsol_
 
edm::EDGetTokenT
< CaloTowerCollection
Towers_
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef
ProductRegistryHelper::TypeLabelList 
TypeLabelList
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 48 of file HLTTauRefProducer.h.

Member Typedef Documentation

Definition at line 57 of file HLTTauRefProducer.h.

Definition at line 58 of file HLTTauRefProducer.h.

Constructor & Destructor Documentation

HLTTauRefProducer::HLTTauRefProducer ( const edm::ParameterSet iConfig)
explicit

Definition at line 33 of file HLTTauRefProducer.cc.

References HI_PhotonSkim_cff::electrons, HLT_25ns14e33_v1_cff::etaMax, edm::ParameterSet::getUntrackedParameter(), fwrapper::jets, patZpeak::muons, interactiveExample::photons, GlobalPosition_Frontier_DevDB_cff::tag, and HLT_25ns14e33_v1_cff::towers.

34 {
35 
36 
37  //One Parameter Set per Collection
38 
39  ParameterSet pfTau = iConfig.getUntrackedParameter<edm::ParameterSet>("PFTaus");
40  PFTaus_ = consumes<reco::PFTauCollection>(pfTau.getUntrackedParameter<InputTag>("PFTauProducer"));
41  auto discs = pfTau.getUntrackedParameter<std::vector<InputTag> >("PFTauDiscriminators");
42  for(edm::InputTag& tag: discs) {
43  PFTauDis_.push_back(consumes<reco::PFTauDiscriminator>(tag));
44  }
45  doPFTaus_ = pfTau.getUntrackedParameter<bool>("doPFTaus",false);
46  ptMinPFTau_= pfTau.getUntrackedParameter<double>("ptMin",15.);
47 
49  Electrons_ = consumes<reco::GsfElectronCollection>(electrons.getUntrackedParameter<InputTag>("ElectronCollection"));
50  doElectrons_ = electrons.getUntrackedParameter<bool>("doElectrons",false);
51  e_doID_ = electrons.getUntrackedParameter<bool>("doID",false);
52  if(e_doID_) {
53  e_idAssocProd_ = consumes<reco::ElectronIDAssociationCollection>(electrons.getUntrackedParameter<InputTag>("IdCollection"));
54  }
55  e_ctfTrackCollectionSrc_ = electrons.getUntrackedParameter<InputTag>("TrackCollection");
56  e_ctfTrackCollection_ = consumes<reco::TrackCollection>(e_ctfTrackCollectionSrc_);
57  ptMinElectron_= electrons.getUntrackedParameter<double>("ptMin",15.);
58  e_doTrackIso_ = electrons.getUntrackedParameter<bool>("doTrackIso",false);
59  e_trackMinPt_= electrons.getUntrackedParameter<double>("ptMinTrack",1.5);
60  e_lipCut_= electrons.getUntrackedParameter<double>("lipMinTrack",1.5);
61  e_minIsoDR_= electrons.getUntrackedParameter<double>("InnerConeDR",0.02);
62  e_maxIsoDR_= electrons.getUntrackedParameter<double>("OuterConeDR",0.6);
63  e_isoMaxSumPt_= electrons.getUntrackedParameter<double>("MaxIsoVar",0.02);
64 
66  Muons_ = consumes<reco::MuonCollection>(muons.getUntrackedParameter<InputTag>("MuonCollection"));
67  doMuons_ = muons.getUntrackedParameter<bool>("doMuons",false);
68  ptMinMuon_= muons.getUntrackedParameter<double>("ptMin",15.);
69 
71  Jets_ = consumes<reco::CaloJetCollection>(jets.getUntrackedParameter<InputTag>("JetCollection"));
72  doJets_ = jets.getUntrackedParameter<bool>("doJets");
73  ptMinJet_= jets.getUntrackedParameter<double>("etMin");
74 
76  Towers_ = consumes<CaloTowerCollection>(towers.getUntrackedParameter<InputTag>("TowerCollection"));
77  doTowers_ = towers.getUntrackedParameter<bool>("doTowers");
78  ptMinTower_= towers.getUntrackedParameter<double>("etMin");
79  towerIsol_= towers.getUntrackedParameter<double>("towerIsolation");
80 
82  Photons_ = consumes<reco::PhotonCollection>(photons.getUntrackedParameter<InputTag>("PhotonCollection"));
83  doPhotons_ = photons.getUntrackedParameter<bool>("doPhotons");
84  ptMinPhoton_= photons.getUntrackedParameter<double>("etMin");
85  photonEcalIso_= photons.getUntrackedParameter<double>("ECALIso");
86 
87 
88  etaMax = iConfig.getUntrackedParameter<double>("EtaMax",2.5);
89 
90 
91  //recoCollections
92  produces<LorentzVectorCollection>("PFTaus");
93  produces<LorentzVectorCollection>("Electrons");
94  produces<LorentzVectorCollection>("Muons");
95  produces<LorentzVectorCollection>("Jets");
96  produces<LorentzVectorCollection>("Photons");
97  produces<LorentzVectorCollection>("Towers");
98 
99 }
T getUntrackedParameter(std::string const &, T const &) const
edm::EDGetTokenT< reco::MuonCollection > Muons_
edm::EDGetTokenT< reco::GsfElectronCollection > Electrons_
edm::EDGetTokenT< reco::TrackCollection > e_ctfTrackCollection_
vector< PseudoJet > jets
edm::EDGetTokenT< CaloTowerCollection > Towers_
edm::EDGetTokenT< reco::PhotonCollection > Photons_
edm::EDGetTokenT< reco::PFTauCollection > PFTaus_
edm::EDGetTokenT< reco::ElectronIDAssociationCollection > e_idAssocProd_
std::vector< edm::EDGetTokenT< reco::PFTauDiscriminator > > PFTauDis_
edm::InputTag e_ctfTrackCollectionSrc_
tuple muons
Definition: patZpeak.py:38
edm::EDGetTokenT< reco::CaloJetCollection > Jets_
HLTTauRefProducer::~HLTTauRefProducer ( )

Definition at line 101 of file HLTTauRefProducer.cc.

101 { }

Member Function Documentation

void HLTTauRefProducer::doElectrons ( edm::Event iEvent,
const edm::EventSetup iES 
)
private

Definition at line 153 of file HLTTauRefProducer.cc.

References edm::AssociationMap< edm::OneToOne< GsfElectronCollection, ElectronIDCollection > >::const_iterator, HI_PhotonSkim_cff::electrons, relval_parameters_module::energy, eta(), HLT_25ns14e33_v1_cff::etaMax, edm::Event::getByToken(), i, xuti::id_tag(), edm::HandleBase::isValid(), HLT_25ns14e33_v1_cff::lip, edm::Handle< T >::product(), EnergyCorrector::pt, edm::Event::put(), and mathSSE::sqrt().

154 {
155  auto_ptr<LorentzVectorCollection> product_Electrons(new LorentzVectorCollection);
156  //Retrieve the collections
157 
159  if(e_doID_){//UGLY HACK UNTIL GET ELETRON ID WORKING IN 210
160 
161  iEvent.getByToken(e_idAssocProd_, pEleID);
162 
163  if (!pEleID.isValid()){
164  edm::LogInfo("")<< "Error! Can't get electronIDAssocProducer by label. ";
165  e_doID_ = false;
166  }
167  }
169  iEvent.getByToken(e_ctfTrackCollection_, pCtfTracks);
170  if (!pCtfTracks.isValid()) {
171  edm::LogInfo("")<< "Error! Can't get " << e_ctfTrackCollectionSrc_.label() << " by label. ";
172  iEvent.put(product_Electrons,"Electrons");
173  return;
174  }
175  const reco::TrackCollection * ctfTracks = pCtfTracks.product();
177  if(iEvent.getByToken(Electrons_,electrons))
178  for(size_t i=0;i<electrons->size();++i)
179  {
180  edm::Ref<reco::GsfElectronCollection> electronRef(electrons,i);
181  bool idDec=false;
182  if(e_doID_){
184  tagIDAssocItr = pEleID->find(electronRef);
185  const reco::ElectronIDRef& id_tag = tagIDAssocItr->val;
186  idDec=id_tag->cutBasedDecision();
187  }else idDec=true;
188  if((*electrons)[i].pt()>ptMinElectron_&&fabs((*electrons)[i].eta())<etaMax&&idDec)
189  {
190  if(e_doTrackIso_){
191  reco::TrackCollection::const_iterator tr = ctfTracks->begin();
192  double sum_of_pt_ele=0;
193  for(;tr != ctfTracks->end();++tr)
194  {
195  double lip = (*electrons)[i].gsfTrack()->dz() - tr->dz();
196  if(tr->pt() > e_trackMinPt_ && fabs(lip) < e_lipCut_){
197  double dphi=fabs(tr->phi()-(*electrons)[i].trackMomentumAtVtx().phi());
198  if(dphi>acos(-1.))dphi=2*acos(-1.)-dphi;
199  double deta=fabs(tr->eta()-(*electrons)[i].trackMomentumAtVtx().eta());
200  double dr_ctf_ele = sqrt(deta*deta+dphi*dphi);
201  if((dr_ctf_ele>e_minIsoDR_) && (dr_ctf_ele<e_maxIsoDR_)){
202  double cft_pt_2 = (tr->pt())*(tr->pt());
203  sum_of_pt_ele += cft_pt_2;
204  }
205  }
206  }
207  double isolation_value_ele = sum_of_pt_ele/((*electrons)[i].trackMomentumAtVtx().Rho()*(*electrons)[i].trackMomentumAtVtx().Rho());
208  if(isolation_value_ele<e_isoMaxSumPt_){
209  LorentzVector vec((*electrons)[i].px(),(*electrons)[i].py(),(*electrons)[i].pz(),(*electrons)[i].energy());
210  product_Electrons->push_back(vec);
211  }
212 
213  }
214  else{
215  LorentzVector vec((*electrons)[i].px(),(*electrons)[i].py(),(*electrons)[i].pz(),(*electrons)[i].energy());
216  product_Electrons->push_back(vec);
217  }
218  }
219  }
220 
221  iEvent.put(product_Electrons,"Electrons");
222 }
int i
Definition: DBlmapReader.cc:9
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
const std::string id_tag("id")
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:13
T eta() const
edm::EDGetTokenT< reco::GsfElectronCollection > Electrons_
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
edm::EDGetTokenT< reco::TrackCollection > e_ctfTrackCollection_
T sqrt(T t)
Definition: SSEVec.h:48
bool isValid() const
Definition: HandleBase.h:75
std::vector< LorentzVector > LorentzVectorCollection
edm::EDGetTokenT< reco::ElectronIDAssociationCollection > e_idAssocProd_
edm::InputTag e_ctfTrackCollectionSrc_
T const * product() const
Definition: Handle.h:81
std::string const & label() const
Definition: InputTag.h:42
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
void HLTTauRefProducer::doElectronsFromZ ( edm::Event ,
const edm::EventSetup ,
std::auto_ptr< LorentzVectorCollection > &   
)
private
void HLTTauRefProducer::doJets ( edm::Event iEvent,
const edm::EventSetup iES 
)
private

Definition at line 249 of file HLTTauRefProducer.cc.

References relval_parameters_module::energy, eta(), HLT_25ns14e33_v1_cff::etaMax, edm::Event::getByToken(), i, fwrapper::jets, and edm::Event::put().

250 {
251  auto_ptr<LorentzVectorCollection> product_Jets(new LorentzVectorCollection);
252  //Retrieve the collection
254  if(iEvent.getByToken(Jets_,jets))
255  for(size_t i = 0 ;i<jets->size();++i)
256  {
257  if((*jets)[i].et()>ptMinJet_&&fabs((*jets)[i].eta())<etaMax)
258  {
259  LorentzVector vec((*jets)[i].px(),(*jets)[i].py(),(*jets)[i].pz(),(*jets)[i].energy());
260  product_Jets->push_back(vec);
261  }
262  }
263  iEvent.put(product_Jets,"Jets");
264 }
int i
Definition: DBlmapReader.cc:9
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
T eta() const
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
vector< PseudoJet > jets
std::vector< LorentzVector > LorentzVectorCollection
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
edm::EDGetTokenT< reco::CaloJetCollection > Jets_
void HLTTauRefProducer::doMuons ( edm::Event iEvent,
const edm::EventSetup iES 
)
private

Definition at line 225 of file HLTTauRefProducer.cc.

References relval_parameters_module::energy, eta(), HLT_25ns14e33_v1_cff::etaMax, edm::Event::getByToken(), i, patZpeak::muons, and edm::Event::put().

226 {
227  auto_ptr<LorentzVectorCollection> product_Muons(new LorentzVectorCollection);
228  //Retrieve the collection
230  if(iEvent.getByToken(Muons_,muons))
231 
232  for(size_t i = 0 ;i<muons->size();++i)
233  {
234 
235  if((*muons)[i].pt()>ptMinMuon_&&fabs((*muons)[i].eta())<etaMax)
236  {
237  LorentzVector vec((*muons)[i].px(),(*muons)[i].py(),(*muons)[i].pz(),(*muons)[i].energy());
238  product_Muons->push_back(vec);
239  }
240  }
241 
242 
243  iEvent.put(product_Muons,"Muons");
244 
245 }
int i
Definition: DBlmapReader.cc:9
edm::EDGetTokenT< reco::MuonCollection > Muons_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
T eta() const
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
std::vector< LorentzVector > LorentzVectorCollection
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
tuple muons
Definition: patZpeak.py:38
void HLTTauRefProducer::doPFTaus ( edm::Event iEvent,
const edm::EventSetup iES 
)
private

Definition at line 121 of file HLTTauRefProducer.cc.

References relval_parameters_module::energy, eta(), HLT_25ns14e33_v1_cff::etaMax, edm::Event::getByToken(), i, and edm::Event::put().

122 {
123  auto_ptr<LorentzVectorCollection> product_PFTaus(new LorentzVectorCollection);
124  //Retrieve the collection
126  if(iEvent.getByToken(PFTaus_,pftaus)) {
127  for(unsigned int i=0; i<pftaus->size(); ++i) {
128  if((*pftaus)[i].pt()>ptMinPFTau_&&fabs((*pftaus)[i].eta())<etaMax)
129  {
130  reco::PFTauRef thePFTau(pftaus,i);
131  bool passAll = true;
134  if(iEvent.getByToken(token, pftaudis)) {
135  if((*pftaudis)[thePFTau] < 0.5) {
136  passAll = false;
137  break;
138  }
139  }
140  }
141  if(passAll) {
142  LorentzVector vec((*pftaus)[i].px(),(*pftaus)[i].py(),(*pftaus)[i].pz(),(*pftaus)[i].energy());
143  product_PFTaus->push_back(vec);
144  }
145  }
146  }
147  }
148  iEvent.put(product_PFTaus,"PFTaus");
149 }
int i
Definition: DBlmapReader.cc:9
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
T eta() const
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
edm::EDGetTokenT< reco::PFTauCollection > PFTaus_
std::vector< LorentzVector > LorentzVectorCollection
std::vector< edm::EDGetTokenT< reco::PFTauDiscriminator > > PFTauDis_
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
void HLTTauRefProducer::doPhotons ( edm::Event iEvent,
const edm::EventSetup iES 
)
private

Definition at line 297 of file HLTTauRefProducer.cc.

References relval_parameters_module::energy, eta(), HLT_25ns14e33_v1_cff::etaMax, edm::Event::getByToken(), i, interactiveExample::photons, and edm::Event::put().

298 {
299  auto_ptr<LorentzVectorCollection> product_Gammas(new LorentzVectorCollection);
300  //Retrieve the collection
302  if(iEvent.getByToken(Photons_,photons))
303  for(size_t i = 0 ;i<photons->size();++i)
304  if((*photons)[i].ecalRecHitSumEtConeDR04()<photonEcalIso_)
305  {
306  if((*photons)[i].et()>ptMinPhoton_&&fabs((*photons)[i].eta())<etaMax)
307  {
308  LorentzVector vec((*photons)[i].px(),(*photons)[i].py(),(*photons)[i].pz(),(*photons)[i].energy());
309  product_Gammas->push_back(vec);
310  }
311  }
312  iEvent.put(product_Gammas,"Photons");
313 }
int i
Definition: DBlmapReader.cc:9
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
T eta() const
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
edm::EDGetTokenT< reco::PhotonCollection > Photons_
std::vector< LorentzVector > LorentzVectorCollection
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
void HLTTauRefProducer::doTowers ( edm::Event iEvent,
const edm::EventSetup iES 
)
private

Definition at line 267 of file HLTTauRefProducer.cc.

References HLT_25ns14e33_v1_cff::DeltaR, relval_parameters_module::energy, eta(), HLT_25ns14e33_v1_cff::etaMax, edm::Event::getByToken(), i, j, p4, edm::Event::put(), and HLT_25ns14e33_v1_cff::towers.

268 {
269  auto_ptr<LorentzVectorCollection> product_Towers(new LorentzVectorCollection);
270  //Retrieve the collection
272  if(iEvent.getByToken(Towers_,towers))
273  for(size_t i = 0 ;i<towers->size();++i)
274  {
275  if((*towers)[i].pt()>ptMinTower_&&fabs((*towers)[i].eta())<etaMax)
276  {
277  //calculate isolation
278  double isolET=0;
279  for(unsigned int j=0;j<towers->size();++j)
280  {
281  if(ROOT::Math::VectorUtil::DeltaR((*towers)[i].p4(),(*towers)[j].p4())<0.5)
282  isolET+=(*towers)[j].pt();
283  }
284  isolET-=(*towers)[i].pt();
285  if(isolET<towerIsol_)
286  {
287  LorentzVector vec((*towers)[i].px(),(*towers)[i].py(),(*towers)[i].pz(),(*towers)[i].energy());
288  product_Towers->push_back(vec);
289  }
290  }
291  }
292  iEvent.put(product_Towers,"Towers");
293 }
int i
Definition: DBlmapReader.cc:9
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
T eta() const
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
double p4[4]
Definition: TauolaWrapper.h:92
int j
Definition: DBlmapReader.cc:9
edm::EDGetTokenT< CaloTowerCollection > Towers_
std::vector< LorentzVector > LorentzVectorCollection
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< float > > LorentzVector
Definition: analysisEnums.h:9
double HLTTauRefProducer::ElectronTrkIsolation ( const reco::TrackCollection ,
const reco::GsfElectron  
)
private
void HLTTauRefProducer::produce ( edm::Event iEvent,
const edm::EventSetup iES 
)
virtual

Implements edm::EDProducer.

Definition at line 103 of file HLTTauRefProducer.cc.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

104 {
105  if(doPFTaus_)
106  doPFTaus(iEvent,iES);
107  if(doElectrons_)
108  doElectrons(iEvent,iES);
109  if(doMuons_)
110  doMuons(iEvent,iES);
111  if(doJets_)
112  doJets(iEvent,iES);
113  if(doPhotons_)
114  doPhotons(iEvent,iES);
115  if(doTowers_)
116  doTowers(iEvent,iES);
117 
118 }
void doPFTaus(edm::Event &, const edm::EventSetup &)
void doTowers(edm::Event &, const edm::EventSetup &)
void doJets(edm::Event &, const edm::EventSetup &)
void doMuons(edm::Event &, const edm::EventSetup &)
void doPhotons(edm::Event &, const edm::EventSetup &)
void doElectrons(edm::Event &, const edm::EventSetup &)

Member Data Documentation

bool HLTTauRefProducer::doElecFromZ_
private

Definition at line 79 of file HLTTauRefProducer.h.

bool HLTTauRefProducer::doElectrons_
private

Definition at line 67 of file HLTTauRefProducer.h.

bool HLTTauRefProducer::doJets_
private

Definition at line 95 of file HLTTauRefProducer.h.

bool HLTTauRefProducer::doMuons_
private

Definition at line 90 of file HLTTauRefProducer.h.

bool HLTTauRefProducer::doPFTaus_
private

Definition at line 62 of file HLTTauRefProducer.h.

bool HLTTauRefProducer::doPhotons_
private

Definition at line 84 of file HLTTauRefProducer.h.

bool HLTTauRefProducer::doTowers_
private

Definition at line 99 of file HLTTauRefProducer.h.

edm::EDGetTokenT<reco::TrackCollection> HLTTauRefProducer::e_ctfTrackCollection_
private

Definition at line 69 of file HLTTauRefProducer.h.

edm::InputTag HLTTauRefProducer::e_ctfTrackCollectionSrc_
private

Definition at line 70 of file HLTTauRefProducer.h.

bool HLTTauRefProducer::e_doID_
private

Definition at line 72 of file HLTTauRefProducer.h.

bool HLTTauRefProducer::e_doTrackIso_
private

Definition at line 73 of file HLTTauRefProducer.h.

double HLTTauRefProducer::e_FromZet_
private

Definition at line 81 of file HLTTauRefProducer.h.

edm::EDGetTokenT<reco::ElectronIDAssociationCollection> HLTTauRefProducer::e_idAssocProd_
private

Definition at line 68 of file HLTTauRefProducer.h.

double HLTTauRefProducer::e_isoMaxSumPt_
private

Definition at line 78 of file HLTTauRefProducer.h.

double HLTTauRefProducer::e_lipCut_
private

Definition at line 75 of file HLTTauRefProducer.h.

double HLTTauRefProducer::e_maxIsoDR_
private

Definition at line 77 of file HLTTauRefProducer.h.

double HLTTauRefProducer::e_minIsoDR_
private

Definition at line 76 of file HLTTauRefProducer.h.

double HLTTauRefProducer::e_trackMinPt_
private

Definition at line 74 of file HLTTauRefProducer.h.

double HLTTauRefProducer::e_zMmax_
private

Definition at line 80 of file HLTTauRefProducer.h.

double HLTTauRefProducer::e_zMmin_
private

Definition at line 80 of file HLTTauRefProducer.h.

edm::EDGetTokenT<reco::GsfElectronCollection> HLTTauRefProducer::Electrons_
private

Definition at line 66 of file HLTTauRefProducer.h.

double HLTTauRefProducer::etaMax
private

Definition at line 103 of file HLTTauRefProducer.h.

edm::EDGetTokenT<reco::CaloJetCollection> HLTTauRefProducer::Jets_
private

Definition at line 94 of file HLTTauRefProducer.h.

edm::EDGetTokenT<reco::MuonCollection> HLTTauRefProducer::Muons_
private

Definition at line 89 of file HLTTauRefProducer.h.

std::vector<edm::EDGetTokenT<reco::PFTauDiscriminator> > HLTTauRefProducer::PFTauDis_
private

Definition at line 61 of file HLTTauRefProducer.h.

edm::EDGetTokenT<reco::PFTauCollection> HLTTauRefProducer::PFTaus_
private

Definition at line 60 of file HLTTauRefProducer.h.

double HLTTauRefProducer::photonEcalIso_
private

Definition at line 85 of file HLTTauRefProducer.h.

edm::EDGetTokenT<reco::PhotonCollection> HLTTauRefProducer::Photons_
private

Definition at line 83 of file HLTTauRefProducer.h.

double HLTTauRefProducer::ptMinElectron_
private

Definition at line 71 of file HLTTauRefProducer.h.

double HLTTauRefProducer::ptMinJet_
private

Definition at line 96 of file HLTTauRefProducer.h.

double HLTTauRefProducer::ptMinMuon_
private

Definition at line 91 of file HLTTauRefProducer.h.

double HLTTauRefProducer::ptMinPFTau_
private

Definition at line 63 of file HLTTauRefProducer.h.

double HLTTauRefProducer::ptMinPhoton_
private

Definition at line 86 of file HLTTauRefProducer.h.

double HLTTauRefProducer::ptMinTower_
private

Definition at line 100 of file HLTTauRefProducer.h.

double HLTTauRefProducer::towerIsol_
private

Definition at line 101 of file HLTTauRefProducer.h.

edm::EDGetTokenT<CaloTowerCollection> HLTTauRefProducer::Towers_
private

Definition at line 98 of file HLTTauRefProducer.h.