statusOld = -1;
}
Int_t igstatus = 0;
- Int_t trnotemp = trackno;
//------------------modified by Mriganka ----------------------
if(ks==1||(imo = mparticle->GetFirstMother())<0 ){
vx = mparticle->Vx();
while(((imo = mparticle->GetFirstMother()) >= 0)&&
(ks = mparticle->GetStatusCode() <1) )
- {
+ {
mparticle = gAlice->GetMCApp()->Particle(imo);
trackpid = mparticle->GetPdgCode();
ks = mparticle->GetStatusCode();
vx = mparticle->Vx();
vy = mparticle->Vy();
vz = mparticle->Vz();
- trnotemp = trackno;
- if(trackpid == 111)
- {
- trackno = trnotemp;
- }
- if(trackpid != 111)
- {
- trackno = imo;
- }
- }
+ trackno=imo;
+
+ }
- if(trackpid == kGamma || trackpid == 11 ||
- trackpid == -11 || trackpid == kPi0) igstatus = 1;
- mtrackpid = trackpid;
- mtrackno = trackno;
- trackpid = trackpidOld;
- trackno = tracknoOld;
+ if(trackpid==kGamma||trackpid==11||trackpid==-11||
+ trackpid==kPi0)igstatus=1;
+ mtrackpid=trackpid;
+ mtrackno=trackno;
+ trackpid=trackpidOld;
+ trackno=tracknoOld;
//-----------------end of modification----------------
xPos = fPMDHit->X();
}
AliDebug(2,Form("Zposition = %f Edeposition = %f",zPos,edep));
-
+ //Float_t zposition = TMath::Abs(zPos);
if (zPos < fZPos)
{
// CPV
//PS Test in September 2003 and 2006
// KeV - ADC conversion for 12bit ADC
// Modified by Ajay
- const Float_t kConstant = 9.0809;
- const Float_t kErConstant = 1.6763;
- const Float_t kSlope = 128.348;
- const Float_t kErSlope = 0.4703;
+ //const Float_t kConstant = 9.0809;
+ //const Float_t kErConstant = 1.6763;
+ //const Float_t kSlope = 128.348;
+ //const Float_t kErSlope = 0.4703;
+
+ const Float_t kConstant = -4.5;
+ const Float_t kErConstant = 14.3;
+ const Float_t kSlope = 68.2;
+ const Float_t kErSlope = 5.21;
Float_t cons = gRandom->Gaus(kConstant,kErConstant);
Float_t slop = gRandom->Gaus(kSlope,kErSlope);
#include "AliPMDddldata.h"
#include "AliPMDUtility.h"
#include "AliESDPmdTrack.h"
-#include "AliRawReader.h"
-
//#include "AliPMDRecoParam.h"
ClassImp(AliPMDQADataMakerRec)
{
// create Raws histograms in Raws subdir
-
+ const Bool_t saveForCorr = kTRUE;
// Preshower plane
TH1I * h0 = new TH1I("hPreEdepM0","ADC Distribution PRE - Module 0", 100, 0, 2000);
// Y vs. X for PRE and CPV planes
TH2F * h48 = new TH2F("hPreXY","PRE plane",200,-100.,100.,200,-100.,100.);
- Add2RawsList(h48, 48);
+ Add2RawsList(h48, 48,saveForCorr);//Ajay
+ //Add2RawsList(h48, 48);//Ajay
TH2F * h49 = new TH2F("hCpvXY","CPV plane",200,-100.,100.,200,-100.,100.);
- Add2RawsList(h49, 49);
+ Add2RawsList(h49, 49,saveForCorr);//Ajay
+ //Add2RawsList(h49, 49);//Ajay
}
//____________________________________________________________________________
{
//Fill prepared histograms with Raw digit properties
- rawReader->Reset() ;
- TObjArray *pmdddlcont = 0x0;
+ TObjArray *pmdddlcont = 0x0;
pmdddlcont = new TObjArray();
AliPMDRawStream stream(rawReader);