QwGeant4
Public Member Functions | Private Attributes
QweakSimCerenkov_RadiatorSD Class Reference

Storing hit information of the Cerenkov radiator. More...

#include <QweakSimCerenkov_RadiatorSD.hh>

Inherits G4VSensitiveDetector.

Public Member Functions

 QweakSimCerenkov_RadiatorSD (G4String name)
 
virtual ~QweakSimCerenkov_RadiatorSD ()
 
void Initialize (G4HCofThisEvent *HCE)
 
G4bool ProcessHits (G4Step *aStep, G4TouchableHistory *ROhist)
 
void EndOfEvent (G4HCofThisEvent *HCE)
 

Private Attributes

QweakSimCerenkovRadiatorHitsCollection * CerenkovRadiatorHitsCollection
 
G4int CerenkovRadiatorCollectionID
 

Detailed Description

Storing hit information of the Cerenkov radiator.

Definition at line 15 of file QweakSimCerenkov_RadiatorSD.hh.

Constructor & Destructor Documentation

QweakSimCerenkov_RadiatorSD::QweakSimCerenkov_RadiatorSD ( G4String  name)

Definition at line 8 of file QweakSimCerenkov_RadiatorSD.cc.

9 : G4VSensitiveDetector(name),
10  CerenkovRadiatorHitsCollection(0),
11  CerenkovRadiatorCollectionID(-1)
12 {
13  collectionName.insert("CerenkovRadiatorCollection");
14 }
QweakSimCerenkov_RadiatorSD::CerenkovRadiatorHitsCollection
QweakSimCerenkovRadiatorHitsCollection * CerenkovRadiatorHitsCollection
Definition: QweakSimCerenkov_RadiatorSD.hh:28
QweakSimCerenkov_RadiatorSD::~QweakSimCerenkov_RadiatorSD ( )
virtual

Definition at line 17 of file QweakSimCerenkov_RadiatorSD.cc.

17 { }

Member Function Documentation

void QweakSimCerenkov_RadiatorSD::EndOfEvent ( G4HCofThisEvent *  HCE)

Definition at line 194 of file QweakSimCerenkov_RadiatorSD.cc.

195 {
196  //G4cout << G4endl << "###### Calling QweakSimCerenkov_RadiatorSD::EndOfEvent() " << G4endl << G4endl;
197 
198 // G4int NbDCHits = DC_hitsCollection->entries();
199 
200 // G4cout << "\n-------->Hits Collection: in this event they are " << NbDCHits
201 // << " hits in the Drift Cells : " << G4endl;
202 // for (G4int i=0;i<NbDCHits;i++) (*DC_hitsCollection)[i]->Print();
203 
204  //G4cout << G4endl << "###### Leaving QweakSimCerenkov_RadiatorSD::EndOfEvent() " << G4endl << G4endl;
205 }
void QweakSimCerenkov_RadiatorSD::Initialize ( G4HCofThisEvent *  HCE)

Definition at line 20 of file QweakSimCerenkov_RadiatorSD.cc.

References CerenkovRadiatorCollectionID, and CerenkovRadiatorHitsCollection.

21 {
22  CerenkovRadiatorHitsCollection = new QweakSimCerenkovRadiatorHitsCollection(SensitiveDetectorName,collectionName[0]);
23 
24  CerenkovRadiatorCollectionID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
25 
26  HCE->AddHitsCollection(CerenkovRadiatorCollectionID , CerenkovRadiatorHitsCollection);
27 }
QweakSimCerenkovRadiatorHitsCollection
G4THitsCollection< QweakSimCerenkov_RadiatorHit > QweakSimCerenkovRadiatorHitsCollection
Definition: QweakSimCerenkov_RadiatorHit.hh:160
QweakSimCerenkov_RadiatorSD::CerenkovRadiatorHitsCollection
QweakSimCerenkovRadiatorHitsCollection * CerenkovRadiatorHitsCollection
Definition: QweakSimCerenkov_RadiatorSD.hh:28
G4bool QweakSimCerenkov_RadiatorSD::ProcessHits ( G4Step *  aStep,
G4TouchableHistory *  ROhist 
)

Definition at line 30 of file QweakSimCerenkov_RadiatorSD.cc.

References CerenkovRadiatorHitsCollection, QweakSimCerenkov_RadiatorHit::GetLogVolume(), QweakSimCerenkov_RadiatorHit::StoreCellPosition(), QweakSimCerenkov_RadiatorHit::StoreCellRotation(), QweakSimCerenkov_RadiatorHit::StoreCreatorProcessName(), QweakSimCerenkov_RadiatorHit::StoreDetectorID(), QweakSimCerenkov_RadiatorHit::StoreGlobalTime(), QweakSimCerenkov_RadiatorHit::StoreKineticEnergy(), QweakSimCerenkov_RadiatorHit::StoreLocalMomentum(), QweakSimCerenkov_RadiatorHit::StoreLocalPosition(), QweakSimCerenkov_RadiatorHit::StoreLogVolume(), QweakSimCerenkov_RadiatorHit::StoreMomentumDirection(), QweakSimCerenkov_RadiatorHit::StoreOriginVertexKineticEnergy(), QweakSimCerenkov_RadiatorHit::StoreOriginVertexMomentumDirection(), QweakSimCerenkov_RadiatorHit::StoreOriginVertexPosition(), QweakSimCerenkov_RadiatorHit::StoreOriginVertexTotalEnergy(), QweakSimCerenkov_RadiatorHit::StoreParentID(), QweakSimCerenkov_RadiatorHit::StoreParticleName(), QweakSimCerenkov_RadiatorHit::StoreParticleType(), QweakSimCerenkov_RadiatorHit::StorePolarization(), QweakSimCerenkov_RadiatorHit::StoreTotalEnergy(), QweakSimCerenkov_RadiatorHit::StoreTrackID(), QweakSimCerenkov_RadiatorHit::StoreWorldMomentum(), and QweakSimCerenkov_RadiatorHit::StoreWorldPosition().

31 {
32 
33 //G4cout << G4endl << "###### Calling QweakSimCerenkov_RadiatorSD::ProcessHits() " << G4endl << G4endl;
34 
35  G4double charge = aStep->GetTrack()->GetDefinition()->GetPDGCharge();
36  if(charge==0. && aStep->GetTrack()->GetTotalEnergy()/MeV < 0.1) {
37 
38  //G4cout << "Return on charge = 0" << G4endl;
39  return false;
40  }
41 
42  if( aStep->GetTrack()->GetDefinition()->GetParticleName() == "opticalphoton") {
43 
44  //std::cout << "Return on optical photon" << std::endl;
45  return false;
46  }
47 
48  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
49 // G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
50  G4TouchableHistory* theTouchable = (G4TouchableHistory*)(preStepPoint->GetTouchable());
51  G4VPhysicalVolume* physVol = theTouchable->GetVolume();
52 
53 // if(strcmp(physVol->GetName(),"CerenkovDetector_Physical")) {
54  if(strcmp(physVol->GetName(),"Radiator_Physical")) {
55  return false;
56  }
57 
58  G4ThreeVector worldPos = preStepPoint->GetPosition();
59  G4ThreeVector localPos = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);
60 
61  //======================================================================================
62  // for test purpose only : print out all the source and primary info of this track
63 
64 // G4cout << "Printing Track Info inside QweakSimCerenkov_RadiatorSD " << G4endl;
65 //
66 // QweakSimTrackInformation* trackInfo;
67 //
68 // trackInfo = (QweakSimTrackInformation*)(aStep->GetTrack()->GetUserInformation());
69 // trackInfo->PrintSourceTrackInfo();
70 // trackInfo->PrintPrimaryTrackInfo();
71 
72  //======================================================================================
73 
74 
75  //G4cout << "PreStep Process in the Cerenkov = " << aStep->GetPreStepPoint()->GetProcessDefinedStep()->GetProcessName() << G4endl;
76  //G4cout << "PostStep Process in the Cerenkov = " << aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName() << G4endl;
77 
78  // here we are only interested if the particle is crossing a boundary
79  // If you want to identify the first step in a volume: pick fGeomBoundary status in preStepPoint
80  // If you want to identify a step out of a volume : pick fGeomBoundary status in postStepPoint
81 
82  if (preStepPoint->GetStepStatus() != fGeomBoundary) {
83 
84  // G4cout << "=======================================================================" << G4endl;
85  // G4cout << ">>>> QweakSimCerenkov_RadiatorSD: We are NOT crossing a boundary <<<<<<" << G4endl;
86  // G4cout << ">>>> Aborting QweakSimCerenkov_RadiatorSD::ProcessHits <<<<<<" << G4endl;
87  // G4cout << "=======================================================================" << G4endl;
88  return false;
89  }
90 
91  G4ThreeVector worldMomentum = preStepPoint->GetMomentum();
92  G4ThreeVector localMomentum = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldMomentum);
93 
94  // Vertex: where this track was generated
95  G4ThreeVector originVertexPosition = aStep->GetTrack()->GetVertexPosition();
96  G4ThreeVector originVertexMomentumDirection = aStep->GetTrack()->GetVertexMomentumDirection();
97  G4double originVertexKineticEnergy = aStep->GetTrack()->GetVertexKineticEnergy();
98 
99  G4double currentKineticEnergy = aStep->GetTrack()->GetKineticEnergy();
100  G4double currentTotalEnergy = aStep->GetTrack()->GetTotalEnergy();
101  G4ThreeVector currentMomentumDirection = aStep->GetTrack()->GetMomentumDirection();
102  G4ThreeVector currentPolarization = aStep->GetTrack()->GetPolarization();
103 
104  G4int trackID = aStep->GetTrack()->GetTrackID();
105  G4int parentID = aStep->GetTrack()->GetParentID();
106 // G4String processName = aStep->GetTrack()->GetCreatorProcess()->GetProcessName();
107  G4String processName = aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName();
108 
109  G4ParticleDefinition* fpParticleDefinition = aStep->GetTrack()->GetDefinition();
110  G4String ParticleName = fpParticleDefinition->GetParticleName();
111 // G4double PDGCharge = fpParticleDefinition->GetPDGCharge();
112  G4int PDGEncoding = fpParticleDefinition->GetPDGEncoding();
113  G4double particlemass = fpParticleDefinition->GetPDGMass()/MeV;
114 
115  G4int MotherCopyNo = theTouchable->GetVolume(1)->GetCopyNo(); // one Mother Volume
116 // G4int DetectorCopyNo = theTouchable->GetVolume()->GetCopyNo(); // but several detectors per MV
117 // G4int DetectorReplicaNo = theTouchable->GetReplicaNumber(); // but several detectors per MV
118 // G4int MotherReplicaNo = theTouchable->GetReplicaNumber(1); // Several MotherVolumes
119 // G4int MotherCopyNo2 = theTouchable->GetVolume(2)->GetCopyNo(); // Several MotherVolumes
120 // G4int MotherReplicaNo2 = theTouchable->GetReplicaNumber(2); // Several MotherVolumes
121 
122 
123 // G4cout << "%%%%%%%%%%%%%%%%%%% Cerenkov MV CopyNumber :" << MotherCopyNo << G4endl;
124 // G4cout << "%%%%%%%%%%%%%%%%%%% Cerenkov CopyNumber :" << DetectorCopyNo << G4endl;
125 // G4cout << "%%%%%%%%%%%%%%%%%%% Cerenkov ReplicaNumber :" << DetectorReplicaNo << G4endl;
126 // G4cout << "%%%%%%%%%%%%%%%%%%% Cerenkov MV ReplicaNumber :" << MotherReplicaNo << G4endl;
127 // G4cout << "%%%%%%%%%%%%%%%%%%% Cerenkov MV CopyNumber2 :" << MotherCopyNo2 << G4endl;
128 // G4cout << "%%%%%%%%%%%%%%%%%%% Cerenkov ReplicaNumber2 :" << MotherReplicaNo2 << G4endl;
129 
130 
131  // get User Track Info
132 // QweakSimTrackInformation* info = (QweakSimTrackInformation*)(aStep->GetTrack()->GetUserInformation());
133 
134 // G4int primaryTrackID = info->GetPrimaryTrackID();
135 // G4double primaryQ2 = info->GetPrimaryQ2();
136 // G4double crossSection = info->GetCrossSection();
137 // G4double crossSectionWeight = info->GetCrossSectionWeight();
138 
139  QweakSimCerenkov_RadiatorHit* aHit = new QweakSimCerenkov_RadiatorHit(MotherCopyNo);
140 
141  aHit->StoreTrackID(trackID);
142  aHit->StoreParentID(parentID);
143  aHit->StoreCreatorProcessName(processName);
144 
145  aHit->StoreDetectorID(MotherCopyNo);
146 
147  aHit->StoreGlobalTime(preStepPoint->GetGlobalTime());
148 
149  aHit->StoreWorldPosition(worldPos);
150  aHit->StoreLocalPosition(localPos);
151 
152  aHit->StoreWorldMomentum(worldMomentum);
153  aHit->StoreLocalMomentum(localMomentum);
154 
155  aHit->StoreOriginVertexPosition(originVertexPosition);
156  aHit->StoreOriginVertexKineticEnergy(originVertexKineticEnergy);
157  //aHit->StoreOriginVertexTotalEnergy(originVertexKineticEnergy); /// \todo beware: total.neq.kinetic (testing only)
158  aHit->StoreOriginVertexTotalEnergy(originVertexKineticEnergy + particlemass);
159  aHit->StoreOriginVertexMomentumDirection(originVertexMomentumDirection);
160 
161  aHit->StoreMomentumDirection(currentMomentumDirection);
162  aHit->StoreKineticEnergy(currentKineticEnergy);
163  aHit->StoreTotalEnergy(currentTotalEnergy);
164  aHit->StorePolarization(currentPolarization);
165 
166  aHit->StoreParticleName(ParticleName);
167  aHit->StoreParticleType(PDGEncoding);
168 
169  // check if it is first touch
170  if(!(aHit->GetLogVolume()))
171  {
172  // store translation and rotation matrix of the drift cell
173  // for the sake of drawing the hit
174  aHit->StoreLogVolume(physVol->GetLogicalVolume());
175  G4AffineTransform aTrans = theTouchable->GetHistory()->GetTopTransform();
176  aTrans.Invert();
177  aHit->StoreCellRotation(aTrans.NetRotation());
178  aHit->StoreCellPosition(aTrans.NetTranslation());
179  }
180 
181 // G4cout << "%%%%%%%%%%%%%%%%%%%%%%% before inserting the hit" << G4endl;
182 
183  CerenkovRadiatorHitsCollection->insert(aHit);
184 
185 // G4cout << G4endl << "###### Leaving QweakSimCerenkovRadiatorSD::ProcessHits() " << G4endl << G4endl;
186 
187 // Peiqing: comment this disable geting secondary info
188 // info->StoreCerenkovHitEnergy(info->GetParticleHistoryLength()-1,currentTotalEnergy);
189 
190  return true;
191 }
QweakSimCerenkov_RadiatorHit::StoreWorldMomentum
void StoreWorldMomentum(G4ThreeVector gpxyz)
Definition: QweakSimCerenkov_RadiatorHit.hh:96
QweakSimCerenkov_RadiatorHit::StoreCellRotation
void StoreCellRotation(G4RotationMatrix rmat)
Definition: QweakSimCerenkov_RadiatorHit.hh:103
QweakSimCerenkov_RadiatorSD::CerenkovRadiatorHitsCollection
QweakSimCerenkovRadiatorHitsCollection * CerenkovRadiatorHitsCollection
Definition: QweakSimCerenkov_RadiatorSD.hh:28
QweakSimCerenkov_RadiatorHit
Handling of a hit in the Cerenkov radiator.
Definition: QweakSimCerenkov_RadiatorHit.hh:14
QweakSimCerenkov_RadiatorHit::StoreLocalMomentum
void StoreLocalMomentum(G4ThreeVector lpxyz)
Definition: QweakSimCerenkov_RadiatorHit.hh:93
QweakSimCerenkov_RadiatorHit::StoreLogVolume
void StoreLogVolume(G4LogicalVolume *val)
Definition: QweakSimCerenkov_RadiatorHit.hh:106
QweakSimCerenkov_RadiatorHit::StoreOriginVertexMomentumDirection
void StoreOriginVertexMomentumDirection(G4ThreeVector opxyz)
Definition: QweakSimCerenkov_RadiatorHit.hh:112
QweakSimCerenkov_RadiatorHit::GetLogVolume
const G4LogicalVolume * GetLogVolume() const
Definition: QweakSimCerenkov_RadiatorHit.hh:107
QweakSimCerenkov_RadiatorHit::StoreMomentumDirection
void StoreMomentumDirection(G4ThreeVector pxyz)
Definition: QweakSimCerenkov_RadiatorHit.hh:123
QweakSimCerenkov_RadiatorHit::StoreOriginVertexPosition
void StoreOriginVertexPosition(G4ThreeVector oxyz)
Definition: QweakSimCerenkov_RadiatorHit.hh:109

+ Here is the call graph for this function:

Field Documentation

G4int QweakSimCerenkov_RadiatorSD::CerenkovRadiatorCollectionID
private

Definition at line 30 of file QweakSimCerenkov_RadiatorSD.hh.

Referenced by Initialize().

QweakSimCerenkovRadiatorHitsCollection* QweakSimCerenkov_RadiatorSD::CerenkovRadiatorHitsCollection
private

Definition at line 28 of file QweakSimCerenkov_RadiatorSD.hh.

Referenced by Initialize(), and ProcessHits().

The documentation for this class was generated from the following files: