Start Presentation

Slide 1: Used CSC data

Slide 2: Measuring ttbar spin correlation code for these slides

  • Measuring method:
    • Atlfast data: distribution of $cos\theta1*cos\theta2$ in reconstruction level divide that in parton level to obtain a correction function (pol(n) fit).
    • Full sim data: the distribution in reconstruction level divide the correction function to recover to its parton level.
    • The mean of the recovered distribution multiplying (-9/0.51) is taken as the estimator of the ttbar correlation parameter A
  • The expected value of A=0.42 in SM and measured value is 0.82+/-0.16
AcerMC_measuring.gif

Slide 3: Why it looks like bias

  • Two analysis with the same technique and entirely different samples
    • Analysis 1: "rome" TopRex, full simulation sample 730/pb and early fast simulation high statistics (10/fb?).
    • Analysis 2: "CSC" AcerMC, full simulation sample 220/pb and atlfast1 fast simulation 220/pb.
  • The measurement results (only statistical error)
    • Analysis 1: 0.67 +/- 0.10 deviation>$2\sigma$
    • Analysis 2: 0.82 +/- 0.16 deviation>$2\sigma$
  • SAME and LARGE deviation from the expected value on the same side. What does it mean:
    • Bias is introduced due to discrepancy between fast and full simulation Question
    • In fact, statistical error Question

Slide 4: Bias comes from full sim electron in crack? NO!

  • remove events with electrons reconstructed in calorimeter crack (bad reconstruction) (1.37 < abs(etaele) < 1.52 ), nothing changed for the results
  • this is not the reason

Slide 5: Bias comes from error of correction function? looks not.

estimation method Change $1\sigma$ for each parameter of the fitted correction function. 5 parameters$\rightarrow$many measurements. take $\sigma$ as error reset each bin with randomly gaussian according to its error bar and re-fit the correction histogram. do 5000 times, take $1\sigma$ as error Statistical error Systematic error total error exist bias
Previous analysis 1 0.09 plot2 --- 0.10 0.08* 0.16 0.22
now analysis 2. Which method is reliable? 0.39 plot1 --- 0.19 0.15** 0.45 0.40
--- 0.18 plot3 0.30
* the variation of some variable in error estimation takes 3%, ** 5% is taken

AcerMC_correctionfunc_error.gif MCNLO_correctionfunc_error_oldfast.gif AcerMC_estimate_error_from_correctiofunction2.gif

Slide 6: Compare correction functions from different samples

  • Difference are clear.
    • This difference is considered as one source of the bias.
    • We still don't understand why and what make the difference, except statisticsTIP.
inversed_correctionfunc_black_oldfast.gif

  • As stated below, 2 independent aflfast data give different measurement, here compare their correction functions
    • different bins are used to see the bins effect (code)
    • fast data 2 is the reprocessing of the same fast data, but change random seed.
6bins 7bins 9bins 13bins 19bins 25bins
AcerMC_Correction_function_6_bins.gif AcerMC_Correction_function_7_bins.gif AcerMC_Correction_function_9_bins.gif AcerMC_Correction_function_13_bins.gif AcerMC_Correction_function_19_bins.gif AcerMC_Correction_function_25_bins.gif

More fast data are produced, A are deduced from them (same data set, different random seed for b tagging)

datasets with trigger, correct with correction function with trigger, correction with histogram 25bins with trigger, correction with histogram 13bins with trigger, correction with histogram 7bins
fastdata1 0.710314 0.705644 0.695296 0.762789
fastdata2 0.452789 0.363568 0.333444 0.464422
fastdata3 0.923959 0.994478 0.613127 0.75341
fastdata4 0.437812 0.471541 0.386576 0.4666
fastdata5 0.354353 0.462082 0.374049 0.397047
fastdata6 0.540373 0.58203 0.496647 0.642595
fastdata7 0.443663 0.521005 0.278286 0.436368
fastdata8 0.616136 0.575533 0.527783 0.630022
fastdata9 0.452907 0.459436 0.368635 0.448882
average 0.548034 0.570590778 0.452649222 0.555792778
RMS 0.167008402 0.175541261 0.130990219 0.134222448

  • the correction functions:
    AcerMC_Correction_function2_9_bins.gif

Slide 7: Pure measuring resolution effect on correction function

  • Fast simulation: according to measurement resolution, we add a small value to the truth as its reconstructed value
    • We can do many times and make many different correction function.
  • The measurement distribution is as blow
    • Take $1\sigma$ as error: 0.15 (left plot) for analysis2. 0.07 for MCNLO (right plot) means it's smaller with higher statistics,
    • so <<0.07 for Analysis1.
  • These error are actually part of error above, can't explain the deviation No
    AcerMC_anglereso_effect_on_correctionfunction2.gif MCNLO_anglereso_effect_on_correctionfunction2.gif

Slide 8: Evidence 1 for not-bias

  • Use MCNLO events, same statistics as "rome" data, the measurement shows no obvious bias.
MCNLO_Measurement.gif

log

open new data...$HOME/analysisResults/forplots_3_jcali_Moore_MCALO_atlfast_btagemulation_jetprecalibration/ANA.ntuple.root
open the 3nd data...$HOME/analysisResults/forplots_3_jcali_Moore_MCALO_btag60_HLT_AODlevelbugfix/ANA.ntuple.root
In CalculateAndSetCorrectionFunction()
 in recons level, events=6104    A=-0.0178893
 in parton level, events=253416  A=0.00158521
Correction function is 1.119092+(0.007044)*x+(-0.004921)*x*x+(-0.000065)*x*x*x+(0.000012)*x*x*x*x
full A:   mean: -0.0238551+/-0.0845253

Slide 9: Evidence 2 for not-bias

  • Use the 340 AcerMC events, where fast and full simulation both reconstructed ttbar events correctly (pure signal).
  • Plot the difference of $-9/0.51*cos\theta1*cos\theta2$ between fast and full simulation, the mean is 0.013 that bias doesn't exist in these events. (0.097 for MCNLO 1675 events)
AcerMC_samepuresignal_fastfull.gif
  • Use the 340 events to repeat the measuring method, got A=0.07. Even larger error from correction function?
AcerMC_Measure_with_bothreco_events.gif

Slide 10: Evidence 3 for not-bias

  • 220/pb AcerMC full simulation samples is split into 2 parts, according to its event number even or odd.
    • Calculate correction function with events of even event number
    • Apply on the events with event number odd.
  • The measured A=-0.84 and 1.10 when exchange the 2 samples . (-0.78 and 1.16 if change correction function bins to 9)
AcerMC_fullsim_2parts_measurement1.gif

AcerMC_fullsim_2parts_measurement2.gif
  • Although no bias in the samples, the deviation is more than $2\sigma$. Do same to MCNLO events, deviation larger than $2\sigma$ is observed too: 0.296882+/-0.118093 and -0.275076+/-0.115798.
  • The correction function is from very low statistics, with large uncertainties around 0.6 and 0.8, which can explain the very large deviation of above. The estimation of the MCNLO sample shows the error are around 0.3 too, to be able to explain the deviation in MCNLO sample.
AcerMC_correctionfunc_error_halffulldata.gif AcerMC_correctionfunc_error_halffulldata2.gif

In CalculateAndSetCorrectionFunction()
 in recons level, events=2523    A=-0.180958
 in parton level, events=253416  A=0.00152187
Correction function is 1.109897+(-0.009322)*x+(-0.004660)*x*x+(0.000038)*x*x*x+(0.000012)*x*x*x*x
full A:   mean: 0.296882+/-0.118093

In CalculateAndSetCorrectionFunction()
 in recons level, events=2534    A=0.0872601
 in parton level, events=253416  A=0.00152187
Correction function is 1.067744+(0.002200)*x+(-0.002326)*x*x+(-0.000007)*x*x*x+(0.000000)*x*x*x*x
full A:   mean: -0.275076+/-0.115798

Slide 11: Cut effect

  • We use quality cut to remove non-well-reconstructed events to improve S/B, use cut on ttbar mass to enhance the asymmetry A.
    • There are obvious difference on these quantities between fast and full simulation. Using same quality cut on both will mean different operation
    • As below, the green is the measurement difference between full simulation and fast simulation (correction function from fast simulation), with the same quality cuts, and the blue use specific quality cut for fast and full respectively.
  • Specific cuts helps to reduce the deviation.

AcerMC sample and MCNLO sample

AcerMC_qualitycut_effects.gif MCNLO_qualitycut_effects.gif

Meaning of number

6 7 8 9 10
no-quality-cut njet=4 pt_e25GeV Mttbar<550G 8+9+pt_n/etmiss>0.99
11 12 13 14
8+9+|Mw-80.4<20 8+9+|Mthad-175|<35 8+9+|Mtlep-175|<35 8+9+10+11+12+13
  • Measure with different quality cut for fast and full simulation: w and top mass cut $2\sigma$ around peak value.
  • A=0.557148+/-0.176043, But can't believe it so much because the best correction is from cut on leptonic top mass, but actually this is the least difference between fast and full.

AcerMC_measuring_diff_qualitycut_fufa.gif

Slide 12: Continue with above

  • With the same selection and quality cuts on the reconstructed quantities, but this time, use $-9/0.51*cos\theta_{1true}*cos\theta_{2true}$ to remove the resolution effect, the discrepancy between full and fast simulation is much smaller, means resolution is part of the reason of the deviation.

AcerMC and MCNLO

AcerMC_qualitycut_effects_truthinfo.gif MCNLO_qualitycut_effects_truthinfo.gif
  • However, we see that
    • The deviation is smaller but still larger than the error bar
    • MCNLO sample: deviation larger than error bar emerges, don't understand! does it mean bias is there introduced by cut and reco procedure effect but is compensated by resolution uncertainty in MCNLO and extended in AcerMC?
  • Take 1300 truth value as measured one (removing reconstruction effects) randomly for 45 times to calculate the correction function and applied on the same sample containing 1192 events of truth value (no reco info, too). The resulted 45 times of measurement (left plot) shows a error 0.005. If we try the $1\sigma$ shift on parameters of correction function, we have 243 times of measurement(right plot), shows a error of 0.35, which is consistent with above: $0.35^2+0.15^2=0.39^2$. The question is which error is reliable? For the first error estimation, 0.005, it's clear that bias is from cuts and reconstruction procedure. for the second, the statistics conceal any possible bias in it.
    AcerMC_onlystatistics1.gif AcerMC_correctionfunc_onlytruthinfo.gif
  • ttbar mass is lower measured in atlfast than full simulation, this makes small bias as 0.07 (farther from truth value). With this cut, we cut 30% of ttbar events, 34%(33%) reconstructed events and 44%(41%) well-reconstructed events for atlfast and full simulation respectively. High-mass ttbar events are easier be reconstructed, may give up this cut.

Slide 13: New method

  • From above:
    • We didn't found any clear sign that a bias exist in the measurement
    • We should find out all potential errors which may resuts in this large deviation
  • A clear sign is statistical fluctuation, especially on the tail of the distribution, contributes much to the deviation, and may be amplified by correction function obtained from other sample.
  • What we learned from the w polarization study is the parameters extracted from fit is not very sensitive to these occasionally happened fluctuations. So we can try to use fit to measure it.
  • Look at the ttbar correlation function
\[ \frac{1}{N}\frac{d^2N}{dcos\theta_1dcos\theta_2}=\frac{1}{4}(1-Ccos\theta_1cos\theta_2)\]
. If we integral on one of cosine of angle, for example $cos\theta_1$, we have
\[ \frac{1}{N}\int{d\frac{dN}{dcos\theta_2}}=(\int^0_{-1}+\int^1_0)\frac{1}{4}(1-Ccos\theta_1cos\theta_2)dcos\theta_1\]
\[ \frac{1}{N}\frac{dN}{dcos\theta_2}=\frac{1}{4}(1+\frac{1}{2}Ccos\theta_2) + \frac{1}{4}(1-\frac{1}{2}Ccos\theta_2) \]
this is not what we want, because it became a constant and we lost information of C. But here we give a tranformation $\theta_2\rightarrow\pi-\theta_2$ when $cos\theta_1>0$ So the integral above became
\[ \frac{1}{N}\frac{dN}{dcos\theta_2}=\frac{1}{4}(1+\frac{1}{2}Ccos\theta_2) + \frac{1}{4}(1+\frac{1}{2}Ccos\theta_2) = \frac{1}{4}(2+ Ccos\theta_2) \]
In the measurement, we can use that tranformation to fill a one-dimension histogram of $cos\theta_2$ and fit it with function "pol1" to extract C. and A=C/0.51
  • Measurement, * code :
    • The Y axis is scaled with (-4) * bins/2 /0.51 after normalization, because first of all, the last formula above is normalized, second, (-4) is the proportional constant; thirdly, the histogram is divided to bins/2 bins each x-axis unit; last, 1/0.51 is the spin analyzer power.
    • integral $cos\theta_2$, of the hadronic decayed top side: A=0.39+/-0.19. The plots are fast truth, fast reco, fast correction function and full measurement.

AcerMC_fitA_measurement5.gif

    • integral $cos\theta_1$, of the leptonic decayed top side: 0.31+/-0.19. The plots are fast correction function and full measurement.

AcerMC_fitA_measurement6.gif

    • test it with MCNLO sample, we got 0.05+/-0.11 and 0.09+/-0.11 for each side integral.
    • Test the sensitivity to "different physics", correction function from AcerMC is used on MCNLO samples, we got 0.17+/-0.11, exchange the samples we have 0.30/-0.18 , it means when the physics changed, we can get the sign
      fitA_MCNLOcorrfunc_AcerMCmeasure.gif fitA_AcerMCcorrfunc_MCNLOmeasure.gif

Slide 14: Use fitting method for Ad code

  • Although the fitting method works well for A measurement, how doest it work for Ad?
  • Ad measurement has been taken as the Mean of 1-D distribution of $cos\Phi$. now use the fit of the theoretical funcition
    \[\frac{1}{N}\frac{dN}{dcos\Phi}=\frac{1}{2}(1-Dcos\Phi) \]
  • We got for AcerMc sample Ad=-0.39+/-0.10 ( Mean method gives -0.45+/-0.10 ), the plot below are fast correction function, full measurement and fast measurement (correct itself). Fitting method works well too.
    AcerMC_fitAd_measurement5.gif
  • Try on the MCNLO sample, Ad=-0.11+/-0.06 (Mean method gives -0.07+/-0.05), the plot below are fast correction function, full measurement and fast measurement (correct itself). The Fitting method is even worse, one fact is the fitting method gives fast measurement 0.08+/-0.05, it means that the error from the correction function is large
    MCNLO_fitAd_measurement5.gif

Slide 15: bug found for the new method.

  • Bug:
    • The distribution of the above is plotted with $cos\theta_1 \frac{cos\theta_2}{|cos\theta_2|}$ and scale factor for both reconstruction and parton levels.
    • But in the correction function application, only $cos\theta_1$ is used, other than $cos\theta_1 \frac{cos\theta_2}{|cos\theta_2|}$
  • After bug fix. the measurement A=0.66, similar as the exist method. (plot 4)
    • AcerMC_fitA_measurement5_bugfix_08_03_24.gif
  • Summary:
    • This method doesn't solve the bias problem, but prove that old method is reliable.
    • The possible reason for bias is explored, see coslep>0 suppressed in full sim.
    • Coshad_lej is slightly affected, so the new method with coshad_lej work better, it is ~0.15 closer to SM value than on coslep.

open new data...$HOME/analysisResults/5205_12605_jcali_atlfast/ANA.ntuple.root
open the 3nd data...$HOME/analysisResults/5205_12605_jcali_fullsim_truthHardProc/ANA.ntuple.root
In CalculateAndSetCorrectionFunction()
Correction function is 0.958208+(-0.196741)*x+(0.183672)*x*x+(0.070694)*x*x*x+(-0.159684)*x*x*x*x
the slope is: -0.391749+/-0.0914739
full A:   mean: -0.446723+/-0.0966451
the slope is: -0.263819+/-0.0863646
fast A:   mean: -0.281153+/-0.0906166
Info in <TCanvas::Print>: GIF file AcerMC_fitAd_measurement5.gif has been created

open new data...$HOME/analysisResults/forplots_3_jcali_Moore_MCALO_atlfast_btagemulation_jetprecalibration/ANA.ntuple.root
open the 3nd data...$HOME/analysisResults/forplots_3_jcali_Moore_MCALO_btag60_HLT_AODlevelbugfix/ANA.ntuple.root
In CalculateAndSetCorrectionFunction()
Correction function is 0.975217+(0.039795)*x+(0.051875)*x*x+(-0.176321)*x*x*x+(0.018757)*x*x*x*x
the slope is: -0.108361+/-0.05298
full A:   mean: -0.0732559+/-0.0472294
the slope is: 0.0759533+/-0.0486825
fast A:   mean: -0.00752075+/-0.0432925
Info in <TCanvas::Print>: GIF file MCNLO_fitAd_measurement5.gif has been created

// the old sample which has duplication problem
open new data...$HOME/analysisResults/5205_jcali_atlfast/ANA.ntuple.root
open the 3nd data...$HOME/analysisResults/5205_forplots_3_jcali_Moore_MCALO_btag60_HLT_truthHardProc/ANA.ntuple.root
In CalculateAndSetCorrectionFunction()
Correction function is 0.908255+(-0.257376)*x+(0.423997)*x*x+(0.301705)*x*x*x+(-0.342519)*x*x*x*x
the slope is: -0.330671+/-0.0951342
full A:   mean: -0.337228+/-0.0999412
the slope is: -0.13286+/-0.0882929
fast A:   mean: -0.15723+/-0.0931127
Info in <TCanvas::Print>: GIF file AcerMC_fitAd_measurement5.gif has been created

Slide 16: Atlfast B tagging effect

  • The tagging weights for b-jets and other jets in ATLFAST events are randomly set, emulating the (IP3D+SV1) weight distribution in full simulation as a function of the pT and eta. Because of this, ∼70% of the selected ATLFAST events and full simulation events are not identical events, although they are read out from same data sample. They are regarded as independent.
  • I re-process once again the data (same data, same program) and produce a new ntuple file. Between this run and old run , they are statistically same if in large scale, now the difference is
    • weight is re-set randomly again to each jet.
    • for individual event, whether it passes n(bjet)>=2 is random.
  • So it is same as a new independent atlfast data is built. To repeat the measurement: run log
  without trigger with trigger (statistical error ~0.19)
method 1: A 0.58 0.45  
mehtod 2: A 0.59 0.49 see plot
method 1: Ad -0.33 -0.28  

  • Question:
    • the atlfast random b tagging are not fully compatible with the b tagging in full simulation? b-tagging in full simulation correlated with some variables?
    • or only statistically contribute large error? but in previous analysis 1, there were enough statistics.
  • Summary
    • Trigger affect the measurement. in t7 notes, trigger is not considered because the effect is much smaller than the bias.
  • Btaging in atlfast is random, means all b jets has the same possibility to be tagged, which is not the case for full data. We compare the truth angle coslep, coehad_lej and their product distribution with btag effect. From the table, we see, that the effect is A_mesur ~= 0.42/0.316*0.394 =0.55, which is biased around 0.13.
A_true njet>=4 njet>=4 && nbjet>=2
fast data 0.307 0.316
full data 0.306 0.394
  • red: fast data, and blue full data
    AcerMC_btag_effect.gif
  • with MC@NLO data (CSC 5200). Question: is the bias propertional to the A?
    MCNLO_measurement.gif

Slide 17: Use Correction function extracted from MCNLO full simulation

  • In MCNLO events, although there is no spin correlation infomation, the physics and detector effect on the objects are same to AcerMC events and .... all other CSC samples.
    • So in principle, we can use the correction function extracted from MCNLO samples and apply on the AcerMC events.
  • Do it without trigger simulation
  SM value measured value $\sigma$
A 0.42 0.31 0.16
Ad -0.29 -0.17 0.10
  • Do it with trigger simulation
  SM value measured value $\sigma$
A 0.42 0.29 0.18
Ad -0.29 -0.12 0.11
  • Extract correction function from AcerMC and MCNLO events:
    Correction_function_fulldata.gif

Slide 18: weighted distribution

  • A is the mean of the distribution of -9/0.51*coslep*coshad_lej. The contribution to A depend both on the value of coslep*coshad_lej and number of events (bin content). So try to weight the bin content with the value of abs(coslep*coshad_lej).
    • AcerMC (red atlfast, blue full sim, left: truth, right: reco)
    • AcerMC_fastfull_contributions.gif

    • MCNLO
    • MCNLO_fastfull_contributions.gif

Slide 19: Summary and Conclusion

  • For sure that cut on ttbar mass make difference between fast and full simulation, but the induced visible bias can't be clearly calculated out because of the low statistics, but expected too small to be responsible for the bias in the analysis.
  • In CSC analysis, the bias can be explained by its large uncertainty from correction function.
    • If the error calculated by the method 1 is more reliable (change $1sigma$ of the parameters of the correction function)
  • In the ROME analysis, we take all possible error 0.09, 0.10, 0.08 from correction function, statistics of full simulation and systematics respectively, which is 0.16 still smaller than the bias 0.25. Concerning the bias induced by cut on ttbar mass, it can be considered as acceptable result.
  • The quality cuts are the possible final source of the bias, reconfiguring quality cut for fast and full simulation respectively, we have a acceptable result A=0.55, but it can't be for sure in the low statistics, either now.
  • With the new measuring procedure, which is not sensitive to the occasionally happened statistical fluctuation, we have good results until now.
  • To do: produce fast simulation data with higher statistics.

Some discoveries

  • ttbar correlation: 0 means flat in cosine distribution of the 2 angles. but combk means flat in angle distribution, so that not flat in cosine plne, but symmetry

Consideration 4 for non-bias

  • The signal and the background mix with a fixed ratio in fast and full simulation in AcerMC 73%(576/788) and 74.8%(510/682). But not the same ratio in MCNLO 113.8%(2364./2078) and 105.3% (1883/1788)
  • With MCNLO sample, the background effect doesn't make bias. The reason is background (combk and non-ttbar events) doesn't have correlation information and should have the same shape as MCNLO sample, so the different ratio of S/B should not affect
  • A similar S/B in AcerMC should not bring effects from background.

Bias on measurement of A code

  • As seen, using the correction function obtained in atlfast to apply on full simulation, the measurement on A is biased of around $2\sigma$.
  • There are 3 type of calculation:
    • dataI, based on truth infomation: no detector and cut effects
    • dataII, based on atlfast measurement: ideal detector and cuts effects
    • dataIII, based on atlfull measurement: realistic detector and cuts effects
  • How to understand the bias and what makes the bias.
    • cuts introduced
    • resolution introduced
    • statistics
  • Considerations
    • correlation between 2 measurements needs more statistics?
    • other method for measuring?
    • at least using correction function from full simulation for real measurement
    • why effects are larger for product of 2 variables?
    • correlation measurement need how good a resolution? When the resolution is worsen than a limit, the correlation will not be measured? Can I have $\sigma_{A}(\sigma_{cos\theta1*cos\theta2})$
    • If a large $\sigma_{cos\theta1*cos\theta2}$ need a large statistics for precise measurement of A. Can I have have $N(\sigma_{A})$
  • Checks and conclusions
    • divide the atlfast data into 2 part. Use the correction function got from part 1 on part 2. Check the bias, if no bias, means not statistical problem; if bias, means statistical problem.
the truth A=0.405,  error 0.03 consistent with SM  (with Mttba<550G)
use data "eventnumber even" to get the correction function and measurement with data "eventnumber odd":
     A=0.48, error 0.23     (change to quality cut specially for atlfast,  A=0.29, error 0.24 )
change the 2 parts of data:
     A=0.43, error 0.22    
means the correction function inside atlfast works well, the bias is not from statisticas only.
    • the same can be done to atlfull
It's fairly bad, before correction, "eventnumber odd" A=-0.01,  "eventnumber even" A=0.51
use data "eventnumber even" to get the correction function and measurement with data "eventnumber odd": 
      A=-0.84, error 0.23
change the 2 parts of data:
      A=1.10 error 0.23 
Means statistical problem in full simulation data.....?
More checks is done to distribution of reconstructed -9/0.51*cos\theta1*cos\theta2 for 2 parts of atlfull data directly without any correction.
The means are -0.02 and 0.51 respectively, error 0.21 for both.
If checks is done to their corresponding truth: 
The means are 0.46 and 0.37 respectively, error are similar.
That means: for the small statistics, .
    • Use dataI , for the truth cuts and atlfast reconstructed cuts, to check whether the resolution bias the cuts much. If yes, change to atlfull reconstructed cuts to see whether larger bias?
    • Artificially put different resolution to truth (dataI), and different statistics, to see whether the correlation disappear without cuts and detector effects? how the angle resolution propogate to the final measurement, and makes needing more statistics. How the resolution bring errors on the final measurement? the resolution may be affect much collaberating with correction function.
    • The correction function from atlfast may bring large uncertainty due to low statistics. The fact is the old correction function is from old atlfast with good statistics, means small uncertainty due to statistics. So the bias should not be from atlfast statistics.
    • The atlfull data may have different composition in least energy jet (different ratio of quarks flavors), so different spin analyzer power , but the fact is there is no effect on Ad, so the effect should be not the reason of the bias.
    • The fact that there is a slope on the ratio between the distribution of full and fast for different data, means the background is not the reason of the bias. but the slope is the reason, what makes the slope?
    • the e event and muon events ratio are different for atlfast and atlfull. As known the e event and muon event have different results(not know why), this may be the reason of the bias. also atlfull has more jets means more combk?
this is not the reason. 
1. we give a weight to e and mu events, to adjust them to the same ratio in atlfast as in atlfull. the measure still gives in atlfull A=0.84.  
2. We have tried to correct efficiency function using the fact that reco efficiency for atlfast and atlfull are different, and the correction doesn't work at all to A.
* atlfull has more jets means more combk? * should use different quality cut for atlfast and atlfull, to decrease the bias, because the event reonstruction (mass distribution) are different for them

Use MCNLO data,

1. which has no ttbar spin correlation, to see the $cos\theta_1$ and $cos\theta_2$ in 2D histogram, which has 2x2 bins. The results is the bin closed to (-1,-1) has events 45200 and the other 3 bins have 45900 events each.

2. To look at $\langle cos\theta_1 cos\theta_2 \rangle$, which is the mean of the distribution, in truth, it gives 0.005.

3. Use the correction function obtained from atlfast in the reconstruction level, the bias is 0.26.

4. If I take the cuts officially defined for T7 note and corresponding correction function, the bias is 0.36. The bias is same as study with AcerMC data.

  • MCNLO data (T7CSCnot selection) corrected with AcerMC atlfast corr_func red atlfast:
    mcnlo12_redfast_bluefull_AcerMC12fastcorrection.gif

  • correction with 2D efficiency::
    MCNLO_biasA_tt_correlation_A_corrected_with_2Defficiency.gif

  • correction with from 2 angle efficiencies:
    MCNLO_biasA_tt_correlation_A_corrected_twice_with_2angleefficiency.gif

  • efficiencies of the 2 anlges:
    MCNLO_biasA_fastfull_efficiency_2angles.gif

  • 2D efficiency for full and fast:
    MCNLO_biasA_fastfull_efficiency_2D.gif

  • resoultion effects on distribution of angle 1D:
    MCNLO_true_resolution_effect1.gif

  • resoultion effects on distribution of angle 2D:
    MCNLO_true_resolution_effect2.gif

Use AcerMC data:

* AcerMC data (selection T7CSCnote) corrected with it's corr_func from atlfast, red atlfast:
AcerMC12_redfast_bluefull_AcerMC12fastcorrection.gif

* to look at the ratio of distribution of $cos\theta_1 cos\theta_2 $ between full sim data and atlfast data. the data used are 3 types: semiemu; semiemu*goodlnb; semiemu*goodlnb*goodjjb and the distribution of ratios are similar and the bias on measurement of A is even larger. It means the background (tau events and cmbk) is not the reason of the bias.

    • full over fast for costheta1costheta2 for 3 different data : semiemu. semiemu*goodlng, semiemu*goodlnb*goodjjb:
      fullfast_full_over_fast_costheta1costheta2_no_correction.gif

    • without quality cut, full over fast for costheta1costheta2 for 3 different data : semiemu. semiemu*goodlng, semiemu*goodlnb*goodjjb::
      fullfast_full_over_fast_costheta1costheta2_no_correction_noqualitycut.gif
  • for truth $cos\theta_1 cos\theta_2 $, the events are reserved only if on the reconstruction level, the event is reconstructed. This is to remove the resolution effect. we see that
    • the bias is very small <0.05
    • if we cut the \|Mwjj-Mw\|<20000, we got bias of 0.15. To change Mw to as measured 80.5 for atlfast and 83.5 for full sim, then the bias is smaller as 0.08. which means the cut result in different events for atlfast ans atlfull, resulting in bias.
    • if we cut the \|Mlnb-Mt\|<35000, we got bias of 0.12. To change Mt as measured 172 for atlfast and 173.3 for atlfull. bias doesn't change.
    • if we cut the \|Mjjb-Mt\|<35000, we got bias of 0.07. To change Mt as measured 172 for atlfast and 177 for atlfull. bias is 0.06
    • if we cut the pt_n/pt_miss>0.99, we got bias of 0.26. if cut pt_n/pt_miss>0.9, bias 0.29, if cut pt_n/pt_miss>0.6, bias 0.03. This cut makes big difference as in the following plot
    • et shrink when reconstructing wln:
      fullfast_et_shrink_in_wln_reco.gif

  • cut effect to costheta1*costheta2:
    biasA_cut_effects_and_correction_recolevel.gif

  • cut effects on costheta1_true*costheta2_true:
    biasA_cut_effects_and_correction_partonlevel.gif

  • correction with 2D efficiency:
    AcerMC_biasA_tt_correlation_A_corrected_with_2Defficiency.gif

  • correction with from 2 angle efficiencies:
    AcerMC_biasA_tt_correlation_A_corrected_twice_with_2angleefficiency.gif

  • efficiencies of the 2 anlges:
    AcerMC_biasA_fastfull_efficiency_2angles.gif

  • 2D efficiency for full and fast:
    AcerMC_biasA_fastfull_efficiency_2D.gif

  • resoultion effects on distribution of angle 1D:
    AcerMC_true_resolution_effect1.gif

  • resoultion effects on distribution of angle 2D:
    AcerMC_true_resolution_effect2.gif

Topic attachments
I Attachment History Action Size Date Who Comment
Unknown file formatcxx ABIAS.cxx r2 r1 manage 60.8 K 2008-01-08 - 14:11 ChengguangZhu  
Unknown file formatcxx Abias.cxx r1 manage 56.8 K 2007-12-23 - 14:43 ChengguangZhu code
GIFgif AcerMC12_redfast_bluefull_AcerMC12fastcorrection.gif r1 manage 14.9 K 2007-12-17 - 19:42 ChengguangZhu AcerMC data (selection T7CSCnote) corrected with it's corr_func from atlfast
GIFgif AcerMC_Correction_function2_9_bins.gif r4 r3 r2 r1 manage 10.3 K 2008-04-06 - 19:36 ChengguangZhu  
GIFgif AcerMC_Correction_function_13_bins.gif r2 r1 manage 6.7 K 2008-04-06 - 02:16 ChengguangZhu  
GIFgif AcerMC_Correction_function_19_bins.gif r2 r1 manage 7.7 K 2008-04-06 - 02:16 ChengguangZhu  
GIFgif AcerMC_Correction_function_25_bins.gif r2 r1 manage 8.3 K 2008-04-06 - 02:17 ChengguangZhu  
GIFgif AcerMC_Correction_function_6_bins.gif r2 r1 manage 5.6 K 2008-04-06 - 02:14 ChengguangZhu  
GIFgif AcerMC_Correction_function_7_bins.gif r2 r1 manage 5.8 K 2008-04-06 - 02:15 ChengguangZhu  
GIFgif AcerMC_Correction_function_9_bins.gif r2 r1 manage 6.1 K 2008-04-06 - 02:15 ChengguangZhu  
GIFgif AcerMC_Measure_with_bothreco_events.gif r1 manage 9.1 K 2008-01-04 - 19:36 ChengguangZhu  
GIFgif AcerMC_anglereso_effect_on_correctionfunction2.gif r1 manage 5.7 K 2008-01-07 - 11:05 ChengguangZhu  
GIFgif AcerMC_biasA_fastfull_efficiency_2D.gif r1 manage 9.1 K 2007-12-23 - 19:57 ChengguangZhu 2D efficiency for full and fast
GIFgif AcerMC_biasA_fastfull_efficiency_2angles.gif r2 r1 manage 12.1 K 2007-12-23 - 20:09 ChengguangZhu  
GIFgif AcerMC_biasA_tt_correlation_A_corrected_twice_with_2angleefficiency.gif r1 manage 38.6 K 2007-12-23 - 19:56 ChengguangZhu correction with from 2 angle efficiencies
GIFgif AcerMC_biasA_tt_correlation_A_corrected_with_2Defficiency.gif r1 manage 42.1 K 2007-12-23 - 19:55 ChengguangZhu correction with 2D efficiency
GIFgif AcerMC_btag_effect.gif r1 manage 13.6 K 2008-04-07 - 01:22 ChengguangZhu  
GIFgif AcerMC_correctionfunc_error.gif r2 r1 manage 6.9 K 2008-01-07 - 13:28 ChengguangZhu  
GIFgif AcerMC_correctionfunc_error_halffulldata.gif r1 manage 5.2 K 2008-01-08 - 18:41 ChengguangZhu  
GIFgif AcerMC_correctionfunc_error_halffulldata2.gif r1 manage 5.4 K 2008-01-08 - 18:41 ChengguangZhu  
GIFgif AcerMC_correctionfunc_onlytruthinfo.gif r3 r2 r1 manage 5.7 K 2008-01-09 - 15:47 ChengguangZhu  
GIFgif AcerMC_estimate_error_from_correctiofunction2.gif r1 manage 7.0 K 2008-01-09 - 18:13 ChengguangZhu  
GIFgif AcerMC_fastfull_contributions.gif r1 manage 12.2 K 2008-04-16 - 23:04 ChengguangZhu  
GIFgif AcerMC_fitA_measurement5.gif r1 manage 23.0 K 2008-01-08 - 11:07 ChengguangZhu  
GIFgif AcerMC_fitA_measurement5_bugfix_08_03_24.gif r1 manage 23.6 K 2008-03-24 - 16:50 ChengguangZhu  
GIFgif AcerMC_fitA_measurement6.gif r2 r1 manage 12.1 K 2008-01-08 - 11:22 ChengguangZhu  
GIFgif AcerMC_fitAd_measurement5.gif r1 manage 18.1 K 2008-01-15 - 11:54 ChengguangZhu  
GIFgif AcerMC_fitAd_measurement52.gif r1 manage 22.0 K 2008-04-04 - 17:20 ChengguangZhu  
GIFgif AcerMC_fullsim_2parts_measurement1.gif r1 manage 9.5 K 2008-01-04 - 19:55 ChengguangZhu  
GIFgif AcerMC_fullsim_2parts_measurement2.gif r1 manage 8.0 K 2008-01-04 - 19:55 ChengguangZhu  
GIFgif AcerMC_measuring.gif r1 manage 8.3 K 2008-01-06 - 18:36 ChengguangZhu  
GIFgif AcerMC_measuring_diff_qualitycut_fufa.gif r1 manage 8.1 K 2008-01-06 - 20:28 ChengguangZhu  
GIFgif AcerMC_onlystatistics1.gif r1 manage 5.4 K 2008-01-09 - 15:26 ChengguangZhu  
GIFgif AcerMC_qualitycut_effects.gif r1 manage 5.3 K 2008-01-05 - 18:45 ChengguangZhu  
GIFgif AcerMC_qualitycut_effects_truthinfo.gif r1 manage 5.0 K 2008-01-05 - 20:25 ChengguangZhu  
GIFgif AcerMC_resolution.gif r1 manage 11.1 K 2008-01-05 - 15:36 ChengguangZhu  
GIFgif AcerMC_samepuresignal_fastfull.gif r2 r1 manage 4.1 K 2008-01-04 - 19:19 ChengguangZhu  
GIFgif AcerMC_true_resolution_effect1.gif r1 manage 41.5 K 2007-12-24 - 18:47 ChengguangZhu resoultion effects on distribution of angle 1D
GIFgif AcerMC_true_resolution_effect2.gif r1 manage 59.0 K 2007-12-24 - 18:47 ChengguangZhu resoultion effects on distribution of angle 2D
GIFgif Correction_function_fulldata.gif r1 manage 5.0 K 2008-04-14 - 00:32 ChengguangZhu  
Unknown file formatcxx ImA.cxx r3 r2 r1 manage 30.7 K 2008-03-24 - 16:53 ChengguangZhu  
Unknown file formatcxx ImAd.cxx r1 manage 21.8 K 2008-01-15 - 12:27 ChengguangZhu  
GIFgif MCNLO_Measurement.gif r1 manage 9.6 K 2008-01-04 - 18:36 ChengguangZhu  
GIFgif MCNLO_anglereso_effect_on_correctionfunction2.gif r1 manage 5.7 K 2008-01-08 - 10:57 ChengguangZhu  
GIFgif MCNLO_biasA_fastfull_efficiency_2D.gif r1 manage 9.3 K 2007-12-23 - 20:03 ChengguangZhu 2D efficiency for full and fast
GIFgif MCNLO_biasA_fastfull_efficiency_2angles.gif r2 r1 manage 11.9 K 2007-12-23 - 20:09 ChengguangZhu  
GIFgif MCNLO_biasA_tt_correlation_A_corrected_twice_with_2angleefficiency.gif r1 manage 37.6 K 2007-12-23 - 20:00 ChengguangZhu correction with from 2 angle efficiencies
GIFgif MCNLO_biasA_tt_correlation_A_corrected_with_2Defficiency.gif r1 manage 43.1 K 2007-12-23 - 19:59 ChengguangZhu correction with 2D efficiency:
GIFgif MCNLO_correctionfunc_error_oldfast.gif r1 manage 6.1 K 2008-01-08 - 13:07 ChengguangZhu  
GIFgif MCNLO_fastfull_contributions.gif r1 manage 13.2 K 2008-04-16 - 23:05 ChengguangZhu  
GIFgif MCNLO_fitAd_measurement5.gif r1 manage 18.1 K 2008-01-15 - 11:54 ChengguangZhu  
GIFgif MCNLO_measurement.gif r1 manage 14.1 K 2008-04-08 - 20:37 ChengguangZhu  
GIFgif MCNLO_qualitycut_effects.gif r1 manage 5.2 K 2008-01-05 - 18:46 ChengguangZhu  
GIFgif MCNLO_qualitycut_effects_truthinfo.gif r1 manage 5.3 K 2008-01-05 - 20:25 ChengguangZhu  
GIFgif MCNLO_true_resolution_effect1.gif r1 manage 40.2 K 2007-12-24 - 18:48 ChengguangZhu resoultion effects on distribution of angle 1D
GIFgif MCNLO_true_resolution_effect2.gif r1 manage 60.4 K 2007-12-24 - 18:53 ChengguangZhu resoultion effects on distribution of angle 2D
GIFgif biasA_cut_effects_and_correction_partonlevel.gif r2 r1 manage 9.3 K 2007-12-23 - 14:59 ChengguangZhu  
GIFgif biasA_cut_effects_and_correction_recolevel.gif r2 r1 manage 10.0 K 2007-12-23 - 14:59 ChengguangZhu  
Unknown file formatcxx com_correction_function.cxx r3 r2 r1 manage 18.8 K 2008-04-06 - 19:37 ChengguangZhu  
GIFgif fitA_AcerMCcorrfunc_MCNLOmeasure.gif r1 manage 6.2 K 2008-01-09 - 13:18 ChengguangZhu  
GIFgif fitA_MCNLOcorrfunc_AcerMCmeasure.gif r1 manage 6.3 K 2008-01-09 - 13:17 ChengguangZhu  
GIFgif fullfast_et_shrink_in_wln_reco.gif r2 r1 manage 4.6 K 2007-12-19 - 14:43 ChengguangZhu  
GIFgif fullfast_full_over_fast_costheta1costheta2_no_correction.gif r1 manage 7.8 K 2007-12-19 - 09:25 ChengguangZhu full over fast for costheta1costheta2 for 3 different data : semiemu. semiemu*goodlng, semiemu*goodlnb*goodjjb
GIFgif fullfast_full_over_fast_costheta1costheta2_no_correction_noqualitycut.gif r1 manage 7.8 K 2007-12-19 - 11:02 ChengguangZhu without quality cut, full over fast for costheta1costheta2 for 3 different data : semiemu. semiemu*goodlng, semiemu*goodlnb*goodjjb:
GIFgif inversed_correctionfunc_black_oldfast.gif r1 manage 5.3 K 2008-01-08 - 13:08 ChengguangZhu  
GIFgif mcnlo12_redfast_bluefull_AcerMC12fastcorrection.gif r1 manage 14.7 K 2007-12-17 - 19:41 ChengguangZhu MCNLO data (T7CSCnot selection) corrected with AcerMC atlfast corr_func
Texttxt runlog1.txt r1 manage 1.2 K 2008-04-04 - 18:02 ChengguangZhu  
Edit | Attach | Watch | Print version | History: r58 < r57 < r56 < r55 < r54 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r58 - 2008-04-16 - ChengguangZhu
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

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