Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/Efficiency.cc

    ra1c9c16 r95aa610  
    1717 */
    1818
     19
    1920/** \class Efficiency
    2021 *
     
    3132#include "classes/DelphesFormula.h"
    3233
     34#include "ExRootAnalysis/ExRootResult.h"
     35#include "ExRootAnalysis/ExRootFilter.h"
    3336#include "ExRootAnalysis/ExRootClassifier.h"
    34 #include "ExRootAnalysis/ExRootFilter.h"
    35 #include "ExRootAnalysis/ExRootResult.h"
    3637
     38#include "TMath.h"
     39#include "TString.h"
     40#include "TFormula.h"
     41#include "TRandom3.h"
     42#include "TObjArray.h"
    3743#include "TDatabasePDG.h"
    38 #include "TFormula.h"
    3944#include "TLorentzVector.h"
    40 #include "TMath.h"
    41 #include "TObjArray.h"
    42 #include "TRandom3.h"
    43 #include "TString.h"
    4445
    45 #include <algorithm>
     46#include <algorithm>
     47#include <stdexcept>
    4648#include <iostream>
    4749#include <sstream>
    48 #include <stdexcept>
    4950
    5051using namespace std;
     
    9394
    9495void Efficiency::Process()
    95 {
     96{ 
    9697  Candidate *candidate;
    97   Double_t pt, eta, phi, e, d0, dz, ctgTheta;
     98  Double_t pt, eta, phi, e;
    9899
    99100  fItInputArray->Reset();
    100   while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
     101  while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
    101102  {
    102103    const TLorentzVector &candidatePosition = candidate->Position;
     
    106107    pt = candidateMomentum.Pt();
    107108    e = candidateMomentum.E();
    108     d0 = candidate->D0;
    109     dz = candidate->DZ;
    110     ctgTheta = candidate->CtgTheta;
    111109
    112110    // apply an efficency formula
    113     if(gRandom->Uniform() > fFormula->Eval(pt, eta, phi, e, d0, dz, ctgTheta)) continue;
    114 
     111    if(gRandom->Uniform() > fFormula->Eval(pt, eta, phi, e)) continue;
     112   
    115113    fOutputArray->Add(candidate);
    116114  }
Note: See TracChangeset for help on using the changeset viewer.