#include <FiberSD.h>

Inheritance diagram for FiberSD:

Public Member Functions
virtual void	clear ()

virtual void	DrawAll ()

virtual void	EndOfEvent (G4HCofThisEvent *HCE)

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

virtual void	Initialize (G4HCofThisEvent *HCE)

virtual void	PrintAll ()

virtual G4bool	ProcessHits (G4Step aStep, G4TouchableHistory ROhist)

virtual	~FiberSD ()

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 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 Member Functions
virtual void	clearHits ()

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

virtual uint32_t	setDetUnitId (G4Step *)

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

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 Attributes
const SimTrackManager *	m_trackManager

FiberG4HitsCollection *	theHC

G4int	theHCID

std::string	theName

HFShower *	theShower

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

Detailed Description

Definition at line 26 of file FiberSD.h.

Constructor & Destructor Documentation

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

Definition at line 21 of file FiberSD.cc.

References SensitiveDetector::AssignSD(), ecaldqm::collectionName, g4SimHits_cfi::HFShower, LogDebug, SensitiveDetectorCatalog::logicalNames(), SensitiveDetector::Register(), and theShower.

                                          :
   SensitiveCaloDetector(name, cpv, clg, p), theName(name),
   m_trackManager(manager), theHCID(-1), theHC(0) {
 
   collectionName.insert(name);
   LogDebug("FiberSim") << "***************************************************"
                << "\n"
                << "*                                                 *"
                << "\n"
                << "* Constructing a FiberSD  with name " << GetName()
                << "\n"
                << "*                                                 *"
                << "\n"
                << "***************************************************";
   theShower = new HFShower(name, cpv, p, 1);
 
   //
   // Now attach the right detectors (LogicalVolumes) to me
   //
   const std::vector<std::string>& lvNames = clg.logicalNames(name);
   this->Register();
   for (std::vector<std::string>::const_iterator it=lvNames.begin();
        it !=lvNames.end(); it++){
     this->AssignSD(*it);
     LogDebug("FiberSim") << "FiberSD : Assigns SD to LV " << (*it);
   }
 }

FiberSD::~FiberSD ( )

virtual

Definition at line 51 of file FiberSD.cc.

References theHC, and theShower.

                   {
  
   if (theShower) delete theShower;
   if (theHC)     delete theHC;
 }

Member Function Documentation

void FiberSD::clear ( void )

virtual

Definition at line 125 of file FiberSD.cc.

Referenced by EndOfEvent().

125 {}

void FiberSD::clearHits ( )

protectedvirtual

Implements SensitiveDetector.

Definition at line 152 of file FiberSD.cc.

152 {}

void FiberSD::DrawAll ( )

virtual

Definition at line 127 of file FiberSD.cc.

127 {}

void FiberSD::EndOfEvent ( G4HCofThisEvent * HCE )

virtual

Reimplemented from SensitiveDetector.

Definition at line 118 of file FiberSD.cc.

References clear(), gather_cfg::cout, LogDebug, and theHC.

                                               {
  
   LogDebug("FiberSim") << "FiberSD: Sees" << theHC->entries() << " hits";
   clear();
   std::cout<<"theHC entries = "<<theHC->entries()<<std::endl;
 }

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

protectedvirtual

Implements SensitiveCaloDetector.

Definition at line 162 of file FiberSD.cc.

162 {}

void FiberSD::Initialize ( G4HCofThisEvent * HCE )

virtual

Reimplemented from SensitiveDetector.

Definition at line 57 of file FiberSD.cc.

References ecaldqm::collectionName, LogDebug, theHC, and theHCID.

                                               {
 
   LogDebug("FiberSim") << "FiberSD : Initialize called for " << GetName();
   theHC = new FiberG4HitsCollection(GetName(), collectionName[0]);
   if (theHCID<0)
     theHCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(theHCID, theHC);
   
 }

void FiberSD::PrintAll ( )

virtual

Definition at line 129 of file FiberSD.cc.

129 {}

G4bool FiberSD::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *	ROhist
	)

virtual

Implements SensitiveDetector.

