#include <HCalSD.h>

Inheritance diagram for HCalSD:

Public Member Functions
virtual double	getEnergyDeposit (G4Step *)

	HCalSD (G4String, const DDCompactView &, SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

virtual bool	ProcessHits (G4Step , G4TouchableHistory )

virtual uint32_t	setDetUnitId (G4Step *step)

void	setNumberingScheme (HcalNumberingScheme *)

virtual	~HCalSD ()

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

Protected Member Functions
virtual bool	filterHit (CaloG4Hit *, double)

virtual void	initRun ()

Protected Member Functions inherited from CaloSD
G4bool	checkHit ()

virtual void	clearHits ()

CaloG4Hit *	createNewHit ()

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

virtual uint16_t	getDepth (G4Step *)

double	getResponseWt (G4Track *)

virtual G4bool	getStepInfo (G4Step *aStep)

virtual int	getTrackID (G4Track *)

G4bool	hitExists ()

void	resetForNewPrimary (G4ThreeVector, double)

G4ThreeVector	setToLocal (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...

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

void	getFromLibrary (G4Step *step)

void	getFromParam (G4Step *step)

void	getHitFibreBundle (G4Step *step, bool type)

void	getHitPMT (G4Step *step)

std::vector< G4String >	getNames (DDFilteredView &)

void	hitForFibre (G4Step *step)

bool	isItConicalBundle (G4LogicalVolume *)

bool	isItFibre (G4LogicalVolume *)

bool	isItFibre (G4String)

bool	isItHF (G4Step *)

bool	isItHF (G4String)

bool	isItPMT (G4LogicalVolume *)

bool	isItScintillator (G4Material *)

bool	isItStraightBundle (G4LogicalVolume *)

double	layerWeight (int, G4ThreeVector, int, int)

void	plotProfile (G4Step *step, G4ThreeVector pos, double edep, double time, int id)

void	readWeightFromFile (std::string)

uint32_t	setDetUnitId (int, G4ThreeVector, int, int)

int	setTrackID (G4Step *step)

Private Attributes
double	betaThr

double	birk1

double	birk2

double	birk3

TH1F *	dist_ [9]

double	eminHitHB

double	eminHitHE

double	eminHitHF

double	eminHitHO

std::vector< G4LogicalVolume * >	fibre1LV

std::vector< G4LogicalVolume * >	fibre2LV

std::vector< G4LogicalVolume * >	fibreLV

std::vector< G4String >	fibreNames

std::vector< int >	hfLevels

std::vector< G4LogicalVolume * >	hfLV

std::vector< G4String >	hfNames

HFShower *	hfshower

TH1F *	hit_ [9]

std::vector< double >	layer0wt

std::map< uint32_t, double >	layerWeights

std::vector< G4Material * >	materials

std::vector< G4String >	matNames

G4int	mumPDG

G4int	mupPDG

HcalNumberingFromDDD *	numberingFromDDD

HcalNumberingScheme *	numberingScheme

std::vector< G4LogicalVolume * >	pmtLV

HFShowerFibreBundle *	showerBundle

HFShowerLibrary *	showerLibrary

HFShowerParam *	showerParam

HFShowerPMT *	showerPMT

TH1F *	time_ [9]

bool	useBirk

bool	useFibreBundle

bool	useHF

bool	useLayerWt

bool	useParam

bool	usePMTHit

bool	useShowerLibrary

Additional Inherited Members
Public 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 31 of file HCalSD.h.

Constructor & Destructor Documentation

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

Definition at line 32 of file HCalSD.cc.

References DDFilteredView::addFilter(), betaThr, birk1, birk2, birk3, dist_, eminHitHB, eminHitHE, eminHitHF, eminHitHO, DDSpecificsFilter::equals, fibre1LV, fibre2LV, fibreLV, fibreNames, mergeVDriftHistosByStation::file, DDFilteredView::firstChild(), g, getDDDArray(), SensitiveDetector::getNames(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hfLevels, hfLV, hfNames, hfshower, hit_, i, edm::Service< T >::isAvailable(), isItFibre(), isItHF(), CaloSD::kmaxIon, CaloSD::kmaxNeutron, CaloSD::kmaxProton, L1TDQM_cfg::labels, layer0wt, testEve_cfg::level, funct::log(), LogDebug, DDFilteredView::logicalPart(), TFileDirectory::make(), DDLogicalPart::material(), materials, matNames, DDFilteredView::mergedSpecifics(), mumPDG, mupPDG, DDName::name(), SensitiveDetector::name, DDBase< N, C >::name(), DDFilteredView::next(), numberingFromDDD, pmtLV, readWeightFromFile(), DDSpecificsFilter::setCriteria(), setNumberingScheme(), showerBundle, showerLibrary, showerParam, showerPMT, CaloSD::suppressHeavy, cond::rpcobtemp::temp, time_, indexGen::title, useBirk, useFibreBundle, useHF, useLayerWt, useParam, usePMTHit, useShowerLibrary, and relativeConstraints::value.

                                                                           : 
   CaloSD(name, cpv, clg, p, manager, 
      p.getParameter<edm::ParameterSet>("HCalSD").getParameter<int>("TimeSliceUnit"),
      p.getParameter<edm::ParameterSet>("HCalSD").getParameter<bool>("IgnoreTrackID")), 
   numberingFromDDD(0), numberingScheme(0), showerLibrary(0), hfshower(0), 
   showerParam(0), showerPMT(0), showerBundle(0){
 
   //static SimpleConfigurable<bool>   on1(false, "HCalSD:UseBirkLaw");
   //static SimpleConfigurable<double> bk1(0.013, "HCalSD:BirkC1");
   //static SimpleConfigurable<double> bk2(0.0568,"HCalSD:BirkC2");
   //static SimpleConfigurable<double> bk3(1.75,  "HCalSD:BirkC3");
   // Values from NIM 80 (1970) 239-244: as implemented in Geant3
   //static SimpleConfigurable<bool> on2(true,"HCalSD:UseShowerLibrary");
 
   edm::ParameterSet m_HC = p.getParameter<edm::ParameterSet>("HCalSD");
   useBirk    = m_HC.getParameter<bool>("UseBirkLaw");
   birk1      = m_HC.getParameter<double>("BirkC1")*(g/(MeV*cm2));
   birk2      = m_HC.getParameter<double>("BirkC2");
   birk3      = m_HC.getParameter<double>("BirkC3");
   useShowerLibrary = m_HC.getParameter<bool>("UseShowerLibrary");
   useParam         = m_HC.getParameter<bool>("UseParametrize");
   bool testNumber  = m_HC.getParameter<bool>("TestNumberingScheme");
   usePMTHit        = m_HC.getParameter<bool>("UsePMTHits");
   betaThr          = m_HC.getParameter<double>("BetaThreshold");
   eminHitHB        = m_HC.getParameter<double>("EminHitHB")*MeV;
   eminHitHE        = m_HC.getParameter<double>("EminHitHE")*MeV;
   eminHitHO        = m_HC.getParameter<double>("EminHitHO")*MeV;
   eminHitHF        = m_HC.getParameter<double>("EminHitHF")*MeV;
   useFibreBundle   = m_HC.getParameter<bool>("UseFibreBundleHits");
   useHF            = m_HC.getUntrackedParameter<bool>("UseHF",true);
   bool forTBH2     = m_HC.getUntrackedParameter<bool>("ForTBH2",false);
   useLayerWt       = m_HC.getUntrackedParameter<bool>("UseLayerWt",false);
   std::string file = m_HC.getUntrackedParameter<std::string>("WtFile","None");
 
 #ifdef DebugLog
   LogDebug("HcalSim") << "***************************************************" 
               << "\n"
               << "*                                                 *"
               << "\n"
               << "* Constructing a HCalSD  with name " << name << "\n"
               << "*                                                 *"
               << "\n"
               << "***************************************************";
 #endif
   edm::LogInfo("HcalSim") << "HCalSD:: Use of HF code is set to " << useHF
               << "\nUse of shower parametrization set to "
               << useParam << "\nUse of shower library is set to " 
               << useShowerLibrary << "\nUse PMT Hit is set to "
               << usePMTHit << " with beta Threshold "<< betaThr
               << "\nUSe of FibreBundle Hit set to "<<useFibreBundle
               << "\n         Use of Birks law is set to      " 
               << useBirk << "  with three constants kB = "
               << birk1 << ", C1 = " << birk2 << ", C2 = " << birk3;
   edm::LogInfo("HcalSim") << "HCalSD:: Suppression Flag " << suppressHeavy
               << " protons below " << kmaxProton << " MeV,"
               << " neutrons below " << kmaxNeutron << " MeV and"
               << " ions below " << kmaxIon << " MeV\n"
               << "         Threshold for storing hits in HB: "
               << eminHitHB << " HE: " << eminHitHE << " HO: "
               << eminHitHO << " HF: " << eminHitHF;
 
   numberingFromDDD = new HcalNumberingFromDDD(name, cpv);
   HcalNumberingScheme* scheme;
   if (testNumber || forTBH2) 
     scheme = dynamic_cast<HcalNumberingScheme*>(new HcalTestNumberingScheme(forTBH2));
   else            
     scheme = new HcalNumberingScheme();
   setNumberingScheme(scheme);
 
   const G4LogicalVolumeStore * lvs = G4LogicalVolumeStore::GetInstance();
   std::vector<G4LogicalVolume *>::const_iterator lvcite;
   G4LogicalVolume* lv;
   std::string attribute, value;
   if (useHF) {
     if (useParam) showerParam = new HFShowerParam(name, cpv, p);
     else {
       if (useShowerLibrary) showerLibrary = new HFShowerLibrary(name, cpv, p);
       hfshower  = new HFShower(name, cpv, p, 0);
     }
 
     // HF volume names
     attribute = "Volume";
     value     = "HF";
     DDSpecificsFilter filter0;
     DDValue           ddv0(attribute,value,0);
     filter0.setCriteria(ddv0,DDSpecificsFilter::equals);
     DDFilteredView fv0(cpv);
     fv0.addFilter(filter0);
     hfNames = getNames(fv0);
     fv0.firstChild();
     DDsvalues_type sv0(fv0.mergedSpecifics());
     std::vector<double> temp =  getDDDArray("Levels",sv0);
     edm::LogInfo("HcalSim") << "HCalSD: Names to be tested for " << attribute 
                 << " = " << value << " has " << hfNames.size()
                 << " elements";
     for (unsigned int i=0; i<hfNames.size(); i++) {
       G4String namv = hfNames[i];
       lv            = 0;
       for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) 
     if ((*lvcite)->GetName() == namv) {
       lv = (*lvcite);
       break;
     }
        hfLV.push_back(lv);
       int level = static_cast<int>(temp[i]);
       hfLevels.push_back(level);
       edm::LogInfo("HcalSim") << "HCalSD:  HF[" << i << "] = " << hfNames[i]
                   << " LV " << hfLV[i] << " at level " 
                   << hfLevels[i];
     }
   
     // HF Fibre volume names
     value     = "HFFibre";
     DDSpecificsFilter filter1;
     DDValue           ddv1(attribute,value,0);
     filter1.setCriteria(ddv1,DDSpecificsFilter::equals);
     DDFilteredView fv1(cpv);
     fv1.addFilter(filter1);
     fibreNames = getNames(fv1);
     edm::LogInfo("HcalSim") << "HCalSD: Names to be tested for " << attribute 
                 << " = " << value << ":";
     for (unsigned int i=0; i<fibreNames.size(); i++) {
       G4String namv = fibreNames[i];
       lv            = 0;
       for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) 
     if ((*lvcite)->GetName() == namv) {
       lv = (*lvcite);
       break;
     }
       fibreLV.push_back(lv);
       edm::LogInfo("HcalSim") << "HCalSD:  (" << i << ") " << fibreNames[i]
                   << " LV " << fibreLV[i];
     }
   
     // HF PMT volume names
     value     = "HFPMT";
     DDSpecificsFilter filter3;
     DDValue           ddv3(attribute,value,0);
     filter3.setCriteria(ddv3,DDSpecificsFilter::equals);
     DDFilteredView fv3(cpv);
     fv3.addFilter(filter3);
     std::vector<G4String> pmtNames = getNames(fv3);
     edm::LogInfo("HcalSim") << "HCalSD: Names to be tested for " << attribute
                 << " = " << value << " have " << pmtNames.size()
                 << " entries";
     for (unsigned int i=0; i<pmtNames.size(); i++) {
       G4String namv = pmtNames[i];
       lv            = 0;
       for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) 
     if ((*lvcite)->GetName() == namv) {
       lv = (*lvcite);
       break;
     }
       pmtLV.push_back(lv);
       edm::LogInfo("HcalSim") << "HCalSD:  (" << i << ") " << pmtNames[i]
                   << " LV " << pmtLV[i];
     }
     if (pmtNames.size() > 0) showerPMT = new HFShowerPMT (name, cpv, p);
   
     // HF Fibre bundles
     value     = "HFFibreBundleStraight";
     DDSpecificsFilter filter4;
     DDValue           ddv4(attribute,value,0);
     filter4.setCriteria(ddv4,DDSpecificsFilter::equals);
     DDFilteredView fv4(cpv);
     fv4.addFilter(filter4);
     std::vector<G4String> fibreNames = getNames(fv4);
     edm::LogInfo("HcalSim") << "HCalSD: Names to be tested for " << attribute
                 << " = " << value << " have " << fibreNames.size()
                 << " entries";
     for (unsigned int i=0; i<fibreNames.size(); i++) {
       G4String namv = fibreNames[i];
       lv            = 0;
       for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) 
     if ((*lvcite)->GetName() == namv) {
       lv = (*lvcite);
       break;
     }
       fibre1LV.push_back(lv);
       edm::LogInfo("HcalSim") << "HCalSD:  (" << i << ") " << fibreNames[i]
                   << " LV " << fibre1LV[i];
     }
 
     value     = "HFFibreBundleConical";
     DDSpecificsFilter filter5;
     DDValue           ddv5(attribute,value,0);
     filter5.setCriteria(ddv5,DDSpecificsFilter::equals);
     DDFilteredView fv5(cpv);
     fv5.addFilter(filter5);
     fibreNames = getNames(fv5);
     edm::LogInfo("HcalSim") << "HCalSD: Names to be tested for " << attribute
                 << " = " << value << " have " << fibreNames.size()
                 << " entries";
     for (unsigned int i=0; i<fibreNames.size(); i++) {
       G4String namv = fibreNames[i];
       lv            = 0;
       for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) 
     if ((*lvcite)->GetName() == namv) {
       lv = (*lvcite);
       break;
     }
       fibre2LV.push_back(lv);
       edm::LogInfo("HcalSim") << "HCalSD:  (" << i << ") " << fibreNames[i]
                   << " LV " << fibre2LV[i];
     }
     if (fibre1LV.size() > 0 || fibre2LV.size() > 0) 
       showerBundle = new HFShowerFibreBundle (name, cpv, p);
   }
 
   //Material list for HB/HE/HO sensitive detectors
   attribute = "ReadOutName";
   DDSpecificsFilter filter2;
   DDValue           ddv2(attribute,name,0);
   filter2.setCriteria(ddv2,DDSpecificsFilter::equals);
   DDFilteredView fv2(cpv);
   fv2.addFilter(filter2);
   bool dodet = fv2.firstChild();
 
   DDsvalues_type sv(fv2.mergedSpecifics());
   //Layer0 Weight
   layer0wt = getDDDArray("Layer0Wt",sv);
   edm::LogInfo("HcalSim") << "HCalSD: " << layer0wt.size() << " Layer0Wt";
   for (unsigned int it=0; it<layer0wt.size(); it++)
     edm::LogInfo("HcalSim") << "HCalSD: [" << it << "] " << layer0wt[it];
 
   const G4MaterialTable * matTab = G4Material::GetMaterialTable();
   std::vector<G4Material*>::const_iterator matite;
   while (dodet) {
     const DDLogicalPart & log = fv2.logicalPart();
     G4String namx = log.name().name();
     if (!isItHF(namx) && !isItFibre(namx)) {
       bool notIn = true;
       for (unsigned int i=0; i<matNames.size(); i++) {
     if (!strcmp(matNames[i].c_str(),log.material().name().name().c_str())){
       notIn = false;
       break;
     }
       }
       if (notIn) {
     namx = log.material().name().name();
     matNames.push_back(namx);
     G4Material* mat;
     for (matite = matTab->begin(); matite != matTab->end(); matite++) {
       if ((*matite)->GetName() == namx) {
         mat = (*matite);
         break;
       }
     }
     materials.push_back(mat);
       }
     }
     dodet = fv2.next();
   }
 
   edm::LogInfo("HcalSim") << "HCalSD: Material names for " << attribute 
               << " = " << name << ":";
   for (unsigned int i=0; i<matNames.size(); i++)
     edm::LogInfo("HcalSim") << "HCalSD: (" << i << ") " << matNames[i]
                 << " pointer " << materials[i];
 
   mumPDG = mupPDG = 0;
   
   if (useLayerWt) readWeightFromFile(file);
 
   for (int i=0;  i<9; i++) {
     hit_[i] = time_[i]= dist_[i] = 0;
   }
 
 #ifdef DebugLog
   edm::Service<TFileService> tfile;
 
   if ( tfile.isAvailable() ) {
     static std::string labels[9] = {"HB", "HE", "HO", "HF Absorber", "HF PMT", "HF Absorber Long", "HF Absorber Short", "HF PMT Long", "HF PMT Short"};
     char name[20], title[60];
     for (int i=0; i<9; i++) {
       sprintf (title, "Hit energy in %s", labels[i].c_str());
       sprintf (name, "HCalSDHit%d", i);
       hit_[i] = tfile->make<TH1F>(name, title, 2000, 0., 2000.);
       sprintf (title, "Energy (MeV)");
       hit_[i]->GetXaxis()->SetTitle(title);
       hit_[i]->GetYaxis()->SetTitle("Hits");
       sprintf (title, "Time of the hit in %s", labels[i].c_str());
       sprintf (name, "HCalSDTime%d", i);
       time_[i] = tfile->make<TH1F>(name, title, 2000, 0., 2000.);
       sprintf (title, "Time (ns)");
       time_[i]->GetXaxis()->SetTitle(title);
       time_[i]->GetYaxis()->SetTitle("Hits");
       sprintf (title, "Longitudinal profile in %s", labels[i].c_str());
       sprintf (name, "HCalSDDist%d", i);
       dist_[i] = tfile->make<TH1F>(name, title, 2000, 0., 2000.);
       sprintf (title, "Distance (mm)");
       dist_[i]->GetXaxis()->SetTitle(title);
       dist_[i]->GetYaxis()->SetTitle("Hits");
     }
   }
 #endif
 }

