+ // 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) ;
+
+ AliEMCALGeometry* geom = GetGeometry();
+ if(geom->GetILOSS()>=0) {
+ for(int i=1600; i<=1603; i++) gMC->Gstpar(idtmed[i], "LOSS", geom->GetILOSS()) ;
+ }
+ if(geom->GetIHADR()>=0) {
+ for(int i=1600; i<=1603; i++) gMC->Gstpar(idtmed[i], "HADR", geom->GetIHADR()) ;
+ }
+}
+
+//____________________________________________________________________________
+AliDigitizer* AliEMCAL::CreateDigitizer(AliRunDigitizer* manager) const
+{
+ //create and return the digitizer
+ return new AliEMCALDigitizer(manager);