B->ppbartaunu analysis

People working on this:

Mark Smith, Matthew Tilley, Mitesh Patel, Ryan Newcombe

-- MarkElliotSmith - 2017-12-04

Presentations:

Repository:

https://gitlab.cern.ch/lhcb-slb/B2pplnu

Stripping:

https://gitlab.cern.ch/lhcb/Stripping/blob/master/Phys/StrippingArchive/python/StrippingArchive/Stripping28/StrippingSL/StrippingB2PPbarMuForTauMu.py

Stripping Info:

For the restripping campaign the two fake proton lines (opposite and same-sign) and the TOPO line have been restripped. The changes are:

  • Tighten the track $\chi^{2} / DOF$ cut from 6 to 3 on the protons. This cut should be applied offline anyway.
  • Tighten the $p\overline{p}$ vtx $\chi^{2}$ cut from 10 to 8.
  • Change prescale on the fake proton lines from 0.15 to 0.06 to keep within the rate.
  • Fix to include both sign of fake proton (take wrt to the muon) - previously we had only taken one which missed a large chunk.
  • Include 3 and 4 body TOPO HLT2 in the topo line.
Ensure that the vtx $\chi^{2}$ and track $\chi^{2}$ are consistently applied to all the samples!

Trigger

https://gitlab.cern.ch/lhcb/Hlt/tree/2017-reference/Hlt/Hlt2Lines/python/Hlt2Lines/XcMuXForTau

Data Location:

  • 2018 - J/psi K - /LHCb/Collision18/Beam6500GeV-VeloClosed-Mag[Down,Up]/Real Data/Reco18/Stripping34/90000000/LEPTONIC.MDST

DEC files of interest

The GEN level cuts on the background modes should be:

# tightCut.Cuts = {
#'[p+]cc' : " in_range( 0.010 , GTHETA , 0.400 )& ( GPT > 750 * MeV ) & (GP > 14600 * MeV)" ,
#'[mu-]cc' : " in_range( 0.010 , GTHETA , 0.400 ) & (GP > 2900 * MeV)"
# }

Mode No File Comments Outstanding Issues Generated
$ \tau $ signal 12513060   No GEN cuts EvtGen model  
$\mu$ normalisation 12513050   No GEN cuts EvtGen model  
$\Delta$ cocktail 12813411 Bu_Delpbarmunu,pX=TightCut2.dec GEN cuts EvtGen model  
$N^{*}$ cocktail 12813401 Bu_pNstmunu,pX=TightCut2.dec GEN cuts EvtGen model  
$B^{0}\to p\bar{p}DX$ 11774004 - needs committing Bd_D0ppbarX,Xmunu=TightCut.dec GEN cuts    
$B^{+}\to p\bar{p}DX$ 12773002 - needs committing Bu_D+ppbarX,Xmunu=TightCut.dec GEN cuts    
$p\Lambda_{c}^{*}\mu\nu$ 12877401 - needs committing Bu_Lcpipipbarmunu,pX=TightCut2.dec GEN cuts $p_{T}$ cut on muon  
$p\Lambda_{c}\mu\nu$, $\Lambda_{c}\to pX$ 12875412 Bu_Lcpbarmunu,pX=TightCut2.dec GEN cuts    
$p\Lambda_{c}\mu\nu$, $\Lambda_{c}\to pK\pi$ 12875010 Bu_Lcpbarmunu,pKpi=TightCut.dec only has acceptance $\theta$ cuts    
$p\Lambda_{c}\mu\nu$, $\Lambda_{c}\to pK^{0}\pi^{0}$ ----        
$ppD\mu\nu$ 12572000   GEN cuts need fixing plus more D modes $p_{T}$ cut on muon  

HistFactory extras

There is a nice set of slides introducing HistFactory here: https://indico.cern.ch/event/509900/contributions/2031935/

That link also includes a patch for HistFactory to implement Barlow-Beeston in more than 1 dimension. There are some more patches and speed gains to be found at https://twiki.cern.ch/twiki/bin/viewauth/LHCbPhysics/HistFactoryInfo

I have only tested the BB fix (HistFactoryModelUtils.cxx) - I couldn't get the other speed patches to work.

Compiling ROOT

Compiling ROOT can be a bit of a pain. The compilers on slc6 are ancient so you need to point CMake in the direction of an LCG installation to use instead. Also make sure you enable all the options.

  • Extract the ROOT source code to some directory (i.e. ~/root_source)
  • In the source directory you need to alter the CMakeLists.txt to use a reasonably modern version of gcc and g++ by adding the following:
    • SET(CMAKE_C_COMPILER /cvmfs/lhcb.cern.ch/lib/lcg/releases/gcc/6.2.0binutils/x86_64-slc6/bin/gcc)
      SET(CMAKE_CXX_COMPILER /cvmfs/lhcb.cern.ch/lib/lcg/releases/gcc/6.2.0binutils/x86_64-slc6/bin/g++)
  • You need to alter your PATH and LD_LIBRARY_PATH to pick up the correct libraries etc (note this is for C shell but I am sure you can figure it out for bash):
    • setenv PATH /cvmfs/lhcb.cern.ch/lib/lcg/releases/gcc/6.2.0binutils/x86_64-slc6/bin\:$PATH
      setenv LD_LIBRARY_PATH $LD_LIBRARY_PATH\:/cvmfs/lhcb.cern.ch/lib/lcg/releases/gcc/6.2.0/x86_64-slc6/lib64/\:/cvmfs/lhcb.cern.ch/lib/lcg/releases/gcc/6.2.0/x86_64-slc6/lib/\:/cvmfs/lhcb.cern.ch/lib/lcg/releases/gcc/6.2.0/x86_64-slc6/lib64/\:/cvmfs/lhcb.cern.ch/lib/lcg/releases/LCG_92/xrootd/4.7.0/x86_64-slc6-gcc62-opt/lib64
  • Now copy in the code patches you wish to use (i.e. for HistFactory)
  • Make a new directory into which ROOT will be compiled (i.e. ~/myROOT) and cd into it.
  • Configure cmake
    • cmake -DAll=ON ~/root_source
  • Now you need to alter CMakeCache.txt so that it picks up the correct libraries for some things
    • //Path to a library.
      XROOTD_XrdCl_LIBRARY:FILEPATH=/cvmfs/lhcb.cern.ch/lib/lcg/releases/LCG_92/xrootd/4.7.0/x86_64-slc6-gcc62-opt/lib64/libXrdCl.so

      //Path to a library.
      XROOTD_XrdClient_LIBRARY:FILEPATH=/cvmfs/lhcb.cern.ch/lib/lcg/releases/LCG_92/xrootd/4.7.0/x86_64-slc6-gcc62-opt/lib64/libXrdClient.so

      //Path to a library.
      XROOTD_XrdMain_LIBRARY:FILEPATH=XROOTD_XrdMain_LIBRARY-NOTFOUND

      //Path to a library.
      XROOTD_XrdUtils_LIBRARY:FILEPATH=/cvmfs/lhcb.cern.ch/lib/lcg/releases/LCG_92/xrootd/4.7.0/x86_64-slc6-gcc62-opt/lib64/libXrdUtils.so

      //Dependencies for the target
      XrdProofd_LIB_DEPENDS:STATIC=general;/cvmfs/lhcb.cern.ch/lib/lcg/releases/LCG_92/xrootd/4.7.0/x86_64-slc6-gcc62-opt/lib64/libXrdUtils.so;general;/cvmfs/lhcb.cern.ch/lib/lcg/releases/LCG_92/xrootd/4.7.0/x86_64-slc6-gcc62-opt/lib64/libXrdClient.so;general;/cvmfs/lhcb.cern.ch/lib/lcg/releases/LCG_92/xrootd/4.7.0/x86_64-slc6-gcc62-opt/lib64/libXrdCl.so;general;rpdutil;general;-lpthread;general;dl;
  • Now compile!
    • cmake --build .
  • At some point the compilation may segfault or throw some other fit. Don't worry - just run the last command again and it should pick up where it left off.
  • Now make sure you are using your newly installed ROOT version:
    • source ~/myROOT/bin/thisroot.csh
Edit | Attach | Watch | Print version | History: r28 | r15 < r14 < r13 < r12 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r13 - 2019-02-05 - MarkElliotSmith
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main 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.
Ideas, requests, problems regarding TWiki? Send feedback