#include <ECalSD.h>

Inheritance diagram for ECalSD:

Public Member Functions
	ECalSD (G4String, const DDCompactView &, SensitiveDetectorCatalog &, edm::ParameterSet const &p, const SimTrackManager *)

virtual uint16_t	getDepth (G4Step *)

virtual double	getEnergyDeposit (G4Step *)

virtual uint16_t	getRadiationLength (G4Step *)

virtual int	getTrackID (G4Track *)

virtual uint32_t	setDetUnitId (G4Step *)

void	setNumberingScheme (EcalNumberingScheme *)

virtual	~ECalSD ()

Public Member Functions inherited from CaloSD
	CaloSD (G4String aSDname, const DDCompactView &cpv, SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, int tSlice=1, bool ignoreTkID=false)

virtual void	clear ()

virtual void	DrawAll ()

virtual void	EndOfEvent (G4HCofThisEvent *eventHC)

void	fillHits (edm::PCaloHitContainer &, std::string n)

virtual void	Initialize (G4HCofThisEvent *HCE)

virtual void	PrintAll ()

virtual bool	ProcessHits (G4Step step, G4TouchableHistory tHistory)

virtual bool	ProcessHits (G4GFlashSpot aSpot, G4TouchableHistory )

virtual	~CaloSD ()

Public Member Functions inherited from SensitiveCaloDetector
	SensitiveCaloDetector (std::string &iname, const DDCompactView &cpv, SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)

Public Member Functions inherited from SensitiveDetector
virtual void	AssignSD (std::string &vname)

Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point)

Local3DPoint	FinalStepPosition (G4Step *s, coordinates)

virtual std::vector< std::string >	getNames ()

Local3DPoint	InitialStepPosition (G4Step *s, coordinates)

std::string	nameOfSD ()

void	NaNTrap (G4Step *step)

void	Register ()

	SensitiveDetector (std::string &iname, const DDCompactView &cpv, SensitiveDetectorCatalog &, edm::ParameterSet const &p)

virtual	~SensitiveDetector ()

Public Member Functions inherited from Observer< const BeginOfRun * >
	Observer ()

void	slotForUpdate (const BeginOfRun *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfEvent * >
	Observer ()

void	slotForUpdate (const BeginOfEvent *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfTrack * >
	Observer ()

void	slotForUpdate (const BeginOfTrack *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfTrack * >
	Observer ()

void	slotForUpdate (const EndOfTrack *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfEvent * >
	Observer ()

void	slotForUpdate (const EndOfEvent *iT)

virtual	~Observer ()

Private Member Functions
double	crystalDepth (G4LogicalVolume *, const G4ThreeVector &)

double	crystalLength (G4LogicalVolume *)

double	curve_LY (G4Step *)

void	getBaseNumber (const G4Step *)

double	getBirkL3 (G4Step *)

std::vector< double >	getDDDArray (const std::string &, const DDsvalues_type &)

std::vector< std::string >	getStringArray (const std::string &, const DDsvalues_type &)

void	initMap (G4String, const DDCompactView &)

Private Attributes
EnergyResolutionVsLumi	ageing

bool	ageingWithSlopeLY

double	birk1

double	birk2

double	birk3

double	birkCut

double	birkSlope

std::string	crystalMat

std::string	depth1Name

std::string	depth2Name

bool	ignoreDepCorr

std::vector< G4LogicalVolume * >	noWeight

EcalNumberingScheme *	numberingScheme

double	slopeLY

bool	storeRL

bool	storeTrack

EcalBaseNumber	theBaseNumber

bool	useBirk

bool	useBirkL3

std::vector< G4LogicalVolume * >	useDepth1

std::vector< G4LogicalVolume * >	useDepth2

bool	useWeight

std::map< G4LogicalVolume *, double >	xtalLMap

Additional Inherited Members
Public Types inherited from SensitiveDetector
enum	coordinates { WorldCoordinates, LocalCoordinates }

Protected Member Functions inherited from CaloSD
G4bool	checkHit ()

virtual void	clearHits ()

CaloG4Hit *	createNewHit ()

virtual bool	filterHit (CaloG4Hit *, double)

double	getAttenuation (G4Step *aStep, double birk1, double birk2, double birk3)

int	getNumberOfHits ()

double	getResponseWt (G4Track *)

virtual G4bool	getStepInfo (G4Step *aStep)

G4bool	hitExists ()

virtual void	initRun ()

void	resetForNewPrimary (const G4ThreeVector &, double)

G4ThreeVector	setToGlobal (const G4ThreeVector &, const G4VTouchable *)

G4ThreeVector	setToLocal (const G4ThreeVector &, const G4VTouchable *)

virtual void	update (const BeginOfRun *)
	This routine will be called when the appropriate signal arrives. More...

virtual void	update (const BeginOfEvent *)
	This routine will be called when the appropriate signal arrives. More...

virtual void	update (const BeginOfTrack *trk)
	This routine will be called when the appropriate signal arrives. More...

virtual void	update (const EndOfTrack *trk)
	This routine will be called when the appropriate signal arrives. More...

virtual void	update (const ::EndOfEvent *)

void	updateHit (CaloG4Hit *)

Protected Member Functions inherited from Observer< const EndOfEvent * >
virtual void	update (const EndOfEvent *)=0
	This routine will be called when the appropriate signal arrives. More...

Protected Attributes inherited from CaloSD
int	checkHits

double	correctT

bool	corrTOFBeam

CaloG4Hit *	currentHit

CaloHitID	currentID

float	edepositEM

float	edepositHAD

double	eminHit

double	eminHitD

G4int	emPDG

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

G4int	epPDG

bool	forceSave

G4int	gammaPDG

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

const SimTrackManager *	m_trackManager

G4ThreeVector	posGlobal

G4StepPoint *	preStepPoint

CaloHitID	previousID

int	primIDSaved

bool	runInit

bool	suppressHeavy

G4Track *	theTrack

double	tmaxHit

bool	useMap

Detailed Description

Definition at line 26 of file ECalSD.h.

Constructor & Destructor Documentation

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

Definition at line 39 of file ECalSD.cc.

References DDFilteredView::addFilter(), ageing, ageingWithSlopeLY, birk1, birk2, birk3, birkCut, birkSlope, crystalMat, depth1Name, depth2Name, DDSpecificsFilter::equals, alcazmumu_cfi::filter, DDFilteredView::firstChild(), g, getDDDArray(), edm::ParameterSet::getParameter(), getStringArray(), edm::ParameterSet::getUntrackedParameter(), ignoreDepCorr, initMap(), CaloSD::kmaxIon, CaloSD::kmaxNeutron, CaloSD::kmaxProton, LogDebug, DDFilteredView::mergedSpecifics(), MeV, DDSpecificsFilter::setCriteria(), EnergyResolutionVsLumi::setLumies(), setNumberingScheme(), slopeLY, storeRL, storeTrack, AlCaHLTBitMon_QueryRunRegistry::string, CaloSD::suppressHeavy, useBirk, useBirkL3, and useWeight.

                                                                       : 
   CaloSD(name, cpv, clg, p, manager, 
      p.getParameter<edm::ParameterSet>("ECalSD").getParameter<int>("TimeSliceUnit"),
      p.getParameter<edm::ParameterSet>("ECalSD").getParameter<bool>("IgnoreTrackID")), 
   numberingScheme(0){
 
   //   static SimpleConfigurable<bool>   on1(false,  "ECalSD:UseBirkLaw");
   //   static SimpleConfigurable<double> bk1(0.00463,"ECalSD:BirkC1");
   //   static SimpleConfigurable<double> bk2(-0.03,  "ECalSD:BirkC2");
   //   static SimpleConfigurable<double> bk3(1.0,    "ECalSD:BirkC3");
   // Values from NIM A484 (2002) 239-244: as implemented in Geant3
   //   useBirk          = on1.value();
   //   birk1            = bk1.value()*(g/(MeV*cm2));
   //   birk2            = bk2.value()*(g/(MeV*cm2))*(g/(MeV*cm2));
 
   edm::ParameterSet m_EC = p.getParameter<edm::ParameterSet>("ECalSD");
   useBirk      = m_EC.getParameter<bool>("UseBirkLaw");
   useBirkL3    = m_EC.getParameter<bool>("BirkL3Parametrization");
   birk1        = m_EC.getParameter<double>("BirkC1")*(g/(MeV*cm2));
   birk2        = m_EC.getParameter<double>("BirkC2");
   birk3        = m_EC.getParameter<double>("BirkC3");
   birkSlope    = m_EC.getParameter<double>("BirkSlope");
   birkCut      = m_EC.getParameter<double>("BirkCut");
   slopeLY      = m_EC.getParameter<double>("SlopeLightYield");
   storeTrack   = m_EC.getParameter<bool>("StoreSecondary");
   ignoreDepCorr= m_EC.getParameter<bool>("IgnoreDepthCorr");
   crystalMat   = m_EC.getUntrackedParameter<std::string>("XtalMat","E_PbWO4");
   bool isItTB  = m_EC.getUntrackedParameter<bool>("TestBeam", false);
   bool nullNS  = m_EC.getUntrackedParameter<bool>("NullNumbering", false);
   storeRL      = m_EC.getUntrackedParameter<bool>("StoreRadLength", false);
   
   ageingWithSlopeLY   = m_EC.getUntrackedParameter<bool>("AgeingWithSlopeLY", false);
   if(ageingWithSlopeLY) ageing.setLumies(p.getParameter<edm::ParameterSet>("ECalSD").getParameter<double>("DelivLuminosity"),
                      p.getParameter<edm::ParameterSet>("ECalSD").getParameter<double>("InstLuminosity"));
 
   //Material list for HB/HE/HO sensitive detectors
   std::string attribute = "ReadOutName";
   DDSpecificsFilter filter;
   DDValue           ddv(attribute,name,0);
   filter.setCriteria(ddv,DDSpecificsFilter::equals);
   DDFilteredView fv(cpv);
   fv.addFilter(filter);
   fv.firstChild();
   DDsvalues_type sv(fv.mergedSpecifics());
   // Use of Weight
   useWeight= true;
   std::vector<double>      tempD = getDDDArray("EnergyWeight",sv);
   if (tempD.size() > 0) { if (tempD[0] < 0.1) useWeight = false; }
   std::vector<std::string> tempS = getStringArray("Depth1Name",sv);
   if (tempS.size() > 0) depth1Name = tempS[0];
   else                  depth1Name = " ";
   tempS = getStringArray("Depth2Name",sv);
   if (tempS.size() > 0) depth2Name = tempS[0];
   else                  depth2Name = " ";
 
   EcalNumberingScheme* scheme=0;
   if (nullNS)                    scheme = 0;
   else if (name == "EcalHitsEB") scheme = dynamic_cast<EcalNumberingScheme*>(new EcalBarrelNumberingScheme());
   else if (name == "EcalHitsEE") scheme = dynamic_cast<EcalNumberingScheme*>(new EcalEndcapNumberingScheme());
   else if (name == "EcalHitsES") {
     if (isItTB) scheme = dynamic_cast<EcalNumberingScheme*>(new ESTBNumberingScheme());
     else        scheme = dynamic_cast<EcalNumberingScheme*>(new EcalPreshowerNumberingScheme());
     useWeight = false;
   } else {edm::LogWarning("EcalSim") << "ECalSD: ReadoutName not supported\n";}
 
   if (scheme)  setNumberingScheme(scheme);
 #ifdef DebugLog
   LogDebug("EcalSim") << "Constructing a ECalSD  with name " << GetName();
 #endif
   if (useWeight) {
     edm::LogInfo("EcalSim")  << "ECalSD:: Use of Birks law is set to      " 
                  << useBirk << "        with three constants kB = "
                  << birk1 << ", C1 = " << birk2 << ", C2 = " 
                  << birk3 <<"\n         Use of L3 parametrization "
                  << useBirkL3 << " with slope " << birkSlope
                  << " and cut off " << birkCut << "\n"
                  << "         Slope for Light yield is set to "
                  << slopeLY;
   } else {
     edm::LogInfo("EcalSim")  << "ECalSD:: energy deposit is not corrected "
                  << " by Birk or light yield curve";
   }
 
   edm::LogInfo("EcalSim") << "ECalSD:: Suppression Flag " << suppressHeavy
               << " protons below " << kmaxProton << " MeV,"
               << " neutrons below " << kmaxNeutron << " MeV and"
               << " ions below " << kmaxIon << " MeV\n"
               << "         Depth1 Name = " << depth1Name
               << " and Depth2 Name = " << depth2Name;
 
   if (useWeight) initMap(name,cpv);
 
 }

ECalSD::~ECalSD ( )

virtual

Definition at line 135 of file ECalSD.cc.

References numberingScheme.

                 {
   if (numberingScheme) delete numberingScheme;
 }

Member Function Documentation

double ECalSD::crystalDepth	(	G4LogicalVolume *	lv,
		const G4ThreeVector &	localPoint
	)

private

Definition at line 440 of file ECalSD.cc.

References funct::abs(), ignoreDepCorr, and xtalLMap.

Referenced by curve_LY(), and getRadiationLength().

                                                              {
 
   auto ite = xtalLMap.find(lv);
   double depth(0);
   if (ite != xtalLMap.end()) {
     if (ignoreDepCorr) depth = (0.5*std::abs(ite->second)+localPoint.z());
     else               depth = std::abs(0.5*(ite->second)+localPoint.z());
   }
   return depth;
 }

double ECalSD::crystalLength ( G4LogicalVolume * lv )

private

Definition at line 434 of file ECalSD.cc.

References funct::abs(), and xtalLMap.

Referenced by curve_LY().

                                                 {
 
   auto ite = xtalLMap.find(lv);
   double length = (ite == xtalLMap.end()) ? 230.0 : std::abs(ite->second);
   return length;
 }

double ECalSD::curve_LY ( G4Step * aStep )

private

Definition at line 394 of file ECalSD.cc.

References ageing, ageingWithSlopeLY, EnergyResolutionVsLumi::calcLightCollectionEfficiencyWeighted(), crystalDepth(), crystalLength(), CaloSD::currentID, LogDebug, CaloSD::setToLocal(), slopeLY, CaloHitID::unitID(), and histoStyle::weight.

Referenced by getEnergyDeposit().

                                      {
 
   G4StepPoint*     stepPoint = aStep->GetPreStepPoint();
   G4LogicalVolume* lv        = stepPoint->GetTouchable()->GetVolume(0)->GetLogicalVolume();
 
   double weight = 1.;
   G4ThreeVector  localPoint = setToLocal(stepPoint->GetPosition(),
                      stepPoint->GetTouchable());
 
   double crlength = crystalLength(lv);
   double depth    = crystalDepth(lv,localPoint);
 
   if(ageingWithSlopeLY){
     //position along the crystal in mm from 0 to 230 (in EB)
     if (depth >= -0.1 || depth <= crlength+0.1)
       weight = ageing.calcLightCollectionEfficiencyWeighted(currentID.unitID(), depth/crlength);
   }
   else{
     double dapd = crlength - depth;
     if (dapd >= -0.1 || dapd <= crlength+0.1) {
       if (dapd <= 100.)
     weight = 1.0 + slopeLY - dapd * 0.01 * slopeLY;
     } else {
       edm::LogWarning("EcalSim") << "ECalSD: light coll curve : wrong distance "
                  << "to APD " << dapd << " crlength = " 
                  << crlength <<" crystal name = " <<lv->GetName()
                  << " z of localPoint = " << localPoint.z() 
                  << " take weight = " << weight;
     }
   }
 #ifdef DebugLog
   LogDebug("EcalSim") << "ECalSD, light coll curve : " << dapd 
               << " crlength = " << crlength
               << " crystal name = " << lv->GetName()
               << " z of localPoint = " << localPoint.z() 
               << " take weight = " << weight;
 #endif
   return weight;
 }

void ECalSD::getBaseNumber ( const G4Step * aStep )

private

Definition at line 452 of file ECalSD.cc.

References EcalBaseNumber::addLevel(), EcalBaseNumber::getCapacity(), cuy::ii, LogDebug, EcalBaseNumber::reset(), EcalBaseNumber::setSize(), and theBaseNumber.

Referenced by setDetUnitId().

                                               {
 
   theBaseNumber.reset();
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   int theSize = touch->GetHistoryDepth()+1;
   if ( theBaseNumber.getCapacity() < theSize ) theBaseNumber.setSize(theSize);
   //Get name and copy numbers
   if ( theSize > 1 ) {
     for (int ii = 0; ii < theSize ; ii++) {
       theBaseNumber.addLevel(touch->GetVolume(ii)->GetName(),touch->GetReplicaNumber(ii));
 #ifdef DebugLog
       LogDebug("EcalSim") << "ECalSD::getBaseNumber(): Adding level " << ii
                           << ": " << touch->GetVolume(ii)->GetName() << "["
                           << touch->GetReplicaNumber(ii) << "]";
 #endif
     }
   }
 
 }

double ECalSD::getBirkL3 ( G4Step * aStep )

private

Definition at line 472 of file ECalSD.cc.

References birk1, birkCut, birkSlope, DeDxDiscriminatorTools::charge(), log, LogDebug, and histoStyle::weight.

Referenced by getEnergyDeposit().

                                       {
 
   double weight = 1.;
   double charge = aStep->GetPreStepPoint()->GetCharge();
 
   if (charge != 0. && aStep->GetStepLength() > 0) {
     G4Material* mat = aStep->GetPreStepPoint()->GetMaterial();
     double density = mat->GetDensity();
     double dedx    = aStep->GetTotalEnergyDeposit()/aStep->GetStepLength();
     double rkb     = birk1/density;
     if (dedx > 0) {
       weight         = 1. - birkSlope*log(rkb*dedx);
       if (weight < birkCut) weight = birkCut;
       else if (weight > 1.) weight = 1.;
     }
 #ifdef DebugLog
     LogDebug("EcalSim") << "ECalSD::getBirkL3 in " << mat->GetName()
                         << " Charge " << charge << " dE/dx " << dedx
                         << " Birk Const " << rkb << " Weight = " << weight 
             << " dE " << aStep->GetTotalEnergyDeposit();
 #endif
   }
   return weight;
 
 }

std::vector< double > ECalSD::getDDDArray	(	const std::string &	str,
		const DDsvalues_type &	sv
	)

private

Definition at line 498 of file ECalSD.cc.

References DDfetch(), DDValue::doubles(), LogDebug, and relativeConstraints::value.

Referenced by ECalSD().

                                                {
 
 #ifdef DebugLog
   LogDebug("EcalSim") << "ECalSD:getDDDArray called for " << str;
 #endif
   DDValue value(str);
   if (DDfetch(&sv,value)) {
 #ifdef DebugLog
     LogDebug("EcalSim") << value;
 #endif
     const std::vector<double> & fvec = value.doubles();
     return fvec;
   } else {
     std::vector<double> fvec;
     return fvec;
   }
 }

uint16_t ECalSD::getDepth ( G4Step * aStep )

virtual

Reimplemented from CaloSD.

Definition at line 227 of file ECalSD.cc.

References any(), getRadiationLength(), PCaloHit::kEcalDepthMask, PCaloHit::kEcalDepthOffset, PCaloHit::kEcalDepthRefz, storeRL, useDepth1, useDepth2, and xtalLMap.

                                         {
   G4LogicalVolume* lv   = aStep->GetPreStepPoint()->GetTouchable()->GetVolume(0)->GetLogicalVolume();
   uint16_t depth  = any(useDepth1,lv) ? 1 : (any(useDepth2,lv) ? 2 : 0);
   if (storeRL) {
     auto ite        = xtalLMap.find(lv);
     uint16_t depth1 = (ite == xtalLMap.end()) ? 0 : (((ite->second) >= 0) ? 0 :
                                                      PCaloHit::kEcalDepthRefz);
     uint16_t depth2 = getRadiationLength(aStep);
     depth          |= (((depth2&PCaloHit::kEcalDepthMask) << PCaloHit::kEcalDepthOffset) | depth1);
   }
   return depth;
 }

double ECalSD::getEnergyDeposit ( G4Step * aStep )

virtual

Reimplemented from CaloSD.

Definition at line 139 of file ECalSD.cc.

References any(), birk1, birk2, birk3, curve_LY(), CaloSD::getAttenuation(), getBirkL3(), CaloSD::getResponseWt(), TrackInformation::isPrimary(), ke, CaloSD::kmaxIon, CaloSD::kmaxNeutron, CaloSD::kmaxProton, LogDebug, MeV, noWeight, NULL, CaloSD::preStepPoint, CaloSD::suppressHeavy, CaloSD::theTrack, useBirk, useBirkL3, useWeight, and histoStyle::weight.

                                               {
   
   if (aStep == NULL) {
     return 0;
   } else {
     preStepPoint        = aStep->GetPreStepPoint();
     G4Track* theTrack   = aStep->GetTrack();
     double wt2  = theTrack->GetWeight();
     G4String nameVolume = preStepPoint->GetPhysicalVolume()->GetName();
 
     // take into account light collection curve for crystals
     double weight = 1.;
     if (suppressHeavy) {
       TrackInformation * trkInfo = (TrackInformation *)(theTrack->GetUserInformation());
       if (trkInfo) {
     int pdg = theTrack->GetDefinition()->GetPDGEncoding();
     if (!(trkInfo->isPrimary())) { // Only secondary particles
       double ke = theTrack->GetKineticEnergy()/MeV;
       if (((pdg/1000000000 == 1 && ((pdg/10000)%100) > 0 &&
         ((pdg/10)%100) > 0)) && (ke<kmaxIon)) weight = 0;
       if ((pdg == 2212) && (ke < kmaxProton))     weight = 0;
       if ((pdg == 2112) && (ke < kmaxNeutron))    weight = 0;
 #ifdef DebugLog
       if (weight == 0) 
         LogDebug("EcalSim") << "Ignore Track " << theTrack->GetTrackID()
                 << " Type " << theTrack->GetDefinition()->GetParticleName()
                 << " Kinetic Energy " << ke << " MeV";
 #endif
     }
       }
     }
     G4LogicalVolume* lv   = aStep->GetPreStepPoint()->GetTouchable()->GetVolume(0)->GetLogicalVolume();
     if (useWeight && !any(noWeight,lv)) {
       weight *= curve_LY(aStep);
       if (useBirk) {
     if (useBirkL3) weight *= getBirkL3(aStep);
     else           weight *= getAttenuation(aStep, birk1, birk2, birk3);
       }
     }
     double wt1  = getResponseWt(theTrack);
     double edep = aStep->GetTotalEnergyDeposit()*weight*wt1;
     /*
     if(wt2 != 1.0) { 
       std::cout << "ECalSD:: " << nameVolume
         <<" LightWeight= " <<weight << " wt1= " <<wt1
         << "  wt2= " << wt2 << "  "
         << " Weighted Energy Deposit " << edep/MeV << " MeV" 
         << std::endl;
       std::cout << theTrack->GetDefinition()->GetParticleName()
         << " " << theTrack->GetKineticEnergy()
         << " Id=" << theTrack->GetTrackID()
         << " IdP=" << theTrack->GetParentID();
       const G4VProcess* pr = theTrack->GetCreatorProcess();
       if(pr) std::cout << " from  " << pr->GetProcessName();
       std::cout << std::endl;
     }
     */
     if(wt2 > 0.0) { edep *= wt2; }
 #ifdef DebugLog
     LogDebug("EcalSim") << "ECalSD:: " << nameVolume
             <<" Light Collection Efficiency " <<weight << ":" <<wt1
             << " Weighted Energy Deposit " << edep/MeV << " MeV";
 #endif
     return edep;
   } 
 }

uint16_t ECalSD::getRadiationLength ( G4Step * aStep )

virtual

Definition at line 240 of file ECalSD.cc.

References crystalDepth(), NULL, CaloSD::setToLocal(), and useWeight.

Referenced by getDepth().

                                                   {
   
   uint16_t thisX0 = 0;
   if (aStep != NULL) {
     G4StepPoint* hitPoint = aStep->GetPreStepPoint();
     G4LogicalVolume* lv   = hitPoint->GetTouchable()->GetVolume(0)->GetLogicalVolume();
     
     if (useWeight) {
       G4ThreeVector  localPoint = setToLocal(hitPoint->GetPosition(),
                          hitPoint->GetTouchable());
       double radl     = hitPoint->GetMaterial()->GetRadlen();
       double depth    = crystalDepth(lv,localPoint);
       thisX0 = (uint16_t)floor(depth/radl);
     } 
   }
   return thisX0;
 }

std::vector< std::string > ECalSD::getStringArray	(	const std::string &	str,
		const DDsvalues_type &	sv
	)

private

Definition at line 517 of file ECalSD.cc.

References DDfetch(), LogDebug, DDValue::strings(), and relativeConstraints::value.

Referenced by ECalSD().

                                                    {
 
 #ifdef DebugLog
   LogDebug("EcalSim") << "ECalSD:getStringArray called for " << str;
 #endif
   DDValue value(str);
   if (DDfetch(&sv,value)) {
 #ifdef DebugLog
     LogDebug("EcalSim") << value;
 #endif
     const std::vector<std::string> & fvec = value.strings();
     return fvec;
   } else {
     std::vector<std::string> fvec;
     return fvec;
   }
 }

int ECalSD::getTrackID ( G4Track * aTrack )

virtual

Reimplemented from CaloSD.

Definition at line 206 of file ECalSD.cc.

References any(), archive::flag, CaloSD::forceSave, CaloSD::getTrackID(), CaloSD::preStepPoint, storeTrack, useDepth1, and useDepth2.

                                       {
 
   int  primaryID(0);
   bool flag(false);
   if (storeTrack) {
     G4LogicalVolume* lv  = preStepPoint->GetTouchable()->GetVolume(0)->GetLogicalVolume();
     if (any(useDepth1,lv)) {
       flag = true;
     } else if (any(useDepth2,lv)) {
       flag = true;
     }
   }
   if (flag) {
     forceSave = true;
     primaryID = aTrack->GetTrackID();
   } else {
     primaryID = CaloSD::getTrackID(aTrack);
   }
   return primaryID;
 }

void ECalSD::initMap	(	G4String	sd,
		const DDCompactView &	cpv
	)

private

Definition at line 277 of file ECalSD.cc.

References DDFilteredView::addFilter(), any(), crystalMat, ddtrap, depth1Name, depth2Name, DDSpecificsFilter::equals, alcazmumu_cfi::filter, DDFilteredView::firstChild(), i, LogDebug, DDFilteredView::logicalPart(), DDLogicalPart::material(), DDName::name(), SensitiveDetector::name, DDBase< N, C >::name(), DDFilteredView::next(), noWeight, DDSolid::parameters(), DDSpecificsFilter::setCriteria(), DDSolid::shape(), DDLogicalPart::solid(), AlCaHLTBitMon_QueryRunRegistry::string, useDepth1, useDepth2, and xtalLMap.

Referenced by ECalSD().

                                                            {
 
   G4String attribute = "ReadOutName";
   DDSpecificsFilter filter;
   DDValue           ddv(attribute,sd,0);
   filter.setCriteria(ddv,DDSpecificsFilter::equals);
   DDFilteredView fv(cpv);
   fv.addFilter(filter);
   fv.firstChild();
 
   std::vector<G4LogicalVolume*> lvused;
   const G4LogicalVolumeStore *  lvs = G4LogicalVolumeStore::GetInstance();
   std::vector<G4LogicalVolume *>::const_iterator lvcite;
   std::map<std::string, G4LogicalVolume *> nameMap;
   for (auto lvi = lvs->begin(), lve = lvs->end(); lvi != lve; ++lvi)
     nameMap.insert(std::make_pair((*lvi)->GetName(), *lvi));
 
   bool dodet=true;
   while (dodet) {
     const std::string &matname = fv.logicalPart().material().name().name();
     const std::string &lvname = fv.logicalPart().name().name();
     G4LogicalVolume* lv = nameMap[lvname];
     int ibec = (lvname.find("EFRY") == std::string::npos) ? 0 : 1;
     int iref = (lvname.find("refl") == std::string::npos) ? 0 : 1;
     int type = (ibec+iref == 1) ? 1 : -1;
     if (depth1Name != " ") {
       if (strncmp(lvname.c_str(), depth1Name.c_str(), 4) == 0) {
     if (!any(useDepth1, lv)) {
       useDepth1.push_back(lv);
 #ifdef DebugLog
       LogDebug("EcalSim") << "ECalSD::initMap Logical Volume " << lvname
                   <<" in Depth 1 volume list";
 #endif
     }
     G4LogicalVolume* lvr = nameMap[lvname + "_refl"];
     if (lvr != 0 && !any(useDepth1, lvr)) {
       useDepth1.push_back(lvr);
 #ifdef DebugLog
       LogDebug("EcalSim") << "ECalSD::initMap Logical Volume " << lvname << "_refl"
                   <<" in Depth 1 volume list";
 #endif
     }
       }
     }
     if (depth2Name != " ") {
       if (strncmp(lvname.c_str(), depth2Name.c_str(), 4) == 0) {
     if (!any(useDepth2, lv)) {
       useDepth2.push_back(lv);
 #ifdef DebugLog
       LogDebug("EcalSim") << "ECalSD::initMap Logical Volume " << lvname
                   <<" in Depth 2 volume list";
 #endif
     }
     G4LogicalVolume* lvr = nameMap[lvname + "_refl"];
     if (lvr != 0 && !any(useDepth2,lvr)) {
       useDepth2.push_back(lvr);
 #ifdef DebugLog
       LogDebug("EcalSim") << "ECalSD::initMap Logical Volume " << lvname << "_refl"
                   <<" in Depth 2 volume list";
 #endif
     }
       }
     }
     if (lv != 0) {
       if (crystalMat.size() == matname.size() && !strcmp(crystalMat.c_str(), matname.c_str())) {
     if (!any(lvused,lv)) {
       lvused.push_back(lv);
       const DDSolid & sol  = fv.logicalPart().solid();
       const std::vector<double> & paras = sol.parameters();
 #ifdef DebugLog
       LogDebug("EcalSim") << "ECalSD::initMap (for " << sd << "): Solid " 
                   << lvname << " Shape " << sol.shape() 
                   << " Parameter 0 = "<< paras[0] 
                   << " Logical Volume " << lv;
 #endif
       if (sol.shape() == ddtrap) {
         double dz = 2*paras[0];
         xtalLMap.insert(std::pair<G4LogicalVolume*,double>(lv,dz*type));
         lv = nameMap[lvname + "_refl"];
         if (lv != 0)
           xtalLMap.insert(std::pair<G4LogicalVolume*,double>(lv,-dz*type));
       }
     }
       } else {
     if (!any(noWeight,lv)) {
       noWeight.push_back(lv);
 #ifdef DebugLog
       LogDebug("EcalSim") << "ECalSD::initMap Logical Volume " << lvname
                   << " Material " << matname <<" in noWeight list";
 #endif
     }
     lv = nameMap[lvname];
     if (lv != 0 && !any(noWeight,lv)) {
       noWeight.push_back(lv);
 #ifdef DebugLog
       LogDebug("EcalSim") << "ECalSD::initMap Logical Volume " << lvname
                   << " Material " << matname <<" in noWeight list";
 #endif
     }
       }
     }
     dodet = fv.next();
   }
 #ifdef DebugLog
   LogDebug("EcalSim") << "ECalSD: Length Table for " << attribute << " = " 
               << sd << ":";   
   std::map<G4LogicalVolume*,double>::const_iterator ite = xtalLMap.begin();
   int i=0;
   for (; ite != xtalLMap.end(); ite++, i++) {
     G4String name = "Unknown";
     if (ite->first != 0) name = (ite->first)->GetName();
     LogDebug("EcalSim") << " " << i << " " << ite->first << " " << name 
             << " L = " << ite->second;
   }
 #endif
 }

uint32_t ECalSD::setDetUnitId ( G4Step * aStep )

virtual

Implements CaloSD.

Definition at line 258 of file ECalSD.cc.

References getBaseNumber(), EcalNumberingScheme::getUnitID(), numberingScheme, and theBaseNumber.

                                             { 
   if (numberingScheme == 0) {
     return EBDetId(1,1)();
   } else {
     getBaseNumber(aStep);
     return numberingScheme->getUnitID(theBaseNumber);
   }
 }

void ECalSD::setNumberingScheme ( EcalNumberingScheme * scheme )

Definition at line 267 of file ECalSD.cc.

References numberingScheme.

Referenced by ECalSD(), and HcalTB04Analysis::update().

                                                            {
   if (scheme != 0) {
     edm::LogInfo("EcalSim") << "EcalSD: updates numbering scheme for " 
                 << GetName() << "\n";
     if (numberingScheme) delete numberingScheme;
     numberingScheme = scheme;
   }
 }

Member Data Documentation

EnergyResolutionVsLumi ECalSD::ageing

private

Definition at line 62 of file ECalSD.h.

Referenced by curve_LY(), and ECalSD().

bool ECalSD::ageingWithSlopeLY

private

Definition at line 63 of file ECalSD.h.

Referenced by curve_LY(), and ECalSD().

double ECalSD::birk1

private