Int_t vol[5];
//char *namep;
- if(gMC->GetMedium() == fMedSens && (destep = gMC->Edep())) {
+ if(gMC->CurrentMedium() == fMedSens && (destep = gMC->Edep())) {
gMC->CurrentVolID(copy);
Int_t vol[8];
//const char *namep;
- if(gMC->GetMedium() == fMedSens && (destep = gMC->Edep())) {
+ if(gMC->CurrentMedium() == fMedSens && (destep = gMC->Edep())) {
gMC->CurrentVolID(copy);
//namep=gMC->CurrentVolName();
Int_t mid=gMC->CurrentMaterial(a,z,den,rad,abs);
Info("Material","id=%i a=%7.2f z=%7.2f den=%9.4f rad=%9.2f abs=%9.2f",mid,a,z,den,rad,abs);
- Int_t iTmedId=gMC->GetMedium();
+ Int_t iTmedId=gMC->CurrentMedium();
const char *sTmed;
switch(iTmedId){
case kAir: sTmed="Air" ;break;
Int_t *idtmed = fIdtmed->GetArray()-499;
Float_t incidenceAngle;
- if(gMC->GetMedium()==idtmed[513] &&
+ if(gMC->CurrentMedium()==idtmed[513] &&
gMC->IsTrackEntering() && gMC->TrackCharge()
&& gMC->CurrentVolID(copy)==fIdSens)
{
Int_t *idtmed = fIdtmed->GetArray()-499;
Float_t incidenceAngle;
- if(gMC->GetMedium()==idtmed[513] &&
+ if(gMC->CurrentMedium()==idtmed[513] &&
gMC->IsTrackEntering() && gMC->TrackCharge()
&& gMC->CurrentVolID(copy)==fIdSens)
{
Int_t *idtmed = fIdtmed->GetArray()-499;
Float_t incidenceAngle;
- if(gMC->GetMedium()==idtmed[513] &&
+ if(gMC->CurrentMedium()==idtmed[513] &&
gMC->IsTrackEntering() && gMC->TrackCharge()
&& gMC->CurrentVolID(copy)==fIdSens)
{
Int_t *idtmed = fIdtmed->GetArray()-499;
Float_t incidenceAngle;
- if(gMC->GetMedium()==idtmed[513] &&
+ if(gMC->CurrentMedium()==idtmed[513] &&
gMC->IsTrackEntering() && gMC->TrackCharge()
&& gMC->CurrentVolID(copy)==fIdSens)
{
Int_t *idtmed = fIdtmed->GetArray()-499;
Float_t incidenceAngle;
- if(gMC->GetMedium()==idtmed[513] &&
+ if(gMC->CurrentMedium()==idtmed[513] &&
gMC->IsTrackEntering() && gMC->TrackCharge()
&& gMC->CurrentVolID(copy)==fIdSens)
{
Int_t *idtmed = fIdtmed->GetArray()-499;
Float_t incidenceAngle;
- if(gMC->GetMedium()==idtmed[513] &&
+ if(gMC->CurrentMedium()==idtmed[513] &&
gMC->IsTrackEntering() && gMC->TrackCharge()
&& gMC->CurrentVolID(copy)==fIdSens)
{
for (j=0;j<10;j++) hits[j]=0;
- if((gMC->GetMedium() == fMedSensZN) || (gMC->GetMedium() == fMedSensZP) ||
- (gMC->GetMedium() == fMedSensGR) || (gMC->GetMedium() == fMedSensF1) ||
- (gMC->GetMedium() == fMedSensF2) || (gMC->GetMedium() == fMedSensZEM)){
+ if((gMC->CurrentMedium() == fMedSensZN) || (gMC->CurrentMedium() == fMedSensZP) ||
+ (gMC->CurrentMedium() == fMedSensGR) || (gMC->CurrentMedium() == fMedSensF1) ||
+ (gMC->CurrentMedium() == fMedSensF2) || (gMC->CurrentMedium() == fMedSensZEM)){
// --- This part is for no shower developement in beam pipe and TDI
-// (gMC->GetMedium() == fMedSensPI) || (gMC->GetMedium() == fMedSensTDI)){
+// (gMC->CurrentMedium() == fMedSensPI) || (gMC->CurrentMedium() == fMedSensTDI)){
// If particle interacts with beam pipe -> return
-// if((gMC->GetMedium() == fMedSensPI) || (gMC->GetMedium() == fMedSensTDI)){
+// if((gMC->CurrentMedium() == fMedSensPI) || (gMC->CurrentMedium() == fMedSensTDI)){
// If option NoShower is set -> StopTrack
// if(fNoShower==1) {
-// if(gMC->GetMedium() == fMedSensPI) {
+// if(gMC->CurrentMedium() == fMedSensPI) {
// knamed = gMC->CurrentVolName();
// if((!strncmp(knamed,"MQ",2)) || (!strncmp(knamed,"YM",2))) fpLostIT += 1;
// if((!strncmp(knamed,"MD1",3))|| (!strncmp(knamed,"YD1",2))) fpLostD1 += 1;
// }
-// if(gMC->GetMedium() == fMedSensTDI) fpLostTDI += 1;
+// if(gMC->CurrentMedium() == fMedSensTDI) fpLostTDI += 1;
// gMC->StopTrack();
// printf("\n # of p lost in Inner Triplet = %d\n",fpLostIT);
// printf("\n # of p lost in D1 = %d\n",fpLostD1);
// *** Light production in fibres
- if((gMC->GetMedium() == fMedSensF1) || (gMC->GetMedium() == fMedSensF2)){
+ if((gMC->CurrentMedium() == fMedSensF1) || (gMC->CurrentMedium() == fMedSensF2)){
//Select charged particles
if((destep=gMC->Edep())){
nphe = gRandom->Poisson(out);
// printf("ZN --- ibeta = %d, ialfa = %d, ibe = %d"
// " -> out = %f, nphe = %d\n", ibeta, ialfa, ibe, out, nphe);
- if(gMC->GetMedium() == fMedSensF1){
+ if(gMC->CurrentMedium() == fMedSensF1){
hits[7] = nphe; //fLightPMQ
hits[8] = 0;
hits[9] = 0;
nphe = gRandom->Poisson(out);
// printf("ZP --- ibeta = %d, ialfa = %d, ibe = %d"
// " -> out = %f, nphe = %d\n", ibeta, ialfa, ibe, out, nphe);
- if(gMC->GetMedium() == fMedSensF1){
+ if(gMC->CurrentMedium() == fMedSensF1){
hits[7] = nphe; //fLightPMQ
hits[8] = 0;
hits[9] = 0;
// --- This part is for no shower developement in beam pipe and TDI
// If particle interacts with beam pipe or TDI -> return
- if((gMC->GetMedium() == fMedSensPI) || (gMC->GetMedium() == fMedSensTDI)){
+ if((gMC->CurrentMedium() == fMedSensPI) || (gMC->CurrentMedium() == fMedSensTDI)){
// If option NoShower is set -> StopTrack
if(fNoShower==1) {
- if(gMC->GetMedium() == fMedSensPI) {
+ if(gMC->CurrentMedium() == fMedSensPI) {
knamed = gMC->CurrentVolName();
if(!strncmp(knamed,"YMQ",3)) fpLostIT += 1;
if(!strncmp(knamed,"YD1",3)) fpLostD1 += 1;
}
- else if(gMC->GetMedium() == fMedSensTDI){ // NB->Cu = TDI or D1 vacuum chamber
+ else if(gMC->CurrentMedium() == fMedSensTDI){ // NB->Cu = TDI or D1 vacuum chamber
knamed = gMC->CurrentVolName();
if(!strncmp(knamed,"MD1",3)) fpLostD1 += 1;
if(!strncmp(knamed,"QTD",3)) fpLostTDI += 1;
return;
}
- if((gMC->GetMedium() == fMedSensZN) || (gMC->GetMedium() == fMedSensZP) ||
- (gMC->GetMedium() == fMedSensGR) || (gMC->GetMedium() == fMedSensF1) ||
- (gMC->GetMedium() == fMedSensF2) || (gMC->GetMedium() == fMedSensZEM)){
+ if((gMC->CurrentMedium() == fMedSensZN) || (gMC->CurrentMedium() == fMedSensZP) ||
+ (gMC->CurrentMedium() == fMedSensGR) || (gMC->CurrentMedium() == fMedSensF1) ||
+ (gMC->CurrentMedium() == fMedSensF2) || (gMC->CurrentMedium() == fMedSensZEM)){
//Particle coordinates
// *** Light production in fibres
- if((gMC->GetMedium() == fMedSensF1) || (gMC->GetMedium() == fMedSensF2)){
+ if((gMC->CurrentMedium() == fMedSensF1) || (gMC->CurrentMedium() == fMedSensF2)){
//Select charged particles
if((destep=gMC->Edep())){
// Ch. debug
//if(ibeta==3) printf("\t %f \t %f \t %f\n",alfa, be, out);
//printf("\t ibeta = %d, ialfa = %d, ibe = %d -> nphe = %d\n\n",ibeta,ialfa,ibe,nphe);
- if(gMC->GetMedium() == fMedSensF1){
+ if(gMC->CurrentMedium() == fMedSensF1){
hits[7] = nphe; //fLightPMQ
hits[8] = 0;
hits[9] = 0;
if(ibe>fNbep) ibe=fNbep;
out = charge*charge*fTablep[ibeta][ialfa][ibe];
nphe = gRandom->Poisson(out);
- if(gMC->GetMedium() == fMedSensF1){
+ if(gMC->CurrentMedium() == fMedSensF1){
hits[7] = nphe; //fLightPMQ
hits[8] = 0;
hits[9] = 0;