QwGeant4
Public Member Functions | Private Attributes
QweakSimVDC_WirePlaneSD Class Reference

Checks for a valid VDC U-WirePlane and/or V-Wireplane hit and stores the hit information. More...

#include <QweakSimVDC_WirePlaneSD.hh>

Inherits G4VSensitiveDetector.

Public Member Functions

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

Private Attributes

QweakSimVDC_WirePlane_HitsCollection * hitsCollection
 
G4int HCID
 

Detailed Description

Checks for a valid VDC U-WirePlane and/or V-Wireplane hit and stores the hit information.

Placeholder for a long explaination

Definition at line 59 of file QweakSimVDC_WirePlaneSD.hh.

Constructor & Destructor Documentation

QweakSimVDC_WirePlaneSD::QweakSimVDC_WirePlaneSD ( G4String  name)

Definition at line 47 of file QweakSimVDC_WirePlaneSD.cc.

References HCID.

48 :G4VSensitiveDetector(name)
49 {
50 // G4cout << G4endl << "###### Calling QweakSimVDC_WirePlaneSD::QweakSimVDC_WirePlaneSD() " << G4endl << G4endl;
51 
52  collectionName.insert("VDCWirePlaneCollection"); //collectionName[0]
53  HCID = -1;
54 
55 // G4cout << G4endl << "###### Leaving QweakSimVDC_WirePlaneSD::QweakSimVDC_WirePlaneSD() " << G4endl << G4endl;
56 }
QweakSimVDC_WirePlaneSD::~QweakSimVDC_WirePlaneSD ( )

Definition at line 59 of file QweakSimVDC_WirePlaneSD.cc.

60 {
61  //G4cout << G4endl << "###### Calling QweakSimVDC_WirePlaneSD::~QweakSimVDC_WirePlaneSD() " << G4endl << G4endl;
62  //delete hitsCollection;
63  //G4cout << G4endl << "###### Leaving QweakSimVDC_WirePlaneSD::~QweakSimVDC_WirePlaneSD() " << G4endl << G4endl;
64 }

Member Function Documentation

void QweakSimVDC_WirePlaneSD::EndOfEvent ( G4HCofThisEvent *  HCE)

Definition at line 266 of file QweakSimVDC_WirePlaneSD.cc.

267 {
268  //G4cout << G4endl << "###### Calling QweakSimVDC_WirePlaneSD::EndOfEvent() " << G4endl << G4endl;
269 
270 // G4int NbHits = hitsCollection->entries();
271 // G4cout << "\n-------->Hits Collection: in this event they are " << NbHits
272 // << " hits in the Wire Planes: " << G4endl;
273 // for (G4int i=0;i<NbHits;i++) (*hitsCollection)[i]->Print();
274 
275  //G4cout << G4endl << "###### Leaving QweakSimVDC_WirePlaneSD::EndOfEvent() " << G4endl << G4endl;
276 
277 }
void QweakSimVDC_WirePlaneSD::Initialize ( G4HCofThisEvent *  HCE)

Definition at line 67 of file QweakSimVDC_WirePlaneSD.cc.

References HCID, and hitsCollection.

68 {
69  //G4cout << G4endl << "###### Calling QweakSimVDC_WirePlaneSD::Initialize() " << G4endl << G4endl;
70 
71  hitsCollection = new QweakSimVDC_WirePlane_HitsCollection(SensitiveDetectorName,collectionName[0]);
72 
73  // if(HCID<0)
74  HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
75 
76  HCE->AddHitsCollection( HCID, hitsCollection );
77 
78  //G4cout << G4endl << "###### Leaving QweakSimVDC_WirePlaneSD::Initialize() " << G4endl << G4endl;
79 }
QweakSimVDC_WirePlaneSD::hitsCollection
QweakSimVDC_WirePlane_HitsCollection * hitsCollection
Definition: QweakSimVDC_WirePlaneSD.hh:71
QweakSimVDC_WirePlane_HitsCollection
G4THitsCollection< QweakSimVDC_WirePlaneHit > QweakSimVDC_WirePlane_HitsCollection
Definition: QweakSimVDC_WirePlaneHit.hh:153
G4bool QweakSimVDC_WirePlaneSD::ProcessHits ( G4Step *  aStep,
G4TouchableHistory *  ROhist 
)

Definition at line 82 of file QweakSimVDC_WirePlaneSD.cc.

References hitsCollection, QweakSimVDC_WirePlaneHit::StoreGlobalTime(), QweakSimVDC_WirePlaneHit::StoreKineticEnergy(), QweakSimVDC_WirePlaneHit::StoreLocalMomentum(), QweakSimVDC_WirePlaneHit::StoreLocalPosition(), QweakSimVDC_WirePlaneHit::StoreOriginVertexKineticEnergy(), QweakSimVDC_WirePlaneHit::StoreOriginVertexMomentumDirection(), QweakSimVDC_WirePlaneHit::StoreOriginVertexPosition(), QweakSimVDC_WirePlaneHit::StoreParticleName(), QweakSimVDC_WirePlaneHit::StoreParticleType(), QweakSimVDC_WirePlaneHit::StoreTotalEnergy(), QweakSimVDC_WirePlaneHit::StoreVDCID(), QweakSimVDC_WirePlaneHit::StoreWirePlaneID(), QweakSimVDC_WirePlaneHit::StoreWorldMomentum(), and QweakSimVDC_WirePlaneHit::StoreWorldPosition().

83 {
84  //G4cout << G4endl << "###### Calling QweakSimVDC_WirePlaneSD::ProcessHits() " << G4endl << G4endl;
85 
86 
87  G4double charge = aStep->GetTrack()->GetDefinition()->GetPDGCharge();
88  // reject non-charged particle hits
89  if (fabs(charge)<0.1)
90  return false;
91 
92  G4String particlename = aStep->GetTrack()->GetDefinition()->GetParticleName();
93  G4int particletype = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
94 
95 
96 
97  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
98  G4StepPoint* postStepPoint = aStep->GetPreStepPoint();
99  G4TouchableHistory* theTouchable = (G4TouchableHistory*)(preStepPoint->GetTouchable());
100  G4VPhysicalVolume* physVol = theTouchable->GetVolume();
101 
102  // we are only interested if the particle is crossing a boundary (glass->vaccuum)
103  if (postStepPoint->GetStepStatus() != fGeomBoundary) {
104 // G4cout << "=======================================================================" << G4endl;
105 // G4cout << ">>>> QweakSimVDC_WirePlaneSD: We are NOT crossing a boundary <<<<<<" << G4endl;
106 // G4cout << ">>>> Aborting QweakSimVDC_WirePlaneSD::ProcessHits <<<<<<" << G4endl;
107 // G4cout << "=======================================================================" << G4endl;
108 
109  return false;
110  }
111 
112 // G4cout << "================== We are in volume :" << physVol->GetName() << G4endl << G4endl;
113 
114 
115  if( (strcmp(physVol->GetName(),"VDC_UPlane_Physical")==0) || (strcmp(physVol->GetName(),"VDC_VPlane_Physical")==0) ) {
116 
117 
118 // if( (strcmp(physVol->GetName(),"VDC_UPlane_Physical")==0) )
119 // {
120 // G4cout << "=============================================================" << G4endl;
121 // G4cout << ">>>>>>> Particle crossing : VDC_UPlane_Physical <<<<<<<<<" << G4endl;
122 // G4cout << "=============================================================" << G4endl;
123 // }
124 
125 
126 // if( (strcmp(physVol->GetName(),"VDC_VPlane_Physical")==0) )
127 // {
128 // G4cout << "=============================================================" << G4endl;
129 // G4cout << ">>>>>>> Particle crossing : VDC_VPlane_Physical <<<<<<<<<" << G4endl;
130 // G4cout << "=============================================================" << G4endl;
131 // }
132 
133 
134  //G4double currentKineticEnergy = aStep->GetTrack()->GetKineticEnergy();
135  G4double currentKineticEnergy = preStepPoint->GetKineticEnergy();
136 
137  //G4double currentTotalEnergy = aStep->GetTrack()->GetTotalEnergy();
138  G4double currentTotalEnergy = preStepPoint->GetTotalEnergy();
139 
140 
141  //G4ThreeVector currentMomentumDirection = aStep->GetTrack()->GetMomentumDirection();
142  G4ThreeVector currentMomentumDirection = preStepPoint->GetMomentumDirection();
143 
144 // G4int trackID = aStep->GetTrack()->GetTrackID();
145 // G4cout << "====> Track ID : " << trackID << G4endl;
146 
147 // G4int parentID = aStep->GetTrack()->GetParentID();
148 // G4cout << "====> Parent ID : " << parentID << G4endl;
149 
150 
151 
152  //G4int MotherCopyNo = theTouchable->GetVolume(1)->GetCopyNo(); // Several MotherVolumes
153  G4int WirePlaneCopyNo = theTouchable->GetVolume()->GetCopyNo(); // but only one WirePlane per MV
154  //G4int MotherReplicaNo = theTouchable->GetReplicaNumber(1); // Several MotherVolumes
155 
156  //G4int MotherCopyNo2 = theTouchable->GetVolume(2)->GetCopyNo(); // Several MotherVolumes
157  G4int MotherReplicaNo2 = theTouchable->GetReplicaNumber(2); // Several MotherVolumes
158 
159 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane MV CopyNumber :" << MotherCopyNo << G4endl;
160 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane CopyNumber :" << WirePlaneCopyNo << G4endl;
161 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane ReplicaNumber :" << MotherReplicaNo << G4endl;
162 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane MV CopyNumber2 :" << MotherCopyNo2 << G4endl;
163 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane ReplicaNumber2 :" << MotherReplicaNo2 << G4endl;
164 
165 
166  G4ThreeVector worldPos = preStepPoint->GetPosition();
167 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane WorldPos.X [cm]:" << worldPos.x()/cm << G4endl;
168 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane WorldPos.Y [cm]:" << worldPos.y()/cm << G4endl;
169 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane WorldPos.Z [cm]:" << worldPos.z()/cm << G4endl;
170 
171  G4ThreeVector localPos = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);
172 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane LocalPos.X [cm]:" << localPos.x()/cm << G4endl;
173 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane LocalPos.Y [cm]:" << localPos.y()/cm << G4endl;
174 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane LocalPos.Z [cm]:" << localPos.z()/cm << G4endl;
175 
176  G4ThreeVector worldMomentum = preStepPoint->GetMomentum();
177 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane WorldMomentum.X [GeV]:" << worldMomentum.x()/GeV << G4endl;
178 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane WorldMomentum.Y [GeV]:" << worldMomentum.y()/GeV << G4endl;
179 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane WorldMomentum.Z [GeV]:" << worldMomentum.z()/GeV << G4endl;
180 
181  G4ThreeVector localMomentum = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldMomentum);
182 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane LocalMomentum.X [GeV]:" << localMomentum.x()/GeV << G4endl;
183 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane LocalMomentum.Y [GeV]:" << localMomentum.y()/GeV << G4endl;
184 // G4cout << "%%%%%%%%%%%%%%%%%%% WirePlane LocalMomentum.Z [GeV]:" << localMomentum.z()/GeV << G4endl;
185 
186 
187 
188  // Vertex: where this track was generated
189  G4ThreeVector originVertexPosition = aStep->GetTrack()->GetVertexPosition();
190  G4ThreeVector originVertexMomentumDirection = aStep->GetTrack()->GetVertexMomentumDirection();
191  G4double originVertexKineticEnergy = aStep->GetTrack()->GetVertexKineticEnergy();
192 
193  // G4cout << "====> originVertexPosition_X [cm] : " << originVertexPosition.x()/cm << G4endl;
194 // G4cout << "====> originVertexPosition_Y [cm] : " << originVertexPosition.y()/cm << G4endl;
195 // G4cout << "====> originVertexPosition_Z [cm] : " << originVertexPosition.z()/cm << G4endl;
196 
197 // G4cout << "====> originVertexMomentumDirection_X : " << originVertexMomentumDirection.x() << G4endl;
198 // G4cout << "====> originVertexMomentumDirection_Y : " << originVertexMomentumDirection.y() << G4endl;
199 // G4cout << "====> originVertexMomentumDirection_Z : " << originVertexMomentumDirection.z() << G4endl;
200 
201 // G4cout << "====> originVertexKineticEnergy [Mev] : " << originVertexKineticEnergy/MeV << G4endl;
202 
203  // get User Track Info
204  //QweakSimTrackInformation* info = (QweakSimTrackInformation*) (aStep->GetTrack()->GetUserInformation());
205  //info->Print();
206 
207 // G4double primaryQ2 = info->GetPrimaryQ2();
208 // G4double crossSection = info->GetCrossSection();
209 // G4double crossSectionWeight = info->GetCrossSectionWeight();
210 // G4int primaryEventNumber = info->GetPrimaryEventNumber();
211 
212 // G4cout << " @@@@@@@@@@@@@@@@@@@@@@@ Original Track Q2 : " << primaryQ2 << G4endl;
213 // G4cout << " @@@@@@@@@@@@@@@@@@@@@@@ Original Track Event Number : " << primaryEventNumber << G4endl;
214 // G4cout << " @@@@@@@@@@@@@@@@@@@@@@@ Original CrossSection Weight (ub) : " << crossSectionWeight << G4endl;
215 
216 
217 
218 
219  // QweakSimVDCHit* aHit = new QweakSimVDCHit(MotherCopyNo); // there is only one plane per motherVolume
220  QweakSimVDC_WirePlaneHit* aHit = new QweakSimVDC_WirePlaneHit(); // there is only one plane per motherVolume
221 
222  aHit->StoreVDCID(MotherReplicaNo2); // 0: Front VDC , 1: Back VDC
223 
224  // old VDC version
225  //aHit->StoreWirePlaneID(MotherReplicaNo); // 0: Front WirePlane , 1: Back WirePlane
226 
227  // new VDC version with foils etc
228  aHit->StoreWirePlaneID(WirePlaneCopyNo); // 0: Front WirePlane , 1: Back WirePlane
229 
230  aHit->StoreWorldPosition(worldPos);
231  aHit->StoreLocalPosition(localPos);
232  aHit->StoreGlobalTime(preStepPoint->GetGlobalTime());
233 
234  aHit->StoreOriginVertexPosition(originVertexPosition);
235  aHit->StoreOriginVertexMomentumDirection(originVertexMomentumDirection);
236  aHit->StoreOriginVertexKineticEnergy(originVertexKineticEnergy);
237 
238 // aHit->StorePrimaryQ2(primaryQ2);
239 // aHit->StoreCrossSection(crossSection);
240 // aHit->StoreCrossSectionWeight(crossSectionWeight);
241 // aHit->StorePrimaryEventNumber(primaryEventNumber);
242 
243  aHit->StoreWorldMomentum(worldMomentum);
244  aHit->StoreLocalMomentum(localMomentum);
245 
246  aHit->StoreParticleName(particlename);
247  aHit->StoreParticleType(particletype);
248 
249  //aHit->StoreMomentumDirection(currentMomentumDirection);
250  aHit->StoreKineticEnergy(currentKineticEnergy);
251  aHit->StoreTotalEnergy(currentTotalEnergy);
252 
253  // now store all the hit information into the hit collection
254  hitsCollection->insert(aHit);
255 
256  } // end of if(strcmp(physVol->GetName(),"VDC_WirePlane_Physical")==0)
257 
258 
259  //G4cout << G4endl << "###### Leaving QweakSimVDC_WirePlaneSD::ProcessHits() " << G4endl << G4endl;
260 
261  return true;
262 }
QweakSimVDC_WirePlaneHit::StoreOriginVertexKineticEnergy
void StoreOriginVertexKineticEnergy(G4double ekin)
Definition: QweakSimVDC_WirePlaneHit.hh:117
QweakSimVDC_WirePlaneHit::StoreKineticEnergy
void StoreKineticEnergy(G4double ekin)
Definition: QweakSimVDC_WirePlaneHit.hh:143
QweakSimVDC_WirePlaneSD::hitsCollection
QweakSimVDC_WirePlane_HitsCollection * hitsCollection
Definition: QweakSimVDC_WirePlaneSD.hh:71
QweakSimVDC_WirePlaneHit::StoreLocalMomentum
void StoreLocalMomentum(G4ThreeVector lpxyz)
Definition: QweakSimVDC_WirePlaneHit.hh:105
QweakSimVDC_WirePlaneHit::StoreLocalPosition
void StoreLocalPosition(G4ThreeVector lxyz)
Definition: QweakSimVDC_WirePlaneHit.hh:99
QweakSimVDC_WirePlaneHit::StoreWirePlaneID
void StoreWirePlaneID(G4int wireplane_ID)
Definition: QweakSimVDC_WirePlaneHit.hh:93
QweakSimVDC_WirePlaneHit::StoreWorldMomentum
void StoreWorldMomentum(G4ThreeVector gpxyz)
Definition: QweakSimVDC_WirePlaneHit.hh:108
QweakSimVDC_WirePlaneHit::StoreWorldPosition
void StoreWorldPosition(G4ThreeVector gxyz)
Definition: QweakSimVDC_WirePlaneHit.hh:102
QweakSimVDC_WirePlaneHit::StoreOriginVertexMomentumDirection
void StoreOriginVertexMomentumDirection(G4ThreeVector pxyz)
Definition: QweakSimVDC_WirePlaneHit.hh:114
QweakSimVDC_WirePlaneHit
Handling of a U-WirePlane and/or V-WirePlane Hit of the VDC.
Definition: QweakSimVDC_WirePlaneHit.hh:42
QweakSimVDC_WirePlaneHit::StoreOriginVertexPosition
void StoreOriginVertexPosition(G4ThreeVector xyz)
Definition: QweakSimVDC_WirePlaneHit.hh:111

+ Here is the call graph for this function:

Field Documentation

G4int QweakSimVDC_WirePlaneSD::HCID
private

Definition at line 73 of file QweakSimVDC_WirePlaneSD.hh.

Referenced by Initialize(), and QweakSimVDC_WirePlaneSD().

QweakSimVDC_WirePlane_HitsCollection* QweakSimVDC_WirePlaneSD::hitsCollection
private

Definition at line 71 of file QweakSimVDC_WirePlaneSD.hh.

Referenced by Initialize(), and ProcessHits().

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