Multiplicity Status

Aus HERMESwiki
Zur Navigation springen Zur Suche springen

Page maintainer: Sylvester and Achim

Workingpage.jpg This is a working page which will continuously change. The content is not reviewed and is intended for people working on this specific topic.

DC19 Status and TODO

Charged Hadron Multiplicities

VMD

  • Cross check VMD fractions.
  • Decide whether we want a real correction, or just to quote the fractions.

eca: I think we definetely want like we always did show that if we do a correction how things change this gives information on what underlying subprocess is dominating

SJ: I agree, and we certainly plan to commit both corrected and uncorrected multiplicity sets to the database, so the fitter can choose what dataset to use. This is more about what dataset we prefer to show as center values in the actual publication - so not the most important issue to discuss right now :-).

Trigger efficiencies

  • Cross check trigger efficiency correction (SJ, in progress)
  • Using efficiencies for different event topologies.

Binnings

  • Put an overview of all current binnings and integrated binnings here
  • We need one additional binning for the Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle F_z^\text{sea}} and Fehler beim Parsen (Konvertierungsfehler. Der Server („https://wikimedia.org/api/rest_“) hat berichtet: „Cannot get mml. Server problem.“): {\displaystyle F_{z}^{\text{val}}} plots.

Year dependence

  • Rerun smearing matrices using MC with specific 2000, 2004, 2005 geometries
    • Remember, use TMC in 2004 and 2005 MC productions, as we do in data!
  • Introduce same hadron target cuts Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z_\text{vert} < 20} for 2000 as we have to use for 2004 and 2005 because of the TMC.

Finalization

  • Rerun over all new (HTC) productions, use smTrackSave to avoid actually using HTC. (AH)
  • Finalize ordering of all corrections and combination of all data sets. (AH)
  • Cross check all final binnings and integration, using the final setup.

Systematic Studies for Charged Hadron Multiplicities

Unfolding procedure

  • Check effect of finite MC statistics on unfolding.
    • Method: Split a big MC sample in subsamples and extract a smearing matrix from each subsample. Then, unfold the data with the different smearing matrices, and fit a Gaussian distribution in each kinematic bin to the different results. The width of this Gaussian is the systematic uncertainty due to the finite MC statistics. Finally, scale this value comparing the subsample sizes to the actual MC sample used for the unfolding (ideally the union of all used subsamples).
    • Current status: first results with the new big disNG 2005 production. 2D binning versus Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle x} and Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} (integrated over Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle p_{h\perp}} ). These results point to a small contribution to the systematic. (This supersedes the estimate from the older pythia production, that was found to be lacking in statistics in certain kinematic corners).

eca: why is the MC statistics still important, I thought Ed, generated much more statistics so this problems should be gone.

SJ: The new MC sample has 2.5 times the generated events, so that's only a reduction of \sqrt(2.5)\approx 1.58 in systematic uncertainty compared to the old MC. We would really need a lot more if we want to decrease the systematic even more (currently it's of the same order of magnitude as the uncertainty due to the RICH unfolding).

  • Check dependence on the current MC tune: compare current tune with an older tune
  • Do we need another 'model-dependence' (eg. disNG vs PYTHIA) check?

eca: what do you mean by model, pythia and disng have the same fragmentation and most of the underlying process are the same.

SJ: Totally agreed, hence the quotation marks around model-dependence. The only thing we would check with this is the effect of having either the elastic tail or the VMD contribution in the actual smearing matrix. I remember that Alex's results in the past indicated that there wasn't that much change when comparing PYTHIA and disNG unfolding.

RICH

  • Sanity check of PMatrix-derived systematic: (AH, SJ)
    • Analyze the two 2005 data productions with a different RICH background.
    • Check if the resulting multiplicities agree within the PMatrix-derived systematic.
  • Decide whether a PEPSI challenge might be preferable.

Year dependence

  • Use PDG method of unconstrained averaging (sect. 5.5.2 p10). (AH)
    • Test whether the final datasets are statistically compatible
    • If they are not, estimate an additional systematic uncertainty using the weight factor necessary to make the sets compatible.
  • Switch to new productions using new geometry file (00e, 04d,05d)
    • Comparison of 05c1,05d1 (new geometry), 05d2 (new Geom + new RICH BKG)
      • Switch 05c -> 05d reduces the hadron yield (Reason (see page 12): Hadrons only available with z vertex inside target cell (but we used to open the hadron z vertex cut to 100cm, to allow for hadrons from decays)
      • Switch 05d1 -> 05d2 (new RICH BKG): (anti-)proton yield is reduced, pion and kaon yield enhanced
        • expected from lower BKG
        • PMT hits become more significant -> less hadrons are identified as (anti-)protons

Track reconstruction

  • Only accept hadrons in different half from DIS lepton, to avoid possible PID problems of a hadron and lepton track moving close together through the spectrometer.
    • Practically: Split the detector in 2 experiments: one with leptons in top and hadrons in bottom, and one where this is reversed. Pass this through the entire analysis chain.
    • Compare the results with what we obtain using the standard configuration.
    • Also: compare top and bottom :-)

eca: very very bad idea, this will select different kinematics therefore you compare things not very easily. Also which problem are you trying to address with this. I think it will more confuse than help to understand things.

SJ: Gunar and Klaus asked us to look into this. I believe the idea is to check whether we have an influence of the fact that the calo PID fails when the hadron and scattered lepton are to close to each other. I agree with your observation that this will severely alter the average event kinematics and I believe a better way is to compare our current result with what we get when using Rebecca and Francesca's method to fix this issue, explained in the links in the PIDLIB paragraph (by discarding the calo PID in events where the calo fails with code -999).

Thoughts

  • Possible to add a particle dimension to smearing matrix? This would allow for a more realistic smearing between particle types and kinematic bins simultaneously, and remove the problem of the RICH response not being constant versus one of the kinematic variables considered in the analysis but not in the P-matrix. Technical difficulties (memory and CPU time) might make this impossible or at least very "expensive". A feasibility study, and if possible, a comparison with standard RICH unfolded results would be interesting. Did anyone already attempt to do this in the past?
  • Possible to extend multidimensional binning to more relevant variables? We always speak about our limiting statistics making it impossible to bin in additional SIDIS variables like Fehler beim Parsen (Konvertierungsfehler. Der Server („https://wikimedia.org/api/rest_“) hat berichtet: „Cannot get mml. Server problem.“): {\displaystyle \phi _{h}} . Is this really true, or are the MC statistics the limiting factor. This should not be the case!

eca: francesca and rebecca have done the best job people can do on this, and they had severe problems with empty bins, which is what we will have as well. Also we publish multiplicities, this means we have to deal with different holes for the semi-inclusive and the inclusive. If we want to do this we should publish unpolarizied sidis cross sections.

SJ: Again, agreed, both points are just purely speculation.

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} Multiplicities

  • Finish cross check of unfolded Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} multiplicities.
    • Current status: disagreement in unfolded multiplicities due to different integration method and error estimation. Slightly differing methods are cause of these disagreements:
      • Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} reconstruction by fit to invariant mass histogram, see the multiplicity comparison plots before unfolding. Add plot after unfolding
        • AH: numerical integration of gaussian over finite range by subtracting background fit from invariant mass histogram
        • SJ: analytical integration of gaussian fit from Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle -\infty} to Fehler beim Parsen (Konvertierungsfehler. Der Server („https://wikimedia.org/api/rest_“) hat berichtet: „Cannot get mml. Server problem.“): {\displaystyle +\infty } .
      • Estimation of statistical uncertainty on reconstructed Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} yield (need plot).
        • AH: based on number of entries under gaussian and number of entries in background.
        • SJ: based on covariance matrix of the gaussian+background fit.
      • Resolve this last issue using a brute-force check similar to the systematic uncertainty due to the limited MC statistics in the smearing matrix. Split MC production in large number of samples (as much as possible), fit all Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle pi^0} for every production, and use the standard deviation of these results as a brute-force estimate of the uncertainty. Compare this value to the AH and SJ methods.
  • Find reason for apparent discrepancy between Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} and average charged pion multiplicities. Studies are collected in the Systematics section. The problem in short:
    • In a scenario of zero to modest isospin symmetry violation (ISV), the Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} multiplicities are expected to be near the average charged pion multiplicities.
    • When comparing the 1D unfolded multiplicities a big discrepancy is observed, with differences up to 20 %
    • The old 1999 publication shows a very good agreement between neutral and average charged pions.
    • This is too big to be caused by a realistic ISV.
    • Candidates are:
      • Uncontrollable systematic effect.
      • Problem with the and/or charged pion reconstruction.
      • Problem with unfolding or related MC problem.
      • A severe conceptual mistake in the analysis (example?).
      • other?
    • Systematic studies can be found in the next section.
  • More to come (maybe change analysis method based on results of systematic studies)

Systematic Studies for Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} Multiplicities

Stability and accuracy of acceptance correction

  • Varying the fiducial calorimeter cuts should not cause major variations in the unfolded, acceptance corrected results. Comparing the standard fiducial volume with reduced volumes yield variations in the Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} yield of the order of 10%. Calorimeter volumes:
    • Standard: Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle _x_\text{clus}_ < 125 } and Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 33 < _y_\text{clus}_ < 105}
    • Medium reduced: and Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 39 < _y_\text{clus}_ < 97.5}
    • Highly reduced: Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle _x_\text{clus}_ < 100 } and Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 45 < _y_\text{clus}_ < 90}

eca: I think after our discussion on tuesday we agree it is not completely true that this cuts make no affects.

SJ: Indeed, these studies only prove that the unfolding and acceptance correction do their job properly, (and also that we get a tremendous drop in statistics once we start cutting big chunks out of the calo). The new fiducial checks have the potential to really have a major impact on the multiplicities.

  • Comparison of MC and data Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} acceptance.
    • In-acceptance comparison of Fehler beim Parsen (Konvertierungsfehler. Der Server („https://wikimedia.org/api/rest_“) hat berichtet: „Cannot get mml. Server problem.“): {\displaystyle m^{\pi ^{0}}/m^{\pi ^{+/-}}} between raw data and tracked MC:
      • Small but not-negligible discrepancy at low to medium Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} , substantial discrepancy at high Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} .
      • The ratio of these ratios shows an interesting dependence, especially when compared to the next point.
    • Comparison of Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle m^{\pi^0}/m^{\pi^{+/-}}} between unfolded data and Born 4Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi} MC:
      • Both Born MC and unfolded data have results that are noticeably different from unity. Shape of disagreement in unfolded ratio similar to the in-acceptance ratio of ratios.
  • Integration of the multiplicities versus before and after unfolding compared to the same thing with those versus Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle x_B} before and after unfolding can provide a handle on the stability of the extraction and unfolding.
    • The values before unfolding are almost identically the same. This supports the stability of the invariant mass fit with differing background.
    • The values after unfolding show a small disagreement on the 10% level. This plot also shows the same comparison for average charged pions. The difference between the Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} and Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle x_B} unfolding is noticeably smaller than the 10% for Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} .
  • Comparison of 1D and 2D unfolded does not show any significant disagreement between both methods.
    • For 2D binning: added 4 bins in to the standard 9 Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} -bins.
  • Multiplicities corrected with the standard 1D unfolding are in agreement with the results of an old-style acceptance factor correction.
  • Detection efficiencies for Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} are 10% lower than those in the 99 publication (page 6).
    • This is not unexpected, the RICH adds 0.1 radiation length compared to the old setup.

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} Reconstruction.

  • Use modified energy correction formula that gets Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi^0} peak back at PDG mass.
    • Multiplicity distribution is very asymmetric (falls of quickly) versus .
    • Wrong particle energy/momentum will shift this asymmetric distribution having major impact on the binned multiplicity.
    • The hadron momentum in the MC has to agree very well with that data. If this is not the case, the bin-to-bin migration factors will be completely wrong due to the violently asymmetric Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} -dependence.
    • Eduard's correction formula for photons:
      • Gives better peak width.
      • But even lower peak position compared to PDG mass. (Calibrated in different energy range)
      • Try to add multiplicative factor that fixes the peak position.
    • This reduces the issue due to discrepancy between data and MC Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} -distributions. Should give an improvement on the unfolded results.
  • Effect of changing invariant mass window for the fit. (AH)
  • Dependence of the fit result to the background function (POL3 or Weibull).
    • Comparison between POL3 and Weibull background, before and after unfolding:
      • Changes appear to be minimal when fitting in same range.
    • Check influence of background window on stability of Weibull fit.
  • other?

Thoughts

  • Possible to do a 2D unfolding of the with sufficient MC data?

Presentation

  • Compile a full and consistent set of plots.
  • Use wiki to dynamically write the release report. The working draft can be found here.

Interpretation: SIDIS Factorization and Extraction of LO PDFs

Other Studies

eca: are these plots for the acceptance or for 4pi?

SJ: all these plots are unfolded to 4pi, the idea to do these studies was prompted by an upcoming Hal C paper, where they use their data in an attempt to fit <pT> and <kT>. There appears to be a lot of interest from the theoretical community in our SIDIS 2D average Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle p_{h\perp}} distributions.

Notes

  • [2010-06-29, SJ]: Francesca asked me if it would be possible to provide them with a fully unfolded multiplicity/SIDIS X-section sample in their 4D binning (5D - 1, after fit). They would like to try using this for their different integrations instead of the PYTHIA cross sections they normally use
    • Not sure whether feasible, they'd like a response before September/October, but preferably sooner.
    • Need to change hardcoded Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle (Q^2, x_B, z, p_{h\perp})} to one in
    • Limiting factor: smearing matrix - do we have enough MC statistics to do the 4D unfolding for the multiplicities? Try to answer before August.
    • If not limited by MC: should be feasible. No systematics etc. needed, just the mults/x-secions. If feasible, then no real work required.