#include <ECalSD.h>

Inheritance diagram for ECalSD:

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

double	getEnergyDeposit (G4Step *) override

virtual uint16_t	getLayerIDForTimeSim (const G4Step *)

virtual uint16_t	getRadiationLength (const G4Step *)

uint32_t	setDetUnitId (const G4Step *) override

void	setNumberingScheme (EcalNumberingScheme *)

	~ECalSD () override

Public Member Functions inherited from CaloSD
	CaloSD (const std::string &aSDname, const DDCompactView &cpv, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false)

void	clear () override

void	clearHits () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

void	fillHits (edm::PCaloHitContainer &, const std::string &) override

void	Initialize (G4HCofThisEvent *HCE) override

void	PrintAll () override

bool	ProcessHits (G4Step step, G4TouchableHistory tHistory) override

bool	ProcessHits (G4GFlashSpot aSpot, G4TouchableHistory ) override

	~CaloSD () override

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

Public Member Functions inherited from SensitiveDetector
void	EndOfEvent (G4HCofThisEvent *eventHC) override

const std::vector< std::string > &	getNames () const

void	Initialize (G4HCofThisEvent *eventHC) override

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

	~SensitiveDetector () override

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 ()

Protected Member Functions
uint16_t	getDepth (const G4Step *) override

int	getTrackID (const G4Track *) override

Protected Member Functions inherited from CaloSD
G4bool	checkHit ()

CaloG4Hit *	createNewHit ()

virtual bool	filterHit (CaloG4Hit *, double)

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

int	getNumberOfHits ()

double	getResponseWt (const 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 *)

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

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

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

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

void	update (const ::EndOfEvent *) override

void	updateHit (CaloG4Hit *)

Protected Member Functions inherited from SensitiveDetector
Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point) const

Local3DPoint	FinalStepPosition (const G4Step *step, coordinates) const

Local3DPoint	InitialStepPosition (const G4Step *step, coordinates) const

Local3DPoint	LocalPostStepPosition (const G4Step *step) const

Local3DPoint	LocalPreStepPosition (const G4Step *step) const

void	NaNTrap (const G4Step *step) const

void	setNames (const std::vector< std::string > &)

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...

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

double	crystalLength (const G4LogicalVolume *)

double	curve_LY (const G4Step *)

void	getBaseNumber (const G4Step *)

double	getBirkL3 (const G4Step *)

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

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

