---

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 <pair.h>
#include "AtlfastCode/ISmearer.h"
#include "AtlfastCode/DefaultSmearer.h"

class HepLorentzVector; //forward declaration

namespace Atlfast {
  class MuonSmearer : virtual public ISmearer, 
                          virtual private DefaultSmearer 
    {
    public:
      //--------------------------------------------------------------------------
      // constructors and destructors- we can leave everything to DefaultSmearer
      //--------------------------------------------------------------------------
      MuonSmearer(const int aseed, const int lumi, const int keymuo): 
        ISmearer(),       DefaultSmearer(aseed), m_lumi(lumi),    m_keymuo(keymuo), 
        m_magic1(0.0005), m_magic2(10./7000.),   m_magic3(0.012), m_magic4(0.02),
        m_percent(0.01) { }
      virtual ~MuonSmearer() { };

      //--------------------------------------------------------------------------
      // the only thing we actually need to provide is the smear() method
      //--------------------------------------------------------------------------
      virtual HepLorentzVector smear(const HepLorentzVector& avec);

    private:
      //      float resolumu_(float* pt, float* eta, float* phi, float* resol);
      pair<float, float> resolms(float pt, float eta, float phi);
      float resolid1(float pt, float eta);
      pair<float, float> resolid2(float pt, float eta);
      float resol1 (float sigmams, float sigmamuon,
                    float sigmaid, float sigmatrack);
      float resol2 (float sigmamuon, float sigmatrack);
      float resol3  ();
      float resol4  ();


      //--------------------------------------------------------------------------
      // 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;
      float m_magic1;
      float m_magic2;
      float m_magic3;
      float m_magic4;
      float m_percent;  //to convert resolumu resolutions which are percerntages
    };

} // end of namespace bracket

#endif

---