TauFinder.h

Go to the documentation of this file.

#ifndef TauFinder_h
#define TauFinder_h 1

#include "marlin/Processor.h"
#include "lcio.h"
#include <string>
#include "TTree.h"
#include "TNtuple.h"
#include "TH1D.h"
#include "TFile.h"
#include "TMath.h"
#include "TSystem.h"
#include <EVENT/MCParticle.h>
#include <EVENT/ReconstructedParticle.h>

using namespace lcio ;
using namespace marlin ;


/** TauFinder processor for marlin.
 * 
 * @author A. Muennich, CERN
 * 
 */

class TauFinder : public Processor {
  
 public:
  
  virtual Processor*  newProcessor() { return new TauFinder ; }
  
  
  TauFinder() ;
  TauFinder(const TauFinder&) = delete;
  TauFinder& operator=(const TauFinder&) = delete;

  /** Called at the begin of the job before anything is read.
   * Use to initialize the processor, e.g. book histograms.
   */
  virtual void init() ;
  
  /** Called for every run.
   */
  virtual void processRunHeader( LCRunHeader* run ) ;
  
  /** Called for every event - the working horse.
   */
  virtual void processEvent( LCEvent * evt ) ; 
  
  
  virtual void check( LCEvent * evt ) ; 
  
  
  /** Called after data processing for clean up.
   */
  virtual void end() ;
  
  
 protected:

  /** Input collection name.
   */
  std::string _colNameMC{}, _colNameRECO{}, _incol{};
  std::string _colNameMCTruth{}, _colNameTauRecLink{};
  std::string _outcol{}, _outcolRest{};
  std::string _OutputFile_Signal{};
  int _nRun=-1;
  int _nEvt=-1;

  float _ptcut=0.0,_ptseed=0.0,_cosTcut=0.0;
  float _coneAngle=0.0,_isoAngle=0.0,_isoE=0.0;
  float _D0seedmin=0.0, _D0seedmax=0.0,_minv=0.0;

  int _fail_minv=0.0,_fail_minv_neg=0.0,_fail_Qtr=0.0,_fail_isoE=0.0,_mergeTries=0.0;
  
  TFile *rootfile=NULL;
  TNtuple *failtuple=NULL;
  
  bool FindTau(std::vector<ReconstructedParticle*> &Qvec,std::vector<ReconstructedParticle*> &Nvec,
           std::vector<std::vector<ReconstructedParticle*> > &tauvec);
  
} ;

#endif