CMS 3D CMS Logo

HcalTB02SD Class Reference

Description: Stores hits of Test Beam 2002 calorimeters. More...

#include <SimG4CMS/HcalTestBeam/interface/HcalTB02SD.h>

Inheritance diagram for HcalTB02SD:

CaloSD SensitiveCaloDetector Observer< const BeginOfRun * > Observer< const BeginOfEvent * > Observer< const EndOfTrack * > Observer< const EndOfEvent * > SensitiveDetector

List of all members.

Public Member Functions

virtual double getEnergyDeposit (G4Step *)
 HcalTB02SD (G4String, const DDCompactView &, SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
virtual uint32_t setDetUnitId (G4Step *step)
void setNumberingScheme (HcalTB02NumberingScheme *scheme)
virtual ~HcalTB02SD ()

Private Member Functions

double crystalLength (G4String)
double curve_LY (G4String &, G4StepPoint *)
void initMap (G4String, const DDCompactView &)

Private Attributes

double birk1
double birk2
double birk3
std::map< G4String, double > lengthMap
HcalTB02NumberingScheme * numberingScheme
bool useBirk
bool useWeight

Detailed Description

Description: Stores hits of Test Beam 2002 calorimeters.

Usage: Activation is done using the XML file by choosing HcalTB02SensitiveDetector

Definition at line 30 of file HcalTB02SD.h.

Constructor & Destructor Documentation

HcalTB02SD::HcalTB02SD ( G4String  name,
const DDCompactView &  cpv,
SensitiveDetectorCatalog &  clg,
edm::ParameterSet const &  p,
const SimTrackManager *  manager 
)

Definition at line 35 of file HcalTB02SD.cc.

References birk1, birk2, birk3, g, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), initMap(), LogDebug, setNumberingScheme(), useBirk, and useWeight.

00038                                                        : 
00039   CaloSD(name, cpv, clg, p, manager), numberingScheme(0) {
00040   
00041   edm::ParameterSet m_SD = p.getParameter<edm::ParameterSet>("HcalTB02SD");
00042   useBirk= m_SD.getUntrackedParameter<bool>("UseBirkLaw",false);
00043   birk1  = m_SD.getUntrackedParameter<double>("BirkC1",0.013)*(g/(MeV*cm2));
00044   birk2  = m_SD.getUntrackedParameter<double>("BirkC2",0.0568);
00045   birk3  = m_SD.getUntrackedParameter<double>("BirkC3",1.75);
00046   useWeight= true;
00047 
00048   HcalTB02NumberingScheme* scheme=0;
00049   if      (name == "EcalHitsEB") {
00050     scheme = dynamic_cast<HcalTB02NumberingScheme*>(new HcalTB02XtalNumberingScheme());
00051     useBirk = false;
00052   } else if (name == "HcalHits") {
00053     scheme = dynamic_cast<HcalTB02NumberingScheme*>(new HcalTB02HcalNumberingScheme());
00054       useWeight= false;
00055   } else {edm::LogWarning("HcalTBSim") << "HcalTB02SD: ReadoutName " << name
00056                                        << " not supported\n";}
00057 
00058   if (scheme)  setNumberingScheme(scheme);
00059   LogDebug("HcalTBSim") 
00060     << "***************************************************" 
00061     << "\n"
00062     << "*                                                 *" 
00063     << "\n"
00064     << "* Constructing a HcalTB02SD  with name " << GetName()
00065     << "\n"
00066     << "*                                                 *"
00067     << "\n"
00068     << "***************************************************" ;
00069   edm::LogInfo("HcalTBSim")  << "HcalTB02SD:: Use of Birks law is set to      "
00070                              << useBirk << "        with three constants kB = "
00071                              << birk1 << ", C1 = " << birk2 << ", C2 = "
00072                              << birk3;
00073 
00074   if (useWeight) initMap(name,cpv);
00075 
00076 }

HcalTB02SD::~HcalTB02SD (  )  [virtual]

Definition at line 78 of file HcalTB02SD.cc.

References numberingScheme.

00078                         {
00079   if (numberingScheme) delete numberingScheme;
00080 }

Member Function Documentation

double HcalTB02SD::crystalLength ( G4String  name  )  [private]

Definition at line 178 of file HcalTB02SD.cc.

References it, and lengthMap.

Referenced by curve_LY().

00178                                               {
00179 
00180   double length = 230.;
00181   std::map<G4String,double>::const_iterator it = lengthMap.find(name);
00182   if (it != lengthMap.end()) length = it->second;
00183   return length;
00184 }

double HcalTB02SD::curve_LY ( G4String &  nameVolume,
G4StepPoint *  stepPoint 
) [private]

Definition at line 153 of file HcalTB02SD.cc.

References crystalLength(), LogDebug, CaloSD::setToLocal(), and weight.

Referenced by getEnergyDeposit().

00153                                                                         {
00154 
00155   double weight = 1.;
00156   G4ThreeVector  localPoint = setToLocal(stepPoint->GetPosition(),
00157                                          stepPoint->GetTouchable());
00158   double crlength = crystalLength(nameVolume);
00159   double dapd = 0.5 * crlength - localPoint.z();
00160   if (dapd >= -0.1 || dapd <= crlength+0.1) {
00161     if (dapd <= 100.)
00162       weight = 1.05 - dapd * 0.0005;
00163   } else {
00164     edm::LogWarning("HcalTBSim") << "HcalTB02SD: light coll curve : wrong "
00165                                  << "distance to APD " << dapd <<" crlength = "
00166                                  << crlength << " crystal name = " <<nameVolume
00167                                  << " z of localPoint = " << localPoint.z() 
00168                                  << " take weight = " << weight;
00169   }
00170   LogDebug("HcalTBSim") << "HcalTB02SD, light coll curve : " << dapd 
00171                         << " crlength = " << crlength
00172                         << " crystal name = " << nameVolume 
00173                         << " z of localPoint = " << localPoint.z() 
00174                         << " take weight = " << weight;
00175   return weight;
00176 }

double HcalTB02SD::getEnergyDeposit ( G4Step *  aStep  )  [virtual]

Reimplemented from CaloSD.

Definition at line 86 of file HcalTB02SD.cc.

References birk1, birk2, birk3, curve_LY(), CaloSD::getAttenuation(), LogDebug, NULL, CaloSD::preStepPoint, useBirk, useWeight, and weight.

00086                                                   {
00087   
00088   if (aStep == NULL) {
00089     return 0;
00090   } else {
00091     preStepPoint        = aStep->GetPreStepPoint();
00092     G4String nameVolume = preStepPoint->GetPhysicalVolume()->GetName();
00093 
00094     // take into account light collection curve for crystals
00095     double weight = 1.;
00096     if (useWeight) weight *= curve_LY(nameVolume, preStepPoint);
00097     if (useBirk)   weight *= getAttenuation(aStep, birk1, birk2, birk3);
00098     double edep   = aStep->GetTotalEnergyDeposit() * weight;
00099     LogDebug("HcalTBSim") << "HcalTB02SD:: " << nameVolume
00100                           <<" Light Collection Efficiency " << weight 
00101                           << " Weighted Energy Deposit " << edep/MeV << " MeV";
00102     return edep;
00103   } 
00104 }

void HcalTB02SD::initMap ( G4String  sd,
const DDCompactView &  cpv 
) [private]

Definition at line 119 of file HcalTB02SD.cc.

References DDFilteredView::addFilter(), DDSplit(), ddtrap, DDSpecificsFilter::equals, filter, first, DDFilteredView::firstChild(), i, it, lengthMap, LogDebug, DDFilteredView::logicalPart(), SensitiveDetector::name, DDBase< N, C >::name(), DDFilteredView::next(), DDSolid::parameters(), DDSpecificsFilter::setCriteria(), DDSolid::shape(), and DDLogicalPart::solid().

Referenced by HcalTB02SD().

00119                                                                {
00120 
00121   G4String attribute = "ReadOutName";
00122   DDSpecificsFilter filter;
00123   DDValue           ddv(attribute,sd,0);
00124   filter.setCriteria(ddv,DDSpecificsFilter::equals);
00125   DDFilteredView fv(cpv);
00126   fv.addFilter(filter);
00127   fv.firstChild();
00128 
00129   bool dodet=true;
00130   while (dodet) {
00131     const DDSolid & sol  = fv.logicalPart().solid();
00132     const std::vector<double> & paras = sol.parameters();
00133     G4String name = DDSplit(sol.name()).first;
00134     LogDebug("HcalTBSim") << "HcalTB02SD::initMap (for " << sd << "): Solid " 
00135                           << name << " Shape " << sol.shape() 
00136                           << " Parameter 0 = " << paras[0];
00137     if (sol.shape() == ddtrap) {
00138       double dz = 2*paras[0];
00139       lengthMap.insert(std::pair<G4String,double>(name,dz));
00140     }
00141     dodet = fv.next();
00142   }
00143   LogDebug("HcalTBSim") << "HcalTB02SD: Length Table for " << attribute 
00144                         << " = " << sd << ":";   
00145   std::map<G4String,double>::const_iterator it = lengthMap.begin();
00146   int i=0;
00147   for (; it != lengthMap.end(); it++, i++) {
00148     LogDebug("HcalTBSim") << " " << i << " " << it->first << " L = " 
00149                           << it->second;
00150   }
00151 }

uint32_t HcalTB02SD::setDetUnitId ( G4Step *  step  )  [virtual]

Implements CaloSD.

Definition at line 106 of file HcalTB02SD.cc.

References HcalTB02NumberingScheme::getUnitID(), and numberingScheme.

00106                                                 { 
00107   return (numberingScheme == 0 ? 0 : (uint32_t)(numberingScheme->getUnitID(aStep)));
00108 }

void HcalTB02SD::setNumberingScheme ( HcalTB02NumberingScheme *  scheme  ) 

Definition at line 110 of file HcalTB02SD.cc.

References numberingScheme.

Referenced by HcalTB02SD().

00110                                                                    {
00111   if (scheme != 0) {
00112     edm::LogInfo("HcalTBSim") << "HcalTB02SD: updates numbering scheme for " 
00113                               << GetName();
00114     if (numberingScheme) delete numberingScheme;
00115     numberingScheme = scheme;
00116   }
00117 }

Member Data Documentation

double HcalTB02SD::birk1 [private]

Definition at line 51 of file HcalTB02SD.h.

Referenced by getEnergyDeposit(), and HcalTB02SD().

double HcalTB02SD::birk2 [private]

Definition at line 51 of file HcalTB02SD.h.

Referenced by getEnergyDeposit(), and HcalTB02SD().

double HcalTB02SD::birk3 [private]

Definition at line 51 of file HcalTB02SD.h.

Referenced by getEnergyDeposit(), and HcalTB02SD().

std::map<G4String,double> HcalTB02SD::lengthMap [private]

Definition at line 52 of file HcalTB02SD.h.

Referenced by crystalLength(), and initMap().

HcalTB02NumberingScheme* HcalTB02SD::numberingScheme [private]

Definition at line 48 of file HcalTB02SD.h.

Referenced by setDetUnitId(), setNumberingScheme(), and ~HcalTB02SD().

bool HcalTB02SD::useBirk [private]

Definition at line 50 of file HcalTB02SD.h.

Referenced by getEnergyDeposit(), and HcalTB02SD().

bool HcalTB02SD::useWeight [private]

Definition at line 49 of file HcalTB02SD.h.

Referenced by getEnergyDeposit(), and HcalTB02SD().

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:23:57 2009 for CMSSW by  doxygen 1.5.4