#include <HcalTB06BeamSD.h>

Inheritance diagram for HcalTB06BeamSD:

Public Member Functions
virtual double	getEnergyDeposit (G4Step *)

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

virtual uint32_t	setDetUnitId (G4Step *step)

virtual	~HcalTB06BeamSD ()

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

virtual void	clear ()

virtual void	DrawAll ()

virtual void	EndOfEvent (G4HCofThisEvent *eventHC)

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

virtual void	Initialize (G4HCofThisEvent *HCE)

virtual void	PrintAll ()

virtual bool	ProcessHits (G4Step step, G4TouchableHistory tHistory)

virtual bool	ProcessHits (G4GFlashSpot aSpot, G4TouchableHistory )

virtual	~CaloSD ()

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

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

Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point)

Local3DPoint	FinalStepPosition (G4Step *s, coordinates)

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

Local3DPoint	InitialStepPosition (G4Step *s, coordinates)

std::string	nameOfSD ()

void	NaNTrap (G4Step *step)

void	Register ()

	SensitiveDetector (std::string &iname, const DDCompactView &cpv, const 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 ()

Static Public Member Functions
static uint32_t	packIndex (int det, int lay, int x, int y)

static void	unpackIndex (const uint32_t &idx, int &det, int &lay, int &x, int &y)

Private Member Functions
std::vector< G4String >	getNames (DDFilteredView &)

bool	isItWireChamber (G4String)

Private Attributes
double	birk1

double	birk2

double	birk3

G4String	matName

const G4Material *	matScin

bool	useBirk

std::vector< G4String >	wcNames

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

Protected Member Functions inherited from CaloSD
G4bool	checkHit ()

virtual void	clearHits ()

CaloG4Hit *	createNewHit ()

virtual bool	filterHit (CaloG4Hit *, double)

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

virtual uint16_t	getDepth (G4Step *)

int	getNumberOfHits ()

double	getResponseWt (G4Track *)

virtual G4bool	getStepInfo (G4Step *aStep)

virtual int	getTrackID (G4Track *)

G4bool	hitExists ()

virtual void	initRun ()

void	resetForNewPrimary (const G4ThreeVector &, double)

G4ThreeVector	setToGlobal (const G4ThreeVector &, const G4VTouchable *)

G4ThreeVector	setToLocal (const G4ThreeVector &, const G4VTouchable *)

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

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

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

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

virtual void	update (const ::EndOfEvent *)

void	updateHit (CaloG4Hit *)

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

Protected Attributes inherited from CaloSD
int	checkHits

double	correctT

bool	corrTOFBeam

CaloG4Hit *	currentHit

CaloHitID	currentID

float	edepositEM

float	edepositHAD

double	eminHit

double	eminHitD

G4int	emPDG

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

G4int	epPDG

bool	forceSave

G4int	gammaPDG

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

const SimTrackManager *	m_trackManager

G4ThreeVector	posGlobal

G4StepPoint *	preStepPoint

CaloHitID	previousID

int	primIDSaved

bool	runInit

bool	suppressHeavy

G4Track *	theTrack

double	tmaxHit

bool	useMap

Detailed Description

Definition at line 20 of file HcalTB06BeamSD.h.

Constructor & Destructor Documentation

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

Definition at line 22 of file HcalTB06BeamSD.cc.

References DDFilteredView::addFilter(), birk1, birk2, birk3, equals, DDFilteredView::firstChild(), g, SensitiveDetector::getNames(), edm::ParameterSet::getParameter(), i, dqm-mbProfile::log, DDFilteredView::logicalPart(), DDLogicalPart::material(), matName, matScin, MeV, DDName::name(), DDBase< N, C >::name(), DDFilteredView::next(), DDSpecificsFilter::setCriteria(), 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";
   DDSpecificsFilter filter1;
   DDValue           ddv1(attribute,value,0);
   filter1.setCriteria(ddv1,DDCompOp::equals);
   DDFilteredView fv1(cpv);
   fv1.addFilter(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";
   DDSpecificsFilter filter2;
   DDValue           ddv2(attribute,name,0);
   filter2.setCriteria(ddv2,DDCompOp::equals);
   DDFilteredView fv2(cpv);
   fv2.addFilter(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.size() > 0) {
     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;
   matScin = G4Material::GetMaterial(matName);
 }

HcalTB06BeamSD::~HcalTB06BeamSD ( )

virtual

Definition at line 100 of file HcalTB06BeamSD.cc.

100 {}

Member Function Documentation

double HcalTB06BeamSD::getEnergyDeposit ( G4Step * aStep )

virtual

Reimplemented from CaloSD.

Definition at line 102 of file HcalTB06BeamSD.cc.

References birk1, birk2, birk3, CaloSD::getAttenuation(), LogDebug, matScin, useBirk, and puppiForMET_cff::weight.

                                                      {
 
   double destep = aStep->GetTotalEnergyDeposit();
   double weight = 1;
   if (useBirk && matScin == aStep->GetPreStepPoint()->GetMaterial()) {
     weight *= getAttenuation(aStep, birk1, birk2, birk3);
   }
   LogDebug("HcalTB06BeamSD") 
     << "HcalTB06BeamSD: Detector " 
     << aStep->GetPreStepPoint()->GetTouchable()->GetVolume()->GetName()
     << " weight " << weight;
   return weight*destep;
 }

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

private

Definition at line 168 of file HcalTB06BeamSD.cc.

References DDFilteredView::firstChild(), i, dqm-mbProfile::log, DDFilteredView::logicalPart(), DDName::name(), SensitiveDetector::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 ( G4String name )

private

Definition at line 184 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::packIndex	(	int	det,
		int	lay,
		int	x,
		int	y
	)

static

Definition at line 138 of file HcalTB06BeamSD.cc.

References TauDecayModes::dec, customizeTrackingMonitorSeedNumber::idx, ixx, iyy, LogDebug, x, and y.

Referenced by setDetUnitId().

                                                                  {
 
   int ix = 0, ixx = x;
   if (x < 0) { ix = 1; ixx =-x;}
   int iy = 0, iyy = y;
   if (y < 0) { iy = 1; iyy =-y;}
   uint32_t idx = (det&15)<<28;      //bits 28-31
   idx         += (lay&127)<<21;     //bits 21-27
   idx         += (iy&1)<<19;        //bit  19
   idx         += (iyy&511)<<10;     //bits 10-18
   idx         += (ix&1)<<9;         //bit   9
   idx         += (ixx&511);         //bits  0-8
 
   LogDebug("HcalTB06BeamSD") << "HcalTB06BeamSD: Detector " << det << " Layer "
                  << lay << " x " << x << " " << ix << " " << ixx 
                  << " y " << y << " " << iy << " " << iyy << " ID " 
                  << std::hex << idx << std::dec;
   return idx;
 }

uint32_t HcalTB06BeamSD::setDetUnitId ( G4Step * step )

virtual

Implements CaloSD.

Definition at line 116 of file HcalTB06BeamSD.cc.

References isItWireChamber(), SensitiveDetector::name, packIndex(), CaloSD::preStepPoint, CaloSD::setToLocal(), x, and y.

                                                     { 
 
   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 hitPoint    = preStepPoint->GetPosition();
     G4ThreeVector localPoint  = setToLocal(hitPoint, touch);
     x   = (int)(localPoint.x()/(0.2*mm));
     y   = (int)(localPoint.y()/(0.2*mm));
   }
 
   return packIndex (det, lay, x, y);
 }

void HcalTB06BeamSD::unpackIndex	(	const uint32_t &	idx,
		int &	det,
		int &	lay,
		int &	x,
		int &	y
	)

static

Definition at line 158 of file HcalTB06BeamSD.cc.

References customizeTrackingMonitorSeedNumber::idx, x, and y.

                                  {
 
   det  = (idx>>28)&15;
   lay  = (idx>>21)&127;
   y    = (idx>>10)&511; if (((idx>>19)&1) == 1) y = -y;
   x    = (idx)&511;     if (((idx>>9)&1)  == 1) x = -x;
 
 }