**************************************************************************/
/*
$Log$
+Revision 1.13 2001/05/21 10:59:09 morsch
+Printouts in debug mode only.
+
+Revision 1.12 2001/05/21 09:39:28 morsch
+Minor modifications on the geometry. (Tapan Nayak)
+
+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");
{
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]);
-
- printf("%f %f \n", xpos[i], ypos[i]);
+ if(fDebug)
+ printf("%s: %f %f \n", ClassName(), xpos[i], ypos[i]);
num_mod = num_mod+1;
- printf("\nNum_mod %d\n",num_mod);
+ if(fDebug)
+ printf("\n%s: Num_mod %d\n",ClassName(),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;
//
- 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 (v0) 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");
-
+ 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 (v0) 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");
+ }
Int_t *idtmed = fIdtmed->GetArray()-599;
fMedSens=idtmed[605-1];
}