LAr Reconstruction: Data vs MC (parabola)

(1)

LAr LAr Reconstruction: Reconstruction:

Data

Data vs vs MC (parabola) MC (parabola)

S.Paganis

S.Paganis (Wisconsin)(Wisconsin) With

With

Manuel,Isabelle,Martin,Karina,Walter Manuel,Isabelle,Martin,Karina,Walter,…,…

LArLAr H8 Meeting, CERN,H8 Meeting, CERN, 55--AprilApril--20052005

(2)

Motivation Motivation

Go through a full MC vs Go through a full MC vs Data analysis within Data analysis within Athena (preparation exercise)

Athena (preparation exercise)

Find differences and understand what is Find differences and understand what is missing from the MC before we have a missing from the MC before we have a

realistic comparison with data (complicated) realistic comparison with data (complicated)

To avoid “apples and oranges” comparisons

To avoid inducing systematics in the Data Analysis due to incomplete MC simulation (profiles, Charge coll, ADC2MeV)

Ambitious: try to make some Data Ambitious: try to make some Data vs vs MC MC comparison plots and draw some prelim

comparison plots and draw some prelim conclusions

conclusions

Example: the parabola reconstruction of the ADCpeak.

(3)

Program Flow (

Program Flow ( release 10.0.1 release 10.0.1 ): ):

Analysis C++

Package

MC: ADC2MEV happens here

Thanks to:

Manuel Galas

Final Physics Plots

jobOptions.G4Ctb_Dig.py

Reconstruction

ESD and CBNT Data: ADC2MEV here

Total Raw Energy (MeV) 180019002000210022002300240025002600

Entries/25.0MeV

0 0.02 0.04 0.06 0.08 0.1 0.12 0.14

CTB04 Data jobOptions.G4Ctb_Sim.py

TBAnalysis on ESD

miniCBNT

+G4Apps

(4)

Analysis Analysis

Run: 1000947 50GeV electronsRun: 1000947 50GeV electrons

Aug.Material scan, but this run has no extra material

Parabola Energy reconstructionParabola Energy reconstruction

50MeV “cubicADCcut” in LArRawChannelSimpleBuilder.cxx

µA2MEV numbers from EMTB

EMTB 3x3 clusteringEMTB 3x3 clustering

No cluster corrections, No Long. weigthsNo cluster corrections, No Long. weigths

No shower cuts yet. Cut only on Tile Energy and phiNo shower cuts yet. Cut only on Tile Energy and phi

MCMC

New “pythonized” version (powerful)

Charge collection corrections

Tried to get “correct” beam profile …

ADC2MEV in Digitization step (careful to choose parabola)

(5)

ADC2MEV (Data

ADC2MEV (Data vs vs MC) MC)

(t(t_drift_drift*W)/e *W)/e ∗∗ 1/SF1/SF

Injection Injection Resistor Resistor 38.147

38.147 uA/VoltuA/Volt

Ramps Ramps How:How:

Volts2 Volts2µµAA

370.3703 370.3703 37.58

37.58 nAnA

S3 (EMB1) S3 (EMB1)

370.3703 370.3703 37.58

37.58 nAnA

S2 (EMB1) S2 (EMB1)

370.3703 370.3703 12.62

12.62 nAnA

S1 (EMB1) S1 (EMB1)

12501250 38.147/R=0.114

38.147/R=0.114 nAnA

PS (EMB1) PS (EMB1)

µµA2MeVA2MeV DAC2Volt

DAC2Volt ADC2DAC

ADC2DAC

MC ADC2MEV(PS) = 7.0 MC ADC2MEV(S1) = 2.5 MC ADC2MEV(S2) = 12.0

Data ADC2MEV(PS) ~ 7.2 Data ADC2MEV(S1) ~ 2.4

Data ADC2MEV(S2) ~ vary 10.0,12.0,16.0

[ ^ADCpeak ^PEDESTAL ]

ADC2MEV ⋅ −

rec

=

E

(6)

How to go to the visible energy for Data:

an example (S.P.) an example (S.P.)

Accordion A

MeV/

947 .

328

Presampler A

MeV/

47 . 1176

8 . 0

|

Accordion A

MeV/

37 . 370

Presampler A

MeV/

1250

8 . 0

|

: tion Reconstruc

EMTBeam

µ µ η

>

<

C eV

ns

C eV

ns

19 19

10 6 . 1 / 6

. 23 470

A MeVper

: Accordion

10 6 . 1 / 6

. 23 420

A MeVper

: Presampler

−

⋅

=

⋅

= µ µ

SF(Presampler η<0.8)=t*W/e/1250 = 0.0496 SF(Accordion η<0.8)=t*W/e/370.37 = 0.18718

(7)

50 50 GeV GeV data data vs vs e and e and π π MC MC

3x3 cluster Energy (GeV)

0 10 20 30 40 50 60

0 100 200 300 400 500 600

Shapes between Data and π MC do not agree

Electron MC energy shifted by -2.5%

(8)

0 10 20 30 40 50 60

10 102

103

50 50 GeV GeV data data vs vs e and e and π π MC MC

Shapes between Data and π MC do not agree

(9)

Energy in Tile0

Energy in Tile0 vs vs Erec Erec in in LAr LAr (Data) (Data)

(GeV) ELAr

0 10 20 30 40 50

(GeV) 0Tile E

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Tail?

(10)

Beam Profile 0.320 0.33 0.34 0.35 0.36 0.37 0.38 0.39 η0.4 0.41 0.42 200

400 600 800 1000 1200 1400

Beam Profile

-0.020 -0.01 0 0.01 0.02 φ0.03 0.04

100 200 300 400 500 600

Beam profile: we need a match with MC Beam profile: we need a match with MC

Dip is not appearing in MC Data

(11)

Erec Erec after only Tile cut after only Tile cut (E (E

₀₀

<300MeV <300MeV ) )

40 42 44 46 48 50 52 54 56 58 60

0 100 200 300 400

500 ^RMS ^1.145

Underflow 3.252 Integral 2680

/ ndf

χ2 7.035 / 5

Prob 0.218

Constant 486.9 ±16.7 Mean 49.92 ± 0.03

Sigma 1.091 ± 0.035

RMS 1.145

/ ndf

χ2 7.035 / 5

Prob 0.218

Constant 486.9 ±16.7 Mean 49.92 ± 0.03

Sigma 1.091 ± 0.035

RMS 1.145

/ ndf

χ2 7.035 / 5

Prob 0.218

Constant 486.9 ±16.7 Mean 49.92 ± 0.03

Sigma 1.091 ± 0.035

Fit on the MC only!

Tail: Early Brem?

ADCpeak? Other?

Electron MC energy shifted by -2.5%

(12)

Erec Erec per per LAr LAr Layer for 50GeV e Layer for 50GeV e (we need the visible Energy)

(we need the visible Energy)

20 40 60 80 100 120 140 160 180

100 200 300 400 500 PreSampler: 3x3 cluster Energy (GeV)

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

0 50 100 150 200 250 300

Strips: 3x3 cluster Energy (GeV)

5 10 15 20 25 30

0 50 100 150 200

250 Strips: inconsistency (ADC2MEV lower than in MC)

Data

(13)

Strips: for

Strips: for Erec Erec >45GeV, Data show large >45GeV, Data show large fraction of Energy outside core

fraction of Energy outside core

3 - 1

7/E Strips: E

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

10-3

10-2

10-1

hfracm

1 1

3 −

±

strip E

strips

E

(14)

Check a strip (

Check a strip ( η,φ η,φ )=(117,2): )=(117,2):

usually part of the electron cluster usually part of the electron cluster

ID Entries Mean RMS

1 4005 24.14 31.02 51.01 / 13

Constant 104.5

Mean 1.910

Sigma 7.433

ID Entries Mean RMS

2 2000 58.14 65.33 22.05 / 13

Constant 89.20

Mean 19.99

Sigma 17.71

(15)

Summary/To check Summary/To check

(before cluster calibrations/corrections) (before cluster calibrations/corrections)

Need to check the agreement of the MC geometry Need to check the agreement of the MC geometry with the runs I am analyzing (true for everybody) with the runs I am analyzing (true for everybody)

Need to run MC with correct beam profileNeed to run MC with correct beam profile

Need to compare layer response visible energyNeed to compare layer response visible energy

Strip response needs to be understoodStrip response needs to be understood

Probably MC is not reflecting correctly the noise (including cross- talk)

Electron tail needs to be understoodElectron tail needs to be understood

Study different beam energiesStudy different beam energies

On the positive side:On the positive side: the nonthe non--tail region of the tail region of the electron energy is in agreement with MC giving a electron energy is in agreement with MC giving a

~1.1GeV Gaussian

~1.1GeV Gaussian σσ for 50GeVfor 50GeV

LAr Reconstruction: Data vs MC (parabola)

LAr LAr Reconstruction: Reconstruction:

Data

Data vs vs MC (parabola) MC (parabola)

Motivation Motivation

Go through a full MC vs Go through a full MC vs Data analysis within Data analysis within Athena (preparation exercise)

Athena (preparation exercise)

Find differences and understand what is Find differences and understand what is missing from the MC before we have a missing from the MC before we have a

realistic comparison with data (complicated) realistic comparison with data (complicated)

Ambitious: try to make some Data Ambitious: try to make some Data vs vs MC MC comparison plots and draw some prelim

comparison plots and draw some prelim conclusions

conclusions

Program Flow (

Program Flow ( release 10.0.1 release 10.0.1 ): ):

Analysis Analysis

ADC2MEV (Data

ADC2MEV (Data vs vs MC) MC)

[ ADCpeak PEDESTAL ]

ADC2MEV ⋅ −

=

E

How to go to the visible energy for Data:

How to go to the visible energy for Data:

an example (S.P.) an example (S.P.)

µ µ η

µ µ η

50 50 GeV GeV data data vs vs e and e and π π MC MC

50 50 GeV GeV data data vs vs e and e and π π MC MC

Energy in Tile0

Energy in Tile0 vs vs Erec Erec in in LAr LAr (Data) (Data)

Beam profile: we need a match with MC Beam profile: we need a match with MC

Erec Erec after only Tile cut after only Tile cut (E (E

<300MeV <300MeV ) )

Erec Erec per per LAr LAr Layer for 50GeV e Layer for 50GeV e (we need the visible Energy)

(we need the visible Energy)

Strips: for

Strips: for Erec Erec >45GeV, Data show large >45GeV, Data show large fraction of Energy outside core

fraction of Energy outside core

1 1

3 −

±

±

strip E

strips

E

Check a strip (

Check a strip ( η,φ η,φ )=(117,2): )=(117,2):

usually part of the electron cluster usually part of the electron cluster

Summary/To check Summary/To check

(before cluster calibrations/corrections) (before cluster calibrations/corrections)

[ ^ADCpeak ^PEDESTAL ]