Quakonia Physics for ALICE Muon Spectrometer

Efficiency correction framework

Steps...

  1. Have fun with simulation & reconstruction
    1. Simulate and reconstruct minimum bias event @ pp 10 TeV within the acceptance (-2.4 < y < -4.0) of ALICE Muon Spectrometer
    2. Get muon tracks in the ESD events
    3. Draw P, Pt, Eta & di-muon invariant mass distribution //2. Reproduce the J/Psi result of efficiency correction //3. Implement the correction framework of J/Psi case to Upsilon
  2. Simulate Upsilon using parameterization
    1. AliGenParam with AliGenMUONlib
    2. "CDF pp 10" parameterization
  3. Study the geometrical acceptance
    1. Acceptance = Nspec / Ngen (Y & pt)
    2. Draw acceptance plots as functions of pt and y like PPR II (Table 6.79)
    3. Compare it with the PPR plots
ALERT! NOTE: The simulation results in PPR II of Muon Spectrometer are obtained by "fast simulation" with parameterization

Simulation & Reconstruction (local)

  • Local production
    Date Generation description Event count Run type Comment
    12 Aug 2009 Muon cocktail pp (AliGenMUONCocktailpp class)
    Y (rapidity) range: -4, -2.4
    1000 ppMBias 23 muon tracks out of 1000 ESD tracks: needs huge event generation for sufficient statistics
    14 Aug 2009 1000 pptrg1mu 867 muon tracks out of 877 ESD tracks
    20 Aug 2009 5000 Almost muons, but too low di-muon events ~ 1/10
    26 Aug 2009 Muon cocktail pp with AliGenParam & AliGenMUONlib
    MUON acceptance (-4 < y < -2.4)
    1000 ppUpsilon "CDF pp 10" parameterization
    for testing generation time
    27 Aug 2009 Only AliGenParam with AliGenMUONlib
    MUON acceptance (-4 < y < -2.4)

Simulation using parameterization

  • AliGenParam (npart, *library, parameter, tname)
    • npart: maximum number of particles per event
    • library: AliGenMUONlib
    • parameter:
      • kPhi, kOmega, kEta, kJpsiFamily, kPsiP, kJpsi, kJpsiFromB, kUpsilonFamily, kUpsilonP, kUpsilonPP, kUpsilon, kCharm, kBeauty, kPion, kKaon, kChic0, kChic1, kChic2, kChic
    • tname: "Vogt" (= "Vogt PbPb"), "Vogt pp", "CDF scaled", "CDF pp" (14 TeV), "CDF pp 10" (10 TeV), "Flat"
      AliGenParam *param = new AliGenParam(1, new AliGenMUONlib(), AliGenMUONlib::kUpsilon, "CDF pp 10");
      param->SetPhiRange(0., 360.);
      param->SetOrigin(0.,0.,0.);
      param->SetSigma(0.,0.,5.3);
      param->SetVertexSmear(kPerEvent);
      param->SetPtRange(0, 100);
      param->SetYRange(-4, -2.4);   // muon acceptance
      //param->SetYRange(-8, 8);    // full phase space
      param->SetForceDecay(kDiMuon);
      param->SetTrackingFlag(1);
      param->Init();

Using ALICE Production

  • ALICE production (Job details)
    Production Description Event count
    LHC09a4 pp minimum bias, 10 TeV 129,082,803
    LHC09a10 MUON Cocktail pp, di-muon trigger 3,420,400

How to access ALICE production data

  • In the alien shell (Alien link),
    • cd /alice/sim/PDC_XX(year)/LHCXXXXX/XXXXXX(Run number)/
  • It is allowed to copy the files to local machine
    • cp command in alien shell
    • TFile::Cp() class in ROOT
      cp aliensh:[alice] [1] /alice/cern.ch/user/s/sahn/ >cp /alice/sim/PDC_XX/LHCXXXX/XXXXXX/001/AliESDs.root file:.
      TFile::Cp() root [0] TFile::Cp("alien:///alice/sim/PDC_XX/LHCXXXX/XXXXXX/001/AliESDs.root","file:./AliESDs.root")
      ALERT! NOTE: in case of using TFile::Cp() class, you must specify the destination path and the file name
      TIP Copy macro using TFile::Cp() class: gcp.C

Analysis using PROOF

  • LHC09a04 dataset on PROOF: /COMMON/COMMON/LHC09a4_run8101X
    TIP All datasets on PROOF can be displayed by following command in ROOT: gProof->ShowDataSets()


References

Notes

  • Z. Conesa del Valle, L. Aphecetche, C. Finck and G. Martinez-Garcia, Efficiency determination strategy in the ALICE muon spectrometer, ALICE Internal Note 2009-007 (2009) (EDMS link)

Talks

Indico search results

-- SangUnAhn - 2009-08-19

Topic attachments
I Attachment History Action Size Date Who Comment
C source code filec gcp.C r1 manage 2.2 K 2009-08-20 - 05:17 SangUnAhn Grid copy macro using TFile::Cp() class
Edit | Attach | Watch | Print version | History: r7 < r6 < r5 < r4 < r3 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r7 - 2009-09-04 - SangUnAhn
 
    • 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