statusOld = -1;
}
Int_t igstatus = 0;
- //------------------modified by Mriganka ----------------------
+
+ Int_t trnotemp = trackno; // Modified on 25th Nov 2009
if(ks==1||(imo = mparticle->GetFirstMother())<0 ){
vx = mparticle->Vx();
vy = mparticle->Vy();
vy = mparticle->Vy();
vz = mparticle->Vz();
- trackno=imo;
-
+ // Modified on 25th Nov 2009
+
+ trnotemp = trackno;
+ if(trackpid == 111)
+ {
+ trackno = trnotemp;
+ }
+ if(trackpid != 111)
+ {
+ trackno=imo;
}
+ // end of modification on 25th Nov 2009
+ }
if(trackpid==kGamma||trackpid==11||trackpid==-11||
trackpid==kPi0)igstatus=1;
// PMD
fDetNo = 0;
}
- //Int_t smn = smnumber - 1;
+
Int_t smn = smnumber;
Int_t ixx = xpad - 1;
Int_t iyy = ypad - 1;
}
Int_t igstatus = 0;
- //-----------------------modified by Mriganka ------------------
+
+ Int_t trnotemp = trackno; // modified on 25th Nov 2009
if(ks==1||(imo = mparticle->GetFirstMother())<0 ){
vx = mparticle->Vx();
vy = mparticle->Vy();
vy = mparticle->Vy();
vz = mparticle->Vz();
- trackno=imo;
-
+ // Modified on 25th Nov 2009
+
+ trnotemp = trackno;
+ if(trackpid == 111)
+ {
+ trackno = trnotemp;
+ }
+ if(trackpid != 111)
+ {
+ trackno=imo;
}
+ }
if(trackpid==kGamma||trackpid==11||trackpid==-11||trackpid==kPi0)
igstatus=1;
trackpid=trackpidOld;
trackno=tracknoOld;
- //-----------------end of modification----------------
xPos = fPMDHit->X();
yPos = fPMDHit->Y();
zPos = fPMDHit->Z();
}
AliDebug(2,Form("ZPosition = %f Edeposition = %d",zPos,edep));
- //Float_t zposition = TMath::Abs(zPos);
if (zPos < fZPos)
{
fDetNo = 0;
}
- //Int_t smn = smnumber - 1;
Int_t smn = smnumber;
Int_t ixx = xpad - 1;
Int_t iyy = ypad - 1;
{
// This converts the simulated edep to ADC according to the
// Test Beam Data
- //PS Test in May 2009
+ //PS Test in May 2009, Voltage @ 1350 V
// KeV - ADC conversion for 12bit ADC
- // Data provided by Ajay and agreed by Dr. Viyogi
// MPV data used for the fit and taken here
- const Float_t kConstant = -9.22;
- const Float_t kErConstant = 0.602225;
- const Float_t kSlope = 79.1119;
- const Float_t kErSlope = 0.844047;
+ const Float_t kConstant = -0.1602;
+ const Float_t kErConstant = 0.9914;
+ const Float_t kSlope = 77.47;
+ const Float_t kErSlope = 3.16;
Float_t cons = gRandom->Gaus(kConstant,kErConstant);
Float_t slop = gRandom->Gaus(kSlope,kErSlope);