Analysis

  • bkgcatfullrunzoomed.gif:
    bkgcatfullrunzoomed.gif
  • No worries about cutting off B mesons travelling out of the detector range in wrong direction!?Bplus_DIRA.gif
    Bplus_DIRA.gif
  • dll mu-e why is there a strange peak around 0 when the likelihood for electron and muon is the same? muminus_PIDe.gif
  • Counter output of run over Run1 Events 1-30000 with v loose preselection cuts (mass window on B 4GeV and chi2/ndof vertex fit of B < 100 !). Interesting asymmetry in reconstruction of B+ and B-
    "# B+ -> pi+ mu+ mu- " 29192 279088 9.5604 21.417 0.0000 1026.0
    "# B- -> pi- mu- mu+ " 29192 274857 9.4155 20.984 0.0000 917.00

  • To loop over the root file produced by DaVinci use the following commands (should really run on TChain but my tchain only contains the one file) TFile * filesignal = TFile::Open("rfio:///castor/cern.ch/user/j/johndan/DaVinciAnalysis/Bu_pimumu_Lumi2_Run_1.root"); filesignal.cd("DecayTreeBu2PiMuMu") TTree * tree = gDirectory->Get("Bu2PiMuMuTuple") tree.MakeClass("Optimisation")

  • //STRIP REST OF OPTIONS

//STRIP OPTIONS PreselB2DiMuon PreselB2DiMuon // From DVPreselB2DiMuon.opts: //In ->DVPreselB2DiMuon->DoPreselB2DiMuon TopAlg += { "GaudiSequencer/SeqPreselB2DiMuon" } ; SeqPreselB2DiMuon.Members = { "MakeResonances/PreselB2DiMuon", "PrintHeader/PrintPreselB2DiMuon" } ; //In ->DVPreselB2DiMuon->DoPreselB2DiMuon->PreselB2DiMuon InputLocations = { "Phys/StdDC06LooseMuons" } ; DecayDescriptor = "J/psi(1S) -> mu+ mu-" ; DaughterFilter.Selections = { " mu+ : KinFilterCriterion/LKin, PVIPFilterCriterion/LPVIP" } ; MinPt = 800*MeV ; MinMomentum = 8000*MeV ; MinIPsignif = 3.0 ; MotherFilter.Selections = { "J/psi(1S) : VtxFilterCriterion, FlightDistanceFilterCriterion" } ; MaxChi2 = 10 ; MinFSPV = 3.0 ; HistoProduce = false ; DaughterPlots.Variables = { "P", "Pt", "IPs" } ; MotherPlots.Variables = { "M", "Chi2", "Vz", "P", "Pt", "IPs", "FS" } ; MotherPlots.Minima = { 3000.0, 0.0, -150., 0.0, 0.0, 0.0, 0.0 } ; MotherPlots.Maxima = { 3000.0, 50.0, 150.0, 20000.0, 10000.0, 10.0, 10.0 } ; //In ->DVPreselB2DiMuon->DoPreselB2DiMuon OutputLevel = 3 ; OutputLevel = 4 ; //Added: MeasureTime = true ; OutputLevel = 3 ; Preselection criteria which must be matched in the signal and background samples to compare them fairly in order to optimise cuts

  • Mu- momentum for the stripped background sample - evidence of stripping cut at 8GeV but much false reconstruction so I'll apply a 8GeV cut in my DaVinci jobs too.
    hmuminus_P.gif

  • Mu- momentum for stripped background sample after cut put in place at P>8GeV in my DaVinci code hmuminus_PWithCut8GeV.gif
-- DanielJohnson - 27 Aug 2008


This topic: Main > TWikiUsers > DanielJohnson > BPiMuMuAnalysis
Topic revision: r9 - 2008-08-27 - DanielJohnson
 
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