void	initMap (const 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	isEB

bool	isEE

std::vector< const G4LogicalVolume * >	noWeight

EcalNumberingScheme *	numberingScheme

double	scaleRL

double	slopeLY

bool	storeLayerTimeSim

bool	storeRL

bool	storeTrack

EcalBaseNumber	theBaseNumber

bool	useBirk

bool	useBirkL3

std::vector< const G4LogicalVolume * >	useDepth1

std::vector< const G4LogicalVolume * >	useDepth2

bool	useWeight

std::map< const G4LogicalVolume *, double >	xtalLMap

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

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 30 of file ECalSD.h.

Constructor & Destructor Documentation

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

Definition at line 41 of file ECalSD.cc.

References ageing, ageingWithSlopeLY, birk1, birk2, birk3, birkCut, birkSlope, crystalMat, depth1Name, depth2Name, ALCARECOTkAlBeamHalo_cff::filter, DDFilteredView::firstChild(), g, getDDDArray(), edm::ParameterSet::getParameter(), getStringArray(), edm::ParameterSet::getUntrackedParameter(), initMap(), edm::Service< T >::isAvailable(), isEB, isEE, gen::k, CaloSD::kmaxIon, CaloSD::kmaxNeutron, CaloSD::kmaxProton, TFileDirectory::make(), DDFilteredView::mergedSpecifics(), MeV, TFileService::mkdir(), scaleRL, EnergyResolutionVsLumi::setLumies(), setNumberingScheme(), slopeLY, storeLayerTimeSim, storeRL, storeTrack, AlCaHLTBitMon_QueryRunRegistry::string, CaloSD::suppressHeavy, compare::tfile, fftjetcommon_cfi::title, useBirk, useBirkL3, useWeight, and TrackerOfflineValidation_Dqm_cff::xmin.

                                                                           : 
   CaloSD(name, cpv, clg, p, manager, 
          (float)(p.getParameter<edm::ParameterSet>("ECalSD").getParameter<double>("TimeSliceUnit")),
          p.getParameter<edm::ParameterSet>("ECalSD").getParameter<bool>("IgnoreTrackID")), 
   numberingScheme(nullptr){
 
   //   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");
   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);
   scaleRL      = m_EC.getUntrackedParameter<double>("ScaleRadLength",1.0);
 
   //Changes for improved timing simulation
   storeLayerTimeSim = m_EC.getUntrackedParameter<bool>("StoreLayerTimeSim", 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";
   DDSpecificsMatchesValueFilter filter{DDValue(attribute,name,0)};
   DDFilteredView fv(cpv,filter);
   fv.firstChild();
   DDsvalues_type sv(fv.mergedSpecifics());
   // Use of Weight
   useWeight= true;
   std::vector<double>      tempD = getDDDArray("EnergyWeight",sv);
   if (!tempD.empty()) { if (tempD[0] < 0.1) useWeight = false; }
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("EcalSim") << "ECalSD:: useWeight " << tempD.size() << ":" 
                           << useWeight << std::endl;
 #endif
   std::vector<std::string> tempS = getStringArray("Depth1Name",sv);
   if (!tempS.empty()) depth1Name = tempS[0];
   else                  depth1Name = " ";
   tempS = getStringArray("Depth2Name",sv);
   if (!tempS.empty()) depth2Name = tempS[0];
   else                  depth2Name = " ";
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("EcalSim") << "Names (Depth 1):" << depth1Name << " (Depth 2):" 
                           << depth2Name << std::endl;
 #endif
   EcalNumberingScheme* scheme=nullptr;
   if (nullNS) {
     scheme = nullptr;
   } else if (name == "EcalHitsEB") {
     scheme = dynamic_cast<EcalNumberingScheme*>(new EcalBarrelNumberingScheme());
     isEB=true;
   } else if (name == "EcalHitsEE") { 
     scheme = dynamic_cast<EcalNumberingScheme*>(new EcalEndcapNumberingScheme());
     isEE=true;
   } 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";
   }
 
   if (scheme)  setNumberingScheme(scheme);
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("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
                           << "\tprotons below " << kmaxProton << " MeV,"
                           << "\tneutrons below " << kmaxNeutron << " MeV"
                           << "\tions below " << kmaxIon << " MeV"
                           << "\n\tDepth1 Name = " << depth1Name
                           << "\tDepth2 Name = " << depth2Name
                           << "\n\tstoreRL" << storeRL << ":" << scaleRL
                           << "\tstoreLayerTimeSim " << storeLayerTimeSim
                           << "\n\ttime Granularity " << p.getParameter<edm::ParameterSet>("ECalSD").getParameter<double>("TimeSliceUnit") << " ns"; 
   if (useWeight) initMap(name,cpv);
 #ifdef plotDebug
   edm::Service<TFileService> tfile;
   if ( tfile.isAvailable() ) {
     TFileDirectory ecDir = tfile->mkdir("ProfileFromECalSD");
     static const std::string ctype[4] = {"EB","EBref","EE","EERef"};
     for (int k=0; k<4; ++k) {
       std::string name = "ECLL_"+ctype[k];
       std::string title= "Local vs Global for "+ctype[k];
       double xmin = (k > 1) ? 3000.0 : 1000.0;
       g2L_[k] = ecDir.make<TH2F>(name.c_str(),title.c_str(),100,xmin,
                                  xmin+1000.,100,0.0,3000.);
     }
   } else {
     for (int k=0; k<4; ++k) g2L_[k] = 0;
   }
 #endif
 
 }

ECalSD::~ECalSD ( )

override

Definition at line 169 of file ECalSD.cc.

References numberingScheme.

                 {
   if (numberingScheme) delete numberingScheme;
 }

Member Function Documentation

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

private

Definition at line 523 of file ECalSD.cc.

References funct::abs(), particleFlowClusterECALTimeSelected_cfi::depth, and xtalLMap.

Referenced by curve_LY(), getLayerIDForTimeSim(), and getRadiationLength().

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

double ECalSD::crystalLength ( const G4LogicalVolume * lv )

private

Definition at line 516 of file ECalSD.cc.

References funct::abs(), and xtalLMap.

Referenced by curve_LY(), and getRadiationLength().

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

double ECalSD::curve_LY ( const G4Step * aStep )

private

Definition at line 485 of file ECalSD.cc.

References ageing, ageingWithSlopeLY, EnergyResolutionVsLumi::calcLightCollectionEfficiencyWeighted(), crystalDepth(), crystalLength(), CaloSD::currentID, particleFlowClusterECALTimeSelected_cfi::depth, CaloSD::preStepPoint, CaloSD::setToLocal(), slopeLY, CaloHitID::unitID(), and mps_merge::weight.

Referenced by getEnergyDeposit().

                                            {
 
   const G4LogicalVolume* lv = preStepPoint->GetTouchable()->GetVolume(0)->GetLogicalVolume();
 
   double weight = 1.;
   const G4ThreeVector localPoint = setToLocal(preStepPoint->GetPosition(),
                                               preStepPoint->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;
     }
   }
   return weight;
 }

void ECalSD::getBaseNumber ( const G4Step * aStep )

private

Definition at line 533 of file ECalSD.cc.

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

Referenced by setDetUnitId().

                                               {
 
   theBaseNumber.reset();
   const G4VTouchable* touch = preStepPoint->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 EDM_ML_DEBUG
       edm::LogInfo("EcalSim") << "ECalSD::getBaseNumber(): Adding level " << ii
                               << ": " << touch->GetVolume(ii)->GetName() << "["
                               << touch->GetReplicaNumber(ii) << "]";
 #endif
     }
   }
 }

double ECalSD::getBirkL3 ( const G4Step * aStep )

private

Definition at line 552 of file ECalSD.cc.

References birk1, birkCut, birkSlope, ALCARECOTkAlJpsiMuMu_cff::charge, cmsBatch::log, CaloSD::preStepPoint, and mps_merge::weight.

Referenced by getEnergyDeposit().

                                             {
 
   double weight = 1.;
   double charge = preStepPoint->GetCharge();
 
   if (charge != 0. && aStep->GetStepLength() > 0.) {
     const G4Material* mat = preStepPoint->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 EDM_ML_DEBUG
     edm::LogInfo("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 578 of file ECalSD.cc.

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

Referenced by ECalSD().

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

uint16_t ECalSD::getDepth ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 267 of file ECalSD.cc.

References any(), TauDecayModes::dec, particleFlowClusterECALTimeSelected_cfi::depth, getLayerIDForTimeSim(), getRadiationLength(), PCaloHit::kEcalDepthMask, PCaloHit::kEcalDepthOffset, PCaloHit::kEcalDepthRefz, storeLayerTimeSim, storeRL, useDepth1, useDepth2, and xtalLMap.

                                               {
   const G4StepPoint* hitPoint = aStep->GetPreStepPoint();
   const G4LogicalVolume* lv = hitPoint->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);
   } else if (storeLayerTimeSim) {
     uint16_t depth2 = getLayerIDForTimeSim(aStep);
     depth          |= ((depth2&PCaloHit::kEcalDepthMask) << PCaloHit::kEcalDepthOffset);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("EcalSim") << "ECalSD::Depth " << std::hex << depth1 << ":"
                               << depth2 << ":" << depth << std::dec << " L "
                               << (ite == xtalLMap.end()) << ":" <<ite->second;
 #endif
   return depth;
 }

double ECalSD::getEnergyDeposit ( G4Step * aStep )

overridevirtual

Reimplemented from CaloSD.

Definition at line 173 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, MeV, noWeight, CaloSD::preStepPoint, CaloSD::suppressHeavy, CaloSD::theTrack, useBirk, useBirkL3, and useWeight.

                                               {
   
   if (aStep == nullptr) {
     return 0;
   } else {
     preStepPoint = aStep->GetPreStepPoint();
     theTrack     = aStep->GetTrack();
     double wt2   = theTrack->GetWeight();
     const 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 EDM_ML_DEBUG
           if (weight == 0) 
             edm::LogInfo("EcalSim") << "Ignore Track " << theTrack->GetTrackID()
                                     << " Type " << theTrack->GetDefinition()->GetParticleName()
                                     << " Kinetic Energy " << ke << " MeV";
 #endif
         }
       }
     }
     const G4LogicalVolume* lv = preStepPoint->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) { 
       edm::LogInfo("EcalSim") << "ECalSD:: " << nameVolume
                               <<" LightWeight= " <<weight << " wt1= " <<wt1
                               << "  wt2= " << wt2 << "  "
                               << " Weighted Energy Deposit " << edep/MeV
                               << " MeV";
       const G4VProcess* pr = theTrack->GetCreatorProcess();
       if (pr) 
         edm::LogInfo("EcalSim") << theTrack->GetDefinition()->GetParticleName()
                                 << " " << theTrack->GetKineticEnergy()
                                 << " Id=" << theTrack->GetTrackID()
                                 << " IdP=" << theTrack->GetParentID()
                                 << " from  " << pr->GetProcessName();
       else
         edm::LogInfo("EcalSim") << theTrack->GetDefinition()->GetParticleName()
                                 << " " << theTrack->GetKineticEnergy()
                                 << " Id=" << theTrack->GetTrackID()
                                 << " IdP=" << theTrack->GetParentID();
     }
     */
     if(wt2 > 0.0) { edep *= wt2; }
 #ifdef EDM_ML_DEBUG
     edm::LogInfo("EcalSim") << "ECalSD:: " << nameVolume
                             << " Light Collection Efficiency " << weight << ":"
                             << wt1 << " Weighted Energy Deposit " << edep/MeV 
                             << " MeV";
 #endif
     return edep;
   } 
 }

