![]() |
![]() |
00001 #include "FWCore/ParameterSet/interface/ParameterSet.h" 00002 00003 //This class header 00004 #include "FastSimulation/CalorimeterProperties/interface/PreshowerLayer2Properties.h" 00005 #include <cmath> 00006 PreshowerLayer2Properties::PreshowerLayer2Properties(const edm::ParameterSet& fastDet) : PreshowerProperties() 00007 { 00008 // Preshower : mumber of Mips / GeV 00009 mips = fastDet.getParameter<double>("PreshowerLayer2_mipsPerGeV"); 00010 thick = fastDet.getParameter<double>("PreshowerLayer2_thickness"); 00011 } 00012 00013 double PreshowerLayer2Properties::thickness(const double eta) const { 00014 00015 // eta is the pseudorapidity 00016 double e = exp(-eta); 00017 double e2 = e*e; 00018 // 1 / cos theta 00019 double cinv = (1.+e2)/(1.-e2); 00020 // double c = (1.-e2)/(1.+e2); 00021 // double s = 2.*e/(1.+e2); 00022 // double t = 2.*e/(1.-e2); 00023 double feta = fabs(eta); 00024 00025 if ( 1.637 < feta && feta < 2.625 ) 00026 { 00027 return thick * fabs(cinv); 00028 } 00029 else 00030 { 00031 return 0; 00032 } 00033 00034 }