Fork me on GitHub

Ignore:
Timestamp:
Feb 12, 2019, 9:29:17 PM (6 years ago)
Author:
Pavel Demin <pavel-demin@…>
Branches:
ImprovedOutputFile, Timing, llp, master
Children:
6455202
Parents:
45e58be
Message:

apply .clang-format to all .h, .cc and .cpp files

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/TaggingParticlesSkimmer.cc

    r45e58be r341014c  
    3535#include "classes/DelphesFormula.h"
    3636
     37#include "ExRootAnalysis/ExRootClassifier.h"
     38#include "ExRootAnalysis/ExRootFilter.h"
    3739#include "ExRootAnalysis/ExRootResult.h"
    38 #include "ExRootAnalysis/ExRootFilter.h"
    39 #include "ExRootAnalysis/ExRootClassifier.h"
    4040
     41#include "TDatabasePDG.h"
     42#include "TFormula.h"
     43#include "TLorentzVector.h"
    4144#include "TMath.h"
     45#include "TObjArray.h"
     46#include "TRandom3.h"
    4247#include "TString.h"
    43 #include "TFormula.h"
    44 #include "TRandom3.h"
    45 #include "TObjArray.h"
    46 #include "TDatabasePDG.h"
    47 #include "TLorentzVector.h"
    4848
    4949#include <algorithm>
    50 #include <stdexcept>
    5150#include <iostream>
    5251#include <sstream>
     52#include <stdexcept>
    5353
    5454using namespace std;
    55 
    5655
    5756//------------------------------------------------------------------------------
     
    7271void TaggingParticlesSkimmer::Init()
    7372{
    74  
     73
    7574  fPTMin = GetDouble("PTMin", 15.0);
    7675  fEtaMax = GetDouble("EtaMax", 2.5);
    77    
     76
    7877  // import input array
    7978  fPartonInputArray = ImportArray(GetString("PartonInputArray", "Delphes/partons"));
    8079  fItPartonInputArray = fPartonInputArray->MakeIterator();
    81  
     80
    8281  fParticleInputArray = ImportArray(GetString("ParticleInputArray", "Delphes/allParticles"));
    8382
    84  
    8583  fClassifier = new TauTaggingPartonClassifier(fParticleInputArray);
    8684  fClassifier->fPTMin = GetDouble("PTMin", 15.0);
    8785  fClassifier->fEtaMax = GetDouble("EtaMax", 2.5);
    8886
    89  
    9087  fFilter = new ExRootFilter(fPartonInputArray);
    9188
     
    9794
    9895void TaggingParticlesSkimmer::Finish()
    99 { 
     96{
    10097  if(fItPartonInputArray) delete fItPartonInputArray;
    10198  if(fFilter) delete fFilter;
     
    127124    if(tau->D1 < 0) continue;
    128125
    129     if(tau->D1 >= fParticleInputArray->GetEntriesFast() ||
    130        tau->D2 >= fParticleInputArray->GetEntriesFast())
     126    if(tau->D1 >= fParticleInputArray->GetEntriesFast() || tau->D2 >= fParticleInputArray->GetEntriesFast())
    131127    {
    132128      throw runtime_error("tau's daughter index is greater than the ParticleInputArray size");
     
    134130
    135131    tauMomentum.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
    136    
     132
    137133    for(i = tau->D1; i <= tau->D2; ++i)
    138134    {
     
    141137      tauMomentum += daughter->Momentum;
    142138    }
    143  
    144    candidate = static_cast<Candidate*>(tau->Clone());
    145    candidate->Momentum = tauMomentum;
    146139
    147    
    148    fOutputArray->Add(candidate);
     140    candidate = static_cast<Candidate *>(tau->Clone());
     141    candidate->Momentum = tauMomentum;
    149142
     143    fOutputArray->Add(candidate);
    150144  }
    151145
    152146  // then add all other partons (except tau's to avoid double counting)
    153  
     147
    154148  fItPartonInputArray->Reset();
    155   while((candidate = static_cast<Candidate*>(fItPartonInputArray->Next())))
     149  while((candidate = static_cast<Candidate *>(fItPartonInputArray->Next())))
    156150  {
    157151    pdgCode = TMath::Abs(candidate->PID);
    158152    if(pdgCode == 15) continue;
    159    
     153
    160154    pt = candidate->Momentum.Pt();
    161155    if(pt < fPTMin) continue;
    162    
     156
    163157    eta = TMath::Abs(candidate->Momentum.Eta());
    164158    if(eta > fEtaMax) continue;
    165        
     159
    166160    fOutputArray->Add(candidate);
    167161  }
    168 
    169 
    170162}
    171 
Note: See TracChangeset for help on using the changeset viewer.