Reconstruction and classification of tau lepton decays with ILD

Abstract

(ProQuest: ... denotes formulae and/or non-USASCII text omitted; see image)

Tau lepton decays with up to two ......s in the final state - ......, ......, and ...... - are used to study the performance of the barrel region of the silicon-tungsten electromagnetic calorimeter (Si-W ECAL) of the International Large Detector (ILD) at the future ...... International Linear Collider. Correct reconstruction of the tau decay mode is crucial for constraining the spin state of tau lepton and measuring the Higgs boson CP state in ...... decays. We find that about 95 % of ......, and 90 % of ...... and ...... decays produced by the ...... process at an ...... beam energy of 125 GeV are correctly reconstructed. In a smaller ILD detector, with the inner Si-W ECAL radius reduced by about 20 %, these efficiencies are reduced by at most 2 %. The ...... mass resolution remains below 10 %. Since failures in tau lepton reconstruction are mainly due to photons, an increase of the ILD magnetic field from 3.5 to 4 T does not bring any significant improvement.

Eur. Phys. J. C (2016) 76:468DOI 10.1140/epjc/s10052-016-4315-2

Special Article - Tools for Experiment and Theory

http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-016-4315-2&domain=pdf

Web End = http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-016-4315-2&domain=pdf

Web End = Reconstruction and classication of tau lepton decays with ILD

T. H. Trana, V. Balagura, V. Boudry, J.-C. Brient, H. Videau

Laboratoire Leprince-Ringuet, cole polytechnique, CNRS/IN2P3, Universit Paris-Saclay, 91128 Palaiseau, France

Received: 21 October 2015 / Accepted: 10 August 2016 / Published online: 20 August 2016 The Author(s) 2016. This article is published with open access at Springerlink.com

Abstract Tau lepton decays with up to two 0s in the nal state + +

, +(+0)

, and a+1(+00)

The reconstruction of the tau decay photons is therefore particularly suited as a benchmark of the ECAL performance.

In this paper we study the reconstruction and classication of + decay modes1 in the international large detector (ILD), a detector concept for the e+e international linear collider (ILC). We concentrate on the reconstruction of three main tau decays with one charged pion and up to two 0s:

+ +

, +

, and a+1

, followed by + +0,

a+1 +00, and 0 . For simplicity, we study

almost monochromatic tau leptons from the decay of the virtual Z0, e+e Z0 + at an e beam

energy of 125 GeV. This corresponds to a centre-of-mass energy of 250 GeV, modied by beamstrahlung and radiative corrections.

The momenta of the tau leptons from this process are larger than, for example, in e+e H0Z0, H0 + at

the same centre-of-mass energy. Therefore, our results correspond to Higgs boson production at higher ILC energies.

Only the task of the correct classication of the three + decay modes mentioned above is considered in this paper.

Backgrounds are not considered, and no attempt is made to extract the + spin state from the momenta of the reconstructed decay products.

The ILD detector has been optimized for the particle ow algorithm approach (PFA) [3,4] to event reconstruction. It combines very lightweight tracker systems with highly granular, imaging, calorimeters. A complete description of ILD is given in [5,6]. In this study the ILD performance with a baseline design is simulated with the software framework Mokka [7] with a parametric geometry description, based on the GEANT4 package [8]. The charged pion from + decay is reconstructed in the silicon vertex and tracker detectors and the time projection chamber. Photon clusters are reconstructed in the Si-W ECAL. We limit our analysis to the barrel part of ECAL, to avoid complications in the transition region

1 The inclusion of charge-conjugate states is implied in this article.

are used to study the performance of the barrel region of the silicon-tungsten electromagnetic calorimeter (Si-W ECAL) of the International Large Detector (ILD) at the future e+e

International Linear Collider. Correct reconstruction of the tau decay mode is crucial for constraining the spin state of tau lepton and measuring the Higgs boson CP state in H + decays. We nd that about 95 % of +

, and

90 % of +

and a+1

decays produced by the e+e

Z0 + process at an e beam energy of 125 GeV

are correctly reconstructed. In a smaller ILD detector, with the inner Si-W ECAL radius reduced by about 20 %, these efciencies are reduced by at most 2 %. The 0 mass resolution remains below 10 %. Since failures in tau lepton reconstruction are mainly due to photons, an increase of the ILD magnetic eld from 3.5 to 4 T does not bring any signicant improvement.

1 Introduction

Tau lepton decays have been used at LEP [1] and many other experiments for precise tests of the Standard Model. The spin state of the tau lepton can be inferred from its decay mode and the momenta of all reconstructed decay products [2]. This will be used at future lepton colliders, for example, to measure the CP (the product of charge conjugation and parity symmetries) state of the Higgs boson decaying into a tau pair, H0 +. The tau decay

products are usually more collimated than QCD jets of similar energy, and the separation of photon clusters in the electromagnetic calorimeter (ECAL) is more difcult. Incorrect determination of the number of photons in the tau nal state results in the wrong reconstruction of the tau decay mode and degrades the spin state measurement.

a e-mail: mailto:[email protected]

Web End [email protected]

123

468 Page 2 of 5 Eur. Phys. J. C (2016) 76 :468

between the barrel and endcap. The hadronic calorimeter is not used.

The Si-W ECAL has a thickness of about 23 radiation lengths at normal incidence. The silicon sensors have a pixelization of 5 5 mm2 and a thickness of 500 m. Due to

their high cost, and the large instrumented area, Si-W ECAL is one of the most expensive sub-detectors in ILD. Several studies of its cost optimization have been performed. In [6] it has been shown that a reduction of the number of ECAL layers from 30 to 26 degrades the jet energy resolution by less than 5 % for jet energies in the range 45250 GeV. In a more recent paper [9], it was demonstrated that a reduction of the ECAL radius by about 20 % degrades the jet energy resolution by at most 8 % in the same energy range, while the total cost of ILD could be lowered by a factor of 1.5.

The current paper supplements the latter optimization study. In addition to the baseline ILD design with an ECAL inner radius of 1843 mm, we present the results on tau decay mode reconstruction in smaller ILD models with inner ECAL radii of 1615 and 1450 mm. These numbers have been chosen according to the size of the large industrially available silicon sensors in order to simplify the nal ECAL mechanical design. The ECAL length is also reduced in order to leave the radius-to-length aspect ratio unchanged. All other detector parameters, such as the ECAL total and layer thicknesses, pixel size, gap between barrel and endcaps, etc., are unchanged.

The degradation of the jet energy and track momentum resolutions in a smaller ILD may be compensated by a higher magnetic eld. Therefore, we have also simulated ILD performance with the magnetic eld increased from the nominal3.5 to 4 T.

2 Simulation and reconstruction

We used a sample of e+e Z0 +( ) decays at

250 GeV produced for the ILD Detailed Baseline Design report included in the ILC Technical Design Report [6].The tau decay is generated by the TAUOLA library [10].Events with an initial state radiation are removed at a generator level, but events with nal state radiation (FSR), Z0 + , are retained. Events with FSR in the tau

decays, + +

, +

, a+1

, are also retained.

The tau branching fractions of the three decays under study without FSR are listed in Table 1. They cover 45 % of all tau decay modes and about 70 % of hadronic decays.

Each of the two -leptons is reconstructed only inside a hemisphere dened by the axis pointing in the tau direction in the Z0 rest frame. Events in which a photon converts before the ECAL are removed by requiring that only one charged track is reconstructed in each hemisphere.

Table 1 Branching fractions of the three tau decays under study with one charged and up to two neutral pions in the nal state [11]

Decay mode Branching fraction [%]

+

10.83 0.06 +

(+ +0) 25.52 0.09 a+1

(a+1 +00) 9.30 0.11

2.1 The Garlic clustering algorithm

The Gamma Reconstruction at a Linear Collider Experiment (Garlic) package [12,13], version 3.0.3, is used to nd photon clusters. The calibration constants for converting the energy deposited in the silicon sensors to the full photon energy are determined using samples of 10 GeV photons. This is done separately for the baseline and two smaller ILD ECAL models to take into account the possible variation of the overall weight of non-instrumented ECAL zones. The difference between the calibration constants of the three models is less than 1 %.

2.2 Photon identication

After nding clusters in the ECAL with the Garlic package, their nal identication as photon showers is performed using the following algorithm, which was specically developed for the purposes of this analysis.

FSR photons are suppressed by requiring the photon energy to be greater than 0.5 GeV. This eliminates 70 % of FSR photons and only 6 % of photons from + +

and a+1

decays.

To distinguish genuine photons from fake clusters originating, for example, from charged + hadronic showers, a boosted decision tree (BDT) with gradient boost is used with the following input variables:

the distance from the cluster barycentre to the + track

at the ECAL front surface;

the cluster depth computed by the Garlic algorithm; the cluster transverse fractal dimension dened as the

logarithm of the number of hits in the cluster, N1, divided by the number of hits when the pixels are grouped into larger pseudo-cells of 4 4 pixels, N4:

F D4 = log(N1/N4)/ log(4). (1)

The fractal dimension, introduced in ref. [14], reects the shower density which is higher in electromagnetic showers.

the mean and standard deviation of the distribution of hit

energies in the cluster; and

123

Eur. Phys. J. C (2016) 76 :468 Page 3 of 5 468

(a)

(b)

Fig. 1 The photon energy asymmetry A = (Emax Emin)/(Emax + Emin) in the case of three (a) and two (b)

reconstructed photons. Solid and dashed curves are for: a

two decay modes + +

and + a+1

; and b

for the generated (Monte Carlo truth) and measured energies, respectively

the fraction of cluster energy deposited between 5 and 10

radiation lengths after the starting point of the cluster.

25 % of all tau decays are used to train the BDT, which is then applied to the remaining 75 % of tau decays. The photon clusters for which more than 50 % of the hit energy is produced by the true Monte Carlo photon from tau decay are used as a signal for a BDT training. The remaining clusters in the same data sample serve as background.

The number of photons from + decay without FSR should be even. Sometimes, however, three photons are reconstructed. This can happen if one of the photons from + +

, + +0, 0 decays is recon

structed as two clusters, or if the decay is accompanied by an FSR photon. The other possibility is the merging of two photons into a single cluster in + a+1

, a+1 +00

decay, or loss of one photon. In the rst case of + +

decay, the photon with smallest energy is soft and the asymmetry between maximal and minimal energies of three photons, dened as

A =

Emax Emin Emax + Emin

, (2)

peaks at one. The distribution of A for both + +

with one reconstructed photon. In the latter case it is assumed that the mass beyond the kinematical limit is due to an external FSR photon, and the event is accepted.

A small fraction of events is accepted by two classiers, mainly when the nal state with three reconstructed photons is compatible with both +

and a+1

decays. In such cases the nal decision is made based on Mreco: +

is pre-

and a+1

is shown in Fig. 1a. In this special case of three photons, if A > 0.8, the energy of the least energetic photon cluster is added to the most energetic one and the cluster is removed from the following analysis.

As a cross-check, the same asymmetry A is plotted in Fig. 1b for the case of two reconstructed photons. It is dominated by + +

decays. Since 0 has no spin, its decay products should have a at asymmetry distribution, ranging from zero to one since the photon is massless. This

is indeed the case for the true Monte Carlo photon energies, save for effects of detector acceptance at one (the solid curve in Fig. 1b). The dashed curve, obtained after reconstruction, shows that the clustering algorithm tends to favor an asymmetric energy partition between two nearby photons.

2.3 Classication of tau decay modes

To distinguish the three considered + decay modes, a second BDT with gradient boost is used with the following input variables:

the invariant mass of all reconstructed + decay products,

Mreco;

the number of reconstructed photons; their total invariant mass; the energy of the photons; and their distance from the reconstructed track at the ECAL

front surface.

For every decay mode a BDT is trained on 25 % of statistics using this decay as a signal and the other two as a background. The three classiers obtained in this way (one for each mode) are applied to the remaining 75 % of events.

If Mreco is larger than the tau lepton mass, the event is rejected, except if it is classied as +

ferred over a+1

if Mreco < 0.85 GeV, and vice versa. Events selected by other combinations of two classiers (+

and

123

468 Page 4 of 5 Eur. Phys. J. C (2016) 76 :468

+

or +

and a+1

) are assigned to the nal state with larger mass. No event passes all three classiers.

3 Results and conclusions

Figure 2 summarizes the probability of correct reconstruction of the tau decay modes + +

, +

and a+1

Reconstructed decay

B = 4.0 T

0.03

0.03

0.02

4.29

4.29

0.15

89.68

89.68

0.40

0.17

0.17

0.05

4.52

4.52

0.16

89.86

89.86

0.40

3.83

3.83

0.21

90.92

90.92

0.21

6.41

6.41

0.32

3.77

3.77

0.20

90.39

90.39

0.22

6.42

6.42

0.32

95.88

95.88

0.21

2.34

2.34

0.11

0.48

0.48

0.09

95.70

95.70

0.21

2.65

2.65

0.12

0.53

0.53

0.10

R inner

ECAL

R

inner =1615 mm

ECAL m

m

1

=

0

5

4

a

1

, in the

central barrel region of the Si-W ECAL, and in the absence of backgrounds. The results are presented for the ILD baseline design with an inner ECAL radius RinnerECAL = 1843 mm,

and for two smaller ILD models with RinnerECAL = 1615 and

1450 mm. The probability does not include the efciency of charged + reconstruction, which is taken to be 100 %.

The off-diagonal values correspond to the probabilities of incorrect reconstruction. For some events the decay mode could not be determined (either Mreco is larger than the tau lepton mass or the event is not accepted by any classier, as explained in Sect. 2.3), therefore the sums of the probabilities in columns and rows do not reach 100 %. The reconstruction of two tau leptons in Z0 +( ) events is performed

independently, so that a failure to reconstruct one does not affect the other. The quoted uncertainties are purely statistical.

The correct reconstruction probability is close to 90 % for +

and a+1

and above 95 % for +

decays. For the smallest ILD with RinnerECAL = 1450 mm the degradation of

the probabilities compared to the baseline design is less than 2 %.

The slightly degraded jet and track momentum resolution of smaller ILD models can be compensated by a higher magnetic eld. This increases the bending of charged tracks, and therefore the typical separation between charged hadrons pairs, and between charged hadrons and photons. We therefore also simulated the ILD performance with the magnetic

1

a

1

a

Simulated decay

Fig. 3 The same as in Fig. 2, but for the magnetic eld of 4 T

eld increased from the nominal 3.5 to 4 T. The corresponding reconstruction probabilities are shown in Fig. 3. The improvement is marginal, less than 1 %. This demonstrates that the failure rate is dominated by photon reconstruction which is almost independent of the magnetic eld.

The photon clusters from the decay of energetic 0s often overlap, which degrades the 0 mass resolution. The effect is stronger for smaller ILD ECAL models, which tend to have less separation between photons. We study this effect using a sample with two reconstructed photons dominated by +

(+ +0) decays, and therefore typically containing a

rather energetic 0. Figure 4a shows that the reconstructed invariant mass of the 0 increases with its energy. This is explained by the fact that the Garlic algorithm splits the overlapping showers geometrically into two disjoint groups of neighboring hits, neglecting the fact that some part of the cluster energy can propagate into the area of another cluster. This tends to overestimate the distance between the cluster barycenters and the opening angle between the photons. The distribution in Fig. 4a is t to a parabolic function and the corresponding correction is applied to recover the nominal 0 mass at all energies. The resulting mass resolution,

Reconstructed decay

1

a

B = 3.5 T

0.09

0.09

0.03

4.68

4.68

0.16

91.12

91.12

0.38

0.09

0.09

0.03

4.86

4.86

0.16

90.50

90.50

0.38

0.15

0.15

0.04

5.47

5.47

0.17

89.05

89.05

0.41

2.75

2.75

0.17

90.46

90.46

0.22

4.86

4.86

0.28

3.27

3.27

0.19

89.91

89.91

0.23

5.90

5.90

0.31

3.29

3.29

0.19

88.78

88.78

0.24

7.07

7.07

0.34

96.76

96.76

0.19

2.48

2.48

0.12

0.40

0.40

0.08

96.39

96.39

0.20

2.89

2.89

0.13

0.49

0.49

0.09

96.19

96.19

0.21

3.44

3.44

0.14

0.66

0.66

0.11

R inner

ECAL m

m

inner =1843 mm

ECAL m

m

R inner

ECAL

R

=

1

6

1

5

=

1

4

0

5

Simulated decay

1

a

1

a

1

a

Fig. 2 The probabilities of correct and wrong reconstruction of the tau decay modes + +

, +

and a+1

in percent, obtained for models with Si-W ECAL inner radius (from left to right) RinnerECAL = 1843

(baseline), 1615 and 1450 mm, and with the nominal magnetic eld of3.5 T. Only statistical uncertainties are shown

123

Eur. Phys. J. C (2016) 76 :468 Page 5 of 5 468

(a)

Fig. 4 a The dependence of the reconstructed mass of 0 on its energy before the correction. The black curve shows the parabolic t. b The corrected 0 mass and the resolution for ILD with three different ECAL radii

shown in Fig. 4b, has a symmetric shape. The relative 0 mass resolution is below 10 % for all considered ILD models. The difference in mass resolution between the nominal and smallest ILD models is about 1 %.

A proper reconstruction of the tau decay modes together with a good momentum measurement of the decay products is mandatory to reconstruct the spin states of tau lepton in the measurement of the Higgs boson CP state in the decay H0 +. The high probability of the correct decay

mode reconstruction demonstrated in this paper shows the full potential of the ILD, even with reduced size, for such a measurement.

Acknowledgments The authors would like to thank Dr. Daniel Jeans, University of Tokyo for the help with the Garlic package used in this study. This work is funded by the Physics Department of cole polytechnique, Palaiseau, France and the Physique des deux Innis et des Origines (P2IO) program.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecomm http://creativecommons.org/licenses/by/4.0/

Web End =ons.org/licenses/by/4.0/ ), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Funded by SCOAP3.

References

1. A. Heister et al. (ALEPH collaboration), Eur. Phys. J. C 20, 401 (2001). http://arxiv.org/abs/hep-ex/0104038

Web End =arXiv:hep-ex/0104038

2. See, for example, the paper S. Berge, W. Bernreuther, H. Spies-berger, Phys. Lett. B 727 488, (2013). http://arxiv.org/abs/1308.2674

Web End =arXiv:1308.2674 [hep-ph] and many references therein

3. J.-C. Brient, H. Videau, in Proccedings of APS/DPF/DBP Summer Study on the Future of Particle Physics (Snowmass, Colorado, 2002). http://arxiv.org/abs/hep-ex/0202004

Web End =arXiv:hep-ex/0202004

4. M.A. Thomson, Nucl. Instrum. Methods A 611, 25 (2009). http://arxiv.org/abs/0907.3577

Web End =arXiv:0907.3577 [physics.ins-det]

5. The ILD Concept Group, ILD Letter of Intent, DESY 2009-87, KEK 2009-6 (2010). http://arxiv.org/abs/1006.3396

Web End =arXiv:1006.3396 [hep-ex]

6. T. Behnke et al., The ILC Technical Design Report. Int. Linear Collid. 4 (2013). http://arxiv.org/abs/1306.6329

Web End =arXiv:1306.6329 [physics.ins-det]

7. P.M. de Freitas, H. Videau, Detector simulation with MOKKA/GEANT4: Present and future. in Proccedings of International Workshop on Linear Colliders (LCWS, Jeju Island, Korea, 2002), pp. 2630

8. S. Agostinelli et al., Nucl. Instrum. Methods. A 506, 250 (2003). SLAC-PUB-9350, FERMILAB-PUB-03-339

9. T.H. Tran, in Proccedings of International Workshop on Future Linear Colliders (LCWS13, Tokyo, Japan). http://arxiv.org/abs/1404.3173

Web End =arXiv:1404.3173 [physics.ins-det]

10. S. Jadach, J.H. Khn, Z. Was, Comput. Phys. Commun. 64, 275 (1990)

11. K.A. Olive et al., Chin. Phys. C 38, 090001 (2014)12. M. Reinhard, J.-C. Brient, in Proccedings of International Linear Collider Workshop (LCWS08 and ILC08, Chicago, Illinois). http://arxiv.org/abs/0902.3042

Web End =arXiv:0902.3042 [hep-ex]

13. D. Jeans, J.-C. Brient, M. Reinhard, J. Instrum. 7, P06003 (2012). http://arxiv.org/abs/1203.0774

Web End =arXiv:1203.0774 [physics.ins-det]

14. M. Ruan, D. Jeans, V. Boudry, J.-C. Brient, H. Videau, Phys. Rev. Lett. 112, 012001 (2014). http://arxiv.org/abs/1312.7662

Web End =arXiv:1312.7662 [physics.ins-det]

123