[u/mrichter/AliRoot.git] / TEvtGen / Photos / f_Init.h

#ifndef _f_Init_included_
#define _f_Init_included_

namespace Photospp
{

const static int NMXHEP = 10000;
const static double PI    = 3.14159265358979324;
const static double TWOPI = 6.28318530717958648;


// extern "C"
//{

	/** Definition of the PHOEVT common block */
	extern struct HEPEVT
	{
		int    nevhep;
		int    nhep;
		int    isthep[NMXHEP];
		int    idhep[NMXHEP];
		int    jmohep[NMXHEP][2];
		int    jdahep[NMXHEP][2];
		double phep[NMXHEP][5];
		double vhep[NMXHEP][4];
	        int    qedrad[NMXHEP];  // to be bool once compatibility with F77 removed       
	        const static int nmxhep=NMXHEP;
		//      NEVPHO,NPHO,ISTPHO(NMXPHO),IDPHO(NMXPHO),
		//  JMOPHO(2,NMXPHO),JDAPHO(2,NMXPHO),PPHO(5,NMXPHO),VPHO(4,NMXPHO)
	        //   int qedrad[NMXHEP]  was an add up 
                //   for  HEPEVT in F77 times. Separate common PH_PHOQED
	        //   also phoif_.chkif[NMXPHO] was add up for PHOEVT
                //   now it is pho.qedrad
	} hep,pho;
        //ph_hepevt_,phoevt_;


	extern struct PHOCOP
	{
		double alpha;
		double xphcut;
	} phocop_;

	extern struct PHNUM
	{
		double iev;
	} phnum_;


	extern struct PHOKEY
	{
		double fsec;
		double fint;
		double expeps;
		int interf;
		int isec;
		int itre;
		int iexp;
		int iftop;
		int ifw;
	} phokey_;

	extern struct PHOSTA
	{
		int status[10];
	        int ifstop;
	} phosta_;

	extern struct PHOLUN
	{
		int phlun;
	} pholun_;

	extern struct PHOPHS
	{
		double xphmax;
		double xphoto;
		double costhg;
		double sinthg;

	} phophs_;
	extern struct TOFROM
	{
		double QQ[4];
		double XM;
		double th1;
		double fi1;

	} tofrom_;

	extern struct PHOPRO
	{
		double probh;
		double corwt;
		double xf;
		int irep;
	} phopro_;

	extern struct PHOREST
	{
		double fi3;
		double fi1;
		double th1;
        	int irep;     //    provably line to be removed

	} phorest_;

	extern struct PHWT
	{
		double beta;
		double wt1;
		double wt2;
		double wt3;

	} phwt_;
	extern struct PHOCORWT
	{
		double phocorwt3;
		double phocorwt2;
		double phocorwt1;

	} phocorwt_;

	extern struct PHOMOM
	{
		double mchsqr;
		double mnesqr;
		double pneutr[5];
	} phomom_;
	extern struct PHOCMS
	{
		double bet[3];
		double gam;
	} phocms_;

	extern struct PHOEXP
	{
	        const static int NX = 10;
		double pro[NX];
		int nchan;
	        int expini;    // bool
	} phoexp_;

	//debug mode on if ipoin <  1 and ipoinm > 1
	extern struct PHLUPY
	{
		int ipoin;
		int ipoinm;
	} phlupy_;

	/** Initialize kinematic corrections */
	void PHCORK(int modcor);

	/** Single branch processing */
	void PHOTOS_MAKE_C(int id);

	/* Central management routine. Defines what action
	   will be performed at point ID. */
	void PHTYPE(int ID);
// }

} // namespace Photospp
#endif
Commit	Line	Data
0ca57c2f	1	#ifndef _f_Init_included_
	2	#define _f_Init_included_
	3
	4	namespace Photospp
	5	{
	6
	7	const static int NMXHEP = 10000;
	8	const static double PI = 3.14159265358979324;
	9	const static double TWOPI = 6.28318530717958648;
	10
	11
	12	// extern "C"
	13	//{
	14
	15	/** Definition of the PHOEVT common block */
	16	extern struct HEPEVT
	17	{
	18	int nevhep;
	19	int nhep;
	20	int isthep[NMXHEP];
	21	int idhep[NMXHEP];
	22	int jmohep[NMXHEP][2];
	23	int jdahep[NMXHEP][2];
	24	double phep[NMXHEP][5];
	25	double vhep[NMXHEP][4];
	26	int qedrad[NMXHEP]; // to be bool once compatibility with F77 removed
	27	const static int nmxhep=NMXHEP;
	28	// NEVPHO,NPHO,ISTPHO(NMXPHO),IDPHO(NMXPHO),
	29	// JMOPHO(2,NMXPHO),JDAPHO(2,NMXPHO),PPHO(5,NMXPHO),VPHO(4,NMXPHO)
	30	// int qedrad[NMXHEP] was an add up
	31	// for HEPEVT in F77 times. Separate common PH_PHOQED
	32	// also phoif_.chkif[NMXPHO] was add up for PHOEVT
	33	// now it is pho.qedrad
	34	} hep,pho;
	35	//ph_hepevt_,phoevt_;
	36
	37
	38
	39	extern struct PHOCOP
	40	{
	41	double alpha;
	42	double xphcut;
	43	} phocop_;
	44
	45	extern struct PHNUM
	46	{
	47	double iev;
	48	} phnum_;
	49
	50
	51
	52	extern struct PHOKEY
	53	{
	54	double fsec;
	55	double fint;
	56	double expeps;
	57	int interf;
	58	int isec;
	59	int itre;
	60	int iexp;
	61	int iftop;
	62	int ifw;
	63	} phokey_;
	64
65	extern struct PHOSTA
66	{
67	int status[10];
68	int ifstop;
69	} phosta_;
70
71	extern struct PHOLUN
72	{
73	int phlun;
74	} pholun_;
75
76	extern struct PHOPHS
77	{
78	double xphmax;
79	double xphoto;
80	double costhg;
81	double sinthg;
82
83	} phophs_;
84	extern struct TOFROM
85	{
86	double QQ[4];
87	double XM;
88	double th1;
89	double fi1;
90
91	} tofrom_;
92
93	extern struct PHOPRO
94	{
95	double probh;
96	double corwt;
97	double xf;
98	int irep;
99	} phopro_;
100
101	extern struct PHOREST
102	{
103	double fi3;
104	double fi1;
105	double th1;
106	int irep; // provably line to be removed
107
108	} phorest_;
109
110	extern struct PHWT
111	{
112	double beta;
113	double wt1;
114	double wt2;
115	double wt3;
116
117	} phwt_;
118	extern struct PHOCORWT
119	{
120	double phocorwt3;
121	double phocorwt2;
122	double phocorwt1;
123
124	} phocorwt_;
125
126	extern struct PHOMOM
127	{
128	double mchsqr;
129	double mnesqr;
130	double pneutr[5];
131	} phomom_;
132	extern struct PHOCMS
133	{
134	double bet[3];
135	double gam;
136	} phocms_;
137
138	extern struct PHOEXP
139	{
140	const static int NX = 10;
141	double pro[NX];
142	int nchan;
143	int expini; // bool
144	} phoexp_;
145
146	//debug mode on if ipoin < 1 and ipoinm > 1
147	extern struct PHLUPY
148	{
149	int ipoin;
150	int ipoinm;
151	} phlupy_;
152
153	/** Initialize kinematic corrections */
154	void PHCORK(int modcor);
155
156	/** Single branch processing */
157	void PHOTOS_MAKE_C(int id);
158
159	/* Central management routine. Defines what action
160	will be performed at point ID. */
161	void PHTYPE(int ID);
162	// }
163
164	} // namespace Photospp
165	#endif