HCalSD::~HCalSD ( )

virtual

Definition at line 332 of file HCalSD.cc.

References hfshower, numberingFromDDD, numberingScheme, showerBundle, showerLibrary, showerParam, and showerPMT.

                 { 
   if (numberingFromDDD) delete numberingFromDDD;
   if (numberingScheme)  delete numberingScheme;
   if (showerLibrary)    delete showerLibrary;
   if (hfshower)         delete hfshower;
   if (showerParam)      delete showerParam;
   if (showerPMT)        delete showerPMT;
   if (showerBundle)     delete showerBundle;
 }

Member Function Documentation

bool HCalSD::filterHit	(	CaloG4Hit *	aHit,
		double	time
	)

protectedvirtual

Reimplemented from CaloSD.

Definition at line 511 of file HCalSD.cc.

References eminHitHB, eminHitHE, eminHitHF, eminHitHO, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getUnitID(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, dtT0WireCalibration_cfg::threshold, and CaloSD::tmaxHit.

                                                    {
 
   double threshold=0;
   DetId theId(aHit->getUnitID());
   switch (theId.subdetId()) {
   case HcalBarrel:
     threshold = eminHitHB; break;
   case HcalEndcap:
     threshold = eminHitHE; break;
   case HcalOuter:
     threshold = eminHitHO; break;
   case HcalForward:
     threshold = eminHitHF; break;
   default:
     break;
   }
   return ((time <= tmaxHit) && (aHit->getEnergyDeposit() > threshold));
 }

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

private

Definition at line 545 of file HCalSD.cc.

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

Referenced by HCalSD().

                                                {
 
 #ifdef DebugLog
   LogDebug("HcalSim") << "HCalSD:getDDDArray called for " << str;
 #endif
   DDValue value(str);
   if (DDfetch(&sv,value)) {
 #ifdef DebugLog
     LogDebug("HcalSim") << value;
 #endif
     const std::vector<double> & fvec = value.doubles();
     int nval = fvec.size();
     if (nval < 1) {
       edm::LogError("HcalSim") << "HCalSD : # of " << str << " bins " << nval
                    << " < 2 ==> illegal";
       throw cms::Exception("Unknown", "HCalSD")
     << "nval < 2 for array " << str << "\n";
     }
     
     return fvec;
   } else {
     edm::LogError("HcalSim") << "HCalSD :  cannot get array " << str;
     throw cms::Exception("Unknown", "HCalSD")
       << "cannot get array " << str << "\n";
   }
 }

double HCalSD::getEnergyDeposit ( G4Step * aStep )

virtual

Reimplemented from CaloSD.

Definition at line 423 of file HCalSD.cc.

References birk1, birk2, birk3, CaloSD::getAttenuation(), CaloSD::getResponseWt(), isItScintillator(), TrackInformation::isPrimary(), CaloSD::kmaxIon, CaloSD::kmaxNeutron, CaloSD::kmaxProton, layer0wt, layerWeight(), CaloSD::suppressHeavy, CaloSD::theTrack, useBirk, useLayerWt, and CommonMethods::weight().

                                              {
 
   double destep = aStep->GetTotalEnergyDeposit();
   double weight = 1;
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   int depth = (touch->GetReplicaNumber(0))%10;
   int det   = ((touch->GetReplicaNumber(1))/1000)-3;
   if (depth==0 && (det==0 || det==1)) weight = layer0wt[det];
   if (useLayerWt) {
     int lay   = (touch->GetReplicaNumber(0)/10)%100 + 1;
     G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
     weight = layerWeight(det+3, hitPoint, depth, lay);
   }
 
   if (suppressHeavy) {
     G4Track* theTrack = aStep->GetTrack();
     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) 
       edm::LogInfo("HcalSim") << "HCalSD:Ignore Track " 
                   << theTrack->GetTrackID() << " Type " 
                   << theTrack->GetDefinition()->GetParticleName()
                   << " Kinetic Energy " << ke << " MeV";
 #endif
       }
     }
   }
 #ifdef DebugLog
   double weight0 = weight;
 #endif
   if (useBirk) {
     G4Material* mat = aStep->GetPreStepPoint()->GetMaterial();
     if (isItScintillator(mat))
       weight *= getAttenuation(aStep, birk1, birk2, birk3);
   }
   double wt1    = getResponseWt(theTrack);
 #ifdef DebugLog
   edm::LogInfo("HcalSim") << "HCalSD: Detector " << det+3 << " Depth " << depth
               << " weight " << weight0 << " " << weight << " " 
               << wt1; 
 #endif
   return weight*wt1*destep;
 }

void HCalSD::getFromLibrary ( G4Step * step )

private

Definition at line 661 of file HCalSD.cc.

References CaloSD::checkHit(), CaloSD::createNewHit(), CaloSD::currentHit, CaloSD::currentID, CaloSD::edepositEM, CaloSD::edepositHAD, CaloSD::emPDG, CaloSD::epPDG, CaloSD::gammaPDG, HFShowerLibrary::getHits(), GetVolume(), i, convertSQLiteXML::ok, plotProfile(), CaloSD::posGlobal, CaloSD::preStepPoint, CaloSD::previousID, CaloSD::resetForNewPrimary(), setDetUnitId(), CaloHitID::setID(), setTrackID(), showerLibrary, CaloSD::theTrack, cond::rpcobgas::time, and CaloSD::updateHit().

Referenced by ProcessHits().

                                           {
 
   preStepPoint  = aStep->GetPreStepPoint(); 
   theTrack      = aStep->GetTrack();   
   int det       = 5;
   bool ok;
 
   std::vector<HFShowerLibrary::Hit> hits = showerLibrary->getHits(aStep, ok);
 
   double etrack    = preStepPoint->GetKineticEnergy();
   int    primaryID = setTrackID(aStep);
   /*
   int    primaryID = 0;
   if (etrack >= energyCut) {
     primaryID    = theTrack->GetTrackID();
   } else {
     primaryID    = theTrack->GetParentID();
     if (primaryID == 0) primaryID = theTrack->GetTrackID();
   }
   */
 
   // Reset entry point for new primary
   posGlobal = preStepPoint->GetPosition();
   resetForNewPrimary(posGlobal, etrack);
 
   G4int particleCode = theTrack->GetDefinition()->GetPDGEncoding();
   if (particleCode==emPDG || particleCode==epPDG || particleCode==gammaPDG) {
     edepositEM  = 1.*GeV; edepositHAD = 0.;
   } else {
     edepositEM  = 0.; edepositHAD = 1.*GeV;
   }
 
 #ifdef DebugLog
   edm::LogInfo("HcalSim") << "HCalSD::getFromLibrary " <<hits.size() 
               << " hits for " << GetName() << " of " << primaryID 
               << " with " << theTrack->GetDefinition()->GetParticleName() 
               << " of " << preStepPoint->GetKineticEnergy()/GeV 
               << " GeV";
 #endif
   for (unsigned int i=0; i<hits.size(); i++) {
     G4ThreeVector hitPoint = hits[i].position;
     int depth              = hits[i].depth;
     double time            = hits[i].time;
     unsigned int unitID    = setDetUnitId(det, hitPoint, depth);
     currentID.setID(unitID, time, primaryID, 0);
 #ifdef DebugLog
     plotProfile(aStep, hitPoint, 1.0*GeV, time, depth);
 #endif
    
     // check if it is in the same unit and timeslice as the previous one
     if (currentID == previousID) {
       updateHit(currentHit);
     } else {
       if (!checkHit()) currentHit = createNewHit();
     }
   }
 
   //Now kill the current track
   if (ok) {
     theTrack->SetTrackStatus(fStopAndKill);
     G4TrackVector tv = *(aStep->GetSecondary());
     for (unsigned int kk=0; kk<tv.size(); kk++) {
       if (tv[kk]->GetVolume() == preStepPoint->GetPhysicalVolume())
     tv[kk]->SetTrackStatus(fStopAndKill);
     }
   }
 }

