#include <LaserAlignmentSimulation.h>

Inheritance diagram for LaserAlignmentSimulation:

Public Member Functions
	LaserAlignmentSimulation (edm::ParameterSet const &theConf)
	constructor More...

void	update (const BeginOfRun *myRun)
	observer for BeginOfRun More...

void	update (const BeginOfEvent *myEvent)
	observer for BeginOfEvent More...

void	update (const G4Step *myStep)
	observer for G4Step More...

void	update (const EndOfEvent *myEvent)
	observer for EndOfEvent More...

void	update (const EndOfRun *myRun)
	observer for EndOfRun More...

void	update (const BeginOfTrack *myTrack)
	observer for BeginOfTrack More...

void	update (const EndOfTrack *myTrack)
	observer for EndOfTrack More...

virtual	~LaserAlignmentSimulation ()
	destructor More...

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

void	slotForUpdate (const G4Step * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfEvent * >
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfRun * >
	Observer ()

void	slotForUpdate (const EndOfRun * 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 SimWatcher
	SimWatcher ()

virtual	~SimWatcher ()

Private Attributes
int	theBarrelHits

int	theDebugLevel

int	theEndcapHits

double	theEnergyLossScalingFactor

MaterialProperties *	theMaterialProperties

int	theMPDebug

edm::ParameterSet	theParameterSet

LaserPrimaryGeneratorAction *	thePrimaryGenerator

double	theSiAbsLengthScale

LaserSteppingAction *	theSteppingAction

G4Timer *	theTimer

Additional Inherited Members

Detailed Description

SimWatcher for the simulation of the Laser Alignment System of the CMS Tracker

Date: 2007/03/20 12:01:00

Revision: 1.3

Author: Maarten Thomas

Definition at line 42 of file LaserAlignmentSimulation.h.

Constructor & Destructor Documentation

LaserAlignmentSimulation::LaserAlignmentSimulation ( edm::ParameterSet const & theConf )

explicit

constructor

Definition at line 27 of file LaserAlignmentSimulation.cc.

References theDebugLevel, theEnergyLossScalingFactor, theMPDebug, theSiAbsLengthScale, and theTimer.

   : theDebugLevel(theConf.getUntrackedParameter<int>("DebugLevel",0)),
     theEnergyLossScalingFactor(theConf.getUntrackedParameter<double>("EnergyLossScalingFactor",1.0)),
     theMPDebug(theConf.getUntrackedParameter<int>("MaterialPropertiesDebugLevel",0)),
     theSiAbsLengthScale(theConf.getUntrackedParameter<double>("SiAbsorptionLengthScalingFactor",1.0)),
     theTimer(), 
     theMaterialProperties(),
     thePrimaryGenerator(), theSteppingAction(),
     theBarrelHits(0), theEndcapHits(0),
         theParameterSet(theConf)
 {
 
   // make some noise
   edm::LogInfo("SimLaserAlignmentSimulation") << " *****     AC1CMS: Configuration from ParameterSet      ***** " 
                     << "\n  AC1CMS: theDebugLevel               = " << theDebugLevel 
                     << "\n  AC1CMS: theEnergyLossScalingFactor  = " << theEnergyLossScalingFactor 
                     << "\n  AC1CMS: theMPDebugLevel             = " << theMPDebug
                     << "\n  AC1CMS: theSiAbsLengthScalingFactor = " << theSiAbsLengthScale;
 
   // declare timer
   theTimer = new G4Timer;
 }

LaserAlignmentSimulation::~LaserAlignmentSimulation ( )

virtual

destructor

Definition at line 50 of file LaserAlignmentSimulation.cc.

References theMaterialProperties, thePrimaryGenerator, theSteppingAction, and theTimer.

 {
   if ( theMaterialProperties != 0 )        { delete theMaterialProperties; }
   if ( theSteppingAction != 0 )            { delete theSteppingAction; }
   if ( thePrimaryGenerator != 0 )          { delete thePrimaryGenerator; }
   if ( theTimer != 0 )                     { delete theTimer; }
 }

Member Function Documentation

void LaserAlignmentSimulation::update ( const BeginOfRun * myRun )

virtual

observer for BeginOfRun

Implements Observer< const BeginOfRun * >.

Definition at line 58 of file LaserAlignmentSimulation.cc.

References LogDebug, theDebugLevel, theMaterialProperties, theMPDebug, theParameterSet, thePrimaryGenerator, theSiAbsLengthScale, theSteppingAction, and theTimer.

 {
   LogDebug("SimLaserAlignmentSimulation") << "<LaserAlignmentSimulation::update(const BeginOfRun * myRun)>"
                 << "\n *****     AC1CMS: Start of Run: " << (*myRun)()->GetRunID() << "     ***** ";
 
   // start timer
   theTimer->Start();
 
 
   // the PrimaryGeneratorAction: defines the used particlegun for the Laser events
   thePrimaryGenerator = new LaserPrimaryGeneratorAction(theParameterSet);
 
   // the UserSteppingAction: at the moment this prints only some information
   theSteppingAction = new LaserSteppingAction(theParameterSet);
 
   // construct your own material properties for setting refractionindex and so on
   theMaterialProperties = new MaterialProperties(theMPDebug, theSiAbsLengthScale);
 
   // list the tree of sensitive detectors
   if (theDebugLevel >= 1)
     {
       G4SDManager * theSDManager = G4SDManager::GetSDMpointer();
       theSDManager->ListTree();
     }
 }

void LaserAlignmentSimulation::update ( const BeginOfEvent * myEvent )

virtual

observer for BeginOfEvent

Implements Observer< const BeginOfEvent * >.

Definition at line 84 of file LaserAlignmentSimulation.cc.

References LaserPrimaryGeneratorAction::GeneratePrimaries(), LogDebug, theBarrelHits, theEndcapHits, and thePrimaryGenerator.

 {
   LogDebug("SimLaserAlignmentSimulation") << "<LaserAlignmentSimulation::update(const BeginOfEvent * myEvent)>"
                 << "\n AC1CMS: Event number = " << (*myEvent)()->GetEventID();
 
   // some statistics for this event
   theBarrelHits = 0;
   theEndcapHits = 0;
 
   // generate the Primaries
   thePrimaryGenerator->GeneratePrimaries((G4Event*)(*myEvent)());
 }

void LaserAlignmentSimulation::update ( const G4Step * myStep )

virtual

observer for G4Step

Implements Observer< const G4Step * >.

Definition at line 101 of file LaserAlignmentSimulation.cc.

References LogDebug, TkAccumulatingSensitiveDetector::ProcessHits(), theBarrelHits, theEndcapHits, theSteppingAction, and LaserSteppingAction::UserSteppingAction().

 {
   LogDebug("SimLaserAlignmentSimulationStepping") << "<LaserAlignmentSimulation::update(const G4Step * myStep)>";
 
   G4Step * theStep = const_cast<G4Step*>(myStep);
 
   // do the LaserSteppingAction
   theSteppingAction->UserSteppingAction(theStep);
 
   // Trigger sensitive detector manually since photon is absorbed
   if ( ( theStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName()== "OpAbsorption" ) )
     {
       LogDebug("SimLaserAlignmentSimulationStepping") << "<LaserAlignmentSimulation::update(const G4Step*)>: Photon was absorbed! ";
       
       
       if ( theStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetSensitiveDetector() )
     {
       LogDebug("SimLaserAlignmentSimulationStepping") << " AC1CMS: Setting the EnergyLoss to " << theStep->GetTotalEnergyDeposit() 
                         << "\n AC1CMS: The z position is " << theStep->GetPreStepPoint()->GetPosition().z()
                         << "\n AC1CMS: the Sensitive Detector: " 
                         << theStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetSensitiveDetector()->GetName()
                         << "\n AC1CMS: the Material: " << theStep->GetPreStepPoint()->GetMaterial()->GetName()
                         << "\n AC1CMS: the Logical Volume: " 
                         << theStep->GetPostStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName();
 
       if (theStep->GetTotalEnergyDeposit() > 0.0)
         {
           // process a hit
           TkAccumulatingSensitiveDetector * theSD = (TkAccumulatingSensitiveDetector*)
         (theStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetSensitiveDetector());
           
           theSD->ProcessHits(theStep, ((G4TouchableHistory *)(theStep->GetPreStepPoint()->GetTouchable())));
 
 
           // some statistics for this event
           if ( ( theStep->GetPostStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName() == "TECModule3RphiActive" ) || 
            ( theStep->GetPostStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName() == "TECModule5RphiActive" ) )
         {
           theEndcapHits++;
         }
           else if ( ( theStep->GetPostStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName() == "TOBActiveSter0" ) ||
             ( theStep->GetPostStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName() == "TOBActiveRphi0" ) ||
             ( theStep->GetPostStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName() == "TIBActiveRphi2" ) )
         {
           theBarrelHits++;
         }
         }
     }
       else 
     {
       LogDebug("SimLaserAlignmentSimulationStepping") << " AC1CMS: No SensitiveDetector available for this Step ... No Hit created :-( "
                         << "\n AC1CMS: The Material was: " << theStep->GetPreStepPoint()->GetMaterial()->GetName(); 
     }
     }
 }

void LaserAlignmentSimulation::update ( const EndOfEvent * myEvent )

virtual

observer for EndOfEvent

Implements Observer< const EndOfEvent * >.

Definition at line 161 of file LaserAlignmentSimulation.cc.

References LogDebug, theBarrelHits, and theEndcapHits.

 {
   LogDebug("SimLaserAlignmentSimulation") << "<LaserAlignmentSimulation::update(const EndOfEvent * myEvent)>"
                 << "\n AC1CMS: End of Event " << (*myEvent)()->GetEventID();
 
   // some statistics for this event
   edm::LogInfo("SimLaserAlignmentSimulation") << " *** Number of Hits: " << theBarrelHits << " / " << theEndcapHits
                     << " (Barrel / Endcaps) *** ";
 }

void LaserAlignmentSimulation::update ( const EndOfRun * myRun )

virtual

observer for EndOfRun

Implements Observer< const EndOfRun * >.

Definition at line 171 of file LaserAlignmentSimulation.cc.

References DEFINE_SIMWATCHER, LogDebug, and theTimer.

 {
   LogDebug("SimLaserAlignmentSimulation") << "<LaserAlignmentSimulation::update(const EndOfRun * myRun)>";
 
   // stop timer
   theTimer->Stop();
   edm::LogInfo("SimLaserAlignmentSimulation") << " AC1CMS: Number of Events = " << (*myRun)()->GetNumberOfEventToBeProcessed()
                     << " " << *theTimer 
                     << " *****     AC1CMS: End of Run: " << (*myRun)()->GetRunID() << "     ***** ";
 }