Definition at line 67 of file FiberSD.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, HFShower::getHits(), hfClusterShapes_cfi::hits, mps_fire::i, LogDebug, position, setDetUnitId(), theHC, and theShower.

                                                                {
 
   //std::vector<HFShower::Hit> hits = theShower->getHits(aStep);
   double zoffset = 1000;
   std::vector<HFShower::Hit> hits = theShower->getHits(aStep,true,zoffset);
 
   
   if (hits.size() > 0) {
     std::vector<HFShowerPhoton> thePE;
     for (unsigned int i=0; i<hits.size(); i++) {
       //std::cout<<"hit position z "<<hits[i].position.z()<<std::endl;
       HFShowerPhoton pe = HFShowerPhoton(hits[i].position.x(), 
                      hits[i].position.y(), 
                      hits[i].position.z(), 
                      hits[i].wavelength, hits[i].time);
       thePE.push_back(pe);
     }
     int trackID = aStep->GetTrack()->GetTrackID();
     G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
     const G4VTouchable* touch = preStepPoint->GetTouchable();
     G4LogicalVolume* lv = touch->GetVolume(0)->GetLogicalVolume();
     int depth = (touch->GetReplicaNumber(0))%10;
     int detID = setDetUnitId(aStep);
     math::XYZPoint theHitPos(preStepPoint->GetPosition().x(),
                  preStepPoint->GetPosition().y(),
                  preStepPoint->GetPosition().z());
     //std::cout<<"presteppoint position z "<<preStepPoint->GetPosition().z()<<std::endl;
     
     FiberG4Hit *aHit = new FiberG4Hit(lv, detID, depth, trackID);
     std::cout<<"hit size "<<hits.size()<<"  npe"<<aHit->npe()<<std::endl;
     std::cout<<"pre hit position "<<aHit->hitPos()<<std::endl;
     aHit->setNpe(hits.size());
     aHit->setPos(theHitPos);
     aHit->setTime(preStepPoint->GetGlobalTime());
     aHit->setPhoton(thePE);
     std::cout<<"ShowerPhoton position "<<thePE[0].x()<<" "<<thePE[0].y()<<" "<<thePE[0].z()<<std::endl;
 
     LogDebug("FiberSim") << "FiberSD: Hit created at " << lv->GetName()
              << " DetID: " << aHit->towerId() << " Depth: " 
              << aHit->depth() << " Track ID: " << aHit->trackId() 
              << " Nb. of Cerenkov Photons: " << aHit->npe()
              << " Time: " << aHit->time() << " at " 
              << aHit->hitPos();
     for (unsigned int i=0; i<thePE.size(); i++) 
       LogDebug("FiberSim") << "FiberSD: PE[" << i << "] " << thePE[i];
 
     theHC->insert(aHit);
   }
   return true;
 }

uint32_t FiberSD::setDetUnitId ( G4Step * aStep )

protectedvirtual

Implements SensitiveDetector.

Definition at line 154 of file FiberSD.cc.

Referenced by ProcessHits().

                                             {
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   int fibre = (touch->GetReplicaNumber(1))%10;
   int cell  = (touch->GetReplicaNumber(2));
   int tower = (touch->GetReplicaNumber(3));
   return ((tower*1000+cell)*10+fibre);
 }

void FiberSD::update ( const BeginOfJob * )

protectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 131 of file FiberSD.cc.

References Exception, edm::EventSetup::get(), HFShower::initRun(), and theShower.

                                            {
 
   const edm::EventSetup* es = (*job)();
   edm::ESHandle<HcalDDDSimConstants>    hdc;
   es->get<HcalSimNumberingRecord>().get(hdc);
   if (hdc.isValid()) {
     HcalDDDSimConstants *hcalConstants = (HcalDDDSimConstants*)(&(*hdc));
     theShower->initRun(0, hcalConstants);
   } else {
     edm::LogError("HcalSim") << "HCalSD : Cannot find HcalDDDSimConstant";
     throw cms::Exception("Unknown", "HCalSD") << "Cannot find HcalDDDSimConstant" << "\n";
   }
 
 }

void FiberSD::update ( const BeginOfRun * )

protectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 146 of file FiberSD.cc.

146 {}

void FiberSD::update ( const BeginOfEvent * )

protectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 148 of file FiberSD.cc.

148 {}

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

protectedvirtual

Definition at line 150 of file FiberSD.cc.

150 {}

Member Data Documentation

const SimTrackManager* FiberSD::m_trackManager

private

Definition at line 59 of file FiberSD.h.

FiberG4HitsCollection* FiberSD::theHC

private

Definition at line 63 of file FiberSD.h.

Referenced by EndOfEvent(), Initialize(), ProcessHits(), and ~FiberSD().

G4int FiberSD::theHCID

private

Definition at line 62 of file FiberSD.h.

Referenced by Initialize().

std::string FiberSD::theName

private

Definition at line 58 of file FiberSD.h.

HFShower* FiberSD::theShower

private

Definition at line 60 of file FiberSD.h.

Referenced by FiberSD(), ProcessHits(), update(), and ~FiberSD().

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation