MuonSmearer.h

Go to the documentation of this file.

// MuonSmearer class
//
// Implements smearing for muons
//
// Namespace Atlfast::
//
//--------------------------------------------------------------
//
// MuonSmearer replicates the smearing code of ATLFAST++'s
// MuonMaker. It provides the smear() method defined
// in the ISmearer interface. It privately inherits from
// the DefaultSmearer class to use its implementation of
// random number services etc.
//
// Authors: H.T.Phillips, P. Clarke, E. Richter-Was, P. Sherwood, R. Steward
//

#ifndef ATLFAST_MUONSMEARER_H
#define ATLFAST_MUONSMEARER_H


#include <utility>       //contains pair

#include "AtlfastEvent/ChargeService.h"

#include "AtlfastAlgs/ISmearer.h"
#include "AtlfastAlgs/DefaultSmearer.h"
#include "AtlfastAlgs/TrackSmearer.h"
#include "AtlfastAlgs/MuonSpectrometer.h"

#include "CLHEP/Vector/LorentzVector.h"

namespace Atlfast {

  using std::pair;

  class MuonSmearer : virtual public ISmearer, 
                      virtual private DefaultSmearer 
    {
    public:
      // Constructor called from Athena code - instantiates TrackSmearer
      MuonSmearer(const int aseed, const int lumi, const int keymuo, 
                  std::string &resolutionFile, MsgStream&);
      // Constructor that can be called from outside Athena - no TrackSmearer
      MuonSmearer(const int aseed, const int lumi, const int keymuo, 
                  std::string &resolutionFile);
      virtual ~MuonSmearer();
      
      //--------------------------------------------------------------------------
      // the only thing we actually need to provide is the smear() method
      //--------------------------------------------------------------------------
      virtual HepLorentzVector smear(const HepMC::GenParticle&);
      virtual HepLorentzVector smear(const HepLorentzVector& avec);

      //      double resolumu_(double* pt, double* eta, double* phi, double* resol);
      pair<double, double> resolMS(const HepLorentzVector& avec);
      pair<double, double> resolID(const HepMC::GenParticle& particle);
      double               resol(double sigmams, double sigmamuon,
                                 double sigmaid, double sigmatrack);
      
      virtual int setSmearParameters( const std::vector<double>& smearValues );
      virtual int setSmearParamSchema( const int smearSchema );
      
    private:
      
      void makeResolutionCalculator(std::string &resolutionFile);

      //--------------------------------------------------------------------------
      // luminosity option. Currently using an int with 1= low lumi
      // 2 = high lumi. This is not very pleasant, but leave for now
      //
      // flag:
      // 1   use muon sepectrometer alone
      // 2   use ID alone
      // >2  combinations of muon spectrometer and ID
      //--------------------------------------------------------------------------
      int   m_lumi;
      int   m_keymuo;
      
      std::vector<double> m_smearParams;
      int m_smearParamSchema;
      
      TrackSmearer*              m_tracksmearer;
      MuonResolutionCalculator*  m_muonResCalculator;
      ChargeService*             m_chargeService;
      
    };
  
} // end of namespace bracket

#endif

---