Alessio Sarti Muon ID studies Page

These pages contain the Muon ID studies summary

Muon Page and LHCb Twiki Page

Log Jobs

Old Jobs

Old jobs can be found in this page

DaVinci v26r3p1

Alessio

Data Stripping TP studies

  • Output can be found in /afs/lnf.infn.it/project/lhcb/user/asarti/workspace/asarti/LocalXML
  • Jobs submitted with cards /afs/cern.ch/user/a/asarti/cmtuser/run2010/strip/LHCb*_Collision10_429*Merged_90000000__CALIBRATIONDST.py : 155 is for MagDw (strip 11) and 215 for MagUp (strip 10)

Job ID Status Name N subjobs Application Backend Data Type Comments Data sub evts
153 completed JpsiDetached 59 DaVinci Dirac Data 10 Jpsi Detached Strip 10 [21/11/10] Mag Up evts
154 completed JpsiDetached 38 DaVinci Dirac Data 10 Jpsi Detached Strip 11 [21/11/10] Mag Dw evts
156 completed Jpsi Prompt 38 DaVinci Dirac Data 10 Jpsi Prompt Strip 11 [21/11/10] Mag Dw evts
157 completed Jpsi Prompt 59 DaVinci Dirac Data 10 Jpsi Prompt Strip 11 [21/11/10] Mag Up evts

DaVinci v25r4p1

Alessio

Data Stripping TP studies

  • Output can be found in /afs/lnf.infn.it/project/lhcb/user/asarti/workspace/asarti/LocalXML
  • Jobs from 72 submitted with cards /afs/cern.ch/user/a/asarti/cmtuser/run2010/strip/LHCb*_Collision10_429*_RealData+RecoStripping-03_90000000__MINIBIASDST.py : 155 is for MagDw and 215 for MagUp

Job ID Status Name N subjobs Application Backend Data Type Comments Data sub evts
77 archived StripTest 182 DaVinci Dirac Data 10 Stripping Lines [18/5/10] ~130M Mag Dw evts no output
78 archived StripTest 144 DaVinci Dirac Data 10 Stripping Lines [18/5/10] ~90M Mag Up evts no output
81 archived StripTest 182 DaVinci Dirac Data 10 Stripping Lines [20/5/10] ~130M Mag Dw evts
82 archived StripTest 144 DaVinci Dirac Data 10 Stripping Lines [20/5/10] ~90M Mag Up evts

Data Tupling JP studies (for JP and bkg studies)

  • Output can be found in /afs/lnf.infn.it/project/lhcb/user/asarti/workspace/asarti/LocalXML
  • Jobs submitted with cards /afs/cern.ch/user/a/asarti/cmtuser/run2010/data/LHCb*_Collision10_429*_RealData+RecoStripping-03_90000000__MINIBIASDST.py : 155 is for MagDw and 215 for MagUp

Job ID Status Name N subjobs Application Backend Data Type Comments Data sub evts
108 archived JpsiData2b 149 DaVinci Dirac Data 10 Jpsi ++,+-,-- [24/5/10] ~90M Mag Up evts
109 archived JpsiData2b 180 DaVinci Dirac Data 10 Jpsi ++,+-,-- [24/5/10] ~130M Mag Dw evts

Data Tupling Tag and probe studies (for misid)

  • Output can be found in /afs/lnf.infn.it/project/lhcb/user/asarti/workspace/asarti/LocalXML
  • Jobs submitted with cards /afs/cern.ch/user/a/asarti/cmtuser/run2010/data/LHCb*_Collision10_429*_RealData+RecoStripping-03_90000000__MINIBIASDST.py : 155 is for MagDw and 215 for MagUp

Job ID Status Name N subjobs Application Backend Data Type Comments Data sub evts
117 submitted JpsiData2b 180 DaVinci Dirac Data 10 Jpsi +- [26/5/10] ~130M Mag Dw evts
118 submitted JpsiData2b 149 DaVinci Dirac Data 10 Jpsi +- [26/5/10] ~90M Mag Up evts