uint16_t ECalSD::getLayerIDForTimeSim ( const G4Step * aStep )

virtual

Definition at line 331 of file ECalSD.cc.

References crystalDepth(), isEB, isEE, and CaloSD::setToLocal().

Referenced by getDepth(), and getRadiationLength().

 {
   float    layerSize = 1*cm; //layer size in cm
   if (!isEB && !isEE)  return 0;
 
   if (aStep != nullptr ) {
     const G4StepPoint* hitPoint = aStep->GetPreStepPoint();
     const G4LogicalVolume* lv = hitPoint->GetTouchable()->GetVolume(0)->GetLogicalVolume();
     const G4ThreeVector  localPoint = setToLocal(hitPoint->GetPosition(),
                                                  hitPoint->GetTouchable());
     double detz     = crystalDepth(lv,localPoint);
     if (detz<0) detz= 0;
     return (int)detz/layerSize;
   }
   return 0;
 }

uint16_t ECalSD::getRadiationLength ( const G4Step * aStep )

virtual

Definition at line 289 of file ECalSD.cc.

References funct::abs(), crystalDepth(), crystalLength(), particleFlowClusterECALTimeSelected_cfi::depth, getLayerIDForTimeSim(), GetRecoTauVFromDQM_MC_cff::kk, CaloSD::preStepPoint, rho, scaleRL, CaloSD::setToLocal(), AlCaHLTBitMon_QueryRunRegistry::string, useWeight, and z.

Referenced by getDepth().

                                                         {
   
   uint16_t thisX0 = 0;
   if (aStep != nullptr) {
     const G4StepPoint* hitPoint = aStep->GetPreStepPoint();
     const G4LogicalVolume* lv = hitPoint->GetTouchable()->GetVolume(0)->GetLogicalVolume();
 #ifdef EDM_ML_DEBUG
     edm::LogInfo("EcalSim") << lv->GetName() << " useWight " << useWeight;
 #endif
     if (useWeight) {
       const G4ThreeVector  localPoint = setToLocal(hitPoint->GetPosition(),
                                                    hitPoint->GetTouchable());
       double radl     = preStepPoint->GetMaterial()->GetRadlen();
       double depth    = crystalDepth(lv,localPoint);
       thisX0 = (uint16_t)floor(scaleRL*depth/radl);
 #ifdef plotDebug
       std::string lvname = lv->GetName();
       int k1 = (lvname.find("EFRY")!=std::string::npos) ? 2 : 0;
       int k2 = (lvname.find("refl")!=std::string::npos) ? 1 : 0;
       int kk = k1+k2;
       double rz = (k1 == 0) ? (hitPoint->GetPosition()).rho() : 
         std::abs((hitPoint->GetPosition()).z());
       edm::LogVerbatim("EcalSim") << lvname << " # " << k1 << ":" << k2 << ":" 
                                   << kk << " rz " << rz << " D " << thisX0;
       g2L_[kk]->Fill(rz,thisX0);
 #endif
 #ifdef EDM_ML_DEBUG
       double crlength = crystalLength(lv);
       edm::LogVerbatim("EcalSim") << lv->GetName() << " Global " 
                                   << hitPoint->GetPosition() << ":" 
                                   << (hitPoint->GetPosition()).rho() 
                                   << " Local " << localPoint 
                                   << " Crystal Length " << crlength 
                                   << " Radl " << radl << " DetZ " << detz 
                                   << " Index " << thisX0 
                                   << " : " << getLayerIDForTimeSim(aStep);
 #endif
     } 
   }
   return thisX0;
 }

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

private

Definition at line 597 of file ECalSD.cc.

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

Referenced by ECalSD().

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

int ECalSD::getTrackID ( const G4Track * aTrack )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 246 of file ECalSD.cc.

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

                                             {
 
   int  primaryID(0);
   bool flag(false);
   if (storeTrack) {
     const 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	(	const G4String &	sd,
		const DDCompactView &	cpv
	)

private

Definition at line 366 of file ECalSD.cc.

References any(), crystalMat, ddtrap, depth1Name, depth2Name, PVValHelper::dz, ALCARECOTkAlBeamHalo_cff::filter, DDFilteredView::firstChild(), mps_fire::i, DDFilteredView::logicalPart(), DDLogicalPart::material(), DDName::name(), dataset::name, DDBase< N, C >::name(), DDFilteredView::next(), noWeight, DDSolid::parameters(), DDSolid::shape(), DDLogicalPart::solid(), AlCaHLTBitMon_QueryRunRegistry::string, useDepth1, useDepth2, and xtalLMap.

Referenced by ECalSD().

                                                                   {
 
   G4String attribute = "ReadOutName";
   DDSpecificsMatchesValueFilter filter{DDValue(attribute,sd,0)};
   DDFilteredView fv(cpv,filter);
   fv.firstChild();
 
   std::vector<const G4LogicalVolume*> lvused;
   const G4LogicalVolumeStore *  lvs = G4LogicalVolumeStore::GetInstance();
   std::vector<const G4LogicalVolume *>::const_iterator lvcite;
   std::map<const std::string, const 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();
     const 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 EDM_ML_DEBUG
           edm::LogInfo("EcalSim") << "ECalSD::initMap Logical Volume " 
                                   << lvname <<" in Depth 1 volume list";
 #endif
         }
         const G4LogicalVolume* lvr = nameMap[lvname + "_refl"];
         if (lvr != nullptr && !any(useDepth1, lvr)) {
           useDepth1.push_back(lvr);
 #ifdef EDM_ML_DEBUG
           edm::LogInfo("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 EDM_ML_DEBUG
           edm::LogInfo("EcalSim") << "ECalSD::initMap Logical Volume " 
                                   << lvname <<" in Depth 2 volume list";
 #endif
         }
         const G4LogicalVolume* lvr = nameMap[lvname + "_refl"];
         if (lvr != nullptr && !any(useDepth2,lvr)) {
           useDepth2.push_back(lvr);
 #ifdef EDM_ML_DEBUG
           edm::LogInfo("EcalSim") << "ECalSD::initMap Logical Volume " 
                                   << lvname << "_refl"
                                   <<" in Depth 2 volume list";
 #endif
         }
       }
     }
     if (lv != nullptr) {
       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 EDM_ML_DEBUG
           edm::LogInfo("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<const G4LogicalVolume*,double>(lv,dz*type));
             lv = nameMap[lvname + "_refl"];
             if (lv != nullptr) {
               xtalLMap.insert(std::pair<const G4LogicalVolume*,double>(lv,-dz*type));
             }
           }
         }
       } else {
         if (!any(noWeight,lv)) {
           noWeight.push_back(lv);
 #ifdef EDM_ML_DEBUG
           edm::LogInfo("EcalSim") << "ECalSD::initMap Logical Volume "
                                   << lvname << " Material " << matname 
                                   << " in noWeight list";
 #endif
         }
         lv = nameMap[lvname];
         if (lv != nullptr && !any(noWeight,lv)) {
           noWeight.push_back(lv);
 #ifdef EDM_ML_DEBUG
           edm::LogInfo("EcalSim") << "ECalSD::initMap Logical Volume " 
                                   << lvname << " Material " << matname 
                                   << " in noWeight list";
 #endif
         }
       }
     }
     dodet = fv.next();
   }
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("EcalSim") << "ECalSD: Length Table for " << attribute << " = " 
                           << sd << ":";   
   int i=0;
   for (auto ite : xtalLMap) {
     G4String name("Unknown");
     if (ite.first != 0) name = (ite.first)->GetName();
     edm::LogInfo("EcalSim") << " " << i << " " << ite.first << " " << name 
                             << " L = " << ite.second;
     ++i;
   }
 #endif
 }

uint32_t ECalSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements CaloSD.

Definition at line 348 of file ECalSD.cc.

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

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

void ECalSD::setNumberingScheme ( EcalNumberingScheme * scheme )

Definition at line 357 of file ECalSD.cc.

References numberingScheme.

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

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

Member Data Documentation

EnergyResolutionVsLumi ECalSD::ageing

private

Definition at line 72 of file ECalSD.h.

Referenced by curve_LY(), and ECalSD().

bool ECalSD::ageingWithSlopeLY

private

Definition at line 73 of file ECalSD.h.

Referenced by curve_LY(), and ECalSD().

double ECalSD::birk1

private