ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Mods/src/PhotonTreeWriter.cc
(Generate patch)

Comparing UserCode/MitPhysics/Mods/src/PhotonTreeWriter.cc (file contents):
Revision 1.13 by fabstoec, Tue Apr 24 11:45:53 2012 UTC vs.
Revision 1.14 by fabstoec, Wed May 2 16:57:19 2012 UTC

# Line 431 | Line 431 | void PhotonTreeWriter::Process()
431      //printf("r9 = %5f, photon sigieie = %5f, seed sigieie = %5f\n",phHard->R9(),
432      //       phHard->CoviEtaiEta(),sqrt(phHard->SCluster()->Seed()->CoviEtaiEta()));
433  
434 +    // MuonStuff
435 +    fDiphotonEvent-> muonPt  = -99.;
436 +    fDiphotonEvent-> muonEta  = -99.;
437 +    fDiphotonEvent-> muDR1  = -99.;
438 +    fDiphotonEvent-> muDR2  = -99.;
439 +    fDiphotonEvent-> muIso1  = -99.;
440 +    fDiphotonEvent-> muIso2  = -99.;
441 +    fDiphotonEvent-> muIso3  = -99.;
442 +    fDiphotonEvent-> muIso4  = -99.;
443 +    fDiphotonEvent-> muD0  = -99.;
444 +    fDiphotonEvent-> muDZ  = -99.;
445 +    fDiphotonEvent-> muChi2  = -99.;  
446 +    fDiphotonEvent-> muNhits = -99;
447 +    fDiphotonEvent-> muNpixhits = -99;
448 +    fDiphotonEvent-> muNegs = -99;
449 +    fDiphotonEvent-> muNMatch = -99;
450 +
451 +    // Electron Stuff
452 +    fDiphotonEvent-> elePt = -99.;
453 +    fDiphotonEvent-> eleEta = -99.;
454 +    fDiphotonEvent-> eleSCEta = -99.;    
455 +    fDiphotonEvent-> eleIso1 = -99.;
456 +    fDiphotonEvent-> eleIso2 = -99.;
457 +    fDiphotonEvent-> eleIso3 = -99.;
458 +    fDiphotonEvent-> eleIso4 = -99.;
459 +    fDiphotonEvent-> eleDist = -99.;
460 +    fDiphotonEvent-> eleDcot = -99.;
461 +    fDiphotonEvent-> eleCoviee = -99.;
462 +    fDiphotonEvent-> eleDphiin = -99.;
463 +    fDiphotonEvent-> eleDetain = -99.;
464 +    fDiphotonEvent-> eleDR1 = -99.;
465 +    fDiphotonEvent-> eleDR2 = -99.;
466 +    fDiphotonEvent-> eleMass1 = -99.;
467 +    fDiphotonEvent-> eleMass2 = -99.;
468 +    fDiphotonEvent-> eleNinnerHits = -99;    
469 +
470      if( fApplyLeptonTag ) {
471        // perform lepton tagging
472        // the diphoton event record will have one more entry; i.e. leptonTag
# Line 441 | Line 477 | void PhotonTreeWriter::Process()
477        fDiphotonEvent->leptonTag = 0;
478  
479        if ( fLeptonTagMuons->GetEntries() > 0 ) {
480 +
481          // need to have dR > 1 for with respect to both photons
482          for(UInt_t iMuon = 0; iMuon <fLeptonTagMuons->GetEntries(); ++iMuon) {
483            if(MathUtils::DeltaR(fLeptonTagMuons->At(iMuon),phHard) < 1.) continue;
484            if(MathUtils::DeltaR(fLeptonTagMuons->At(iMuon),phSoft) < 1.) continue;
485  
486            fDiphotonEvent->leptonTag = 1;
487 +
488 +          fDiphotonEvent-> muonPt  = fLeptonTagMuons->At(iMuon)->Pt();
489 +          fDiphotonEvent-> muonEta = fLeptonTagMuons->At(iMuon)->Eta();
490 +          fDiphotonEvent-> muDR1   = MathUtils::DeltaR(fLeptonTagMuons->At(iMuon),phHard);
491 +          fDiphotonEvent-> muDR2   = MathUtils::DeltaR(fLeptonTagMuons->At(iMuon),phSoft);
492 +          fDiphotonEvent-> muIso1   = (fLeptonTagMuons->At(iMuon)->IsoR03SumPt() + fLeptonTagMuons->At(iMuon)->IsoR03EmEt() + fLeptonTagMuons->At(iMuon)->IsoR03HadEt() - fPileUpDen->At(0)->RhoRandomLowEta() * TMath::Pi() * 0.3 * 0.3)/ fLeptonTagMuons->At(iMuon)->Pt();
493 +          fDiphotonEvent-> muIso2   = (fLeptonTagMuons->At(iMuon)->IsoR03SumPt() + fLeptonTagMuons->At(iMuon)->IsoR03EmEt() + fLeptonTagMuons->At(iMuon)->IsoR03HadEt() - fPileUpDen->At(0)->RhoRandom() * TMath::Pi() * 0.3 * 0.3)/ fLeptonTagMuons->At(iMuon)->Pt();
494 +          fDiphotonEvent-> muIso3   = (fLeptonTagMuons->At(iMuon)->IsoR03SumPt() + fLeptonTagMuons->At(iMuon)->IsoR03EmEt() + fLeptonTagMuons->At(iMuon)->IsoR03HadEt() - fPileUpDen->At(0)->RhoLowEta() * TMath::Pi() * 0.3 * 0.3)/ fLeptonTagMuons->At(iMuon)->Pt();
495 +          fDiphotonEvent-> muIso4   = (fLeptonTagMuons->At(iMuon)->IsoR03SumPt() + fLeptonTagMuons->At(iMuon)->IsoR03EmEt() + fLeptonTagMuons->At(iMuon)->IsoR03HadEt() - fPileUpDen->At(0)->Rho() * TMath::Pi() * 0.3 * 0.3)/ fLeptonTagMuons->At(iMuon)->Pt();
496 +          fDiphotonEvent-> muD0  = TMath::Abs(fLeptonTagMuons->At(iMuon)->BestTrk()->D0Corrected(*fPV->At(0)));
497 +          fDiphotonEvent-> muDZ  = TMath::Abs(fLeptonTagMuons->At(iMuon)->BestTrk()->DzCorrected(*fPV->At(0)));
498 +          fDiphotonEvent-> muChi2  = fLeptonTagMuons->At(iMuon)->GlobalTrk()->Chi2()/fLeptonTagMuons->At(iMuon)->GlobalTrk()->Ndof();
499 +                  
500 +          fDiphotonEvent-> muNhits = fLeptonTagMuons->At(iMuon)->BestTrk()->NHits();
501 +          fDiphotonEvent-> muNpixhits = fLeptonTagMuons->At(iMuon)->BestTrk()->NPixelHits();
502 +          fDiphotonEvent-> muNegs = fLeptonTagMuons->At(iMuon)->NSegments();
503 +          fDiphotonEvent-> muNMatch = fLeptonTagMuons->At(iMuon)->NMatches();
504 +
505            break;
506          }
507        }
# Line 460 | Line 515 | void PhotonTreeWriter::Process()
515            if ( TMath::Abs( (phSoft->Mom()+fLeptonTagElectrons->At(iEle)->Mom()).M()-91.19 ) < 5. ) continue;
516            
517            fDiphotonEvent->leptonTag = 2;
518 +
519 +          fDiphotonEvent-> elePt = fLeptonTagElectrons->At(iEle)->Pt();
520 +          fDiphotonEvent-> eleEta = fLeptonTagElectrons->At(iEle)->Eta();
521 +          fDiphotonEvent-> eleSCEta = fLeptonTagElectrons->At(iEle)->SCluster()->Eta();
522 +          fDiphotonEvent-> eleIso1 = (fLeptonTagElectrons->At(iEle)->TrackIsolationDr03() + fLeptonTagElectrons->At(iEle)->EcalRecHitIsoDr03() + fLeptonTagElectrons->At(iEle)->HcalTowerSumEtDr03() - fPileUpDen->At(0)->RhoRandomLowEta() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
523 +          fDiphotonEvent-> eleIso2 = (fLeptonTagElectrons->At(iEle)->TrackIsolationDr03() + fLeptonTagElectrons->At(iEle)->EcalRecHitIsoDr03() + fLeptonTagElectrons->At(iEle)->HcalTowerSumEtDr03() - fPileUpDen->At(0)->RhoRandom() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
524 +          fDiphotonEvent-> eleIso3 = (fLeptonTagElectrons->At(iEle)->TrackIsolationDr03() + fLeptonTagElectrons->At(iEle)->EcalRecHitIsoDr03() + fLeptonTagElectrons->At(iEle)->HcalTowerSumEtDr03() - fPileUpDen->At(0)->RhoLowEta() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
525 +          fDiphotonEvent-> eleIso4 = (fLeptonTagElectrons->At(iEle)->TrackIsolationDr03() + fLeptonTagElectrons->At(iEle)->EcalRecHitIsoDr03() + fLeptonTagElectrons->At(iEle)->HcalTowerSumEtDr03() - fPileUpDen->At(0)->Rho() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
526 +          fDiphotonEvent-> eleDist = fLeptonTagElectrons->At(iEle)->ConvPartnerDist();
527 +          fDiphotonEvent-> eleDcot = fLeptonTagElectrons->At(iEle)->ConvPartnerDCotTheta();
528 +          fDiphotonEvent-> eleCoviee = fLeptonTagElectrons->At(iEle)->CoviEtaiEta();
529 +          fDiphotonEvent-> eleDphiin = TMath::Abs(fLeptonTagElectrons->At(iEle)->DeltaPhiSuperClusterTrackAtVtx());
530 +          fDiphotonEvent-> eleDetain = TMath::Abs(fLeptonTagElectrons->At(iEle)->DeltaEtaSuperClusterTrackAtVtx());
531 +          fDiphotonEvent-> eleDR1 = MathUtils::DeltaR(fLeptonTagElectrons->At(iEle),phHard);
532 +          fDiphotonEvent-> eleDR2 = MathUtils::DeltaR(fLeptonTagElectrons->At(iEle),phSoft);
533 +          fDiphotonEvent-> eleMass1 = (phHard->Mom()+fLeptonTagElectrons->At(iEle)->Mom()).M();
534 +          fDiphotonEvent-> eleMass2 = (phSoft->Mom()+fLeptonTagElectrons->At(iEle)->Mom()).M();
535 +          fDiphotonEvent-> eleNinnerHits =      fLeptonTagElectrons->At(iEle)->Trk()->NExpectedHitsInner();
536 +
537            break;
538          }
539        }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines