Fork me on GitHub

Changeset 50d7259 in git


Ignore:
Timestamp:
Feb 27, 2015, 11:16:07 AM (10 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
76ab6be
Parents:
b55dc12
git-author:
Michele Selvaggi <michele.selvaggi@…> (02/27/15 11:11:12)
git-committer:
Michele Selvaggi <michele.selvaggi@…> (02/27/15 11:16:07)
Message:

corrected small bug in TauTaggingPartonClassifier responsible for auto-vetoing taus with explicit W daughter, moved TauTaggingPartonClassifier definition from TauTagging.cc to TauTagging.h

Location:
modules
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • modules/TauTagging.cc

    rb55dc12 r50d7259  
    3333#include "classes/DelphesFactory.h"
    3434#include "classes/DelphesFormula.h"
    35 
    36 #include "ExRootAnalysis/ExRootResult.h"
    37 #include "ExRootAnalysis/ExRootFilter.h"
    38 #include "ExRootAnalysis/ExRootClassifier.h"
    3935
    4036#include "TMath.h"
     
    5349using namespace std;
    5450
    55 //------------------------------------------------------------------------------
    56 
    57 class TauTaggingPartonClassifier : public ExRootClassifier
    58 {
    59 public:
    60 
    61   TauTaggingPartonClassifier(const TObjArray *array);
    62 
    63   Int_t GetCategory(TObject *object);
    64 
    65   Double_t fEtaMax, fPTMin;
    66 
    67   const TObjArray *fParticleInputArray;
    68 };
    6951
    7052//------------------------------------------------------------------------------
     
    7961{
    8062  Candidate *tau = static_cast<Candidate *>(object);
    81   Candidate *daughter = 0;
     63  Candidate *daughter1 = 0;
     64  Candidate *daughter2 = 0;
     65 
    8266  const TLorentzVector &momentum = tau->Momentum;
    83   Int_t pdgCode, i;
     67  Int_t pdgCode, i, j;
    8468
    8569  pdgCode = TMath::Abs(tau->PID);
     
    9882  for(i = tau->D1; i <= tau->D2; ++i)
    9983  {
    100     daughter = static_cast<Candidate *>(fParticleInputArray->At(i));
    101     pdgCode = TMath::Abs(daughter->PID);
    102     if(pdgCode == 11 || pdgCode == 13 || pdgCode == 15 || pdgCode == 24) return -1;
     84    daughter1 = static_cast<Candidate *>(fParticleInputArray->At(i));
     85    pdgCode = TMath::Abs(daughter1->PID);
     86    if(pdgCode == 11 || pdgCode == 13 || pdgCode == 15) return -1;
     87    else if(pdgCode == 24)
     88    {
     89     if(daughter1->D1 < 0) return -1;
     90     for(j = daughter1->D1; j <= daughter1->D2; ++j)
     91     {
     92       daughter2 = static_cast<Candidate*>(fParticleInputArray->At(j));
     93       pdgCode = TMath::Abs(daughter2->PID);
     94       if(pdgCode == 11 || pdgCode == 13) return -1;
     95     }
     96       
     97    }
    10398  }
    10499
  • modules/TauTagging.h

    rb55dc12 r50d7259  
    3131
    3232#include "classes/DelphesModule.h"
     33#include "ExRootAnalysis/ExRootResult.h"
     34#include "ExRootAnalysis/ExRootFilter.h"
     35#include "ExRootAnalysis/ExRootClassifier.h"
    3336
    3437#include <map>
     
    7679};
    7780
     81
     82//------------------------------------------------------------------------------
     83
     84class TauTaggingPartonClassifier : public ExRootClassifier
     85{
     86public:
     87
     88  TauTaggingPartonClassifier(const TObjArray *array);
     89
     90  Int_t GetCategory(TObject *object);
     91
     92  Double_t fEtaMax, fPTMin;
     93
     94  const TObjArray *fParticleInputArray;
     95};
     96
     97
    7898#endif
Note: See TracChangeset for help on using the changeset viewer.