Newer
Older
Presentations / tau23mu_run2 / VRD_slides_04_10_2017 / mchrzasz.tex
@Marcin CHrzaszcz Marcin CHrzaszcz on 3 Oct 2017 18 KB tau23mu vrd
\documentclass[11 pt,xcolor={dvipsnames,svgnames,x11names,table}]{beamer}

\usepackage[english]{babel}
\usepackage{polski}


\usetheme[
	bullet=circle,		% Other option: square
	bigpagenumber,		% circled page number on lower right
	topline=true,			% colored bar at the top of the frame
	shadow=false,			% Shading for beamer blocks
	watermark=BG_lower,	% png file for the watermark
	]{Flip}

%\logo{\kern+1.em\includegraphics[height=1cm]{SHiP-3_LightCharcoal}}


\usepackage[lf]{berenis}
\usepackage[LY1]{fontenc}
\usepackage[utf8]{inputenc}

\usefonttheme{professionalfonts}
\usepackage[no-math]{fontspec}
\defaultfontfeatures{Mapping=tex-text}	% This seems to be important for mapping glyphs properly

\setmainfont{Times}			% Beamer ignores "main font" in favor of sans font
\setsansfont{Times}			% This is the font that beamer will use by default
% \setmainfont{Gill Sans Light}		% Prettier, but harder to read

\setbeamerfont{title}{family=\fontspec{Gillius ADF}}

\input t1augie.fd

%\newcommand{\handwriting}{\fontspec{augie}} % From Emerald City, free font
%\newcommand{\handwriting}{\usefont{T1}{fau}{m}{n}} % From Emerald City, free font
% \newcommand{\handwriting}{}	% If you prefer no special handwriting font or don't have augie

%% Gill Sans doesn't look very nice when boldfaced
%% This is a hack to use Helvetica instead
%% Usage: \textbf{\forbold some stuff}
%\newcommand{\forbold}{\fontspec{Arial}}

\usepackage{graphicx}
\usepackage[export]{adjustbox}

\usepackage{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{bm}
\usepackage{colortbl}
\usepackage{mathrsfs} 			% For Weinberg-esque letters
\usepackage{cancel}				% For "SUSY-breaking" symbol
\usepackage{slashed}            % for slashed characters in math mode
\usepackage{bbm}                % for \mathbbm{1} (unit matrix)
\usepackage{amsthm}				% For theorem environment
\usepackage{multirow}			% For multi row cells in table
\usepackage{arydshln} 			% For dashed lines in arrays and tables
\usepackage{siunitx}
\usepackage{xhfill}
\usepackage{grffile}
\usepackage{textpos}
\usepackage{subfigure}
\usepackage{tikz}

%\usepackage{hepparticles}
\usepackage[italic]{hepparticles}

\usepackage{hepnicenames}

% Drawing a line
\tikzstyle{lw} = [line width=20pt]
\newcommand{\topline}{%
  \tikz[remember picture,overlay] {%
    \draw[crimsonred] ([yshift=-23.5pt]current page.north west)
             -- ([yshift=-23.5pt,xshift=\paperwidth]current page.north west);}}



% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
\usepackage{tikzfeynman}		% For Feynman diagrams
\usetikzlibrary{arrows,shapes}
\usetikzlibrary{trees}
\usetikzlibrary{matrix,arrows} 				% For commutative diagram
% http://www.felixl.de/commu.pdf
\usetikzlibrary{positioning}				% For "above of=" commands
\usetikzlibrary{calc,through}				% For coordinates
\usetikzlibrary{decorations.pathreplacing}  % For curly braces
% http://www.math.ucla.edu/~getreuer/tikz.html
\usepackage{pgffor}							% For repeating patterns

\usetikzlibrary{decorations.pathmorphing}	% For Feynman Diagrams
\usetikzlibrary{decorations.markings}
\tikzset{
	% >=stealth', %%  Uncomment for more conventional arrows
	vector/.style={decorate, decoration={snake}, draw},
	provector/.style={decorate, decoration={snake,amplitude=2.5pt}, draw},
	antivector/.style={decorate, decoration={snake,amplitude=-2.5pt}, draw},
	fermion/.style={draw=gray, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=gray]{>}}}},
	fermionbar/.style={draw=gray, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=gray]{<}}}},
	fermionnoarrow/.style={draw=gray},
	gluon/.style={decorate, draw=black,
		decoration={coil,amplitude=4pt, segment length=5pt}},
	scalar/.style={dashed,draw=black, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=black]{>}}}},
	scalarbar/.style={dashed,draw=black, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=black]{<}}}},
	scalarnoarrow/.style={dashed,draw=black},
	electron/.style={draw=black, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=black]{>}}}},
	bigvector/.style={decorate, decoration={snake,amplitude=4pt}, draw},
}

% TIKZ - for block diagrams,
% from http://www.texample.net/tikz/examples/control-system-principles/
% \usetikzlibrary{shapes,arrows}
\tikzstyle{block} = [draw, rectangle,
minimum height=3em, minimum width=6em]




\usetikzlibrary{backgrounds}
\usetikzlibrary{mindmap,trees}	% For mind map
\newcommand{\degree}{\ensuremath{^\circ}}
\newcommand{\E}{\mathrm{E}}
\newcommand{\Var}{\mathrm{Var}}
\newcommand{\Cov}{\mathrm{Cov}}
\newcommand\Ts{\rule{0pt}{2.6ex}}       % Top strut
\newcommand\Bs{\rule[-1.2ex]{0pt}{0pt}} % Bottom strut

\graphicspath{{images/}}	% Put all images in this directory. Avoids clutter.

% SOME COMMANDS THAT I FIND HANDY
% \renewcommand{\tilde}{\widetilde} % dinky tildes look silly, dosn't work with fontspec
\newcommand{\comment}[1]{\textcolor{comment}{\footnotesize{#1}\normalsize}} % comment mild
\newcommand{\Comment}[1]{\textcolor{Comment}{\footnotesize{#1}\normalsize}} % comment bold
\newcommand{\COMMENT}[1]{\textcolor{COMMENT}{\footnotesize{#1}\normalsize}} % comment crazy bold
\newcommand{\Alert}[1]{\textcolor{Alert}{#1}} % louder alert
\newcommand{\ALERT}[1]{\textcolor{ALERT}{#1}} % loudest alert
%% "\alert" is already a beamer pre-defined
\newcommand*{\Scale}[2][4]{\scalebox{#1}{$#2$}}%

\def\Put(#1,#2)#3{\leavevmode\makebox(0,0){\put(#1,#2){#3}}}

\usepackage{gmp}
\usepackage[final]{feynmp-auto}

\usepackage[backend=bibtex,style=numeric-comp,firstinits=true]{biblatex}
\bibliography{bib}
\setbeamertemplate{bibliography item}[text]

\makeatletter\let\frametextheight\beamer@frametextheight\makeatother

% suppress frame numbering for backup slides
% you always need the appendix for this!
\newcommand{\backupbegin}{
	\newcounter{framenumberappendix}
	\setcounter{framenumberappendix}{\value{framenumber}}
}
\newcommand{\backupend}{
	\addtocounter{framenumberappendix}{-\value{framenumber}}
	\addtocounter{framenumber}{\value{framenumberappendix}}
}


\definecolor{links}{HTML}{2A1B81}
%\hypersetup{colorlinks,linkcolor=,urlcolor=links}

% For shapo's formulas:
\def\lsi{\raise0.3ex\hbox{$<$\kern-0.75em\raise-1.1ex\hbox{$\sim$}}}
\def\gsi{\raise0.3ex\hbox{$>$\kern-0.75em\raise-1.1ex\hbox{$\sim$}}}
\newcommand{\lsim}{\mathop{\lsi}}
\newcommand{\gsim}{\mathop{\gsi}}
\newcommand{\wt}{\widetilde}
%\newcommand{\ol}{\overline}
\newcommand{\Tr}{\rm{Tr}}
\newcommand{\tr}{\rm{tr}}
\newcommand{\eqn}[1]{&\hspace{-0.7em}#1\hspace{-0.7em}&}
\newcommand{\vev}[1]{\rm{$\langle #1 \rangle$}}
\newcommand{\abs}[1]{\rm{$\left| #1 \right|$}}
\newcommand{\eV}{\rm{eV}}
\newcommand{\keV}{\rm{keV}}
\newcommand{\GeV}{\rm{GeV}}
\newcommand{\im}{\rm{Im}}
\newcommand{\disp}{\displaystyle}
\def\be{\begin{equation}}
\def\ee{\end{equation}}
\def\ba{\begin{eqnarray}}
\def\ea{\end{eqnarray}}
\def\d{\partial}
\def\l{\left(}
\def\r{\right)}
\def\la{\langle}
\def\ra{\rangle}
\def\e{{\rm e}}
\def\Br{{\rm Br}}

\def\ARROW{{\color{JungleGreen}{$\Rrightarrow$}}\xspace}
\def\ARROWR{{\color{WildStrawberry}{$\Rrightarrow$}}\xspace}

\author{ {Marcin Chrzaszcz} (CERN)}
\institute{UZH}
\title[$\tau \to \mu \mu\mu$ Run 2 plans]{$\tau \to \mu \mu\mu$ Run 2 plans}
\date{4 October 2017}


\begin{document}
\tikzstyle{every picture}+=[remember picture]

{
\setbeamertemplate{sidebar right}{\llap{\includegraphics[width=\paperwidth,height=\paperheight]{bubble2}}}
\begin{frame}[c]%{\phantom{title page}}
\begin{center}
\begin{center}
	\begin{columns}
		\begin{column}{0.75\textwidth}
	\flushright \bfseries \Large {$\tau \to \mu \mu\mu$ Run 2 plans}
		\end{column}
                \begin{column}{0.02\textwidth}
                  {~}
                  \end{column}
                \begin{column}{0.23\textwidth}
                 % \hspace*{-1.cm}
                  \vspace*{-3mm}
                  \includegraphics[width=0.6\textwidth]{lhcb-logo}
                  \end{column}

	\end{columns}
\end{center}
	\quad
	\vspace{3em}
\begin{columns}
\begin{column}{0.44\textwidth}
\flushright \vspace{-2.8em} {  \fontspec{Zapfino} Marcin Chrzaszcz\\\vspace{-0.1em}\small \href{mailto:mchrzasz@cern.ch}{mchrzasz@cern.ch}}

\end{column}
\begin{column}{0.53\textwidth}
\hspace{2.6cm}
\includegraphics[height=1.6cm]{cern}
\end{column}
\end{columns}


\vspace{1em}
	\footnotesize{\large With M. Palutan (INFN Frascati), F. Teubert (CERN),
	
	 R. Vazquez Gomez (CERN)}
\vspace{0.5em}

	\textcolor{normal text.fg!50!Comment}{VRD meeting, CERN\\October 4, 2017}
\end{center}
\end{frame}
}

\begin{frame}\frametitle{Lepton Flavour/Number Violation}
\begin{small}
 Lepton Flavour Violation(LFV):
\end{small}


\begin{footnotesize}

After $\Pmuon$ was discovered it was logical to think of it as an excited $\Pelectron$.
\begin{columns}
\column{3in}
\begin{itemize}
\item Expected: $B(\mu\to\Pe\gamma) \approx  10^{-4}$
\item Unless another $\Pnu$, in intermediate vector boson loop cancels. 
\end{itemize}
\column{2in}
{~}\includegraphics[width=0.98\textwidth]{rabi.png}

\end{columns}
\begin{columns}
\column{1in}
{~}
\column{3in}
\begin{block}{I.I.Rabi:}
"Who ordered that?"
\end{block}
\column{2in}
{~}{~}{~}\includegraphics[scale=0.08]{II_Rabi.jpg}

\end{columns}


\begin{itemize}
\item Up to this day charged LFV is being searched for in various decay modes.
\item LFV was already found in neutrino sector.
\end{itemize}
\end{footnotesize}


\begin{footnotesize}

\begin{columns}
\column{3.5in}
\begin{small}
 Lepton Number Violation (LNV)
\end{small}

\begin{itemize}
\item Even with LFV, lepton number is a conserved quantity. 
\item Many new thesis predict it violation(Majorana neutrinos)
\item Searched in so called Neutrinoless double $\beta$ decays.
\end{itemize}

\column{1.5in}
\includegraphics[width=0.8\textwidth]{Double_beta_decay_feynman.png}

\end{columns}

\end{footnotesize}
%Double_beta_decay_feynman.png

\end{frame}


 \begin{frame}
        \frametitle{"The Rule of Three"}
        
        \ARROW Since 4 years thee experiments are chasing the $\tau \to \mu \mu \mu$ decay:
        
 \begin{columns}
   %     \column{2.5in}
  
   \begin{column}{2.1in}
            \begin{alertblock}{ $\Ptau \to \Pmu \Pmu \Pmu$ limits ($ \color{white} 90\,\%$ CL)}

              \begin{description}
                \item[BaBar(FC)] $3.3\times 10^{-8}$
                \item[Belle(FC)] $2.1\times 10^{-8}$
                \item[LHCb(CLs)] $4.6\times 10^{-8}~(3 {\rm fb}^{-1})$ 
                \item[HFAG(CLs)] $1.2 \times 10^{-8}$
              \end{description}
            \end{alertblock}   
\end{column}
   \begin{column}{2.5in}    
          \includegraphics[width=1\textwidth]{zom.png}\\
   

\end{column}

 \end{columns}
 {~}\\
    \begin{itemize}
      	\item Thanks to 3 experiments we have a world limit: $\mathcal{B}(\Ptau \to \Pmu \Pmu \Pmu)< 1.2 \times 10^{-8}$ at 90\% CL.
      \item Would be nice to put LHCb ahead of the game!
   	\item We are getting close to B-factories with the Run2 data but still would be nice to have some boost!

         \end{itemize}    
 

      \end{frame}     


\begin{frame}[c]{Extra Gain ? }
	\begin{minipage}{\textwidth}

\ARROW The stripping that we are using takes 3 \texttt{StdLooseMuons} and construct a $\tau$ candidate. \\
\ARROW Now the stripping is loose in other variables.

\begin{center}
\begin{tabular}{| c || c | c |}
\hline
 &  \multicolumn{2}{c|}{Gain wrt. 3 isM $[\%]$ }\\
Category &  2 isM & 2isM + 3AccM \\ \hline
$ B \to \tau$  & $44.53$ & $12.67$  \\        
$ B \to D_s \to \tau$  & $44.37$ & $12.44$ \\
$ D_s \to \tau $    & $43.78$ & $12.34$ \\
$ B \to D \to \tau $    & $46.31$ & $13.97$\\
$ D \to \tau$         & $43.70$ & $11.87$ \\ \hline




\end{tabular}
\end{center}
Legend:
\begin{itemize}
\item isM -IsMuon.
\item AccM - Muon in Acceptance.
\end{itemize}	
	

\ARROW There is extra gain $\mathcal{O}(40~\%)$ that could be included!\\
\ARROW We have already modified the Stripping line and will be run in the S28 restripping.


	\end{minipage}\end{frame}


\begin{frame}[c]{Can we trigger?}


\begin{center}
\begin{Huge}
Initial trigger study
\end{Huge}
\end{center}
\end{frame}


\begin{frame}[c]{IsMuon == 3  WITH TOS BUG FIXED!!}
	\begin{minipage}{\textwidth}
\begin{center}
\begin{tabular}{|c ||  c | c|}
\hline
Name & Eff & Gain wrt. to previous \\ \hline
Hlt2TriMuonTau23Mu &  0.76343 & 0.763439 \\
Hlt2SingleMuonLowPT & 0.764768 & 0.001328\\
Hlt2CharmHadDpToPimPipPipTurbo & 0.76515 & 0.00038\\
Hlt2RareCharmD02KPi & 0.76528 & 0.00012\\
Hlt2SingleMuon & 0.765369 & 0.00008\\
Hlt2LowMultDiMuon & 0.765454 & 0.00006\\
Hlt2XcMuXForTauB2XcMu & 0.765454 & 0.0\\ \hline
\end{tabular}

\end{center}
\begin{alertblock}{How to read table}
Example: Hlt2CharmHadDpToPimPipPipTurbo, 0.76444, 0.00071\\
{\color{green}Eff. of }Hlt2CharmHadDpToPimPipPipTurbo ||  Hlt2SingleMuonLowPT || Hlt2TriMuonTau23Mu is $0.76444$.\\
{\color{green}Gain of adding } Hlt2CharmHadDpToPimPipPipTurbo to previous  Hlt2SingleMuonLowPT || Hlt2TriMuonTau23Mu is $0.00071$.
\end{alertblock}



	\end{minipage}
\end{frame}















\begin{frame}[c]{IsMuon == 2 WITH TOS BUG FIXED!!}
	\begin{minipage}{\textwidth}
\begin{center}
\begin{tabular}{|c ||  c | c|}
\hline
Name & Eff & Gain wrt. to previous \\ \hline
Hlt2CharmHadDpToPimPipPipTurbo & $0.22792$ & $0.22792$ \\
Hlt2SingleMuon & $0.28015$ &  $0.052231$ \\
Hlt2SingleMuonLowPT & $0.28975$ & $0.00960$ \\
Hlt2RareCharmD02KPi & $0.293869$ & $0.004115$ \\
Hlt2CharmHadDspToKpPimPipTurbo& $0.297034$ & $0.003165$\\
Hlt2CharmHadInclSigc2PiLc2HHX  & $0.29946$ & $0.00242$ \\
Hlt2RareCharmD02MuMu & $0.30115$ & $0.00168$\\
Hlt2LowMultDiMuon & $0.302732$ & $0.001582$\\
Hlt2TopoMu3Body & $0.303999$ & $0.001266$ \\
Hlt2CharmHadDspToKmPipPipTurbo & $0.30452$ & $0.00052$ \\ \hline

\hline
\end{tabular}

\end{center}




	\end{minipage}
\end{frame}


\begin{frame}[c]{Trigger conclusion}
	\begin{minipage}{\textwidth}

\begin{itemize}
\item For the 3 \texttt{ISMuon} case we are safe.
\item The 2  \texttt{ISMuon} :
\begin{itemize}
\item Problem will be fixed in 2017 data after bug fix of the Hlt2DiMuonDetached line.
\item For 2015, 2016 we will define our strategy later.
\end{itemize}
\end{itemize}




	\end{minipage}
\end{frame}




\begin{frame}[c]{MC filtering}
	\begin{minipage}{\textwidth}
\ARROW Producing MC becomes more and more a problem.\\
\ARROW Need a smart way to reduce the data set at the production level:
\begin{alertblock}{MC Filtering:}
Select only the events that are needed for the analysis \ARROW Saving disk space.\\
Caveat: once lost events cannot be recovered.
\end{alertblock}
\begin{exampleblock}{What to filter on?}
\ARROWR Stripping?\\
\ARROWR Trigger?\\ 
\ARROWR Reconstruction!
\end{exampleblock}
\ARROW Filtering on MCTRUTHMACHED is the safest thing to do!


	\end{minipage}\end{frame}



\begin{frame}[c]{How does this work?}
	\begin{minipage}{\textwidth}
\ARROW You check that the Particle lists (ex. \texttt{Phys/StdAllNoPIDsMuons/Particles}) contain all particles from the decay (ex. $\tau \to 3 \mu$).\\
\ARROW Very robust! Sensitive only to bugs in reconstruction and TruthMaching.\\
\ARROW Already tested this in $\PB \to \PKstar \mu \mu$ with success.\\
\ARROW Now implemented for $\tau \to \mu\mu\mu$ analysis.


	\end{minipage}\end{frame}


\begin{frame}[c]{Filter efficiency}
	\begin{minipage}{\textwidth}
\begin{center}
\begin{tabular}{|c|c|}
\hline
Channel & Fraction of save events\\ \hline
$\tau \to 3 \mu$ & $0.553$ \\
$D_s \to \phi \pi$ & $0.481$ \\\hline

\end{tabular}
\end{center}

\ARROW FILTER IS IMMUNE TO ANY \texttt{isMuon} REQUIREMENT.\\
\ARROW We save around $50~\%$ of disk space!




	\end{minipage}\end{frame}



\begin{frame}[c]{Cross check}
	\begin{minipage}{\textwidth}

\ARROW As a cross-check is performed:
\begin{itemize}
\item Strip the filtered dst.
\item Strip the non filtered dst.
\end{itemize}
\ARROW After truthmaching they both should be identical.\\
\ARROW And they are:\\
\includegraphics[angle=-90,width=0.49\textwidth]{Ds.pdf}
\includegraphics[angle=-90,width=0.49\textwidth]{tau23mu.pdf}


(Yes there are two IDENTICAL histograms on each plot :) )
	\end{minipage}\end{frame}













\begin{frame}[c]{Conclusions}
	\begin{minipage}{\textwidth}

\ARROW There is a potential $40~\%$ gain of the signal!\\
\ARROW Stripping implemented.\\
\ARROW Filter for MC production developed and tested!\\
\ARROW We gain $50~\%$ on disk space!\\
\ARROW Trigger strategy on the way.\\
\ARROW We would like to ask for 2015,2016,2017 Filtered signal MC.

	\end{minipage}
\end{frame}



\backupbegin




\begin{frame}[c]{IsMuon == 3  WITH TOS BUG!!}
	\begin{minipage}{\textwidth}
\begin{center}
\begin{tabular}{|c ||  c | c|}
\hline
Name & Eff & Gain wrt. to previous \\ \hline
Hlt2TriMuonTau23MuDecision & 0.76343 & 0.76343 \\
Hlt2SingleMuonLowPTDecision& 0.76372 & 0.00136\\
Hlt2CharmHadDpToPimPipPipTurbo & 0.76444 & 0.00071\\
Hlt2SingleMuonDecision & 0.76459 & 0.00015\\
Hlt2RareCharmD02KPiDecision &0.76474 & 0.00015\\
Hlt2LowMultDiMuon & 0.76482 & 0.00008\\
Hlt2XcMuXForTauB2XcMuDecision & 0.76482 & 0.0\\
%Hlt2XcMuXForTauB2XcFakeMuDecision & 0.76482 & 0.0\\
%Hlt2TriMuonDetachedDecision & 0.76482 & 0.0\\
%Hlt2TransparentDecision & 0.76482 & 0.0\\
\hline
\end{tabular}

\end{center}
\begin{alertblock}{How to read table}
Example: Hlt2CharmHadDpToPimPipPipTurbo, 0.76444, 0.00071\\
{\color{green}Eff. of }Hlt2CharmHadDpToPimPipPipTurbo ||  Hlt2SingleMuonLowPT || Hlt2TriMuonTau23Mu is $0.76444$.\\
{\color{green}Gain of adding } Hlt2CharmHadDpToPimPipPipTurbo to previous  Hlt2SingleMuonLowPT || Hlt2TriMuonTau23Mu is $0.00071$.
\end{alertblock}



	\end{minipage}
\end{frame}


\begin{frame}[c]{IsMuon == 2 WITH TOS BUG!!}
	\begin{minipage}{\textwidth}
\begin{center}
\begin{tabular}{|c ||  c | c|}
\hline
Name & Eff & Gain wrt. to previous \\ \hline
Hlt2CharmHadDpToPimPipPip & $0.18769$ &  $0.18769$\\
Hlt2SingleMuonD & $0.25864$ & $0.07095$ \\
Hlt2SingleMuonLowPT & $0.27076$ & $0.01212$\\
Hlt2RareCharmD02KPi & $0.27431$ & $0.003543$\\
Hlt2CharmHadDspToKpPimPip & $0.27738$ & $0.00307$\\
Hlt2CharmHadInclSigc2Pi & $0.279347$ & $0.00195$ \\
Hlt2RareCharmD02MuMu & $0.28121$ & $0.00186$ \\
Hlt2LowMultDiMuon\_PS & $0.28279$ & $0.00158$ \\
Hlt2TopoMu3Body & $0.28400$ & $0.001212$\\
Hlt2PIDKs2PiPiLLTurboCalib & $0.284568$ & $0.00055$\\ 
\hline
\end{tabular}

\end{center}




	\end{minipage}
\end{frame}







\backupend

\end{document}