251 if (calojets.isValid()) {
253 mycalojets=*calojets;
255 typedef reco::CaloJetCollection::const_iterator cjiter;
257 for ( cjiter
i=mycalojets.begin();
i!=mycalojets.end();
i++) {
263 jcale[jcal] =
i->energy();
264 jcalemf[jcal] =
i->emEnergyFraction();
274 if (calocorjets.isValid()) {
276 mycalocorjets=*calocorjets;
278 typedef reco::CaloJetCollection::const_iterator ccorjiter;
280 for ( ccorjiter
i=mycalocorjets.begin();
i!=mycalocorjets.end();
i++) {
301 if(tower->energy() > thresholdForSavingTowers)
303 towet[jtow] = tower->et();
304 toweta[jtow] = tower->eta();
305 towphi[jtow] = tower->phi();
306 towen[jtow] = tower->energy();
307 towem[jtow] = tower->emEnergy();
308 towhd[jtow] = tower->hadEnergy();
309 towoe[jtow] = tower->outerEnergy();
317 if (recmets.isValid()) {
318 typedef reco::CaloMETCollection::const_iterator cmiter;
319 for ( cmiter
i=recmets->begin();
i!=recmets->end();
i++) {
327 typedef reco::METCollection::const_iterator iter;
328 for ( iter
i=ht->begin();
i!=ht->end();
i++) {
337 if (genjets.isValid()) {
341 typedef reco::GenJetCollection::const_iterator gjiter;
343 for ( gjiter
i=mygenjets.begin();
i!=mygenjets.end();
i++) {
349 jgene[jgen] =
i->energy();
358 if (genmets.isValid()) {
359 typedef reco::GenMETCollection::const_iterator gmiter;
360 for ( gmiter
i=genmets->begin();
i!=genmets->end();
i++) {
372 if (taujets.isValid()) {
377 typedef reco::HLTTauCollection::const_iterator tauit;
379 for(tauit
i=mytaujets.begin();
i!=mytaujets.end();
i++){
387 l25tauPt[itau] =
i->getL25LeadTrackPtValue();
404 typedef reco::PFTauCollection::const_iterator pftauit;
408 for(pftauit
i=taus.begin();
i!=taus.end();
i++){
413 ohpfTauJetPt[ipftau] =
i->pfTauTagInfoRef()->pfjetRef()->pt();
415 pfMHTx = pfMHTx +
i->pfTauTagInfoRef()->pfjetRef()->px();
416 pfMHTy = pfMHTy +
i->pfTauTagInfoRef()->pfjetRef()->py();
422 if( (
i->leadPFNeutralCand()).isNonnull())
424 if((
i->leadPFChargedHadrCand()).isNonnull())
426 float maxPtTrkIso = 0;
427 for (
unsigned int iTrk = 0; iTrk <
i->isolationPFChargedHadrCands().size(); iTrk++)
429 if(
i->isolationPFChargedHadrCands()[iTrk]->pt() > maxPtTrkIso) maxPtTrkIso =
i->isolationPFChargedHadrCands()[iTrk]->pt();
433 float maxPtGammaIso = 0;
434 for (
unsigned int iGamma = 0; iGamma <
i->isolationPFGammaCands().size(); iGamma++)
436 if(
i->isolationPFGammaCands()[iGamma]->pt() > maxPtGammaIso) maxPtGammaIso =
i->isolationPFGammaCands()[iGamma]->pt();
441 pfMHT =
sqrt(pfMHTx*pfMHTx + pfMHTy*pfMHTy);
447 if(recoPfTaus.isValid()) {
448 float minTrkPt = minPtCH;
449 float minGammaPt = minPtGamma;
453 typedef reco::PFTauCollection::const_iterator pftauit;
456 for(pftauit
i=taus.begin();
i!=taus.end();
i++){
463 if( (
i->leadPFNeutralCand()).isNonnull())
465 if((
i->leadPFChargedHadrCand()).isNonnull())
468 for (
unsigned int iTrk = 0; iTrk <
i->isolationPFChargedHadrCands().size(); iTrk++)
470 if(
i->isolationPFChargedHadrCands()[iTrk]->pt() > minTrkPt) myTrks++;
475 for (
unsigned int iGamma = 0; iGamma <
i->isolationPFGammaCands().size(); iGamma++)
477 if(
i->isolationPFGammaCands()[iGamma]->pt() > minGammaPt) myGammas++;
483 if(theRecoPFTauDiscrByTanCOnePercent.
isValid()){
485 if(theRecoPFTauDiscrByIsolation.
isValid()){
487 if(theRecoPFTauDiscrAgainstMuon.
isValid()){
489 if(theRecoPFTauDiscrAgainstElec.
isValid()){
491 if(theRecoPFTauDiscrByTanCHalfPercent.
isValid()){
493 if(theRecoPFTauDiscrByTanCQuarterPercent.
isValid()){
495 if(theRecoPFTauDiscrByTanCTenthPercent.
isValid()){
507 typedef reco::PFJetCollection::const_iterator pfJetit;
511 for(pfJetit
i=Jets.begin();
i!=Jets.end();
i++){
517 pfMHTx = pfMHTx +
i->px();
518 pfMHTy = pfMHTy +
i->py();
521 pfMHT =
sqrt(pfMHTx*pfMHTx + pfMHTy*pfMHTy);
528 nrpj = recoPFJets->size();
531 typedef reco::PFJetCollection::const_iterator pfJetit;
533 for(pfJetit
i=Jets.begin();
i!=Jets.end();
i++){
float * recopfTauDiscrAgainstElec
std::vector< PFTau > PFTauCollection
collection of PFTau objects
int * jpfrecochargedMultiplicity
float * ohpfTauLeadTrackPt
float * recopfTauDiscrByTancQuarterPercent
std::vector< GenJet > GenJetCollection
collection of GenJet objects
float * jpfrecochargedHadronFraction
std::vector< CaloTower >::const_iterator const_iterator
float * recopfTauDiscrByIso
float * recopfTauLeadTrackPt
float * recopfTauDiscrByTancHalfPercent
float * jpfreconeutralHadronFraction
float * recopfTauLeadPionPt
float * recopfTauDiscrAgainstMuon
std::vector< HLTTau > HLTTauCollection
collection of HLTTau objects
float * recopfTauDiscrByTancOnePercent
float * ohpfTauLeadPionPt
float * jpfreconeutralEMFraction
int * jpfreconeutralMultiplicity
std::vector< PFJet > PFJetCollection
collection of PFJet objects
float * jpfrecochargedEMFraction
float * recopfTauDiscrByTancTenthPercent
std::vector< CaloJet > CaloJetCollection
collection of CaloJet objects