**************************************************************************/
/*
$Log$
+Revision 1.11 2001/05/14 14:01:04 morsch
+AliPMDv0 coarse geometry and AliPMDv1 detailed simulation, completely revised versions by Tapan Nayak.
*/
//
z_air1= -dpara_sm1[2] + th_base + dpara_air[2];
gMC->Gspos("EAIR", 1, "ESMA", 0., 0., z_air1, 0, "ONLY");
z_gas=z_air1+dpara_air[2]+ th_pcb + dpara1[2];
- gMC->Gspos("EHC1", 1, "ESMA", 0., 0., z_gas, 0, "ONLY");
+ //Line below Commented for version 0 of PMD routine
+ // gMC->Gspos("EHC1", 1, "ESMA", 0., 0., z_gas, 0, "ONLY");
z_air2=z_gas+dpara1[2]+ th_pcb + dpara_air[2];
gMC->Gspos("EAIR", 2, "ESMA", 0., 0., z_air2, 0, "ONLY");
z_air1= -dpara_sm1[2] + th_pcb + dpara_air[2];
gMC->Gspos("EAIR", 3, "ESMB", 0., 0., z_air1, 0, "ONLY");
z_gas=z_air1+dpara_air[2]+ th_pcb + dpara1[2];
- gMC->Gspos("EHC1", 2, "ESMB", 0., 0., z_gas, 0, "ONLY");
+ //Line below Commented for version 0 of PMD routine
+ // gMC->Gspos("EHC1", 2, "ESMB", 0., 0., z_gas, 0, "ONLY");
z_air2=z_gas+dpara1[2]+ th_pcb + dpara_air[2];
gMC->Gspos("EAIR", 4, "ESMB", 0., 0., z_air2, 0, "ONLY");
z_air1= -dpara_sm2[2] + th_base + dpara2_air[2];
gMC->Gspos("EAIX", 1, "ESMX", 0., 0., z_air1, 0, "ONLY");
z_gas=z_air1+dpara2_air[2]+ th_pcb + dpara2[2];
- gMC->Gspos("EHC2", 1, "ESMX", 0., 0., z_gas, 0, "ONLY");
+ //Line below Commented for version 0 of PMD routine
+ // gMC->Gspos("EHC2", 1, "ESMX", 0., 0., z_gas, 0, "ONLY");
z_air2=z_gas+dpara2[2]+ th_pcb + dpara2_air[2];
gMC->Gspos("EAIX", 2, "ESMX", 0., 0., z_air2, 0, "ONLY");
z_air1= -dpara_sm2[2] + th_pcb + dpara2_air[2];
gMC->Gspos("EAIX", 3, "ESMY", 0., 0., z_air1, 0, "ONLY");
z_gas=z_air1+dpara2_air[2]+ th_pcb + dpara2[2];
- gMC->Gspos("EHC2", 2, "ESMY", 0., 0., z_gas, 0, "ONLY");
+ //Line below Commented for version 0 of PMD routine
+ // gMC->Gspos("EHC2", 2, "ESMY", 0., 0., z_gas, 0, "ONLY");
z_air2=z_gas+dpara2[2]+ th_pcb + dpara2_air[2];
gMC->Gspos("EAIX", 4, "ESMY", 0., 0., z_air2, 0, "ONLY");
z_air1= -dpara_sm3[2] + th_base + dpara3_air[2];
gMC->Gspos("EAIP", 1, "ESMP", 0., 0., z_air1, 0, "ONLY");
z_gas=z_air1+dpara3_air[2]+ th_pcb + dpara3[2];
- gMC->Gspos("EHC3", 1, "ESMP", 0., 0., z_gas, 0, "ONLY");
+ //Line below Commented for version 0 of PMD routine
+ // gMC->Gspos("EHC3", 1, "ESMP", 0., 0., z_gas, 0, "ONLY");
z_air2=z_gas+dpara3[2]+ th_pcb + dpara3_air[2];
gMC->Gspos("EAIP", 2, "ESMP", 0., 0., z_air2, 0, "ONLY");
z_air1= -dpara_sm3[2] + th_pcb + dpara3_air[2];
gMC->Gspos("EAIP", 3, "ESMQ", 0., 0., z_air1, 0, "ONLY");
z_gas=z_air1+dpara3_air[2]+ th_pcb + dpara3[2];
- gMC->Gspos("EHC3", 2, "ESMQ", 0., 0., z_gas, 0, "ONLY");
+ //Line below Commented for version 0 of PMD routine
+ // gMC->Gspos("EHC3", 2, "ESMQ", 0., 0., z_gas, 0, "ONLY");
z_air2=z_gas+dpara3[2]+ th_pcb + dpara3_air[2];
gMC->Gspos("EAIP", 4, "ESMQ", 0., 0., z_air2, 0, "ONLY");
Float_t Al_rod[3] ;
Al_rod[0] = sm_length * 3/2. - gaspmd[5]/2 - boundary ;
- Al_rod[1] = boundary;
+ Al_rod[1] = boundary - 0.5*cell_radius*root3_2;
Al_rod[2] = dm_thick/2.;
gMC->Gsvolu("EALM","BOX ", idtmed[698], Al_rod, 3);
xpos[i]=xmod[j] + xsup[i]*TMath::Cos(theta[j]) - ysup[i]*TMath::Sin(theta[j]);
ypos[i]=ymod[j] + xsup[i]*TMath::Sin(theta[j]) + ysup[i]*TMath::Cos(theta[j]);
- if(fDebug)
- printf("%s: %f %f \n", ClassName(), xpos[i], ypos[i]);
+ printf("%f %f \n", xpos[i], ypos[i]);
num_mod = num_mod+1;
- if(fDebug)
- printf("\n%s: Num_mod %d\n",ClassName(),num_mod);
+ printf("\nNum_mod %d\n",num_mod);
gMC->Gsposp("EMM1", num_mod + 6, "EPMD", xpos[i],ypos[i], 0., irotate[j], "ONLY", dpara_emm1, 6);
Int_t i;
kdet=1;
//
- if(fDebug) {
- printf("\n%s: ",ClassName());
- for(i=0;i<35;i++) printf("*");
- printf(" PMD_INIT ");
- for(i=0;i<35;i++) printf("*");
- printf("\n%s: ",ClassName());
- printf(" PMD simulation package (v1) initialised\n");
- printf("%s: parameters of pmd\n",ClassName());
- printf("%s: %10.2f %10.2f %10.2f %10.2f\n",ClassName(),cell_radius,cell_wall,cell_depth,zdist1 );
- printf("%s: ",ClassName());
- for(i=0;i<80;i++) printf("*");
- printf("\n");
- }
+ printf("\n");
+ for(i=0;i<35;i++) printf("*");
+ printf(" PMD_INIT ");
+ for(i=0;i<35;i++) printf("*");
+ printf("\n");
+ printf(" PMD simulation package (v1) initialised\n");
+ printf(" parameters of pmd\n");
+ printf("%10.2f %10.2f %10.2f %10.2f\n", cell_radius,cell_wall,cell_depth,zdist1 );
+
+ for(i=0;i<80;i++) printf("*");
+ printf("\n");
Int_t *idtmed = fIdtmed->GetArray()-599;
fMedSens=idtmed[605-1];