/*
 * DO NOT MANUALLY UPDATE THIS FILE.
 * If needed: make cpp_headers
 *
 * fortran2c.py - automatic conversion
 *       dcdrbs.inc -> dcdrbs.h
 */

#ifndef __DCDRBS__
#define __DCDRBS__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//=== Dcdrbs ===========================================================*
//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//     DeCay DiRection BiaSing:                                         *
//                                                                      *
//     Created on   30 March 1998   by    Alfredo Ferrari & Paola Sala  *
//                                                   Infn - Milan       *
//                                                                      *
//          Aldcdr = Lamda for eps=1-cos(theta). The polar direction    *
//                   around the direction corresponding to the one      *
//                   wished by the user is sample from                  *
//                      P(eps) = exp (-eps/Aldcdr)                      *
//      U,V,Wdcdrb = Lab direction the neutrino should go along         *
//          Sdcdrb = Sin(acos(Wdcdrb))                                  *
//          Spdcdr = Vdcdrb / Sdcdrb                                    *
//          Cpdcdr = Udcdrb / Sdcdrb                                    *
//          Aldcdc = current Lamda for eps=1-cos(theta)                 *
//          Sth0dc = Sin(Theta_0), where (Theta_0,Phi_0) is the direc-  *
//                   tion corresponding to U,V,Wdcdrb in the decaying   *
//                   particle CMS                                       *
//          Cth0dc = Cos(Theta_0)                                       *
//          Sph0dc = Sin(Phi_0)                                         *
//          Cph0dc = Cos(Phi_0)                                         *
//          K12ddb = possible index of (up to two) neutrino types to be *
//                   preferentially biased in multibody decays.         *
//                   K12ddb = k1 + k2 * 1000                            *
//          Kpdcdr = number (1,2,3...) of the decay product whose       *
//                   direction is biased                                *
//       Ldcdrb(i) = flag for decay direction biasing for i_th particle *
//                   type (Paprop numbering scheme)                     *
//          Ldcdbc = current flag for decay direction biasing           *
//          L3bddb = flag for decay direction biasing of 3-body decays  *
//                                                                      *
//----------------------------------------------------------------------*
//
//
extern struct {
	double	aldcdr;
	double	udcdrb;
	double	vdcdrb;
	double	wdcdrb;
	double	sdcdrb;
	double	spdcdr;
	double	cpdcdr;
	double	aldcdc;
	double	sth0dc;
	double	cth0dc;
	double	sph0dc;
	double	cph0dc;
	int	k12ddb;
	int	kpdcdr;
	logical	ldcdrb[NALLWP+6+1];
	logical	ldcdbc;
	logical	l3bddb;
} dcdrbs_;
const int dcdrbs_ldcdrb_base = -6;


#ifdef __cplusplus
}
#endif
#endif