FMGLfromTKEfficiencyHandler Class Reference

#include <FMGLfromTKEfficiencyHandler.h>

List of all members.

Public Member Functions
	FMGLfromTKEfficiencyHandler (const RandomEngine *engine)
bool	kill (const SimTrack &)
	~FMGLfromTKEfficiencyHandler ()
Private Attributes
float	Effic_Eta [nEtaBins]
const RandomEngine *	random
Static Private Attributes
static const int	nEtaBins = 60

---

Detailed Description

Definition at line 13 of file FMGLfromTKEfficiencyHandler.h.

---

Constructor & Destructor Documentation

FMGLfromTKEfficiencyHandler::FMGLfromTKEfficiencyHandler

(

const RandomEngine *

engine

)

Definition at line 11 of file FMGLfromTKEfficiencyHandler.cc.

References gather_cfg::cout, Effic_Eta, connectstrParser::f1, alignmentValidation::fname, i, FileInPath::name(), nEtaBins, and path().

  : random(engine) {;

  std::string fname = "FastSimulation/ParamL3MuonProducer/data/efficiencyGL_TK.data";
  std::string path(getenv("CMSSW_SEARCH_PATH"));
  FileInPath f1(path,fname);
  
  if ( f1() == 0) {
    std::cout << "File efficiencyGL_TK.data not found in " << path << std::endl;
    throw Genexception(" efficiency list not found for FMGLfromTKEfficiencyHandler.");
  } else {
    std::cout << "Reading " << f1.name() << std::endl;
  }
  std::ifstream & listfile = *f1();
  
  double eff=0.;
  int nent;
  listfile >> nent;
  if (nent != nEtaBins) { 
    std::cout << " *** ERROR -> FMGLfromTKEfficiencyHandler : nEta bins " << nent << " instead of " << nEtaBins << std::endl;
  }
  for (int i=0; i<nEtaBins; i++) {
    listfile >> eff;
    Effic_Eta[i]=eff;
  }

}

FMGLfromTKEfficiencyHandler::~FMGLfromTKEfficiencyHandler

(

)

Definition at line 39 of file FMGLfromTKEfficiencyHandler.cc.

References Effic_Eta.

                                                         {
  delete Effic_Eta;
}

---

Member Function Documentation

bool FMGLfromTKEfficiencyHandler::kill

(

const SimTrack &

aTrack

)

Definition at line 43 of file FMGLfromTKEfficiencyHandler.cc.

References Effic_Eta, eta(), funct::exp(), RandomEngine::flatShoot(), CoreSimTrack::momentum(), nEtaBins, ExpressReco_HICollisions_FallBack::pt, random, and mathSSE::sqrt().

Referenced by ParamL3MuonProducer::produce().

                                                              {
  
  // At least eight hit in the tracker : To be tuned !!!
  //if ( aTrack.recHits().size() < 8 ) return false;

  // At least zero reconstructed  Pixel hits : To be tuned !!!
  //int seed = 0;
  //for ( unsigned i=1; i<6; ++i ) 
  //  if ( aTrack.isARecHit(i) ) ++seed;
  //if ( seed < 0 ) return false;

  double myEffEta=0. , myEffPt=0. , myCorrection=1. , myEff;

  // Eta dependence : 
  double eta = fabs(aTrack.momentum().eta());
  if (eta < 2.40) {
    int iEtaBin = (int) ( (eta/2.40) * nEtaBins);
    myEffEta = Effic_Eta[iEtaBin];
  } else return false;

  // Pt dependence (same as L3 muons for "track-only" Global Muons):
  double pt = std::sqrt(aTrack.momentum().perp2());
  if (eta<1.04) {
    if (pt>40.) myEffPt = 0.9583 - pt*5.82726e-05;
    //    else if (pt>4.10) myEffPt = 0.952*(1.-exp(-(pt-4.072)));
    else if (pt>4.0) myEffPt = 0.952*(1.-sqrt(exp(-(pt-3.2))));
    //    myCorrection = 1.045;
    myCorrection = 1.077;
  }
  else if (eta<2.07) {
    if (pt>173.) myEffPt = 0.991 - pt*3.46562e-05;
    //    else if (pt>3.10) myEffPt = 0.985*(1.-exp(-(pt-3.061)));
    else if (pt>3.0) myEffPt = 0.985*(1.-sqrt(exp(-(pt-2.4))));
    //    myCorrection = 1.027;
    myCorrection = 1.028;
  }
  else if (eta<2.40) {
    if (pt>26.) myEffPt = 0.9221 - pt*7.75139e-05;
    else if (pt>3.00) myEffPt = 0.927*(1.-exp(-sqrt(pt-1.617)));
    //    myCorrection = 1.061;
    myCorrection = 1.133;
  }
  else return false;

  /*
  // Special high Pt muons treatment :
  if (pt>400.) {
    float myEffHighPt=1.;
    float pttev = pt/1000.;
    if (eta < 0.3) {
      myEffHighPt = 0.952 - 0.033*pttev;
    } else if ( eta < 0.6) {
      myEffHighPt = 0.973 - 0.033*pttev;
    } else if ( eta < 0.9) {
      myEffHighPt = 0.969 - 0.045*pttev;
    } else if ( eta < 1.2) {
      myEffHighPt = 0.968 - 0.058*pttev;
    } else if ( eta < 1.5) {
      myEffHighPt = 0.966 - 0.074*pttev;
    } else if ( eta < 1.8) {
      myEffHighPt = 0.955 - 0.11*pttev;
    } else if ( eta < 2.1) {
      myEffHighPt = 0.982 - 0.11*pttev;
    } else {
      myEffHighPt = 0.958 - 0.14*pttev;
    }
    myEffPt = (myEffPt>myEffHighPt? myEffHighPt : myEffPt);
  }
  */
  /*
  if (eta<1.04) {
    if (pt>40.) myEffPt = 0.9583 - pt*5.82726e-05;
    else if (pt>4.10) myEffPt = 0.952*(1.-exp(-(pt-4.072)));
    myCorrection = 1.115;
  }
  else if (eta<2.07) {
    if (pt>173.) myEffPt = 0.991 - pt*3.46562e-05;
    else if (pt>3.10) myEffPt = 0.985*(1.-exp(-(pt-3.061)));
    myCorrection = 1.034;
  }
  else if (eta<2.40) {
    if (pt>26.) myEffPt = 0.9221 - pt*7.75139e-05;
    else if (pt>3.00) myEffPt = 0.927*(1.-exp(-sqrt(pt-1.617)));
    myCorrection = 1.157;
  }
  else return false;
  */

  myEff = myEffEta*myEffPt*myCorrection;

  double prob = random->flatShoot();
  return (myEff > prob) ;

}

---

Member Data Documentation

float FMGLfromTKEfficiencyHandler::Effic_Eta[nEtaBins] [private]

Definition at line 25 of file FMGLfromTKEfficiencyHandler.h.

Referenced by FMGLfromTKEfficiencyHandler(), kill(), and ~FMGLfromTKEfficiencyHandler().

const int FMGLfromTKEfficiencyHandler::nEtaBins = 60 [static, private]

Definition at line 24 of file FMGLfromTKEfficiencyHandler.h.

Referenced by FMGLfromTKEfficiencyHandler(), and kill().

const RandomEngine* FMGLfromTKEfficiencyHandler::random [private]

Definition at line 23 of file FMGLfromTKEfficiencyHandler.h.

Referenced by kill().