#include <HcalTB06BeamSD.h>

Inheritance diagram for HcalTB06BeamSD:

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

uint32_t	setDetUnitId (const G4Step *step) override

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

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

virtual bool	getFromLibrary (const G4Step *step)

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
std::vector< G4String >	getNames (DDFilteredView &)

bool	isItWireChamber (const G4String &)

Private Attributes
double	birk1

double	birk2

double	birk3

G4String	matName

bool	useBirk

std::vector< G4String >	wcNames

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

Definition at line 20 of file HcalTB06BeamSD.h.

Constructor & Destructor Documentation

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

Definition at line 25 of file HcalTB06BeamSD.cc.

References birk1, birk2, birk3, DDFilteredView::firstChild(), g, SensitiveDetector::getNames(), edm::ParameterSet::getParameter(), mps_fire::i, cmsBatch::log, DDFilteredView::logicalPart(), DDLogicalPart::material(), matName, MeV, DDName::name(), DDBase< N, C >::name(), DDFilteredView::next(), AlCaHLTBitMon_QueryRunRegistry::string, useBirk, relativeConstraints::value, and wcNames.

                                                    : 
   CaloSD(name, cpv, clg, p, manager) {
 
   // Values from NIM 80 (1970) 239-244: as implemented in Geant3
   edm::ParameterSet m_HC = p.getParameter<edm::ParameterSet>("HcalTB06BeamSD");
   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");
 
   edm::LogInfo("HcalTB06BeamSD") << "HcalTB06BeamSD:: Use of Birks law is set "
                  << "to " << useBirk << "  with three "
                  << "constants kB = " << birk1 << ", C1 = " 
                  << birk2 << ", C2 = " << birk3;
 
   std::string attribute, value;
 
   // Wire Chamber volume names
   attribute = "Volume";
   value     = "WireChamber";
   DDSpecificsMatchesValueFilter filter1{DDValue(attribute,value,0)};
   DDFilteredView fv1(cpv,filter1);
   wcNames = getNames(fv1);
   edm::LogInfo("HcalTB06BeamSD") << "HcalTB06BeamSD:: Names to be tested for " 
                  << attribute << " = " << value << ": " 
                  << wcNames.size() << " paths";
   for (unsigned int i=0; i<wcNames.size(); i++)
     edm::LogInfo("HcalTB06BeamSD") << "HcalTB06BeamSD:: (" << i << ") " 
                    << wcNames[i];
 
   //Material list for scintillator detector
   attribute = "ReadOutName";
   DDSpecificsMatchesValueFilter filter2{DDValue(attribute,name,0)};
   DDFilteredView fv2(cpv,filter2);
   bool dodet = fv2.firstChild();
 
   std::vector<G4String> matNames;
   std::vector<int>      nocc;
   while (dodet) {
     const DDLogicalPart & log = fv2.logicalPart();
     matName = log.material().name().name();
     bool notIn = true;
     for (unsigned int i=0; i<matNames.size(); i++) {
       if (matName == matNames[i]) {notIn = false; nocc[i]++;}
     }
     if (notIn) {
       matNames.push_back(matName);
       nocc.push_back(0);
     }
     dodet = fv2.next();
   }
   if (!matNames.empty()) {
     matName = matNames[0];
     int occ = nocc[0];
     for (unsigned int i = 0; i < matNames.size(); i++) {
       if (nocc[i] > occ) {
     occ     = nocc[i];
     matName = matNames[i];
       }
     }
   } else {
     matName = "Not Found";
   }
   edm::LogInfo("HcalTB06BeamSD") << "HcalTB06BeamSD: Material name for " 
                  << attribute << " = " << name << ":" 
                  << matName;
 }

HcalTB06BeamSD::~HcalTB06BeamSD ( )

override

Definition at line 97 of file HcalTB06BeamSD.cc.

97 {}

Member Function Documentation

double HcalTB06BeamSD::getEnergyDeposit ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 99 of file HcalTB06BeamSD.cc.

References birk1, birk2, birk3, CaloSD::getAttenuation(), matName, useBirk, and mps_merge::weight.

                                                            {
 
   double destep = aStep->GetTotalEnergyDeposit();
   double weight = 1;
   if (useBirk && matName == aStep->GetPreStepPoint()->GetMaterial()->GetName())
     weight *= getAttenuation(aStep, birk1, birk2, birk3);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTB06BeamSD") << "HcalTB06BeamSD: Detector " 
                      << aStep->GetPreStepPoint()->GetTouchable()->GetVolume()->GetName()
                      << " weight " << weight;
 #endif
   return weight*destep;
 }

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

private

Definition at line 134 of file HcalTB06BeamSD.cc.

References DDFilteredView::firstChild(), mps_fire::i, cmsBatch::log, DDFilteredView::logicalPart(), DDName::name(), dataset::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();
     G4String name = log.name().name();
     bool ok = true;
     for (unsigned int i=0; i<tmp.size(); i++)
       if (name == tmp[i]) ok = false;
     if (ok) tmp.push_back(name);
     dodet = fv.next();
   }
   return tmp;
 }

bool HcalTB06BeamSD::isItWireChamber ( const G4String & name )

private

Definition at line 150 of file HcalTB06BeamSD.cc.

References wcNames.

Referenced by setDetUnitId().

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

uint32_t HcalTB06BeamSD::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 113 of file HcalTB06BeamSD.cc.

References createfilelist::int, isItWireChamber(), dataset::name, HcalTestBeamNumbering::packIndex(), CaloSD::setToLocal(), x, and y.

                                                           { 
 
   const G4StepPoint* preStepPoint = aStep->GetPreStepPoint(); 
   const G4VTouchable* touch = preStepPoint->GetTouchable();
   G4String name             = preStepPoint->GetPhysicalVolume()->GetName();
 
   int det = 1;
   int lay = 0, x = 0, y = 0;
   if (!isItWireChamber(name)) {
     lay     = (touch->GetReplicaNumber(0));
   } else {
     det = 2;
     lay = (touch->GetReplicaNumber(1));
     G4ThreeVector localPoint  = setToLocal(preStepPoint->GetPosition(), touch);
     x   = (int)(localPoint.x()/(0.2*mm));
     y   = (int)(localPoint.y()/(0.2*mm));
   }
 
   return HcalTestBeamNumbering::packIndex (det, lay, x, y);
 }