Changes between Version 1 and Version 2 of LowMassSMHiggs
- Timestamp:
- Apr 6, 2012, 4:33:02 PM (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
LowMassSMHiggs
v1 v2 7 7 We want to investigate the following signatures: 8 8 9 WH - > ell+- + 2b<br />ZH -> ell+ ell- + 2b<br />tth ->ell+ ell- 4b9 WH -> ell+- + 2b[[br]]ZH -> ell+ ell- + 2b[[br]]tth -> ell+ ell- 4b 10 10 11 11 <br /> We are now focusing on ttH 12 12 === II. Goal === 13 13 14 The project is to investigate the significance that can be achieved in the search for a Higgs produced in association with top quark pair <br />considering the '''dileptonic channel'''.<br />The main background is %$ t \bar t $% + 2 jets.<br /><br />The significance can be studied in different scenarios, for example<br />one can start to fill the following table14 The project is to investigate the significance that can be achieved in the search for a Higgs produced in association with top quark pair[[br]]considering the '''dileptonic channel'''.[[br]]The main background is $ t \bar t $ + 2 jets.[[br]][[br]]The significance can be studied in different scenarios, for example[[br]]one can start to fill the following table 15 15 16 16 || Statistical significance || "thin" TF || "Broad TF" || … … 18 18 || events with no ISR || S_3 || S_4 || 19 19 20 <br />By estimating S_1,S_2 ... we will provide a reasonable estimate <br />of the maximum significance that can be reached at the LHC, and will<br />also show what are the most important factors controlling this maximum significance.20 <br />By estimating S_1,S_2 ... we will provide a reasonable estimate[[br]]of the maximum significance that can be reached at the LHC, and will[[br]]also show what are the most important factors controlling this maximum significance. 21 21 22 22 The strategy is to organize the work into a validation procedure (first step) and a pheno study (second step). 23 23 === III. Validation procedure === 24 24 25 Idea: apply MEM to reconstruct the mass of the Higgs or reconstruct the fraction of signal <br />events <strong>UNDER CONTROLLED CONDITIONS.</strong><br /><br />"under controlled conditions" means that the events in the prepared samples follow EXACTLY the probability <br />distribution function that is used to evaluate the matrix element weights.<br /><br />In this way, one can set up the calculation of the weights and check that <br />NO BIAS is observed in the final result, and hence validate the whole procedure in the absence of systematic errors. <br />The procedure for the calculation of the weights should be the same as the one used later on for the second step <br />(i.e. one should consider a finite resolution on jet energies, correct for ISR if necessary, ...)<br /><br />Only the samples of events are prepared in an artificial way so that we control exactly how the events are distributed in phase-space.<br />In particular: <br />- the energy of the final-state partons are smeared exaclty according to the shape of the transfer functions<br />- the effect of ISR (if taken into account) is to boost the events in the transverse plan, according to a known distribution in pT25 Idea: apply MEM to reconstruct the mass of the Higgs or reconstruct the fraction of signal[[br]]events **g>UNDER CONTROLLED CONDITIONS.**[[br]][[br]]"under controlled conditions" means that the events in the prepared samples follow EXACTLY the probability [[br]]distribution function that is used to evaluate the matrix element weights.[[br]][[br]]In this way, one can set up the calculation of the weights and check that [[br]]NO BIAS is observed in the final result, and hence validate the whole procedure in the absence of systematic errors. [[br]]The procedure for the calculation of the weights should be the same as the one used later on for the second step [[br]](i.e. one should consider a finite resolution on jet energies, correct for ISR if necessary, ...)[[br]][[br]]Only the samples of events are prepared in an artificial way so that we control exactly how the events are distributed in phase-space.[[br]]In particular: [[br]]- the energy of the final-state partons are smeared exaclty according to the shape of the transfer functions[[br]] - the effect of ISR (if taken into account) is to boost the events in the transverse plan, according to a known distribution in pT 26 26 27 27 The idea is that once this procedure is validated, it can be used in a reliable way for the pheno study. 28 28 29 The subsection below gives the <strong> work plan</strong>for the validation procedure29 The subsection below gives the**g> work plan** for the validation procedure 30 30 ===== A. Reconstruction of m_H (signal events only) ===== 31 31 ====== 1. Parton level + infinite resolution (DONE) ====== 32 32 33 Reconstruction of the mass of the Higgs using a pure sample of parton-level signal events, and considering a narrow transfer function <br />for jet energies <br /><br />This is done (Priscilla): there is no bias in the reconstructed mass of the Higgs ->OK33 Reconstruction of the mass of the Higgs using a pure sample of parton-level signal events, and considering a narrow transfer function[[br]]for jet energies [[br]][[br]]This is done (Priscilla): there is no bias in the reconstructed mass of the Higgs -> OK 34 34 35 35 ====== 2. Parton level + finite resolution (DONE) ====== 36 36 37 generation of a parton-level event sample (no showering), smearing of the partons energies according to a "broad" transfer <br />function, reconstruction of the mass of the Higgs with the same TF that has been used to smear the parton energies.\\<br /><br />Things to keep in mind:38 * when we smear the energy of the partons, we are forced to apply some cuts on jet energies <br />->one should include the acceptance term in the likelihood.39 * to save some time, we are now considering only the gluon-gluon initiated process (for this <br />sanity check phase, it is ok). But then we should also consider '''only''' the gluon-gluon initiated process<br />when generating the parton-level events. Otherwise we may introduce a bias.40 * one should also check the convergence for the evalutation of the matrix element weights in the regime where <br />the resolution on b-jet energies is much worse than the width of the Higgs.<br />This is a delicate point, because the Higgs decay process %$ H \rightarrow b \bar b$% is overconstrained.<br />By default madweight consider 2 integration channels:37 generation of a parton-level event sample (no showering), smearing of the partons energies according to a "broad" transfer[[br]]function, reconstruction of the mass of the Higgs with the same TF that has been used to smear the parton energies.\\[[br]][[br]]Things to keep in mind: 38 * when we smear the energy of the partons, we are forced to apply some cuts on jet energies[[br]]-> one should include the acceptance term in the likelihood. 39 * to save some time, we are now considering only the gluon-gluon initiated process (for this[[br]]sanity check phase, it is ok). But then we should also consider '''only''' the gluon-gluon initiated process[[br]]when generating the parton-level events. Otherwise we may introduce a bias. 40 * one should also check the convergence for the evalutation of the matrix element weights in the regime where[[br]]the resolution on b-jet energies is much worse than the width of the Higgs.[[br]]This is a delicate point, because the Higgs decay process $ H \rightarrow b \bar b$ is overconstrained.[[br]] By default madweight consider 2 integration channels: 41 41 1. in channel 1, the invariant mass of the Higgs is mapped onto one variable of integration, 42 1. in channel 2, the energies of the b-quarks originating from the Higgs are mapped onto 2 variables of integration, <br />but the invariant mass of the Higgs is not.42 1. in channel 2, the energies of the b-quarks originating from the Higgs are mapped onto 2 variables of integration,[[br]] but the invariant mass of the Higgs is not. 43 43 44 44 When the width of the Higgs is orders of magnitude smaller than the resolution in jet energies, we expect that channel 1 is the most appropriate. This is indeed the case. But I also observed that running with \textbf{only} channel 1 makes a big difference: when I compare the values of the weights calculated with the one-channel integrator and with the two-channel integrator, the weights are systematically underestimated in the second case. The difference is quite sizable when we look at the likelihood: the difference is roughtly 4 units of $Log(L)$. So I would suggest to run MadWeight with only one channel of integration (the first one). This can be done by copying the files main_code_one_channel.f and data_one_channel.inc (available in the drop box) in the MW_P1_gg_bbxbmu+vmbxmu-vmx. … … 63 63 === IV. Pheno study === 64 64 65 <br /> - Redo the analysis, but with samples of events that are as realistic as possible. <br /><br />- Evaluate all the systematic uncertainties.65 <br /> - Redo the analysis, but with samples of events that are as realistic as possible.[[br]][[br]] - Evaluate all the systematic uncertainties. 66 66 === V. Inputs of the analysis === 67 67 68 <br /> <br />We will do the analysis for the LHC at '''14 TeV'''.<br />There are several input parameters that need to be fixed right now.<br /><br />Even during the validation procedure, it will be very useful if we <br />consider realistic values the parameters associated with the <br />final-state cuts, the reconstruction efficiencies, <br />the b-taggings and the energy resolution.<br />In such a way, the significance that we will obtain<br />at the end of the validation procedure will not be completely unrealistic,<br />and this will give us some insights to jump into the second part<br />(e.g. if we find that the significance is extremely low for a given <br />signature even under ideal conditions, it may not be worth to push <br />the analysis further for this signature.)68 <br />[[br]]We will do the analysis for the LHC at '''14 TeV'''.[[br]]There are several input parameters that need to be fixed right now.[[br]][[br]]Even during the validation procedure, it will be very useful if we [[br]]consider realistic values the parameters associated with the [[br]]final-state cuts, the reconstruction efficiencies, [[br]]the b-taggings and the energy resolution.[[br]]In such a way, the significance that we will obtain[[br]]at the end of the validation procedure will not be completely unrealistic,[[br]]and this will give us some insights to jump into the second part[[br]](e.g. if we find that the significance is extremely low for a given [[br]]signature even under ideal conditions, it may not be worth to push [[br]]the analysis further for this signature.) 69 69 70 70 For the theoretical parameters, we can stick to the default param_card.dat file on the web. 71 71 ===== A. Cuts ===== 72 72 73 <br />We need to agree on a set of cuts to be applied on the jets and on the leptons. <br />I think a resonable set of cuts are (see http://arxiv.org/pdf/1106.0902.pdf)73 <br />We need to agree on a set of cuts to be applied on the jets and on the leptons.[[br]]I think a resonable set of cuts are (see http://arxiv.org/pdf/1106.0902.pdf) 74 74 75 pT(jets) > 30 GeV, ||eta(jet)||<2.4 delta R (p_i,p_j) >0.3 with p_i, p_j =jet or lepton75 pT(jets) > 30 GeV, ||eta(jet)||<2.4 delta R (p_i,p_j) > 0.3 with p_i, p_j =jet or lepton 76 76 77 pT(e) > 20 GeV, ||eta(e)||<2.5, pT(mu)> 30 GeV, ||eta(mu)||<2.577 pT(e) > 20 GeV, ||eta(e)||<2.5, pT(mu)> 30 GeV, ||eta(mu)||< 2.5 78 78 79 <br /> <br /><em>Parton-level cuts vs reconstructed-level cuts</em>:<br /><br />In the validation procedure, parton-level cuts are different from reconstructed-level cuts<br />because:<br /> - parton-level events are boosted in the transverse plan (if ISR is taken into account) <br /> - final-state parton energies are smeared according to the shape of the transfer function<br /><br />So one need to apply looser cuts at the parton-level, and then apply the correct set of cuts <br />at the "reconstructed level".79 <br />[[br]]<em>Parton-level cuts vs reconstructed-level cuts</em>:[[br]][[br]]In the validation procedure, parton-level cuts are different from reconstructed-level cuts[[br]]because:[[br]] - parton-level events are boosted in the transverse plan (if ISR is taken into account) [[br]] - final-state parton energies are smeared according to the shape of the transfer function[[br]][[br]]So one need to apply looser cuts at the parton-level, and then apply the correct set of cuts [[br]]at the "reconstructed level". 80 80 ===== B. Transfer function ===== 81 81 82 <br />For the parametrization of the transfer functions, <br />we can stick to the usual asumptions: a superposition<br />of two Gaussian distributions for the energy of the jets, <br />a delta function for all other visible quantities.82 <br />For the parametrization of the transfer functions,[[br]]we can stick to the usual asumptions: a superposition[[br]]of two Gaussian distributions for the energy of the jets, [[br]]a delta function for all other visible quantities. 83 83 84 84 The parametrization of the TF for jet energies is given by … … 86 86 87 87 88 with %$ \delta=E_p-E_j $% (parton-level energy minus reconstructed energy).It would be good to choose values for the parameters<br /> The parameters %$p_i$% can be assumed to depend linearly on the '''parton-level''' energy (%$ p_i=a_1+b_i*E_p $%).88 with $ \delta=E_p-E_j $ (parton-level energy minus reconstructed energy).It would be good to choose values for the parameters[[br]] The parameters $p_i$ can be assumed to depend linearly on the '''parton-level''' energy ($ p_i=a_1+b_i*E_p $). 89 89 90 90 || || '''a_i''' || '''b_i''' || … … 95 95 || p_5 || XXXX || XXXX || 96 96 97 It would be good to choose values for the parameters %$a_i, b_i$% in the TF that capture the typical resolution of the <br />CMS detector. Olivier, do you think you could get these values ?97 It would be good to choose values for the parameters $a_i, b_i$ in the TF that capture the typical resolution of the [[br]]CMS detector. Olivier, do you think you could get these values ? 98 98 99 99 ANSWER from Olivier:
