/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/FastSimulation/CalorimeterProperties/src/PreshowerLayer1Properties.cc

Go to the documentation of this file.

#include "FWCore/ParameterSet/interface/ParameterSet.h"
//This class header
#include "FastSimulation/CalorimeterProperties/interface/PreshowerLayer1Properties.h"
#include <cmath>
PreshowerLayer1Properties::PreshowerLayer1Properties(const edm::ParameterSet& fastDet) 
  : PreshowerProperties()
{
  // Preshower : mumber of Mips / GeV 
  mips = fastDet.getParameter<double>("PreshowerLayer1_mipsPerGeV");
  thick = fastDet.getParameter<double>("PreshowerLayer1_thickness");
}

double PreshowerLayer1Properties::thickness(double eta) const { 

  // eta is the pseudorapidity
  double e  = exp(-eta);
  double e2 = e*e;
  // 1 / cos theta
  double cinv  = (1.+e2)/(1.-e2);
  //    double c  = (1.-e2)/(1.+e2);
  //    double s  = 2.*e/(1.+e2);
  //    double t  = 2.*e/(1.-e2);
  double feta = fabs(eta);
  
  if ( 1.623 < feta && feta < 2.611 ) 
    {
      return thick * fabs(cinv); 
    }
  else
    {
      return 0;
    }
  
}