Analysis
- bkgcatfullrunzoomed.gif:
- No worries about cutting off B mesons travelling out of the detector range in wrong direction!?Bplus_DIRA.gif
- dll mu-e why is there a strange peak around 0 when the likelihood for electron and muon is the same?
- 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 |
- Preliminary Background Run information:
- 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 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.
- Mu- momentum for stripped background sample after cut put in place at P>8GeV in my DaVinci code
--
DanielJohnson - 27 Aug 2008