2 #include "TVirtualMCApplication.h"
3 #include "TVirtualMCStack.h"
5 #include "TFlukaCodes.h"
8 #include "Fdimpar.h" //(DIMPAR)
9 #include "Fdblprc.h" //(DBLPRC)
10 #include "Ftrackr.h" //(TRACKR)
11 #include "Fopphst.h" //(OPPHST)
12 #include "Fflkstk.h" //(FLKSTK)
13 #include "Fltclcm.h" //(LTCLCM)
14 #include "Fpaprop.h" //(PAPROP)
15 #include "Falldlt.h" //(ALLDLT)
16 #include "Fdpdxcm.h" //(DPDXCM)
17 #include "Fflkmat.h" //(FLKMAT)
19 #include "TGeoManager.h"
22 # define mgdraw mgdraw_
24 # define mgdraw MGDRAW
28 void mgdraw(Int_t& icode, Int_t& mreg)
30 TFluka* fluka = (TFluka*) gMC;
31 if (mreg == fluka->GetDummyRegion()) return;
33 // Make sure that stack has currrent track Id
36 TVirtualMCStack* cppstack = fluka->GetStack();
38 if (TRACKR.jtrack == -1) {
39 trackId = OPPHST.louopp[OPPHST.lstopp];
41 trackId = FLKSTK.ispark[FLKSTK.npflka][mkbmx2-1];
44 trackId = TRACKR.ispusr[mkbmx2-1];
47 Int_t verbosityLevel = fluka->GetVerbosityLevel();
48 if (TRACKR.jtrack < -6) {
49 // from -7 to -12 = "heavy" fragment
51 // id < -6 was skipped in stuprf => if (kpart < -6) return;
52 if (verbosityLevel >= 3) {
53 cout << "mgdraw: (heavy fragment) jtrack < -6 =" << TRACKR.jtrack
54 << " assign parent pdg=" << fluka->PDGFromId(TRACKR.ispusr[mkbmx2 - 3]) << endl;
56 // TRACKR.jtrack = TRACKR.ispusr[mkbmx2 - 3];
59 cppstack->SetCurrentTrack(trackId);
62 Int_t mlttc = TRACKR.lt1trk; // LTCLCM.mlatm1;
63 fluka->SetMreg(mreg, mlttc);
64 fluka->SetIcode((FlukaProcessCode_t) icode);
65 fluka->SetCaller(kMGDRAW);
68 Int_t volId = fluka->CurrentVolID(nodeId);
69 Int_t crtlttc = gGeoManager->GetCurrentNodeId()+1;
71 // check region lattice consistency (debug Ernesto)
72 // *****************************************************
73 if(verbosityLevel>=3 && mreg != volId && !gGeoManager->IsOutside() ) {
74 cout << " mgdraw: track=" << trackId << " pdg=" << fluka->PDGFromId(TRACKR.jtrack)
75 << " icode=" << icode << " gNstep=" << fluka->GetNstep() << endl
76 << " fluka mreg=" << mreg << " mlttc=" << mlttc << endl
77 << " TGeo volId=" << volId << " crtlttc=" << crtlttc << endl
78 << " common TRACKR lt1trk=" << TRACKR.lt1trk << " lt2trk=" << TRACKR.lt2trk << endl
79 << " common LTCLCM newlat=" << LTCLCM.newlat << " mlatld=" << LTCLCM.mlatld << endl
80 << " mlatm1=" << LTCLCM.mlatm1 << " mltsen=" << LTCLCM.mltsen << endl
81 << " mltsm1=" << LTCLCM.mltsm1 << " mlattc=" << LTCLCM.mlattc << endl;
82 if( mlttc == crtlttc ) cout << " *************************************************************" << endl;
84 // *****************************************************
87 Int_t med = FLKMAT.medium[mreg - 1]; // Medium
88 Int_t msd = DPDXCM.msdpdx[med - 1]; // Iionisation model
90 if (!TRACKR.ispusr[mkbmx2 - 2]) {
91 if (verbosityLevel >= 3) {
92 cout << endl << "mgdraw: energy deposition for:" << trackId
94 << " pdg=" << fluka->PDGFromId(TRACKR.jtrack)
95 << " flukaid="<< TRACKR.jtrack
97 << " np =" << ALLDLT.nalldl
104 // Primary ionsiations
105 Int_t nprim = ALLDLT.nalldl;
106 if (nprim >= mxalld) {
108 Warning("mgdraw", "nprim > mxalld, nprim: %6d pdg: %6d mreg %6d p %13.3f step %13.3f\n",
110 fluka->PDGFromId(TRACKR.jtrack),
115 fluka->PrimaryIonisationStepping(nprim);
118 (TVirtualMCApplication::Instance())->Stepping();
119 fluka->SetTrackIsNew(kFALSE);
124 // Tracking is being resumed after secondary tracking
126 if (verbosityLevel >= 3) {
127 cout << endl << "mgdraw: resuming Stepping(): " << trackId << endl;
130 fluka->SetTrackIsNew(kTRUE);
131 fluka->SetCaller(kMGResumedTrack);
132 (TVirtualMCApplication::Instance())->Stepping();
134 // Reset flag and stored values
135 TRACKR.ispusr[mkbmx2 - 2] = 0;
136 for (Int_t i = 0; i < 9; i++) TRACKR.spausr[i] = -1.;
139 if (verbosityLevel >= 3) {
140 cout << endl << " !!! I am in mgdraw - first Stepping() after resume: " << icode << endl;
141 cout << " Track= " << trackId << " region = " << mreg << endl;
144 fluka->SetTrackIsNew(kFALSE);
145 fluka->SetCaller(kMGDRAW);
147 (TVirtualMCApplication::Instance())->Stepping();
149 Int_t nprim = ALLDLT.nalldl;
150 // Protection against nprim > mxalld
151 if (nprim >= mxalld) {
153 Warning("mgdraw", "nprim > mxalld, nprim: %6d pdg: %6d mreg %6d p %13.3f step %13.3f\n",
155 fluka->PDGFromId(TRACKR.jtrack),
160 fluka->PrimaryIonisationStepping(nprim);
161 } // primary ionisation switched on
162 } // tracking resumed
164 } // end of extern "C"