MC

  • Output can be found in /afs/lnf.infn.it/project/lhcb/user/asarti/workspace/asarti/LocalXML
  • Jobs submitted with cards /afs/cern.ch/user/a/asarti/cmtuser/run2010/mc/MC_2010_*prod*_2010-Sim03Reco03-withTruth_30000000__DST.py : 428935 is mag down data ~10M minbias and 429376 is mag up data ~10M minbias.

Job ID Status Name N subjobs Application Backend Data Type Comments Data sub evts
111 completed MC 17 DaVinci Dirac MC 10 Jpsi ++,+-,-- [25/5/10] ~10M Mag Dw evts + added tistos, TCK 1810
112 completed MC 18 DaVinci Dirac MC 10 Jpsi ++,+-,-- [25/5/10] ~10M Mag Up evts + added tistos, TCK 1810
113 completed MC 17 DaVinci Dirac MC 10 Jpsi ++,+-,-- [26/5/10] ~10M Mag Dw evts + added tistos, TCK 1810 + added Evtinfo
114 completed MC 18 DaVinci Dirac MC 10 Jpsi ++,+-,-- [26/5/10] ~10M Mag Up evts + added tistos, TCK 1810 + added Evtinfo
120 submitted MC 29 DaVinci Dirac MC 10 Jpsi ++,+-,-- [28/5/10] ~10M Mag Dw MinBias dimuon evts + added tistos, TCK 1810, added Evtinfo
121 submitted MC 30 DaVinci Dirac MC 10 Jpsi ++,+-,-- [28/5/10] ~10M Mag Dw Inclusive Jpsi evts + added tistos, TCK 1810, added Evtinfo

Datasets

We have 4 ntuples: 345122, 346100, 354178 and 352165 MB events (Tot: 1397565). Number is taken from bookeeping (valid +-O(100); Davinci log removed). The remaining 190 ntuples do correspond to 3139480 MB events [number taken from DaVinci log grepping for 'processed' events]. The grand total is 4537045 evts. The events list can be found in /afs/cern.ch/user/a/asarti/muon_group/vol0/mucalib .

MC type Ganga id Dirac id Tot events so far
Bsmumu 1398 582597 - 19 78524
BdKpi 1410, 1412 585992 - 024 and 586136 - 586165 203749
Bdpipi 1408,1409 585901 - 585933 and 585941 - 585970 194292
BsKK 1407 585878 -885 22365
BsKpi 1406 585836 - 868 93948

To each job correspond a stat of 3.5kev /file (processed from DaVinci)

Muon ID starting point

Starting tips:

  • use DaVinci v19r8
  • Analyzed 140 k Lambda B -> Jpsi (L -> ppi) and 200 k Jpsi Incl

Code area

  • $BASE/DaVinciUser/v8r3/src/MuonIDCalib_2body.cpp, .h --> codice per storare le info per MuonID per ogni decadimento a due corpi

Few explanations are needed:

  • Main exec is FillHist
    • Check for a MC particle with f_ismcParticle [+ some mass dist]
    • Build isinregion array
    • Check for MC decay [muon decay in flight] and MC muon [coming from B,D and Jpsi].
    • Check for preselection cuts and acceptance. Preselection acts on Ecal, Hcal, P, Pt, IPS, chi2 track... Acceptance checks inner and outer in M1 and M5 (why both?)
  • For 'preselection passing particles'
    • ic == ic_presel (ic1,2 and 3 are on momentum ranges)
    • MIP defined with Ecan and Hcal cut
    • ic_ismuon having required ismuon or isloosemuon given from f_IsMuon(iP,nhFoi,mom,isM);
  • Efficiencies computed using ic_ismuon (ismuon over preselected)

Option files area

Macro area macro

Inside macro/muonID

  • MuonIDCalib_jpsi.C, .h --> main macro per il calcolo di distanze etc sia per J/Psi che per Lambda (da chiamare MuonIDCalib_2body.C)
  • landau_fit.C : fit della distanza per i mu con una doppia Landau
  • landau_fit2.C: fit della distanza per i non-mu con una "funzione di accettanza"
  • plotta_probmu.C: produce i due plot in alto della mia slide n.11 presentata al meeting del bs->mumu.
  • scan_histos.C : produce i due plot in basso della stessa slide.
  • plotta_eff_stations.C e compara_eff_stations.C producono il plot della funzione di efficienza per stazione di slide 3
  • plotta_efficiency.C e compara_eff.C : producono (piu' o meno) i plots della slide 4

Open Problems

Z -> Dimuon decays inefficiency

When translating from tracks to objects inside the Muon detector there's an acceptance that should be considered:

Inclusive Jpsi mu efficiency

See talk from G. Sabatino

Calibration

Proto Particle Filters

Selection of control samples

  • Developed under cmtuser/calibMu
  • Managed to run on (verified with script in script/DaVinciJobUtil) : 4676017 L0 events (Lambda) and 1298384 (Jpsi)
  • Output files from this run on L0 min bias DC06 are in muonID/MuonCalibrator/Lambda(Jpsi)Monitor.root

Jpsi Calibration sample

  • Defined DaVinciJob.py + JpsiCalib.opts to be used for Jpsi calib.
  • Added TupleTool
  • Removed Muon RAW data : How to rebuild muonID for cross checks? We'll see
  • Took from unbiased Jpsi standard the template for applying cuts

Selection
Use StdNoPIDsMuons (from StandardMuons)
Use 2 no pid long tracks: took std loose muons and removed the require det Muon
track chi2/nDoF < 6 (added, using trackquality filter criterion)
pt(mu)1 or pt(mu)2 > 500MeV/c (already in unbiased Jpsi)
p(mu)1 and p(mu)2 > 3 GeV/c (added, using kinfilter criterion)
Eecal(mu)1 and Eecal(mu)2 < 3000MeV/c (added, using CombinedParticleMaker, requires calo)
Ehcal(mu)1 and Ehcal(mu)2 < 10000MeV/c (added, using CombinedParticleMaker, requires calo)
chi2/nDoF(J/Y) < 20 (modified from 15 to 20)
(M(mu1,mu2) - MJ/Y) < 300MeV/c (modified from 400 to 300)

The ecal and hcal cuts are performed on the Xcal energy collected 'along the track path'. The energy of the Ecal cluster associated to the track is also defined ('CaloChargedEcal') in proto->info() but is not yet implemented in the ProtoParticleCALOFilter.

Detector requirements:

Job n. output Selection
19 muongroup, vol11
JpsiCalib.PhysDesktop.InputLocations  = [ "StdNoPIDsMuons" ]
JpsiCalib.DecayDescriptor = "[J/psi(1S) -> mu+ mu-]cc" ;
JpsiCalib.CombinationCut = "(ADAMASS('J/psi(1S)')<300*MeV)"
JpsiCalib.DaughtersCuts = { "mu+" : "(TRCHI2<6) & (P>3*GeV) & (PT>500*MeV) & (ISLONG)",
             "mu-" : "(TRCHI2<6) & (P>3*GeV) & (PT>500*MeV) & (ISLONG)" }
JpsiCalib.MotherCut = " (MM>2.7*GeV) & (MM<10.1*GeV) & (VFASPF(VCHI2)<20) "
25 muongroup, vol11
JpsiCalib.PhysDesktop.InputLocations  = [ "StdNoPIDsMuons" ]
JpsiCalib.DecayDescriptor = "[J/psi(1S) -> mu+ mu-]cc" ;
JpsiCalib.CombinationCut = "(ADAMASS('J/psi(1S)')<300*MeV)"
JpsiCalib.DaughtersCuts = { "mu+" : "(TRCHI2DOF<6) & (P>3*GeV) & (PT>500*MeV) & (ISLONG)",
             "mu-" : "(TRCHI2DOF<6) & (P>3*GeV) & (PT>500*MeV) & (ISLONG)" }
JpsiCalib.MotherCut = " (MM>2.7*GeV) & (MM<10.1*GeV) & (VFASPF(VCHI2/VDOF)<20) "

Lambda Calibration sample

  • Defined DaVinciJob.py + LambdaCalib.opts to be used for Lambda calib.
  • Added TupleTool
  • Took from unbiased Jpsi standard the template for applying cuts

Selection
Use DD and LL tracks (removed LD and U combinations) Using NoPidParticles as input
LL selection
Lambda : Vtx chi2 less than 3. AND ipsmin greather than 5 : implemented with DaVinciFilter
Lambda : Pointing direction (angle btw Particle vertex and Primary vertex) greater than 0.999999 : implemented with LoKi Filter
Z diff less than 600 : implemented with LoKi Filter
Lambda : IM 50 MeV (macro were 6) : implemented with DaVinciFilter
DD selection
Lambda: Vtx chi2 less than 5. : implemented with DaVinciFilter
Lambda : Pointing direction (angle btw Particle vertex and Primary vertex) greater than 0.999995 : implemented with LoKi Filter
z diff btw 500 and 2200 : implemented with LoKi Filter
Lambda : IM 50 MeV (macro were 6) : implemented with DaVinciFilter

Removed lambda min pT cut, changed cut on significance of flight distance from PV into an ipsmin cut

Job n. output Selection
20 muongroup, vol11
 
LambdaCalibLL.PhysDesktop.InputLocations  = [ "StdNoPIDsPions","StdNoPIDsProtons" ]
LambdaCalibLL.DecayDescriptor = "[Lambda0 -> p+ pi-]cc" ;
LambdaCalibLL.CombinationCut = " (ADAMASS('Lambda0')<50*MeV) "
LambdaCalibLL.OutputLevel = INFO
LambdaCalibLL.DaughtersCuts = { "p+" : " (PT>400*MeV) & (ISLONG) & (MIPCHI2DV(PRIMARY)>3) ",
                 "pi-" : " (PT>100*MeV) & (ISLONG) & (MIPCHI2DV(PRIMARY)>3) " }
LambdaCalibLL.MotherCut = " (VFASPF(VCHI2)<3) & (0<BPVVDZ) & (BPVVDZ<600) & (BPVDIRA>0.999999) & (MIPCHI2DV(PRIMARY)>5) " 
LambdaCalibDD.DaughtersCuts = { "p+" : "(PT>1100*MeV) & (ISDOWN)",
                 "pi-" : "(PT>250*MeV) & (ISDOWN) & (MIPCHI2DV(PRIMARY)>2) " }
LambdaCalibDD.MotherCut = " (VFASPF(VCHI2)<5) & (500<BPVVDZ) & (BPVVDZ<2200) & (BPVDIRA>0.999995) "
Still to be implemented: cut on the Ks mass. And the cut on the pt of the daughter particles wrt the mother ones (Arm-Podo).

To do

Data

Run on MinBias L0 events, select Jpsi and Lambda without ID.

Lambda in DaVinci v19r8 with option file DVMuonIDCalib_lambda.opts. Nell'area ~galanf/cmtuser/DaVinci_v19r8/Phys/DaVinciUser/v8r3/options/cards trovi il file: mb.opts che sono 4,618,488 Minimum Bias event L0 yes su cui passare il codice.

Analysis

  • Fit con landau . fit
  • Distanza ripesata in bin di impulso
  • Identificazione dei parametri delle landau con l'impulso (pluggare nel codice)?
  • Codice da usare per ottimizzare tagli in calib
  • Purezza in funzione di p

-- AlessioSarti - 28 May 2008

Edit | Attach | Watch | Print version | History: r48 < r47 < r46 < r45 < r44 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r48 - 2010-11-22 - AlessioSarti
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback