BTaggIng < Main

Main Web>DerSchrecklicheSven>AthenaStuff>BTaggIng (2005-05-30, AndreasWildauer) (raw view)

* NOT VALID ANYMORE *

-----------

---++ How to get weights for the !LifetimeTag and !SecVtxTags

Basics: you are running 10.0.x (x>=1) and have checked out the 3 packages described here: 
<verbatim>
https://twiki.cern.ch/twiki/bin/view/Atlas/BTagging#Running_on_AOD
</verbatim>
which enables you to rerun b-tagging on AOD level.

Due to persistency problems the procedure to obtain weights has to be (almost) entirely shifted to your analysis algorithm.
The only files which had to be updated are:

   * !LikelihoodTool.cxx
   * !LikelihoodTool.h

Since I do not want to branch the packages and the HEAD versions are incompatible with 10.0.x already I did not commit any of these changes into the repository. Please copy by hand the following two files:
<verbatim>
/afs/cern.ch/user/w/wildauer/public/LikelihoodTool.cxx 
/afs/cern.ch/user/w/wildauer/public/LikelihoodTool.h
</verbatim>
to the apropriet locations inside !JetTagTools (src and !JetTagTools dir).

Now recompile the packages.

---+++ Changes in your analysis algorithm
For every tag tool you need a separate likelihood tool in your algorithm. Here is an example for the !LifetimeTag2D tool.
Assume the analysis algorithms are called !myAnalysis.h and !myAnalysis.cxx. In !myAnalysis.h add the following:
   * add a forward declaration
<verbatim>
class LikelihoodTool;
</verbatim>
   * and a datamember for the tool (one per tag, remember ... this on is for the 2D tag)
<verbatim>
LikelihoodTool* m_likelihoodTool;
</verbatim>
In the source files add:
   * the following inludes
<verbatim>
      #include "JetTagTools/ILikelihoodTool.h"
      #include "JetTagTools/LikelihoodTool.h"
</verbatim>
   * in the initialize() method get the tool:
<verbatim>
      /// retrieve the Likelihood tool
      ILikelihoodTool* likelihoodTool(0);
      StatusCode sc = toolSvc->retrieveTool("Analysis::LikelihoodTool", "AnalysisLikelihoodTool2D" , likelihoodTool);
      if (StatusCode::SUCCESS != sc && likelihoodTool == 0)
      {
          mLog << MSG::ERROR << "Can't get handle on AnalysisLikelihoodTool2D" << endreq;
          return StatusCode::FAILURE;
      }
      m_likelihoodTool = dynamic_cast< LikelihoodTool* >(likelihoodTool);
      if (m_likelihoodTool==0)
      {
          mLog << MSG::ERROR << "dynamic cast from ILikelihood to Likelihood did not work!" << endreq;
          return StatusCode::FAILURE;
      }
      m_likelihoodTool->readInHistosFromFile("JetTagRefFileSig/REF/t3333/");
      m_likelihoodTool->readInHistosFromFile("JetTagRefFileBkg/REF/t3333/");
      m_likelihoodTool->printStatus();
</verbatim>
   * in your execute() method add the following lines and you have the weight:
<verbatim>
      std::vector<double> signif = bLifeInfo2D->significance();
      m_likelihoodTool->setLhVariableValue("significance2D", signif);
      
      double weight2D(m_likelihoodTool->calculateWeight());
</verbatim>
*Important:* the variable !bLifeInfo2D is a pointer to the !LifetimeInfo object in the !JetTag object. If you do not have it already in the analysis let me know :-)

   * last step: modify your !jobOptions and add these lines:
<verbatim>
     AnalysisLikelihoodTool2D = Service ("ToolSvc.AnalysisLikelihoodTool2D")
     AnalysisLikelihoodTool2D.allLhVariables      = [ "significance2D" ]
     AnalysisLikelihoodTool2D.useTheseLhVariables = [ "significance2D" ]
</verbatim>
For the 1D and 3D weight repeat everything but change 2D correspondingly. For the !SecVtxTags it is a bit more work (more variables) which I will post in a bit. Please try this one first.

---++ Attached you find very preliminar plots to show that it works ...

-- Main.DerSchrecklicheSven - 23 May 2005

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
ps	Weight1D.ps	r1	manage	10.7 K	2005-05-23 - 09:26	AndreasWildauer	1D weight
ps	Weight2D.ps	r1	manage	10.9 K	2005-05-23 - 09:25	AndreasWildauer	2D weight
ps	Weight3D.ps	r1	manage	10.9 K	2005-05-23 - 09:27	AndreasWildauer	3D weight

Topic revision: r2 - 2005-05-30 - AndreasWildauer

Main

Webs

Public webs

Welcome Guest

- Cern Search
- TWiki Search
- Google Search
Main All webs

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