CellSmearer.cxx

Go to the documentation of this file.

// CellSmearer.cxx
//
// Implementation of the Cell Smearer class
//
// Only the smear() method is overridden
// Based on JetSmearer, lowlumi, 0.1*0.1
//
// Authors: J.Couchman
//
//

#include "AtlfastAlgs/CellSmearer.h"
#include <cmath>
#include <iostream>
#include "CLHEP/Vector/LorentzVector.h"
#include "CLHEP/Random/JamesRandom.h"
#include "CLHEP/Random/RandGauss.h"
#include "CLHEP/Units/SystemOfUnits.h"


namespace Atlfast {

  HepLorentzVector CellSmearer::smear(const HepMC::GenParticle& particle){
    return smear(particle.momentum());
  }

  HepLorentzVector CellSmearer::smear (const HepLorentzVector& vec) {
//........................................................
//.....smear Cells energy
//........................................................
//     parametrizes smearing for hadronic energy deposition
//     parametrization from L. Poggioli
//     no pileup added

    float sigma=0.;
    HepLorentzVector smearedVec(vec);
    float aa, bb;
    float sqrtene = sqrt(vec.e()/GeV);
    float abseta  = fabs(vec.pseudoRapidity());
    while(1) {
      aa=randGauss()->fire();
      bb=randGauss()->fire();
      if(abseta < m_BarrelForwardEta) sigma = aa*0.5/sqrtene + bb*0.03;
      else            sigma = aa*1.0/sqrtene + bb*0.07;
      if(1.+sigma > .0) break;
    }
    smearedVec.setPx(vec.px()*(1.0+sigma));
    smearedVec.setPy(vec.py()*(1.0+sigma));
    smearedVec.setPz(vec.pz()*(1.0+sigma));
    smearedVec.setE(vec.e()*(1.0+sigma));
    return smearedVec;
    
  }

  //cct: implement the setSmearParameters method
  int CellSmearer::setSmearParameters (const std::vector<double>& /*smearValues*/){
    return 0;
  }
  //cct: implement the setSmearParamSchema method
  int CellSmearer::setSmearParamSchema ( const int /*smearSchema*/){
    return 0;
  }
  
}//end of namespace bracket

---