#include <FastTimerSD.h>

Inheritance diagram for FastTimerSD:

Public Member Functions
virtual void	clear ()

virtual void	DrawAll ()

virtual void	EndOfEvent (G4HCofThisEvent *eventHC)

	FastTimerSD (std::string, const DDCompactView &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

void	fillHits (edm::PSimHitContainer &, std::string use)

virtual double	getEnergyDeposit (G4Step *step)

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

virtual void	Initialize (G4HCofThisEvent *HCE)

virtual void	PrintAll ()

virtual bool	ProcessHits (G4Step , G4TouchableHistory )

virtual uint32_t	setDetUnitId (G4Step *)

virtual	~FastTimerSD ()

Public Member Functions inherited from SensitiveTkDetector
	SensitiveTkDetector (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)

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 BeginOfJob * >
	Observer ()

void	slotForUpdate (const BeginOfJob *iT)

virtual	~Observer ()

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 EndOfEvent * >
	Observer ()

void	slotForUpdate (const EndOfEvent *iT)

virtual	~Observer ()

Protected Attributes
float	edepositEM

float	edepositHAD

G4int	emPDG

G4int	epPDG

G4int	gammaPDG

Private Member Functions
virtual void	clearHits ()

void	CreateNewHit ()

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

void	GetStepInfo (G4Step *aStep)

G4bool	HitExists ()

void	ResetForNewPrimary ()

G4ThreeVector	SetToLocal (const G4ThreeVector &global)

G4ThreeVector	SetToLocalExit (const G4ThreeVector &globalPoint)

void	StoreHit (BscG4Hit *)

void	Summarize ()

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

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

void	UpdateHit ()

Private Attributes
BscG4Hit *	currentHit

G4VPhysicalVolume *	currentPV

float	edeposit

float	Eloss

G4ThreeVector	entrancePoint

int	eventno

G4ThreeVector	exitPoint

const FastTimeDDDConstants *	ftcons

G4int	hcID

G4ThreeVector	hitPoint

G4ThreeVector	hitPointExit

G4ThreeVector	hitPointLocal

G4ThreeVector	hitPointLocalExit

float	incidentEnergy

std::string	name

float	Pabs

int	ParentId

short	ParticleType

float	PhiAtEntry

G4StepPoint *	postStepPoint

G4StepPoint *	preStepPoint

uint32_t	previousUnitID

G4int	primaryID

G4int	primID

TrackingSlaveSD *	slave

G4ThreeVector	theEntryPoint

G4ThreeVector	theExitPoint

BscG4HitCollection *	theHC

const SimTrackManager *	theManager

float	ThetaAtEntry

G4Track *	theTrack

float	Tof

G4int	tsID

G4double	tSlice

G4int	tSliceID

int	type_

uint32_t	unitID

float	Vx

float	Vy

float	Vz

float	X

float	Y

float	Z

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

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

Detailed Description

Definition at line 34 of file FastTimerSD.h.

Constructor & Destructor Documentation

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

Definition at line 40 of file FastTimerSD.cc.

References DDFilteredView::addFilter(), SensitiveDetector::AssignSD(), ecaldqm::collectionName, gather_cfg::cout, dbl_to_int(), equals, alcazmumu_cfi::filter, DDFilteredView::firstChild(), getDDDArray(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), SensitiveDetectorCatalog::logicalNames(), DDFilteredView::mergedSpecifics(), SensitiveDetector::Register(), DDSpecificsFilter::setCriteria(), slave, AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, and type_.

                                              :
   SensitiveTkDetector(name, cpv, clg, p), ftcons(0), name(name),
   hcID(-1), theHC(0), theManager(manager), currentHit(0), theTrack(0), 
   currentPV(0), unitID(0),  previousUnitID(0), preStepPoint(0), 
   postStepPoint(0), eventno(0) {
     
   //Add FastTimer Sentitive Detector Name
   collectionName.insert(name);
     
     
   //Parameters
   edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("FastTimerSD");
   int verbn = m_p.getUntrackedParameter<int>("Verbosity");
     
   SetVerboseLevel(verbn);
 #ifdef EDM_ML_DEBUG
   std::cout << "*******************************************************\n"
         << "*                                                     *\n"
         << "* Constructing a FastTimerSD  with name " << name << "\n"
         << "*                                                     *\n"
         << "*******************************************************\n";
 #endif    
     
   slave  = new TrackingSlaveSD(name);
     
   //
   // attach detectors (LogicalVolumes)
   //
   std::vector<std::string> lvNames = clg.logicalNames(name);
 
   this->Register();
 
   for (std::vector<std::string>::iterator it=lvNames.begin();  
        it !=lvNames.end(); it++) {
     this->AssignSD(*it);
     edm::LogInfo("FastTimerSim") << "FastTimerSD : Assigns SD to LV " << (*it);
   }
     
   std::string attribute = "ReadOutName";
   DDSpecificsFilter filter;
   DDValue           ddv(attribute,name,0);
   filter.setCriteria(ddv,DDCompOp::equals);
   DDFilteredView fv(cpv);
   fv.addFilter(filter);
   fv.firstChild();
   DDsvalues_type sv(fv.mergedSpecifics());
   std::vector<int> temp = dbl_to_int(getDDDArray("Type",sv));
   type_  = temp[0];
 
   edm::LogInfo("FastTimerSim") << "FastTimerSD: Instantiation completed for "
                    << name << " of type " << type_;
 }

FastTimerSD::~FastTimerSD ( )

virtual

Definition at line 97 of file FastTimerSD.cc.

References slave.

                           { 
   if (slave)  delete slave; 
 }

Member Function Documentation

void FastTimerSD::clear ( void )

virtual

Definition at line 389 of file FastTimerSD.cc.

389 {}

void FastTimerSD::clearHits ( )

privatevirtual

Implements SensitiveDetector.

Definition at line 440 of file FastTimerSD.cc.

References TrackingSlaveSD::Initialize(), and slave.

Referenced by update().

                            {
   slave->Initialize();
 }

void FastTimerSD::CreateNewHit ( )

private

Definition at line 266 of file FastTimerSD.cc.

References gather_cfg::cout, currentHit, currentPV, Eloss, hitPoint, hitPointLocal, hitPointLocalExit, incidentEnergy, NULL, Pabs, ParentId, ParticleType, PhiAtEntry, primaryID, BscG4Hit::setEnergyLoss(), BscG4Hit::setEntry(), BscG4Hit::setEntryLocalP(), BscG4Hit::setExitLocalP(), BscG4Hit::setIncidentEnergy(), BscG4Hit::setPabs(), BscG4Hit::setParentId(), BscG4Hit::setParticleType(), BscG4Hit::setPhiAtEntry(), BscG4Hit::setThetaAtEntry(), BscG4Hit::setTimeSlice(), BscG4Hit::setTof(), BscG4Hit::setTrackID(), BscG4Hit::setUnitID(), BscG4Hit::setVx(), BscG4Hit::setVy(), BscG4Hit::setVz(), BscG4Hit::setX(), BscG4Hit::setY(), BscG4Hit::setZ(), StoreHit(), ThetaAtEntry, theTrack, Tof, tSlice, tSliceID, unitID, UpdateHit(), Vx, Vy, Vz, X, Y, and Z.

Referenced by ProcessHits().

                                {
 
 #ifdef EDM_ML_DEBUG
   std::cout << "FastTimerSD CreateNewHit for" << " PV "
         << currentPV->GetName() << " PVid = " << currentPV->GetCopyNo()
         << " Unit " << unitID << std::endl;
   std::cout << " primary " << primaryID << " time slice " << tSliceID 
         << " For Track  " << theTrack->GetTrackID() << " which is a "
         << theTrack->GetDefinition()->GetParticleName();
        
   if (theTrack->GetTrackID()==1) {
     std::cout << " of energy "     << theTrack->GetTotalEnergy();
   } else {
     std::cout << " daughter of part. " << theTrack->GetParentID();
   }
 
   std::cout << " and created by " ;
   if (theTrack->GetCreatorProcess()!=NULL)
     std::cout << theTrack->GetCreatorProcess()->GetProcessName() ;
   else 
     std::cout << "NO process";
   std::cout << std::endl;
 #endif          
     
   currentHit = new BscG4Hit;
   currentHit->setTrackID(primaryID);
   currentHit->setTimeSlice(tSlice);
   currentHit->setUnitID(unitID);
   currentHit->setIncidentEnergy(incidentEnergy);
 
   currentHit->setPabs(Pabs);
   currentHit->setTof(Tof);
   currentHit->setEnergyLoss(Eloss);
   currentHit->setParticleType(ParticleType);
   currentHit->setThetaAtEntry(ThetaAtEntry);
   currentHit->setPhiAtEntry(PhiAtEntry);
 
   currentHit->setEntry(hitPoint);
 
   currentHit->setEntryLocalP(hitPointLocal);
   currentHit->setExitLocalP(hitPointLocalExit);
 
   currentHit->setParentId(ParentId);
   currentHit->setVx(Vx);
   currentHit->setVy(Vy);
   currentHit->setVz(Vz);
 
   currentHit->setX(X);
   currentHit->setY(Y);
   currentHit->setZ(Z);
 
   UpdateHit();
   
   StoreHit(currentHit);
 }    

void FastTimerSD::DrawAll ( )

virtual

Definition at line 391 of file FastTimerSD.cc.

391 {}

void FastTimerSD::EndOfEvent ( G4HCofThisEvent * eventHC )

virtual

Reimplemented from SensitiveDetector.

Definition at line 359 of file FastTimerSD.cc.

References gather_cfg::cout, TauDecayModes::dec, BscG4Hit::getEnergyLoss(), BscG4Hit::getEntry(), BscG4Hit::getPabs(), BscG4Hit::getParticleType(), BscG4Hit::getPhiAtEntry(), BscG4Hit::getThetaAtEntry(), BscG4Hit::getTof(), BscG4Hit::getTrackID(), BscG4Hit::getUnitID(), j, TrackingSlaveSD::processHits(), slave, Summarize(), and theHC.

                                               {
 
   // here we loop over transient hits and make them persistent
   for (int j=0; j<theHC->entries(); j++) {
     BscG4Hit* aHit = (*theHC)[j];
 #ifdef EDM_ML_DEBUG
     std::cout << "hit number " << j << " unit ID = " << std::hex 
           << aHit->getUnitID() << std::dec << " entry z "
           << aHit->getEntry().z() << " entry theta "
           << aHit->getThetaAtEntry() << std::endl;
 #endif
     Local3DPoint locExitPoint(0,0,0);
     Local3DPoint locEntryPoint(aHit->getEntry().x(),
                    aHit->getEntry().y(),
                    aHit->getEntry().z());
     slave->processHits(PSimHit(locEntryPoint,locExitPoint,
                    aHit->getPabs(),
                    aHit->getTof(),
                    aHit->getEnergyLoss(),
                    aHit->getParticleType(),
                    aHit->getUnitID(),
                    aHit->getTrackID(),
                    aHit->getThetaAtEntry(),
                    aHit->getPhiAtEntry()));
   }
   Summarize();
 }

void FastTimerSD::fillHits	(	edm::PSimHitContainer &	c,
		std::string	use
	)

virtual

Implements SensitiveTkDetector.

Definition at line 400 of file FastTimerSD.cc.

References TrackingSlaveSD::hits(), gen::n, TrackingSlaveSD::name(), and slave.

                                                               {
   if (slave->name() == n) c=slave->hits();
 }

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

private

Definition at line 450 of file FastTimerSD.cc.

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

Referenced by FastTimerSD().

                                                     {
 
   DDValue value(str);
   if (DDfetch(&sv,value)) {
     const std::vector<double> & fvec = value.doubles();
     int nval = fvec.size();
     if (nval < 1) {
       edm::LogError("FastTimerSim") << "FastTimerSD : # of " << str
                     << " bins " << nval << " < 1 ==> illegal";
       throw cms::Exception("DDException") << "FastTimerSD: cannot get array " << str;
     }
     return fvec;
   } else {
     edm::LogError("FastTimerSim") << "FastTimerSD: cannot get array " << str;
     throw cms::Exception("DDException") << "FastTimerSD: cannot get array " << str;
   }
 }

double FastTimerSD::getEnergyDeposit ( G4Step * step )

virtual

Definition at line 101 of file FastTimerSD.cc.

Referenced by GetStepInfo().

                                                   {
   return aStep->GetTotalEnergyDeposit();
 }

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

virtual

Reimplemented from SensitiveDetector.

Definition at line 444 of file FastTimerSD.cc.

References TrackingSlaveSD::name(), slave, and groupFilesInBlocks::temp.

                                           {
   std::vector<std::string> temp;
   temp.push_back(slave->name());
   return temp;
 }

void FastTimerSD::GetStepInfo ( G4Step * aStep )

private

Definition at line 137 of file FastTimerSD.cc.

References gather_cfg::cout, currentPV, TauDecayModes::dec, edeposit, edepositEM, edepositHAD, Eloss, emPDG, epPDG, gammaPDG, getEnergyDeposit(), GeV, hitPoint, hitPointExit, hitPointLocal, hitPointLocalExit, Pabs, ParentId, ParticleType, PhiAtEntry, postStepPoint, preStepPoint, primaryID, setDetUnitId(), ThetaAtEntry, theTrack, Tof, tSlice, tSliceID, unitID, Vx, Vy, Vz, X, Y, and Z.

Referenced by ProcessHits().

                                            {
   
   preStepPoint = aStep->GetPreStepPoint(); 
   postStepPoint= aStep->GetPostStepPoint(); 
   theTrack     = aStep->GetTrack();   
   hitPoint     = preStepPoint->GetPosition();   
   currentPV    = preStepPoint->GetPhysicalVolume();
   hitPointExit = postStepPoint->GetPosition();  
 
   hitPointLocal = preStepPoint->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
   hitPointLocalExit = preStepPoint->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPointExit);
 
 
   G4int particleCode = theTrack->GetDefinition()->GetPDGEncoding();
 #ifdef EDM_ML_DEBUG
   std::cout << "FastTimerSD :particleType =  " 
         << theTrack->GetDefinition()->GetParticleName() << std::endl;
 #endif
   if (particleCode == emPDG ||
       particleCode == epPDG ||
       particleCode == gammaPDG ) {
     edepositEM  = getEnergyDeposit(aStep); edepositHAD = 0.;
   } else {
     edepositEM  = 0.; edepositHAD = getEnergyDeposit(aStep);
   }
   edeposit = aStep->GetTotalEnergyDeposit();
   tSlice    = (100*postStepPoint->GetGlobalTime() )/CLHEP::nanosecond;
   tSliceID  = (int) tSlice;
   unitID    = setDetUnitId(aStep);
 #ifdef EDM_ML_DEBUG
   std::cout << "FastTimerSD:unitID = " << std::hex << unitID << std::dec<<"\n";
 #endif
   primaryID    = theTrack->GetTrackID();
   //  Position     = hitPoint;
   Pabs         = aStep->GetPreStepPoint()->GetMomentum().mag()/CLHEP::GeV;
   Tof          = aStep->GetPostStepPoint()->GetGlobalTime()/CLHEP::nanosecond;
   Eloss        = aStep->GetTotalEnergyDeposit()/CLHEP::GeV;
   ParticleType = theTrack->GetDefinition()->GetPDGEncoding();      
   ThetaAtEntry = aStep->GetPreStepPoint()->GetPosition().theta()/CLHEP::deg;
   PhiAtEntry   = aStep->GetPreStepPoint()->GetPosition().phi()/CLHEP::deg;
 
   ParentId = theTrack->GetParentID();
   Vx = theTrack->GetVertexPosition().x();
   Vy = theTrack->GetVertexPosition().y();
   Vz = theTrack->GetVertexPosition().z();
   X  = hitPoint.x();
   Y  = hitPoint.y();
   Z  = hitPoint.z();
 }

G4bool FastTimerSD::HitExists ( )

private

Definition at line 208 of file FastTimerSD.cc.

References currentHit, newFWLiteAna::found, BscG4Hit::getTimeSliceID(), BscG4Hit::getTrackID(), BscG4Hit::getUnitID(), j, previousUnitID, primaryID, primID, ResetForNewPrimary(), theHC, tsID, tSliceID, unitID, and UpdateHit().

Referenced by ProcessHits().

                               {
   if (primaryID<1) {
     edm::LogWarning("FastTimerSim") << "***** FastTimerSD error: primaryID = " 
                     << primaryID
                     << " maybe detector name changed";
   }
 
   // Update if in the same detector, time-slice and for same track   
   if (tSliceID == tsID && unitID==previousUnitID) {
     UpdateHit();
     return true;
   }
   // Reset entry point for new primary
   if (primaryID != primID)
     ResetForNewPrimary();
    
   //look in the HitContainer whether a hit with the same primID, unitID,
   //tSliceID already exists:
    
   G4bool found = false;
 
   for (int j=0; j<theHC->entries()&&!found; j++) {
     BscG4Hit* aPreviousHit = (*theHC)[j];
     if (aPreviousHit->getTrackID()     == primaryID &&
     aPreviousHit->getTimeSliceID() == tSliceID  &&
     aPreviousHit->getUnitID()      == unitID       ) {
       currentHit = aPreviousHit;
       found      = true;
     }
   }          
 
   if (found) {
     UpdateHit();
     return true;
   } else {
     return false;
   }    
 }

void FastTimerSD::Initialize ( G4HCofThisEvent * HCE )

virtual

Reimplemented from SensitiveDetector.

Definition at line 105 of file FastTimerSD.cc.

References ecaldqm::collectionName, gather_cfg::cout, hcID, name, primID, theHC, and tsID.

                                                   { 
 #ifdef EDM_ML_DEBUG
   std::cout << "FastTimerSD : Initialize called for " << name << std::endl;
 #endif
 
   theHC = new BscG4HitCollection(name, collectionName[0]);
   if (hcID<0) 
     hcID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(hcID, theHC);
 
   tsID   = -2;
   primID = -2;
 }

void FastTimerSD::PrintAll ( )

virtual

Definition at line 393 of file FastTimerSD.cc.

References gather_cfg::cout, and theHC.

                            {
 #ifdef EDM_ML_DEBUG
   std::cout << "FastTimerSD: Collection " << theHC->GetName() << std::endl;
 #endif
   theHC->PrintAllHits();
 } 

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

virtual

Implements SensitiveDetector.

Definition at line 120 of file FastTimerSD.cc.

References gather_cfg::cout, CreateNewHit(), edeposit, GetStepInfo(), HitExists(), NULL, and theHC.

                                                                    {
 
   if (aStep == NULL) {
     return true;
   } else {
     GetStepInfo(aStep);
 #ifdef EDM_ML_DEBUG
     std::cout << "FastTimerSD :  number of hits = " << theHC->entries() <<"\n";
 #endif
     if (HitExists() == false && edeposit>0. ){ 
       CreateNewHit();
       return true;
     }
   }
   return true;
 } 

void FastTimerSD::ResetForNewPrimary ( )

private

Definition at line 248 of file FastTimerSD.cc.

References entrancePoint, exitPoint, hitPoint, hitPointExit, incidentEnergy, preStepPoint, SetToLocal(), and SetToLocalExit().

Referenced by HitExists().

                                      {
   entrancePoint  = SetToLocal(hitPoint);
   exitPoint      = SetToLocalExit(hitPointExit);
   incidentEnergy = preStepPoint->GetKineticEnergy();
 }

uint32_t FastTimerSD::setDetUnitId ( G4Step * aStep )

virtual

Implements SensitiveDetector.

Definition at line 187 of file FastTimerSD.cc.

References funct::abs(), gather_cfg::cout, TauDecayModes::dec, ftcons, FastTimeDDDConstants::getEtaPhi(), FastTimeDDDConstants::getZPhi(), DetId::rawId(), and type_.

Referenced by GetStepInfo().

                                                  { 
 
   //Find the depth segment
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   G4ThreeVector global = aStep->GetPreStepPoint()->GetPosition();
   G4ThreeVector local  = touch->GetHistory()->GetTopTransform().TransformPoint(global);
   int        iz = (global.z() > 0) ? 1 : -1;
   std::pair<int,int> izphi = ((ftcons) ? ((type_ == 1) ? 
                       (ftcons->getZPhi(std::abs(local.z()),local.phi())) : 
                       (ftcons->getEtaPhi(local.perp(),local.phi()))) :
                   (std::pair<int,int>(0,0)));
   uint32_t id = FastTimeDetId(type_,izphi.first,izphi.second,iz).rawId();
 #ifdef EDM_ML_DEBUG
   std::cout << "Volume " << touch->GetVolume(0)->GetName() << ":" << global.z()
         << " Iz(eta)phi " << izphi.first << ":"  << izphi.second << ":" 
         << iz  << " id " << std::hex << id << std::dec << std::endl;
 #endif
   return id;
 }

G4ThreeVector FastTimerSD::SetToLocal ( const G4ThreeVector & global )

private

Definition at line 343 of file FastTimerSD.cc.

References preStepPoint, and theEntryPoint.

Referenced by ResetForNewPrimary().

                                                                 {
 
   const G4VTouchable* touch= preStepPoint->GetTouchable();
   theEntryPoint = touch->GetHistory()->GetTopTransform().TransformPoint(global);
   return theEntryPoint;  
 }

G4ThreeVector FastTimerSD::SetToLocalExit ( const G4ThreeVector & globalPoint )

private

Definition at line 351 of file FastTimerSD.cc.

References postStepPoint, and theExitPoint.

Referenced by ResetForNewPrimary().

                                                                          {
 
   const G4VTouchable* touch= postStepPoint->GetTouchable();
   theExitPoint = touch->GetHistory()->GetTopTransform().TransformPoint(globalPoint);
   return theExitPoint;  
 }

void FastTimerSD::StoreHit ( BscG4Hit * hit )

private

Definition at line 255 of file FastTimerSD.cc.

References primID, and theHC.

Referenced by CreateNewHit().

                                        {
 
   if (primID<0) return;
   if (hit == 0) {
     edm::LogWarning("FastTimerSim") << "FastTimerSD: hit to be stored is NULL !!";
   } else {
     theHC->insert( hit );
   }
 }

void FastTimerSD::Summarize ( )

private

Definition at line 387 of file FastTimerSD.cc.

Referenced by EndOfEvent().

387 {}

void FastTimerSD::update ( const BeginOfJob * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 404 of file FastTimerSD.cc.

References gather_cfg::cout, Exception, ftcons, and edm::EventSetup::get().

                                                {
 
   const edm::EventSetup* es = (*job)();
   edm::ESHandle<FastTimeDDDConstants> fdc;
   es->get<IdealGeometryRecord>().get(fdc);
   if (fdc.isValid()) {
     ftcons = &(*fdc);
   } else {
     edm::LogError("FastTimerSim") << "FastTimerSD : Cannot find FastTimeDDDConstants";
     throw cms::Exception("Unknown", "FastTimerSD") << "Cannot find FastTimeDDDConstants\n";
   }
 #ifdef EDM_ML_DEBUG
   std::cout << "FastTimerSD::Initialized with FastTimeDDDConstants\n";
 #endif
 }

void FastTimerSD::update ( const BeginOfRun * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 428 of file FastTimerSD.cc.

References emPDG, epPDG, and gammaPDG.

                                            {
 
   G4ParticleTable * theParticleTable = G4ParticleTable::GetParticleTable();
   G4String particleName;
   emPDG = theParticleTable->FindParticle(particleName="e-")->GetPDGEncoding();
   epPDG = theParticleTable->FindParticle(particleName="e+")->GetPDGEncoding();
   gammaPDG = theParticleTable->FindParticle(particleName="gamma")->GetPDGEncoding();
 
 } 

void FastTimerSD::update ( const BeginOfEvent * )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 420 of file FastTimerSD.cc.

References clearHits(), gather_cfg::cout, and eventno.

                                                 {
 #ifdef EDM_ML_DEBUG
   std::cout << "Dispatched BeginOfEvent for " << GetName() << " !\n" ;
 #endif
    clearHits();
    eventno = (*i)()->GetEventID();
 }

void FastTimerSD::update ( const ::EndOfEvent * )

privatevirtual

Definition at line 438 of file FastTimerSD.cc.

438 {}

void FastTimerSD::UpdateHit ( )

private

Definition at line 323 of file FastTimerSD.cc.

References BscG4Hit::addEnergyDeposit(), gather_cfg::cout, currentHit, edepositEM, edepositHAD, Eloss, postStepPoint, previousUnitID, primaryID, primID, BscG4Hit::setEnergyLoss(), tsID, tSliceID, and unitID.

Referenced by CreateNewHit(), and HitExists().

                             {
 
   if (Eloss > 0.) {
     currentHit->addEnergyDeposit(edepositEM,edepositHAD);
 
 #ifdef EDM_ML_DEBUG
     std::cout << "updateHit: add eloss " << Eloss <<std::endl;
     std::cout << "CurrentHit="<< currentHit<< ", PostStepPoint = "
           << postStepPoint->GetPosition() << std::endl;
 #endif
     currentHit->setEnergyLoss(Eloss);
   }  
 
   // buffer for next steps:
   tsID           = tSliceID;
   primID         = primaryID;
   previousUnitID = unitID;
 }