-// --- Set decent energy thresholds for gamma and electron tracking
-
- // Tracking threshold for photons and electrons in Lead
- Float_t cutgam=10.e-5; // 100 kev;
- Float_t cutele=10.e-5; // 100 kev;
- TString ntmp(GetTitle());
- ntmp.ToUpper();
- if(ntmp.Contains("10KEV")) {
- cutele = cutgam = 1.e-5;
- } else if(ntmp.Contains("50KEV")) {
- cutele = cutgam = 5.e-5;
- } else if(ntmp.Contains("100KEV")) {
- cutele = cutgam = 1.e-4;
- } else if(ntmp.Contains("200KEV")) {
- cutele = cutgam = 2.e-4;
- } else if(ntmp.Contains("500KEV")) {
- cutele = cutgam = 5.e-4;
- }
-
- gMC->Gstpar(idtmed[1600],"CUTGAM", cutgam);
- gMC->Gstpar(idtmed[1600],"CUTELE", cutele); // 1MEV -> 0.1MEV; 15-aug-05
- gMC->Gstpar(idtmed[1600],"BCUTE", cutgam); // BCUTE and BCUTM start from GUTGUM
- gMC->Gstpar(idtmed[1600],"BCUTM", cutgam); // BCUTE and BCUTM start from GUTGUM
- // --- Generate explicitly delta rays in Lead ---
- gMC->Gstpar(idtmed[1600], "LOSS",3.) ;
- gMC->Gstpar(idtmed[1600], "DRAY",1.) ;
- gMC->Gstpar(idtmed[1600], "DCUTE", cutele) ;
- gMC->Gstpar(idtmed[1600], "DCUTM", cutele) ;
-
-// --- in aluminium parts ---
- gMC->Gstpar(idtmed[1602],"CUTGAM", cutgam) ;
- gMC->Gstpar(idtmed[1602],"CUTELE", cutele) ;
- gMC->Gstpar(idtmed[1602],"BCUTE", cutgam); // BCUTE and BCUTM start from GUTGUM
- gMC->Gstpar(idtmed[1602],"BCUTM", cutgam); // BCUTE and BCUTM start from GUTGUM
- gMC->Gstpar(idtmed[1602], "LOSS",3.) ;
- gMC->Gstpar(idtmed[1602], "DRAY",1.) ;
- gMC->Gstpar(idtmed[1602], "DCUTE", cutele) ;
- gMC->Gstpar(idtmed[1602], "DCUTM", cutele) ;
-
-// --- and finally thresholds for photons and electrons in the scintillator ---
- gMC->Gstpar(idtmed[1601],"CUTGAM", cutgam) ;
- gMC->Gstpar(idtmed[1601],"CUTELE", cutele) ;// 1MEV -> 0.1MEV; 15-aug-05
- gMC->Gstpar(idtmed[1601],"BCUTE", cutgam); // BCUTE and BCUTM start from GUTGUM
- gMC->Gstpar(idtmed[1601],"BCUTM", cutgam); // BCUTE and BCUTM start from GUTGUM
- gMC->Gstpar(idtmed[1601], "LOSS",3.) ; // generate delta rays
- gMC->Gstpar(idtmed[1601], "DRAY",1.) ;
- gMC->Gstpar(idtmed[1601], "DCUTE", cutele) ;
- gMC->Gstpar(idtmed[1601], "DCUTM", cutele) ;
-
- // S steel -
- gMC->Gstpar(idtmed[1603],"CUTGAM", cutgam);
- gMC->Gstpar(idtmed[1603],"CUTELE", cutele);
- gMC->Gstpar(idtmed[1603],"BCUTE", cutgam); // BCUTE and BCUTM start from GUTGUM
- gMC->Gstpar(idtmed[1603],"BCUTM", cutgam); // BCUTE and BCUTM start from GUTGUM
- // --- Generate explicitly delta rays
- gMC->Gstpar(idtmed[1603], "LOSS",3.);
- gMC->Gstpar(idtmed[1603], "DRAY",1.);
- gMC->Gstpar(idtmed[1603], "DCUTE", cutele) ;
- gMC->Gstpar(idtmed[1603], "DCUTM", cutele) ;