#include <FCALTestbeamSetupSD.hh>

Inheritance diagram for FCALTestbeamSetupSD:

Public Member Functions
	FCALTestbeamSetupSD (G4String)

	~FCALTestbeamSetupSD ()

void	Initialize (G4HCofThisEvent *)

G4bool	ProcessHits (G4Step , G4TouchableHistory )

void	EndOfEvent (G4HCofThisEvent *)

void	clear ()

void	DrawAll ()

void	PrintAll ()

Public Member Functions inherited from G4VSensitiveDetector
	G4VSensitiveDetector (G4String name)

	G4VSensitiveDetector (const G4VSensitiveDetector &right)

virtual	~G4VSensitiveDetector ()

const G4VSensitiveDetector &	operator= (const G4VSensitiveDetector &right)

G4int	operator== (const G4VSensitiveDetector &right) const

G4int	operator!= (const G4VSensitiveDetector &right) const

G4bool	Hit (G4Step *aStep)

void	SetROgeometry (G4VReadOutGeometry *value)

void	SetFilter (G4VSDFilter *value)

G4int	GetNumberOfCollections () const

G4String	GetCollectionName (G4int id) const

void	SetVerboseLevel (G4int vl)

void	Activate (G4bool activeFlag)

G4bool	isActive () const

G4String	GetName () const

G4String	GetPathName () const

G4String	GetFullPathName () const

G4VReadOutGeometry *	GetROgeometry () const

G4VSDFilter *	GetFilter () const

virtual G4VSensitiveDetector *	Clone () const

Data Fields
G4double	EBeamS1

G4double	EBeamS2

G4double	EBeamS3

G4double	EHoleScint

G4double	EBeamHole

G4double	EBeamDead

G4int	TailCatcherID

G4double	ETailVis [8]

G4double	ETailDep [7]

Additional Inherited Members
Protected Member Functions inherited from G4VSensitiveDetector
virtual G4int	GetCollectionID (G4int i)

Protected Attributes inherited from G4VSensitiveDetector
G4CollectionNameVector	collectionName

G4String	SensitiveDetectorName

G4String	thePathName

G4String	fullPathName

G4int	verboseLevel

G4bool	active

G4VReadOutGeometry *	ROgeometry

G4VSDFilter *	filter

Detailed Description

Definition at line 46 of file FCALTestbeamSetupSD.hh.

Constructor & Destructor Documentation

FCALTestbeamSetupSD::FCALTestbeamSetupSD ( G4String name )

Definition at line 54 of file FCALTestbeamSetupSD.cc.

54 : G4VSensitiveDetector(name)

55 {;}

G4VSensitiveDetector::G4VSensitiveDetector

G4VSensitiveDetector(G4String name)

Definition: G4VSensitiveDetector.cc:33

FCALTestbeamSetupSD::~FCALTestbeamSetupSD ( )

Definition at line 59 of file FCALTestbeamSetupSD.cc.

60 {;}

Member Function Documentation

void FCALTestbeamSetupSD::clear ( void )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 175 of file FCALTestbeamSetupSD.cc.

176 {}

void FCALTestbeamSetupSD::DrawAll ( void )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 180 of file FCALTestbeamSetupSD.cc.

181 {}

void FCALTestbeamSetupSD::EndOfEvent ( G4HCofThisEvent * )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 130 of file FCALTestbeamSetupSD.cc.

References testem0::app, EBeamDead, EBeamHole, EBeamS1, EBeamS2, EBeamS3, EHoleScint, ETailDep, ETailVis, G4cout, G4endl, and python.hepunit::MeV.

 {
   G4cout << " Visisble Energy in S1 , S2 , S3 in (MeV)" << G4endl;
   G4cout << EBeamS1/MeV << " " << EBeamS2/MeV << " " << EBeamS3/MeV << " " << G4endl;
 
   G4cout << " Visible Energy in Hole Counter  (MeV) " << G4endl;
   G4cout << EHoleScint/MeV << " " << EBeamHole/MeV << G4endl;
 
   G4cout << " Visible Energy in Upstream Dead Materials " << G4endl;
   G4cout << EBeamDead/MeV << G4endl;
 
   G4cout << " Visible Energy in Tail Catcher Scintillator" << G4endl;
   G4int j;
   for (j=1; j<8 ; j++) {G4cout <<  ETailVis[j]/MeV << " " ;};  G4cout << G4endl;
  
   G4cout << " Visible Energy in Tail Catcher Absorber" << G4endl;
   for (j=1; j<7 ; j++) {G4cout <<  ETailDep[j]/MeV << " " ;};  G4cout << G4endl;
 
 // Write data in File
   //-------------------
   G4String FileName = "Beam_802_1mm.dat";
   std::ios::openmode iostemp;
   if(InitBeam == 0) {
     iostemp = std::ios::out;
     InitBeam++;
   } else {
     iostemp = std::ios::out|std::ios::app;   
   };
   
   std::ofstream BeamDatafile(FileName, iostemp);
   // BeamDatafile.precision(5);
 
   BeamDatafile << EBeamS1/MeV << " " << EBeamS2/MeV << " " << EBeamS3/MeV << std::endl;
   BeamDatafile << EBeamHole/MeV << std::endl;
   BeamDatafile << EBeamDead/MeV << std::endl;
   for (j=1; j<8 ; j++) { BeamDatafile << ETailVis[j]/MeV << " " ;} ; BeamDatafile << std::endl;
   for (j=1; j<7 ; j++) { BeamDatafile << ETailDep[j]/MeV << " " ;} ; BeamDatafile << std::endl;
  
   BeamDatafile.close();
 
        
 }

void FCALTestbeamSetupSD::Initialize ( G4HCofThisEvent * )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 64 of file FCALTestbeamSetupSD.cc.

References EBeamDead, EBeamHole, EBeamS1, EBeamS2, EBeamS3, EHoleScint, ETailDep, ETailVis, and TailCatcherID.

 {
   EBeamS1 = EBeamS2 = EBeamS3 = 0.;
   EHoleScint = EBeamHole = 0.;
   EBeamDead = 0;
   G4int j;
   for (j =0 ; j<8 ; j++) { ETailVis[j] = 0.;};
   for (j =0 ; j<7 ; j++) { ETailDep[j] = 0.;};
   TailCatcherID = 0;
 }

void FCALTestbeamSetupSD::PrintAll ( void )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 185 of file FCALTestbeamSetupSD.cc.

186 {}

G4bool FCALTestbeamSetupSD::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *
	)

virtual

Implements G4VSensitiveDetector.

Definition at line 77 of file FCALTestbeamSetupSD.cc.

References EBeamDead, EBeamHole, EBeamS1, EBeamS2, EBeamS3, EHoleScint, ETailDep, ETailVis, G4VPhysicalVolume::GetCopyNo(), G4VPhysicalVolume::GetName(), G4Step::GetPreStepPoint(), G4Step::GetTotalEnergyDeposit(), G4StepPoint::GetTouchable(), G4TouchableHistory::GetVolume(), and TailCatcherID.

 {
   
   G4double edep = aStep->GetTotalEnergyDeposit();
   if (edep==0.) return false;
 
   G4TouchableHistory* theTouchable
     = (G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable());  
   G4VPhysicalVolume* physVol = theTouchable->GetVolume();
 
   
   if(strcmp(physVol->GetName(),"ScintS1Physical")==0) { 
     EBeamS1 = EBeamS1 + edep;}
   
   if(strcmp(physVol->GetName(),"ScintS2Physical")==0) { 
     EBeamS2 = EBeamS2 + edep;}
   
   if(strcmp(physVol->GetName(),"ScintS3Physical")==0) {  
     EBeamS3 = EBeamS3 + edep;}
 
   if(strcmp(physVol->GetName(),"HoleScintPhysical")==0){ EHoleScint += edep;}
   if(strcmp(physVol->GetName(),"HoleCntrScintPhysical")==0){
     EBeamHole = EBeamHole + edep;}
 
   if(strcmp(physVol->GetName(),"MWPCPhysical")==0) { EBeamDead += edep;};
   if(strcmp(physVol->GetName(),"HoleCntrPbPhysical")==0) { EBeamDead += edep;};
   if(strcmp(physVol->GetName(),"HoleCntrAlPhysical")==0) { EBeamDead += edep;};
   if(strcmp(physVol->GetName(),"LeadWallPhysical")==0) { EBeamDead += edep;};
   if(strcmp(physVol->GetName(),"IronWallPhysical")==0) { EBeamDead += edep;};
 
   if(strcmp(physVol->GetName(),"BigScintPhysical")==0) {
     TailCatcherID = physVol->GetCopyNo();
     if(TailCatcherID > 0)  ETailVis[TailCatcherID] += edep;
   }
   if(strcmp(physVol->GetName(),"SmallScintPhysical")==0) {
     TailCatcherID = physVol->GetCopyNo();
     if(TailCatcherID > 0)  ETailVis[TailCatcherID + 3] += edep;
   }
   
   if(strcmp(physVol->GetName(),"BigIronPhysical")==0) {
     TailCatcherID = physVol->GetCopyNo();
     if(TailCatcherID > 0)  ETailDep[TailCatcherID] += edep;
   }
   if(strcmp(physVol->GetName(),"SmallIronPhysical")==0) {
     TailCatcherID = physVol->GetCopyNo();
     if(TailCatcherID > 0)  ETailDep[TailCatcherID+2] += edep;
   }
 
   return true;
 }