(r19) D0ToHHPi0Page < LHCb

LHCb Web>D0ToHHPi0Page (revision 19) (raw view)~~EditAttachPDF~~

<noautolink>

---+!! D0 to hhpi0 analysis

---++!! Documentation

---++ !! Talks
   * Mika Vesterinen, pi0 in the HLT, Hlt Ops meeting, 21st June 2013 [[http://indico.cern.ch/conferenceDisplay.py?confId=237936][link]]
   * Mika Vesterinen, Update on hhpi0, Charm WG meeting, 15th May 2013 [[https://indico.cern.ch/conferenceOtherViews.py?view=lhcb_meeting&confId=246868][link]]
   * Mika Vesterinen, Status of D0->hhpi0, Joint Charm/SLB WG meeting, 27th March 2013 [[https://indico.cern.ch/conferenceOtherViews.py?view=lhcb_meeting&confId=203974][link]]

---++ !! Stripping-line development

   * Documented in [[D0ToHHPi0StrippingDevelopmentPage][this]] page.

---++!! Location of ntuples

   * LATEST version: Stripping 20, ~180 pb<sup>-1</sup> 
   =/afs/hep.man.ac.uk/d/lhcb-charm_1/will/data_files/merged2.root=
   =/castor/cern.ch/user/m/mvesteri/HHPi0/Stripping20_180pb_251112/merged.root=
   =/castor/cern.ch/user/m/mvesteri/HHPi0/Stripping20_180pb_251112/options_Mika.py=
   =/castor/cern.ch/user/m/mvesteri/HHPi0/Stripping20_180pb_251112/ganga_submit.py=
   =/castor/cern.ch/user/m/mvesteri/HHPi0/Stripping20_180pb_251112/merge.sh=

   * 2012 Data (stripping 20, ~ 681 pb<sup>-1</sup>): =/castor/cern.ch/user/p/pgarosi/NtupleData2012Strip20=

   * Generator level simulation produced by Frederic: <br />100k events of D0--&gt;pi-pi+pi0 according to the !BaBar isobar Dalitz model: =/castor/cern.ch/user/g/gersabec/charm/pipipi0/Dalitz_100000evt.root= <p> </p>
---++!! Relevant codes
   * DaVinci/ganga scripts live in svn here
   <verbatim> 
SetupProject Urania v1r1 --build-env
getpack Phys/D0ToHHPi0Analysis
</verbatim>
   * Original scripts from Paola, to produce ntuples of data and MC =/afs/cern.ch/user/p/pgarosi/public/ntuple_code=. <br />In each sub-directory options and LFNs files.
---++ Stripping lines
   * D<sup>*+</sup> &rarr; D<sup>0</sup> (&rarr; K<sup>-</sup> &pi;<sup>+</sup> &pi;<sup>0</sup>) &pi;<sup>+</sup> <br /><br /> To call the stripping lines: 
      * ="Phys/DstarD2hhpi0PromptDst2D2RSPi0RLine/Particles"= for ="[ D*(2010)+ -&gt; (^D0 -&gt; (^K*(892)0 -&gt; ^K- ^pi+) (^pi0 -&gt; ^gamma ^gamma)) ^pi+]cc "=
      * ="Phys/DstarD2hhpi0PromptDst2D2WSPi0RLine/Particles"= for ="[ D*(2010)+ -&gt; (^D0 -&gt; (^K*(892)0 -&gt; ^K+ ^pi-) (^pi0 -&gt; ^gamma ^gamma)) ^pi+]cc "=
      * ="Phys/DstarD2hhpi0PromptDst2D2KKPi0RLine/Particles"= for ="[ D*(2010)+ -&gt; (^D0 -&gt; (^K*(892)0 -&gt; ^K- ^K+) (^pi0 -&gt; ^gamma ^gamma)) ^pi+]cc "=
      * ="Phys/DstarD2hhpi0PromptDst2D2PiPiPi0RLine/Particles"= for ="[ D*(2010)+ -&gt; (^D0 -&gt; (^K*(892)0 -&gt; ^pi- ^pi+) (^pi0 -&gt; ^gamma ^gamma)) ^pi+]cc "=
   * B<sup>-</sup> &rarr; D<sup>0</sup> (&rarr; K<sup>-</sup> &pi;<sup>+</sup> &pi;<sup>0</sup>) &mu;<sup>-</sup> <br /><br /> To call the stripping lines: <br /> 
      * ="Phys/b2D0MuXKPiPi0ResolvedCharmFromBSemiLine/Particles"= for ="[ B- -&gt; (^D0 -&gt; ^K- ^pi+ (^pi0 -&gt; ^gamma ^gamma)) ^mu- ]cc "=
      * ="Phys/b2D0MuXKKPi0ResolvedCharmFromBSemiLine/Particles"= for ="[ B- -&gt; (^D0 -&gt; ^K- ^K+ (^pi0 -&gt; ^gamma ^gamma)) ^mu- ]cc "=
      * ="Phys/b2D0MuXPiPiPi0ResolvedCharmFromBSemiLine/Particles"= for ="[ B- -&gt; (^D0 -&gt; ^pi- ^pi+ (^pi0 -&gt; ^gamma ^gamma)) ^mu- ]cc "=
   * B<sup>0</sup> &rarr; &mu;<sup>-</sup> D<sup>*+</sup> &rarr; D<sup>0</sup> (&rarr; K<sup>-</sup> &pi;<sup>+</sup> &pi;<sup>0</sup>) &pi;<sup>+</sup> <br /><br /> To call the stripping lines: <br /> 
      * ="Phys/b2DstarMuXKPiPi0ResolvedCharmFromBSemiLine/Particles"= for ="[ B0 -&gt; (^D*(2010)+ -&gt; (^D0 -&gt; ^K- ^pi+ (^pi0 -&gt; ^gamma ^gamma)) ^pi+) ^mu- ]cc "=
      * ="Phys/b2DstarMuXKKPi0ResolvedCharmFromBSemiLine/Particles"= for ="[ B0 -&gt; (^D*(2010)+ -&gt; (^D0 -&gt; ^K- ^K+ (^pi0 -&gt; ^gamma ^gamma)) ^pi+) ^mu- ]cc "=
      * ="Phys/b2DstarMuXPiPiPi0ResolvedCharmFromBSemiLine/Particles"= for ="[ B0 -&gt; (^D*(2010)+ -&gt; (^D0 -&gt; ^pi- ^pi+ (^pi0 -&gt; ^gamma ^gamma)) ^pi+) ^mu- ]cc "=
---++!! MC samples

<verbatim>
dkfiles/Dst_D0pi,Kpipi0=DecProdCut,PHSP.dec
bk_query=BKQuery(path='/MC/MC11a/Beam3500GeV-2011-MagDown-Nu2-EmNoCuts/Sim05a/Trig0x40760037Flagged/Reco12a/Stripping17NoPrescalingFlagged/27163400/ALLSTREAMS.DST')
INPUTDATA = bk_query.getDataset()
print INPUTDATA
</verbatim>

---++!! Related Talks

---++ !! Theory papers
   * Y. Grossman et al., Testing for new physics in singly Cabibbo suppressed D decays, [[http://arxiv.org/abs/1204.3557v1][http://arxiv.org/abs/1204.3557v1]]
   * J. Brod et al., Size of direct CP violation in singly Cabibbo-suppressed D decays, [[http://arxiv.org/abs/1111.5000][http://arxiv.org/abs/1111.5000]]

---++ !! Experimental papers
   * http://arxiv.org/abs/hep-ex/0608009 (BaBar analysis of pipipi0 and KKpi0 with 230 fb^-1, 2006) <p> </p>
   * http://arxiv.org/abs/hep-ex/0703037 (BaBar analysis of pipipi0 in B-&gt;D0 K, 2007) <p> </p>
   * http://arxiv.org/abs/0807.4544 (BaBar mixing analysis with Kpipi0, 2012) <p> </p>
---++!! Trigger studies

* [[HltEfficiencyTable][Hlt efficiency table]] -- now including D*-&gt;pi (D-&gt;Kpipi0) MC11a

---++!! Stripping 20 signal mass plots 

---++++!! Cuts
<verbatim>
  D from B:
  TString cuts = "D_M > 0";       
  cuts += "&& nCandidate == 0";                                                                                                                      
  cuts += "&& BL0Global_TOS == 1";
  cuts += "&& BHlt1Global_TOS == 1";
  cuts += "&& BHlt2Global_TOS == 1";
  cuts += "&& ((pi0_M > 120 && pi0_M < 150))";
  if(path.Contains("Dstr")){
    cuts += "&& Dstr_M-D_M > 142 && Dstr_M-D_M < 148";
  }
  Prompt:
  TString cuts = "D_M > 0";      
  cuts += "&& nCandidate == 0";                                                                                                                                  
  cuts += "&& TMath::Abs(D_M-1865) < 60";
  cuts += "&& (DstrL0Global_TIS == 1 || DstrL0Global_TOS==1)";
  cuts += "&& DstrHlt1Global_TOS == 1";
  //cuts += "&& DstrHlt2Global_TOS == 1";    //not applied                                                                                                                                 
  cuts += "&& ((pi0_M > 120 && pi0_M < 150))";
</verbatim>
---++++!! D from SL B, single tag
     <img src="%ATTACHURLPATH%/KPiPi0R_SLsingletag.png" alt="KPiPi0R_SLsingletag.png" width="330" height="260" />
     <img src="%ATTACHURLPATH%/PiPiPi0R_SLsingletag.png" alt="PiPiPi0R_SLsingletag.png" width="330" height="260" />
     <img src="%ATTACHURLPATH%/KKPi0R_SLsingletag.png" alt="KKPi0R_SLsingletag.png" width="330" height="260" />

---++++!! D from SL B, double tag

     <img src="%ATTACHURLPATH%/KPiPi0R_SLdoubletag.png" alt="KPiPi0R_SLdoubletag.png" width="330" height="260" />
     <img src="%ATTACHURLPATH%/PiPiPi0R_SLdoubletag.png" alt="PiPiPi0R_SLdoubletag.png" width="330" height="260" />
     <img src="%ATTACHURLPATH%/KKPi0R_SLdoubletag.png" alt="KKPi0R_SLdoubletag.png" width="330" height="260" />
  
---++++!! Prompt with D* tag

     <img src="%ATTACHURLPATH%/KPiPi0R_Prompt.png" alt="KPiPi0R_Prompt.png" width="330" height="260" />
     <img src="%ATTACHURLPATH%/PiPiPi0R_Prompt.png" alt="PiPiPi0R_Prompt.png" width="330" height="260" />
     <img src="%ATTACHURLPATH%/KKPi0R_Prompt.png" alt="KKPi0R_Prompt.png" width="330" height="260" />
-- Main.MikaVesterinen

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
png	KKPi0R_Prompt.png	r2 r1	manage	20.1 K	2012-11-26 - 18:08	MikaVesterinen
png	KKPi0R_SLdoubletag.png	r2 r1	manage	19.8 K	2012-11-26 - 18:08	MikaVesterinen
png	KKPi0R_SLsingletag.png	r2 r1	manage	18.8 K	2012-11-26 - 18:08	MikaVesterinen
png	KPiPi0R_Prompt.png	r2 r1	manage	17.7 K	2012-11-26 - 18:08	MikaVesterinen
png	KPiPi0R_SLdoubletag.png	r2 r1	manage	19.4 K	2012-11-26 - 18:08	MikaVesterinen
png	KPiPi0R_SLsingletag.png	r2 r1	manage	20.0 K	2012-11-26 - 18:08	MikaVesterinen
png	PiPiPi0R_Prompt.png	r2 r1	manage	19.6 K	2012-11-26 - 18:08	MikaVesterinen
png	PiPiPi0R_SLdoubletag.png	r2 r1	manage	20.7 K	2012-11-26 - 18:08	MikaVesterinen
png	PiPiPi0R_SLsingletag.png	r2 r1	manage	17.2 K	2012-11-26 - 18:08	MikaVesterinen
pdf	SalvatoreTesoroReport.pdf	r1	manage	795.2 K	2013-01-26 - 17:34	MikaVesterinen	Report by Salvatore Tesoro

Topic revision: r19 - 2013-06-24 - MikaVesterinen

LHCb

Public webs

Welcome Guest

- Cern Search
- TWiki Search
- Google Search
LHCb All webs

Copyright &© 2008-2021 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