#include <SimG4CMS/Forward/interface/CastorSD.h>

Inheritance diagram for CastorSD:

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

uint32_t	setDetUnitId (const G4Step *step) override

void	setNumberingScheme (CastorNumberingScheme *scheme)

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

G4bool	ProcessHits (G4Step step, G4TouchableHistory ) override

bool	ProcessHits (G4GFlashSpot aSpot, G4TouchableHistory ) override

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

bool	isCaloSD () const

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

	~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
double	getEnergyDeposit (const G4Step *) override

bool	getFromLibrary (const G4Step *) override

Protected Member Functions inherited from CaloSD
bool	checkHit ()

CaloG4Hit *	createNewHit (const G4Step *)

virtual void	endEvent ()

virtual bool	filterHit (CaloG4Hit *, double)

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

virtual uint16_t	getDepth (const G4Step *)

int	getNumberOfHits ()

double	getResponseWt (const G4Track *)

virtual int	getTrackID (const G4Track *)

bool	hitExists (const G4Step *)

virtual void	initEvent (const BeginOfEvent *)

virtual void	initRun ()

void	processHit (const G4Step *step)

void	resetForNewPrimary (const G4Step *)

void	setNumberCheckedHits (int val)

void	setParameterized (bool val)

G4ThreeVector	setToGlobal (const G4ThreeVector &, const G4VTouchable *) const

G4ThreeVector	setToLocal (const G4ThreeVector &, const G4VTouchable *) const

virtual int	setTrackID (const G4Step *)

void	setUseMap (bool val)

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
TrackInformation *	cmsTrackInformation (const G4Track *aTrack)

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
uint32_t	rotateUnitID (uint32_t, const G4Track *, const CastorShowerEvent &)

Private Attributes
double	energyThresholdSL

G4LogicalVolume *	lvC3EF

G4LogicalVolume *	lvC3HF

G4LogicalVolume *	lvC4EF

G4LogicalVolume *	lvC4HF

G4LogicalVolume *	lvCAST

double	non_compensation_factor

CastorNumberingScheme *	numberingScheme

CastorShowerLibrary *	showerLibrary

bool	useShowerLibrary

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

Protected Attributes inherited from CaloSD
CaloG4Hit *	currentHit

CaloHitID	currentID

float	edepositEM

float	edepositHAD

double	eminHit

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

bool	forceSave

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

G4ThreeVector	posGlobal

CaloHitID	previousID

bool	suppressHeavy

double	tmaxHit

Detailed Description

Description: Stores hits of Castor in appropriate container

Usage: Used in sensitive detector builder

Definition at line 30 of file CastorSD.h.

Constructor & Destructor Documentation

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

Definition at line 31 of file CastorSD.cc.

References g4SimHits_cfi::CastorShowerLibrary, energyThresholdSL, edm::ParameterSet::getParameter(), GeV, lvC3EF, lvC3HF, lvC4EF, lvC4HF, lvCAST, non_compensation_factor, setNumberingScheme(), CaloSD::setParameterized(), showerLibrary, and useShowerLibrary.

                                                    : 
   CaloSD(name, cpv, clg, p, manager), numberingScheme(nullptr), lvC3EF(nullptr),
   lvC3HF(nullptr), lvC4EF(nullptr), lvC4HF(nullptr), lvCAST(nullptr) {
   
   edm::ParameterSet m_CastorSD = p.getParameter<edm::ParameterSet>("CastorSD");
   useShowerLibrary  = m_CastorSD.getParameter<bool>("useShowerLibrary");
   energyThresholdSL = m_CastorSD.getParameter<double>("minEnergyInGeVforUsingSLibrary");
   energyThresholdSL = energyThresholdSL*GeV;   //  Convert GeV => MeV 
           
   non_compensation_factor = m_CastorSD.getParameter<double>("nonCompensationFactor");
   
   if (useShowerLibrary) {
     showerLibrary = new CastorShowerLibrary(name, p);
     setParameterized(true);
   }
   setNumberingScheme(new CastorNumberingScheme());
   
   edm::LogInfo("ForwardSim") 
     << "********************************************************\n"
     << "* Constructing a CastorSD  with name " << GetName() << "\n"
     << "* Using Castor Shower Library: " << useShowerLibrary << "\n"
     << "********************************************************";
     
   const G4LogicalVolumeStore * lvs = G4LogicalVolumeStore::GetInstance();
   for (auto lv : *lvs) {
     if (strcmp((lv->GetName()).c_str(),"C3EF") == 0) { lvC3EF = lv; }
     if (strcmp((lv->GetName()).c_str(),"C3HF") == 0) { lvC3HF = lv; }
     if (strcmp((lv->GetName()).c_str(),"C4EF") == 0) { lvC4EF = lv; }
     if (strcmp((lv->GetName()).c_str(),"C4HF") == 0) { lvC4HF = lv; }
     if (strcmp((lv->GetName()).c_str(),"CAST") == 0) { lvCAST = lv; }
     if (lvC3EF != nullptr && lvC3HF != nullptr && lvC4EF != nullptr && 
         lvC4HF != nullptr && lvCAST != nullptr) { break; }
   }
   edm::LogInfo("ForwardSim") << "CastorSD:: LogicalVolume pointers\n"
                              << lvC3EF << " for C3EF; " << lvC3HF 
                              << " for C3HF; " << lvC4EF << " for C4EF; " 
                              << lvC4HF << " for C4HF; " 
                              << lvCAST << " for CAST. ";
 }

CastorSD::~CastorSD ( )

override

Definition at line 76 of file CastorSD.cc.

References showerLibrary.

                     {
   delete showerLibrary;
 }

Member Function Documentation

double CastorSD::getEnergyDeposit ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 82 of file CastorSD.cc.

References a, funct::abs(), pfBoostedDoubleSVAK8TagInfos_cfi::beta, ALCARECOTkAlJpsiMuMu_cff::charge, edmIntegrityCheck::d, TrackInformation::getCastorHitPID(), TrackInformation::hasCastorHit(), G4TrackToParticleID::isStableHadronIon(), lvC3EF, lvC3HF, lvC4EF, lvC4HF, SiStripPI::max, min(), non_compensation_factor, pi, alignCSCRings::r, Scenarios_cff::scale, TrackInformation::setCastorHitPID(), mathSSE::sqrt(), and funct::tan().

                                                       {
   
   double NCherPhot = 0.;
 
   // Get theTrack 
   auto const theTrack = aStep->GetTrack();
 
   // preStepPoint information *********************************************
   auto const preStepPoint = aStep->GetPreStepPoint();
   auto const currentPV    = preStepPoint->GetPhysicalVolume();
   auto const currentLV    = currentPV->GetLogicalVolume();
 
 #ifdef debugLog
   edm::LogInfo("ForwardSim") << "CastorSD::getEnergyDeposit:"
                              << "\n CurrentStepNumber , TrackID , ParentID, Particle , VertexPosition ,"
                              << " LogicalVolumeAtVertex , PV, Time"
                              << "\n  TRACKINFO: " 
                              << theTrack->GetCurrentStepNumber() 
                              << " , " 
                              << theTrack->GetTrackID() 
                              << " , " 
                              << theTrack->GetParentID() 
                              << " , "
                              << theTrack->GetDefinition()->GetParticleName() 
                              << " , "
                              << theTrack->GetVertexPosition() 
                              << " , "
                              << theTrack->GetLogicalVolumeAtVertex()->GetName() 
                              << " , "
                              << currentPV->GetName()
                              << " , " 
                              << theTrack->GetGlobalTime();
 #endif
 
   // if particle moves from interaction point or "backwards (halo)
   const G4ThreeVector&  hit_mom  = preStepPoint->GetMomentumDirection();
   const G4ThreeVector&  hitPoint = preStepPoint->GetPosition();        
   double zint = hitPoint.z();
   
   // Check if theTrack is a muon (if so, DO NOT use Shower Library) 
   G4int parCode = theTrack->GetDefinition()->GetPDGEncoding();
   bool isHad = G4TrackToParticleID::isStableHadronIon(theTrack);
 
   double meanNCherPhot=0;
   G4double charge = preStepPoint->GetCharge();
   // VI: no Cerenkov light from neutrals
   if(0.0 == charge) { return meanNCherPhot; }
 
   G4double beta = preStepPoint->GetBeta();
   const double bThreshold = 0.67;
   // VI: no Cerenkov light from non-relativistic particles
   if(beta < bThreshold) { return meanNCherPhot; }
     
   // remember primary particle hitting the CASTOR detector
   TrackInformationExtractor TIextractor;
   TrackInformation& trkInfo = TIextractor(theTrack);
   if (!trkInfo.hasCastorHit()) {
     trkInfo.setCastorHitPID(parCode);
   }
   G4double stepl = aStep->GetStepLength()/cm;
 
   //int castorHitPID = trkInfo.hasCastorHit() ? std::abs(trkInfo.getCastorHitPID())
   //  : std::abs(parCode);
 
   // *************************************************
   // take into account light collection curve for plate
   //      double weight = curve_Castor(nameVolume, preStepPoint);
   //      double edep   = aStep->GetTotalEnergyDeposit() * weight;
   // *************************************************
   
   
   // *************************************************
   /*    comments for sensitive volumes:      
         C001 ...-... CP06,CBU1 ...-...CALM --- > fibres and bundle 
         for first release of CASTOR
         CASF  --- > quartz plate  for first and second releases of CASTOR  
         GF2Q, GFNQ, GR2Q, GRNQ 
         for tests with my own test geometry of HF (on ask of Gavrilov)
         C3TF, C4TF - for third release of CASTOR
   */  
 #ifdef debugLog
   edm::LogInfo("ForwardSim") << "CastorSD::getEnergyDeposit: for ID=" 
                              << theTrack->GetTrackID() << " LV: " << currentLV->GetName() 
                              << " isHad:" << isHad << " pdg=" << parCode 
                              << " castorPID=" << trkInfo.getCastorHitPID()
                              << " sl= " << stepl << " Edep= " << aStep->GetTotalEnergyDeposit(); 
 #endif
   if (currentLV == lvC3EF || currentLV == lvC4EF || currentLV == lvC3HF ||
       currentLV == lvC4HF) {
     
     const double nMedium = 1.4925;
     //     double photEnSpectrDL = (1./400.nm-1./700.nm)*10000000.cm/nm; /* cm-1  */
     //     double photEnSpectrDL = 10714.285714;
     
     const double photEnSpectrDE = 1.24;    /* see below   */
     /*     E = 2pi*(1./137.)*(eV*cm/370.)/lambda  =     */
     /*       = 12.389184*(eV*cm)/lambda                 */
     /*     Emax = 12.389184*(eV*cm)/400nm*10-7cm/nm  = 3.01 eV     */
     /*     Emin = 12.389184*(eV*cm)/700nm*10-7cm/nm  = 1.77 eV     */
     /*     delE = Emax - Emin = 1.24 eV  --> */
     /*   */
     /* default for Castor nameVolume  == "CASF" or (C3TF & C4TF)  */
     
     const double thFullRefl = 23.;  /* 23.dergee */
     const double thFullReflRad = thFullRefl*pi/180.;
     
     /* default for Castor nameVolume  == "CASF" or (C3TF & C4TF)  */
     const double thFibDir = 45.;  /* .dergee */
     /* for test HF geometry volumes:   
        if(nameVolume == "GF2Q" || nameVolume == "GFNQ" ||
        nameVolume == "GR2Q" || nameVolume == "GRNQ")
        thFibDir = 0.0; // .dergee
     */
     const double thFibDirRad = thFibDir*pi/180.;
     
     // theta of charged particle in LabRF(hit momentum direction):
     double costh =hit_mom.z()/sqrt(hit_mom.x()*hit_mom.x()+
                                    hit_mom.y()*hit_mom.y()+
                                    hit_mom.z()*hit_mom.z());
     if (zint < 0) costh = -costh;
     double th = acos(std::min(std::max(costh,double(-1.)),double(1.)));
     
     // just in case (can do bot use):
     if (th < 0.) th += twopi;
     
     // theta of cone with Cherenkov photons w.r.t.direction of charged part.:
     double costhcher =1./(nMedium*beta);
     double thcher = acos(std::min(std::max(costhcher,double(-1.)),double(1.)));
     
     // diff thetas of charged part. and quartz direction in LabRF:
     double DelFibPart = std::abs(th - thFibDirRad);
     
     // define real distances:
     double d = std::abs(tan(th)-tan(thFibDirRad));   
     
     double a = tan(thFibDirRad)+tan(std::abs(thFibDirRad-thFullReflRad));   
     double r = tan(th)+tan(std::abs(th-thcher));       
     
     // define losses d_qz in cone of full reflection inside quartz direction
     double d_qz;
 #ifdef debugLog
     int variant(0);
 #endif
     if(DelFibPart > (thFullReflRad + thcher) ) {
       d_qz = 0.; 
     }
     else {
       if((th + thcher) < (thFibDirRad+thFullReflRad) && 
          (th - thcher) > (thFibDirRad-thFullReflRad)) {
         d_qz = 1.; 
 #ifdef debugLog
         variant=1;
 #endif
       }
       else {
         if((thFibDirRad + thFullReflRad) < (th + thcher) && 
            (thFibDirRad - thFullReflRad) > (th - thcher) ) {
           d_qz = 0.; 
 #ifdef debugLog
           variant=2;
 #endif
         } else {
 #ifdef debugLog
           variant=3;
 #endif          
           // use crossed length of circles(cone projection)
           // dC1/dC2 : 
           double arg_arcos = 0.;
           double tan_arcos = 2.*a*d;
           if(tan_arcos != 0.) arg_arcos =(r*r-a*a-d*d)/tan_arcos; 
           arg_arcos = std::abs(arg_arcos);
           double th_arcos = acos(std::min(std::max(arg_arcos,-1.),1.));
           d_qz = std::abs(th_arcos/CLHEP::twopi);
         }
       }
     }
         
     // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     
     meanNCherPhot = 370.*charge*charge*
       ( 1. - 1./(nMedium*nMedium*beta*beta) )*photEnSpectrDE*stepl;
       
     double scale = isHad ? non_compensation_factor : 1.0;
 
     int poissNCherPhot = std::max(0, (int)G4Poisson(meanNCherPhot * scale));
       
     const double effPMTandTransport = 0.19;
     const double ReflPower = 0.1;
     double proba = d_qz + (1-d_qz)*ReflPower;
     NCherPhot = poissNCherPhot*effPMTandTransport*proba*0.307;
 #ifdef debugLog
       edm::LogInfo("ForwardSim") << " Nph= " << NCherPhot << " Np= " << poissNCherPhot
                                  << " eff= " << effPMTandTransport << " pb= " << proba
                                  << " Nmean= " << meanNCherPhot
                                  << " q=" << charge << " beta=" << beta 
                                  << " nMedium= " << nMedium << " sl= " << stepl
                                  << " Nde=" << photEnSpectrDE;
 #endif
   }
   return NCherPhot;
 }

bool CastorSD::getFromLibrary ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 363 of file CastorSD.cc.

References funct::abs(), pfBoostedDoubleSVAK8TagInfos_cfi::beta, ALCARECOTkAlJpsiMuMu_cff::charge, CaloSD::currentID, dot(), CaloSD::edepositEM, CaloSD::edepositHAD, energyThresholdSL, CastorShowerEvent::getDetID(), CastorShowerEvent::getNhit(), CastorShowerEvent::getNphotons(), CastorShowerEvent::getPrimE(), CastorShowerLibrary::getShowerHits(), CastorShowerEvent::getTime(), GeV, TrackInformation::hasCastorHit(), hfClusterShapes_cfi::hits, mps_fire::i, G4TrackToParticleID::isGammaElectronPositron(), G4TrackToParticleID::isStableHadronIon(), lvC3EF, lvC3HF, lvC4EF, lvC4HF, lvCAST, M_PI, non_compensation_factor, CaloSD::processHit(), dttmaxenums::R, CaloSD::resetForNewPrimary(), rotateUnitID(), Scenarios_cff::scale, TrackInformation::setCastorHitPID(), CaloHitID::setID(), CaloSD::setTrackID(), showerLibrary, mathSSE::sqrt(), theta(), ntuplemaker::time, and useShowerLibrary.

                                                  {
 
 //
 //   Method to get hits from the Shower Library
 //    
 //   "updateHit" save the Hits to a CaloG4Hit container
 //
 
   auto const theTrack = aStep->GetTrack();
   auto parCode = theTrack->GetDefinition()->GetPDGEncoding();
 
   // remember primary particle hitting the CASTOR detector    
   TrackInformationExtractor TIextractor;
   TrackInformation& trkInfo = TIextractor(theTrack);
   if (!trkInfo.hasCastorHit()) {
     trkInfo.setCastorHitPID(parCode);
   }
 
   if (!useShowerLibrary) { return false; }
 
   // preStepPoint information *********************************************
   
   auto const preStepPoint = aStep->GetPreStepPoint();
   auto const currentPV    = preStepPoint->GetPhysicalVolume();
   auto const currentLV    = currentPV->GetLogicalVolume();
 
 #ifdef debugLog
     edm::LogInfo("ForwardSim") << "CastorSD::getFromLibrary: for ID=" << theTrack->GetTrackID() 
                                << " parentID= " << theTrack->GetParentID() << " " 
                                << theTrack->GetDefinition()->GetParticleName() 
                                << " LV: " << currentLV->GetName() << " PV: " << currentPV->GetName()
                                << "\n eta= " << theTrack->GetPosition().eta() 
                                << " phi= " << theTrack->GetPosition().phi()
                                << " z(cm)= " << theTrack->GetPosition().z()/cm 
                                << " time(ns)= " << theTrack->GetGlobalTime() 
                                << " E(GeV)= " << theTrack->GetTotalEnergy()/GeV;
   
 #endif
   
   const G4ThreeVector&  hitPoint = preStepPoint->GetPosition();        
   const G4ThreeVector&  hit_mom  = preStepPoint->GetMomentumDirection();
   double zint = hitPoint.z();
   double pz   = hit_mom.z();
   
   // Check if theTrack moves backward 
   bool backward = (pz * zint < 0.) ? true : false;
   
   // Check that theTrack is above the energy threshold to use Shower Library 
   bool aboveThreshold = (theTrack->GetKineticEnergy() > energyThresholdSL) 
     ? true : false;
     
   // Check if theTrack is a muon (if so, DO NOT use Shower Library) 
   bool isEM = G4TrackToParticleID::isGammaElectronPositron(parCode);
   bool isHad = G4TrackToParticleID::isStableHadronIon(theTrack);
   
   // angle condition
   double theta_max = M_PI - 3.1305; // angle in radians corresponding to -5.2 eta
   double R_mom=sqrt(hit_mom.x()*hit_mom.x() + hit_mom.y()*hit_mom.y());
   double theta = atan2(R_mom,std::abs(pz));
   bool angleok = (theta < theta_max) ? true : false;
   
   // OkToUse
   double R = sqrt(hitPoint.x()*hitPoint.x() + hitPoint.y()*hitPoint.y());
   bool dot = (zint < -14450. && R < 45.) ? true : false;
   bool inRange = (zint < -14700. || R > 193.) ? false : true;
   //bool OkToUse = (inRange && !dot) ? true : false;
   
   bool particleWithinShowerLibrary = aboveThreshold && (isEM || isHad)
     && (!backward) && inRange && !dot && angleok && currentLV == lvCAST;
 
 #ifdef debugLog
   edm::LogInfo("ForwardSim") << "CastorSD::getFromLibrary: ID= " << theTrack->GetTrackID()
                              << " E>E0 " << aboveThreshold
                              << " isEM " << isEM << " isHad " << isHad 
                  << " backword " << backward
                              << " Ok " << (inRange && !dot) << " angle " << angleok << " LV: " 
                              << currentLV->GetName() << "  " << (currentLV == lvCAST) 
                              << " " << particleWithinShowerLibrary 
                              << " Edep= " << aStep->GetTotalEnergyDeposit();
 #endif
 
   // Use Castor shower library if energy is above threshold, is not a muon 
   // and is not moving backward 
   if (!particleWithinShowerLibrary) {
     //---VI: This code is for backward compatibility and should be removed 
 
     if (currentLV == lvC3EF || currentLV == lvC4EF || currentLV == lvC3HF ||
         currentLV == lvC4HF) {
       G4double edep     = aStep->GetTotalEnergyDeposit();
       G4double beta     = preStepPoint->GetBeta();
       G4double charge   = preStepPoint->GetCharge();
       double bThreshold = 0.67;
       if(edep == 0.0 && charge != 0.0 && beta > bThreshold) { 
         G4Poisson(0.0); 
       }
     }
     //---
     return false;
   }
   
   // ****    Call method to retrieve hits from the ShowerLibrary   ****
   // always kill primary
   bool isKilled(true);
   CastorShowerEvent hits = showerLibrary->getShowerHits(aStep, isKilled);
   
   //---VI:
   //if(!isKilled) { return false; }
   
   int primaryID = setTrackID(aStep);
 
   // Reset entry point for new primary
   resetForNewPrimary(aStep);
 
 #ifdef debugLog
   edm::LogInfo("ForwardSim") 
     << "CastorSD::getFromLibrary:  " 
     << hits.getNhit() << " hits for " << GetName() << " from " 
     << theTrack->GetDefinition()->GetParticleName() << " of "
     << preStepPoint->GetKineticEnergy()/GeV << " GeV and trackID " 
     << theTrack->GetTrackID() << " isHad: " << isHad;
 #endif
 
   // Scale to correct energy
   double E_track = preStepPoint->GetTotalEnergy() ;
   double E_SLhit = hits.getPrimE() * CLHEP::GeV ;
   double scale = E_track/E_SLhit ;
 
   //Non compensation 
   if (isHad){
     scale *= non_compensation_factor; // if hadronic extend the scale with the non-compensation factor
   }
   //  Loop over hits retrieved from the library
   for (unsigned int i=0; i<hits.getNhit(); ++i) {
     
     // Get nPhotoElectrons and set edepositEM / edepositHAD accordingly
     double nPhotoElectrons = hits.getNphotons(i)*scale;
 
     if(isEM)  {
        edepositEM  = nPhotoElectrons; 
        edepositHAD = 0.;                 
     } else {
        edepositEM  = 0.;                  
        edepositHAD = nPhotoElectrons;
     }
     
     // Get hit position and time
     double                time = hits.getTime(i);
     
     // Get hit detID
     unsigned int        unitID = hits.getDetID(i);
     
     // Make the detID "rotation" from one sector to another taking into account the 
     // sectors of the impinging particle (theTrack) and of the particle that produced 
     // the 'hits' retrieved from shower library   
     unsigned int rotatedUnitID = rotateUnitID(unitID , theTrack , hits);
     currentID.setID(rotatedUnitID, time, primaryID, 0);
     processHit(aStep);
   }  //  End of loop over hits
   return isKilled;
 }

