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
      

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")

    NOTE: in case of using TFile::Cp() class, you must specify the destination path and the file name
    Copy macro using TFile::Cp() class: gcp.C

Analysis using PROOF

• LHC09a04 dataset on PROOF: /COMMON/COMMON/LHC09a4_run8101X
  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

• Philippe's PWG3-MUON-EVO talk: CORRFW: Quarkonia

Indico search results

• keyword: CORRFW: what is Correction Framework?
• keyword: Xavier: X. Lopez's talks about efficiency correction framework on PWG3-MUON-EVO

-- SangUnAhn - 2009-08-19