Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/IdentificationMap.cc

    r95aa610 r341014c  
    1717 */
    1818
    19 
    2019/** \class IdentificationMap
    2120 *
     
    3332#include "classes/DelphesFormula.h"
    3433
     34#include "ExRootAnalysis/ExRootClassifier.h"
     35#include "ExRootAnalysis/ExRootFilter.h"
    3536#include "ExRootAnalysis/ExRootResult.h"
    36 #include "ExRootAnalysis/ExRootFilter.h"
    37 #include "ExRootAnalysis/ExRootClassifier.h"
    3837
     38#include "TDatabasePDG.h"
     39#include "TFormula.h"
     40#include "TLorentzVector.h"
    3941#include "TMath.h"
     42#include "TObjArray.h"
     43#include "TRandom3.h"
    4044#include "TString.h"
    41 #include "TFormula.h"
    42 #include "TRandom3.h"
    43 #include "TObjArray.h"
    44 #include "TDatabasePDG.h"
    45 #include "TLorentzVector.h"
    4645
    4746#include <algorithm>
    48 #include <stdexcept>
    4947#include <iostream>
    5048#include <sstream>
     49#include <stdexcept>
    5150
    5251using namespace std;
     
    7978
    8079  fEfficiencyMap.clear();
    81   for(i = 0; i < size/3; ++i)
     80  for(i = 0; i < size / 3; ++i)
    8281  {
    8382    formula = new DelphesFormula;
    84     formula->Compile(param[i*3 + 2].GetString());
    85     pdg = param[i*3].GetInt();
    86     fEfficiencyMap.insert(make_pair(pdg, make_pair(param[i*3 + 1].GetInt(), formula)));
     83    formula->Compile(param[i * 3 + 2].GetString());
     84    pdg = param[i * 3].GetInt();
     85    fEfficiencyMap.insert(make_pair(pdg, make_pair(param[i * 3 + 1].GetInt(), formula)));
    8786  }
    8887
     
    129128  Double_t pt, eta, phi, e;
    130129  TMisIDMap::iterator itEfficiencyMap;
    131   pair <TMisIDMap::iterator, TMisIDMap::iterator> range;
     130  pair<TMisIDMap::iterator, TMisIDMap::iterator> range;
    132131  DelphesFormula *formula;
    133132  Int_t pdgCodeIn, pdgCodeOut, charge;
     
    136135
    137136  fItInputArray->Reset();
    138   while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
     137  while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
    139138  {
    140139    const TLorentzVector &candidatePosition = candidate->Position;
     
    144143    pt = candidateMomentum.Pt();
    145144    e = candidateMomentum.E();
    146    
     145
    147146    pdgCodeIn = candidate->PID;
    148147    charge = candidate->Charge;
     
    171170      {
    172171        // change PID of particle
    173         if(pdgCodeOut != 0) candidate->PID = charge*pdgCodeOut;
     172        if(pdgCodeOut != 0) candidate->PID = charge * pdgCodeOut;
    174173        fOutputArray->Add(candidate);
    175174        break;
Note: See TracChangeset for help on using the changeset viewer.