uint32_t CastorSD::rotateUnitID	(	uint32_t	unitID,
		const G4Track *	track,
		const CastorShowerEvent &	shower
	)

private

Definition at line 304 of file CastorSD.cc.

References printConversionInfo::aux, CastorShowerEvent::getPrimPhi(), createfilelist::int, LogDebug, M_PI, and reco::btau::trackPhi.

Referenced by getFromLibrary().

                                                                                                       {
 // ==============================================================
 //
 //   o   Exploit Castor phi symmetry to return newUnitID for  
 //       shower hits based on track phi coordinate 
 //
 // ==============================================================
   
   // Get 'track' phi:
   double   trackPhi = track->GetPosition().phi(); 
   if(trackPhi<0) trackPhi += 2*M_PI ;
   // Get phi from primary that gave rise to SL 'shower':
   double  showerPhi = shower.getPrimPhi(); 
   if(showerPhi<0) showerPhi += 2*M_PI ;
   // Delta phi:
   
   //  Find the OctSector for which 'track' and 'shower' belong
   int  trackOctSector = (int) (  trackPhi / (M_PI/4) ) ;
   int showerOctSector = (int) ( showerPhi / (M_PI/4) ) ;
   
   uint32_t  newUnitID;
   uint32_t         sec = ( ( unitID>>4 ) & 0xF ) ;
   uint32_t  complement = ( unitID & 0xFFFFFF0F ) ;
   
   // Get 'track' z:
   double   trackZ = track->GetPosition().z();
   
   int aux ;
   int dSec = 2*(trackOctSector - showerOctSector) ;
   if(trackZ>0)  // Good for revision 1.9 of CastorNumberingScheme
   {
     int sec1 = sec-dSec;
     //    sec -= dSec ;
     if(sec1<0) sec1  += 16;
     if(sec1>15) sec1 -= 16;
     sec = (uint32_t)(sec1);
   } else {
     if( dSec<0 ) sec += 16 ;
     sec += dSec ;
     aux  = (int) (sec/16) ;
     sec -= aux*16 ;
   }
   sec  = sec<<4 ;
   newUnitID = complement | sec ;
   
 #ifdef debugLog
   if(newUnitID != unitID) {
     LogDebug("ForwardSim") << "\n CastorSD::rotateUnitID:  " 
                            << "\n     unitID = " << unitID 
                            << "\n  newUnitID = " << newUnitID ; 
   }
 #endif
   
   return newUnitID ;
 
 }

uint32_t CastorSD::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 286 of file CastorSD.cc.

References CastorNumberingScheme::getUnitID(), and numberingScheme.

                                                    {
   return (numberingScheme == nullptr ? 0 : numberingScheme->getUnitID(aStep));
 }

void CastorSD::setNumberingScheme ( CastorNumberingScheme * scheme )

Definition at line 292 of file CastorSD.cc.

References numberingScheme.

Referenced by CastorSD().

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