Table of contents:

Mission

Want to study the E over p bump as a function of correction methods, different datasets and pT.

People

Tomas Aronsson

Plan

The plan is to use the official framework, PWG4 in my case, to do a simple E over p plot. I use this in my B-tag class when I do track matching to find electrons.

Cuts and conditions

  • Use the PWG4 PartCorr framework
  • Run over LHC10c/d/e pass 2 data, 50M events each (in total 150M events, turned out to be ~170M events)
  • In order to assume compatibility, the AliAnaPhoton.cxx was cloned and used with the accompanying AddTaskPi0.C macro to add task with certain parameters. It's all taken from trunk.
  • So far 4 different cases of non-linearity corrections has been used no-correction, MC, GammaGamma and Test beam
  • For each track exactly ONE cluster is matched: the closest one in dphi/deta.
    • Require track to point within roughly the EMCal
    • Require track pT>1.0 GeV
    • Cut at dE/dx = 70 from the central detectors
    • Match CLOSEST cluster, one cluster for each track. If no clusters, then no E/p.
    • Cut on 0.3 GeV cluster energy
    • Cut on NCells>1




The different non-linearity correction is from AliEMCALRecoUtils.cxx:

  • kPi0MC: energy /= (fNonLinearityParams[0]*exp(-fNonLinearityParams[1]/energy))+((fNonLinearityParams[2]/(fNonLinearityParams[3]*2.*TMath::Pi())*exp(-(energy-fNonLinearityParams[4])*(energy-fNonLinearityParams[4])/(2.*fNonLinearityParams[3]*fNonLinearityParams[3]))));
    • //Non-Linearity correction (from MC with function ([0]*exp(-[1]/E))+(([2]/([3]*2.*TMath::Pi())*exp(-(E-[4])^2/(2.*[3]^2)))))
    • //Double_t fNonLinearityParams[0] = 1.001;
    • //Double_t fNonLinearityParams[1] = -0.01264;
    • //Double_t fNonLinearityParams[2] = -0.03632;
    • //Double_t fNonLinearityParams[3] = 0.1798;
    • //Double_t fNonLinearityParams[4] = -0.522;

  • kPi0GammaGamma: energy /= fNonLinearityParams[0]*(1-fNonLinearityParams[1]*exp(-fNonLinearityParams[2]*energy));
    • //Non-Linearity correction (Nicolas from Dimitri Data with function C*[1-a*exp(-b*E)])
    • //fNonLinearityParams[0] = 0.139393/0.1349766;
    • //fNonLinearityParams[1] = 0.0566186;
    • //fNonLinearityParams[2] = 0.982133;

  • kPi0BeamTest: energy /= fNonLinearityParams[0]/(1+fNonLinearityParams[1]*exp(-energy/fNonLinearityParams[2]));
    • //From beam test, Alexei's results, for different ZS thresholds
    • // th=30 MeV; th = 45 MeV; th = 75 MeV
    • //fNonLinearityParams[0] = 0.107; 1.003; 1.002
    • //fNonLinearityParams[1] = 0.894; 0.719; 0.797
    • //fNonLinearityParams[2] = 0.246; 0.334; 0.358

  • kEvi: there is a plan to run over the correction function from Evi: f=new TF1("f","0.983*(([1]*exp(-x/[2]))+(([2]/([3]*2.*TMath::Pi()))*exp(-(x-[4])^2/(2.*[3]^2))))",0.1,30.0); f->SetParameters(1.01, 9.83529e-01, -1.84235e+02, -2.05019e+00, -5.89423e+00); But this is to be done.

Results

Here are all the plots that are sliced according to pT and then fitted with a gaussian. The mean and sigma is then plotted on three different graphs, one for each data set.

Gaussian fit of the electron bump

  • LHC10c pass2: No correction:
    LHC10cpass2_nocorr.gif

  • LHC10c pass2: MC correction:
    LHC10cpass2_MCcorr.gif

  • LHC10c pass2: Test beam correction:
    LHC10cpass2_Testbeamcorr.gif

  • LHC10d pass2: No correction:
    LHC10dpass2_nocorr.gif

  • LHC10d pass2: MC correction:
    LHC10dpass2_MCcorr.gif

  • LHC10d pass2: Test beam correction:
    LHC10dpass2_Testbeamcorr.gif

  • LHC10e pass2: No correction:
    LHC10epass2_nocorr.gif

  • LHC10e pass2: MC correction:
    LHC10epass2_MCcorr.gif

  • LHC10e pass2: Test beam correction:
    LHC10epass2_Testbeamcorr.gif

Gaussian mean/sigma vs pT

  • LHC10c pass2: Gaussian mean vs pT:
    LHC10cpass2_mean.gif

  • LHC10d pass2: Gaussian mean vs pT:
    LHC10dpass2_mean.gif

  • LHC10e pass2: Gaussian mean vs pT:
    LHC10epass2_mean.gif

  • LHC10c pass2: Gaussian sigma vs pT:
    LHC10cpass2_sigma.gif

  • LHC10d pass2: Gaussian sigma vs pT:
    LHC10dpass2_sigma.gif

  • LHC10e pass2: Gaussian sigma vs pT:
    LHC10epass2_sigma.gif
Topic attachments
I Attachment History Action Size Date Who Comment
GIFgif LHC10cpass2_GammaGammacorr.gif r1 manage 21.2 K 2011-02-15 - 05:07 TomasAronsson LHC10c pass2: GammaGamma correction
GIFgif LHC10cpass2_MCcorr.gif r1 manage 21.3 K 2011-02-15 - 05:07 TomasAronsson LHC10c pass2: MC correction
GIFgif LHC10cpass2_Testbeamcorr.gif r1 manage 21.0 K 2011-02-15 - 05:08 TomasAronsson LHC10c pass2: Test beam correction
GIFgif LHC10cpass2_gaus.gif r1 manage 9.1 K 2011-02-15 - 05:13 TomasAronsson LHC10c pass2: Gaussian vs pT
GIFgif LHC10cpass2_mean.gif r1 manage 9.0 K 2011-02-15 - 18:21 TomasAronsson LHC10c pass2: Gaussian mean vs pT
GIFgif LHC10cpass2_nocorr.gif r1 manage 21.1 K 2011-02-15 - 05:07 TomasAronsson LHC10c pass2: No correction
GIFgif LHC10cpass2_sigma.gif r1 manage 11.3 K 2011-02-15 - 18:22 TomasAronsson LHC10c pass2: Gaussian sigma vs pT
GIFgif LHC10dpass2_GammaGammacorr.gif r1 manage 20.8 K 2011-02-15 - 05:09 TomasAronsson LHC10d pass2: GammaGamma correction
GIFgif LHC10dpass2_MCcorr.gif r1 manage 20.9 K 2011-02-15 - 05:09 TomasAronsson LHC10d pass2: MC correction
GIFgif LHC10dpass2_Testbeamcorr.gif r1 manage 21.1 K 2011-02-15 - 05:10 TomasAronsson LHC10d pass2: Test beam correction
GIFgif LHC10dpass2_gaus.gif r1 manage 8.7 K 2011-02-15 - 05:13 TomasAronsson LHC10d pass2: Gaussian vs pT
GIFgif LHC10dpass2_mean.gif r1 manage 8.3 K 2011-02-15 - 18:22 TomasAronsson LHC10d pass2: Gaussian mean vs pT
GIFgif LHC10dpass2_nocorr.gif r1 manage 20.7 K 2011-02-15 - 05:09 TomasAronsson LHC10d pass2: No correction
GIFgif LHC10dpass2_sigma.gif r1 manage 9.5 K 2011-02-15 - 18:23 TomasAronsson LHC10d pass2: Gaussian sigma vs pT
GIFgif LHC10epass2_GammaGammacorr.gif r1 manage 20.8 K 2011-02-15 - 05:11 TomasAronsson LHC10e pass2: GammaGamma correction
GIFgif LHC10epass2_MCcorr.gif r1 manage 20.6 K 2011-02-15 - 05:11 TomasAronsson LHC10e pass2: MC correction
GIFgif LHC10epass2_Testbeamcorr.gif r1 manage 21.2 K 2011-02-15 - 05:12 TomasAronsson LHC10e pass2: Test beam correction
GIFgif LHC10epass2_gaus.gif r1 manage 8.7 K 2011-02-15 - 05:14 TomasAronsson LHC10e pass2: Gaussian vs pT
GIFgif LHC10epass2_mean.gif r1 manage 8.4 K 2011-02-15 - 18:22 TomasAronsson LHC10e pass2: Gaussian mean vs pT
GIFgif LHC10epass2_nocorr.gif r1 manage 20.8 K 2011-02-15 - 05:10 TomasAronsson LHC10e pass2: No correction
GIFgif LHC10epass2_sigma.gif r1 manage 9.9 K 2011-02-15 - 18:24 TomasAronsson LHC10e pass2: Gaussian sigma vs pT
Edit | Attach | Watch | Print version | History: r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r2 - 2011-02-15 - TomasAronsson
 
    • 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