Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/TimeSmearing.cc

    r341014c r043e5b3  
    1717 */
    1818
     19
    1920/** \class TimeSmearing
    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
    4546#include <algorithm>
     47#include <stdexcept>
    4648#include <iostream>
    4749#include <sstream>
    48 #include <stdexcept>
    4950
    5051using namespace std;
     
    5354
    5455TimeSmearing::TimeSmearing() :
    55   fItInputArray(0)
     56fItInputArray(0)
    5657{
    5758}
     
    9697
    9798  fItInputArray->Reset();
    98   while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
     99  while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
    99100  {
    100101    const TLorentzVector &candidateInitialPosition = candidate->InitialPosition;
    101102    const TLorentzVector &candidateFinalPosition = candidate->Position;
    102103
    103     ti = candidateInitialPosition.T() * 1.0E-3 / c_light;
    104     tf = candidateFinalPosition.T() * 1.0E-3 / c_light;
     104    ti = candidateInitialPosition.T()*1.0E-3/c_light;
     105    tf = candidateFinalPosition.T()*1.0E-3/c_light;
    105106
    106107    // apply smearing formula
     
    108109    ti = ti + tf_smeared - tf;
    109110    tf = tf_smeared;
     111   
     112    mother = candidate;
     113    candidate = static_cast<Candidate*>(candidate->Clone());
     114    candidate->InitialPosition.SetT(ti*1.0E3*c_light);
     115    candidate->Position.SetT(tf*1.0E3*c_light);
    110116
    111     mother = candidate;
    112     candidate = static_cast<Candidate *>(candidate->Clone());
    113     candidate->InitialPosition.SetT(ti * 1.0E3 * c_light);
    114     candidate->Position.SetT(tf * 1.0E3 * c_light);
    115 
    116     candidate->ErrorT = fTimeResolution * 1.0E3 * c_light;
     117    candidate->ErrorT = fTimeResolution*1.0E3*c_light;
    117118
    118119    candidate->AddCandidate(mother);
Note: See TracChangeset for help on using the changeset viewer.