void HCalSD::getFromParam ( G4Step * step )

private

Definition at line 777 of file HCalSD.cc.

References CaloSD::checkHit(), CaloSD::createNewHit(), CaloSD::currentHit, CaloSD::currentID, CaloSD::edepositEM, CaloSD::edepositHAD, HFShowerParam::getHits(), i, plotProfile(), CaloSD::posGlobal, CaloSD::preStepPoint, CaloSD::previousID, setDetUnitId(), CaloHitID::setID(), setTrackID(), showerParam, cond::rpcobgas::time, and CaloSD::updateHit().

Referenced by ProcessHits().

                                         {
 
   std::vector<HFShowerParam::Hit> hits = showerParam->getHits(aStep);
   int nHit = static_cast<int>(hits.size());
 
   if (nHit > 0) {
     preStepPoint  = aStep->GetPreStepPoint();
     int primaryID = setTrackID(aStep);
    
     int det   = 5;
 #ifdef DebugLog
     edm::LogInfo("HcalSim") << "HCalSD::getFromParam " << nHit << " hits for " 
                 << GetName() << " of " << primaryID << " with " 
                 <<  aStep->GetTrack()->GetDefinition()->GetParticleName()
                 << " of " << preStepPoint->GetKineticEnergy()/GeV 
                 << " GeV in detector type " << det;
 #endif
     for (int i = 0; i<nHit; i++) {
       G4ThreeVector hitPoint = hits[i].position;
       int depth              = hits[i].depth;
       double time            = hits[i].time;
       unsigned int unitID    = setDetUnitId(det, hitPoint, depth);
       currentID.setID(unitID, time, primaryID, 0);
       edepositEM             = hits[i].edep*GeV; 
       edepositHAD            = 0;
 #ifdef DebugLog
       plotProfile(aStep, hitPoint, edepositEM, time, depth);
 #endif
 
       // check if it is in the same unit and timeslice as the previous one
       if (currentID == previousID) {
         updateHit(currentHit);
       } else {
         posGlobal = preStepPoint->GetPosition();
         if (!checkHit()) currentHit = createNewHit();
       }
     }
   }
 }

void HCalSD::getHitFibreBundle	(	G4Step *	step,
		bool	type
	)

private

Definition at line 882 of file HCalSD.cc.

References beta, CaloSD::checkHit(), CaloSD::createNewHit(), CaloSD::currentHit, CaloSD::currentID, CaloSD::edepositEM, CaloSD::edepositHAD, CaloSD::energyCut, HFShowerFibreBundle::getHits(), HFShowerFibreBundle::getRadius(), HcalNumberingScheme::getUnitID(), HcalForward, LogDebug, numberingFromDDD, numberingScheme, phi, plotProfile(), CaloSD::posGlobal, CaloSD::preStepPoint, CaloSD::previousID, CaloSD::resetForNewPrimary(), CaloHitID::setID(), showerBundle, CaloSD::theTrack, cond::rpcobgas::time, tmp, HcalNumberingFromDDD::unitID(), and CaloSD::updateHit().

Referenced by ProcessHits().

                                                         {
 
   preStepPoint = aStep->GetPreStepPoint();
   theTrack     = aStep->GetTrack();
   double edep  = showerBundle->getHits(aStep, type);
 
   if (edep >= 0) {
     double etrack    = preStepPoint->GetKineticEnergy();
     int    primaryID = 0;
     if (etrack >= energyCut) {
       primaryID    = theTrack->GetTrackID();
     } else {
       primaryID    = theTrack->GetParentID();
       if (primaryID == 0) primaryID = theTrack->GetTrackID();
     }
     // Reset entry point for new primary
     posGlobal = preStepPoint->GetPosition();
     resetForNewPrimary(posGlobal, etrack);
 
     //
     int    det      = static_cast<int>(HcalForward);
     G4ThreeVector hitPoint = preStepPoint->GetPosition();   
     double rr       = (hitPoint.x()*hitPoint.x() + hitPoint.y()*hitPoint.y());
     double phi      = (rr == 0. ? 0. :atan2(hitPoint.y(),hitPoint.x()));
     double etaR     = showerBundle->getRadius();
     int depth       = 1;
     if (etaR < 0) {
       depth         = 2;
       etaR          =-etaR;
     }
     if (hitPoint.z() < 0) etaR =-etaR;
 #ifdef DebugLog
     edm::LogInfo("HcalSim") << "HCalSD::Hit for Detector " << det << " etaR "
                 << etaR << " phi " << phi/deg << " depth " <<depth;
 #endif
     double time = (aStep->GetPostStepPoint()->GetGlobalTime());
     uint32_t unitID = 0;
     if (numberingFromDDD) {
       HcalNumberingFromDDD::HcalID tmp = numberingFromDDD->unitID(det,etaR,phi,
                                   depth,1);
       if (numberingScheme) unitID = numberingScheme->getUnitID(tmp);
     }
     if (type) currentID.setID(unitID, time, primaryID, 3);
     else      currentID.setID(unitID, time, primaryID, 2);
 
     edepositHAD = aStep->GetTotalEnergyDeposit();
     edepositEM  =-edepositHAD + (edep*GeV);
 #ifdef DebugLog
     plotProfile(aStep, hitPoint, edep*GeV, time, depth);
     double beta = preStepPoint->GetBeta();
     LogDebug("HcalSim") << "HCalSD::getHitFibreBundle 1 hit for " << GetName() 
                         << " of " << primaryID << " with " 
             << theTrack->GetDefinition()->GetParticleName()
                         << " of " << preStepPoint->GetKineticEnergy()/GeV 
                         << " GeV with velocity " << beta << " UnitID "
                         << std::hex << unitID << std::dec;
 #endif
     // check if it is in the same unit and timeslice as the previous one
     if (currentID == previousID) {
       updateHit(currentHit);
     } else {
       if (!checkHit()) currentHit = createNewHit();
     }
   } // non-zero energy deposit
 }

void HCalSD::getHitPMT ( G4Step * step )

private

Definition at line 817 of file HCalSD.cc.

References beta, CaloSD::checkHit(), CaloSD::createNewHit(), CaloSD::currentHit, CaloSD::currentID, CaloSD::edepositEM, CaloSD::edepositHAD, CaloSD::energyCut, HFShowerPMT::getHits(), HFShowerPMT::getRadius(), HcalNumberingScheme::getUnitID(), HcalForward, LogDebug, numberingFromDDD, numberingScheme, phi, plotProfile(), CaloSD::posGlobal, CaloSD::preStepPoint, CaloSD::previousID, CaloSD::resetForNewPrimary(), CaloHitID::setID(), showerPMT, CaloSD::theTrack, cond::rpcobgas::time, tmp, HcalNumberingFromDDD::unitID(), and CaloSD::updateHit().

Referenced by ProcessHits().

                                      {
 
   preStepPoint = aStep->GetPreStepPoint();
   theTrack     = aStep->GetTrack();
   double edep  = showerPMT->getHits(aStep);
 
   if (edep >= 0) {
     double etrack    = preStepPoint->GetKineticEnergy();
     int    primaryID = 0;
     if (etrack >= energyCut) {
       primaryID    = theTrack->GetTrackID();
     } else {
       primaryID    = theTrack->GetParentID();
       if (primaryID == 0) primaryID = theTrack->GetTrackID();
     }
     // Reset entry point for new primary
     posGlobal = preStepPoint->GetPosition();
     resetForNewPrimary(posGlobal, etrack);
 
     //
     int    det      = static_cast<int>(HcalForward);
     G4ThreeVector hitPoint = preStepPoint->GetPosition();   
     double rr       = (hitPoint.x()*hitPoint.x() + hitPoint.y()*hitPoint.y());
     double phi      = (rr == 0. ? 0. :atan2(hitPoint.y(),hitPoint.x()));
     double etaR     = showerPMT->getRadius();
     int depth       = 1;
     if (etaR < 0) {
       depth         = 2;
       etaR          =-etaR;
     }
     if (hitPoint.z() < 0) etaR =-etaR;
 #ifdef DebugLog
     edm::LogInfo("HcalSim") << "HCalSD::Hit for Detector " << det << " etaR "
                 << etaR << " phi " << phi/deg << " depth " <<depth;
 #endif
     double time = (aStep->GetPostStepPoint()->GetGlobalTime());
     uint32_t unitID = 0;
     if (numberingFromDDD) {
       HcalNumberingFromDDD::HcalID tmp = numberingFromDDD->unitID(det,etaR,phi,
                                   depth,1);
       if (numberingScheme) unitID = numberingScheme->getUnitID(tmp);
     }
     currentID.setID(unitID, time, primaryID, 1);
 
     edepositHAD = aStep->GetTotalEnergyDeposit();
     edepositEM  =-edepositHAD + (edep*GeV);
 #ifdef DebugLog
     plotProfile(aStep, hitPoint, edep*GeV, time, depth);
     double beta = preStepPoint->GetBeta();
     LogDebug("HcalSim") << "HCalSD::getHitPMT 1 hit for " << GetName() 
                         << " of " << primaryID << " with " 
             << theTrack->GetDefinition()->GetParticleName()
                         << " of " << preStepPoint->GetKineticEnergy()/GeV 
                         << " GeV with velocity " << beta << " UnitID "
                         << std::hex << unitID << std::dec;
 #endif
     // check if it is in the same unit and timeslice as the previous one
     if (currentID == previousID) {
       updateHit(currentHit);
     } else {
       if (!checkHit()) currentHit = createNewHit();
     }
   }
 }

std::vector< G4String > HCalSD::getNames ( DDFilteredView & fv )

private

Definition at line 573 of file HCalSD.cc.

References DDFilteredView::firstChild(), i, funct::log(), DDFilteredView::logicalPart(), DDName::name(), DDBase< N, C >::name(), DDFilteredView::next(), convertSQLiteXML::ok, and tmp.

                                                        {
 
   std::vector<G4String> tmp;
   bool dodet = fv.firstChild();
   while (dodet) {
     const DDLogicalPart & log = fv.logicalPart();
     bool ok = true;
     for (unsigned int i=0; i<tmp.size(); i++) {
       if (!strcmp(tmp[i].c_str(),log.name().name().c_str())) {
     ok = false;
     break;
       }
     }
     if (ok) tmp.push_back(log.name().name());
     dodet = fv.next();
   }
   return tmp;
 }

void HCalSD::hitForFibre ( G4Step * step )

private

Definition at line 729 of file HCalSD.cc.

References CaloSD::checkHit(), CaloSD::createNewHit(), CaloSD::currentHit, CaloSD::currentID, CaloSD::edepositEM, CaloSD::edepositHAD, CaloSD::emPDG, CaloSD::epPDG, CaloSD::gammaPDG, HFShower::getHits(), hfshower, i, plotProfile(), CaloSD::posGlobal, CaloSD::preStepPoint, CaloSD::previousID, setDetUnitId(), CaloHitID::setID(), setTrackID(), CaloSD::theTrack, cond::rpcobgas::time, and CaloSD::updateHit().

Referenced by ProcessHits().

                                        {
 
   preStepPoint  = aStep->GetPreStepPoint();
   theTrack      = aStep->GetTrack();
   int primaryID = setTrackID(aStep);
 
   int det   = 5;
   std::vector<HFShower::Hit> hits = hfshower->getHits(aStep);
 
   G4int particleCode = theTrack->GetDefinition()->GetPDGEncoding();
   if (particleCode==emPDG || particleCode==epPDG || particleCode==gammaPDG) {
     edepositEM  = 1.*GeV; edepositHAD = 0.;
   } else {
     edepositEM  = 0.; edepositHAD = 1.*GeV;
   }
  
 #ifdef DebugLog
   edm::LogInfo("HcalSim") << "HCalSD::hitForFibre " << hits.size() 
               << " hits for " << GetName() << " of " << primaryID 
               << " with " << theTrack->GetDefinition()->GetParticleName() 
               << " of " << preStepPoint->GetKineticEnergy()/GeV 
               << " GeV in detector type " << det;
 #endif
   if (hits.size() > 0) {
 
     G4ThreeVector hitPoint = preStepPoint->GetPosition();
     int           depth    = 
       (preStepPoint->GetTouchable()->GetReplicaNumber(0))%10;
     unsigned int unitID    = setDetUnitId(det, hitPoint, depth);
 
     for (unsigned int i=0; i<hits.size(); i++) {
       double time            = hits[i].time;
       currentID.setID(unitID, time, primaryID, 0);
 #ifdef DebugLog
       plotProfile(aStep, hitPoint, 1.0*GeV, time, depth);
 #endif
       // check if it is in the same unit and timeslice as the previous one
       if (currentID == previousID) {
         updateHit(currentHit);
       } else {
         posGlobal = preStepPoint->GetPosition();
         if (!checkHit()) currentHit = createNewHit();
       }
     }
   }
 
 }

void HCalSD::initRun ( )

protectedvirtual

Reimplemented from CaloSD.

Definition at line 497 of file HCalSD.cc.

References HFShowerParam::initRun(), HFShowerLibrary::initRun(), mumPDG, mupPDG, showerLibrary, and showerParam.

                      {
 
   G4ParticleTable * theParticleTable = G4ParticleTable::GetParticleTable();
   G4String          particleName;
   mumPDG = theParticleTable->FindParticle(particleName="mu-")->GetPDGEncoding();
   mupPDG = theParticleTable->FindParticle(particleName="mu+")->GetPDGEncoding();
 #ifdef DebugLog
   edm::LogInfo("HcalSim") << "HCalSD: Particle code for mu- = " << mumPDG
               << " for mu+ = " << mupPDG;
 #endif
   if (showerLibrary) showerLibrary->initRun(theParticleTable);
   if (showerParam)   showerParam->initRun(theParticleTable);
 }

bool HCalSD::isItConicalBundle ( G4LogicalVolume * lv )

private

Definition at line 645 of file HCalSD.cc.

References fibre2LV.

Referenced by ProcessHits().

                                                    {
 
   std::vector<G4LogicalVolume*>::const_iterator ite = fibre2LV.begin();
   for (; ite != fibre2LV.end(); ite++) 
     if (lv == *ite) return true;
   return false;
 }

bool HCalSD::isItFibre ( G4LogicalVolume * lv )

private

Definition at line 613 of file HCalSD.cc.

References fibreLV.

Referenced by HCalSD(), and ProcessHits().

                                            {
 
   std::vector<G4LogicalVolume*>::const_iterator ite = fibreLV.begin();
   for (; ite != fibreLV.end(); ite++) 
     if (lv == *ite) return true;
   return false;
 }

bool HCalSD::isItFibre ( G4String name )

private

Definition at line 621 of file HCalSD.cc.

References fibreNames.

                                      {
 
   std::vector<G4String>::const_iterator it = fibreNames.begin();
   for (; it != fibreNames.end(); it++) 
     if (name == *it) return true;
   return false;
 }

bool HCalSD::isItHF ( G4Step * aStep )

private

Definition at line 592 of file HCalSD.cc.

References hfLevels, hfLV, and hfNames.

Referenced by HCalSD(), and ProcessHits().

                                    {
 
   const G4VTouchable* touch    = aStep->GetPreStepPoint()->GetTouchable();
   int levels = ((touch->GetHistoryDepth())+1);
   for (unsigned int it=0; it < hfNames.size(); it++) {
     if (levels >= hfLevels[it]) {
       G4LogicalVolume* lv = touch->GetVolume(levels-hfLevels[it])->GetLogicalVolume();
       if (lv == hfLV[it]) return true;
     }
   }
   return false;
 }

bool HCalSD::isItHF ( G4String name )

private

Definition at line 605 of file HCalSD.cc.

References hfNames.

                                   {
 
   std::vector<G4String>::const_iterator it = hfNames.begin();
   for (; it != hfNames.end(); it++) 
     if (name == *it) return true;
   return false;
 }

bool HCalSD::isItPMT ( G4LogicalVolume * lv )

private

Definition at line 629 of file HCalSD.cc.

References pmtLV.

Referenced by ProcessHits().

                                          {
 
   std::vector<G4LogicalVolume*>::const_iterator ite = pmtLV.begin();
   for (; ite != pmtLV.end(); ite++) 
     if (lv == *ite) return true;
   return false;
 }

bool HCalSD::isItScintillator ( G4Material * mat )

private

Definition at line 653 of file HCalSD.cc.

References materials.

Referenced by getEnergyDeposit().

                                               {
 
   std::vector<G4Material*>::const_iterator ite = materials.begin();
   for (; ite != materials.end(); ite++) 
     if (mat == *ite) return true;
   return false;
 }

bool HCalSD::isItStraightBundle ( G4LogicalVolume * lv )

private

Definition at line 637 of file HCalSD.cc.

References fibre1LV.

Referenced by ProcessHits().

                                                     {
 
   std::vector<G4LogicalVolume*>::const_iterator ite = fibre1LV.begin();
   for (; ite != fibre1LV.end(); ite++) 
     if (lv == *ite) return true;
   return false;
 }

double HCalSD::layerWeight	(	int	det,
		G4ThreeVector	pos,
		int	depth,
		int	lay
	)

private

Definition at line 995 of file HCalSD.cc.

References HcalNumberingFromDDD::HcalID::etaR, HcalNumberingFromDDD::HcalID::lay, layerWeights, numberingFromDDD, HcalTestNumbering::packHcalIndex(), HcalNumberingFromDDD::HcalID::phis, HcalNumberingFromDDD::HcalID::subdet, tmp, HcalNumberingFromDDD::unitID(), and HcalNumberingFromDDD::HcalID::zside.

Referenced by getEnergyDeposit().

                                                                          { 
 
   double wt = 1.;
   if (numberingFromDDD) {
     //get the ID's as eta, phi, depth, ... indices
     HcalNumberingFromDDD::HcalID tmp = numberingFromDDD->unitID(det, pos, 
                                 depth, lay);
     uint32_t id = HcalTestNumbering::packHcalIndex(tmp.subdet, tmp.zside,
                            1, tmp.etaR, tmp.phis,
                            tmp.lay);
     std::map<uint32_t,double>::const_iterator ite = layerWeights.find(id);
     if (ite != layerWeights.end()) wt = ite->second;
 #ifdef DebugLog
     edm::LogInfo("HcalSim") << "HCalSD::layerWeight: ID " << std::hex << id 
                 << std::dec << " (" << tmp.subdet << "/"  
                 << tmp.zside << "/1/" << tmp.etaR << "/" 
                 << tmp.phis << "/"  << tmp.lay << ") Weight " <<wt;
 #endif
   }
   return wt;
 }

void HCalSD::plotProfile	(	G4Step *	step,
		G4ThreeVector	pos,
		double	edep,
		double	time,
		int	id
	)

private

Definition at line 1017 of file HCalSD.cc.

References abs, dist_, newFWLiteAna::found, hit_, LogDebug, n, SensitiveDetector::name, cond::rpcobgas::time, and time_.

Referenced by getFromLibrary(), getFromParam(), getHitFibreBundle(), getHitPMT(), hitForFibre(), and ProcessHits().

                                   { 
 
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   static G4String modName[8] = {"HEModule", "HVQF", "HBModule", "MBAT", "MBBT", "MBBTC", "MBBT_R1P", "MBBT_R1M"};
   G4ThreeVector local;
   bool found=false;
   double depth=-2000;
   int idx = 4;
   for (int n=0; n<touch->GetHistoryDepth(); n++) {
     G4String name = touch->GetVolume(n)->GetName();
 #ifdef DebugLog
     LogDebug("HcalSim") << "plotProfile Depth " << n << " Name " << name;
 #endif
     for (unsigned int ii=0; ii<8; ii++) {
       if (name == modName[ii]) {
     found = true;
     int dn = touch->GetHistoryDepth() - n;
     local = touch->GetHistory()->GetTransform(dn).TransformPoint(global);
     if      (ii == 0) {depth = local.z() - 4006.5; idx = 1;}
     else if (ii == 1) {depth = local.z() + 825.0;  idx = 3;}
     else if (ii == 2) {depth = local.x() - 1775.;  idx = 0;}
     else              {depth = local.y() + 15.;    idx = 2;}
     break;
       }
     }
     if (found) break;
   }
   if (!found) depth = std::abs(global.z()) - 11500;
 #ifdef DebugLog
   LogDebug("HcalSim") << "plotProfile Found " << found << " Global " << global << " Local " << local << " depth " << depth << " ID " << id << " EDEP " << edep << " Time " << time;
 #endif
   if (hit_[idx]  != 0) hit_[idx]->Fill(edep);
   if (time_[idx] != 0) time_[idx]->Fill(time,edep);
   if (dist_[idx] != 0) dist_[idx]->Fill(depth,edep);
   int jd = 2*idx + id - 7;
   if (jd >= 0 && jd < 4) {
     jd += 5;
     if (hit_[jd]  != 0) hit_[jd]->Fill(edep);
     if (time_[jd] != 0) time_[jd]->Fill(time,edep);
     if (dist_[jd] != 0) dist_[jd]->Fill(depth,edep);
   }
 }

bool HCalSD::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *
	)

virtual

Reimplemented from CaloSD.

Definition at line 342 of file HCalSD.cc.

References CaloSD::createNewHit(), CaloSD::currentHit, CaloSD::edepositEM, CaloSD::edepositHAD, getFromLibrary(), getFromParam(), getHitFibreBundle(), getHitPMT(), CaloSD::getStepInfo(), CaloSD::hitExists(), hitForFibre(), isItConicalBundle(), isItFibre(), isItHF(), isItPMT(), isItStraightBundle(), LogDebug, mumPDG, mupPDG, SensitiveDetector::NaNTrap(), NULL, plotProfile(), showerBundle, showerPMT, useFibreBundle, useParam, usePMTHit, and useShowerLibrary.

                                                               {
   //  TimeMe t1( theHitTimer, false);
   
   NaNTrap( aStep ) ;
   
   if (aStep == NULL) {
     return true;
   } else {
     G4LogicalVolume* lv = 
       aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume();
     G4String nameVolume = lv->GetName();
     if (isItHF(aStep)) {
       G4int parCode =aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
       if (useParam) {
 #ifdef DebugLog
         LogDebug("HcalSim") << "HCalSD: Hit from parametrization in " 
                             << nameVolume << " for Track " 
                             << aStep->GetTrack()->GetTrackID()
                             <<" (" << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
         getFromParam(aStep);
       } else {
         bool notaMuon = true;
         if (parCode == mupPDG || parCode == mumPDG ) notaMuon = false;
         if (useShowerLibrary && notaMuon) {
 #ifdef DebugLog
           LogDebug("HcalSim") << "HCalSD: Starts shower library from " 
                               << nameVolume << " for Track " 
                               << aStep->GetTrack()->GetTrackID()
                               <<" (" << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
           getFromLibrary(aStep);
         } else if (isItFibre(lv)) {
 #ifdef DebugLog
           LogDebug("HcalSim") << "HCalSD: Hit at Fibre in " << nameVolume 
                               << " for Track " 
                               << aStep->GetTrack()->GetTrackID()
                               <<" ("  << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
           hitForFibre(aStep);
         }
       }
     } else if (isItPMT(lv)) {
 #ifdef DebugLog
       LogDebug("HcalSim") << "HCalSD: Hit from PMT parametrization from " 
                           <<  nameVolume << " for Track " 
                           << aStep->GetTrack()->GetTrackID() << " ("
                           << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
       if (usePMTHit && showerPMT) getHitPMT(aStep);
     } else if (isItStraightBundle(lv) || isItConicalBundle(lv)) {
 #ifdef DebugLog
       LogDebug("HcalSim") << "HCalSD: Hit from FibreBundle from "
                           << nameVolume << " for Track " 
                           << aStep->GetTrack()->GetTrackID() << " ("
                           << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
       if (useFibreBundle && showerBundle) 
     getHitFibreBundle(aStep, isItConicalBundle(lv));
     } else {
 #ifdef DebugLog
       LogDebug("HcalSim") << "HCalSD: Hit from standard path from " 
                           <<  nameVolume << " for Track " 
                           << aStep->GetTrack()->GetTrackID() << " ("
                           << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
       if (getStepInfo(aStep)) {
 #ifdef DebugLog
     if (edepositEM+edepositHAD > 0) 
       plotProfile(aStep, aStep->GetPreStepPoint()->GetPosition(), 
               edepositEM+edepositHAD, 
               aStep->GetPostStepPoint()->GetGlobalTime(), 0);
 #endif
         if (hitExists() == false && edepositEM+edepositHAD>0.) 
           currentHit = createNewHit();
       }
     }
     return true;
   }
 } 

void HCalSD::readWeightFromFile ( std::string fName )

private

Definition at line 969 of file HCalSD.cc.

References recoMuon::in, EdgesToViz::infile, layerWeights, HcalTestNumbering::packHcalIndex(), phi, and useLayerWt.

Referenced by HCalSD().

                                                {
 
   std::ifstream infile;
   int entry=0;
   infile.open(fName.c_str(), std::ios::in);
   if (infile) {
     int    det, zside, etaR, phi, lay;
     double wt;
     while (infile >> det >> zside >> etaR >> phi >> lay >> wt) {
       uint32_t id = HcalTestNumbering::packHcalIndex(det,zside,1,etaR,phi,lay);
       layerWeights.insert(std::pair<uint32_t,double>(id,wt));
       entry++;
 #ifdef DebugLog
       edm::LogInfo("HcalSim") << "HCalSD::readWeightFromFile:Entry " << entry
                   << " ID " << std::hex << id << std::dec << " ("
                   << det << "/" << zside << "/1/" << etaR << "/"
                   << phi << "/" << lay << ") Weight " << wt;
 #endif
     }
     infile.close();
   }
   edm::LogInfo("HcalSim") << "HCalSD::readWeightFromFile: reads " << entry
               << " weights from " << fName;
   if (entry <= 0) useLayerWt = false;
 }

uint32_t HCalSD::setDetUnitId ( G4Step * step )

virtual

Implements CaloSD.

Definition at line 475 of file HCalSD.cc.

References CaloSD::preStepPoint.

Referenced by getFromLibrary(), getFromParam(), and hitForFibre().

                                             { 
 
   G4StepPoint* preStepPoint = aStep->GetPreStepPoint(); 
   const G4VTouchable* touch = preStepPoint->GetTouchable();
   G4ThreeVector hitPoint    = preStepPoint->GetPosition();
 
   int depth = (touch->GetReplicaNumber(0))%10 + 1;
   int lay   = (touch->GetReplicaNumber(0)/10)%100 + 1;
   int det   = (touch->GetReplicaNumber(1))/1000;
 
   return setDetUnitId (det, hitPoint, depth, lay);
 }

uint32_t HCalSD::setDetUnitId	(	int	det,
		G4ThreeVector	pos,
		int	depth,
		int	lay = `1`
	)

private

Definition at line 531 of file HCalSD.cc.

References HcalNumberingScheme::getUnitID(), numberingFromDDD, numberingScheme, tmp, and HcalNumberingFromDDD::unitID().

                               { 
 
   uint32_t id = 0;
   if (numberingFromDDD) {
     //get the ID's as eta, phi, depth, ... indices
     HcalNumberingFromDDD::HcalID tmp = numberingFromDDD->unitID(det, pos, 
                                 depth, lay);
     //get the ID
     if (numberingScheme) id = numberingScheme->getUnitID(tmp);
   }
   return id;
 }

void HCalSD::setNumberingScheme ( HcalNumberingScheme * scheme )

Definition at line 488 of file HCalSD.cc.

References numberingScheme.

Referenced by HCalSD(), HcalTestAnalysis::update(), SimG4HcalValidation::update(), and HcalTB04Analysis::update().

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

int HCalSD::setTrackID ( G4Step * step )

private

Definition at line 948 of file HCalSD.cc.

References TrackInformation::getIDonCaloSurface(), CaloSD::preStepPoint, CaloSD::previousID, CaloSD::resetForNewPrimary(), CaloSD::theTrack, and CaloHitID::trackID().

Referenced by getFromLibrary(), getFromParam(), and hitForFibre().

                                      {
 
   theTrack     = aStep->GetTrack();
 
   double etrack = preStepPoint->GetKineticEnergy();
   TrackInformation * trkInfo = (TrackInformation *)(theTrack->GetUserInformation());
   int      primaryID = trkInfo->getIDonCaloSurface();
   if (primaryID == 0) {
 #ifdef DebugLog
     edm::LogInfo("HcalSim") << "HCalSD: Problem with primaryID **** set by "
                 << "force to TkID **** " <<theTrack->GetTrackID();
 #endif
     primaryID = theTrack->GetTrackID();
   }
 
   if (primaryID != previousID.trackID())
     resetForNewPrimary(preStepPoint->GetPosition(), etrack);
 
   return primaryID;
 }