[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
#include "Fopphst.h"  //(OPPHST) fluka common
#include "Fstack.h"   //(STACK)  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 = -1;
    TVirtualMCStack* cppstack = fluka->GetStack();
    
    if (TRACKR.jtrack == -1) {
	trackId = OPPHST.louopp[OPPHST.lstopp];
	if (trackId == 0) {
	    trackId = STACK.ispark[STACK.lstack][mkbmx2-1];
	}
    } else {
	trackId = TRACKR.ispusr[mkbmx2-1];
    }
    
    cppstack->SetCurrentTrack(trackId);
//
//    
    fluka->SetMreg(mreg);
    fluka->SetNewreg(mreg);
    fluka->SetIcode(icode);
    fluka->SetCaller(4);

    if (!TRACKR.ispusr[mkbmx2 - 2]) {
	//
	// Single step
	if (verbosityLevel >= 3) {
	    cout << endl << " !!! I am in mgdraw - calling Stepping(): " << icode << endl;
	    cout << endl << " Track Id = " << trackId << " region = " << mreg << endl;
	}
	(TVirtualMCApplication::Instance())->Stepping();
	fluka->SetTrackIsNew(kFALSE);
    } else {
	//
	// Tracking is being resumed after secondary tracking
	//
	// Reset flag
	TRACKR.ispusr[mkbmx2 - 2] = 0;
//
	fluka->SetTrackIsNew(kTRUE);

	if (verbosityLevel >= 3) {
	    cout << endl << " !!! I am in mgdraw - resuming Stepping(): " << trackId << endl;
	}

	(TVirtualMCApplication::Instance())->Stepping();
	fluka->SetTrackIsNew(kFALSE);

	if (verbosityLevel >= 3) {
	    cout << endl << " !!! I am in mgdraw - first Stepping() after resume: " << icode << endl;
	    cout << endl << " Track Id = " << trackId << " region = " << mreg << endl;
	}
	(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"
a7bb59a2	4
fa3d1cc7	5	#include "TFluka.h"
a7bb59a2	6
b1b2005d	7	// Fluka include
	8	#include "Fdimpar.h" //(DIMPAR) fluka include
	9	#include "Fdblprc.h" //(DBLPRC) fluka common
	10	#include "Ftrackr.h" //(TRACKR) fluka common
3a625972	11	#include "Fopphst.h" //(OPPHST) fluka common
3ad1c4d5	12	#include "Fstack.h" //(STACK) fluka common
b1b2005d	13
fa3d1cc7	14	#ifndef WIN32
	15	# define mgdraw mgdraw_
	16	#else
	17	# define mgdraw MGDRAW
	18	#endif
	19
	20	extern "C" {
	21	void mgdraw(Int_t& icode, Int_t& mreg)
	22	{
b1b2005d	23	TFluka* fluka = (TFluka*) gMC;
57dc5a4a	24	Int_t verbosityLevel = fluka->GetVerbosityLevel();
b1b2005d	25	//
b1b2005d	26	// Make sure that stack has currrent track Id
57dc5a4a	27	//
3a625972	28	Int_t trackId = -1;
b1b2005d	29	TVirtualMCStack* cppstack = fluka->GetStack();
3a625972	30
3a625972	31	if (TRACKR.jtrack == -1) {
e45bbee7	32	trackId = OPPHST.louopp[OPPHST.lstopp];
3ad1c4d5	33	if (trackId == 0) {
3ad1c4d5	34	trackId = STACK.ispark[STACK.lstack][mkbmx2-1];
3ad1c4d5	35	}
3a625972	36	} else {
	37	trackId = TRACKR.ispusr[mkbmx2-1];
	38	}
	39
b1b2005d	40	cppstack->SetCurrentTrack(trackId);
	41	//
	42	//
55d6cb0f	43	fluka->SetMreg(mreg);
55d6cb0f	44	fluka->SetNewreg(mreg);
b1b2005d	45	fluka->SetIcode(icode);
a1f42b9c	46	fluka->SetCaller(4);
57dc5a4a	47
	48	if (!TRACKR.ispusr[mkbmx2 - 2]) {
	49	//
	50	// Single step
	51	if (verbosityLevel >= 3) {
	52	cout << endl << " !!! I am in mgdraw - calling Stepping(): " << icode << endl;
	53	cout << endl << " Track Id = " << trackId << " region = " << mreg << endl;
	54	}
	55	(TVirtualMCApplication::Instance())->Stepping();
	56	fluka->SetTrackIsNew(kFALSE);
	57	} else {
	58	//
	59	// Tracking is being resumed after secondary tracking
	60	//
	61	// Reset flag
	62	TRACKR.ispusr[mkbmx2 - 2] = 0;
	63	//
	64	fluka->SetTrackIsNew(kTRUE);
	65
	66	if (verbosityLevel >= 3) {
	67	cout << endl << " !!! I am in mgdraw - resuming Stepping(): " << trackId << endl;
	68	}
	69
	70	(TVirtualMCApplication::Instance())->Stepping();
	71	fluka->SetTrackIsNew(kFALSE);
	72
	73	if (verbosityLevel >= 3) {
	74	cout << endl << " !!! I am in mgdraw - first Stepping() after resume: " << icode << endl;
	75	cout << endl << " Track Id = " << trackId << " region = " << mreg << endl;
	76	}
	77	(TVirtualMCApplication::Instance())->Stepping();
	78	}
3a625972	79
cad96957	80
fa3d1cc7	81	} // end of mgdraw
	82	} // end of extern "C"
	83