[u/mrichter/AliRoot.git] / TFluka / mgdraw.cxx

#include <Riostream.h>
#include "TVirtualMCApplication.h"
#include "TVirtualMCStack.h"
#include "TFluka.h"
// Fluka include
#include "Fdimpar.h"  //(DIMPAR) fluka include
#include "Fdblprc.h"  //(DBLPRC) fluka common
#include "Ftrackr.h"  //(TRACKR) fluka common

#ifndef WIN32
# define mgdraw mgdraw_
#else
# define mgdraw MGDRAW
#endif

extern "C" {
void mgdraw(Int_t& icode, Int_t& mreg)
{
    TFluka* fluka =  (TFluka*) gMC;
    Int_t verbosityLevel = fluka->GetVerbosityLevel();
//
//  Make sure that stack has currrent track Id
    Int_t trackId = TRACKR.ispusr[mkbmx2-1];
    TVirtualMCStack* cppstack = fluka->GetStack();
    cppstack->SetCurrentTrack(trackId);
//
//    
    Int_t oldreg = ((TFluka*) gMC)->GetMreg();
    if (oldreg != mreg) {
//
//  Boundary Crossing
//
	fluka->SetNewreg(mreg);
	if (oldreg == -1) fluka->SetMreg(mreg);
	if (verbosityLevel >= 3)
	    printf("Boundary Crossing %d %d \n", oldreg, mreg);
    } else {
	fluka->SetMreg(mreg);
	fluka->SetNewreg(mreg);
	if (verbosityLevel >= 3)
	    printf("Normal step %d %d \n", oldreg, mreg);
    }
    fluka->SetIcode(icode);

    cout << endl << " !!! I am in mgdraw - calling Stepping()" << endl;
    cout << endl << " Track Id =" << trackId << endl;
    
    fluka->FutoTest();

    if (oldreg != mreg) {
//
//  Double step for boundary crossing
//
	fluka->SetTrackIsExiting();
	(TVirtualMCApplication::Instance())->Stepping();
	fluka->SetMreg(mreg);
	fluka->SetTrackIsEntering();
	(TVirtualMCApplication::Instance())->Stepping();
	fluka->SetTrackIsInside();
    } else {
	(TVirtualMCApplication::Instance())->Stepping();
    }
} // end of mgdraw
} // end of extern "C"
Commit	Line	Data
fa3d1cc7	1	#include <Riostream.h>
8fd1d27e	2	#include "TVirtualMCApplication.h"
b1b2005d	3	#include "TVirtualMCStack.h"
fa3d1cc7	4	#include "TFluka.h"
b1b2005d	5	// Fluka include
	6	#include "Fdimpar.h" //(DIMPAR) fluka include
	7	#include "Fdblprc.h" //(DBLPRC) fluka common
	8	#include "Ftrackr.h" //(TRACKR) fluka common
	9
fa3d1cc7	10	#ifndef WIN32
	11	# define mgdraw mgdraw_
	12	#else
	13	# define mgdraw MGDRAW
	14	#endif
	15
	16	extern "C" {
	17	void mgdraw(Int_t& icode, Int_t& mreg)
	18	{
b1b2005d	19	TFluka* fluka = (TFluka*) gMC;
	20	Int_t verbosityLevel = fluka->GetVerbosityLevel();
	21	//
	22	// Make sure that stack has currrent track Id
	23	Int_t trackId = TRACKR.ispusr[mkbmx2-1];
	24	TVirtualMCStack* cppstack = fluka->GetStack();
	25	cppstack->SetCurrentTrack(trackId);
	26	//
	27	//
8fd1d27e	28	Int_t oldreg = ((TFluka*) gMC)->GetMreg();
	29	if (oldreg != mreg) {
	30	//
	31	// Boundary Crossing
	32	//
b1b2005d	33	fluka->SetNewreg(mreg);
	34	if (oldreg == -1) fluka->SetMreg(mreg);
	35	if (verbosityLevel >= 3)
	36	printf("Boundary Crossing %d %d \n", oldreg, mreg);
8fd1d27e	37	} else {
b1b2005d	38	fluka->SetMreg(mreg);
	39	fluka->SetNewreg(mreg);
	40	if (verbosityLevel >= 3)
	41	printf("Normal step %d %d \n", oldreg, mreg);
8fd1d27e	42	}
b1b2005d	43	fluka->SetIcode(icode);
b1b2005d	44
8fd1d27e	45	cout << endl << " !!! I am in mgdraw - calling Stepping()" << endl;
b1b2005d	46	cout << endl << " Track Id =" << trackId << endl;
	47
	48	fluka->FutoTest();
8fd1d27e	49
	50	if (oldreg != mreg) {
	51	//
	52	// Double step for boundary crossing
	53	//
b1b2005d	54	fluka->SetTrackIsExiting();
8fd1d27e	55	(TVirtualMCApplication::Instance())->Stepping();
b1b2005d	56	fluka->SetMreg(mreg);
b1b2005d	57	fluka->SetTrackIsEntering();
8fd1d27e	58	(TVirtualMCApplication::Instance())->Stepping();
b1b2005d	59	fluka->SetTrackIsInside();
8fd1d27e	60	} else {
	61	(TVirtualMCApplication::Instance())->Stepping();
	62	}
fa3d1cc7	63	} // end of mgdraw
	64	} // end of extern "C"
	65