Tau Polarization in T1 & Charged Higgs Samples

AOD Samples

Top D3PDs are produced from these AODs:
mc11_7TeV.109851.Pythia_ttbar_bWincbHminus_Hminus130TauInc.merge.AOD.e825_s1310_s1300_r3043_r2993
mc11_7TeV.105200.T1_McAtNlo_Jimmy.merge.AOD.e835_s1272_s1274_r3043_r2993
They are located on OMEGA, here:
/group/atlas/data/tauD3PD/topD3PD/MC/makingD3PD

D3PD Maker on lxplus

Top D3PDs are produced on lxplus, with Athena release 17.0.5.6.1, in
/afs/cern.ch/user/s/sbedikia/testarea/AtlasPhysics-17.0.5.6.1
To set up the environment:
 export AtlasSetup=/afs/cern.ch/atlas/software/dist/AtlasSetup
 alias asetup='source $AtlasSetup/scripts/asetup.sh'
 cd /afs/cern.ch/user/s/sbedikia/testarea/AtlasPhysics-17.0.5.6.1
 asetup AtlasPhysics,17.0.5.6.1,here
And run the command:
Reco_trf.py \
autoConfiguration=everything --ignoreerrors=True RunNumber=999407 \
 preExec='from D3PDMakerConfig.D3PDMakerFlags import D3PDMakerFlags;D3PDMakerFlags.TruthWriteHadrons.set_Value_and_Lock(True)' \
 inputAODFile=XXX \
outputNTUP_TOPFile=XXX \
  maxEvents=XXX
The TopInputsD3PDMaker TWiki is used.

D3PD Maker on the grid

Top D3PDs from the AODs have also been produced on the grid. Use dq2-get
 
user.sbedikia.mc11_7TeV.105200.T1_McAtNlo_Jimmy.merge.AOD.e835_s1272_s1274_r3043_r2993.NTUP_TOP_24032012/
user.sbedikia.mc11_7TeV.109851.Pythia_ttbar_bWincbHminus_Hminus130TauInc.merge.AOD.e825_s1310_s1300_r3043_r2993.NTUP_TOP_2403/
They are also located on OMEGA, here:
/group/atlas/data/tauD3PD/topD3PD/MC/fromGrid

Top D3PDs

Branching Fractions

T1 sample represents all ttbar, except fully hadronic, decays

T1 sample = 58% of all ttbar decays

'hadronic decays' include W-> jets and W->tau_had nu

The fraction of all ttbar decays that have a lepton + hadronically decaying tau = 3.7%. The fraction of lepton+hadronically decaying tau in T1 = 3.7%/58% = 5%.

The Branching Ratio of Ws is: BR(W->e) : BR(W->mu): BR(W->tau): BR(W->jets) is 1:1:1:6 (W->ud & W->sc, 3 colors each)

Assumptions made in D3PD: in T1: BR(t->W b) = 100% in Charged Higgs sample: BR(t -> H b) = 1% BR(ChH -> tau nu ) = 100%

Selection of Top Quarks in mc_ container

In the T1 sample, exactly two top quarks in each event are selected, with status code =155.

In the Single Charged Higgs sample, exactly two top quarks in each event are selected, with status code =3

Selection of W and H- bosons

In the both samples, only bosons from the selected top quarks were considered. Bosons can radiate a number of times before decaying into non-boson children. Only the last boson in the decay chain which decays to non-boson children is saved. The others are discarded by the selection.

The T1 sample contains exactly two non-intermediate W bosons from top quarks(status code 155) per event.

The single charged Higgs sample contains events with exactly one W boson(either + and - charged) and one negatively-charged Higgs boson. Bosons with status code = 2 from top quarks with status code =3 are selected.

Selection of Taus & Tau Daughters

Final State Particles

A 1 Prong hadronically decaying tau, as defined in my code, has at least 1 final-state particle with pdgId >16, with exactly 1 final-state particle of charge = +/- 1.

Distributions in the T1 sample: event quantities

And now, with the W-> tb processes included:

Warning: The distribution below has a known error, but still provides partially useful information:

Observed in the T1 sample: 1969 taus from Ws in 4995 T1 events 350 tau decays to electrons of 1968 taus = 18% 354 tau decays to muons of 1968 taus =18% 976 1Prong tau decays of 1968 taus = 50% 1 2Prong tau decay, leaked from 3P bin; unphysical 287 3Prong tau decays of 1968 taus = 15% 976+1+287=1264 hadronic tau decays of 1968 taus = 64%

Distributions in the Charged Higgs sample: event quantities

Observed in the Charged Higgs sample: 9034 taus from Charged Higgs decays in 9040 Charged Higgs events 1603 tau decays to electrons of 9034 taus = 17.7% 1648 tau decays to muons of 9034 taus =18.2% 4458 1Prong tau decays of 9034 taus = 49.3% 1297 3Prong tau decays of 9034 taus = 14% 7 5Prong tau decays of 9034 < 0.1%

Distributions in the T1 sample (left) and Charged Higgs sample (right) for 1 Prong Truth-Level Taus

Charged:

1. The ChargedEnergy distribution shows the mc_E of the charged final-state particles( charged pions). T1 on left, Charged Higgs on right.

2. The ChargedPionPt distribution shows the mc_pt of the charged final-state particles(charged pions). T1 on left, Charged Higgs on right.

3. The ChargedPionEt distribution shows the transverse energy calculated from the 4-momentum of the charged final-state particle( charged pion). T1 on left, Charged Higgs on right.

Neutral(non-neutrino):

4. The NeutralEnergy distribution shows the mc_E of all final-state non-charged non-neutrino tau daughters. T1 on left, Charged Higgs on right.

5. The NeutralPt distribution shows the transverse momentum calculated from the sum of the 4-momenta of all final-state non-charged, non-neutrino tau daughters. T1 on left, Charged Higgs on right.

6. The NeutralEt distribution shows the transverse energy calculated from the sum of the 4-momenta of all final-state non-charged, non-neutrino tau daughters. T1 on left, Charged Higgs on right.

The Neutrino:

7. The NeutrinoEnergy distribution shows the mc_E of the tau neutrino tau daughter. T1 on left, Charged Higgs on right.

8. The NeutrinoPt distribution shows the mc_pt of the tau neutrino tau daughter. T1 on left, Charged Higgs on right.

Visible:

9. The VisiblePt distribution shows the the transverse momentum calculated from the sum of the 4-momenta of neutral(non-neutrino) and charged tau daughters. T1 on left, Charged Higgs on right.

10. The VisibleEt distribution shows the the transverse energy calculated from the sum of the 4-momenta of neutral(non-neutrino) and charged tau daughters. T1 on left, Charged Higgs on right.

11. The VisibleE distribution shows the the energy calculated from the sum of the 4-momenta of neutral(non-neutrino) and charged tau daughters. T1 on left, Charged Higgs on right.

Tau:

12. The TauEt distribution shows the transverse energy calculated from the 4-momentum of the generator-level tau. T1 on left, Charged Higgs on right. T1 on left, Charged Higgs on right.

13. The TauPt distribution shows the transverse momentum calculated from the 4-momentum of the generator-level tau. T1 on left, Charged Higgs on right.

14. The TauE distribution shows the energy calculated from the 4-momentum of the generator-level tau. T1 on left, Charged Higgs on right.

SumAll:

15. The SumAllEt distribution shows the transverse energy calculated from the sum of the 4-momenta of all tau daughters( visible and neutrino). T1 on left, Charged Higgs on right.

16. The SumAllE distribution shows the energy calculated from the sum of the 4-momenta of all tau daughters( visible and neutrino). T1 on left, Charged Higgs on right.

17. The SumAllPt distribution shows the transverse momentum calculated from the sum of the 4-momenta of all tau daughters( visible and neutrino). T1 on left, Charged Higgs on right.

Diff:

18. Per Tau, the difference:

 TauE - SumAllE
, in GeV. In the T1 sample:

In the Charged Higgs sample:

Charged Asymmetry in T1 sample (left) and Charged Higgs sample (right) for 1P truth-level taus

The UpsilonHisto distribution shows the charged asymmetry for 1P truth-level taus from the T1 sample. It was previously calculated incorrectly as

 2* chargedPt/TauEt -1 

The UpsilonHisto distributions (T1 on left, Charged Higgs on right) show the charged asymmetry for 1P truth-level taus calculated as

 2* chargedPt/VisiblePt -1 

The distributions from both samples are plotted below. Each is normalized to unity. The plot on the right is on log scale.

Reconstruction Level Plots: Charged Higgs sample

The charged asymmetry distributions are shown for 1Prong generator- level & reconstructed-level taus truth-matched(dR<0.2) to generator-level taus, from the Charged Higgs sample. The peak at upsilon=1 is smeared out by the momentum resolution.

The charged asymmetry for reconstruction-level truth-matched taus with tau_seedCalo_numTrack==1 is given as :

 Charged Asymmetry := 2( tau_leadtrkpt/ tau_pt) -1 

And now, normalized:

The reconstruction efficiency as a function of the charged asymmetry: For the following distribution below, the horizontal axis is charged asymmetry. The vertical axis is the fraction of 1 track reconstruction-level taus matched within

 dR<0.2 
to a 1 prong generator-level tau, in bins of 0.04 in upsilon.

The Identification efficiency as a function of charged asymmetry is given below. The horizontal axis is charged asymmetry. The vertical axis is the fraction of 1 track taus with

 BDTj>0.7 
matched within
dR<0.2 
to a 1 prong generator-level tau.

The Upsilon resolution function is calculated, per tau, to be

 := Upsilon_reco/Upsilon_gen  -1 
. All taus are included except those with Upsilon_gen=1, =0 & Upsilon_reco=1. The taus with no neutral pions, i.e. those with charged asymetry equal to unity are not included.

Collaboration

Discussions, technical help, ideas, etc. : C. Cuenca, Z. Czyczuła, S. Demers, M. Scherzer

Next Step: Quantifying the Sensitivity

About 1200 ttbar -> lepton+1Ptau_had events are expected in 5 inv fb of data. The ATLAS limit , at m_chHiggs =130GeV, of the BR(t->chHiggs b) ~1 - 2%.
 https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2012-011/ 
, Figure 7d. Create pseudo-data with ROOT TF1::GetRandom & use RooFit with the truth-level templates.

-- SusieBedikian - 22-Feb-2012

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng ChargedEnergy.png r2 r1 manage 8.7 K 2012-03-06 - 12:15 SusieBedikian  
PNGpng ChargedEnergy2.png r1 manage 8.7 K 2012-03-06 - 12:20 SusieBedikian  
PNGpng ChargedEnergyChH.png r1 manage 8.3 K 2012-03-15 - 12:09 SusieBedikian  
PDFpdf ChargedPionEt.pdf r2 r1 manage 1.7 K 2012-03-06 - 02:26 SusieBedikian  
PNGpng ChargedPionEt.png r2 r1 manage 8.2 K 2012-03-06 - 12:15 SusieBedikian  
PNGpng ChargedPionEt2.png r1 manage 8.2 K 2012-03-06 - 12:20 SusieBedikian  
PNGpng ChargedPionEtChH.png r1 manage 9.6 K 2012-03-15 - 12:09 SusieBedikian  
PDFpdf ChargedPionPt.pdf r1 manage 14.6 K 2012-03-01 - 16:40 SusieBedikian  
PNGpng ChargedPionPt.png r1 manage 8.2 K 2012-03-06 - 12:14 SusieBedikian  
PNGpng ChargedPionPt2.png r1 manage 8.2 K 2012-03-06 - 12:20 SusieBedikian  
PNGpng ChargedPionPtChH.png r1 manage 9.6 K 2012-03-15 - 12:11 SusieBedikian  
PDFpdf DecayMatrix.pdf r1 manage 13.7 K 2012-02-29 - 14:00 SusieBedikian  
PNGpng DecayMatrix2.png r1 manage 41.2 K 2012-03-06 - 12:40 SusieBedikian  
PNGpng DecayMatrix_wtb.png r1 manage 32.1 K 2012-03-08 - 13:02 SusieBedikian  
PDFpdf DecayPair.pdf r1 manage 15.5 K 2012-03-01 - 17:20 SusieBedikian  
PNGpng DecayPair.png r1 manage 11.5 K 2012-03-06 - 12:22 SusieBedikian  
PNGpng DecayPair2.png r1 manage 11.5 K 2012-03-06 - 12:24 SusieBedikian  
PNGpng DecayPaircopy.png r1 manage 47.6 K 2012-03-02 - 13:36 SusieBedikian  
PNGpng Diff.png r1 manage 6.4 K 2012-03-15 - 13:04 SusieBedikian  
PNGpng DiffChH.png r1 manage 5.6 K 2012-03-15 - 13:03 SusieBedikian  
PNGpng GenTauReco_rebinned.png r1 manage 22.9 K 2012-03-16 - 15:49 SusieBedikian  
PNGpng GeneratorLevelUpsilonSame.png r1 manage 15.0 K 2012-03-26 - 13:57 SusieBedikian  
PNGpng GeneratorLevelUpsilonSameLOG.png r1 manage 12.4 K 2012-03-26 - 14:07 SusieBedikian  
PNGpng IDEff12.png r1 manage 9.2 K 2012-03-19 - 00:02 SusieBedikian  
PDFpdf LepHadStable.pdf r1 manage 13.2 K 2012-03-01 - 16:39 SusieBedikian  
PNGpng LepHadStable.png r1 manage 8.5 K 2012-03-06 - 02:51 SusieBedikian  
PNGpng LepHadStable2.png r1 manage 8.5 K 2012-03-06 - 12:25 SusieBedikian  
PNGpng LepHadStableChH.png r1 manage 10.3 K 2012-03-15 - 12:12 SusieBedikian  
PDFpdf LeptonicHadronic.pdf r1 manage 13.1 K 2012-03-01 - 16:39 SusieBedikian  
PDFpdf NeutralEnergy.pdf r1 manage 14.2 K 2012-03-01 - 16:43 SusieBedikian  
PNGpng NeutralEnergy.png r1 manage 8.0 K 2012-03-06 - 02:52 SusieBedikian  
PNGpng NeutralEnergyChH.png r1 manage 8.4 K 2012-03-15 - 12:12 SusieBedikian  
PNGpng NeutralEt.png r1 manage 7.7 K 2012-03-06 - 02:52 SusieBedikian  
PNGpng NeutralEtChH.png r1 manage 7.9 K 2012-03-15 - 12:13 SusieBedikian  
PNGpng NeutralPt.png r1 manage 7.7 K 2012-03-06 - 02:52 SusieBedikian  
PNGpng NeutralPtChH.png r1 manage 8.0 K 2012-03-15 - 12:13 SusieBedikian  
PDFpdf NeutrinoEnergy.pdf r1 manage 14.5 K 2012-03-01 - 16:43 SusieBedikian  
PNGpng NeutrinoEnergy.png r1 manage 7.9 K 2012-03-06 - 12:08 SusieBedikian  
PNGpng NeutrinoEnergy2.png r1 manage 7.9 K 2012-03-06 - 12:18 SusieBedikian  
PNGpng NeutrinoEnergyChH.png r1 manage 8.7 K 2012-03-15 - 12:15 SusieBedikian  
PDFpdf NeutrinoPt.pdf r1 manage 14.7 K 2012-03-01 - 16:43 SusieBedikian  
PNGpng NeutrinoPt.png r1 manage 8.5 K 2012-03-06 - 12:11 SusieBedikian  
PNGpng NeutrinoPt2.png r1 manage 8.5 K 2012-03-06 - 12:19 SusieBedikian  
PNGpng NeutrinoPtChH.png r1 manage 9.7 K 2012-03-15 - 12:35 SusieBedikian  
PDFpdf OtherEnergy.pdf r1 manage 14.4 K 2012-03-01 - 16:42 SusieBedikian  
PDFpdf OtherPt.pdf r1 manage 19.0 K 2012-03-01 - 16:42 SusieBedikian  
PNGpng ReconstructionEff12.png r1 manage 9.9 K 2012-03-19 - 00:02 SusieBedikian  
PNGpng Reconstruction_Level_Charged_Asymmetry.png r1 manage 24.3 K 2012-03-16 - 17:27 SusieBedikian  
PNGpng ResolutionFRecoTM.png r1 manage 7.0 K 2012-03-19 - 00:01 SusieBedikian  
PDFpdf SumAllE.pdf r1 manage 14.6 K 2012-03-01 - 16:38 SusieBedikian  
PNGpng SumAllE.png r1 manage 8.7 K 2012-03-06 - 02:54 SusieBedikian  
PNGpng SumAllEChH.png r1 manage 9.2 K 2012-03-15 - 12:42 SusieBedikian  
PDFpdf SumAllEt.pdf r1 manage 14.8 K 2012-02-29 - 18:17 SusieBedikian  
PNGpng SumAllEt.png r1 manage 9.6 K 2012-03-06 - 02:55 SusieBedikian  
PNGpng SumAllEtChH.png r1 manage 8.5 K 2012-03-15 - 12:42 SusieBedikian  
PDFpdf SumAllPt.pdf r1 manage 14.8 K 2012-03-01 - 16:38 SusieBedikian  
PNGpng SumAllPt.png r1 manage 9.6 K 2012-03-06 - 02:55 SusieBedikian  
PNGpng SumAllPtChH.png r1 manage 8.6 K 2012-03-15 - 12:46 SusieBedikian  
PDFpdf TauDecayTypes_mesons.pdf r1 manage 14.4 K 2012-03-01 - 16:38 SusieBedikian  
PNGpng TauDecayTypes_mesons.png r1 manage 12.3 K 2012-03-06 - 02:55 SusieBedikian  
PNGpng TauDecayTypes_mesons2.png r1 manage 12.3 K 2012-03-06 - 12:25 SusieBedikian  
PNGpng TauE.png r1 manage 8.4 K 2012-03-06 - 02:56 SusieBedikian  
PNGpng TauEChH.png r1 manage 9.1 K 2012-03-15 - 12:40 SusieBedikian  
PDFpdf TauEt.pdf r1 manage 14.8 K 2012-03-01 - 16:49 SusieBedikian  
PNGpng TauEt.png r1 manage 9.4 K 2012-03-06 - 12:11 SusieBedikian  
PNGpng TauEt2.png r1 manage 9.4 K 2012-03-06 - 12:19 SusieBedikian  
PNGpng TauEtChH.png r1 manage 8.4 K 2012-03-15 - 12:40 SusieBedikian  
PDFpdf TauEt_andSumAllEt.pdf r1 manage 15.2 K 2012-03-01 - 16:50 SusieBedikian  
PNGpng TauEt_andSumAllEt.png r1 manage 9.9 K 2012-03-06 - 02:56 SusieBedikian  
PNGpng TauPt.png r1 manage 9.4 K 2012-03-06 - 02:57 SusieBedikian  
PNGpng TauPtChH.png r1 manage 8.4 K 2012-03-15 - 12:43 SusieBedikian  
PDFpdf UpsilonHisto.pdf r1 manage 14.1 K 2012-02-29 - 16:41 SusieBedikian  
PNGpng UpsilonHisto.png r1 manage 7.2 K 2012-03-06 - 12:17 SusieBedikian  
PDFpdf UpsilonHistoChH.pdf r1 manage 14.1 K 2012-02-29 - 16:41 SusieBedikian  
PNGpng UpsilonHisto_LOGScale.png r1 manage 20.2 K 2012-03-15 - 11:48 SusieBedikian  
PNGpng UpsilonHisto_SamePlot.png r1 manage 23.8 K 2012-03-15 - 11:48 SusieBedikian  
PNGpng UpsilonHisto_visibleptDenom.png r1 manage 7.2 K 2012-03-15 - 12:00 SusieBedikian  
PNGpng UpsilonHisto_visibleptDenomChH.png r1 manage 9.2 K 2012-03-15 - 11:50 SusieBedikian  
PDFpdf VisibleE.pdf r1 manage 14.5 K 2012-03-01 - 16:37 SusieBedikian  
PNGpng VisibleE.png r1 manage 7.8 K 2012-03-06 - 02:57 SusieBedikian  
PNGpng VisibleEChH.png r1 manage 8.4 K 2012-03-15 - 12:39 SusieBedikian  
PDFpdf VisibleEt.pdf r1 manage 14.7 K 2012-03-01 - 16:37 SusieBedikian  
PNGpng VisibleEt.png r1 manage 8.8 K 2012-03-06 - 02:58 SusieBedikian  
PNGpng VisibleEtChH.png r1 manage 9.2 K 2012-03-15 - 12:36 SusieBedikian  
PDFpdf VisiblePt.pdf r1 manage 14.7 K 2012-03-01 - 16:36 SusieBedikian  
PNGpng VisiblePt.png r1 manage 8.7 K 2012-03-06 - 02:58 SusieBedikian  
PNGpng VisiblePtChH.png r1 manage 9.2 K 2012-03-15 - 11:00 SusieBedikian  
Edit | Attach | Watch | Print version | History: r24 < r23 < r22 < r21 < r20 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r24 - 2012-03-27 - SusieBedikian
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback