Loading...
Searching...
No Matches
BesCgemDigitizer Class Reference
#include <BesCgemDigitizer.hh>
Inheritance diagram for BesCgemDigitizer:
Public Member Functions | |
| BesCgemDigitizer (G4String modName) | |
| ~BesCgemDigitizer () | |
| virtual void | Digitize () |
| void | Digitize_v1 () |
| void | Digitize_v2 () |
| void | Digitize_v3 () |
| void | DoLorentzDiffusion (const G4int f_ID_layer, G4ThreeVector &f_XYZ_pre, G4ThreeVector &f_XYZ_post, const int f_F_print=0) |
| void | Smear (void) |
| void | AddNoise (void) |
| void | GetIDFromXYZ (const G4int f_ID_layer, const G4ThreeVector f_XYZ, G4double &f_x, G4double &f_v, G4int &f_ID_sheet, G4int &f_ID_x, G4int &f_ID_v) |
| void | GetIDFromSegmentInSameSheet (const G4int f_ID_x_pre, const G4int f_ID_v_pre, const G4int f_ID_x_post, const G4int f_ID_v_post, G4int &f_N_hit_strip_x, G4int &f_N_hit_strip_v, G4int &f_ID_x_start, G4int &f_ID_v_start) |
| void | GetMiddleIDFromSegmentCrossedSheet (const G4int f_ID_layer, const G4int f_ID_sheet_pre, const G4double f_x_pre, const G4double f_v_pre, const G4int f_ID_sheet_post, const G4double f_x_post, const G4double f_v_post, G4int &f_ID_sheet_mid_pre, G4int &f_ID_x_mid_pre, G4int &f_ID_v_mid_pre, G4int &f_ID_sheet_mid_post, G4int &f_ID_x_mid_post, G4int &f_ID_v_mid_post) |
| void | GetXYZFromID () |
Detailed Description
Definition at line 61 of file BesCgemDigitizer.hh.
Constructor & Destructor Documentation
◆ BesCgemDigitizer()
| BesCgemDigitizer::BesCgemDigitizer | ( | G4String | modName | ) |
Definition at line 71 of file BesCgemDigitizer.cc.
72:G4VDigitizerModule(modName)
73{
74 /* Digi collection name handled by BesCgemDigitizer */
75 collectionName.push_back("BesCgemDigisCollection");
76
77 ISvcLocator* svcLocator = Gaudi::svcLocator();
78 ICgemGeomSvc* ISvc;
79 StatusCode sc=svcLocator->service("CgemGeomSvc", ISvc);
80
82
83 ISvcLocator* svcLocator2 = Gaudi::svcLocator();
84 ICgemDigitizerSvc* ISvc2;
85 StatusCode sc2=svcLocator2->service("CgemDigitizerSvc", ISvc2);
86
88 m_E_threshold = 0.*eV; /* Energy threshold */
89
90 m_F_lorentz = 1; /* Flag of LorentzDiffusion(): 0-OFF; 1-ON */
91 m_F_smear = 0; /* Flag of Smear(): 0-OFF; 1-ON */
92 m_F_noise = 0; /* Flag of AddNoise(): 0-OFF; 1-ON */
93
94 m_F_printStrip = 0; /* Flag of Print Strip information: 0-OFF; 1-ON */
95 m_F_printHitStrip = 0; /* Flag of Print Hit Strip information: 0-OFF; 1-ON */
96 m_F_printDigi = 0; /* Flag of Print Digi information: 0-OFF; 1-ON */
97 m_F_ntuple = 0; /* Flag of fill ntuple: 0-OFF; 1-ON */
98
99 m_DigitizerVer = 2; /* Flag of Digitizer version: 2: Geo projection, 3: full digitization by CgemDigitizerSvc */
100 static IJobOptionsSvc* jobSvc = 0;
101 sc = svcLocator->service("JobOptionsSvc", jobSvc);
102 if ( sc.isFailure() ) {
103 std::cout << "BesCgemDigitizer: Can't get the JobOptionsSvc " << std::endl;
104 //return sc;
105 }
106 else
107 {
108 const vector<const Property*>* properties = jobSvc->getProperties("BesSim");
109 if ( properties != NULL ) {
110 for ( unsigned int i = 0; i < properties->size(); ++i ) {
111 if ( properties->at(i)->name() == "CgemDigitizer" ) {
112 string strRnd = properties->at(i)->toString();
113 sscanf(strRnd.c_str(), "%i", &m_DigitizerVer);
114 cout<<"BesCgemDigitizer: DigitizerVer is set to "<<m_DigitizerVer<<" by jobOptions"<<endl;
115 break;
116 }
117 }
118 }
119 }
120
121 ISvcLocator* iface = Gaudi::svcLocator();
122
123 if(m_F_ntuple){
125 SmartIF<IDataProviderSvc> eds(iface->service("EventDataSvc"));
126
128 std::cout << "###################" << __LINE__ << " ntupleSvc is not valid"<< std::endl;
129 }
130
132 std::cout << "###################" << __LINE__ << std::endl;
134 if( m_nt1){
135
136 m_nt1->addItem("nevt",m_evt);
137 m_nt1->addItem("nhit",m_nhit,0,1000000);
138 m_nt1->addIndexedItem("layer",m_nhit,m_layer);
139 m_nt1->addIndexedItem("sheet",m_nhit,m_sheet);
140 m_nt1->addIndexedItem("phi",m_nhit,m_phi);
141 m_nt1->addIndexedItem("v",m_nhit,m_v);
142
143 }
144 std::cout << "###################" << __LINE__ << std::endl;
145 }
146 }
147 std::cout << "###################" << __LINE__ << std::endl;
148}
Definition: CgemDigitizerSvc.h:44
Definition: CgemGeomSvc.h:25
Definition: ICgemDigitizerSvc.h:9
Definition: ICgemGeomSvc.h:22
◆ ~BesCgemDigitizer()
| BesCgemDigitizer::~BesCgemDigitizer | ( | ) |
Definition at line 151 of file BesCgemDigitizer.cc.
152{
153}
Member Function Documentation
◆ AddNoise()
| void BesCgemDigitizer::AddNoise | ( | void | ) |
Definition at line 1227 of file BesCgemDigitizer.cc.
1228{
1229}
Referenced by Digitize_v1(), and Digitize_v2().
◆ Digitize()
|
virtual |
Definition at line 156 of file BesCgemDigitizer.cc.
157{
160}
◆ Digitize_v1()
| void BesCgemDigitizer::Digitize_v1 | ( | ) |
Definition at line 161 of file BesCgemDigitizer.cc.
162{
163 IMessageSvc* msgSvc;
166 bool printFlag=false;
167 ISvcLocator* iface = Gaudi::svcLocator();
168 SmartIF<IDataProviderSvc> eds(iface->service("EventDataSvc"));
169 SmartDataPtr<Event::EventHeader> eventHeader(eds,"/Event/EventHeader");
170 if(m_F_ntuple) m_evt = eventHeader->eventNumber();
171 if(printFlag) {
172 log<< MSG::INFO << "|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||" << endreq;
173 log<< MSG::INFO << "INFO : CgemSim::BesCgemDigitizer::Digitize(), Begin !!!" << endreq;
174 }
175
176 /* Flag to check over whether the strip hit before */
177 for (G4int i=0; i<3; i++) /* Layer: 3 */
178 {
179 for (G4int j=0; j<2; j++) /* Sheet: 2 */
180 {
181 for (G4int k=0; k<2; k++) /* Strip Type: 2 */
182 {
183 for (G4int l=0; l<1500; l++) /* StripID: 1419 */
184 {
185 m_F_hit[i][j][k][l] = -1;
186 }
187 }
188 }
189 } /* End of 'for (int i=0; i<3; i++)' */
190
191 /* Digi manager, singleton */
192 G4DigiManager *gv_digi_manager = G4DigiManager::GetDMpointer();
193
194 /* Retrieve Hit Collection ID */
195 G4int lvi_ID_HC = -1;
196 lvi_ID_HC = gv_digi_manager->GetHitsCollectionID("BesCgemHitsCollection");
197
198 /* Retrieve Hit Collection by ID */
199 BesCgemHitsCollection *lv_HC = 0;
200 if (lvi_ID_HC >= 0)
201 {
202 lv_HC = (BesCgemHitsCollection*) (gv_digi_manager->GetHitsCollection(lvi_ID_HC));
203 }
204 else
205 {
206 log<< MSG::ERROR << "ERROR : CgemSim::BesCgemDigitizer::Digitize(), Fail to retrieve HitsCollection!"
207 << endreq;
208 return ;
209 }
210
211 /* Local Variables used by hit and digi in the following */
212 /* Local Variables by hit */
213 G4int lvi_N_hit = 0; /* Number of hits in hit collection */
214 G4int lvi_ID_track = 0; /* Track ID of hit (track, event) */
215 G4int lvi_ID_layer = 0; /* CgemLayer ID of hit */
216 G4double lvd_global_time = 0.; /* Time since the event is created */
217 G4double lvd_E_deposit = 0.; /* Energy deposited */
218 G4double lvd_E_average = 0.; /* Energy average */
219 G4ThreeVector lv3_XYZ_pre (0., 0., 0.); /* Position of pre step point */
220 G4ThreeVector lv3_XYZ_post(0., 0., 0.); /* Position of post step point */
221 //Add by sunxh
222 //G4int lvi_ID_hit = 0;
223 G4int lvi_RdtElectron = 0;
224 //end
225
226 /* Local Variables by CgemLayer */
228 G4int lvi_N_sheet = 0; /* Number of readout sheet */
229 G4int lvi_N_strip[2]; /* 0-X; 1-V */
230 G4double lvd_R_layer = 0.;
231 G4double lvd_L_layer = 0.;
232 G4double lvd_A_stero = 0.; /* Angle of readout V-X */
233 G4double lvd_W_sheet = 0.; /* Width of readout sheet */
234 G4double lvd_W_pitch_x = 0.; /* Width of readout pitch for x*/
235 G4double lvd_W_pitch_v = 0.; /* Width of readout pitch for v*/
236
237 /* Local Variables by Digi */
238 G4double lvd_x_pre, lvd_x_post;
239 G4double lvd_v_pre, lvd_v_post;
240 G4int lvi_ID_sheet_pre , lvi_ID_x_pre, lvi_ID_x_post;
241 G4int lvi_ID_sheet_post, lvi_ID_v_pre, lvi_ID_v_post;
242 G4int lvi_N_hit_strip_x, lvi_N_hit_strip_v;
243 G4int lvi_ID_x_start , lvi_ID_v_start ;
244 G4int lvi_ID_sheet_mid_pre , lvi_ID_x_mid_pre , lvi_ID_v_mid_pre ;
245 G4int lvi_ID_sheet_mid_post, lvi_ID_x_mid_post, lvi_ID_v_mid_post;
246 G4int lvi_N_hit_strip_x_pre , lvi_N_hit_strip_v_pre , lvi_ID_x_start_pre , lvi_ID_v_start_pre ;
247 G4int lvi_N_hit_strip_x_post, lvi_N_hit_strip_v_post, lvi_ID_x_start_post, lvi_ID_v_start_post;
248 G4int lvi_N_hit_strip[2], lvi_N_hit_strip_pre[2];
249 G4int lvi_ID_strip_start[2], lvi_ID_strip_start_pre[2], lvi_ID_strip_start_post[2];
250
251 G4int lvi_ID_sheet = 0;
252 G4int lvi_ID_strip = 0;
253 G4double lvd_E_strip = 0.;
254
255 G4int lvi_case = 0;
256 /* Create digi collection by hit collection */
257 if(lv_HC)
258 {
260
261 /* Process hits and get digi */
262 lvi_N_hit = lv_HC->entries();
263 for (G4int ii=0; ii<lvi_N_hit; ii++)
264 {
265 /* Get information from the hit */
266 BesCgemHit *lv_hit = (*lv_HC)[ii];
267
268 if(printFlag)
269 {
270 log<< MSG::INFO << "..............................................................................."
271 << endreq;
272 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Digitize hit " << ii
273 << " of " << lvi_N_hit << endreq;
274 }
275
276 //Add by sunxh
277 lvi_RdtElectron = lv_hit->GetRdtEl ();
278 //lvi_ID_hit = lv_hit->GetHitID ();
279 //end
280
281 lvi_ID_track = lv_hit->GetTrackID ();
282 lvi_ID_layer = lv_hit->GetLayerID ();
283 lvd_global_time = lv_hit->GetGlobalTime ();
284 lvd_E_deposit = lv_hit->GetTotalEnergyDeposit ();
285 lv3_XYZ_pre = lv_hit->GetPositionOfPrePoint ();
286 lv3_XYZ_post = lv_hit->GetPositionOfPostPoint ();
287
288 lv_cgem_layer = m_cgem_geomsvc->getCgemLayer(lvi_ID_layer);
289 lvd_R_layer = lv_cgem_layer->getInnerROfAnode()*mm;
290 lvd_L_layer = lv_cgem_layer->getLengthOfCgemLayer()*mm;
291 lvi_N_sheet = lv_cgem_layer->getNumberOfSheet();
292 lvd_A_stero = lv_cgem_layer->getAngleOfStereo()*deg;
293 lvd_W_sheet = lv_cgem_layer->getWidthOfSheet()*mm;
294 lvd_W_pitch_x = lv_cgem_layer->getWidthOfPitchX()*mm;
295 lvd_W_pitch_v = lv_cgem_layer->getWidthOfPitchV()*mm;
296
297 lvd_W_sheet = lvd_R_layer * twopi / lvi_N_sheet;
298
299 lvi_N_strip[0] = ceil(lvd_W_sheet / lvd_W_pitch_x);
300 lvi_N_strip[1] = ceil((lvd_W_sheet*cos(lvd_A_stero)+lvd_L_layer*fabs(sin(lvd_A_stero))) / lvd_W_pitch_v);
301
302 //Add by sunxh
303 if(lvi_RdtElectron ==1) continue;// 1: close the effect of delta electron
304 //End Add
305
306 /* Print Strip Information! */
307 if (m_F_printStrip ==1)
308 {
309 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Strip information:" << endreq;
310 log<< MSG::INFO << "ID_layer "
311 << "R "
312 << "N_sheet "
313 << "W_sheet "
314 << "N_X_strip "
315 << "N_V_strip "
316 << "A_stero "
317 << endreq;
318 log<< MSG::INFO << left << setw( 9) << lvi_ID_layer
319 << left << setw( 8) << lvd_R_layer
320 << left << setw( 8) << lvi_N_sheet
321 << left << setw( 8) << lvd_W_sheet
322 << left << setw(10) << lvi_N_strip[0]
323 << left << setw(10) << lvi_N_strip[1]
324 << left << setw( 8) << lvd_A_stero
325 << endreq;
326 log<< MSG::INFO <<"lv3_XYZ_pre, lv3_XYZ_post"<<lv3_XYZ_pre<<", "<<lv3_XYZ_post<<endreq;
327 log<< MSG::INFO << " " << endreq;
328 }
329
330 /* Lorentz diffusion */
331 if (m_F_lorentz == 1)
332 {
333 if(printFlag)
334 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Begin to Lorentz Diffuse!" << endreq;
335 DoLorentzDiffusion(lvi_ID_layer, lv3_XYZ_pre, lv3_XYZ_post);
336 }
337
338 /* Smear */
339 if (m_F_smear == 1)
340 {
341 if(printFlag)
342 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Begin to Smear!" << endreq;
343 Smear();
344 }
345
346 /* Add noise */
347 if (m_F_noise == 1)
348 {
349 if(printFlag)
350 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Begin to Add Noise!" << endreq;
351 AddNoise();
352 }
353
354 /* Get ID from position XYZ of readout layer */
355 /* Pre Point of the step */
356 GetIDFromXYZ(lvi_ID_layer, lv3_XYZ_pre,
357 lvd_x_pre, lvd_v_pre,
358 lvi_ID_sheet_pre, lvi_ID_x_pre, lvi_ID_v_pre);
359
360 /* Post Point of the step */
361 GetIDFromXYZ(lvi_ID_layer, lv3_XYZ_post,
362 lvd_x_post, lvd_v_post,
363 lvi_ID_sheet_post, lvi_ID_x_post, lvi_ID_v_post);
364
365 /* Suppose one step cannot cross over one total sheet */
366 /* Segment in the same sheet */
367 if (lvi_ID_sheet_pre == lvi_ID_sheet_post and
370 {
371 lvi_case = 0;
372 if(printFlag)
373 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Step in same sheet!" << endreq;
374
376 lvi_ID_x_post , lvi_ID_v_post , /* Input */
377 lvi_N_hit_strip_x , lvi_N_hit_strip_v , /* Output */
378 lvi_ID_x_start , lvi_ID_v_start ); /* Output */
379
380 lvi_N_hit_strip[0] = lvi_N_hit_strip_x;
381 lvi_N_hit_strip[1] = lvi_N_hit_strip_v;
382
383 lvi_ID_strip_start[0] = lvi_ID_x_start;
384 lvi_ID_strip_start[1] = lvi_ID_v_start;
385
386 m_phi[ii]=(lvi_ID_x_pre+lvi_ID_x_post)/2*lvd_W_pitch_x;
387 m_v[ii] =(lvi_ID_v_pre+lvi_ID_v_post)/2*lvd_W_pitch_v;
388 m_layer[ii]=lvi_ID_layer;
389 m_sheet[ii]=lvi_ID_sheet_pre;
390 }
391 /* Segment in different sheets, split it into two segemnts which are in same sheet */
392 else
393 {
394 if(printFlag)
395 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Step in different sheet!" << endreq;
396
397 lvi_case = 1;
398 GetMiddleIDFromSegmentCrossedSheet(lvi_ID_layer,
399 lvi_ID_sheet_pre , lvd_x_pre , lvd_v_pre ,
400 lvi_ID_sheet_post, lvd_x_post, lvd_v_post,
401 lvi_ID_sheet_mid_pre , lvi_ID_x_mid_pre , lvi_ID_v_mid_pre,
402 lvi_ID_sheet_mid_post , lvi_ID_x_mid_post , lvi_ID_v_mid_post);
403
405 lvi_ID_x_mid_pre , lvi_ID_v_mid_pre , /* Input */
406 lvi_N_hit_strip_x_pre , lvi_N_hit_strip_v_pre , /* Output */
407 lvi_ID_x_start_pre , lvi_ID_v_start_pre ); /* Output */
408
410 lvi_ID_x_mid_post , lvi_ID_v_mid_post , /* Input */
411 lvi_N_hit_strip_x_post , lvi_N_hit_strip_v_post , /* Output */
412 lvi_ID_x_start_post , lvi_ID_v_start_post ); /* Output */
413
414 if(m_F_ntuple){
415 m_phi[ii]=(lvi_ID_x_pre+lvi_ID_x_mid_pre)/2*0.65;
416 m_v[ii] =(lvi_ID_v_pre+lvi_ID_v_mid_pre)/2*0.65;
417 m_layer[ii]=lvi_ID_layer;
418 m_sheet[ii]=lvi_ID_sheet_mid_pre;
419 }
420
421 lvi_N_hit_strip[0] = lvi_N_hit_strip_x_pre + lvi_N_hit_strip_x_post;
422 lvi_N_hit_strip[1] = lvi_N_hit_strip_v_pre + lvi_N_hit_strip_v_post;
423 lvi_N_hit_strip_pre[0] = lvi_N_hit_strip_x_pre;
424 lvi_N_hit_strip_pre[1] = lvi_N_hit_strip_v_pre;
425 lvi_ID_strip_start_pre[0] = lvi_ID_x_start_pre;
426 lvi_ID_strip_start_pre[1] = lvi_ID_v_start_pre;
427 lvi_ID_strip_start_post[0] = lvi_ID_x_start_post;
428 lvi_ID_strip_start_post[1] = lvi_ID_v_start_post;
429 }
430
431
432 /* Print hit strip information */
433 if (m_F_printHitStrip == 1)
434 {
435 log<< MSG::INFO << " " << endreq;
436 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Hit strip information:"
437 << endreq;
438 log<< MSG::INFO << "PreSheet "
439 << "PostSheet "
440 << "PreX "
441 << "PostX "
442 << "PreV "
443 << "PostV "
444 << "N_X "
445 << "N_V "
446 << "StartX "
447 << "StartV "
448 << endreq;
449 log<< MSG::INFO << left << setw( 9) << lvi_ID_sheet_pre
450 << left << setw(10) << lvi_ID_sheet_post
451 << left << setw( 5) << lvi_ID_x_pre
452 << left << setw( 6) << lvi_ID_x_post
453 << left << setw( 5) << lvi_ID_v_pre
454 << left << setw( 6) << lvi_ID_v_post
455 << left << setw( 4) << lvi_N_hit_strip[0]
456 << left << setw( 4) << lvi_N_hit_strip[1]
457 << left << setw( 7) << lvi_ID_strip_start[0]
458 << left << setw( 7) << lvi_ID_strip_start[1]
459 << endreq;
460 log<< MSG::INFO << " " << endreq;
461 }
462
463 /* uniform energy through sgement pre-post */
464 lvd_E_average = lvd_E_deposit / (lvi_N_hit_strip[0] + lvi_N_hit_strip[1]);
465
466 /* Print Information of Digi */
467 if (m_F_printDigi == 1)
468 {
469 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Digi information:" << endreq;
470 log<< MSG::INFO << left << setw( 8) << "TrackID "
471 << left << setw( 8) << "LayerID "
472 << left << setw( 8) << "SheetID "
473 << left << setw( 8) << "X-0,V-1 "
474 << left << setw( 8) << "StripID "
475 << left << setw(16) << "GlobalTime "
476 << left << setw(16) << "EnergyDeposit "
477 << endreq;
478 }
479
480 /* Add by sunxh */
481 vector<int> ident_tmp;
482 //End Add
483
484 /* Create Digi */
485 for (G4int jj=0; jj<2; jj++) /* 0-X; 1-V */
486 {
487 for (G4int kk=0; kk<lvi_N_hit_strip[jj]; kk++)
488 {
489 if (lvi_case==0)
490 {
491 //G4cout << "### step in same sheet" << G4endl;
492 lvi_ID_sheet = lvi_ID_sheet_pre;
493 lvi_ID_strip = lvi_ID_strip_start[jj] + kk;
494 }
495 else
496 {
497 if (kk < lvi_N_hit_strip_pre[jj])
498 {
499 //G4cout << "*** step in different sheet" << G4endl;
500 lvi_ID_sheet = lvi_ID_sheet_pre;
501 lvi_ID_strip = lvi_ID_strip_start_pre[jj] + kk;
502 }
503 else
504 {
505 //G4cout << "*** step in different sheet" << G4endl;
506 lvi_ID_sheet = lvi_ID_sheet_post;
507 lvi_ID_strip = lvi_ID_strip_start_post[jj] + kk - lvi_N_hit_strip_pre[jj];
508 }
509 }
510
511 /* Compute energy of the hit strip */
512 G4int &idx_digi = m_F_hit[lvi_ID_layer][lvi_ID_sheet][jj][lvi_ID_strip];
513 if (idx_digi == -1)
514 {
515 lvd_E_strip = lvd_E_average ;
516
517 /* Signal should be larger than threshold */
518 if (lvd_E_strip >= m_E_threshold)
519 {
521
522 lv_new_digi->SetTrackID ( lvi_ID_track );
523 lv_new_digi->SetLayerID ( lvi_ID_layer );
524 lv_new_digi->SetSheetID ( lvi_ID_sheet );
525 lv_new_digi->SetStripType ( jj );
526 lv_new_digi->SetStripID ( lvi_ID_strip );
527 lv_new_digi->SetGlobalTime ( lvd_global_time );
528 lv_new_digi->SetEnergyDeposit ( lvd_E_strip );
529
530 m_digi_collection->insert(lv_new_digi);
531
532 idx_digi = m_digi_collection->entries() - 1;
533
534 //Add by sunxh
536 ident_tmp.push_back(ident);
537 //End Add
538
539 /* Print Digi information! */
540 if (m_F_printDigi == 1)
541 {
542 G4cout << left << setw( 8) << lvi_ID_track
543 << left << setw( 8) << lvi_ID_layer
544 << left << setw( 8) << lvi_ID_sheet
545 << left << setw( 8) << jj
546 << left << setw( 8) << lvi_ID_strip
547 << left << setw(16) << lvd_global_time
548 << left << setw(16) << lvd_E_strip
549 << G4endl;
550 }
551 }
552 }
553 else
554 {
555 lvd_E_strip = (*m_digi_collection)[idx_digi]->GetEnergyDeposit() + lvd_E_average;
556
557
558 //Add by sunxh
560 ident_tmp.push_back(Identifier(ident));
561 //End Add
562
563 (*m_digi_collection)[idx_digi]->SetEnergyDeposit(lvd_E_strip);
564 }
565
566 } /* End of 'for (G4int kk=0; kk<lvi_N_hit_strip[jj]; kk++)' */
567 } /* End of 'for (G4int jj=0; jj<2; jj++)' */
568
569 /* Add by sunxh */
570 int array_tmp[2000];
571 for(unsigned int hh =0; hh < ident_tmp.size(); hh++){
572 array_tmp[hh] = ident_tmp[hh];
573 }
574 lv_hit->AddIdentifier(array_tmp,ident_tmp.size());
575 /*End Add*/
576
577 } /* End of 'for (G4int ii=0; ii<lvi_N_hit; ii++)' */
578
579 if(m_F_ntuple) m_nhit = lvi_N_hit;
580 StoreDigiCollection(m_digi_collection); /* */
581
582 } /* End of 'if(lv_HC)' */
583 if(m_F_ntuple) m_nt1->write();
584}
G4TDigiCollection< BesCgemDigi > BesCgemDigisCollection
Definition: BesCgemDigi.hh:91
G4THitsCollection< BesCgemHit > BesCgemHitsCollection
Definition: BesCgemHit.hh:177
Definition: BesCgemDigi.hh:47
void SetEnergyDeposit(G4double f_E_deposit)
Definition: BesCgemDigi.hh:68
void GetIDFromXYZ(const G4int f_ID_layer, const G4ThreeVector f_XYZ, G4double &f_x, G4double &f_v, G4int &f_ID_sheet, G4int &f_ID_x, G4int &f_ID_v)
Definition: BesCgemDigitizer.cc:1233
void DoLorentzDiffusion(const G4int f_ID_layer, G4ThreeVector &f_XYZ_pre, G4ThreeVector &f_XYZ_post, const int f_F_print=0)
Definition: BesCgemDigitizer.cc:1157
void GetIDFromSegmentInSameSheet(const G4int f_ID_x_pre, const G4int f_ID_v_pre, const G4int f_ID_x_post, const G4int f_ID_v_post, G4int &f_N_hit_strip_x, G4int &f_N_hit_strip_v, G4int &f_ID_x_start, G4int &f_ID_v_start)
Definition: BesCgemDigitizer.cc:1297
void GetMiddleIDFromSegmentCrossedSheet(const G4int f_ID_layer, const G4int f_ID_sheet_pre, const G4double f_x_pre, const G4double f_v_pre, const G4int f_ID_sheet_post, const G4double f_x_post, const G4double f_v_post, G4int &f_ID_sheet_mid_pre, G4int &f_ID_x_mid_pre, G4int &f_ID_v_mid_pre, G4int &f_ID_sheet_mid_post, G4int &f_ID_x_mid_post, G4int &f_ID_v_mid_post)
Definition: BesCgemDigitizer.cc:1313
Definition: BesCgemHit.hh:54
G4ThreeVector GetPositionOfPrePoint() const
Definition: BesCgemHit.hh:118
void AddIdentifier(G4int f_ID_Identifier[2000], G4int N_dim)
Definition: BesCgemHit.hh:95
G4ThreeVector GetPositionOfPostPoint() const
Definition: BesCgemHit.hh:119
G4double GetTotalEnergyDeposit() const
Definition: BesCgemHit.hh:116
Definition: CgemGeoLayer.h:9
static value_type getIntID(unsigned int f_layer, unsigned int f_sheet, unsigned int f_strip_type, unsigned int f_strip)
Definition: CgemID.cxx:104
Definition: Identifier.h:20
◆ Digitize_v2()
| void BesCgemDigitizer::Digitize_v2 | ( | ) |
Definition at line 587 of file BesCgemDigitizer.cc.
588{
589 IMessageSvc* msgSvc;
592 bool printFlag=true;
593
594 /* Flag to check over whether the strip hit before */
595 for (G4int i=0; i<3; i++) /* Layer: 3 */
596 {
597 for (G4int j=0; j<2; j++) /* Sheet: 2 */
598 {
599 for (G4int k=0; k<2; k++) /* Strip Type: 2 */
600 {
601 for (G4int l=0; l<1500; l++) /* StripID: 1419 */
602 {
603 m_F_hit[i][j][k][l] = -1;
604 }
605 }
606 }
607 } /* End of 'for (int i=0; i<3; i++)' */
608
609 /* Digi manager, singleton */
610 G4DigiManager *gv_digi_manager = G4DigiManager::GetDMpointer();
611
612 /* Retrieve Hit Collection ID */
613 G4int lvi_ID_HC = -1;
614 lvi_ID_HC = gv_digi_manager->GetHitsCollectionID("BesCgemHitsCollection");
615
616 /* Retrieve Hit Collection by ID */
617 BesCgemHitsCollection *lv_HC = 0;
618 if (lvi_ID_HC >= 0)
619 {
620 lv_HC = (BesCgemHitsCollection*) (gv_digi_manager->GetHitsCollection(lvi_ID_HC));
621 }
622 else
623 {
624 log<< MSG::ERROR << "CgemSim::BesCgemDigitizer::Digitize(), Fail to retrieve HitsCollection!"
625 << endreq;
626 return ;
627 }
628
629 /* Create digi collection by hit collection */
630 if(lv_HC)
631 {
633
634 /* Process hits and get digi */
635 G4int lvi_N_hit = lv_HC->entries();
636 for (G4int ii=0; ii<lvi_N_hit; ii++)
637 {
638 /* Get information from the hit */
639 BesCgemHit *lv_hit = (*lv_HC)[ii];
640
641 //Add by sunxh
642 G4int lvi_RdtElectron = lv_hit->GetRdtEl ();
643 if(lvi_RdtElectron ==1) continue;// 1: close the effect of delta electron
644 //End Add
645
646 G4int lvi_ID_track = lv_hit->GetTrackID ();
647 G4int lvi_ID_layer = lv_hit->GetLayerID ();
648 G4int lvi_pdg = lv_hit->GetPDGCode ();
649 //G4int lvi_ID_parent = lv_hit->GetParentID ();
650 //if(lvi_pdg!=13) continue;
651 G4double lvd_global_time = lv_hit->GetGlobalTime ();
652 G4double lvd_E_deposit = lv_hit->GetTotalEnergyDeposit ();
653 G4ThreeVector lv3_XYZ_pre = lv_hit->GetPositionOfPrePoint ();
654 G4ThreeVector lv3_XYZ_post = lv_hit->GetPositionOfPostPoint ();
656 G4int lvi_N_sheet = lv_cgem_layer->getNumberOfSheet();
657 if(printFlag) {
658 log<< MSG::INFO << "..............................................................................."
659 << endreq;
660 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Digitize hit " << ii
661 << " of " << lvi_N_hit <<", E="<<lvd_E_deposit<<", layer "<<lvi_ID_layer<<", nsheet="<<lvi_N_sheet<< endreq;
662 log<< MSG::INFO << " lvi_pdg "<<lvi_pdg<<", from "<<lv3_XYZ_pre<<" to "<<lv3_XYZ_post<<endreq;
663 }
664
665 /* Lorentz diffusion */
666 if (m_F_lorentz == 1)
667 {
668 if(printFlag)
669 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Begin to Lorentz Diffuse!" << endreq;
670 DoLorentzDiffusion(lvi_ID_layer, lv3_XYZ_pre, lv3_XYZ_post);
671 if(printFlag)
672 log<< MSG::INFO << " after Lorentz Diffuse: from "<<lv3_XYZ_pre<<" to "<<lv3_XYZ_post << endreq;
673 }// lv3_XYZ_pre, lv3_XYZ_post already on readout plane
674
675 /* Smear */
676 if (m_F_smear == 1)
677 {
678 if(printFlag)
679 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Begin to Smear!" << endreq;
680 Smear();
681 }
682
683 /* Add noise */
684 if (m_F_noise == 1)
685 {
686 if(printFlag)
687 log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Begin to Add Noise!" << endreq;
688 AddNoise();
689 }
690
691 vector<int> vecID[2][2];// [iSheet][isXV]
692 int Ntot[2]={0,0};// [isXV]
693 for(int i=0; i<lvi_N_sheet; i++)
694 {
696 readoutPlane->getFiredStripID(lv3_XYZ_pre,lv3_XYZ_post,vecID[i][0],vecID[i][1]);
697 Ntot[0]+=vecID[i][0].size();
698 Ntot[1]+=vecID[i][1].size();
699 }
700 if(printFlag)
701 log<< MSG::INFO << " fired "<<Ntot[0]<<" X strips, "<<Ntot[1]<<" V strips" << endreq;
702 double lvd_E_average=lvd_E_deposit/(Ntot[0]+Ntot[1]);
703
704 /* Add by sunxh */
705 vector<int> ident_tmp;
706 //End Add
707
708 for(int i=0; i<lvi_N_sheet; i++)
709 {
710 //cout<<"* sheet "<<i<<endl;
711 for(int j=0;j<2;j++)// isXV
712 {
713 //cout<<" * XV "<<j<<": "<<endl;
714 for(unsigned int iStrip=0; iStrip<vecID[i][j].size(); iStrip++)
715 {
716 int ID=vecID[i][j][iStrip];
717 G4int &idx_digi=m_F_hit[lvi_ID_layer][i][j][ID];
718 BesCgemDigi *lv_new_digi;
719 double lvd_E_strip=lvd_E_average;
720 //cout<<" ID, E="<<ID<<", "<<lvd_E_strip<<endl;
721 if(idx_digi == -1)
722 {
724 lv_new_digi->SetTrackID ( lvi_ID_track );
725 lv_new_digi->SetLayerID ( lvi_ID_layer );
726 lv_new_digi->SetSheetID ( i);
727 lv_new_digi->SetStripType ( j );
728 lv_new_digi->SetStripID ( ID );
729 lv_new_digi->SetGlobalTime ( lvd_global_time );
730 lv_new_digi->SetEnergyDeposit ( lvd_E_strip );
731 m_digi_collection->insert(lv_new_digi);
732 idx_digi = m_digi_collection->entries() - 1;
733
734 //Add by sunxh
736 ident_tmp.push_back(ident);
737 //End Add
738
739 }
740 else
741 {
742 lvd_E_strip += (*m_digi_collection)[idx_digi]->GetEnergyDeposit();
743
744 //Add by sunxh
746 ident_tmp.push_back(Identifier(ident));
747 //End Add
748
749 (*m_digi_collection)[idx_digi]->SetEnergyDeposit(lvd_E_strip);
750 }
751 }// loop strip
752 }// loop XV
753 }// loop sheet
754
755 /* Add by sunxh */
756 int array_tmp[2000];
757 for(unsigned int hh =0; hh < ident_tmp.size(); hh++){
758 array_tmp[hh] = ident_tmp[hh];
759 }
760 lv_hit->AddIdentifier(array_tmp,ident_tmp.size());
761
762 // TArrayI tmp = lv_hit->GetIdentifier();
763
764 /*End Add*/
765
766 }// loop hits ii
767 StoreDigiCollection(m_digi_collection); /* */
768 }//if(lv_HC)
769}
Definition: CgemGeoReadoutPlane.h:14
void getFiredStripID(G4ThreeVector pos1, G4ThreeVector pos2, vector< int > &vecXID, vector< int > &vecVID) const
Definition: CgemGeoReadoutPlane.cxx:290
CgemGeoReadoutPlane * getReadoutPlane(int iLayer, int iSheet) const
Definition: CgemGeomSvc.h:140
Referenced by Digitize().
◆ Digitize_v3()
| void BesCgemDigitizer::Digitize_v3 | ( | ) |
Definition at line 771 of file BesCgemDigitizer.cc.
772{
773 IMessageSvc* msgSvc;
776 bool printFlag=false; //printFlag=true;
777
778 /* Flag to check over whether the strip hit before */
779 for (G4int i=0; i<3; i++) /* Layer: 3 */
780 {
781 for (G4int j=0; j<2; j++) /* Sheet: 2 */
782 {
783 for (G4int k=0; k<2; k++) /* Strip Type: 2 */
784 {
785 for (G4int l=0; l<1500; l++) /* max StripID: 1395=2790/2 */
786 {
787 m_F_hit[i][j][k][l] = -1;
788 m_cgemHitTimePnt[i][j][k][l] = NULL;
789 }
790 }
791 }
792 } /* End of 'for (int i=0; i<3; i++)' */
793
794 /* Digi manager, singleton */
795 G4DigiManager *gv_digi_manager = G4DigiManager::GetDMpointer();
796
797 /* Retrieve Hit Collection ID */
798 G4int lvi_ID_HC = -1;
799 lvi_ID_HC = gv_digi_manager->GetHitsCollectionID("BesCgemHitsCollection");
800
801 /* Retrieve Hit Collection by ID */
802 BesCgemHitsCollection *lv_HC = 0;
803 if (lvi_ID_HC >= 0)
804 {
805 lv_HC = (BesCgemHitsCollection*) (gv_digi_manager->GetHitsCollection(lvi_ID_HC));
806 }
807 else
808 {
809 log<< MSG::ERROR << "CgemSim::BesCgemDigitizer::Digitize(), Fail to retrieve HitsCollection!"
810 << endreq;
811 return ;
812 }
813
814 /* Create digi collection by hit collection */
815 if(lv_HC)
816 {
818
819 /* Process hits and get digi */
820 G4int lvi_N_hit = lv_HC->entries();
821 for (G4int ii=0; ii<lvi_N_hit; ii++)
822 {
823 /* Get information from the hit */
824 BesCgemHit *lv_hit = (*lv_HC)[ii];
825
826 //Add by sunxh
827 G4int lvi_RdtElectron = lv_hit->GetRdtEl ();
828 if(lvi_RdtElectron ==1) continue;// 1: close the effect of delta electron
829 //End Add
830
831 G4int lvi_ID_track = lv_hit->GetTrackID ();
832 G4int lvi_ID_layer = lv_hit->GetLayerID ();
833 G4int lvi_pdg = lv_hit->GetPDGCode ();
834 G4ThreeVector lvi_momentum = lv_hit->GetMomentumOfPrePoint();
835 //G4int lvi_ID_parent = lv_hit->GetParentID ();
836 //if(lvi_pdg!=13) continue;
837 //G4double lvd_global_time = lv_hit->GetGlobalTime ();
838 G4double lvd_E_deposit = lv_hit->GetTotalEnergyDeposit ();
839 G4ThreeVector lv3_XYZ_pre = lv_hit->GetPositionOfPrePoint ();
840 G4ThreeVector lv3_XYZ_post = lv_hit->GetPositionOfPostPoint ();
842 G4int lvi_N_sheet = lv_cgem_layer->getNumberOfSheet();
843 if(printFlag) {
844 //log<< MSG::INFO << "..............................................................................."
845 // << endreq;
846 //log<< MSG::INFO << "CgemSim::BesCgemDigitizer::Digitize(), Digitize hit " << ii
847 // << " of " << lvi_N_hit <<", E="<<lvd_E_deposit<<", layer "<<lvi_ID_layer<<", nsheet="<<lvi_N_sheet<< endreq;
848 //log<< MSG::INFO << " lvi_pdg "<<lvi_pdg<<", from "<<lv3_XYZ_pre<<" to "<<lv3_XYZ_post<<endreq;
849 cout << "..............................................................................."
850 << endl;
851 cout<< "CgemSim::BesCgemDigitizer::Digitize(), Digitize hit " << ii
852 << " of " << lvi_N_hit <<", E="<<lvd_E_deposit<<", layer "<<lvi_ID_layer<<", nsheet="<<lvi_N_sheet<< endl;
853 cout<< " lvi_pdg "<<lvi_pdg<<", from "<<lv3_XYZ_pre<<" to "<<lv3_XYZ_post<<endl;
854 }
855
856 // if((0==particle)&&(-1==charge)) {particleName = "electron";} add p Lepton e- 11
857 // else if((0==particle)&&(1==charge)) {particleName = "e+";} add p Lepton e+ -11
858 // else if((1==particle)&&(1==charge)) {particleName = "mu+";} add p Lepton mu+ -13
859 // else if((1==particle)&&(-1==charge)) {particleName = "mu-";} add p Lepton mu- 13
860 // else if((2==particle)&&(1==charge)) {particleName = "pi+";} add p Meson pi+ 211
861 // else if((2==particle)&&(-1==charge)) {particleName = "pi-";} add p Meson pi- -211
862 // else if((3==particle)&&(1==charge)) {particleName = "K+";} add p Meson K+ 321
863 // else if((3==particle)&&(-1==charge)) {particleName = "K-";} add p Meson K- -321
864 // else if((4==particle)&&(1==charge)) {particleName = "p";} add p Baryon p+ 2212
865 // else if((4==particle)&&(-1==charge)) {particleName = "p-bar";} add p Baryon anti-p- -2212
866
867 int particle_temp = -1;
868 int charge_temp = 0;
869 if(lvi_pdg == 11){particle_temp=0;charge_temp=-1;}
870 if(lvi_pdg == -11){particle_temp=0;charge_temp=1;}
871 if(lvi_pdg == -13){particle_temp=1;charge_temp=1;}
872 if(lvi_pdg == 13){particle_temp=1;charge_temp=-1;}
873 if(lvi_pdg == 211){particle_temp=2;charge_temp=1;}
874 if(lvi_pdg == -211){particle_temp=2;charge_temp=-1;}
875 if(lvi_pdg == 321){particle_temp=3;charge_temp=1;}
876 if(lvi_pdg == -321){particle_temp=3;charge_temp=-1;}
877 if(lvi_pdg == 2212){particle_temp=4;charge_temp=1;}
878 if(lvi_pdg == -2212){particle_temp=4;charge_temp=-1;}
879 if(lvi_pdg == 1011){particle_temp=5;charge_temp=1;}
880 if(lvi_pdg == 1022){particle_temp=6;charge_temp=2;}
881
882 if(particle_temp==-1) {
883 cout<<"BesCgemDigitizer::Digitize_v3(): no digitization for pdg "<<lvi_pdg<<endl;
884 continue;
885 }
886
887 double posIn[3];
888 double posOut[3];
889 posIn[0] = lv3_XYZ_pre.x();
890 posIn[1] = lv3_XYZ_pre.y();
891 posIn[2] = lv3_XYZ_pre.z();
892 posOut[0] = lv3_XYZ_post.x();
893 posOut[1] = lv3_XYZ_post.y();
894 posOut[2] = lv3_XYZ_post.z();
895
896 //if(sqrt(posIn[0]*posIn[0]+posIn[1]*posIn[1]) < sqrt(posOut[0]*posOut[0]+posOut[1]*posOut[1])){
897
898 /* Get information from digitization service */
899 if(printFlag) { std::cout << "just before m_cgemDigiSvc->setTrack "<<endl;}
900 m_cgemDigiSvc->setTrack(lvi_ID_layer, particle_temp, charge_temp, lvi_momentum.mag()/1000., posIn, posOut);
901 if(printFlag) {
902 std::cout << "layer "<<lvi_ID_layer << ", particle " << particle_temp << ", Q= " << charge_temp << ", p= " << lvi_momentum.mag()/1000.
903 << " GeV/c, (" <<posIn[0] << ", " <<posIn[1]<< ", " <<posIn[2]<< ") to (" <<posOut[0] << ", " <<posOut[1]<< ", " <<posOut[2] <<")"<< std::endl;
904 std::cout << "R_in, R_out = "<< sqrt(posIn[0]*posIn[0]+posIn[1]*posIn[1]) << " " << sqrt(posOut[0]*posOut[0]+posOut[1]*posOut[1]) << std::endl;
905 std::cout << "check--> NXStrips, NVStrips = " << m_cgemDigiSvc->getNXstrips() << m_cgemDigiSvc->getNVstrips()<< std::endl;
906 }
907
908 vector<int> vecID[2][2];// [iSheet][isXV]
909 vector<double> vecQ[2][2]; // [iSheet][isXV]
910 vector<double> vecT[2][2]; // [iSheet][isXV]
911 int Ntot[2]={0,0};// [isXV]
912 //for(int i=0; i<2; i++){
913 // for(int j=0; j<2; j++){
914 // vecID[i][j].clear();
915 // vecQ[i][j].clear();
916 // vecT[i][j].clear();
917 // }
918 //}
919
921 //if(m_cgemDigiSvc->getXstripSheet(i)==0) vecID[0][0].push_back(m_cgemDigiSvc->getXstripID(i));
922 //if(m_cgemDigiSvc->getXstripSheet(i)==1) vecID[1][0].push_back(m_cgemDigiSvc->getXstripID(i));
923 //if(m_cgemDigiSvc->getXstripSheet(i)==0) vecQ[0][0].push_back(m_cgemDigiSvc->getXstripQ(i));
924 //if(m_cgemDigiSvc->getXstripSheet(i)==1) vecQ[1][0].push_back(m_cgemDigiSvc->getXstripQ(i));
925 //if(m_cgemDigiSvc->getXstripSheet(i)==0) vecT[0][0].push_back(m_cgemDigiSvc->getXstripT(i));
926 //if(m_cgemDigiSvc->getXstripSheet(i)==1) vecT[1][0].push_back(m_cgemDigiSvc->getXstripT(i));
927 //std::cout << "zhaojy check sheet" << m_cgemDigiSvc->getXstripSheet(i) << std::endl;
929 vecID[sheetID][0].push_back(m_cgemDigiSvc->getXstripID(i));
930 vecQ[sheetID][0].push_back(m_cgemDigiSvc->getXstripQ(i));
931 vecT[sheetID][0].push_back(m_cgemDigiSvc->getXstripT(i));
932 }
934 //if(m_cgemDigiSvc->getVstripSheet(i)==0) vecID[0][1].push_back(m_cgemDigiSvc->getVstripID(i));
935 //if(m_cgemDigiSvc->getVstripSheet(i)==1) vecID[1][1].push_back(m_cgemDigiSvc->getVstripID(i));
936 //if(m_cgemDigiSvc->getVstripSheet(i)==0) vecQ[0][1].push_back(m_cgemDigiSvc->getVstripQ(i));
937 //if(m_cgemDigiSvc->getVstripSheet(i)==1) vecQ[1][1].push_back(m_cgemDigiSvc->getVstripQ(i));
938 //if(m_cgemDigiSvc->getVstripSheet(i)==0) vecT[0][1].push_back(m_cgemDigiSvc->getVstripT(i));
939 //if(m_cgemDigiSvc->getVstripSheet(i)==1) vecT[1][1].push_back(m_cgemDigiSvc->getVstripT(i));
941 vecID[sheetID][1].push_back(m_cgemDigiSvc->getVstripID(i));
942 vecQ[sheetID][1].push_back(m_cgemDigiSvc->getVstripQ(i));
943 vecT[sheetID][1].push_back(m_cgemDigiSvc->getVstripT(i));
944 }
945
946 for(int i=0; i<lvi_N_sheet; i++)
947 {
948 Ntot[0]+=vecID[i][0].size();
949 Ntot[1]+=vecID[i][1].size();
950 }
951
952 if(printFlag)
953 cout << " fired "<<Ntot[0]<<" X strips, "<<Ntot[1]<<" V strips" << endl;
954 //double lvd_E_average=lvd_E_deposit/(Ntot[0]+Ntot[1]);
955
956 /* Add by sunxh */
957 vector<int> ident_tmp;
958 //End Add
959
960 double Q_saturation= 45.0; //-----Q Satruation
961 for(int i=0; i<lvi_N_sheet; i++)
962 {
963 //cout<<"* sheet "<<i<<endl;
964 for(int j=0;j<2;j++)// isXV
965 {
966 //cout<<" * XV "<<j<<": "<<endl;
967 for(unsigned int iStrip=0; iStrip<vecID[i][j].size(); iStrip++)
968 {
969 int ID=vecID[i][j][iStrip];
970 G4int &idx_digi=m_F_hit[lvi_ID_layer][i][j][ID];
971 BesCgemDigi *lv_new_digi;
972 //double lvd_E_strip=lvd_E_average;
973 //cout<<" ID, E="<<ID<<", "<<lvd_E_strip<<endl;
974 if(idx_digi == -1)
975 {
977 lv_new_digi->SetTrackID ( lvi_ID_track );
978 lv_new_digi->SetLayerID ( lvi_ID_layer );
979 lv_new_digi->SetSheetID ( i);
980 lv_new_digi->SetStripType ( j );
981 lv_new_digi->SetStripID ( ID );
982
983 //if(j==0){ // X strip
984 // if(i==0){ // Sheet 0
985 // for(int k=0; k<vecID[0][0].size(); k++){
986 // if(vecID[0][0][k]==ID){
987 // lv_new_digi->SetGlobalTime ( vecT[0][0][k] );
988 // lv_new_digi->SetEnergyDeposit ( vecQ[0][0][k] );
989 // }
990 // }
991 // }
992 // else if(i==1){ // Sheet 1
993 // for(int k=0; k<vecID[1][0].size(); k++){
994 // if(vecID[1][0][k]==ID){
995 // lv_new_digi->SetGlobalTime ( vecT[1][0][k] );
996 // lv_new_digi->SetEnergyDeposit ( vecQ[1][0][k] );
997 // }
998 // }
999 // }
1000 //}
1001 //else if(j==1){ // V strip
1002 // if(i==0){ // Sheet 0
1003 // for(int k=0; k<vecID[0][1].size(); k++){
1004 // if(vecID[0][1][k]==ID){
1005 // lv_new_digi->SetGlobalTime ( vecT[0][1][k] );
1006 // lv_new_digi->SetEnergyDeposit ( vecQ[0][1][k] );
1007 // }
1008 // }
1009 // }
1010 // else if(i==1){ // Sheet 1
1011 // for(int k=0; k<vecID[1][1].size(); k++){
1012 // if(vecID[1][1][k]==ID){
1013 // lv_new_digi->SetGlobalTime ( vecT[1][1][k] );
1014 // lv_new_digi->SetEnergyDeposit ( vecQ[1][1][k] );
1015 // }
1016 // }
1017 // }
1018 //}
1019 lv_new_digi->SetGlobalTime ( vecT[i][j][iStrip] );
1020 double EnergyDeposit_temp=vecQ[i][j][iStrip];
1021 Q_saturation=m_cgemDigiSvc->getQsaturation(lvi_ID_layer,i,j,ID);
1022 if(EnergyDeposit_temp > Q_saturation) EnergyDeposit_temp = Q_saturation;
1023 lv_new_digi->SetEnergyDeposit ( EnergyDeposit_temp );
1024
1025 m_digi_collection->insert(lv_new_digi);
1026 idx_digi = m_digi_collection->entries() - 1;
1027
1028 // -- update the association between Cgem true hits and digis
1029 lv_hit->AddDigiIdxQ(idx_digi);
1030 lv_hit->SetDigiIdxT(idx_digi);
1031 m_cgemHitTimePnt[lvi_ID_layer][i][j][ID]=lv_hit;
1032
1033 //Add by sunxh
1035 ident_tmp.push_back(ident);
1036 //End Add
1037
1038 } // create a new digi
1039 else
1040 {
1041 //if(j==0){ // X strip
1042 // if(i==0){ // Sheet 0
1043 // for(int k=0; k<vecID[0][0].size(); k++){
1044 // if(vecID[0][0][k]==ID){
1045 // double EnergyDeposit_temp = (*m_digi_collection)[idx_digi]->GetEnergyDeposit()+vecQ[0][0][k];
1046 // double GlobalTime_temp = (*m_digi_collection)[idx_digi]->GetGlobalTime();
1047 // if(GlobalTime_temp > Q_saturation) GlobalTime_temp = Q_saturation;
1048 // (*m_digi_collection)[idx_digi]->SetEnergyDeposit(EnergyDeposit_temp);
1049 // if(vecT[0][0][k]<GlobalTime_temp) (*m_digi_collection)[idx_digi]->SetGlobalTime(vecT[0][0][k]);
1050 // }
1051 // }
1052 // }
1053 // else if(i==1){ // Sheet 1
1054 // for(int k=0; k<vecID[1][0].size(); k++){
1055 // if(vecID[1][0][k]==ID){
1056 // double EnergyDeposit_temp = (*m_digi_collection)[idx_digi]->GetEnergyDeposit()+vecQ[1][0][k];
1057 // double GlobalTime_temp = (*m_digi_collection)[idx_digi]->GetGlobalTime();
1058 // if(GlobalTime_temp > Q_saturation) GlobalTime_temp = Q_saturation;
1059 // (*m_digi_collection)[idx_digi]->SetEnergyDeposit(EnergyDeposit_temp);
1060 // if(vecT[1][0][k]<GlobalTime_temp) (*m_digi_collection)[idx_digi]->SetGlobalTime(vecT[1][0][k]);
1061 // }
1062 // }
1063 // }
1064 //}
1065 //else if(j==1){ // V strip
1066 // if(i==0){ // Sheet 0
1067 // for(int k=0; k<vecID[0][1].size(); k++){
1068 // if(vecID[0][1][k]==ID){
1069 // double EnergyDeposit_temp = (*m_digi_collection)[idx_digi]->GetEnergyDeposit()+vecQ[0][1][k];
1070 // double GlobalTime_temp = (*m_digi_collection)[idx_digi]->GetGlobalTime();
1071 // if(GlobalTime_temp > Q_saturation) GlobalTime_temp = Q_saturation;
1072 // (*m_digi_collection)[idx_digi]->SetEnergyDeposit(EnergyDeposit_temp);
1073 // if(vecT[0][1][k]<GlobalTime_temp) (*m_digi_collection)[idx_digi]->SetGlobalTime(vecT[0][1][k]);
1074 // }
1075 // }
1076 // }
1077 // else if(i==1){ // Sheet 1
1078 // for(int k=0; k<vecID[1][1].size(); k++){
1079 // if(vecID[1][1][k]==ID){
1080 // double EnergyDeposit_temp = (*m_digi_collection)[idx_digi]->GetEnergyDeposit()+vecQ[1][1][k];
1081 // double GlobalTime_temp = (*m_digi_collection)[idx_digi]->GetGlobalTime();
1082 // if(GlobalTime_temp > Q_saturation) GlobalTime_temp = Q_saturation;
1083 // (*m_digi_collection)[idx_digi]->SetEnergyDeposit(EnergyDeposit_temp);
1084 // if(vecT[1][1][k]<GlobalTime_temp) (*m_digi_collection)[idx_digi]->SetGlobalTime(vecT[1][1][k]);
1085 // }
1086 // }
1087 // }
1088 //}
1089
1090 // update energy
1091 double EnergyDeposit_temp = (*m_digi_collection)[idx_digi]->GetEnergyDeposit()+vecQ[i][j][iStrip];
1092 Q_saturation=m_cgemDigiSvc->getQsaturation(lvi_ID_layer,i,j,ID);
1093 if(EnergyDeposit_temp > Q_saturation) EnergyDeposit_temp = Q_saturation;
1094 (*m_digi_collection)[idx_digi]->SetEnergyDeposit(EnergyDeposit_temp);
1095
1096 // update time
1097 double GlobalTime_temp = (*m_digi_collection)[idx_digi]->GetGlobalTime();
1098 if(vecT[i][j][iStrip]<GlobalTime_temp)
1099 {
1100 (*m_digi_collection)[idx_digi]->SetGlobalTime(vecT[i][j][iStrip]);
1101 lv_hit->SetDigiIdxT(idx_digi);
1102 BesCgemHit* lastTCgemHit = m_cgemHitTimePnt[lvi_ID_layer][i][j][ID];
1103 lastTCgemHit->SetDigiIdxT(-1);
1104 m_cgemHitTimePnt[lvi_ID_layer][i][j][ID]=lv_hit;
1105 }
1106 //-----check zhaojy---------------------
1107 //std::cout << "zhaojy check simulation --> layer:" << lvi_ID_layer << " sheet:" << i << " x/v:" << j << " ID:" << ID << " Q:" << (*m_digi_collection)[idx_digi]->GetEnergyDeposit() << " T:" << (*m_digi_collection)[idx_digi]->GetGlobalTime() << std::endl;
1108 //--------------------------------------
1109
1110 // -- update Cgem true Hit
1111 lv_hit->AddDigiIdxQ(idx_digi);
1112
1113 //Add by sunxh
1115 ident_tmp.push_back(Identifier(ident));
1116 //End Add
1117 }// else update an existing digi
1118 }// loop strip
1119 }// loop XV
1120 }// loop sheet
1121
1122 /* Add by sunxh */
1123 int array_tmp[2000];
1124 int size_tmp=ident_tmp.size();
1125 if(size_tmp>2000) size_tmp=2000;
1126 for(unsigned int hh =0; hh < size_tmp; hh++){
1127 array_tmp[hh] = ident_tmp[hh];
1128 }
1129 lv_hit->AddIdentifier(array_tmp,size_tmp);
1130 // TArrayI tmp = lv_hit->GetIdentifier();
1131 /*End Add*/
1132
1133 }// loop hits ii
1134
1135 StoreDigiCollection(m_digi_collection); /* */
1136
1137 //-----check zhaojy---------------------
1138 if(printFlag) {
1139 for(int i=0; i<m_digi_collection->entries(); i++){
1140 std::cout << "check CGEM digitization in simulation --> layer:" << (*m_digi_collection)[i]->GetLayerID()
1141 << " sheet:" << (*m_digi_collection)[i]->GetSheetID()
1142 << " x/v:" << (*m_digi_collection)[i]->GetStripType()
1143 << " ID:" << (*m_digi_collection)[i]->GetStripID()
1144 << " Q:" << (*m_digi_collection)[i]->GetEnergyDeposit()
1145 << " T:" << (*m_digi_collection)[i]->GetGlobalTime() << std::endl;
1146 }
1147 }
1148 //--------------------------------------
1149
1150 }//if(lv_HC)
1151}
G4ThreeVector GetMomentumOfPrePoint() const
Definition: BesCgemHit.hh:122
double getQsaturation(int layer, int sheet, int view, int id)
Definition: CgemDigitizerSvc.h:75
StatusCode setTrack(int layer, int particle, int charge, double p, double trkPosIn[], double trkPosOut[])
Definition: CgemDigitizerSvc.cxx:282
Referenced by Digitize().
◆ DoLorentzDiffusion()
| void BesCgemDigitizer::DoLorentzDiffusion | ( | const G4int | f_ID_layer, |
| G4ThreeVector & | f_XYZ_pre, | ||
| G4ThreeVector & | f_XYZ_post, | ||
| const int | f_F_print = 0 |
||
| ) |
Definition at line 1157 of file BesCgemDigitizer.cc.
1158{
1159 IMessageSvc* msgSvc;
1162 /* !!! Attention! JUST a SIMPLE-TEST-VERSION, it is not a real one!!! */
1164 G4double lvd_R_layer = f_cgem_layer->getInnerROfAnode()*mm;
1165
1166 G4double lvd_A_phi_pre = f_XYZ_pre.phi();
1167 //G4double lvd_A_theta_pre = f_XYZ_pre.theta();
1168 G4double lvd_A_phi_post = f_XYZ_post.phi();
1169 //G4double lvd_A_theta_post = f_XYZ_post.theta();
1170
1171 G4double lvd_XYZ_x_pre = f_XYZ_pre. x() ;
1172 G4double lvd_XYZ_x_post = f_XYZ_post. x() ;
1173 G4double lvd_XYZ_z_pre = f_XYZ_pre. z() ;
1174 G4double lvd_XYZ_z_post = f_XYZ_post. z() ;
1175
1176 /* Print information */
1177 if (f_F_print == 1)
1178 {
1179 log<< MSG::INFO << "Position before Lorentz Diffuse: " << endreq;
1180 log<< MSG::INFO << "PreX=" << f_XYZ_pre.x()
1181 << " PreY=" << f_XYZ_pre.y()
1182 << " PreZ=" << f_XYZ_pre.z()
1183 << " PostX=" << f_XYZ_post.x()
1184 << " PostY=" << f_XYZ_post.y()
1185 << " PostZ=" << f_XYZ_post.z()
1186 << endreq;
1187 log<< MSG::INFO << "" << endreq;
1188 }
1189
1190 f_XYZ_pre.setX ( lvd_R_layer * cos ( lvd_A_phi_pre ));
1191 f_XYZ_pre.setY ( lvd_R_layer * sin ( lvd_A_phi_pre ));
1192 // f_XYZ_pre.setZ ( lvd_R_layer / tan ( lvd_A_theta_pre ));
1193
1194 f_XYZ_post.setX ( lvd_R_layer * cos ( lvd_A_phi_post ));
1195 f_XYZ_post.setY ( lvd_R_layer * sin ( lvd_A_phi_post ));
1196 // f_XYZ_post.setZ ( lvd_R_layer / tan ( lvd_A_theta_post ));
1197
1198 /* Print information */
1199 if (f_F_print == 1)
1200 {
1201 log<< MSG::INFO << "Position after Lorentz Diffuse: " << endreq;
1202 log<< MSG::INFO << "PreX=" << f_XYZ_pre.x()
1203 << " PreY=" << f_XYZ_pre.y()
1204 << " PreZ=" << f_XYZ_pre.z()
1205 << " PostX=" << f_XYZ_post.x()
1206 << " PostY=" << f_XYZ_post.y()
1207 << " PostZ=" << f_XYZ_post.z()
1208 << endreq;
1209 log<< MSG::INFO << "" << endreq;
1210
1211 log<< MSG::INFO << "Postition Delta before and after Lorentz Diffuse: " << endreq;
1216 << endreq;
1217 log<< MSG::INFO << "" << endreq;
1218 }
1219}
Referenced by Digitize_v1(), and Digitize_v2().
◆ GetIDFromSegmentInSameSheet()
| void BesCgemDigitizer::GetIDFromSegmentInSameSheet | ( | const G4int | f_ID_x_pre, |
| const G4int | f_ID_v_pre, | ||
| const G4int | f_ID_x_post, | ||
| const G4int | f_ID_v_post, | ||
| G4int & | f_N_hit_strip_x, | ||
| G4int & | f_N_hit_strip_v, | ||
| G4int & | f_ID_x_start, | ||
| G4int & | f_ID_v_start | ||
| ) |
Definition at line 1297 of file BesCgemDigitizer.cc.
1301{
1302 /* Number of hit strips */
1303 f_N_hit_strip_x = abs(f_ID_x_post - f_ID_x_pre) + 1;
1304 f_N_hit_strip_v = abs(f_ID_v_post - f_ID_v_pre) + 1;
1305
1306 /* Starting strip of hit strip */
1307 f_ID_x_start = (f_ID_x_pre <= f_ID_x_post) ? f_ID_x_pre : f_ID_x_post;
1308 f_ID_v_start = (f_ID_v_pre <= f_ID_v_post) ? f_ID_v_pre : f_ID_v_post;
1309}
Referenced by Digitize_v1().
◆ GetIDFromXYZ()
| void BesCgemDigitizer::GetIDFromXYZ | ( | const G4int | f_ID_layer, |
| const G4ThreeVector | f_XYZ, | ||
| G4double & | f_x, | ||
| G4double & | f_v, | ||
| G4int & | f_ID_sheet, | ||
| G4int & | f_ID_x, | ||
| G4int & | f_ID_v | ||
| ) |
Definition at line 1233 of file BesCgemDigitizer.cc.
1235{
1236 IMessageSvc* msgSvc;
1239 /* CgemLayer information about readout strip */
1241 G4double lvd_R_layer = lv_cgem_layer->getInnerROfAnode()*mm;
1242 G4double lvd_L_layer = lv_cgem_layer->getLengthOfCgemLayer()*mm;
1243 G4double lvd_A_stero = lv_cgem_layer->getAngleOfStereo()*deg;
1244 G4double lvd_W_sheet = lv_cgem_layer->getWidthOfSheet()*mm;
1245 G4double lvd_W_pitch_x = lv_cgem_layer->getWidthOfPitchX()*mm;
1246 G4double lvd_W_pitch_v = lv_cgem_layer->getWidthOfPitchV()*mm;
1247 G4int lvi_N_sheet = lv_cgem_layer->getNumberOfSheet();
1248
1249 lvd_W_sheet = lvd_R_layer * twopi / lvi_N_sheet;
1250
1251 /* lvd_A_phi should be in (0,twopi) */
1252 G4double lvd_A_phi = f_XYZ.phi();
1253 if (lvd_A_phi < 0)
1254 {
1255 lvd_A_phi += twopi;
1256 }
1257 if (lvd_A_phi >= twopi)
1258 {
1259 lvd_A_phi -= twopi;
1260 }
1261
1262 /* Transfer Cylinder to Rectangle */
1263 G4double lvd_x = f_XYZ.rho() * lvd_A_phi;
1264 G4double lvd_v = f_XYZ.z() + lvd_L_layer/2.;
1265
1266
1267 /* ID of Sheet */
1268 if (lvd_x < 0 or lvd_x > lvi_N_sheet*lvd_W_sheet)
1269 {
1270 log<< MSG::ERROR << "ERROR : CgemSim::BesCgemDigitizer::GetIDFromXYZ(), X overlaps range!" << endreq;
1271 log<< MSG::ERROR << "The overlaps X=" << lvd_x
1272 << " Its range should be < " << lvi_N_sheet*lvd_W_sheet
1273 << " at layerID=" << f_ID_layer
1274 << endreq;
1275 return ;
1276 }
1277 else
1278 {
1279 for (G4int i=0; i<lvi_N_sheet; i++)
1280 {
1281 if (lvd_x >= i*lvd_W_sheet and lvd_x < (i+1)*lvd_W_sheet)
1282 {
1283 f_ID_sheet = i;
1284 }
1285 }
1286 }
1287
1288 /* ID of Strip X and V */
1289 f_x = lvd_x - f_ID_sheet * lvd_W_sheet;
1290 f_v = lvd_v;
1291 f_ID_x = floor( f_x / lvd_W_pitch_x );
1293}
Referenced by Digitize_v1().
◆ GetMiddleIDFromSegmentCrossedSheet()
| void BesCgemDigitizer::GetMiddleIDFromSegmentCrossedSheet | ( | const G4int | f_ID_layer, |
| const G4int | f_ID_sheet_pre, | ||
| const G4double | f_x_pre, | ||
| const G4double | f_v_pre, | ||
| const G4int | f_ID_sheet_post, | ||
| const G4double | f_x_post, | ||
| const G4double | f_v_post, | ||
| G4int & | f_ID_sheet_mid_pre, | ||
| G4int & | f_ID_x_mid_pre, | ||
| G4int & | f_ID_v_mid_pre, | ||
| G4int & | f_ID_sheet_mid_post, | ||
| G4int & | f_ID_x_mid_post, | ||
| G4int & | f_ID_v_mid_post | ||
| ) |
Definition at line 1313 of file BesCgemDigitizer.cc.
1318{
1320 G4double lvd_R_layer = lv_cgem_layer->getInnerROfAnode()*mm;
1321 G4double lvd_W_sheet = lv_cgem_layer->getWidthOfSheet()*mm;
1322 G4int lvi_N_sheet = lv_cgem_layer->getNumberOfSheet();
1323 G4double lvd_W_pitch_x = lv_cgem_layer->getWidthOfPitchX()*mm;
1324 G4double lvd_W_pitch_v = lv_cgem_layer->getWidthOfPitchV()*mm;
1325 G4double lvd_A_stero = lv_cgem_layer->getAngleOfStereo()*deg;
1326 G4double lvd_v_mid = 0;
1327
1328 lvd_W_sheet = lvd_R_layer * twopi / lvi_N_sheet;
1329
1330 G4int lvi_N_strip_x = ceil(lvd_W_sheet / lvd_W_pitch_x);
1331
1332 /* Sheet ID of mid point */
1333 f_ID_sheet_mid_pre = f_ID_sheet_pre;
1334 f_ID_sheet_mid_post = f_ID_sheet_post;
1335
1336 /* Strip ID of mid point */
1337 if (f_x_pre > f_x_post)
1338 {
1339 /* X Strip ID of mid point */
1340 f_ID_x_mid_pre = lvi_N_strip_x - 1;
1341 f_ID_x_mid_post = 0;
1342
1343 /* V Strip ID of mid point */
1344 lvd_v_mid = f_x_post*(f_v_pre-f_v_post)/(lvd_W_sheet-f_x_pre+f_x_post) + f_v_post;
1345
1346 f_ID_v_mid_pre = floor((lvd_W_sheet*cos(lvd_A_stero)+lvd_v_mid*sin(lvd_A_stero)) / lvd_W_pitch_v);
1347 f_ID_v_mid_post = floor(lvd_v_mid*sin(lvd_A_stero) / lvd_W_pitch_v);
1348 }
1349 else
1350 {
1351 /* X Strip ID of mid point */
1352 f_ID_x_mid_pre = 0;
1353 f_ID_x_mid_post = lvi_N_strip_x - 1;
1354
1355 /* V Strip ID of mid point */
1356 lvd_v_mid = (lvd_W_sheet-f_x_post)*(f_v_pre-f_v_post)/(lvd_W_sheet-f_x_post+f_x_pre) + f_v_post;
1357
1358 f_ID_v_mid_pre = floor(lvd_v_mid*sin(lvd_A_stero) / lvd_W_pitch_v);
1359 f_ID_v_mid_post = floor((lvd_W_sheet*cos(lvd_A_stero)+lvd_v_mid*sin(lvd_A_stero)) / lvd_W_pitch_v);
1360 }
1361}
Referenced by Digitize_v1().
◆ GetXYZFromID()
| void BesCgemDigitizer::GetXYZFromID | ( | ) |
Definition at line 1365 of file BesCgemDigitizer.cc.
1366{
1367}
◆ Smear()
| void BesCgemDigitizer::Smear | ( | void | ) |
Definition at line 1222 of file BesCgemDigitizer.cc.
1223{
1224}
Referenced by Digitize_v1(), and Digitize_v2().
The documentation for this class was generated from the following files:
