- ((TFluka*) gMC)->SetIcode(icode);
- ((TFluka*) gMC)->SetMreg(mreg);
- ((TFluka*) gMC)->SetNewreg(newreg);
- ((TFluka*) gMC)->SetXsco(xsco);
- ((TFluka*) gMC)->SetYsco(ysco);
- ((TFluka*) gMC)->SetZsco(zsco);
-// cout << endl << " !!! I am in bxdraw - calling Stepping()" << mreg << endl;
- ((TFluka*) gMC)->FutoTest();
- (TVirtualMCApplication::Instance())->Stepping();
+ TFluka* fluka = (TFluka*) gMC;
+ Int_t oldlttc = LTCLCM.mlatm1;
+ Int_t newlttc = LTCLCM.newlat;
+ fluka->SetIcode((FlukaProcessCode_t)icode);
+ fluka->SetNewreg(newreg,newlttc);
+ fluka->SetXsco(xsco);
+ fluka->SetYsco(ysco);
+ fluka->SetZsco(zsco);
+ Int_t verbosityLevel = fluka->GetVerbosityLevel();
+ Bool_t debug = (verbosityLevel>=3)?kTRUE:kFALSE;
+//
+// Double step for boundary crossing
+//
+ fluka->SetTrackIsNew(kFALSE); // has to be called BEFORE Stepping()
+ if (mreg != fluka->GetDummyRegion()) {
+ if (debug) printf("bxdraw (ex) \n");
+ fluka->SetTrackIsExiting();
+ fluka->SetCaller(kBXExiting);
+ fluka->SetMreg(mreg,oldlttc);
+ (TVirtualMCApplication::Instance())->Stepping();
+ }
+ if (newreg != fluka->GetDummyRegion()) {
+ if (debug) printf("bxdraw (en) \n");
+ fluka->SetCaller(kBXEntering);
+ fluka->SetTrackIsEntering();
+ if (fluka->GetDummyBoundary() == 1) fluka->SetDummyBoundary(2);
+ fluka->SetMreg(newreg,newlttc);
+ (TVirtualMCApplication::Instance())->Stepping();
+ }
+