[u/mrichter/AliRoot.git] / STRUCT / AliALIFE.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

//----------------------------------------------------------------------
// Helper class to write geometry in ALIFE format
// in parallel with Geant geometry definition.
// Author: A.Morsch
//----------------------------------------------------------------------

#include <AliALIFE.h>

ClassImp(AliALIFE)


AliALIFE::AliALIFE(const char* file1, const char* file2)
{  
// Constructor
    fNBodies = 0; 
    fNVolumes= 0; 
    fBodyFile   = file1;   // File for Fluka bodies
    fVolumeFile = file2;
    fFile1=fopen(fBodyFile,"w");
    fFile2=fopen(fVolumeFile,"w");
    BodyHeader();
    VolumeHeader();
    fDefaultVolume1 = "DEFAU1";
    fDefaultVolume2 = "DEFAU2";    
}

AliALIFE::AliALIFE()
{
// Default constructor
    fBodyFile   = "FlukaBody.inp";  
    fVolumeFile = "FlukaVolume.inp";
    fFile1=fopen(fBodyFile,"w");
    fFile2=fopen(fVolumeFile,"w");
    BodyHeader();
    VolumeHeader();
}

void AliALIFE::BodyHeader()
{
// Write header for body definitions
    fprintf(fFile1,"*\n");        
    fprintf(fFile1,"GEOBEGIN                                                              COMBINAT\n");
    fprintf(fFile1,"    0    0         AliRoot Generated\n");
    fprintf(fFile1,"*\n");
    fprintf(fFile1,"*== Body Definitions ================================================\n");
    fprintf(fFile1,"*\n");
                            
}


void AliALIFE::VolumeHeader()
{
// Write header for region (volume)  definitions
    fprintf(fFile2,"REGIONS\n");
    fprintf(fFile2,"*\n");
    fprintf(fFile2,"*== Region Definitions =================================================\n");
    fprintf(fFile2,"*\n");
}


void AliALIFE:: Cylinder(Float_t rmin, Float_t rmax, 
			 Float_t zmin, Float_t zmax, 
			 Float_t pos[3], 
			 const char* Material,
			 const char* Field,
			 const char* Cuts) 
{
// Simple cylinder
//
// Bodies
// ^^^^^^
    char name1[5], name2[5], name3[5], name4[5];

//  outer radius
    sprintf(name1, "R%4.4d", fNBodies++);
    fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",name1, pos[0], pos[1], rmax); // inner radius
    sprintf(name2, "R%4.4d", fNBodies++);
    fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",name2, pos[0], pos[1], rmin); 
// z-max
    sprintf(name3, "Z%4.4d", fNBodies++);
    fprintf(fFile1,"%5s XYP%10.3f\n",name3, zmax); 
// z-min
    sprintf(name4, "Z%4.4d", fNBodies++);
    fprintf(fFile1,"%5s XYP%10.3f\n",name4, zmin); 
//
// Volumes
// ^^^^^^^

    fprintf(fFile2,">%s:%s\n", Material, Field);
    fprintf(fFile2,"EMFCUT=%s\n", Cuts);
    fprintf(fFile2,"WW-FACTOR=%s\n", Cuts);
    if (rmin >0) {
	fprintf(fFile2,"+%s-%s+%s-%s\n", name1, name2, name3, name4);
    } else {
	fprintf(fFile2,"+%s+%s-%s\n", name1, name3, name4);
    }
    
    fprintf(fFile2,"\n");    
}

void AliALIFE::OnionCylinder(Float_t* r, Int_t nr, Float_t zmin, Float_t zmax,
			     Float_t pos[3], const char** Materials, 
			     const char** Fields, const char** Cuts) 
{
//
// Concentric cylinders
//

// Bodies
// ^^^^^^
    char nameRin[6], nameRou[6], nameZin[6], nameZou[6];
// z-limits
// z-max
    sprintf(nameZou, "Z%4.4d", fNBodies++);
    nameZou[5]='\0';
    
    fprintf(fFile1,"%5s XYP%10.3f\n",nameZou, zmax); 
// z-min
    sprintf(nameZin, "Z%4.4d", fNBodies++);
    nameZin[5]='\0';
    fprintf(fFile1,"%5s XYP%10.3f\n",nameZin, zmin); 

// inner radius
    sprintf(nameRin, "R%4.4d", fNBodies++);
    nameRin[5]='\0';
    fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",nameRin, pos[0], pos[1], r[0]); 

    Int_t i;
    for (i=1; i<nr; i++) {
//  outer radius
	sprintf(nameRou, "R%4.4d", fNBodies++);
	nameRou[5]='\0';
	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
		nameRou, pos[0], pos[1], r[i]); 	
//
// Volumes
// ^^^^^^^
	if (Fields && Cuts) {
	    fprintf(fFile2,">%s:%s\n", Materials[i-1], Fields[i-1]);
	    fprintf(fFile2,"EMFCUT=%s\n", Cuts[i-1]);
	    fprintf(fFile2,"WW-FACTOR=%s\n", Cuts[i-1]);
	} else {
	    fprintf(fFile2,">%s:%s\n", Materials[i-1], "MF");
	    fprintf(fFile2,"EMFCUT=%s\n", "$UNSHIELDED");
	    fprintf(fFile2,"WW-FACTOR=%s\n", "$UNSHIELDED");
	}
	if (r[i-1] != 0.) {
	    fprintf(fFile2,"+%5s-%5s+%5s-%5s\n", nameZou, nameZin, nameRou, nameRin);
	} else {
	    fprintf(fFile2,"+%5s-%5s+%5s\n", nameZou, nameZin, nameRou);
	}
	fprintf(fFile2,"\n");
	strcpy(nameRin,nameRou);
    }
}


void AliALIFE::Cone(Float_t rmin1, Float_t rmin2, 
		     Float_t rmax1, Float_t rmax2,
		     Float_t zmin, Float_t zmax, 
		     Float_t pos[3], 
		     const char* Material,
		     const char* Field,
		     const char* Cuts) 
{
// Simple cone 

//
// Bodies
// ^^^^^^
    char nameCou[6], nameCin[6];
    Float_t d, r1, r2;
    char nameZin[6], nameZou[6];
// z-max
    sprintf(nameZou, "Z%4.4d", fNBodies++);
    fprintf(fFile1,"%5s XYP%10.3f\n",nameZou, zmax); 
// z-min
    sprintf(nameZin, "Z%4.4d", fNBodies++);
    fprintf(fFile1,"%5s XYP%10.3f\n",nameZin, zmin); 
    
//  outer radius
    d=zmax-zmin;
    sprintf(nameCou, "C%4.4d", fNBodies++);
    if (rmax1 >= 0. && rmax2 >= 0.) {
	if (rmax1!=rmax2) {
	    if (rmax1 > rmax2) {
		pos[2]=zmin;
		r1=rmax1;
		r2=rmax2;
	    } else {
		d=-d;
		pos[2]=zmax;
		r1=rmax2;
		r2=rmax1;
	    }
	    fprintf(fFile1,"%5s TRC%10.3f%10.3f%10.3f%10.3f%10.3f%10.3f\n",
		    nameCou, pos[0], pos[1], pos[2], 0., 0., d); 
	    fprintf(fFile1,"         %10.3f%10.3f\n",r1,r2);
	} else {
	    fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
		    nameCou, pos[0], pos[1], rmax1); 
	} 
    }else {
      fDefaultVolume1 = nameCou;
    }

    
// inner radius
    sprintf(nameCin, "C%4.4d", fNBodies++);
    if (rmin1 >= 0. && rmin2 >= 0.) {
	if (rmin1!=rmin2) {
	    if (rmin1 != 0 && rmin2 !=0) {
		d=zmax-zmin;
		if (rmin1 > rmin2) {
		    pos[2]=zmin;
		    r1=rmin1;
		    r2=rmin2;
		} else {
		    pos[2]=zmax;
		    r1=rmin2;
		    r2=rmin1;
		    d=-d;
		}
		fprintf(fFile1,"%5s TRC%10.3f%10.3f%10.3f%10.3f%10.3f%10.3f\n",
			nameCin, pos[0], pos[1], pos[2], 0., 0., d); 
		fprintf(fFile1,"         %10.3f%10.3f\n",r1,r2);
	    } 
	} else {
	    fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
		    nameCin, pos[0], pos[1], rmin1); 
	} 
    }else {
      fDefaultVolume2 = nameCin;
    }


//
// Volumes
// ^^^^^^^
    fprintf(fFile2,">%s:%s\n", Material, Field);
    fprintf(fFile2,"EMFCUT=%s\n", Cuts);
    fprintf(fFile2,"WW-FACTOR=%s\n", Cuts);
    if (rmin1 != 0 && rmin2 !=0) {
	fprintf(fFile2,"+%s-%s+%s-%s\n", nameCou, nameCin, nameZou, nameZin);
    } else {
	fprintf(fFile2,"+%s+%s-%s\n", nameCou, nameZou, nameZin);
    }
    fprintf(fFile2,"\n");        
}

void AliALIFE::OnionCone (Float_t* r1, Float_t* r2, Int_t nr, 
			  Float_t zmin, Float_t zmax,
			  Float_t pos[3], const char** Materials,
			  const char** Fields, const char** Cuts) 
{
// Concentric cones
//
    char nameCin[6], nameCou[6], nameZin[6], nameZou[6];
//
// Bodies
// ^^^^^^
    Float_t ri1, ri2, ro1, ro2, d;

// z-max
    sprintf(nameZou, "Z%4.4d", fNBodies++);
    fprintf(fFile1,"%5s XYP%10.3f\n",nameZou, zmax); 
// z-min
    sprintf(nameZin, "Z%4.4d", fNBodies++);
    fprintf(fFile1,"%5s XYP%10.3f\n",nameZin, zmin); 
     
//  inner radius
    d=zmax-zmin;
    Bool_t hasInner=kFALSE;
    Bool_t isCylinder=kFALSE;
    if (r1[0]>=0 && r2[0]>=0) {
	if (r1[0] != 0. &&  r2[0] !=0.) {
	    if (r1[0]!=r2[0]) {
		hasInner=kTRUE;
		if (r1[0] > r2[0]) {
		    pos[2]=zmin;
		    ri1=r1[0];
		    ri2=r2[0];
		} else {
		    d=-d;
		    pos[2]=zmax;
		    ri1=r2[0];
		    ri2=r1[0];
		}
		sprintf(nameCin, "C%4.4d", fNBodies++);
		nameCin[5]='\0';
		fprintf(fFile1,"%5s TRC%10.3f%10.3f%10.3f%10.3f%10.3f%10.3f\n",
			nameCin, pos[0], pos[1], pos[2], 0., 0., d); 
		fprintf(fFile1,"         %10.3f%10.3f\n",ri1,ri2);
	    } else {
		isCylinder=kTRUE;
		sprintf(nameCin, "C%4.4d", fNBodies++);
		nameCin[5]='\0';
		fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
			nameCin, pos[0], pos[1], r1[0]); 
	    }
	}
    } else {
      fDefaultVolume1 = nameCin;
    }
    
    
// outer radius
    Int_t i;
    for (i=1; i<nr; i++) {
	if (r1[i] >= 0. && r2[i] >=0) {
	    if (r1[i]!=r2[i]) {
		d=zmax-zmin;
		if (r1[i] > r2[i]) {
		    pos[2]=zmin;
		    ro1=r1[i];
		    ro2=r2[i];
		} else {
		    pos[2]=zmax;
		    ro1=r2[i];
		    ro2=r1[i];
		    d=-d;
		}
		sprintf(nameCou, "C%4.4d", fNBodies++);
		nameCou[5]='\0';
		fprintf(fFile1,"%5s TRC%10.3f%10.3f%10.3f%10.3f%10.3f%10.3f\n",
			nameCou, pos[0], pos[1], pos[2], 0., 0., d); 
		fprintf(fFile1,"         %10.3f%10.3f\n",ro1,ro2);
	    } else {
		isCylinder=kTRUE;
		sprintf(nameCou, "C%4.4d", fNBodies++);
		nameCou[5]='\0';
		fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
			nameCou, pos[0], pos[1], r1[i]); 
	    }
	} else {
	  fDefaultVolume1 = nameCou;
	}
	
// Volumes
// ^^^^^^^
	if (Fields && Cuts) {
	    fprintf(fFile2,">%s:%s\n", Materials[i-1], Fields[i-1]);
	    fprintf(fFile2,"EMFCUT=%s\n", Cuts[i-1]);
	    fprintf(fFile2,"WW-FACTOR=%s\n", Cuts[i-1]);
	} else {
	    fprintf(fFile2,">%s:%s\n", Materials[i-1], "MF");
	    fprintf(fFile2,"EMFCUT=%s\n", "$UNSHIELDED");
	    fprintf(fFile2,"WW-FACTOR=%s\n", "$UNSHIELDED");
	}
	if (hasInner || isCylinder) {
	    if (!isCylinder) fprintf(fFile2,"+%5s-%5s+%5s-%5s\n", 
		    nameCou, nameCin, nameZou, nameZin);
	    if (isCylinder)  fprintf(fFile2,"+%5s+%5s-%5s\n",
		    nameCou, nameZou, nameZin); 
	} else {
	    fprintf(fFile2,"+%5s\n", nameCou);
	    hasInner=kTRUE;
	}
	fprintf(fFile2,"\n");        
	strcpy(nameCin,nameCou);
    }
}


void AliALIFE::PolyCone(Float_t* rmin, Float_t* rmax, Float_t* z, 
			Int_t nz,
			Float_t pos[3], 
			const char* Material,
			const char* Field,
			const char* Cuts) 
{
//
// Equivalent to the Geant3 PCON
    Int_t i;
    for (i=0; i<nz-1; i++) {
// skip (z_i = z_i+1)
	if (z[i] == z[i+1]) continue;
	Cone(rmin[i], rmin[i+1], rmax[i], rmax[i+1], z[i], z[i+1], pos,
	     Material, Field, Cuts);
    }
}

void AliALIFE::OnionPolyCone(Float_t** r, Float_t* z,
			     Int_t nr, Int_t nz,
			     Float_t pos[3], 
			     const char** Materials,
			     const char** Fields,const  char** Cuts)
{
//
// Concentric PCONS
    Int_t i, j;
    for (i=0; i<nz-1; i++) {
// skip (z_i = z_i+1)
	if (z[i] == z[i+1]) continue;
	for (j=0; j<nr-1; j++) {
	    if (Fields && Cuts) {
		Cone(r[i][j], r[i+1][j], r[i][j+1], r[i+1][j+1], 
		     z[i], z[i+1], pos,
		     Materials[i], Fields[i], Cuts[i]);
	    } else {
		Cone(r[i][j], r[i+1][j], r[i][j+1], r[i+1][j+1], 
		     z[i], z[i+1], pos,
		     Materials[i]);
	    }
	}
    }
}

void AliALIFE::Comment(const char* Comment)
{
// Insert comment line
    fprintf(fFile1,"*%s\n", Comment);        
    fprintf(fFile2,"*%s\n", Comment);        
}


void AliALIFE::Finish(Bool_t iremove)
{
// Finish geometry definition
    char s[BUFSIZ];
    fprintf(fFile1,"      END\n");
    fclose(fFile2);
    fFile2=fopen(fVolumeFile,"r");
    while (fgets(s, BUFSIZ, fFile2)) {
	fputs(s,fFile1);
    }
    
    fclose(fFile1);
    fclose(fFile2);    
    if (iremove) {
	remove(fVolumeFile);
	remove(fBodyFile);
    }
}

void AliALIFE::Copy(TObject&) const
{
    //
    // Copy 
    //
    Fatal("Copy","Not implemented!\n");
}
Commit	Line	Data
660eb9a2	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
eeacf08b	16	/* $Id$ */
660eb9a2	17
eeacf08b	18	//----------------------------------------------------------------------
	19	// Helper class to write geometry in ALIFE format
	20	// in parallel with Geant geometry definition.
	21	// Author: A.Morsch
	22	//----------------------------------------------------------------------
88cb7938	23
660eb9a2	24	#include <AliALIFE.h>
	25
	26	ClassImp(AliALIFE)
	27
9b62c596	28
9b62c596	29	AliALIFE::AliALIFE(const char* file1, const char* file2)
660eb9a2	30	{
	31	// Constructor
	32	fNBodies = 0;
	33	fNVolumes= 0;
	34	fBodyFile = file1; // File for Fluka bodies
	35	fVolumeFile = file2;
	36	fFile1=fopen(fBodyFile,"w");
	37	fFile2=fopen(fVolumeFile,"w");
	38	BodyHeader();
	39	VolumeHeader();
	40	fDefaultVolume1 = "DEFAU1";
	41	fDefaultVolume2 = "DEFAU2";
	42	}
	43
9b62c596	44	AliALIFE::AliALIFE()
660eb9a2	45	{
	46	// Default constructor
	47	fBodyFile = "FlukaBody.inp";
	48	fVolumeFile = "FlukaVolume.inp";
	49	fFile1=fopen(fBodyFile,"w");
	50	fFile2=fopen(fVolumeFile,"w");
	51	BodyHeader();
	52	VolumeHeader();
	53	}
	54
	55	void AliALIFE::BodyHeader()
	56	{
	57	// Write header for body definitions
	58	fprintf(fFile1,"*\n");
	59	fprintf(fFile1,"GEOBEGIN COMBINAT\n");
	60	fprintf(fFile1," 0 0 AliRoot Generated\n");
	61	fprintf(fFile1,"*\n");
	62	fprintf(fFile1,"*== Body Definitions ================================================\n");
	63	fprintf(fFile1,"*\n");
	64
	65	}
	66
	67
	68	void AliALIFE::VolumeHeader()
	69	{
	70	// Write header for region (volume) definitions
	71	fprintf(fFile2,"REGIONS\n");
	72	fprintf(fFile2,"*\n");
	73	fprintf(fFile2,"*== Region Definitions =================================================\n");
	74	fprintf(fFile2,"*\n");
	75	}
	76
	77
	78	void AliALIFE:: Cylinder(Float_t rmin, Float_t rmax,
	79	Float_t zmin, Float_t zmax,
	80	Float_t pos[3],
01a62e80	81	const char* Material,
	82	const char* Field,
	83	const char* Cuts)
660eb9a2	84	{
	85	// Simple cylinder
	86	//
	87	// Bodies
	88	// ^^^^^^
	89	char name1[5], name2[5], name3[5], name4[5];
	90
	91	// outer radius
	92	sprintf(name1, "R%4.4d", fNBodies++);
	93	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",name1, pos[0], pos[1], rmax); // inner radius
	94	sprintf(name2, "R%4.4d", fNBodies++);
	95	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",name2, pos[0], pos[1], rmin);
	96	// z-max
	97	sprintf(name3, "Z%4.4d", fNBodies++);
	98	fprintf(fFile1,"%5s XYP%10.3f\n",name3, zmax);
	99	// z-min
	100	sprintf(name4, "Z%4.4d", fNBodies++);
	101	fprintf(fFile1,"%5s XYP%10.3f\n",name4, zmin);
	102	//
	103	// Volumes
	104	// ^^^^^^^
	105
	106	fprintf(fFile2,">%s:%s\n", Material, Field);
	107	fprintf(fFile2,"EMFCUT=%s\n", Cuts);
	108	fprintf(fFile2,"WW-FACTOR=%s\n", Cuts);
	109	if (rmin >0) {
	110	fprintf(fFile2,"+%s-%s+%s-%s\n", name1, name2, name3, name4);
	111	} else {
	112	fprintf(fFile2,"+%s+%s-%s\n", name1, name3, name4);
	113	}
	114
	115	fprintf(fFile2,"\n");
	116	}
	117
	118	void AliALIFE::OnionCylinder(Float_t* r, Int_t nr, Float_t zmin, Float_t zmax,
01a62e80	119	Float_t pos[3], const char** Materials,
01a62e80	120	const char Fields, const char Cuts)
660eb9a2	121	{
	122	//
	123	// Concentric cylinders
	124	//
	125
	126	// Bodies
	127	// ^^^^^^
	128	char nameRin[6], nameRou[6], nameZin[6], nameZou[6];
	129	// z-limits
	130	// z-max
	131	sprintf(nameZou, "Z%4.4d", fNBodies++);
	132	nameZou[5]='\0';
	133
	134	fprintf(fFile1,"%5s XYP%10.3f\n",nameZou, zmax);
	135	// z-min
	136	sprintf(nameZin, "Z%4.4d", fNBodies++);
	137	nameZin[5]='\0';
	138	fprintf(fFile1,"%5s XYP%10.3f\n",nameZin, zmin);
	139
	140	// inner radius
	141	sprintf(nameRin, "R%4.4d", fNBodies++);
	142	nameRin[5]='\0';
	143	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",nameRin, pos[0], pos[1], r[0]);
	144
	145	Int_t i;
	146	for (i=1; i<nr; i++) {
	147	// outer radius
	148	sprintf(nameRou, "R%4.4d", fNBodies++);
	149	nameRou[5]='\0';
	150	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
	151	nameRou, pos[0], pos[1], r[i]);
	152	//
	153	// Volumes
	154	// ^^^^^^^
	155	if (Fields && Cuts) {
	156	fprintf(fFile2,">%s:%s\n", Materials[i-1], Fields[i-1]);
	157	fprintf(fFile2,"EMFCUT=%s\n", Cuts[i-1]);
	158	fprintf(fFile2,"WW-FACTOR=%s\n", Cuts[i-1]);
	159	} else {
	160	fprintf(fFile2,">%s:%s\n", Materials[i-1], "MF");
	161	fprintf(fFile2,"EMFCUT=%s\n", "$UNSHIELDED");
	162	fprintf(fFile2,"WW-FACTOR=%s\n", "$UNSHIELDED");
	163	}
	164	if (r[i-1] != 0.) {
	165	fprintf(fFile2,"+%5s-%5s+%5s-%5s\n", nameZou, nameZin, nameRou, nameRin);
	166	} else {
	167	fprintf(fFile2,"+%5s-%5s+%5s\n", nameZou, nameZin, nameRou);
	168	}
	169	fprintf(fFile2,"\n");
886e839a	170	strcpy(nameRin,nameRou);
660eb9a2	171	}
	172	}
	173
	174
	175	void AliALIFE::Cone(Float_t rmin1, Float_t rmin2,
	176	Float_t rmax1, Float_t rmax2,
	177	Float_t zmin, Float_t zmax,
	178	Float_t pos[3],
01a62e80	179	const char* Material,
	180	const char* Field,
	181	const char* Cuts)
660eb9a2	182	{
	183	// Simple cone
	184
	185	//
	186	// Bodies
	187	// ^^^^^^
	188	char nameCou[6], nameCin[6];
	189	Float_t d, r1, r2;
	190	char nameZin[6], nameZou[6];
	191	// z-max
	192	sprintf(nameZou, "Z%4.4d", fNBodies++);
	193	fprintf(fFile1,"%5s XYP%10.3f\n",nameZou, zmax);
	194	// z-min
	195	sprintf(nameZin, "Z%4.4d", fNBodies++);
	196	fprintf(fFile1,"%5s XYP%10.3f\n",nameZin, zmin);
	197
	198	// outer radius
	199	d=zmax-zmin;
eeacf08b	200	sprintf(nameCou, "C%4.4d", fNBodies++);
660eb9a2	201	if (rmax1 >= 0. && rmax2 >= 0.) {
	202	if (rmax1!=rmax2) {
	203	if (rmax1 > rmax2) {
	204	pos[2]=zmin;
	205	r1=rmax1;
	206	r2=rmax2;
	207	} else {
	208	d=-d;
	209	pos[2]=zmax;
	210	r1=rmax2;
	211	r2=rmax1;
	212	}
660eb9a2	213	fprintf(fFile1,"%5s TRC%10.3f%10.3f%10.3f%10.3f%10.3f%10.3f\n",
	214	nameCou, pos[0], pos[1], pos[2], 0., 0., d);
	215	fprintf(fFile1," %10.3f%10.3f\n",r1,r2);
	216	} else {
660eb9a2	217	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
	218	nameCou, pos[0], pos[1], rmax1);
	219	}
	220	}else {
b7fab4b5	221	fDefaultVolume1 = nameCou;
660eb9a2	222	}
	223
	224
	225
	226	// inner radius
eeacf08b	227	sprintf(nameCin, "C%4.4d", fNBodies++);
660eb9a2	228	if (rmin1 >= 0. && rmin2 >= 0.) {
	229	if (rmin1!=rmin2) {
	230	if (rmin1 != 0 && rmin2 !=0) {
	231	d=zmax-zmin;
	232	if (rmin1 > rmin2) {
	233	pos[2]=zmin;
	234	r1=rmin1;
	235	r2=rmin2;
	236	} else {
	237	pos[2]=zmax;
	238	r1=rmin2;
	239	r2=rmin1;
	240	d=-d;
	241	}
660eb9a2	242	fprintf(fFile1,"%5s TRC%10.3f%10.3f%10.3f%10.3f%10.3f%10.3f\n",
	243	nameCin, pos[0], pos[1], pos[2], 0., 0., d);
	244	fprintf(fFile1," %10.3f%10.3f\n",r1,r2);
	245	}
	246	} else {
660eb9a2	247	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
	248	nameCin, pos[0], pos[1], rmin1);
	249	}
	250	}else {
b7fab4b5	251	fDefaultVolume2 = nameCin;
660eb9a2	252	}
	253
	254
	255
	256	//
	257	// Volumes
	258	// ^^^^^^^
	259	fprintf(fFile2,">%s:%s\n", Material, Field);
	260	fprintf(fFile2,"EMFCUT=%s\n", Cuts);
	261	fprintf(fFile2,"WW-FACTOR=%s\n", Cuts);
	262	if (rmin1 != 0 && rmin2 !=0) {
	263	fprintf(fFile2,"+%s-%s+%s-%s\n", nameCou, nameCin, nameZou, nameZin);
	264	} else {
	265	fprintf(fFile2,"+%s+%s-%s\n", nameCou, nameZou, nameZin);
	266	}
	267	fprintf(fFile2,"\n");
	268	}
	269
	270	void AliALIFE::OnionCone (Float_t* r1, Float_t* r2, Int_t nr,
	271	Float_t zmin, Float_t zmax,
01a62e80	272	Float_t pos[3], const char** Materials,
01a62e80	273	const char Fields, const char Cuts)
660eb9a2	274	{
	275	// Concentric cones
	276	//
	277	char nameCin[6], nameCou[6], nameZin[6], nameZou[6];
	278	//
	279	// Bodies
	280	// ^^^^^^
	281	Float_t ri1, ri2, ro1, ro2, d;
	282
	283	// z-max
	284	sprintf(nameZou, "Z%4.4d", fNBodies++);
	285	fprintf(fFile1,"%5s XYP%10.3f\n",nameZou, zmax);
	286	// z-min
	287	sprintf(nameZin, "Z%4.4d", fNBodies++);
	288	fprintf(fFile1,"%5s XYP%10.3f\n",nameZin, zmin);
	289
	290	// inner radius
	291	d=zmax-zmin;
	292	Bool_t hasInner=kFALSE;
	293	Bool_t isCylinder=kFALSE;
	294	if (r1[0]>=0 && r2[0]>=0) {
	295	if (r1[0] != 0. && r2[0] !=0.) {
	296	if (r1[0]!=r2[0]) {
	297	hasInner=kTRUE;
	298	if (r1[0] > r2[0]) {
	299	pos[2]=zmin;
	300	ri1=r1[0];
	301	ri2=r2[0];
	302	} else {
	303	d=-d;
	304	pos[2]=zmax;
	305	ri1=r2[0];
	306	ri2=r1[0];
	307	}
	308	sprintf(nameCin, "C%4.4d", fNBodies++);
	309	nameCin[5]='\0';
	310	fprintf(fFile1,"%5s TRC%10.3f%10.3f%10.3f%10.3f%10.3f%10.3f\n",
	311	nameCin, pos[0], pos[1], pos[2], 0., 0., d);
	312	fprintf(fFile1," %10.3f%10.3f\n",ri1,ri2);
	313	} else {
	314	isCylinder=kTRUE;
	315	sprintf(nameCin, "C%4.4d", fNBodies++);
	316	nameCin[5]='\0';
	317	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
	318	nameCin, pos[0], pos[1], r1[0]);
	319	}
	320	}
	321	} else {
b7fab4b5	322	fDefaultVolume1 = nameCin;
660eb9a2	323	}
	324
	325
	326
	327	// outer radius
	328	Int_t i;
	329	for (i=1; i<nr; i++) {
	330	if (r1[i] >= 0. && r2[i] >=0) {
	331	if (r1[i]!=r2[i]) {
	332	d=zmax-zmin;
	333	if (r1[i] > r2[i]) {
	334	pos[2]=zmin;
	335	ro1=r1[i];
	336	ro2=r2[i];
	337	} else {
	338	pos[2]=zmax;
	339	ro1=r2[i];
	340	ro2=r1[i];
	341	d=-d;
	342	}
	343	sprintf(nameCou, "C%4.4d", fNBodies++);
	344	nameCou[5]='\0';
	345	fprintf(fFile1,"%5s TRC%10.3f%10.3f%10.3f%10.3f%10.3f%10.3f\n",
	346	nameCou, pos[0], pos[1], pos[2], 0., 0., d);
	347	fprintf(fFile1," %10.3f%10.3f\n",ro1,ro2);
	348	} else {
	349	isCylinder=kTRUE;
	350	sprintf(nameCou, "C%4.4d", fNBodies++);
	351	nameCou[5]='\0';
	352	fprintf(fFile1,"%5s ZCC%10.3f%10.3f%10.3f\n",
	353	nameCou, pos[0], pos[1], r1[i]);
	354	}
	355	} else {
b7fab4b5	356	fDefaultVolume1 = nameCou;
660eb9a2	357	}
	358
	359	// Volumes
	360	// ^^^^^^^
	361	if (Fields && Cuts) {
	362	fprintf(fFile2,">%s:%s\n", Materials[i-1], Fields[i-1]);
	363	fprintf(fFile2,"EMFCUT=%s\n", Cuts[i-1]);
	364	fprintf(fFile2,"WW-FACTOR=%s\n", Cuts[i-1]);
	365	} else {
	366	fprintf(fFile2,">%s:%s\n", Materials[i-1], "MF");
	367	fprintf(fFile2,"EMFCUT=%s\n", "$UNSHIELDED");
	368	fprintf(fFile2,"WW-FACTOR=%s\n", "$UNSHIELDED");
	369	}
38e826ad	370	if (hasInner \|\| isCylinder) {
38e826ad	371	if (!isCylinder) fprintf(fFile2,"+%5s-%5s+%5s-%5s\n",
660eb9a2	372	nameCou, nameCin, nameZou, nameZin);
38e826ad	373	if (isCylinder) fprintf(fFile2,"+%5s+%5s-%5s\n",
38e826ad	374	nameCou, nameZou, nameZin);
660eb9a2	375	} else {
	376	fprintf(fFile2,"+%5s\n", nameCou);
	377	hasInner=kTRUE;
	378	}
	379	fprintf(fFile2,"\n");
886e839a	380	strcpy(nameCin,nameCou);
660eb9a2	381	}
	382	}
	383
	384
	385	void AliALIFE::PolyCone(Float_t* rmin, Float_t* rmax, Float_t* z,
	386	Int_t nz,
	387	Float_t pos[3],
01a62e80	388	const char* Material,
	389	const char* Field,
	390	const char* Cuts)
660eb9a2	391	{
	392	//
	393	// Equivalent to the Geant3 PCON
	394	Int_t i;
	395	for (i=0; i<nz-1; i++) {
	396	// skip (z_i = z_i+1)
	397	if (z[i] == z[i+1]) continue;
	398	Cone(rmin[i], rmin[i+1], rmax[i], rmax[i+1], z[i], z[i+1], pos,
	399	Material, Field, Cuts);
	400	}
	401	}
	402
	403	void AliALIFE::OnionPolyCone(Float_t** r, Float_t* z,
	404	Int_t nr, Int_t nz,
	405	Float_t pos[3],
01a62e80	406	const char** Materials,
01a62e80	407	const char Fields,const char Cuts)
660eb9a2	408	{
	409	//
	410	// Concentric PCONS
	411	Int_t i, j;
	412	for (i=0; i<nz-1; i++) {
	413	// skip (z_i = z_i+1)
	414	if (z[i] == z[i+1]) continue;
	415	for (j=0; j<nr-1; j++) {
	416	if (Fields && Cuts) {
	417	Cone(r[i][j], r[i+1][j], r[i][j+1], r[i+1][j+1],
	418	z[i], z[i+1], pos,
	419	Materials[i], Fields[i], Cuts[i]);
	420	} else {
	421	Cone(r[i][j], r[i+1][j], r[i][j+1], r[i+1][j+1],
	422	z[i], z[i+1], pos,
	423	Materials[i]);
	424	}
	425	}
	426	}
	427	}
	428
01a62e80	429	void AliALIFE::Comment(const char* Comment)
660eb9a2	430	{
	431	// Insert comment line
	432	fprintf(fFile1,"*%s\n", Comment);
	433	fprintf(fFile2,"*%s\n", Comment);
	434	}
	435
	436
e59713af	437	void AliALIFE::Finish(Bool_t iremove)
660eb9a2	438	{
	439	// Finish geometry definition
	440	char s[BUFSIZ];
	441	fprintf(fFile1," END\n");
	442	fclose(fFile2);
	443	fFile2=fopen(fVolumeFile,"r");
	444	while (fgets(s, BUFSIZ, fFile2)) {
	445	fputs(s,fFile1);
	446	}
	447
	448	fclose(fFile1);
	449	fclose(fFile2);
e59713af	450	if (iremove) {
	451	remove(fVolumeFile);
	452	remove(fBodyFile);
	453	}
660eb9a2	454	}
660eb9a2	455
cbd70aad	456	void AliALIFE::Copy(TObject&) const
9b62c596	457	{
	458	//
	459	// Copy
	460	//
	461	Fatal("Copy","Not implemented!\n");
	462	}
660eb9a2	463
	464
	465
	466
	467
	468
	469
	470
	471
	472