//don't use masked channels;
if (mROC ->GetValue(irow,ipad)) continue;
if (mRefROC->GetValue(irow,ipad)) continue;
- if (nRefROC->GetValue(irow,ipad)==0) continue;
+ if (nRefROC->GetValue(irow,ipad)==0) continue;
Float_t deviation=(nROC->GetValue(irow,ipad)/nRefROC->GetValue(irow,ipad))-1;
for (Int_t i=0;i<npar;++i){
if (deviation>vThres[i])
(*valArray[ipar])[naccept]=state[ipar];
naccept++;
}
+ if (naccept<2) return 0;
TMatrixD *pstat=new TMatrixD(9,3);
TMatrixD &stat=*pstat;
for (Int_t ipar=0; ipar<9; ipar++){
// sigmaCut - maximal allowed deviation from mean in terms of RMS
if (!array) return 0;
if (array->GetEntries()<=0) return 0;
+ if (!(&stat)) return 0;
// error increase in 1 hour
const Double_t kerrsTime[9]={
0.00001, 0.00001, 0.00001,
fVdriftBins = 500;
fVdriftStart= -0.1;
fVdriftEnd = 0.1;
- fRunBins = 100001;
+ fRunBins = 1000001;
fRunStart = -1.5;
- fRunEnd = 99999.5;
+ fRunEnd = 999999.5;
Int_t binsVdriftLaser[4] = {fTimeBins , fPtBins , fVdriftBins*20, fRunBins };
Double_t xminVdriftLaser[4] = {fTimeStart, fPtStart, fVdriftStart , fRunStart};
Double_t ptrelative1 = AliTPCcalibDB::GetPTRelative(tstamp,fRun,1);
Double_t temp0 = AliTPCcalibDB::GetTemperature(tstamp,fRun,0);
Double_t temp1 = AliTPCcalibDB::GetTemperature(tstamp,fRun,1);
+ Double_t vdcorr = AliTPCcalibDB::Instance()->GetVDriftCorrectionTime(tstamp,fRun,0,1);
TVectorD vecGoofie(20);
AliDCSSensorArray* goofieArray = AliTPCcalibDB::Instance()->GetGoofieSensors(fRun);
if (goofieArray){
"temp0="<<temp0<<
"temp1="<<temp1<<
"vecGoofie.="<<&vecGoofie<<
+ "vdcorr="<<vdcorr<<
//laser
"rejectA="<<isReject[0]<<
"rejectC="<<isReject[1]<<
Double_t ptrelative1 = AliTPCcalibDB::GetPTRelative(tstamp,fRun,1);
Double_t temp0 = AliTPCcalibDB::GetTemperature(tstamp,fRun,0);
Double_t temp1 = AliTPCcalibDB::GetTemperature(tstamp,fRun,1);
+ Double_t vdcorr = AliTPCcalibDB::Instance()->GetVDriftCorrectionTime(tstamp,fRun,0,1);
TVectorD vecGoofie(20);
AliDCSSensorArray* goofieArray = AliTPCcalibDB::Instance()->GetGoofieSensors(fRun);
if (goofieArray){
"temp0="<<temp0<<
"temp1="<<temp1<<
"vecGoofie.=<<"<<&vecGoofie<<
+ "vdcorr="<<vdcorr<<
//
// accumulated values
//
TTimeStamp tstamp(fTime);
Double_t ptrelative0 = AliTPCcalibDB::GetPTRelative(tstamp,fRun,0);
Double_t ptrelative1 = AliTPCcalibDB::GetPTRelative(tstamp,fRun,1);
+ Double_t vdcorr = AliTPCcalibDB::Instance()->GetVDriftCorrectionTime(tstamp,fRun,0,1);
(*cstream)<<"tpctpc"<<
"run="<<fRun<< // run number
"event="<<fEvent<< // event number
"mag="<<fMagF<< // magnetic field
"ptrel0.="<<ptrelative0<<
"ptrel1.="<<ptrelative1<<
+ "vdcorr="<<vdcorr<< // drift correction applied
//
"xyz.="<<&gxyz<< // global position
"tIn.="<<&trackIn<< // refitterd track in
pTPC.GetDirection(gdTPC.GetMatrixArray());
pITS.GetXYZ(gpITS.GetMatrixArray());
pITS.GetDirection(gdITS.GetMatrixArray());
+ Double_t vdcorr = AliTPCcalibDB::Instance()->GetVDriftCorrectionTime(tstamp,fRun,0,1);
(*cstream)<<"itstpc"<<
"run="<<fRun<< // run number
"event="<<fEvent<< // event number
"temp0="<<temp0<<
"temp1="<<temp1<<
"vecGoofie.="<<&vecGoofie<<
+ "vdcorr="<<vdcorr<< // drift correction applied
//
"nmed="<<kglast<< // number of entries to define median and RMS
"vMed.="<<&vecMedian<< // median of deltas
pTPC.GetDirection(gdTPC.GetMatrixArray());
pTRD.GetXYZ(gpTRD.GetMatrixArray());
pTRD.GetDirection(gdTRD.GetMatrixArray());
+ Double_t vdcorr = AliTPCcalibDB::Instance()->GetVDriftCorrectionTime(tstamp,fRun,0,1);
(*cstream)<<"trdtpc"<<
"run="<<fRun<< // run number
"event="<<fEvent<< // event number
"temp0="<<temp0<<
"temp1="<<temp1<<
"vecGoofie.="<<&vecGoofie<<
+ "vdcorr="<<vdcorr<< // drift correction applied
//
"nmed="<<kglast<< // number of entries to define median and RMS
"vMed.="<<&vecMedian<< // median of deltas
pTPC.GetDirection(gdTPC.GetMatrixArray());
pTOF.GetXYZ(gpTOF.GetMatrixArray());
pTOF.GetDirection(gdTOF.GetMatrixArray());
+ Double_t vdcorr = AliTPCcalibDB::Instance()->GetVDriftCorrectionTime(tstamp,fRun,0,1);
(*cstream)<<"toftpc"<<
"run="<<fRun<< // run number
"event="<<fEvent<< // event number
"temp0="<<temp0<<
"temp1="<<temp1<<
"vecGoofie.="<<&vecGoofie<<
+ "vdcorr="<<vdcorr<< // drift correction applied
//
"nmed="<<kglast<< // number of entries to define median and RMS
"vMed.="<<&vecMedian<< // median of deltas