* *
* Permission to use, copy, modify and distribute this software and its *
* documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
+ * without fee, provided that the above copyright notice appears in all *
* copies and that both the copyright notice and this permission notice *
* appear in the supporting documentation. The authors make no claims *
* about the suitability of this software for any purpose. It is *
////
/////////////////////////////////////////////////////////////////////////
-
+
#include "AliLog.h"
#include "AliTPCMonitor.h"
#include "AliTPCMonitorMappingHandler.h"
-#include "AliTPCMonitorDateFile.h"
-#include "AliTPCMonitorDateFormat.h"
-#include "AliTPCMonitorAltro.h"
#include "AliTPCMonitorFFT.h"
-#include "AliRawReaderRoot.h"
#include "AliRawReader.h"
+#include "AliRawReaderRoot.h"
#include "AliRawEventHeaderBase.h"
+#include "AliAltroRawStreamV3.h"
#include "TH2F.h"
#include "TF1.h"
#include "TMath.h"
#include "TPaveText.h"
#include "TFile.h"
#include <Riostream.h>
-#include "AliTPCMonitorDateMonitor.h"
ClassImp(AliTPCMonitor)
+const Int_t AliTPCMonitor::fgkHwMaskFEC = 0x0780;
+const Int_t AliTPCMonitor::fgkHwMaskBranch = 0x0800;
+const Int_t AliTPCMonitor::fgkHwMaskFECChannel = 0x007f;
+const Int_t AliTPCMonitor::fgkHwMaskAltroChannel = 0x000f;
+const Int_t AliTPCMonitor::fgkHwMaskAltroChip = 0x0070;
+const Int_t AliTPCMonitor::fgkHwMaskRCU = 0x7000;
+
//____________________________________________________________________________
AliTPCMonitor::AliTPCMonitor(char* name, char* title) :
- AliTPCMonitorConfig(name,title),
- fPad(new Int_t*[GetMaxHwAddr()]),
- fPadMapHw(new Float_t[GetMaxHwAddr()]),
- fPadMapRCU(new Int_t*[GetMaxHwAddr()]),
- fHistIROC(0),
- fHistOROC(0),
- fHistIROCIndex(0),
- fHistOROCIndex(0),
- fHistIROCTime(0),
- fHistOROCTime(0),
- fHistIROCClone(0),
- fHistOROCClone(0),
- fHistIROCRMS(0),
- fHistOROCRMS(0),
- fHistIROCBASE(0),
- fHistOROCBASE(0),
- fHistIROCSUM(0),
- fHistOROCSUM(0),
- fHistChannelTime(0),
- fHistAddrMapIndex(0),
- fHistAddrMaxAdc(0),
- fHistAddrBaseMean(0),
- fHistAddrMaxAdcX(0),
- fHistAddrAdcSum(0),
- fHistAddrBaseRms(0),
- fHistDistrSumIROC(0),
- fHistDistrMaxIROC(0),
- fHistDistrSumOROC(0),
- fHistDistrMaxOROC(0),
- fHistDistrBase2dIROC(0),
- fHistDistrBase2dOROC(0),
- fHistDistrBaseRmsIROC(0),
- fHistDistrBaseMeanIROC(0),
- fHistDistrBaseRmsOROC(0),
- fHistDistrBaseMeanOROC(0),
- fHistGlobalMaxA(0),
- fHistGlobalMaxC(0),
- fHistList(new TObjArray()),
-
- fkNRowsIroc(63),
- fkNRowsOroc(96),
-
- fkNPadsIroc(110),
- fkNPadsOroc(140),
- fkNPadMinIroc(-55),
- fkNPadMinOroc(-70),
- fkNPadMaxIroc(55),
- fkNPadMaxOroc(70),
- fVerb(0),
- fLastEv(0),
- fEventNumber(0),
- fEventNumberOld(0),
- fDisableFit(0),
- fExecGlob(0),
- fExecPlaneMax(0),
- fExecPadIrocRms(0),
- fExecPadOrocRms(0),
- fRunId(0),
- fEqId(0),
- fPadUsedRoc(-1),
- fPadUsedHwAddr(-1),
- fGdcId(0),
- fLdcId(0),
- fLdcIdOld(1),
- fMapEqidsSec(new Int_t*[36]),
- fMapEqidsRcu(new Int_t[1000]),
- fMirror(1),
- fChannelIter(0),
- fMapHand(0),
- fReaderROOT(0),
- fReaderDATE(0),
- fReaderDATEMon(0)
+AliTPCMonitorConfig(name,title),
+fPad(new Int_t*[GetMaxHwAddr()]),
+fPadMapHw(new Float_t[GetMaxHwAddr()]),
+fPadMapRCU(new Int_t*[GetMaxHwAddr()]),
+fHistIROC(0),
+fHistOROC(0),
+fHistIROCIndex(0),
+fHistOROCIndex(0),
+fHistIROCTime(0),
+fHistOROCTime(0),
+fHistIROCClone(0),
+fHistOROCClone(0),
+fHistIROCRMS(0),
+fHistOROCRMS(0),
+fHistIROCBASE(0),
+fHistOROCBASE(0),
+fHistIROCSUM(0),
+fHistOROCSUM(0),
+fHistChannelTime(0),
+fHistAddrMapIndex(0),
+fHistAddrMaxAdc(0),
+fHistAddrBaseMean(0),
+fHistAddrMaxAdcX(0),
+fHistAddrAdcSum(0),
+fHistAddrBaseRms(0),
+fHistDistrSumIROC(0),
+fHistDistrMaxIROC(0),
+fHistDistrSumOROC(0),
+fHistDistrMaxOROC(0),
+fHistDistrBase2dIROC(0),
+fHistDistrBase2dOROC(0),
+fHistDistrBaseRmsIROC(0),
+fHistDistrBaseMeanIROC(0),
+fHistDistrBaseRmsOROC(0),
+fHistDistrBaseMeanOROC(0),
+fHistGlobalMaxA(0),
+fHistGlobalMaxC(0),
+fHistList(new TObjArray()),
+
+fkNRowsIroc(63),
+fkNRowsOroc(96),
+
+fkNPadsIroc(110),
+fkNPadsOroc(140),
+fkNPadMinIroc(-55),
+fkNPadMinOroc(-70),
+fkNPadMaxIroc(55),
+fkNPadMaxOroc(70),
+fVerb(0),
+fLastEv(0),
+fEventNumber(0),
+fEventNumberOld(0),
+fDisableFit(0),
+fExecGlob(0),
+fExecPlaneMax(0),
+fExecPadIrocRms(0),
+fExecPadOrocRms(0),
+fRunId(0),
+fEqId(0),
+fPadUsedRoc(-1),
+fPadUsedHwAddr(-1),
+fGdcId(0),
+fLdcId(0),
+fLdcIdOld(1),
+fMapEqidsSec(new Int_t*[36]),
+fMapEqidsRcu(new Int_t[1000]),
+fMirror(1),
+fChannelIter(0),
+fMapHand(0),
+fRawReader(0)
{
- // Constructor
+ // Constructor
for(Int_t i = 0; i<GetMaxHwAddr(); i++) { fPad[i] = new Int_t[GetTimeBins()];}
for(Int_t i = 0; i<GetMaxHwAddr(); i++) { fPadMapRCU[i] = new Int_t[6];}
for(Int_t i = 0; i<36; i++) { fMapEqidsSec[i] = new Int_t[6];}
-
- if (gDirectory)
+
+ if (gDirectory)
+ {
+ if (!gDirectory->GetList())
{
- if (!gDirectory->GetList())
- {
- Warning("Build","Current directory is not a valid directory");
- return;
- }
- AliTPCMonitor *hold = (AliTPCMonitor*)gDirectory->GetList()->FindObject(GetName());
- if(hold)
- {
- Warning("Build","Replacing existing histogram: %s (Potential memory leak).",GetName());
- gDirectory->GetList()->Remove(hold);
- }
- gDirectory->Append(this);
+ Warning("Build","Current directory is not a valid directory");
+ return;
}
+ AliTPCMonitor *hold = (AliTPCMonitor*)gDirectory->GetList()->FindObject(GetName());
+ if(hold)
+ {
+ Warning("Build","Replacing existing histogram: %s (Potential memory leak).",GetName());
+ gDirectory->GetList()->Remove(hold);
+ }
+ gDirectory->Append(this);
+ }
CreateHistos();
SetEqIds();
-
+
}
//____________________________________________________________________________
AliTPCMonitor::AliTPCMonitor(const AliTPCMonitor &monitor):
- AliTPCMonitorConfig(monitor.GetName(),monitor.GetTitle()),
- fPad(new Int_t*[GetMaxHwAddr()]),
- fPadMapHw(new Float_t[GetMaxHwAddr()]),
- fPadMapRCU(new Int_t*[GetMaxHwAddr()]),
- fHistIROC(0),
- fHistOROC(0),
- fHistIROCIndex(0),
- fHistOROCIndex(0),
- fHistIROCTime(0),
- fHistOROCTime(0),
- fHistIROCClone(0),
- fHistOROCClone(0),
- fHistIROCRMS(0),
- fHistOROCRMS(0),
- fHistIROCBASE(0),
- fHistOROCBASE(0),
- fHistIROCSUM(0),
- fHistOROCSUM(0),
- fHistChannelTime(0),
- fHistAddrMapIndex(0),
- fHistAddrMaxAdc(0),
- fHistAddrBaseMean(0),
- fHistAddrMaxAdcX(0),
- fHistAddrAdcSum(0),
- fHistAddrBaseRms(0),
- fHistDistrSumIROC(0),
- fHistDistrMaxIROC(0),
- fHistDistrSumOROC(0),
- fHistDistrMaxOROC(0),
- fHistDistrBase2dIROC(0),
- fHistDistrBase2dOROC(0),
- fHistDistrBaseRmsIROC(0),
- fHistDistrBaseMeanIROC(0),
- fHistDistrBaseRmsOROC(0),
- fHistDistrBaseMeanOROC(0),
- fHistGlobalMaxA(0),
- fHistGlobalMaxC(0),
- fHistList(new TObjArray()),
- fkNRowsIroc(monitor.fkNRowsIroc),
- fkNRowsOroc(monitor.fkNRowsOroc),
- fkNPadsIroc(monitor.fkNPadsIroc),
- fkNPadsOroc(monitor.fkNPadsOroc),
- fkNPadMinIroc(monitor.fkNPadMinIroc),
- fkNPadMinOroc(monitor.fkNPadMinOroc),
- fkNPadMaxIroc(monitor.fkNPadMaxIroc),
- fkNPadMaxOroc(monitor.fkNPadMaxOroc),
- fVerb(monitor.fVerb),
- fLastEv(monitor.fLastEv),
- fEventNumber(monitor.fEventNumber),
- fEventNumberOld(monitor.fEventNumberOld),
- fDisableFit(monitor.fDisableFit),
- fExecGlob(monitor.fExecGlob),
- fExecPlaneMax(monitor.fExecPlaneMax),
- fExecPadIrocRms(monitor.fExecPadIrocRms),
- fExecPadOrocRms(monitor.fExecPadOrocRms),
- fRunId(monitor.fRunId),
- fEqId(monitor.fEqId),
- fPadUsedRoc(monitor.fPadUsedRoc),
- fPadUsedHwAddr(monitor.fPadUsedHwAddr),
- fGdcId(monitor.fGdcId),
- fLdcId(monitor.fLdcId),
- fLdcIdOld(monitor.fLdcIdOld),
- fMapEqidsSec(new Int_t*[36]),
- fMapEqidsRcu(new Int_t[1000]),
- fMirror(monitor.fMirror),
- fChannelIter(monitor.fChannelIter),
- fMapHand(monitor.fMapHand),
- fReaderROOT(monitor.fReaderROOT),
- fReaderDATE(monitor.fReaderDATE),
- fReaderDATEMon(monitor.fReaderDATEMon)
+AliTPCMonitorConfig(monitor.GetName(),monitor.GetTitle()),
+fPad(new Int_t*[GetMaxHwAddr()]),
+fPadMapHw(new Float_t[GetMaxHwAddr()]),
+fPadMapRCU(new Int_t*[GetMaxHwAddr()]),
+fHistIROC(0),
+fHistOROC(0),
+fHistIROCIndex(0),
+fHistOROCIndex(0),
+fHistIROCTime(0),
+fHistOROCTime(0),
+fHistIROCClone(0),
+fHistOROCClone(0),
+fHistIROCRMS(0),
+fHistOROCRMS(0),
+fHistIROCBASE(0),
+fHistOROCBASE(0),
+fHistIROCSUM(0),
+fHistOROCSUM(0),
+fHistChannelTime(0),
+fHistAddrMapIndex(0),
+fHistAddrMaxAdc(0),
+fHistAddrBaseMean(0),
+fHistAddrMaxAdcX(0),
+fHistAddrAdcSum(0),
+fHistAddrBaseRms(0),
+fHistDistrSumIROC(0),
+fHistDistrMaxIROC(0),
+fHistDistrSumOROC(0),
+fHistDistrMaxOROC(0),
+fHistDistrBase2dIROC(0),
+fHistDistrBase2dOROC(0),
+fHistDistrBaseRmsIROC(0),
+fHistDistrBaseMeanIROC(0),
+fHistDistrBaseRmsOROC(0),
+fHistDistrBaseMeanOROC(0),
+fHistGlobalMaxA(0),
+fHistGlobalMaxC(0),
+fHistList(new TObjArray()),
+fkNRowsIroc(monitor.fkNRowsIroc),
+fkNRowsOroc(monitor.fkNRowsOroc),
+fkNPadsIroc(monitor.fkNPadsIroc),
+fkNPadsOroc(monitor.fkNPadsOroc),
+fkNPadMinIroc(monitor.fkNPadMinIroc),
+fkNPadMinOroc(monitor.fkNPadMinOroc),
+fkNPadMaxIroc(monitor.fkNPadMaxIroc),
+fkNPadMaxOroc(monitor.fkNPadMaxOroc),
+fVerb(monitor.fVerb),
+fLastEv(monitor.fLastEv),
+fEventNumber(monitor.fEventNumber),
+fEventNumberOld(monitor.fEventNumberOld),
+fDisableFit(monitor.fDisableFit),
+fExecGlob(monitor.fExecGlob),
+fExecPlaneMax(monitor.fExecPlaneMax),
+fExecPadIrocRms(monitor.fExecPadIrocRms),
+fExecPadOrocRms(monitor.fExecPadOrocRms),
+fRunId(monitor.fRunId),
+fEqId(monitor.fEqId),
+fPadUsedRoc(monitor.fPadUsedRoc),
+fPadUsedHwAddr(monitor.fPadUsedHwAddr),
+fGdcId(monitor.fGdcId),
+fLdcId(monitor.fLdcId),
+fLdcIdOld(monitor.fLdcIdOld),
+fMapEqidsSec(new Int_t*[36]),
+fMapEqidsRcu(new Int_t[1000]),
+fMirror(monitor.fMirror),
+fChannelIter(monitor.fChannelIter),
+fMapHand(monitor.fMapHand),
+fRawReader(monitor.fRawReader)
{
// copy constructor
fHistAddrMaxAdcX=(TH1F*)monitor.fHistAddrMaxAdcX->Clone(); fHistList->Add(fHistAddrMaxAdcX);
fHistAddrAdcSum=(TH1F*)monitor.fHistAddrAdcSum->Clone(); fHistList->Add(fHistAddrAdcSum);
fHistAddrBaseRms=(TH1F*)monitor.fHistAddrBaseRms->Clone(); fHistList->Add(fHistAddrBaseRms);
-
+
fHistDistrSumIROC=(TH1F*)monitor.fHistDistrSumIROC->Clone(); fHistList->Add(fHistDistrSumIROC);
fHistDistrMaxIROC=(TH1F*)monitor.fHistDistrMaxIROC->Clone(); fHistList->Add(fHistDistrMaxIROC);
fHistDistrSumOROC=(TH1F*)monitor.fHistDistrSumOROC->Clone(); fHistList->Add(fHistDistrSumOROC);
fHistDistrMaxOROC=(TH1F*)monitor.fHistDistrMaxOROC->Clone(); fHistList->Add(fHistDistrMaxOROC);
-
+
fHistDistrBase2dIROC=(TH2F*)monitor.fHistDistrBase2dIROC->Clone(); fHistList->Add(fHistDistrBase2dIROC);
fHistDistrBase2dOROC=(TH2F*)monitor.fHistDistrBase2dOROC->Clone(); fHistList->Add(fHistDistrBase2dOROC);
fHistDistrBaseRmsIROC=(TH1D*)monitor.fHistDistrBaseRmsIROC->Clone(); fHistList->Add(fHistDistrBaseRmsIROC);
fHistGlobalMaxA=(TH2S*)monitor.fHistGlobalMaxA->Clone(); fHistList->Add(fHistGlobalMaxA);
fHistGlobalMaxC=(TH2S*)monitor.fHistGlobalMaxC->Clone(); fHistList->Add(fHistGlobalMaxC);
- // fPad = new Int_t*[monitor.GetMaxHwAddr()];
+ // fPad = new Int_t*[monitor.GetMaxHwAddr()];
for(Int_t i = 0; i<GetMaxHwAddr(); i++) { fPad[i] = new Int_t[monitor.GetTimeBins()];}
- //fPadMapRCU = new Int_t*[monitor.GetMaxHwAddr()];
+ //fPadMapRCU = new Int_t*[monitor.GetMaxHwAddr()];
for(Int_t i = 0; i<monitor.GetMaxHwAddr(); i++) { fPadMapRCU[i] = new Int_t[6];}
//fPadMapHw = new Float_t[monitor.GetMaxHwAddr()];
for(Int_t i = 0; i<36; i++) { fMapEqidsSec[i] = new Int_t[6];}
SetEqIds();
- if (gDirectory)
+ if (gDirectory)
+ {
+ if (!gDirectory->GetList())
{
- if (!gDirectory->GetList())
- {
- Warning("Build","Current directory is not a valid directory");
-
- }
- else
- {
- AliTPCMonitor *hold = (AliTPCMonitor*)gDirectory->GetList()->FindObject(monitor.GetName());
- if(hold)
- {
- Warning("Build","Replacing existing histogram: %s (Potential memory leak).",monitor.GetName());
- gDirectory->GetList()->Remove(hold);
- }
- gDirectory->Append(this);
- }
- }
+ Warning("Build","Current directory is not a valid directory");
+
+ }
+ else
+ {
+ AliTPCMonitor *hold = (AliTPCMonitor*)gDirectory->GetList()->FindObject(monitor.GetName());
+ if(hold)
+ {
+ Warning("Build","Replacing existing histogram: %s (Potential memory leak).",monitor.GetName());
+ gDirectory->GetList()->Remove(hold);
+ }
+ gDirectory->Append(this);
+ }
+ }
}
{
// assigment operator
if(this!=&monitor)
+ {
+ ((AliTPCMonitorConfig *)this)->operator=(monitor);
+ fkNRowsIroc=monitor.fkNRowsIroc;
+ fkNRowsOroc=monitor.fkNRowsOroc;
+ fkNPadsIroc=monitor.fkNPadsIroc;
+ fkNPadsOroc=monitor.fkNPadsOroc;
+ fkNPadMinIroc=monitor.fkNPadMinIroc;
+ fkNPadMinOroc=monitor.fkNPadMinOroc;
+ fkNPadMaxIroc=monitor.fkNPadMaxIroc;
+ fkNPadMaxOroc=monitor.fkNPadMaxOroc;
+ fVerb=monitor.fVerb;
+ fLastEv=monitor.fLastEv;
+ fEventNumber=monitor.fEventNumber;
+ fEventNumberOld=monitor.fEventNumberOld;
+ fDisableFit=monitor.fDisableFit;
+ fExecGlob=monitor.fExecGlob;
+ fExecPlaneMax=monitor.fExecPlaneMax;
+ fExecPadIrocRms=monitor.fExecPadIrocRms;
+ fExecPadOrocRms=monitor.fExecPadOrocRms;
+ fRunId=monitor.fRunId;
+ fEqId=monitor.fEqId;
+ fPadUsedRoc=monitor.fPadUsedRoc;
+ fPadUsedHwAddr=monitor.fPadUsedHwAddr;
+ fGdcId=monitor.fGdcId;
+ fLdcId=monitor.fLdcId;
+ fLdcIdOld=monitor.fLdcIdOld;
+ fMapHand=monitor.fMapHand;
+ fRawReader=monitor.fRawReader;
+
+
+ fHistList = new TObjArray();
+ fHistIROC=(TH2F*)monitor.fHistIROC->Clone(); fHistList->Add(fHistIROC);
+ fHistOROC=(TH2F*)monitor.fHistOROC->Clone(); fHistList->Add(fHistOROC);
+ fHistIROCIndex=(TH2S*)monitor.fHistIROCIndex->Clone(); fHistList->Add(fHistIROCIndex);
+ fHistOROCIndex=(TH2S*)monitor.fHistOROCIndex->Clone(); fHistList->Add(fHistOROCIndex);
+ fHistIROCTime=(TH2F*)monitor.fHistIROCTime->Clone(); fHistList->Add(fHistIROCTime);
+ fHistOROCTime=(TH2F*)monitor.fHistOROCTime->Clone(); fHistList->Add(fHistOROCTime);
+ fHistIROCClone=(TH2F*)monitor.fHistIROCClone->Clone(); fHistList->Add(fHistIROCClone);
+ fHistOROCClone=(TH2F*)monitor.fHistOROCClone->Clone(); fHistList->Add(fHistOROCClone);
+ fHistIROCRMS=(TH2F*)monitor.fHistIROCRMS->Clone(); fHistList->Add(fHistIROCRMS);
+ fHistOROCRMS=(TH2F*)monitor.fHistOROCRMS->Clone(); fHistList->Add(fHistOROCRMS);
+ fHistIROCBASE=(TH2F*)monitor.fHistIROCBASE->Clone(); fHistList->Add(fHistIROCBASE);
+ fHistOROCBASE=(TH2F*)monitor.fHistOROCBASE->Clone(); fHistList->Add(fHistOROCBASE);
+ fHistIROCSUM=(TH2F*)monitor.fHistIROCSUM->Clone(); fHistList->Add(fHistIROCSUM);
+ fHistOROCSUM=(TH2F*)monitor.fHistOROCSUM->Clone(); fHistList->Add(fHistOROCSUM);
+
+ fHistChannelTime=(TH2F*)monitor.fHistChannelTime->Clone(); fHistList->Add(fHistChannelTime);
+
+ fHistAddrMapIndex=(TH1F*)monitor.fHistAddrMapIndex->Clone(); fHistList->Add(fHistAddrMapIndex);
+ fHistAddrMaxAdc=(TH1F*)monitor.fHistAddrMaxAdc->Clone(); fHistList->Add(fHistAddrMaxAdc);
+ fHistAddrBaseMean=(TH1F*)monitor.fHistAddrBaseMean->Clone(); fHistList->Add(fHistAddrBaseMean);
+ fHistAddrMaxAdcX=(TH1F*)monitor.fHistAddrMaxAdcX->Clone(); fHistList->Add(fHistAddrMaxAdcX);
+ fHistAddrAdcSum=(TH1F*)monitor.fHistAddrAdcSum->Clone(); fHistList->Add(fHistAddrAdcSum);
+ fHistAddrBaseRms=(TH1F*)monitor.fHistAddrBaseRms->Clone(); fHistList->Add(fHistAddrBaseRms);
+
+
+ fHistDistrSumIROC=(TH1F*)monitor.fHistDistrSumIROC->Clone(); fHistList->Add(fHistDistrSumIROC);
+ fHistDistrMaxIROC=(TH1F*)monitor.fHistDistrMaxIROC->Clone(); fHistList->Add(fHistDistrMaxIROC);
+ fHistDistrSumOROC=(TH1F*)monitor.fHistDistrSumOROC->Clone(); fHistList->Add(fHistDistrSumOROC);
+ fHistDistrMaxOROC=(TH1F*)monitor.fHistDistrMaxOROC->Clone(); fHistList->Add(fHistDistrMaxOROC);
+
+ fHistDistrBase2dIROC=(TH2F*)monitor.fHistDistrBase2dIROC->Clone(); fHistList->Add(fHistDistrBase2dIROC);
+ fHistDistrBase2dOROC=(TH2F*)monitor.fHistDistrBase2dOROC->Clone(); fHistList->Add(fHistDistrBase2dOROC);
+ fHistDistrBaseRmsIROC=(TH1D*)monitor.fHistDistrBaseRmsIROC->Clone(); fHistList->Add(fHistDistrBaseRmsIROC);
+ fHistDistrBaseMeanIROC=(TH1D*)monitor.fHistDistrBaseMeanIROC->Clone(); fHistList->Add(fHistDistrBaseMeanIROC);
+ fHistDistrBaseRmsOROC=(TH1D*)monitor.fHistDistrBaseRmsOROC->Clone(); fHistList->Add(fHistDistrBaseRmsOROC);
+ fHistDistrBaseMeanOROC=(TH1D*)monitor.fHistDistrBaseMeanOROC->Clone(); fHistList->Add(fHistDistrBaseMeanOROC);
+
+ fHistGlobalMaxA=(TH2S*)monitor.fHistGlobalMaxA->Clone(); fHistList->Add(fHistGlobalMaxA);
+ fHistGlobalMaxC=(TH2S*)monitor.fHistGlobalMaxC->Clone(); fHistList->Add(fHistGlobalMaxC);
+
+ fPad = new Int_t*[monitor.GetMaxHwAddr()];
+ for(Int_t i = 0; i<GetMaxHwAddr(); i++) { fPad[i] = new Int_t[monitor.GetTimeBins()];}
+
+ fPadMapRCU = new Int_t*[monitor.GetMaxHwAddr()];
+ for(Int_t i = 0; i<monitor.GetMaxHwAddr(); i++) { fPadMapRCU[i] = new Int_t[6];}
+
+ fPadMapHw = new Float_t[monitor.GetMaxHwAddr()];
+
+ fMapEqidsRcu = new Int_t[1000];
+ fMapEqidsSec = new Int_t*[36];
+ for(Int_t i = 0; i<36; i++) { fMapEqidsSec[i] = new Int_t[6];}
+ SetEqIds();
+
+ if (gDirectory)
{
- ((AliTPCMonitorConfig *)this)->operator=(monitor);
- fkNRowsIroc=monitor.fkNRowsIroc;
- fkNRowsOroc=monitor.fkNRowsOroc;
- fkNPadsIroc=monitor.fkNPadsIroc;
- fkNPadsOroc=monitor.fkNPadsOroc;
- fkNPadMinIroc=monitor.fkNPadMinIroc;
- fkNPadMinOroc=monitor.fkNPadMinOroc;
- fkNPadMaxIroc=monitor.fkNPadMaxIroc;
- fkNPadMaxOroc=monitor.fkNPadMaxOroc;
- fVerb=monitor.fVerb;
- fLastEv=monitor.fLastEv;
- fEventNumber=monitor.fEventNumber;
- fEventNumberOld=monitor.fEventNumberOld;
- fDisableFit=monitor.fDisableFit;
- fExecGlob=monitor.fExecGlob;
- fExecPlaneMax=monitor.fExecPlaneMax;
- fExecPadIrocRms=monitor.fExecPadIrocRms;
- fExecPadOrocRms=monitor.fExecPadOrocRms;
- fRunId=monitor.fRunId;
- fEqId=monitor.fEqId;
- fPadUsedRoc=monitor.fPadUsedRoc;
- fPadUsedHwAddr=monitor.fPadUsedHwAddr;
- fGdcId=monitor.fGdcId;
- fLdcId=monitor.fLdcId;
- fLdcIdOld=monitor.fLdcIdOld;
- fMapHand=monitor.fMapHand;
- fReaderROOT=monitor.fReaderROOT;
- fReaderDATE=monitor.fReaderDATE;
-
-
- fHistList = new TObjArray();
- fHistIROC=(TH2F*)monitor.fHistIROC->Clone(); fHistList->Add(fHistIROC);
- fHistOROC=(TH2F*)monitor.fHistOROC->Clone(); fHistList->Add(fHistOROC);
- fHistIROCIndex=(TH2S*)monitor.fHistIROCIndex->Clone(); fHistList->Add(fHistIROCIndex);
- fHistOROCIndex=(TH2S*)monitor.fHistOROCIndex->Clone(); fHistList->Add(fHistOROCIndex);
- fHistIROCTime=(TH2F*)monitor.fHistIROCTime->Clone(); fHistList->Add(fHistIROCTime);
- fHistOROCTime=(TH2F*)monitor.fHistOROCTime->Clone(); fHistList->Add(fHistOROCTime);
- fHistIROCClone=(TH2F*)monitor.fHistIROCClone->Clone(); fHistList->Add(fHistIROCClone);
- fHistOROCClone=(TH2F*)monitor.fHistOROCClone->Clone(); fHistList->Add(fHistOROCClone);
- fHistIROCRMS=(TH2F*)monitor.fHistIROCRMS->Clone(); fHistList->Add(fHistIROCRMS);
- fHistOROCRMS=(TH2F*)monitor.fHistOROCRMS->Clone(); fHistList->Add(fHistOROCRMS);
- fHistIROCBASE=(TH2F*)monitor.fHistIROCBASE->Clone(); fHistList->Add(fHistIROCBASE);
- fHistOROCBASE=(TH2F*)monitor.fHistOROCBASE->Clone(); fHistList->Add(fHistOROCBASE);
- fHistIROCSUM=(TH2F*)monitor.fHistIROCSUM->Clone(); fHistList->Add(fHistIROCSUM);
- fHistOROCSUM=(TH2F*)monitor.fHistOROCSUM->Clone(); fHistList->Add(fHistOROCSUM);
-
- fHistChannelTime=(TH2F*)monitor.fHistChannelTime->Clone(); fHistList->Add(fHistChannelTime);
-
- fHistAddrMapIndex=(TH1F*)monitor.fHistAddrMapIndex->Clone(); fHistList->Add(fHistAddrMapIndex);
- fHistAddrMaxAdc=(TH1F*)monitor.fHistAddrMaxAdc->Clone(); fHistList->Add(fHistAddrMaxAdc);
- fHistAddrBaseMean=(TH1F*)monitor.fHistAddrBaseMean->Clone(); fHistList->Add(fHistAddrBaseMean);
- fHistAddrMaxAdcX=(TH1F*)monitor.fHistAddrMaxAdcX->Clone(); fHistList->Add(fHistAddrMaxAdcX);
- fHistAddrAdcSum=(TH1F*)monitor.fHistAddrAdcSum->Clone(); fHistList->Add(fHistAddrAdcSum);
- fHistAddrBaseRms=(TH1F*)monitor.fHistAddrBaseRms->Clone(); fHistList->Add(fHistAddrBaseRms);
-
-
- fHistDistrSumIROC=(TH1F*)monitor.fHistDistrSumIROC->Clone(); fHistList->Add(fHistDistrSumIROC);
- fHistDistrMaxIROC=(TH1F*)monitor.fHistDistrMaxIROC->Clone(); fHistList->Add(fHistDistrMaxIROC);
- fHistDistrSumOROC=(TH1F*)monitor.fHistDistrSumOROC->Clone(); fHistList->Add(fHistDistrSumOROC);
- fHistDistrMaxOROC=(TH1F*)monitor.fHistDistrMaxOROC->Clone(); fHistList->Add(fHistDistrMaxOROC);
-
- fHistDistrBase2dIROC=(TH2F*)monitor.fHistDistrBase2dIROC->Clone(); fHistList->Add(fHistDistrBase2dIROC);
- fHistDistrBase2dOROC=(TH2F*)monitor.fHistDistrBase2dOROC->Clone(); fHistList->Add(fHistDistrBase2dOROC);
- fHistDistrBaseRmsIROC=(TH1D*)monitor.fHistDistrBaseRmsIROC->Clone(); fHistList->Add(fHistDistrBaseRmsIROC);
- fHistDistrBaseMeanIROC=(TH1D*)monitor.fHistDistrBaseMeanIROC->Clone(); fHistList->Add(fHistDistrBaseMeanIROC);
- fHistDistrBaseRmsOROC=(TH1D*)monitor.fHistDistrBaseRmsOROC->Clone(); fHistList->Add(fHistDistrBaseRmsOROC);
- fHistDistrBaseMeanOROC=(TH1D*)monitor.fHistDistrBaseMeanOROC->Clone(); fHistList->Add(fHistDistrBaseMeanOROC);
-
- fHistGlobalMaxA=(TH2S*)monitor.fHistGlobalMaxA->Clone(); fHistList->Add(fHistGlobalMaxA);
- fHistGlobalMaxC=(TH2S*)monitor.fHistGlobalMaxC->Clone(); fHistList->Add(fHistGlobalMaxC);
-
- fPad = new Int_t*[monitor.GetMaxHwAddr()];
- for(Int_t i = 0; i<GetMaxHwAddr(); i++) { fPad[i] = new Int_t[monitor.GetTimeBins()];}
-
- fPadMapRCU = new Int_t*[monitor.GetMaxHwAddr()];
- for(Int_t i = 0; i<monitor.GetMaxHwAddr(); i++) { fPadMapRCU[i] = new Int_t[6];}
-
- fPadMapHw = new Float_t[monitor.GetMaxHwAddr()];
-
- fMapEqidsRcu = new Int_t[1000];
- fMapEqidsSec = new Int_t*[36];
- for(Int_t i = 0; i<36; i++) { fMapEqidsSec[i] = new Int_t[6];}
- SetEqIds();
-
- if (gDirectory)
- {
- if (!gDirectory->GetList())
- {
- Warning("Build","Current directory is not a valid directory");
- return *this;
- }
- AliTPCMonitor *hold = (AliTPCMonitor*)gDirectory->GetList()->FindObject(monitor.GetName());
- if(hold)
- {
- Warning("Build","Replacing existing histogram: %s (Potential memory leak).",monitor.GetName());
- gDirectory->GetList()->Remove(hold);
- }
- gDirectory->Append(this);
- }
- fReaderDATEMon=monitor.fReaderDATEMon;
+ if (!gDirectory->GetList())
+ {
+ Warning("Build","Current directory is not a valid directory");
+ return *this;
+ }
+ AliTPCMonitor *hold = (AliTPCMonitor*)gDirectory->GetList()->FindObject(monitor.GetName());
+ if(hold)
+ {
+ Warning("Build","Replacing existing histogram: %s (Potential memory leak).",monitor.GetName());
+ gDirectory->GetList()->Remove(hold);
+ }
+ gDirectory->Append(this);
}
+ }
return *this;
}
//____________________________________________________________________________
void AliTPCMonitor::CreateHistos()
{
- // Create histograms to be displayed
-
+ // Create histograms to be displayed
+
if(fVerb) cout << " create new ones " << endl;
fHistIROC = new TH2F("fHistIROC" ,"fHistIROC" ,fkNRowsIroc,0,fkNRowsIroc,fkNPadsIroc, fkNPadMinIroc, fkNPadMaxIroc); fHistList->Add(fHistIROC);
- fHistOROC = new TH2F("fHistOROC" ,"fHistOROC" ,fkNRowsOroc,0,fkNRowsOroc,fkNPadsOroc, fkNPadMinOroc, fkNPadMaxOroc); fHistList->Add(fHistOROC);
-
+ fHistOROC = new TH2F("fHistOROC" ,"fHistOROC" ,fkNRowsOroc,0,fkNRowsOroc,fkNPadsOroc, fkNPadMinOroc, fkNPadMaxOroc); fHistList->Add(fHistOROC);
+
fHistIROCIndex = new TH2S("fHistIROCIndex" ,"fHistIROCIndex" ,fkNRowsIroc,0,fkNRowsIroc,fkNPadsIroc, fkNPadMinIroc, fkNPadMaxIroc); fHistList->Add(fHistIROCIndex);
fHistOROCIndex = new TH2S("fHistOROCIndex" ,"fHistOROCIndex" ,fkNRowsOroc,0,fkNRowsOroc,fkNPadsOroc, fkNPadMinOroc, fkNPadMaxOroc); fHistList->Add(fHistOROCIndex);
-
+
fHistIROCTime = new TH2F("fHistIROCTime" ,"fHistIROCTime" ,fkNRowsIroc,0,fkNRowsIroc,fkNPadsIroc, fkNPadMinIroc, fkNPadMaxIroc); fHistList->Add(fHistIROCTime);
fHistOROCTime = new TH2F("fHistOROCTime" ,"fHistOROCTime" ,fkNRowsOroc,0,fkNRowsOroc,fkNPadsOroc, fkNPadMinOroc, fkNPadMaxOroc); fHistList->Add(fHistOROCTime);
fHistIROCRMS = new TH2F("fHistIROCRMS" ,"fHistIROCRMS" ,fkNRowsIroc,0,fkNRowsIroc,fkNPadsIroc, fkNPadMinIroc, fkNPadMaxIroc); fHistList->Add(fHistIROCRMS);
fHistOROCRMS = new TH2F("fHistOROCRMS" ,"fHistOROCRMS" ,fkNRowsOroc,0,fkNRowsOroc,fkNPadsOroc, fkNPadMinOroc, fkNPadMaxOroc); fHistList->Add(fHistOROCRMS);
-
+
fHistIROCSUM = new TH2F("fHistIROCSUM" ,"fHistIROCSUM" ,fkNRowsIroc,0,fkNRowsIroc,fkNPadsIroc, fkNPadMinIroc, fkNPadMaxIroc); fHistList->Add(fHistIROCSUM);
fHistOROCSUM = new TH2F("fHistOROCSUM" ,"fHistOROCSUM" ,fkNRowsOroc,0,fkNRowsOroc,fkNPadsOroc, fkNPadMinOroc, fkNPadMaxOroc); fHistList->Add(fHistOROCSUM);
-
+
fHistIROCBASE = new TH2F("fHistIROCBASE" ,"fHistIROCBASE" ,fkNRowsIroc,0,fkNRowsIroc,fkNPadsIroc, fkNPadMinIroc, fkNPadMaxIroc); fHistList->Add(fHistIROCBASE);
fHistOROCBASE = new TH2F("fHistOROCBASE" ,"fHistOROCBASE" ,fkNRowsOroc,0,fkNRowsOroc,fkNPadsOroc, fkNPadMinOroc, fkNPadMaxOroc); fHistList->Add(fHistOROCBASE);
-
-
+
+
fHistChannelTime = new TH2F("fHistChannelTime" ,"fHistChannelTime" ,GetNumOfChannels(),0,GetNumOfChannels(),GetTimeBins(),0,GetTimeBins());fHistList->Add(fHistChannelTime);
fHistAddrMapIndex = new TH1F("fHistAddrMapIndex" ,"fHistAddrMapIndex" ,GetMaxHwAddr() ,0,GetMaxHwAddr()); fHistList->Add(fHistAddrMapIndex);
fHistAddrMaxAdc = new TH1F("fHistAddrMaxAdc" ,"fHistAddrMaxAdc" ,GetMaxHwAddr(),0,GetMaxHwAddr()); fHistList->Add(fHistAddrMaxAdc);
fHistDistrMaxIROC = new TH1F("fHistDistrMaxIROC" ,"fHistDistrMaxIROC" ,500,0.0,1000.0); fHistList->Add(fHistDistrMaxIROC);
fHistDistrSumOROC = new TH1F("fHistDistrSumOROC" ,"fHistDistrSumOROC" ,400,0.0,4000.0); fHistList->Add(fHistDistrSumOROC);
fHistDistrMaxOROC = new TH1F("fHistDistrMaxOROC" ,"fHistDistrMaxOROC" ,500,0.0,1000.0); fHistList->Add(fHistDistrMaxOROC);
-
+
fHistDistrBase2dIROC = new TH2F("fHistDistrBase2dIROC" ,"fHistDistrBase2dIROC",100,0.0,100.0,100,0.0,10.0); fHistList->Add(fHistDistrBase2dIROC);
fHistDistrBase2dOROC = new TH2F("fHistDistrBase2dOROC" ,"fHistDistrBase2dOROC",100,0.0,100.0,100,0.0,10.0); fHistList->Add(fHistDistrBase2dOROC);
// Depending on the value of the sector id all sectors (sectorid == -1) are processed.
//
// In this case only the maximum values are calculated per pad and filled to the global histograms
- // In a second loop the last processed(displayed) sector will be processed (sectorid!=-1)
+ // In a second loop the last processed(displayed) sector will be processed (sectorid!=-1)
// again and the baseline rms and further quantities are calculated
//
// If only one sector should be processed SetProcOneSector(1) should be set.
// In this case only the specified (last/last displayed) sector will be processed.
//
// If GetProcNextEvent()==0 the same event will be processed again
-
-
+
+
Int_t sectorid = 0;
Int_t retflag = 0; // id of last sector + 1000, or error flag
if(GetProcNextEvent()==1 && fLastEv) { AliInfo("Last event already processed"); }
if(GetProcNextEvent()==0 || GetProcOneSector()==1 ) sectorid = GetLastSector();
else sectorid = -1;
- // first iteration
- retflag = ReadData(sectorid);
+ // first iteration
+ retflag = ReadDataNew(sectorid);
SetLastProcFile(GetFile());
DrawHists(3);
- // second iteration
- if(sectorid==-1 && retflag >1000)
- {
- AliInfo("Second read cycle");
- SetProcNextEvent(0);
- if(GetLastSectorDisplayed()==-1) {sectorid = GetLastSector() ; }
- else {sectorid = GetLastSectorDisplayed(); SetLastSector(sectorid) ; }
- retflag = ReadData(sectorid);
- }
+ // second iteration
+ if(sectorid==-1 && retflag >1000)
+ {
+ AliInfo("Second read cycle");
+ SetProcNextEvent(0);
+ if(GetLastSectorDisplayed()==-1) {sectorid = GetLastSector() ; }
+ else {sectorid = GetLastSectorDisplayed(); SetLastSector(sectorid) ; }
+ retflag = ReadDataNew(sectorid);
+ }
SetLastSectorDisplayed(sectorid) ;
fMapHand->ReadfecHwMap(GetLastSector());
SetEventProcessed(1);
return retflag;
}
-
//__________________________________________________________________
-Int_t AliTPCMonitor::ReadData(Int_t secid)
+Int_t AliTPCMonitor::ReadDataNew(Int_t secid)
{
- // Read Data File/Stream for specified Format.
- // Payload will be extracted from either ROOT or DATE format
- // and passed to FillHistsDecode for decoding of the adc information
-
- Int_t format = GetFormat();
-
-
- //if(format==2 && !gSystem->Getenv("DATE_ROOT")) { AliError("DATE not initialized on this system"); return 11;}
-
- if( format==0) {return ReadDataROOT(secid);}
- else if(format==1) {return ReadDataDATEFile(secid);}
- else if(format==2)
- {
- return ReadDataDATEStream(secid);
- }
-
- AliWarning("Function should already be left");
- return 11;
-}
-//__________________________________________________________________
-Int_t AliTPCMonitor::ReadDataDATEFile(Int_t secid)
-{
- // Read Data in Date format.
- // Date file and monitor will be handled in this function.
-
- if(fReaderROOT) { delete fReaderROOT ; fReaderROOT=0;}
-
-#ifdef ALI_DATE
- if(fReaderDATEMon) { delete fReaderDATEMon ; fReaderDATEMon=0;}
-#endif
-
- Char_t* eventPtr = 0;
- AliTPCMonitorDateFormat* dateform = 0;
-
- // Create objects //
- if( fReaderDATE==0 || ( fReaderDATE!=0 && (strcmp(GetLastProcFile(),GetFile())!=0) ) )
- {
- cout << " Create new DATE file "<< endl;
- if(fReaderDATE) { delete fReaderDATE ; fReaderDATE=0; }
- fReaderDATE = new AliTPCMonitorDateFile() ;
- fReaderDATE->SetName("fReaderDATE") ;
- fReaderDATE->OpenDateFile(GetFile()) ;
- fEventNumber =0;
- }
-
- // Get Event pointer ///
- if(fReaderDATE->IsDateFileOpen() ==false ) { AliError("Could not open Date File"); return 11 ; }
- // Rewind event if new event number is smaller than old one
- if(fEventNumberOld>fEventNumber ) { fReaderDATE->ResetFilePos(); }
- while(GetProcNextEvent() || fEventNumber==0)
- {
- fReaderDATE->ReadEvent();
- eventPtr = fReaderDATE->GetMemoryPointer();
- dateform = new AliTPCMonitorDateFormat(eventPtr);
- Int_t currentev = dateform->GetEventID();
- if(fEventNumber <= currentev ){ break; }
- }
-
- eventPtr = fReaderDATE->GetMemoryPointer();
-
- if(dateform==0) dateform = new AliTPCMonitorDateFormat(eventPtr);
- fEventNumber = dateform->GetEventID();
- fEventNumberOld = dateform->GetEventID();
-
- if(fReaderDATE->IsLastEvent()) fLastEv =1;
- if(dateform->IsEventEndOfRun()) { AliInfo("Event is end of run: eventType END_OF_RUN "); }
-
- if(fReaderDATE->IsEventValid() == false ) {AliInfo("Skipping Event because invalid"); return 10; }
- if(dateform->IsEventPhysicsEvent() == false ) {AliInfo("Skipping Header/event because not Physics event"); return 12; }
-
- ResetHistos() ;
-
- Int_t lastrcu = ReadDataDATESubEventLoop(dateform,secid);
-
- delete dateform;
- if(fVerb) cout << " last rcu " << lastrcu << endl;
- return lastrcu;
-}
-
-
-#ifdef ALI_DATE
-//__________________________________________________________________
-Int_t AliTPCMonitor::ReadDataDATEStream(Int_t secid)
-{
- // Read Data from DATE stream.
- // Can also be used for DATE file.
-
-
- if(fReaderROOT) { delete fReaderROOT ; fReaderROOT=0;}
- if(fReaderDATE) { delete fReaderDATE ; fReaderDATE =0; }
-
- Char_t* eventPtr = 0;
-
- AliTPCMonitorDateFormat* dateform =0 ;
-
- // Create objects ///
- if((fReaderDATEMon==0 || ( fReaderDATEMon!=0 && (strcmp(GetLastProcFile(),GetFile())!=0) )))
- {
- if(fReaderDATEMon!=0)
- {
- fReaderDATEMon->Logout();
- delete fReaderDATEMon;
- }
- fReaderDATEMon = new AliTPCMonitorDateMonitor();
- fReaderDATEMon->SetName("DMon");
- Int_t status = fReaderDATEMon->OpenMonitoring(GetFile());
- if(status) { AliError(Form("Could not read event online: Error: %s",fReaderDATEMon->DecodeError(status))); return 11; }
- }
-
- if(GetProcNextEvent())
- {
- fReaderDATEMon->Free();
- Int_t status = fReaderDATEMon->GetEvent();
- if(status) { AliError(Form("Could not read event online: Error: %s",fReaderDATEMon->DecodeError(status))); return 11 ;}
- //fReaderDATEMon->Logout();
- }
-
- eventPtr = fReaderDATEMon->GetEventPointerasChar();
-
- if(dateform==0) dateform = new AliTPCMonitorDateFormat(eventPtr);
- fEventNumber = dateform->GetEventID();
- fEventNumberOld = dateform->GetEventID();
-
- if(dateform->IsEventEndOfRun()) { AliInfo("Event is end of run: eventType END_OF_RUN "); }
- if(dateform->IsEventPhysicsEvent() == false ) {AliInfo("Skipping Header/event because not Physics event"); return 12; }
-
- ResetHistos() ;
-
- Int_t lastrcu = ReadDataDATESubEventLoop(dateform,secid);
-
- delete dateform;
- if(fVerb) cout << " last rcu " << lastrcu << endl;
- return lastrcu;
-}
-#else
-//__________________________________________________________________
-Int_t AliTPCMonitor::ReadDataDATEStream(Int_t /*secid*/)
-{
- // Read Data from DATE stream.
- // Can also be used for DATE file.
- // In case DATE is not install
- // this method is dummy
-
- AliError("DATE not initialized on this system");
- return 11;
-}
-#endif
-
-Int_t AliTPCMonitor::ReadDataDATESubEventLoop(AliTPCMonitorDateFormat* dateform, Int_t secid)
-{
- // Loop over DATE Subevents
-
- Bool_t exitSubEventLoop = false;
- Int_t lastrcuid = 0;
- Int_t lasteq = 0;
- Int_t start = 1;
- Int_t rcupatch = 0;
- Char_t* eventPtr = 0;
- UInt_t* eventPtrUI = 0;
-
- fChannelIter =0;
-
- while(exitSubEventLoop == false)
- {
- if(start==1) { dateform->GotoSubEventHeader() ; start=0;}
- else { dateform->GotoNextSubEventHeader(); }
- if(dateform->IsLastSubEventHeader()) exitSubEventLoop = true;
-
- if(fVerb) cout << " next subevent LDC " << (Int_t)dateform->GetSubEventLDC() << endl;
-
- lasteq = 0 ;
- Int_t neq = 0 ;
-
- while(lasteq==0)
- {
- if(neq ==0){ dateform->GotoFirstEquipment();}
- else { dateform->GotoNextEquipment();}
-
- fGdcId = dateform->GetSubEventGDC() ;
- fLdcId = dateform->GetSubEventLDC();
- fRunId = dateform->GetEventRunID();
- fEqId = dateform->GetEquipmentID();
- rcupatch = GetRCUPatch(fRunId, fEqId);
- lastrcuid = (rcupatch+1000);
- neq++;
-
- if(fLdcIdOld!=fLdcId && fChannelIter!=0) {
- if(secid==-1)
- {
- FillGlobal(GetLastSector());
- ResetArrays();
- fChannelIter =0;
- }
- else
- {
- return lastrcuid;
- }
- }
- fLdcIdOld = fLdcId ;
-
-
- if(dateform->IsLastEquipment() != false ) lasteq = 1;
-
- eventPtr = dateform->GetFirstDataPointer();
- eventPtrUI = (UInt_t *) eventPtr;
- Int_t payload = dateform->GetPayloadSize(); // 40Bit words
-
- if(fVerb)DumpHeader(dateform);
- if(fVerb) cout << "Check sector and fEqId " << endl;
-
- if(fVerb && secid >0 )
- {
- cout << "secid : "<< secid << " fEqId "<< fEqId << " equipment 0 form map " << endl;
- cout << fMapEqidsSec[secid][0] << endl;
- cout << " showed map_eqids " << endl;
- }
-
- if(CheckEqId(secid,fEqId))
- {
- if(fVerb) cout << " init altro " << endl;
- AliTPCMonitorAltro* altro = new AliTPCMonitorAltro((UInt_t *)eventPtrUI,(payload/4),1); //hier
- altro->SetWrite10Bit(GetWrite10Bit());
- altro->SetActFilename(GetFile());
- if(fVerb) cout << " allocated 10bit " << endl;
- altro->Allocate10BitArray();
- altro->Decodeto10Bit(fEqId);
- AliInfo(Form("Process eqid %i , patch %i ",fEqId,rcupatch%6));
- FillHistsDecode(altro,(rcupatch%6),secid);
- if(fChannelIter!=0 && secid==-1 ) SetSectorFilled(rcupatch/6);
- delete altro;
- }// if(CheckID)
- }// while last eq
- SetLastSector(rcupatch/6);
- }
- if(fChannelIter!=0 && secid==-1){ FillGlobal(GetLastSector()); }
- return lastrcuid;
-
-}
-//__________________________________________________________________
-Int_t AliTPCMonitor::ReadDataROOT(Int_t secid)
-{
- // Read in data in ROOT format
- if(fReaderDATE){ delete fReaderDATE ; fReaderDATE=0;}
- fChannelIter =0;
- Int_t lastrcuid = 0;
- Int_t rcupatch = 0;
- Int_t equipmentSize = 0;
- UChar_t* eventPtr = 0 ;
- UInt_t* eventPtrUI = 0 ;
- Int_t evtype = 0;
-
- if(fReaderROOT==0 || ( fReaderROOT!=0 && (strcmp(GetLastProcFile(),GetFile())!=0) ) )
- {
- if(fVerb) cout << "AliTPCMonitor::ReadDataROOT create new reader " << endl;
- if(fReaderROOT) { delete fReaderROOT ; fReaderROOT=0; }
- fReaderROOT = new AliRawReaderRoot(GetFile());
- if(!fReaderROOT) { AliError("Could not initiate AliRawReaderRoo "); return 10;}
- fEventNumber=0;
- }
-
- // reset to beginning of the event
- fReaderROOT->Reset();
+ // Read Data File/Stream for specified Format.
+ // Payload will be extracted from either ROOT or DATE format
+ // and passed to FillHistsDecode for decoding of the adc information
+
+ if (!fRawReader){
+ fRawReader = AliRawReader::Create(GetFile());
+ SetLastProcFile(GetFile());
+ } else if (strcmp(GetLastProcFile(),GetFile())!=0){
+ delete fRawReader;
+ fRawReader = AliRawReader::Create(GetFile());
+// printf("New file!!!\n");
+ SetLastProcFile(GetFile());
+ }
- // Rewind event if new event number is smaller than old one
- if(fEventNumberOld>fEventNumber) fReaderROOT->RewindEvents();
+ if (!fRawReader){
+ AliWarning("Coult not initialize raw reader");
+ return 11;
+ }
- while(GetProcNextEvent() || fEventNumber==0)
+ Bool_t skip=kTRUE;
+ while(skip && GetProcNextEvent())
{
- if(fVerb) cout << "AliTPCMonitor::ReadDataROOT get event " << endl;
- if(!fReaderROOT->NextEvent()) { AliError("Could not get next Event"); return 11 ;}
-// Int_t currentev = *(fReaderROOT->GetEventId());
- //!! think about
- Int_t currentev=fReaderROOT->GetEventIndex();
-// break;
- // skip all events but physics, calibration and software trigger events!
- UInt_t eventType=fReaderROOT->GetType();
+ if(fVerb) cout << "AliTPCMonitor::ReadDataNew get event " << endl;
+ if(fRawReader->IsA()==AliRawReaderRoot::Class()){
+ printf("Root, NextEvent: %d\n",GetEventID());
+ if (!fRawReader->GotoEvent(GetEventID())){AliError("Could not get next Event"); return 11 ;}
+ } else if(!fRawReader->NextEvent()) { AliError("Could not get next Event"); return 11 ;}
+ // skip all events but physics, calibration and software trigger events!
+ UInt_t eventType=fRawReader->GetType();
if ( !(eventType==AliRawEventHeaderBase::kPhysicsEvent ||
eventType==AliRawEventHeaderBase::kCalibrationEvent ||
eventType==AliRawEventHeaderBase::kSystemSoftwareTriggerEvent ||
eventType==AliRawEventHeaderBase::kDetectorSoftwareTriggerEvent) ) {
- if (fVerb) cout<< "Skipping event! Its neither of 'physics, calibration and software trigger event'" << endl;
- continue;
+ if (fVerb) cout<< "Skipping event! Its neither of 'physics, calibration and software trigger event'" << endl;
+ if(fRawReader->IsA()==AliRawReaderRoot::Class()){
+ if (fEventNumber<fRawReader->GetNumberOfEvents()-1) ++fEventNumber;
+ else {AliError("No more events");return 11;}
+ }
+ continue;
}
- if(fEventNumber <= currentev ){ break; }
+ skip=kFALSE;
+ //test if the TPC has data
+ UChar_t *data=0;
+ fRawReader->Select("TPC");
+ Int_t eventNr=fRawReader->GetEventIndex();
+ if (!fRawReader->ReadNextData(data)) {
+ skip=kTRUE;
+ printf("%d - %d\n",fEventNumber,fRawReader->GetNumberOfEvents());
+ if(fRawReader->IsA()==AliRawReaderRoot::Class()){
+ if (fEventNumber<fRawReader->GetNumberOfEvents()-1){
+ ++eventNr;
+ printf("inc conter\n");
+ }
+ else {
+ AliError("No more events");
+ return 11;
+ }
+ }
+ }
+ fEventNumber = eventNr;
+ fEventNumberOld = eventNr;
}
-// fEventNumber = *(fReaderROOT->GetEventId());
-// fEventNumberOld = *(fReaderROOT->GetEventId());
- fEventNumber = fReaderROOT->GetEventIndex();
- fEventNumberOld = fReaderROOT->GetEventIndex();
+// printf("secid: %d\n",secid);
- ResetHistos() ;
+ //========================== Histogram filling ======================
+ ResetHistos() ;
+ AliAltroRawStreamV3 *altro=new AliAltroRawStreamV3(fRawReader);
+ altro->SelectRawData("TPC");
+ altro->Reset();
+ fChannelIter=0;
- while(fReaderROOT->ReadHeader())
- {
- if(fVerb) cout << "AliTPCMonitor::ReadDataROOT read header " << endl;
- fGdcId = fReaderROOT->GetGDCId() ;
- fLdcId = fReaderROOT->GetLDCId() ;
- fRunId = fReaderROOT->GetRunNumber() ;
- equipmentSize = fReaderROOT->GetEquipmentSize();
- fEqId = fReaderROOT->GetEquipmentId();
- evtype = fReaderROOT->GetType();
- rcupatch = GetRCUPatch(fRunId, fEqId);
- lastrcuid = (rcupatch+1000);
-
- if(evtype==1) { AliWarning(Form("EventType==1 in event %i ",fEventNumber)); return 10; }
- if(equipmentSize==0) { AliWarning(Form("Equipmentsize ==0 in event %i ",fEventNumber)); return 10; }
-
- if(fVerb) DumpHeader(fReaderROOT);
-
- if(fLdcIdOld!=fLdcId && fChannelIter!=0) {
- if(secid==-1)
- {
- FillGlobal(GetLastSector());
- ResetArrays();
- fChannelIter =0;
- }
- else
- {
- return lastrcuid;
- }
- }
- fLdcIdOld = fLdcId ;
-
- if(CheckEqId(secid,fEqId))
- {
- fReaderROOT->ReadNextData(eventPtr);
- eventPtrUI = (UInt_t *) eventPtr;
- Int_t offset = 16;
- Int_t fsize = (Int_t)((equipmentSize/4) -offset) +1;
- AliTPCMonitorAltro* altro = new AliTPCMonitorAltro((UInt_t *)eventPtrUI,fsize,2);
- altro->SetWrite10Bit(GetWrite10Bit());
- if(fVerb) cout << "AliTPCMonitor::ReadDataROOT Set Write10bit to " << GetWrite10Bit() << endl;
- altro->SetActFilename(GetFile());
- altro->Allocate10BitArray();
- altro->Decodeto10Bit(fEqId);
- AliInfo(Form("Process sector %i eqid %i , patch %i ",rcupatch/6,fEqId,rcupatch%6));
- FillHistsDecode(altro,(rcupatch%6),secid);
- if(fChannelIter!=0 && secid==-1 ) SetSectorFilled(rcupatch/6);
- delete altro;
- }
- SetLastSector(rcupatch/6);
- }
- if(fChannelIter!=0 && secid==-1) { FillGlobal(GetLastSector()); }
- return lastrcuid;
-}
-
-
-//____________________________________________________________________________
-void AliTPCMonitor::FillHistsDecode(AliTPCMonitorAltro* altro ,Int_t rcupatch, Int_t secid)
-{
- // Decode Channels, calculate base mean and rms depending on the
- // value of secid (see ProcessEvent) and fill histograms
-
-
- if(fVerb) cout << "FillHistsDecode : rcupatch " << rcupatch << " id " << secid << endl;
- Int_t timestamp = 0;
- Int_t sampleiter = 0;
- Int_t samplelength = 0;
- Int_t samplebins = 0;
- Float_t max = 0;
- Float_t maxx = 0;
- Float_t sum = 0.0;
- Int_t sumn = 0;
- Int_t blockpos = 0;
Int_t hw = 0;
- Int_t nwords = 0;
Int_t nextHwAddress = 0;
- Int_t ntime = 0;
- Int_t adc = 0;
- Int_t nextpos = 0;
- Int_t lastpos = altro->Get10BitArraySize()-1;
- Short_t* entries = altro->Get10BitArray();
- Double_t hrms = 0.0;
- Double_t hmean = 0.0;
- Int_t supnextpos = 0;
-// TH1D* hbase = new TH1D("hbase","hbase",GetTimeBins(),0.5,(GetTimeBins()+0.5));
- Float_t fbase[1000];
-
-
- while(lastpos>0)
- {
- nextpos = altro->DecodeTrailer(lastpos);
- supnextpos = altro->GetNextTrailerPos();
- if(nextpos==-1) { break; }
- Int_t itimebin=0; //timebins in this pad
-
- lastpos = nextpos;
- blockpos = altro->GetTrailerBlockPos();
- hw = altro->GetTrailerHwAddress();
- nwords = altro->GetTrailerNWords();
- nextHwAddress = ( hw + (rcupatch<<12) );
- fPad[fChannelIter][0] = nextHwAddress ;
-
- if(fPadMapHw[nextHwAddress]!=-1 )
+ Int_t rcupatch = 0;
+ Int_t maxADC = 0;
+ Int_t maxx = 0;
+ Int_t sum = 0;
+ Int_t sumn = 0;
+ Int_t lastrcuid = 0;
+
+ while ( altro->NextDDL() ){
+ fGdcId = fRawReader->GetGDCId() ;
+ fLdcId = fRawReader->GetLDCId() ;
+ fRunId = fRawReader->GetRunNumber() ;
+ fEqId = fRawReader->GetEquipmentId();
+ rcupatch = GetRCUPatch(fRunId, fEqId);
+ Int_t rcupatchSector=rcupatch%6;
+ lastrcuid = (rcupatch+1000);
+// printf("RCU patch: %d, LDC: %d, EqId: %d\n",rcupatch, fLdcId, fEqId);
+ if(fLdcIdOld!=fLdcId && fChannelIter!=0) {
+ if(secid==-1)
{
- //Int_t hw_before1 = fPad[fChannelIter-2][0];
- //Int_t hw_before2 = fPad[fChannelIter-3][0];
-
- if(fVerb){ cout <<"\n //// Ambiguous hwaddress " << nextHwAddress << " write 10bit files and check file for eqid : " << fEqId << " /////// " << endl;}
- continue;
-
- if( TMath::Abs(fPadMapRCU[nextHwAddress][4] - fChannelIter)==1)
- {
- if(fVerb) cout << "//// Difference to previous channel==1 : reset branch bit of hw from last channel to 1 " << endl;
- Int_t hwreset = (nextHwAddress + (1<<11));
- fPad[fChannelIter-1][0] = hwreset;
-
- fPadMapHw[hwreset] = fChannelIter-1 ;
- fPadMapRCU[hwreset][0]= rcupatch ;
- fPadMapRCU[hwreset][1]= ((hwreset & AliTPCMonitorAltro::GetHwMaskBranch()) >> 11);
- fPadMapRCU[hwreset][2]= ((hwreset & AliTPCMonitorAltro::GetHwMaskFEC()) >>7 );
- fPadMapRCU[hwreset][3]= ( hwreset & AliTPCMonitorAltro::GetHwMaskFECChannel() );
- fPadMapRCU[hwreset][4]= fChannelIter-1;
- fPadMapRCU[hwreset][5]= altro->GetTrailerPos();
- }
- }
-
+ FillGlobal(GetLastSector());
+ ResetArrays();
+// printf("filled sector: %d\n",GetLastSector());
+ fChannelIter =0;
+ }
+ else
+ {
+ return lastrcuid;
+// if (rcupatch/6!=secid) continue;
+ }
+// fChannelIter=0;
+ }
+ if (!CheckEqId(secid,fEqId)) continue;
+ while ( altro->NextChannel() ){
+ hw=altro->GetHWAddress();
+ nextHwAddress = ( hw + (rcupatchSector<<12) );
+ fPad[fChannelIter][0] = nextHwAddress ;
fPadMapHw[nextHwAddress] = fChannelIter ;
- fPadMapRCU[nextHwAddress][0]= rcupatch ;
- fPadMapRCU[nextHwAddress][1]= ((nextHwAddress & AliTPCMonitorAltro::GetHwMaskBranch())>> 11) ;
- fPadMapRCU[nextHwAddress][2]= ((nextHwAddress & AliTPCMonitorAltro::GetHwMaskFEC()) >>7);
- fPadMapRCU[nextHwAddress][3]= (nextHwAddress & AliTPCMonitorAltro::GetHwMaskFECChannel());
- fPadMapRCU[nextHwAddress][4]= fChannelIter;
- fPadMapRCU[nextHwAddress][5]= altro->GetTrailerPos();
-
- timestamp = 0;
- sampleiter = 0;
- samplelength = 0;
- samplebins = 0;
-
- max = 0.0;
- maxx = 0.0;
- sum = 0.0;
- sumn = 0;
-
-// hbase->Reset();
-
- for(Int_t iterwords = 0 ; iterwords< nwords ; iterwords++)
- {
- if(entries[blockpos-iterwords]==682) { continue; }
- if(entries[blockpos-iterwords]!=682 && sampleiter==0)
- {
- samplelength = entries[blockpos-iterwords];
- sampleiter = samplelength;
- samplebins = 0;
- timestamp = entries[blockpos-iterwords-1];
- iterwords++;
- sampleiter-=2;
- }
- else
- {
- ntime = timestamp-samplebins;
- adc = entries[blockpos-iterwords];
- fPad[fChannelIter][ntime] = adc;
-//if( (adc!=0) && (ntime>=GetRangeBaseMin() ) && (ntime<GetRangeBaseMax() )) {hbase->Fill(adc) ;}
- if( (adc!=0) && (ntime>=GetRangeBaseMin() ) && (ntime<GetRangeBaseMax() )) {fbase[itimebin]=adc;itimebin++;}
- if( (adc>max) && (ntime>=GetRangeMaxAdcMin()) && (ntime<GetRangeMaxAdcMax() )) {max = adc;maxx = ntime ;}
- if( (ntime>=GetRangeSumMin()) && (ntime<GetRangeSumMax() )) {sum+=adc; sumn++;}
- samplebins++;
- sampleiter--;
- }
- }
-
-// hmean = hbase->GetMean();
-// hbase->GetXaxis()->SetRangeUser(hmean- hmean/3 , hmean + hmean/3);
-// hmean = hbase->GetMean();
-// hrms = hbase->GetRMS();
- if (itimebin>0){
- hmean = TMath::Mean(itimebin, fbase);
- hrms = TMath::RMS(itimebin, fbase);
- }
-
- if( GetPedestals()==1) fHistAddrMaxAdc->SetBinContent( nextHwAddress,max- hmean);
- else fHistAddrMaxAdc->SetBinContent( nextHwAddress,max );
+ maxADC=0;
+ while ( altro->NextBunch() ){
+ Int_t startTbin = (Int_t)altro->GetStartTimeBin();
+ Int_t bunchlength = (Int_t)altro->GetBunchLength();
+ const UShort_t *sig = altro->GetSignals();
+ for (Int_t iTimeBin = 0; iTimeBin<bunchlength; iTimeBin++){
+ Int_t adc=(Int_t)sig[iTimeBin];
+ Int_t ntime=startTbin-iTimeBin;
+ //fill channel information
+ fPad[fChannelIter][ntime] = adc;
+ if( (adc>maxADC) && (ntime>=GetRangeMaxAdcMin()) && (ntime<GetRangeMaxAdcMax() )) {maxADC = adc;maxx = ntime ;}
+ if( (ntime>=GetRangeSumMin()) && (ntime<GetRangeSumMax() )) {sum+=adc; sumn++;}
+
+ }
+ }
+ //get pedestal, noise
+ Float_t pedestal=TMath::Mean(GetRangeBaseMax()-GetRangeBaseMin(),fPad[fChannelIter]+GetRangeBaseMin());
+ Float_t noise =TMath::RMS(GetRangeBaseMax()-GetRangeBaseMin(),fPad[fChannelIter]+GetRangeBaseMin());
+ fHistAddrMaxAdc->SetBinContent(nextHwAddress,maxADC-GetPedestals()*pedestal);
if(secid!=-1)
- {
- if(rcupatch<2)
- {
- fHistDistrBase2dIROC->Fill(hmean,hrms);
- fHistDistrSumIROC->Fill(sum);
- if( GetPedestals()==1 ) { fHistDistrMaxIROC->Fill(max-hmean); fHistDistrSumIROC->Fill(sum -sumn*hmean);}
- else { fHistDistrMaxIROC->Fill(max); fHistDistrSumIROC->Fill(sum );}
- }
- else
- {
- fHistDistrBase2dOROC->Fill(hmean,hrms);
- fHistDistrSumOROC->Fill(sum);
- if( GetPedestals()==1 ){ fHistDistrMaxOROC->Fill(max-hmean); fHistDistrSumOROC->Fill(sum -sumn*hmean);}
- else { fHistDistrMaxOROC->Fill(max); fHistDistrSumOROC->Fill(sum) ;}
- }
-
- fHistAddrAdcSum->SetBinContent( nextHwAddress,sum);
- fHistAddrMapIndex->SetBinContent(nextHwAddress,fChannelIter);
- fHistAddrBaseMean->SetBinContent(nextHwAddress,hmean);
- fHistAddrMaxAdcX->SetBinContent( nextHwAddress,maxx);
- fHistAddrBaseRms->SetBinContent( nextHwAddress,hrms);
- }
- fChannelIter++;
- if(nextpos<0) { AliError("Error : next pos < 0 "); break ;}
+ {
+ if(rcupatchSector<2)
+ {
+ fHistDistrBase2dIROC->Fill(pedestal,noise);
+ fHistDistrSumIROC->Fill(sum);
+ fHistDistrMaxIROC->Fill(maxADC-pedestal*GetPedestals());
+ fHistDistrSumIROC->Fill(sum -sumn*pedestal*GetPedestals());
+ }
+ else
+ {
+ fHistDistrBase2dOROC->Fill(pedestal,noise);
+ fHistDistrSumOROC->Fill(sum);
+ fHistDistrMaxOROC->Fill(maxADC-pedestal*GetPedestals());
+ fHistDistrSumOROC->Fill(sum -sumn*pedestal*GetPedestals());
+ }
+
+ fHistAddrAdcSum->SetBinContent( nextHwAddress,sum);
+ fHistAddrMapIndex->SetBinContent(nextHwAddress,fChannelIter);
+ fHistAddrBaseMean->SetBinContent(nextHwAddress,pedestal);
+ fHistAddrMaxAdcX->SetBinContent( nextHwAddress,maxx);
+ fHistAddrBaseRms->SetBinContent( nextHwAddress,noise);
+ }
+
+ ++fChannelIter;
}
-// delete hbase;
- return ;
+ SetLastSector(rcupatch/6);
+ if(fChannelIter!=0 && secid==-1 ) SetSectorFilled(rcupatch/6);
+ fLdcIdOld = fLdcId ;
+ }
+ delete altro;
+ if(fChannelIter!=0 && secid==-1) { FillGlobal(GetLastSector());}
+ return lastrcuid;
}
+
+
+
+
+
+
+
+
+
//____________________________________________________________________________
void AliTPCMonitor::FillHistsPadPlane()
{
if(fVerb)cout << "AliTPCMonitor::FillHistsPadPlane() Start " << endl;
if(fVerb)PrintConfig();
-
+
Int_t pad = 0;
Int_t row = 0;
Int_t padmax = 0;
Int_t hwadd = 0;
- for(Int_t ch = 0; ch<fChannelIter; ch++)
+ for(Int_t ch = 0; ch<fChannelIter; ch++)
+ {
+ hwadd= fPad[ch][0];
+ fHistChannelTime->SetCellContent(ch,0,hwadd);
+
+ for(Int_t bin = 1; bin <GetTimeBins(); bin++)
{
- hwadd= fPad[ch][0];
- fHistChannelTime->SetCellContent(ch,0,hwadd);
-
- for(Int_t bin = 1; bin <GetTimeBins(); bin++)
- {
- if( fHistChannelTime->GetCellContent(ch,bin)!=0) cout << " cellcontent already set " << endl;
- if( GetPedestals()==1 ) fHistChannelTime->SetCellContent(ch,bin,(fPad[ch][bin]- fHistAddrBaseMean->GetBinContent(hwadd)));
- else fHistChannelTime->SetCellContent(ch,bin,(fPad[ch][bin]));
- }
-
- pad = fMapHand->GetPad( hwadd);
- row = fMapHand->GetPadRow(hwadd);
- padmax = fMapHand->GetNumofPads(row);
-
- if(row<63)
- {
- fHistIROC->SetCellContent( row +1 ,pad +55 -padmax/2 +1,fHistAddrMaxAdc->GetBinContent( hwadd));
- fHistIROCIndex->SetCellContent(row +1 ,pad +55 -padmax/2 +1,ch);
- fHistIROCRMS->SetCellContent( row +1 ,pad +55 -padmax/2 +1,fHistAddrBaseRms->GetBinContent( hwadd));
- fHistIROCBASE->SetCellContent( row +1 ,pad +55 -padmax/2 +1,fHistAddrBaseMean->GetBinContent(hwadd));
- fHistIROCSUM->SetCellContent( row +1 ,pad +55 -padmax/2 +1,fHistAddrAdcSum->GetBinContent( hwadd));
- }
- else
- {
- fHistOROC->SetCellContent( row-63 +1 ,pad +70 -padmax/2 +1,fHistAddrMaxAdc->GetBinContent( hwadd));
- fHistOROCIndex->SetCellContent(row-63 +1 ,pad +70 -padmax/2 +1,ch);
- fHistOROCRMS->SetCellContent( row-63 +1 ,pad +70 -padmax/2 +1,fHistAddrBaseRms->GetBinContent( hwadd));
- fHistOROCBASE->SetCellContent( row-63 +1 ,pad +70 -padmax/2 +1,fHistAddrBaseMean->GetBinContent(hwadd));
- fHistOROCSUM->SetCellContent( row-63 +1 ,pad +70 -padmax/2 +1,fHistAddrAdcSum->GetBinContent( hwadd));
- }
+ if( fHistChannelTime->GetCellContent(ch,bin)!=0) cout << " cellcontent already set " << endl;
+ if( GetPedestals()==1 ) fHistChannelTime->SetCellContent(ch,bin,(fPad[ch][bin]- fHistAddrBaseMean->GetBinContent(hwadd)));
+ else fHistChannelTime->SetCellContent(ch,bin,(fPad[ch][bin]));
}
+
+ pad = fMapHand->GetPad( hwadd);
+ row = fMapHand->GetPadRow(hwadd);
+ padmax = fMapHand->GetNumofPads(row);
+
+ if(row<63)
+ {
+ fHistIROC->SetCellContent( row +1 ,pad +55 -padmax/2 +1,fHistAddrMaxAdc->GetBinContent( hwadd));
+ fHistIROCIndex->SetCellContent(row +1 ,pad +55 -padmax/2 +1,ch);
+ fHistIROCRMS->SetCellContent( row +1 ,pad +55 -padmax/2 +1,fHistAddrBaseRms->GetBinContent( hwadd));
+ fHistIROCBASE->SetCellContent( row +1 ,pad +55 -padmax/2 +1,fHistAddrBaseMean->GetBinContent(hwadd));
+ fHistIROCSUM->SetCellContent( row +1 ,pad +55 -padmax/2 +1,fHistAddrAdcSum->GetBinContent( hwadd));
+ }
+ else
+ {
+ fHistOROC->SetCellContent( row-63 +1 ,pad +70 -padmax/2 +1,fHistAddrMaxAdc->GetBinContent( hwadd));
+ fHistOROCIndex->SetCellContent(row-63 +1 ,pad +70 -padmax/2 +1,ch);
+ fHistOROCRMS->SetCellContent( row-63 +1 ,pad +70 -padmax/2 +1,fHistAddrBaseRms->GetBinContent( hwadd));
+ fHistOROCBASE->SetCellContent( row-63 +1 ,pad +70 -padmax/2 +1,fHistAddrBaseMean->GetBinContent(hwadd));
+ fHistOROCSUM->SetCellContent( row-63 +1 ,pad +70 -padmax/2 +1,fHistAddrAdcSum->GetBinContent( hwadd));
+ }
+ }
fHistChannelTime->GetXaxis()->SetRange(0,fChannelIter);
fHistChannelTime->GetYaxis()->SetRange(0,GetTimeBins());
//____________________________________________________________________________
void AliTPCMonitor::ResetArrays()
{
- // Reset data arrays
- for(Int_t row = 0 ; row < fkNRowsIroc; row++)
+ // Reset data arrays
+ for(Int_t row = 0 ; row < fkNRowsIroc; row++)
+ {
+ for(Int_t pad = 0 ; pad < fkNPadsIroc ; pad++)
{
- for(Int_t pad = 0 ; pad < fkNPadsIroc ; pad++)
- {
- fHistIROCIndex->SetCellContent(row+1,pad+1,-1);
- }
+ fHistIROCIndex->SetCellContent(row+1,pad+1,-1);
}
- for(Int_t row = 0 ; row < fkNRowsOroc; row++)
+ }
+ for(Int_t row = 0 ; row < fkNRowsOroc; row++)
+ {
+ for(Int_t pad = 0 ; pad < fkNPadsOroc ; pad++)
{
- for(Int_t pad = 0 ; pad < fkNPadsOroc ; pad++)
- {
- fHistOROCIndex->SetCellContent(row+1,pad+1,-1);
- }
+ fHistOROCIndex->SetCellContent(row+1,pad+1,-1);
}
+ }
- for(Int_t ch= 0; ch<GetMaxHwAddr(); ch++)
- {
- fHistAddrMaxAdcX->SetBinContent(ch,-1);
- fHistAddrMapIndex->SetBinContent(ch,-1);
- fHistAddrMaxAdc->SetBinContent( ch, 0);
- fHistAddrBaseMean->SetBinContent( ch, 0);
- fHistAddrAdcSum->SetBinContent( ch, 0);
- fHistAddrBaseRms->SetBinContent(ch, 0);
- }
+ for(Int_t ch= 0; ch<GetMaxHwAddr(); ch++)
+ {
+ fHistAddrMaxAdcX->SetBinContent(ch,-1);
+ fHistAddrMapIndex->SetBinContent(ch,-1);
+ fHistAddrMaxAdc->SetBinContent( ch, 0);
+ fHistAddrBaseMean->SetBinContent( ch, 0);
+ fHistAddrAdcSum->SetBinContent( ch, 0);
+ fHistAddrBaseRms->SetBinContent(ch, 0);
+ }
for(Int_t ch = 0; ch< GetNumOfChannels(); ch++)
+ {
+ for(Int_t bin = 0; bin< GetTimeBins(); bin++)
{
- for(Int_t bin = 0; bin< GetTimeBins(); bin++)
- {
- fPad[ch][bin] = 0;
- }
- }
- for(Int_t ch = 0; ch< GetMaxHwAddr(); ch++)
- {
- fPadMapHw[ch]=-1;
- fPadMapRCU[ch][0]=-1;
- fPadMapRCU[ch][1]=-1;
- fPadMapRCU[ch][2]=-1;
- fPadMapRCU[ch][3]=-1;
- fPadMapRCU[ch][4]=-1;
- fPadMapRCU[ch][5]=-1;
+ fPad[ch][bin] = 0;
}
+ }
+ for(Int_t ch = 0; ch< GetMaxHwAddr(); ch++)
+ {
+ fPadMapHw[ch]=-1;
+ fPadMapRCU[ch][0]=-1;
+ fPadMapRCU[ch][1]=-1;
+ fPadMapRCU[ch][2]=-1;
+ fPadMapRCU[ch][3]=-1;
+ fPadMapRCU[ch][4]=-1;
+ fPadMapRCU[ch][5]=-1;
+ }
fPadUsedHwAddr=-1;
}
//____________________________________________________________________________
void AliTPCMonitor::ResetHistos()
{
- // Reset all but
+ // Reset all but
for(Int_t i =0; i<fHistList->GetEntries(); i++)
- {
- if(GetProcNextEvent()==0 && strcmp(((TH1*)fHistList->At(i))->GetName(),"SIDE A")==0) continue;
- if(GetProcNextEvent()==0 && strcmp(((TH1*)fHistList->At(i))->GetName(),"SIDE C")==0) continue;
- ((TH1*)fHistList->At(i))->Reset();
- }
+ {
+ if(GetProcNextEvent()==0 && strcmp(((TH1*)fHistList->At(i))->GetName(),"SIDE A")==0) continue;
+ if(GetProcNextEvent()==0 && strcmp(((TH1*)fHistList->At(i))->GetName(),"SIDE C")==0) continue;
+ ((TH1*)fHistList->At(i))->Reset();
+ }
ResetArrays();
}
//____________________________________________________________________________
void AliTPCMonitor::DeleteHistos()
{
- // Delete histograms
+ // Delete histograms
for(Int_t i =0; i<fHistList->GetEntries(); i++)
- {
- delete (TH1*)fHistList->At(i);
- }
+ {
+ delete (TH1*)fHistList->At(i);
+ }
}
//__________________________________________________________________
Int_t AliTPCMonitor::CheckEqId(Int_t secid,Int_t eqid)
{
- // Check if equipment id corresponds to any rcu patch in sector
+ // Check if equipment id corresponds to any rcu patch in sector
// Equipment ids changed during commisioning in 2006 (starting from run 704)
// However Runids started from 0 again in 2007
// Now only runids from commissioning in 2006 after runid 704 and all new ones are supported.
// Comment in equipment check for runids < 704 if old runs should be processed
-
+
if(fVerb) cout << "AliTPCMonitor::CheckEqId : SectorId " << secid << " EquipmentId " << eqid << " runid " << fRunId << endl;
Int_t passed =1;
//skip all eqids which do not belong to the TPC
if ( eqid<768||eqid>983 ) return 0;
//
if(fRunId<704 && 0) // commented out --> runs with runid < 704 in 2006 are not recognized anymore
+ {
+ if( (secid>-1) && (secid<36) ) // if ( secid is in range) { take only specific eqids} else { take all }
{
- if( (secid>-1) && (secid<36) ) // if ( secid is in range) { take only specific eqids} else { take all }
- {
- if( (secid==13) && ( eqid!=408 && eqid!=409 && eqid!=509 && eqid!=512 && eqid!=513 && eqid!=517 )) {passed=0;}
- else if( (secid==4) && ( eqid!=404 && eqid!=504 && eqid!=407 && eqid!=503 && eqid!=508 && eqid!=506 )) {passed=0;}
- }
- else {if(fVerb) cout << "passed check "<< endl; }
+ if( (secid==13) && ( eqid!=408 && eqid!=409 && eqid!=509 && eqid!=512 && eqid!=513 && eqid!=517 )) {passed=0;}
+ else if( (secid==4) && ( eqid!=404 && eqid!=504 && eqid!=407 && eqid!=503 && eqid!=508 && eqid!=506 )) {passed=0;}
}
+ else {if(fVerb) cout << "passed check "<< endl; }
+ }
else
+ {
+ if( (secid>-1) && (secid<36) ) // if ( secid is in range) { take only specific eqids} else { take all }
{
- if( (secid>-1) && (secid<36) ) // if ( secid is in range) { take only specific eqids} else { take all }
- {
- if(eqid!=fMapEqidsSec[secid][0] && eqid!= fMapEqidsSec[secid][1] && eqid!=fMapEqidsSec[secid][2] &&
- eqid!=fMapEqidsSec[secid][3] && eqid!= fMapEqidsSec[secid][4] && eqid!=fMapEqidsSec[secid][5] ) {passed=0;}
- }
- else {if(fVerb) cout << "passed check "<< endl;}
+ if(eqid!=fMapEqidsSec[secid][0] && eqid!= fMapEqidsSec[secid][1] && eqid!=fMapEqidsSec[secid][2] &&
+ eqid!=fMapEqidsSec[secid][3] && eqid!= fMapEqidsSec[secid][4] && eqid!=fMapEqidsSec[secid][5] ) {passed=0;}
}
-
+ else {if(fVerb) cout << "passed check "<< endl;}
+ }
+
return passed;
}
void AliTPCMonitor::SetEqIds()
{
// Set mapping for equipment ids
- for(Int_t i = 0; i<36 ; i++)
- {
- for(Int_t j = 0; j<6; j++)
- {
- if(j<2) fMapEqidsSec[i][j]= 768+i*2+j;
- else fMapEqidsSec[i][j]= 840+i*4+j-2;
- }
+ for(Int_t i = 0; i<36 ; i++)
+ {
+ for(Int_t j = 0; j<6; j++)
+ {
+ if(j<2) fMapEqidsSec[i][j]= 768+i*2+j;
+ else fMapEqidsSec[i][j]= 840+i*4+j-2;
}
+ }
for(Int_t i = 0; i<36 ; i++)
- {
- for(Int_t j = 0; j<6; j++)
- {
- if(j<2) fMapEqidsRcu[768+i*2+j] = i*6 +j;
- else fMapEqidsRcu[840+i*4+j-2] = i*6 +j;
- }
+ {
+ for(Int_t j = 0; j<6; j++)
+ {
+ if(j<2) fMapEqidsRcu[768+i*2+j] = i*6 +j;
+ else fMapEqidsRcu[840+i*4+j-2] = i*6 +j;
}
+ }
}
//__________________________________________________________________
void AliTPCMonitor::FillGlobal(Int_t sector)
{
-
- // Fill global histograms with max adc for each channel
+
+ // Fill global histograms with max adc for each channel
TH2S* hglob =0;
if((sector/18) ==0) hglob = fHistGlobalMaxA;
Float_t xdr = 0;
Float_t ydr = 0;
- for(Int_t hw = 0; hw<fHistAddrMaxAdc->GetNbinsX(); hw++)
+ for(Int_t hw = 0; hw<fHistAddrMaxAdc->GetNbinsX(); hw++)
+ {
+ max = (Int_t)fHistAddrMaxAdc->GetBinContent(hw);
+ if(max!=-1)
{
- max = (Int_t)fHistAddrMaxAdc->GetBinContent(hw);
- if(max!=-1)
- {
- pad = fMapHand->GetPad( hw);
- row = fMapHand->GetPadRow(hw);
- padmax = fMapHand->GetNumofPads(row);
- if (sector%36>17) fMirror=-1;
- else fMirror=1;
- GetXY(xval ,yval , padmax,row ,pad);
- xdr = xval*m11 +yval*m12;
- ydr = xval*m21 +yval*m22;
- if(hglob->GetBinContent(hglob->GetXaxis()->FindBin(xdr),hglob->GetYaxis()->FindBin(ydr))==0) hglob->Fill(xdr,ydr,(Int_t)max);
- }
+ pad = fMapHand->GetPad( hw);
+ row = fMapHand->GetPadRow(hw);
+ padmax = fMapHand->GetNumofPads(row);
+ if (sector%36>17) fMirror=-1;
+ else fMirror=1;
+ GetXY(xval ,yval , padmax,row ,pad);
+ xdr = xval*m11 +yval*m12;
+ ydr = xval*m21 +yval*m22;
+ if(hglob->GetBinContent(hglob->GetXaxis()->FindBin(xdr),hglob->GetYaxis()->FindBin(ydr))==0) hglob->Fill(xdr,ydr,(Int_t)max);
}
+ }
}
{
// Get x and y position of pad
- if(row<63)
- {
- xval = fMirror*( 2*padmax -4*pad -2);
- yval = 852.25 +7.5*row;
- }
- else
- {
- xval = fMirror*( 3*padmax -6*pad -3);
- if((row-63)<63) { yval = 10*(row-63) +1351; }
- else { yval = 15*(row-63-64)+1993.5; }
- }
+ if(row<63)
+ {
+ xval = fMirror*( 2*padmax -4*pad -2);
+ yval = 852.25 +7.5*row;
+ }
+ else
+ {
+ xval = fMirror*( 3*padmax -6*pad -3);
+ if((row-63)<63) { yval = 10*(row-63) +1351; }
+ else { yval = 15*(row-63-64)+1993.5; }
+ }
}
//__________________________________________________________________
Int_t AliTPCMonitor::GetPadAtX(Float_t xval, Int_t row) const
{
-
+
// Get pad number at given position in x
-
+
Int_t padmax = fMapHand->GetNumofPads(row);
Int_t pad = 0;
//__________________________________________________________________
void AliTPCMonitor::DrawHists(Int_t histos)
{
-
+
// Draw sets of histograms
// histos==1 : 2Dim histos for MAX adc and add executables
- // histos==2 : distributions max/rms/sum
+ // histos==2 : distributions max/rms/sum
// histos==3 : global max adc for specified SideA/C
-
-
- if(fVerb) cout << " Draw histos " << endl;
+
+
+ if(fVerb) cout << " Draw histos " << endl;
Char_t cside[10];
if(GetLastSector()/18==0 ) sprintf(cside,"A");
else sprintf(cside,"C");
Char_t titleBASE[256]; sprintf(titleBASE ,"Baseline RMS<->Mean (timebin: %i-%i) %s" ,GetRangeBaseMin() ,GetRangeBaseMax() ,titleSEC);
Char_t titleMEAN[256]; sprintf(titleMEAN ,"Baseline Mean (timebin: %i-%i) %s" ,GetRangeBaseMin() ,GetRangeBaseMax() ,titleSEC);
Char_t titleRMS[256] ; sprintf(titleRMS ,"Baseline RMS (timebin: %i-%i) %s" ,GetRangeBaseMin() ,GetRangeBaseMax() ,titleSEC);
-
- if(histos==1)
- {
+
+ if(histos==1)
+ {
// IROC _______________________________________________________________
- TCanvas* ciroc = 0;
- ciroc = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("ciroc");
- if(!ciroc)
- {
- ciroc = CreateCanvas("ciroc");
- fExecPlaneMax=0;
- }
- ciroc->cd();
-
- fHistIROC->SetXTitle("row");
- fHistIROC->SetYTitle("pad");
- if(GetPedestals()) fHistIROC->SetZTitle("max ADC (baseline sub)");
- else fHistIROC->SetZTitle("max ADC ");
- fHistIROC->SetTitle(titleIROC);
- fHistIROC->SetMinimum(0.01);
- fHistIROC->Draw("COLZ");
- ciroc->UseCurrentStyle();
-
-
- fHistIROCTime->SetXTitle("row"); fHistIROCTime->SetZTitle("peak time (fit)"); fHistIROCTime->SetYTitle("pad"); fHistIROCTime->SetTitle(titleIROC);
- fHistIROCRMS->SetXTitle("row"); fHistIROCRMS->SetZTitle( "baseline rms (ADC)"); fHistIROCRMS->SetYTitle("pad"); fHistIROCRMS->SetTitle(titleIROC);
-
+ TCanvas* ciroc = 0;
+ ciroc = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("ciroc");
+ if(!ciroc)
+ {
+ ciroc = CreateCanvas("ciroc");
+ fExecPlaneMax=0;
+ }
+ ciroc->cd();
+
+ fHistIROC->SetXTitle("row");
+ fHistIROC->SetYTitle("pad");
+ if(GetPedestals()) fHistIROC->SetZTitle("max ADC (baseline sub)");
+ else fHistIROC->SetZTitle("max ADC ");
+ fHistIROC->SetTitle(titleIROC);
+ fHistIROC->SetMinimum(0.01);
+ fHistIROC->Draw("COLZ");
+ ciroc->UseCurrentStyle();
+
+
+ fHistIROCTime->SetXTitle("row"); fHistIROCTime->SetZTitle("peak time (fit)"); fHistIROCTime->SetYTitle("pad"); fHistIROCTime->SetTitle(titleIROC);
+ fHistIROCRMS->SetXTitle("row"); fHistIROCRMS->SetZTitle( "baseline rms (ADC)"); fHistIROCRMS->SetYTitle("pad"); fHistIROCRMS->SetTitle(titleIROC);
+
// OROC
- TCanvas* coroc = 0;
- coroc =(TCanvas*)gROOT->GetListOfCanvases()->FindObject("coroc");
- if(!coroc) {
- coroc = CreateCanvas("coroc");
- fExecPlaneMax=0;
- }
- coroc->cd();
-
- fHistOROC->SetXTitle("row");
- fHistOROC->SetYTitle("pad");
- if(GetPedestals()) fHistOROC->SetZTitle("max ADC (baseline sub)");
- else fHistOROC->SetZTitle("max ADC ");
- fHistOROC->SetTitle(titleOROC);
- fHistOROC->SetMinimum(0.01);
- fHistOROC->Draw("COLZ");
- coroc->UseCurrentStyle();
-
-
- fHistOROCTime->SetXTitle("row"); fHistOROCTime->SetZTitle("peak time (fit) (timebins)"); fHistOROCTime->SetYTitle("pad"); fHistOROCTime->SetTitle(titleOROC);
- fHistOROCRMS->SetXTitle("row"); fHistOROCRMS->SetZTitle("baseline rms (ADC)"); fHistOROCRMS->SetYTitle("pad"); fHistOROCRMS->SetTitle(titleOROC);
-
- // SUM
- Char_t namesum[256] ; sprintf(namesum,"ADC sum (bins: %i, %i)",GetRangeSumMin() ,GetRangeSumMax() );
- fHistIROCSUM->SetXTitle("row"); fHistIROCSUM->SetZTitle(namesum); fHistIROCSUM->SetYTitle("pad"); fHistIROCSUM->SetTitle(titleIROC);
- fHistOROCSUM->SetXTitle("row"); fHistOROCSUM->SetZTitle(namesum); fHistOROCSUM->SetYTitle("pad"); fHistOROCSUM->SetTitle(titleOROC);
+ TCanvas* coroc = 0;
+ coroc =(TCanvas*)gROOT->GetListOfCanvases()->FindObject("coroc");
+ if(!coroc) {
+ coroc = CreateCanvas("coroc");
+ fExecPlaneMax=0;
+ }
+ coroc->cd();
+
+ fHistOROC->SetXTitle("row");
+ fHistOROC->SetYTitle("pad");
+ if(GetPedestals()) fHistOROC->SetZTitle("max ADC (baseline sub)");
+ else fHistOROC->SetZTitle("max ADC ");
+ fHistOROC->SetTitle(titleOROC);
+ fHistOROC->SetMinimum(0.01);
+ fHistOROC->Draw("COLZ");
+ coroc->UseCurrentStyle();
+
+
+ fHistOROCTime->SetXTitle("row"); fHistOROCTime->SetZTitle("peak time (fit) (timebins)"); fHistOROCTime->SetYTitle("pad"); fHistOROCTime->SetTitle(titleOROC);
+ fHistOROCRMS->SetXTitle("row"); fHistOROCRMS->SetZTitle("baseline rms (ADC)"); fHistOROCRMS->SetYTitle("pad"); fHistOROCRMS->SetTitle(titleOROC);
+
+ // SUM
+ Char_t namesum[256] ; sprintf(namesum,"ADC sum (bins: %i, %i)",GetRangeSumMin() ,GetRangeSumMax() );
+ fHistIROCSUM->SetXTitle("row"); fHistIROCSUM->SetZTitle(namesum); fHistIROCSUM->SetYTitle("pad"); fHistIROCSUM->SetTitle(titleIROC);
+ fHistOROCSUM->SetXTitle("row"); fHistOROCSUM->SetZTitle(namesum); fHistOROCSUM->SetYTitle("pad"); fHistOROCSUM->SetTitle(titleOROC);
// BASE
- Char_t namebase[256] ; sprintf(namebase ,"base mean (timbebin: %i, %i )",GetRangeBaseMin(),GetRangeBaseMax());
- fHistIROCBASE->SetXTitle("row"); fHistIROCBASE->SetZTitle(namebase); fHistIROCBASE->SetYTitle("pad"); fHistIROCBASE->SetTitle(titleIROC);
- fHistOROCBASE->SetXTitle("row"); fHistOROCBASE->SetZTitle(namebase); fHistOROCBASE->SetYTitle("pad"); fHistOROCBASE->SetTitle(titleOROC);
-
- if(fHistIROCClone) fHistIROCClone->Delete();
- if(fHistOROCClone) fHistOROCClone->Delete();
- fHistIROCClone = (TH2F*)fHistIROC->Clone("fHistIROCClone");
- fHistOROCClone = (TH2F*)fHistOROC->Clone("fHistOROCClone");
-
+ Char_t namebase[256] ; sprintf(namebase ,"base mean (timbebin: %i, %i )",GetRangeBaseMin(),GetRangeBaseMax());
+ fHistIROCBASE->SetXTitle("row"); fHistIROCBASE->SetZTitle(namebase); fHistIROCBASE->SetYTitle("pad"); fHistIROCBASE->SetTitle(titleIROC);
+ fHistOROCBASE->SetXTitle("row"); fHistOROCBASE->SetZTitle(namebase); fHistOROCBASE->SetYTitle("pad"); fHistOROCBASE->SetTitle(titleOROC);
+
+ if(fHistIROCClone) fHistIROCClone->Delete();
+ if(fHistOROCClone) fHistOROCClone->Delete();
+ fHistIROCClone = (TH2F*)fHistIROC->Clone("fHistIROCClone");
+ fHistOROCClone = (TH2F*)fHistOROC->Clone("fHistOROCClone");
+
// Executables
- if(fExecPlaneMax==0)
- {
- Char_t carry1[100];
- sprintf(carry1,".x %s/TPC/AliTPCMonitorExec.C(1)",gSystem->Getenv("ALICE_ROOT"));
- ciroc->AddExec("pad",carry1);
- coroc->AddExec("pad",carry1);
-
- Char_t carry2[100];
- sprintf(carry2,".x %s/TPC/AliTPCMonitorExec.C(2)",gSystem->Getenv("ALICE_ROOT"));
- ciroc->AddExec("row",carry2);
- coroc->AddExec("row",carry2);
- fExecPlaneMax=1;
- }
- coroc->Update();
- ciroc->Update();
- }
- else if(histos==2)
+ if(fExecPlaneMax==0)
{
+ Char_t carry1[100];
+ sprintf(carry1,".x %s/TPC/AliTPCMonitorExec.C(1)",gSystem->Getenv("ALICE_ROOT"));
+ ciroc->AddExec("pad",carry1);
+ coroc->AddExec("pad",carry1);
+
+ Char_t carry2[100];
+ sprintf(carry2,".x %s/TPC/AliTPCMonitorExec.C(2)",gSystem->Getenv("ALICE_ROOT"));
+ ciroc->AddExec("row",carry2);
+ coroc->AddExec("row",carry2);
+ fExecPlaneMax=1;
+ }
+ coroc->Update();
+ ciroc->Update();
+ }
+ else if(histos==2)
+ {
// MAX ADC distribution ____________________________________________
- TCanvas* cmax = 0;
- cmax = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("cmax");
- if(!cmax) cmax = CreateCanvas("cmax");
+ TCanvas* cmax = 0;
+ cmax = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("cmax");
+ if(!cmax) cmax = CreateCanvas("cmax");
- cmax->cd();
- fHistDistrMaxIROC->GetXaxis()->SetRangeUser(0.0,1000.0);
- fHistDistrMaxIROC->SetXTitle("max ADC (ADC)");
- fHistDistrMaxIROC->SetYTitle("counts");
- fHistDistrMaxIROC->SetTitle(titleMAX);
- fHistDistrMaxIROC->Draw("");
- fHistDistrMaxOROC->SetLineColor(2);
- fHistDistrMaxOROC->Draw("same");
+ cmax->cd();
+ fHistDistrMaxIROC->GetXaxis()->SetRangeUser(0.0,1000.0);
+ fHistDistrMaxIROC->SetXTitle("max ADC (ADC)");
+ fHistDistrMaxIROC->SetYTitle("counts");
+ fHistDistrMaxIROC->SetTitle(titleMAX);
+ fHistDistrMaxIROC->Draw("");
+ fHistDistrMaxOROC->SetLineColor(2);
+ fHistDistrMaxOROC->Draw("same");
- if(fHistDistrMaxOROC->GetMaximum()> fHistDistrMaxIROC->GetMaximum()) fHistDistrMaxIROC->SetMaximum(fHistDistrMaxOROC->GetMaximum()*1.1);
-
- TLegend* legio = new TLegend(0.6,0.6,0.8,0.8);
- legio->SetFillColor(0);
- legio->AddEntry(fHistDistrMaxIROC,"IROC","l");
- legio->AddEntry(fHistDistrMaxOROC,"OROC","l");
- legio->Draw("same");
+ if(fHistDistrMaxOROC->GetMaximum()> fHistDistrMaxIROC->GetMaximum()) fHistDistrMaxIROC->SetMaximum(fHistDistrMaxOROC->GetMaximum()*1.1);
+
+ TLegend* legio = new TLegend(0.6,0.6,0.8,0.8);
+ legio->SetFillColor(0);
+ legio->AddEntry(fHistDistrMaxIROC,"IROC","l");
+ legio->AddEntry(fHistDistrMaxOROC,"OROC","l");
+ legio->Draw("same");
// ADC sum distribution
- TCanvas* csum = 0;
- csum = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("csum");
- if(!csum) csum = CreateCanvas("csum") ;
- csum->cd();
+ TCanvas* csum = 0;
+ csum = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("csum");
+ if(!csum) csum = CreateCanvas("csum") ;
+ csum->cd();
- fHistDistrSumIROC->SetXTitle("sum ADC (ADC)");
- fHistDistrSumIROC->SetYTitle("counts");
- fHistDistrSumIROC->SetTitle(titleSUM);
- fHistDistrSumIROC->Draw("");
- fHistDistrSumOROC->SetLineColor(2);
- fHistDistrSumOROC->Draw("same");
- if(fHistDistrSumOROC->GetMaximum()> fHistDistrSumIROC->GetMaximum()) fHistDistrSumIROC->SetMaximum(fHistDistrSumOROC->GetMaximum()*1.1);
- legio->Draw("same");
+ fHistDistrSumIROC->SetXTitle("sum ADC (ADC)");
+ fHistDistrSumIROC->SetYTitle("counts");
+ fHistDistrSumIROC->SetTitle(titleSUM);
+ fHistDistrSumIROC->Draw("");
+ fHistDistrSumOROC->SetLineColor(2);
+ fHistDistrSumOROC->Draw("same");
+ if(fHistDistrSumOROC->GetMaximum()> fHistDistrSumIROC->GetMaximum()) fHistDistrSumIROC->SetMaximum(fHistDistrSumOROC->GetMaximum()*1.1);
+ legio->Draw("same");
// BASELINE MEAN distribution
- TCanvas* cbasemean = 0;
- cbasemean = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("cbasemean");
- if(!cbasemean) cbasemean = CreateCanvas("cbasemean");
- cbasemean->cd();
+ TCanvas* cbasemean = 0;
+ cbasemean = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("cbasemean");
+ if(!cbasemean) cbasemean = CreateCanvas("cbasemean");
+ cbasemean->cd();
- fHistDistrBaseMeanIROC = fHistDistrBase2dIROC->ProjectionX("fHistDistrBaseMeanIROC");
- fHistDistrBaseMeanIROC->SetXTitle("base mean (ADC)");
- fHistDistrBaseMeanIROC->SetYTitle("counts");
- fHistDistrBaseMeanIROC->SetTitle(titleMEAN);
- fHistDistrBaseMeanIROC->Draw("");
+ fHistDistrBaseMeanIROC = fHistDistrBase2dIROC->ProjectionX("fHistDistrBaseMeanIROC");
+ fHistDistrBaseMeanIROC->SetXTitle("base mean (ADC)");
+ fHistDistrBaseMeanIROC->SetYTitle("counts");
+ fHistDistrBaseMeanIROC->SetTitle(titleMEAN);
+ fHistDistrBaseMeanIROC->Draw("");
- fHistDistrBaseMeanOROC = fHistDistrBase2dOROC->ProjectionX("fHistDistrBaseMeanOROC");
- fHistDistrBaseMeanOROC->SetLineColor(2);
- fHistDistrBaseMeanOROC->Draw("same");
- if(fHistDistrBaseMeanOROC->GetMaximum()>fHistDistrBaseMeanIROC->GetMaximum()) fHistDistrBaseMeanIROC->SetMaximum(fHistDistrBaseMeanOROC->GetMaximum()*1.1);
- legio->Draw("same");
-
- TCanvas* cbaserms = 0;
- cbaserms = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("cbaserms");
- if(!cbaserms) cbaserms = CreateCanvas("cbaserms") ;
- cbaserms->cd();
+ fHistDistrBaseMeanOROC = fHistDistrBase2dOROC->ProjectionX("fHistDistrBaseMeanOROC");
+ fHistDistrBaseMeanOROC->SetLineColor(2);
+ fHistDistrBaseMeanOROC->Draw("same");
+ if(fHistDistrBaseMeanOROC->GetMaximum()>fHistDistrBaseMeanIROC->GetMaximum()) fHistDistrBaseMeanIROC->SetMaximum(fHistDistrBaseMeanOROC->GetMaximum()*1.1);
+ legio->Draw("same");
+
+ TCanvas* cbaserms = 0;
+ cbaserms = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("cbaserms");
+ if(!cbaserms) cbaserms = CreateCanvas("cbaserms") ;
+ cbaserms->cd();
// BASELINE RMS distribution
- fHistDistrBaseRmsIROC = fHistDistrBase2dIROC->ProjectionY("fHistDistrBaseRmsIROC");
- fHistDistrBaseRmsIROC->SetXTitle("base rms (ADC)");
- fHistDistrBaseRmsIROC->SetYTitle("counts");
- fHistDistrBaseRmsIROC->SetTitle(titleRMS);
- fHistDistrBaseRmsIROC->Draw("");
+ fHistDistrBaseRmsIROC = fHistDistrBase2dIROC->ProjectionY("fHistDistrBaseRmsIROC");
+ fHistDistrBaseRmsIROC->SetXTitle("base rms (ADC)");
+ fHistDistrBaseRmsIROC->SetYTitle("counts");
+ fHistDistrBaseRmsIROC->SetTitle(titleRMS);
+ fHistDistrBaseRmsIROC->Draw("");
- fHistDistrBaseRmsOROC = fHistDistrBase2dOROC->ProjectionY("fHistDistrBaseRmsOROC");
- fHistDistrBaseRmsOROC->SetLineColor(2);
- fHistDistrBaseRmsOROC->Draw("same");
- if(fHistDistrBaseRmsOROC->GetMaximum()>fHistDistrBaseRmsIROC->GetMaximum()) fHistDistrBaseRmsIROC->SetMaximum(fHistDistrBaseRmsOROC->GetMaximum()*1.1);
- legio->Draw("same");
+ fHistDistrBaseRmsOROC = fHistDistrBase2dOROC->ProjectionY("fHistDistrBaseRmsOROC");
+ fHistDistrBaseRmsOROC->SetLineColor(2);
+ fHistDistrBaseRmsOROC->Draw("same");
+ if(fHistDistrBaseRmsOROC->GetMaximum()>fHistDistrBaseRmsIROC->GetMaximum()) fHistDistrBaseRmsIROC->SetMaximum(fHistDistrBaseRmsOROC->GetMaximum()*1.1);
+ legio->Draw("same");
- cmax->Update();
- csum->Update();
- cbasemean->Update();
- cbaserms->Update();
- }
+ cmax->Update();
+ csum->Update();
+ cbasemean->Update();
+ cbaserms->Update();
+ }
else if(histos==3)
- {
+ {
// GLOBAL MAX ADC _________________________________
- if(GetProcNextEvent()==1)
- {
- TCanvas* cglobA =0;
- TCanvas* cglobC =0;
-
- if(!(cglobC=(TCanvas*)gROOT->GetListOfCanvases()->FindObject("SIDE C all"))) cglobC = CreateCanvas("SIDE C all");
- if(!(cglobA=(TCanvas*)gROOT->GetListOfCanvases()->FindObject("SIDE A all"))) cglobA = CreateCanvas("SIDE A all");
-
- Char_t globtitle1[256]; sprintf(globtitle1,"SIDE A Run %05i (EventID %i)",fRunId,fEventNumber);
- Char_t globtitle2[256]; sprintf(globtitle2,"SIDE C Run %05i (EventID %i)",fRunId,fEventNumber);
-
- fHistGlobalMaxA->SetTitle(globtitle1);
- fHistGlobalMaxC->SetTitle(globtitle2);
- fHistGlobalMaxA->SetXTitle("x/mm");
- fHistGlobalMaxA->SetYTitle("y/mm");
- fHistGlobalMaxC->SetXTitle("x/mm");
- fHistGlobalMaxC->SetYTitle("y/mm");
-
- if(GetPedestals()==0) { fHistGlobalMaxA->SetZTitle("max adc (not baseline sub)"); fHistGlobalMaxC->SetZTitle("max adc (not baseline sub)"); }
- else { fHistGlobalMaxA->SetZTitle("max adc "); fHistGlobalMaxC->SetZTitle("max adc "); }
-
- fHistGlobalMaxA->SetMinimum(0.01);
- fHistGlobalMaxC->SetMinimum(0.01);
-
- cglobC->cd() ; fHistGlobalMaxC->Draw("COLZ");
- cglobA->cd() ; fHistGlobalMaxA->Draw("COLZ");
-
- Char_t nameom[256];
- sprintf(nameom,".x %s/TPC/AliTPCMonitorExec.C(3)",gSystem->Getenv("ALICE_ROOT"));
-
- if(fExecGlob==0)
- {
- if(fVerb)cout << " set exec " << nameom << endl;
- cglobC->AddExec("glob",nameom);
- cglobA->AddExec("glob",nameom);
- fExecGlob = 1;
- }
- else
- {
- cglobC->DeleteExec("glob");
- cglobA->DeleteExec("glob");
-
- if(fVerb) cout << " set exec " << nameom << endl;
- cglobC->AddExec("glob",nameom);
- cglobA->AddExec("glob",nameom);
-
- }
- cglobC->Update();
- cglobA->Update();
- }
+ if(GetProcNextEvent()==1)
+ {
+ TCanvas* cglobA =0;
+ TCanvas* cglobC =0;
+
+ if(!(cglobC=(TCanvas*)gROOT->GetListOfCanvases()->FindObject("SIDE C all"))) cglobC = CreateCanvas("SIDE C all");
+ if(!(cglobA=(TCanvas*)gROOT->GetListOfCanvases()->FindObject("SIDE A all"))) cglobA = CreateCanvas("SIDE A all");
+
+ Char_t globtitle1[256]; sprintf(globtitle1,"SIDE A Run %05i (EventID %i)",fRunId,fEventNumber);
+ Char_t globtitle2[256]; sprintf(globtitle2,"SIDE C Run %05i (EventID %i)",fRunId,fEventNumber);
+
+ fHistGlobalMaxA->SetTitle(globtitle1);
+ fHistGlobalMaxC->SetTitle(globtitle2);
+ fHistGlobalMaxA->SetXTitle("x/mm");
+ fHistGlobalMaxA->SetYTitle("y/mm");
+ fHistGlobalMaxC->SetXTitle("x/mm");
+ fHistGlobalMaxC->SetYTitle("y/mm");
+
+ if(GetPedestals()==0) { fHistGlobalMaxA->SetZTitle("max adc (not baseline sub)"); fHistGlobalMaxC->SetZTitle("max adc (not baseline sub)"); }
+ else { fHistGlobalMaxA->SetZTitle("max adc "); fHistGlobalMaxC->SetZTitle("max adc "); }
+
+ fHistGlobalMaxA->SetMinimum(0.01);
+ fHistGlobalMaxC->SetMinimum(0.01);
+
+ cglobC->cd() ; fHistGlobalMaxC->Draw("COLZ");
+ cglobA->cd() ; fHistGlobalMaxA->Draw("COLZ");
+ Char_t nameom[256];
+ sprintf(nameom,".x %s/TPC/AliTPCMonitorExec.C(3)",gSystem->Getenv("ALICE_ROOT"));
+
+ if(fExecGlob==0)
+ {
+ if(fVerb)cout << " set exec " << nameom << endl;
+ cglobC->AddExec("glob",nameom);
+ cglobA->AddExec("glob",nameom);
+ fExecGlob = 1;
+ }
+ else
+ {
+ cglobC->DeleteExec("glob");
+ cglobA->DeleteExec("glob");
+
+ if(fVerb) cout << " set exec " << nameom << endl;
+ cglobC->AddExec("glob",nameom);
+ cglobA->AddExec("glob",nameom);
+
+ }
+ cglobC->Update();
+ cglobA->Update();
}
+
+ }
}
//__________________________________________________________________
void AliTPCMonitor::DrawRMSMap()
{
- // Draw 2Dim rms histos for IROC and OROC
+ // Draw 2Dim rms histos for IROC and OROC
// and set executables for canvases
TCanvas* crmsoroc =0;
TCanvas* crmsiroc =0;
if(!(crmsoroc = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("crmsoroc"))) crmsoroc = CreateCanvas("crmsoroc");
if(!(crmsiroc = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("crmsiroc"))) crmsiroc = CreateCanvas("crmsiroc");
-
+
crmsiroc->cd(); fHistIROCRMS->Draw("COLZ");
crmsoroc->cd(); fHistOROCRMS->Draw("COLZ");
-
+
Char_t carry1[100]; sprintf(carry1,".x %s/TPC/AliTPCMonitorExec.C(1)",gSystem->Getenv("ALICE_ROOT"));
Char_t carry2[100]; sprintf(carry2,".x %s/TPC/AliTPCMonitorExec.C(2)",gSystem->Getenv("ALICE_ROOT"));
-
- if(fExecPadIrocRms==0)
- {
- crmsiroc->AddExec("pad",carry1);
- crmsiroc->AddExec("row",carry2);
- fExecPadIrocRms=1;
- }
- if(fExecPadOrocRms==0)
- {
- crmsoroc->AddExec("pad",carry1);
- crmsoroc->AddExec("row",carry2);
- fExecPadOrocRms=1;
- }
+ if(fExecPadIrocRms==0)
+ {
+ crmsiroc->AddExec("pad",carry1);
+ crmsiroc->AddExec("row",carry2);
+ fExecPadIrocRms=1;
+ }
+
+ if(fExecPadOrocRms==0)
+ {
+ crmsoroc->AddExec("pad",carry1);
+ crmsoroc->AddExec("row",carry2);
+ fExecPadOrocRms=1;
+ }
crmsiroc->Update();
crmsoroc->Update();
- DrawHists(2);
-
+ DrawHists(2);
+
}
//__________________________________________________________________
void AliTPCMonitor::ExecPad()
{
-
- // Executable for Pad
- // Show time profile for channel the mouse is pointing at
+
+ // Executable for Pad
+ // Show time profile for channel the mouse is pointing at
Int_t event = gPad->GetEvent();
if (event != 51) return;
Float_t upx = gPad->AbsPixeltoX(px);
Float_t y = gPad->PadtoY(upy);
Float_t x = gPad->PadtoX(upx);
-
+
Int_t setrange = 0;
TCanvas* cpad = 0;
Char_t projhist[60];
Char_t namesel[256];
Char_t namecanv[256];
-
+
Int_t xbinmin = 0;
Int_t xbinmax = 0;
Float_t ybinmin = 0;
Float_t ybinmax = 0;
Int_t rocid = 0;
-
- // Check wich Canvas executed the event
+ // Check wich Canvas executed the event
TH2S* fHistIndex=0;
sprintf(namesel,select->GetName());
- if(strcmp(namesel,"fHistOROC")==0 || strcmp(namesel,"fHistOROCRMS")==0 || strcmp(namesel,"fHistOROCTime")==0 )
- {
- rocid = 1;
- fPadUsedRoc =1;
- sprintf(projhist,"ProjectionOROC");
- sprintf(namecanv,"coroc_ch");
- fHistIndex = fHistOROCIndex;
- }
- if(strcmp(namesel,"fHistIROC")==0 || strcmp(namesel,"fHistIROCRMS")==0 || strcmp(namesel,"fHistIROCTime")==0 )
- {
- rocid = 0;
- fPadUsedRoc=0;
- sprintf(projhist,"ProjectionIROC");
- sprintf(namecanv,"ciroc_ch");
- fHistIndex = fHistIROCIndex;
- }
+ if(strcmp(namesel,"fHistOROC")==0 || strcmp(namesel,"fHistOROCRMS")==0 || strcmp(namesel,"fHistOROCTime")==0 )
+ {
+ rocid = 1;
+ fPadUsedRoc =1;
+ sprintf(projhist,"ProjectionOROC");
+ sprintf(namecanv,"coroc_ch");
+ fHistIndex = fHistOROCIndex;
+ }
+ if(strcmp(namesel,"fHistIROC")==0 || strcmp(namesel,"fHistIROCRMS")==0 || strcmp(namesel,"fHistIROCTime")==0 )
+ {
+ rocid = 0;
+ fPadUsedRoc=0;
+ sprintf(projhist,"ProjectionIROC");
+ sprintf(namecanv,"ciroc_ch");
+ fHistIndex = fHistIROCIndex;
+ }
- // Check if Canvas already existed and get Ranges from former Prjection histogram
- if((cpad=(TCanvas*)gROOT->GetListOfCanvases()->FindObject(namecanv)))
- {
- cpad->cd();
- if(gROOT->Get(projhist))
- {
- setrange = 1;
- xbinmin = ((TH1D*)gROOT->Get(projhist))->GetXaxis()->GetFirst();
- xbinmax = ((TH1D*)gROOT->Get(projhist))->GetXaxis()->GetLast();
- ybinmin = ((TH1D*)gROOT->Get(projhist))->GetMinimum();
- ybinmax = ((TH1D*)gROOT->Get(projhist))->GetMaximum();
- delete gROOT->Get("legfit");
- delete gROOT->Get("fg");
- }
- }
- else
+ // Check if Canvas already existed and get Ranges from former Prjection histogram
+ if((cpad=(TCanvas*)gROOT->GetListOfCanvases()->FindObject(namecanv)))
+ {
+ cpad->cd();
+ if(gROOT->Get(projhist))
{
- cpad = CreateCanvas(namecanv); cpad->cd();
+ setrange = 1;
+ xbinmin = ((TH1D*)gROOT->Get(projhist))->GetXaxis()->GetFirst();
+ xbinmax = ((TH1D*)gROOT->Get(projhist))->GetXaxis()->GetLast();
+ ybinmin = ((TH1D*)gROOT->Get(projhist))->GetMinimum();
+ ybinmax = ((TH1D*)gROOT->Get(projhist))->GetMaximum();
+ delete gROOT->Get("legfit");
+ delete gROOT->Get("fg");
}
+ }
+ else
+ {
+ cpad = CreateCanvas(namecanv); cpad->cd();
+ }
- // Get Bin
+ // Get Bin
Int_t testy = fHistIndex->GetYaxis()->FindBin(y);
Int_t testx = fHistIndex->GetXaxis()->FindBin(x);
Int_t binchannel = (Int_t)fHistIndex->GetCellContent(testx,testy);
if(binchannel>30000 || binchannel<0) return;
-
- if(gROOT->Get(projhist)) delete gROOT->Get(projhist);
- // Get Projection
- TH1D *hp = (TH1D*)(((TH1D*)fHistChannelTime->ProjectionY("hp",binchannel,binchannel))->Clone(projhist));
+
+ TH1D *hp=(TH1D*)gROOT->Get(projhist);
+ if(!hp) hp=new TH1D(projhist,projhist,GetTimeBins(),0,GetTimeBins());//delete gROOT->Get(projhist);
+ // Get Projection
+// TH1D *hp = (TH1D*)(((TH1D*)fHistChannelTime->ProjectionY("hp",binchannel,binchannel))->Clone(projhist));
+// TH1D *hp = fHistChannelTime->ProjectionY(projhist,binchannel,binchannel);
+// hp->Reset();
+ Int_t nbinsx=fHistChannelTime->GetNbinsX();
+ for (Int_t itb=0;itb<GetTimeBins();++itb){
+ hp->GetArray()[itb]=fHistChannelTime->GetArray()[binchannel+itb*(nbinsx+2)];
+ }
+// printf("hi\n");
+// return;
// Make title and Pave for channel Info
Char_t title[256];
nhw = hwadd;
nmax = (Int_t)hp->GetMaximum();
-
+
TPaveText* legstat = new TPaveText(0.18,0.65,0.3,0.8,"NDC");
-
+
Int_t connector = fMapHand->GetFECconnector(hwadd);
Int_t fecnr = fMapHand->GetFECfromHw(hwadd);
Int_t fecch = fMapHand->GetFECchannel(hwadd);
Int_t fecloc = fMapHand->U2fGetFECinRCU(fecnr) ;
Int_t feclocbran = fMapHand->U2fGetFECinBranch(fecnr);
Int_t branch = fMapHand->U2fGetBranch(fecnr);
-
- Short_t fecget = (hwadd & AliTPCMonitorAltro::GetHwMaskFEC()) >> 7;
- Short_t branchget = (hwadd & AliTPCMonitorAltro::GetHwMaskBranch())>> 11;
-
-
- Char_t nstat1[100]; Char_t nstat2[100]; Char_t nstat3[100]; Char_t nstat4[100];
+
+ Short_t fecget = (hwadd & fgkHwMaskFEC) >> 7;
+ Short_t branchget = (hwadd & fgkHwMaskBranch)>> 11;
+
+
+ Char_t nstat1[100]; Char_t nstat2[100]; Char_t nstat3[100]; Char_t nstat4[100];
Char_t nstat5[100]; Char_t nstat6[100]; Char_t nstat7[100]; Char_t nstat8[100];
sprintf(nstat1,"Branch (map) \t %i (%i) \n",branchget,branch);
sprintf(nstat5,"Altro chip\t %i \n",altrochip);
sprintf(nstat6,"Altro chan\t %i \n",altrochannel);
- legstat->AddText(nstat1); legstat->AddText(nstat2); legstat->AddText(nstat8); legstat->AddText(nstat7);
+ legstat->AddText(nstat1); legstat->AddText(nstat2); legstat->AddText(nstat8); legstat->AddText(nstat7);
legstat->AddText(nstat3); legstat->AddText(nstat4); legstat->AddText(nstat5); legstat->AddText(nstat6);
-
+
sprintf(title,"Row=%d Pad=%d Hw =%d maxADC =%d count =%d",npadRow,npad,nhw,nmax,binchannel);
- hp->SetName(projhist);
+// hp->SetName(projhist);
hp->SetTitleSize(0.04);
hp->SetTitle(title);
hp->SetYTitle("ADC");
hp->SetXTitle("Timebin");
hp->GetXaxis()->SetTitleColor(1);
- if(setrange)
- {
- hp->GetXaxis()->SetRange(xbinmin,xbinmax);
- hp->SetMinimum(ybinmin);
- hp->SetMaximum(ybinmax);
- }
- else
- {
- hp->SetMinimum(0.0);
- hp->SetMaximum(1000.0);
- }
+ if(setrange)
+ {
+ hp->GetXaxis()->SetRange(xbinmin,xbinmax);
+ hp->SetMinimum(ybinmin);
+ hp->SetMaximum(ybinmax);
+ }
+ else
+ {
+ hp->SetMinimum(0.0);
+ hp->SetMaximum(1000.0);
+ }
cpad->cd();
hp->Draw();
// Make Fit to peak
- if(GetPedestals() && fDisableFit==0)
+ if(GetPedestals() && fDisableFit==0)
+ {
+ Int_t maxx = (Int_t)fHistAddrMaxAdcX->GetBinContent(hwadd);
+ Float_t max = (Float_t)fHistAddrMaxAdc->GetBinContent(hwadd);
+ Float_t base = (Float_t)fHistAddrBaseMean->GetBinContent(hwadd);
+ if(base!=0)
{
- Int_t maxx = (Int_t)fHistAddrMaxAdcX->GetBinContent(hwadd);
- Float_t max = (Float_t)fHistAddrMaxAdc->GetBinContent(hwadd);
- Float_t base = (Float_t)fHistAddrBaseMean->GetBinContent(hwadd);
- if(base!=0)
- {
- if( ((max+base)/base)>1.2)
- {
- TF1* fg = new TF1("fg",AliTPCMonitor::Gamma4,maxx-5,maxx+5,4);
- fg->SetParName(0,"Normalisation");
- fg->SetParName(1,"Minimum");
- fg->SetParName(2,"Width");
- fg->SetParName(3,"Base");
- fg->SetParameter(0,max);
- fg->SetParameter(1,maxx-2);
- fg->SetParameter(2,1.5);
- fg->FixParameter(3,0);
- fg->SetLineColor(4);
- fg->SetLineWidth(1);
- hp->Fit("fg","RQ");
-
- TLegend* legfit = new TLegend(0.6,0.7,0.7,0.8);
- legfit->AddEntry("fg","#Gamma 4 fit","l");
- legfit->SetFillColor(0);
- legfit->SetName("legfit");
- legfit->Draw("same");
- }
- }
+ if( ((max+base)/base)>1.2)
+ {
+ TF1* fg = new TF1("fg",AliTPCMonitor::Gamma4,maxx-5,maxx+5,4);
+ fg->SetParName(0,"Normalisation");
+ fg->SetParName(1,"Minimum");
+ fg->SetParName(2,"Width");
+ fg->SetParName(3,"Base");
+ fg->SetParameter(0,max);
+ fg->SetParameter(1,maxx-2);
+ fg->SetParameter(2,1.5);
+ fg->FixParameter(3,0);
+ fg->SetLineColor(4);
+ fg->SetLineWidth(1);
+ hp->Fit("fg","RQ");
+
+ TLegend* legfit = new TLegend(0.6,0.7,0.7,0.8);
+ legfit->AddEntry("fg","#Gamma 4 fit","l");
+ legfit->SetFillColor(0);
+ legfit->SetName("legfit");
+ legfit->Draw("same");
+ }
}
+ }
legstat->SetFillColor(0);
legstat->Draw("same");
cpad->Update();
//__________________________________________________________________
void AliTPCMonitor::ExecRow()
{
-
- // Executable for Pad
- // Show profile of max adc over given pad row
- // and 2dim histo adc(pad-in-row,time bin)
-
+
+ // Executable for Pad
+ // Show profile of max adc over given pad row
+ // and 2dim histo adc(pad-in-row,time bin)
+
Int_t event = gPad->GetEvent();
if (event != 61) return;
gPad->cd();
TObject *select = gPad->GetSelected();
if(!select) return;
if(!select->InheritsFrom("TH2")) { return; }
-
+
Int_t rocid = 0;
// Char_t namehist[50];
Char_t rowhist[60];
sprintf(rowhistxmax ,"hxmax");
sprintf(rowhistmax , "hrowmax");
- // get position
+ // get position
Int_t px = gPad->GetEventX();
Int_t py = gPad->GetEventY();
Float_t upy = gPad->AbsPixeltoY(py);
TCanvas*cxmax = 0;
TH2S* fHistIndex = 0;
-
- // ranges from already existing histos
+
+ // ranges from already existing histos
Int_t rowtimexmin = 0;
Int_t rowtimexmax = 0;
Int_t rowtimeymin = 0;
Int_t profrowxmax = 0;
Double_t profrowymin = 0;
Double_t profrowymax = 0;
-
+
Int_t profxxmin = 0;
Int_t profxxmax = 0;
Double_t profxymin = 0;
Double_t profxymax = 0;
-
-
+
+
Int_t setrange = 0;
-
-
- if( strcmp(select->GetName(),"fHistIROC")==0 || strcmp(select->GetName(),"fHistIROCRMS")==0 ) { fHistIndex = fHistIROCIndex; rocid =1; }
+
+
+ if( strcmp(select->GetName(),"fHistIROC")==0 || strcmp(select->GetName(),"fHistIROCRMS")==0 ) { fHistIndex = fHistIROCIndex; rocid =1; }
else if(strcmp(select->GetName(),"fHistOROC")==0 || strcmp(select->GetName(),"fHistOROCRMS")==0 ) { fHistIndex = fHistOROCIndex; rocid =2; }
else { cout << " not implemented for this histo " << endl; return; }
gPad->GetCanvas()->FeedbackMode(kTRUE);
-
-
+
+
// check if canvases exist //
crowtime = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("crowtime");
crowmax = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("crowmax");
if(!crowmax) crowmax = CreateCanvas("crowmax") ;
if(!cxmax ) cxmax = CreateCanvas("cxmax") ;
- // check ranges of already existing histos
- if(gROOT->Get(rowhist))
- {
- rowtimexmin = ((TH2F*)gROOT->Get(rowhist))->GetXaxis()->GetFirst();
- rowtimexmax = ((TH2F*)gROOT->Get(rowhist))->GetXaxis()->GetLast();
- rowtimeymin = ((TH2F*)gROOT->Get(rowhist))->GetYaxis()->GetFirst();
- rowtimeymax = ((TH2F*)gROOT->Get(rowhist))->GetYaxis()->GetLast();
- rowtimezmin = ((TH2F*)gROOT->Get(rowhist))->GetMinimum();
- rowtimezmax = ((TH2F*)gROOT->Get(rowhist))->GetMaximum();
-
- profrowxmin = ((TH1F*)gROOT->Get(rowhistmax))->GetXaxis()->GetFirst();
- profrowxmax = ((TH1F*)gROOT->Get(rowhistmax))->GetXaxis()->GetLast();
- profrowymin = ((TH1F*)gROOT->Get(rowhistmax))->GetMinimum();
- profrowymax = ((TH1F*)gROOT->Get(rowhistmax))->GetMaximum();
-
- profxxmin = ((TH1F*)gROOT->Get(rowhistxmax))->GetXaxis()->GetFirst();
- profxxmax = ((TH1F*)gROOT->Get(rowhistxmax))->GetXaxis()->GetLast();
- profxymin = ((TH1F*)gROOT->Get(rowhistxmax))->GetMinimum();
- profxymax = ((TH1F*)gROOT->Get(rowhistxmax))->GetMaximum();
-
- setrange =1;
-
- delete gROOT->Get(rowhist);
- delete gROOT->Get(rowhistmax);
- delete gROOT->Get(rowhistsum);
- delete gROOT->Get("hxmax");
- delete gROOT->Get("legrow");
- }
+ // check ranges of already existing histos
+ if(gROOT->Get(rowhist))
+ {
+ rowtimexmin = ((TH2F*)gROOT->Get(rowhist))->GetXaxis()->GetFirst();
+ rowtimexmax = ((TH2F*)gROOT->Get(rowhist))->GetXaxis()->GetLast();
+ rowtimeymin = ((TH2F*)gROOT->Get(rowhist))->GetYaxis()->GetFirst();
+ rowtimeymax = ((TH2F*)gROOT->Get(rowhist))->GetYaxis()->GetLast();
+ rowtimezmin = ((TH2F*)gROOT->Get(rowhist))->GetMinimum();
+ rowtimezmax = ((TH2F*)gROOT->Get(rowhist))->GetMaximum();
+
+ profrowxmin = ((TH1F*)gROOT->Get(rowhistmax))->GetXaxis()->GetFirst();
+ profrowxmax = ((TH1F*)gROOT->Get(rowhistmax))->GetXaxis()->GetLast();
+ profrowymin = ((TH1F*)gROOT->Get(rowhistmax))->GetMinimum();
+ profrowymax = ((TH1F*)gROOT->Get(rowhistmax))->GetMaximum();
+
+ profxxmin = ((TH1F*)gROOT->Get(rowhistxmax))->GetXaxis()->GetFirst();
+ profxxmax = ((TH1F*)gROOT->Get(rowhistxmax))->GetXaxis()->GetLast();
+ profxymin = ((TH1F*)gROOT->Get(rowhistxmax))->GetMinimum();
+ profxymax = ((TH1F*)gROOT->Get(rowhistxmax))->GetMaximum();
+
+ setrange =1;
+
+ delete gROOT->Get(rowhist);
+ delete gROOT->Get(rowhistmax);
+ delete gROOT->Get(rowhistsum);
+ delete gROOT->Get("hxmax");
+ delete gROOT->Get("legrow");
+ }
// get channel for xy bin -> getRow -> getNrows -> getHw for each Pad in Row -> get channel for each hw -> make proj
Int_t testy = fHistIndex->GetYaxis()->FindBin(y);
if(binchannel>30000) return;
if(binchannel<=0 ) { crowtime->Update() ; crowmax->Update() ; return ; }
- // get hwaddress
+ // get hwaddress
Int_t hwadd = (Int_t)fHistChannelTime->GetCellContent(binchannel,0);
Int_t row = fMapHand->GetPadRow(hwadd);
Int_t pad = fMapHand->GetPad(hwadd) ;
Int_t numofpads = fMapHand->GetNumofPads(row);
- // create histos
+ // create histos
TH2F *hrowtime = new TH2F(rowhist , "" ,numofpads,0,numofpads,GetTimeBins(),0.0,GetTimeBins());
TH1F *hrowmax = new TH1F(rowhistmax , "" ,numofpads,0,numofpads);
TH1F *hxmax = new TH1F(rowhistxmax, "" ,159,0,159 );
// Row profile ///////////
if(fVerb) cout << " Number of pads " << numofpads << endl;
- for(Int_t padnr = 0; padnr<numofpads;padnr++)
- {
- Int_t addrinrow = fMapHand->GetPadAddInRow(row,padnr );
- Int_t channel = (Int_t)fHistAddrMapIndex->GetBinContent(addrinrow);
- if(channel==-1) continue;
+ for(Int_t padnr = 0; padnr<numofpads;padnr++)
+ {
+ Int_t addrinrow = fMapHand->GetPadAddInRow(row,padnr );
+ Int_t channel = (Int_t)fHistAddrMapIndex->GetBinContent(addrinrow);
+ if(channel==-1) continue;
+
+ hrowmax->SetBinContent(padnr+1,fHistAddrMaxAdc->GetBinContent(addrinrow));
+ TH1D *hp = fHistChannelTime->ProjectionY("hp",channel,channel);
+ for(Int_t time = 0;time<GetTimeBins();time++) {
- hrowmax->SetBinContent(padnr+1,fHistAddrMaxAdc->GetBinContent(addrinrow));
- TH1D *hp = fHistChannelTime->ProjectionY("hp",channel,channel);
- for(Int_t time = 0;time<GetTimeBins();time++) {
-
- Float_t val = hp->GetBinContent(time);
- hrowtime->SetCellContent(padnr+1,time+1,val);
- }
+ Float_t val = hp->GetBinContent(time);
+ hrowtime->SetCellContent(padnr+1,time+1,val);
}
+ }
// X profile /////////////
Double_t xval = 0.0;
Double_t yval = 0.0;
GetXY(xval,yval,numofpads,row,pad);
- Int_t padnr = 0;
+ Int_t padnr = 0;
Int_t hw = 0;
for(Int_t nrow = 0; nrow<159; nrow++)
+ {
+ padnr = GetPadAtX(xval,nrow);
+ if(padnr>=0)
{
- padnr = GetPadAtX(xval,nrow);
- if(padnr>=0)
- {
- hw = fMapHand->GetPadAddInRow(nrow,padnr);
- if(fPadMapHw[hw]==-1){ continue ; }
- else { hxmax->SetBinContent(nrow+1,fHistAddrMaxAdc->GetBinContent(hw)) ; }
- }
+ hw = fMapHand->GetPadAddInRow(nrow,padnr);
+ if(fPadMapHw[hw]==-1){ continue ; }
+ else { hxmax->SetBinContent(nrow+1,fHistAddrMaxAdc->GetBinContent(hw)) ; }
}
+ }
cxmax->cd();
Char_t hxtitle[50] ; sprintf(hxtitle,"max adc in pads at x=%5.1f mm",xval);
hxmax->Draw("l");
if(setrange)
- {
- hxmax->GetXaxis()->SetRange(profxxmin,profxxmax);
- hxmax->SetMinimum(profxymin);
- hxmax->SetMaximum(profxymax);
- }
+ {
+ hxmax->GetXaxis()->SetRange(profxxmin,profxxmax);
+ hxmax->SetMinimum(profxymin);
+ hxmax->SetMaximum(profxymax);
+ }
cxmax->Update();
-
+
crowtime->cd();
Char_t title[256];
Char_t titlemax[256];
- if(rocid==1) {sprintf(title,"%s Row=%d",((TH2*)select)->GetTitle(),row) ; sprintf(titlemax,"IROC max/sum Row=%d",row );}
+ if(rocid==1) {sprintf(title,"%s Row=%d",((TH2*)select)->GetTitle(),row) ; sprintf(titlemax,"IROC max/sum Row=%d",row );}
else {sprintf(title,"%s Row=%d",((TH2*)select)->GetTitle(),row-63); sprintf(titlemax,"OROC max/sum Row=%d",row-63);}
if(fVerb) cout << " set name " << endl;
-
+
// row vs time
crowtime->cd();
hrowtime->SetTitleSize(0.04);
hrowtime->SetYTitle("timbin");
hrowtime->SetXTitle("pad in row");
hrowtime->SetZTitle("signal (ADC)");
-
+
hrowtime->GetXaxis()->SetTitleColor(1);
hrowtime->SetMaximum(1000.0);
hrowtime->SetMinimum(0.0);
- if(setrange)
- {
- hrowtime->GetXaxis()->SetRange(rowtimexmin,rowtimexmax);
- hrowtime->GetYaxis()->SetRange(rowtimeymin,rowtimeymax);
- hrowtime->SetMinimum(rowtimezmin);
- hrowtime->SetMaximum(rowtimezmax);
- }
-
+ if(setrange)
+ {
+ hrowtime->GetXaxis()->SetRange(rowtimexmin,rowtimexmax);
+ hrowtime->GetYaxis()->SetRange(rowtimeymin,rowtimeymax);
+ hrowtime->SetMinimum(rowtimezmin);
+ hrowtime->SetMaximum(rowtimezmax);
+ }
+
hrowtime->Draw("COLZ");
crowtime->UseCurrentStyle();
crowtime->Update();
-
+
// max and sum /////////////////////////
crowmax->cd();
if(setrange) {
hrowmax->SetYTitle("max adc");
hrowmax->SetXTitle("pad in row");
hrowmax->GetXaxis()->SetTitleColor(1);
-
+
hrowmax->SetLineColor(2);
- hrowmax->Draw("l");
+ hrowmax->Draw("l");
crowmax->Update();
return;
void AliTPCMonitor::Write10bitChannel()
{
- // Write 10 bit words form histogram for active(last pointed) channel
+ // Write 10 bit words form histogram for active(last pointed) channel
if(fPadUsedHwAddr==-1){ AliWarning(" No active pad "); return ;}
- Int_t pad = (Int_t)fMapHand->GetPad( fPadUsedHwAddr);
- Int_t row = (Int_t)fMapHand->GetPadRow(fPadUsedHwAddr);
+ Int_t pad = (Int_t)fMapHand->GetPad( fPadUsedHwAddr);
+ Int_t row = (Int_t)fMapHand->GetPadRow(fPadUsedHwAddr);
Int_t channel = (Int_t)fPadMapHw[fPadUsedHwAddr];
Char_t filenameroot[256];
Char_t filenamedat[256];
Char_t projhist[256];
-
+
if(fPadUsedRoc==1) { sprintf(projhist,"ProjectionOROC"); }
if(fPadUsedRoc==0) { sprintf(projhist,"ProjectionIROC"); }
sprintf(filenamedat, "Channel_Run%05i_EventID_%i_Pad_%i_Row_%i.dat" ,fRunId,fEventNumber,pad,row);
sprintf(filenameroot,"Channel_Run%05i_EventID_%i_Pad_%i_Row_%i.root" ,fRunId,fEventNumber,pad,row);
- TH1D* hpr = 0;
- if((hpr=(TH1D*)gROOT->Get(projhist)))
+ TH1D* hpr = 0;
+ if((hpr=(TH1D*)gROOT->Get(projhist)))
+ {
+ // root histo
+ TFile f(filenameroot,"recreate");
+ hpr->Write();
+ f.Close();
+
+ // raw singal
+ ofstream datout(filenamedat,ios::out);
+ datout <<"Timebin \t ADC value " << endl;
+ for(Int_t i = 1; i <GetTimeBins(); i++)
{
- // root histo
- TFile f(filenameroot,"recreate");
- hpr->Write();
- f.Close();
-
- // raw singal
- ofstream datout(filenamedat,ios::out);
- datout <<"Timebin \t ADC value " << endl;
- for(Int_t i = 1; i <GetTimeBins(); i++)
- {
- datout << i << " \t \t " << fPad[channel][i] << endl;
- }
- datout.close();
+ datout << i << " \t \t " << fPad[channel][i] << endl;
}
+ datout.close();
+ }
else
- {
- AliWarning("No projection histo found ");
- }
+ {
+ AliWarning("No projection histo found ");
+ }
}
//__________________________________________________________________
void AliTPCMonitor::ExecTransform()
{
-
+
// Make Fast Fourier Transformation for active pad
// fft is only performed for a data sample of size 2^n
- // reduce window according to largest power of 2 which is smaller than the viewing range
-
+ // reduce window according to largest power of 2 which is smaller than the viewing range
+
Char_t namecanv[256];
Char_t namecanv2[256];
Char_t projhist[256];
if(fPadUsedRoc==1) { sprintf(namecanv,"coroc_ch_trans") ; sprintf(namecanv2,"coroc_ch_trans2") ; sprintf(projhist,"ProjectionOROC"); }
if(fPadUsedRoc==0) { sprintf(namecanv,"ciroc_ch_trans") ; sprintf(namecanv2,"ciroc_ch_trans2") ; sprintf(projhist,"ProjectionIROC"); }
-
+
TH1D* hproj = 0;
if((TH1D*)gROOT->Get(projhist)==0){AliWarning("Proj histo does not exist \n Move mouse over 2d histo choose channel \n and drag mouse form histo again!"); return ;}
else hproj = (TH1D*)gROOT->Get(projhist) ;
-
+
if(fPadUsedRoc==1) { sprintf(namehtrimag,"htransimagfreq_oroc"); sprintf(namehtrreal,"htransrealfreq_oroc"); sprintf(namehtrmag,"htransmagfreq_oroc"); }
else { sprintf(namehtrimag,"htransimagfreq_iroc"); sprintf(namehtrreal,"htransrealfreq_iroc"); sprintf(namehtrmag,"htransmagfreq_iroc"); }
-
+
if( gROOT->Get(namehtrimag)) delete gROOT->Get(namehtrimag);
if( gROOT->Get(namehtrreal)) delete gROOT->Get(namehtrreal);
if( gROOT->Get(namehtrmag)) delete gROOT->Get(namehtrmag);
-
+
TCanvas *ctrans = 0;
if(!(ctrans = (TCanvas*)gROOT->GetListOfCanvases()->FindObject(namecanv)))
{
ctrans2 = CreateCanvas(namecanv2);
// ctrans2->Divide(1,2);
}
-
+
Int_t binfirst = hproj->GetXaxis()->GetFirst();
Int_t binlast = hproj->GetXaxis()->GetLast();
Int_t bins = binlast -binfirst +1;
-
+
Int_t power = 0;
for(Int_t pot = 0; pot<=10 ; pot++)
{
Int_t comp = (Int_t)TMath::Power(2,pot);
if(bins>=comp)power = pot;
}
-
+
bins = (Int_t)TMath::Power(2,power);
-
+
// sampling frequency ;
Double_t deltat = 1.0/(Float_t)GetSamplingFrequency();
-
+
// output histo
TH1D* htransrealfreq = new TH1D(namehtrreal,namehtrreal,10000,-1/(2*deltat),1/(2*deltat));
TH1D* htransimagfreq = new TH1D(namehtrimag,namehtrimag,10000,-1/(2*deltat),1/(2*deltat));
TH1D* htransmag = new TH1D(namehtrmag,namehtrmag,10000,-1/(2*deltat),1/(2*deltat));
-
+
Char_t titlereal[256];
Char_t titleimag[256];
Char_t titlemag[256];
if(fPadUsedRoc==1) { sprintf(titlereal,"OROC DFT real part"); sprintf(titleimag,"OROC DFT imag part"); sprintf(titlemag,"OROC DFT magnitude"); }
else { sprintf(titlereal,"IROC DFT real part"); sprintf(titleimag,"IROC DFT imag part"); sprintf(titlemag,"IROC DFT magnitude"); }
-
+
htransrealfreq->SetTitle(titlereal); htransrealfreq->SetXTitle("f/hz"); htransrealfreq->SetYTitle("z_{real}(f)");
htransimagfreq->SetTitle(titleimag); htransimagfreq->SetXTitle("f/hz"); htransimagfreq->SetYTitle("z_{imag}(f)");
htransmag->SetTitle(titlemag); htransmag->SetXTitle("f/hz"); htransmag->SetYTitle("mag(f)");
-
- // create complex packed data array
+
+ // create complex packed data array
const Int_t kdatasiz = 2*bins;
Double_t* data = new Double_t[kdatasiz];
for(Int_t i=0;i<2*bins;i++) { data[i] = 0.0;}
for(Int_t i=0;i<bins;i++) { data[2*i] = (Double_t)hproj->GetBinContent(binfirst+i); }
-
+
// make fourier transformation
AliTPCMonitorFFT* four = new AliTPCMonitorFFT();
four->ComplexRadix2ForwardWrap(data,1,bins);
-
- // write output and fill histos forward
+
+ // write output and fill histos forward
Double_t freq = 0.0;
for(Int_t i=0;i<2*bins;i++)
{
htransmag->Fill( freq, TMath::Sqrt(data[2*i]*data[2*i]+data[2*i+1]*data[2*i+1]) );
}
}
-
+
ctrans->cd(1);
htransrealfreq->Draw();
ctrans->cd(2);
htransmag->Draw();
ctrans2->Update();
delete four;
- delete data;
+ delete [] data;
}
//__________________________________________________________________
void AliTPCMonitor::ShowSel(Int_t* compval)
{
- // Show only selected components
- // First restore original histogram from clone
- // Than remove all not matching pads form histos
+ // Show only selected components
+ // First restore original histogram from clone
+ // Than remove all not matching pads form histos
Int_t connector = 0;
Int_t fecnr = 0;
Int_t emptyI = 1;
Int_t index = -1;
Int_t hwadd = 0;
-
+
Float_t maxiroc = fHistIROCClone->GetMaximum();
Float_t maxoroc = fHistOROCClone->GetMaximum();
-
- // restore original histos
- for(Int_t row = 0; row<fkNRowsIroc; row++)
+
+ // restore original histos
+ for(Int_t row = 0; row<fkNRowsIroc; row++)
+ {
+ for(Int_t pad = 0; pad<fkNPadsIroc; pad++)
{
- for(Int_t pad = 0; pad<fkNPadsIroc; pad++)
- {
index = (Int_t)fHistIROCIndex->GetCellContent(row+1,pad+1);
if(index==-1)continue;
else fHistIROC->SetCellContent(row+1,pad+1,fHistIROCClone->GetCellContent(row+1,pad+1));
- }
}
- for(Int_t row = 0; row<fkNRowsOroc; row++)
+ }
+ for(Int_t row = 0; row<fkNRowsOroc; row++)
+ {
+ for(Int_t pad = 0; pad<fkNPadsOroc; pad++)
{
- for(Int_t pad = 0; pad<fkNPadsOroc; pad++)
- {
- index = (Int_t)fHistOROCIndex->GetCellContent(row+1,pad+1);
- if(index==-1)continue;
- else fHistOROC->SetCellContent(row+1,pad+1,fHistOROCClone->GetCellContent(row+1,pad+1));
- }
+ index = (Int_t)fHistOROCIndex->GetCellContent(row+1,pad+1);
+ if(index==-1)continue;
+ else fHistOROC->SetCellContent(row+1,pad+1,fHistOROCClone->GetCellContent(row+1,pad+1));
}
+ }
// remove not matching entries from fHistIROC/fHistOROC
Int_t npads =0;
Int_t subrows =0;
- for(Int_t row = 0; row< (fkNRowsIroc + fkNRowsOroc); row++)
+ for(Int_t row = 0; row< (fkNRowsIroc + fkNRowsOroc); row++)
+ {
+ if(row<fkNRowsIroc) { fHist=fHistIROC ; fHistIndex = fHistIROCIndex; npads = fkNPadsIroc; subrows =0 ;}
+ else { fHist=fHistOROC ; fHistIndex = fHistOROCIndex; npads = fkNPadsOroc; subrows =fkNRowsIroc;}
+
+ for(Int_t pad = 0; pad<npads; pad++)
{
- if(row<fkNRowsIroc) { fHist=fHistIROC ; fHistIndex = fHistIROCIndex; npads = fkNPadsIroc; subrows =0 ;}
- else { fHist=fHistOROC ; fHistIndex = fHistOROCIndex; npads = fkNPadsOroc; subrows =fkNRowsIroc;}
+ index = (Int_t)fHistIndex->GetCellContent(row -subrows +1,pad+1);
+ if(index==-1) continue ;
+ hwadd = (Int_t)fHistChannelTime->GetCellContent(index,0);
+
+ // global fecnr
+ fecnr = fMapHand->GetFECfromHw(hwadd);
+ if(compval[0]!=-1 && fecnr!=compval[0]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
+
+ // rcu
+ rcuget = (hwadd & fgkHwMaskRCU)>> 12;
+ if(compval[1]!=-1 && rcuget!=compval[1]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
+
+ // branch
+ branch = fMapHand->U2fGetBranch(fecnr) ;
+ if(compval[2]!=-1 && branch!=compval[2]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
- for(Int_t pad = 0; pad<npads; pad++)
- {
- index = (Int_t)fHistIndex->GetCellContent(row -subrows +1,pad+1);
- if(index==-1) continue ;
- hwadd = (Int_t)fHistChannelTime->GetCellContent(index,0);
-
- // global fecnr
- fecnr = fMapHand->GetFECfromHw(hwadd);
- if(compval[0]!=-1 && fecnr!=compval[0]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
-
- // rcu
- rcuget = (hwadd & AliTPCMonitorAltro::GetHwMaskRCU())>> 12;
- if(compval[1]!=-1 && rcuget!=compval[1]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
-
- // branch
- branch = fMapHand->U2fGetBranch(fecnr) ;
- if(compval[2]!=-1 && branch!=compval[2]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
-
- // local fec
- feclocbran= fMapHand->U2fGetFECinBranch(fecnr) ;
- if(compval[3]!=-1 && feclocbran!=compval[3]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
-
- // connector
- connector = fMapHand->GetFECconnector(hwadd);
- if(compval[4]!=-1 && connector!=compval[4]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
-
- // Altro chip
- altrochip = fMapHand->GetAltro(hwadd);
- if(compval[5]!=-1 && altrochip!=compval[5]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
- emptyI =0;
- }
+ // local fec
+ feclocbran= fMapHand->U2fGetFECinBranch(fecnr) ;
+ if(compval[3]!=-1 && feclocbran!=compval[3]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
+
+ // connector
+ connector = fMapHand->GetFECconnector(hwadd);
+ if(compval[4]!=-1 && connector!=compval[4]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
+
+ // Altro chip
+ altrochip = fMapHand->GetAltro(hwadd);
+ if(compval[5]!=-1 && altrochip!=compval[5]) { fHist->SetCellContent(row-subrows+1,pad+1,0); continue; }
+ emptyI =0;
}
+ }
TCanvas* c1 =0;
TCanvas* c2 =0;
- if(gROOT->GetListOfCanvases()->FindObject("ciroc"))
- {
- c1 = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("ciroc");
- c1->cd() ;
- fHistIROC->Draw("COLZ");
- fHistIROC->SetMaximum(maxiroc);
- fHistIROC->SetMinimum(0.0);
- c1->Update();
- }
- if(gROOT->GetListOfCanvases()->FindObject("coroc"))
- {
- c2 = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("coroc");
- c2->cd() ;
- fHistOROC->Draw("COLZ");
- fHistOROC->SetMaximum(maxoroc);
- fHistOROC->SetMinimum(0.0);
- c2->Update();
- }
+ if(gROOT->GetListOfCanvases()->FindObject("ciroc"))
+ {
+ c1 = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("ciroc");
+ c1->cd() ;
+ fHistIROC->Draw("COLZ");
+ fHistIROC->SetMaximum(maxiroc);
+ fHistIROC->SetMinimum(0.0);
+ c1->Update();
+ }
+ if(gROOT->GetListOfCanvases()->FindObject("coroc"))
+ {
+ c2 = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("coroc");
+ c2->cd() ;
+ fHistOROC->Draw("COLZ");
+ fHistOROC->SetMaximum(maxoroc);
+ fHistOROC->SetMinimum(0.0);
+ c2->Update();
+ }
return ;
}
{
// Resize canvases and delete some of them
- Char_t carry1[100];
+ Char_t carry1[100];
sprintf(carry1,".x %s/TPC/AliTPCMonitorExec.C(1)",gSystem->Getenv("ALICE_ROOT"));
- Char_t carry3[100];
+ Char_t carry3[100];
sprintf(carry3,".x %s/TPC/AliTPCMonitorExec.C(2)",gSystem->Getenv("ALICE_ROOT"));
if(fVerb) cout << " canv 1 " << endl;
if(gROOT->GetListOfCanvases()->FindObject( "ciroc_ch")) { delete gROOT->GetListOfCanvases()->FindObject("ciroc_ch") ; }
// for 2dim plots delete create and draw again
- if(gROOT->GetListOfCanvases()->FindObject("ciroc"))
- {
- delete gROOT->GetListOfCanvases()->FindObject("ciroc");
- TCanvas* ciroc = CreateCanvas("ciroc");
- ciroc->cd();
- fHistIROC->Draw("COLZ");
- ciroc->AddExec("pad",carry1);
- ciroc->AddExec("row",carry3);
- fExecPlaneMax=1;
- ciroc->Update();
- }
+ if(gROOT->GetListOfCanvases()->FindObject("ciroc"))
+ {
+ delete gROOT->GetListOfCanvases()->FindObject("ciroc");
+ TCanvas* ciroc = CreateCanvas("ciroc");
+ ciroc->cd();
+ fHistIROC->Draw("COLZ");
+ ciroc->AddExec("pad",carry1);
+ ciroc->AddExec("row",carry3);
+ fExecPlaneMax=1;
+ ciroc->Update();
+ }
// for 2dim plots delete create and draw again
- if(gROOT->GetListOfCanvases()->FindObject("coroc"))
- {
- delete gROOT->GetListOfCanvases()->FindObject("coroc");
- TCanvas* coroc = CreateCanvas("coroc");
- coroc->cd();
- fHistOROC->Draw("COLZ");
-
- coroc->AddExec("pad",carry1);
- coroc->AddExec("row",carry3);
- coroc->Update();
- fExecPlaneMax=1;
- }
+ if(gROOT->GetListOfCanvases()->FindObject("coroc"))
+ {
+ delete gROOT->GetListOfCanvases()->FindObject("coroc");
+ TCanvas* coroc = CreateCanvas("coroc");
+ coroc->cd();
+ fHistOROC->Draw("COLZ");
+
+ coroc->AddExec("pad",carry1);
+ coroc->AddExec("row",carry3);
+ coroc->Update();
+ fExecPlaneMax=1;
+ }
if(gROOT->GetListOfCanvases()->FindObject( "cbasemean")) { delete gROOT->GetListOfCanvases()->FindObject("cbasemean"); }
if(gROOT->GetListOfCanvases()->FindObject( "cbase")) { delete gROOT->GetListOfCanvases()->FindObject("cbase");}
if(gROOT->GetListOfCanvases()->FindObject( "crmsiroc")) { delete gROOT->GetListOfCanvases()->FindObject("crmsiroc"); fExecPadIrocRms = 0; }
if(gROOT->GetListOfCanvases()->FindObject( "crowmax")) { delete gROOT->GetListOfCanvases()->FindObject("crowmax"); }
if(gROOT->GetListOfCanvases()->FindObject( "crowmax")) { delete gROOT->GetListOfCanvases()->FindObject("crowmax"); }
-
+
}
{
// Executable for global Histogram
// Will be called from /TPC/AliTPCMonitorExec.C(3)
- // Call ProcessEvent for same event and sector pointed at
+ // Call ProcessEvent for same event and sector pointed at
Int_t side = 0;
Int_t sector = 0;
if(!select->InheritsFrom("TH2")) {gPad->SetUniqueID(0); return -1; }
if( strcmp(select->GetName(),"hglobal" )==0 || ( strcmp(select->GetName(),"SIDE A" )==0) ) side = 0;
else if( strcmp(select->GetName(),"hglobal2")==0 || ( strcmp(select->GetName(),"SIDE C" )==0) ) side = 1;
-
+
// get position
Int_t px = gPad->GetEventX();
Int_t py = gPad->GetEventY();
Float_t upx = gPad->AbsPixeltoX(px);
Float_t y = gPad->PadtoY(upy);
Float_t x = gPad->PadtoX(upx);
-
+
Int_t testy = ((TH2*)select)->GetYaxis()->FindBin(y);
Int_t testx = ((TH2*)select)->GetXaxis()->FindBin(x);
if(((TH2*)select)->GetCellContent(testx,testy)==0) return -1 ;
-
+
Float_t alpha = 0.0;
if(x>0.0 && y > 0.0) alpha = TMath::Abs(TMath::ATan(TMath::Abs(x/y)));
if(x>0.0 && y < 0.0) alpha = TMath::Abs(TMath::ATan(TMath::Abs(y/x)));
if(fVerb) cout << "AliTPCMonitor::ExecProcess() next side " << side << " next sector " << sector << endl;
return (Int_t)ProcessEvent();
-
+
}
//__________________________________________________________________
Int_t AliTPCMonitor::GetRCUPatch(Int_t runid, Int_t eqid) const
{
- // Return RCU patch index for given equipment id eqid
+ // Return RCU patch index for given equipment id eqid
Int_t patch = 0;
//if(runid>=704)
if ( eqid<768 || eqid>983 ) return 0; //no TPC eqid
if(runid>=0)
- {
- if(eqid>=1000) return 0;
- patch = fMapEqidsRcu[eqid] ;
- }
+ {
+ if(eqid>=1000) return 0;
+ patch = fMapEqidsRcu[eqid] ;
+ }
else
- {
- if(eqid==408) {patch = 13*6+4 +0; }
- if(eqid==409) {patch = 13*6+5 +0; }
- if(eqid==509) {patch = 13*6+0 +0; }
- if(eqid==512) {patch = 13*6+3 +0; }
- if(eqid==513) {patch = 13*6+1 +0; }
- if(eqid==517) {patch = 13*6+2 +0; }
-
- if(eqid==404) {patch = 4*6+5 +0; }
- if(eqid==504) {patch = 4*6+4 +0; }
- if(eqid==407) {patch = 4*6+3 +0; }
- if(eqid==503) {patch = 4*6+2 +0; }
- if(eqid==508) {patch = 4*6+1 +0; }
- if(eqid==506) {patch = 4*6+0 +0; }
- }
+ {
+ if(eqid==408) {patch = 13*6+4 +0; }
+ if(eqid==409) {patch = 13*6+5 +0; }
+ if(eqid==509) {patch = 13*6+0 +0; }
+ if(eqid==512) {patch = 13*6+3 +0; }
+ if(eqid==513) {patch = 13*6+1 +0; }
+ if(eqid==517) {patch = 13*6+2 +0; }
+
+ if(eqid==404) {patch = 4*6+5 +0; }
+ if(eqid==504) {patch = 4*6+4 +0; }
+ if(eqid==407) {patch = 4*6+3 +0; }
+ if(eqid==503) {patch = 4*6+2 +0; }
+ if(eqid==508) {patch = 4*6+1 +0; }
+ if(eqid==506) {patch = 4*6+0 +0; }
+ }
return patch;
}
//__________________________________________________________________
-void AliTPCMonitor::DumpHeader(AliRawReaderRoot * readerROOT) const
+void AliTPCMonitor::DumpHeader(AliRawReader * reader) const
{
- // Dump Event header for ROOT format
-
- cout << "EventHeader : fReaderROOT->GetEquipmentSize() :" << readerROOT->GetEquipmentSize() << endl;
- cout << "EventHeader : fReaderROOT->GetType() :" << readerROOT->GetType() << endl;
- cout << "EventHeader : fReaderROOT->GetRunNumber() :" << readerROOT->GetRunNumber() << endl;
- cout << "EventHeader : fReaderROOT->GetEventId() :" << *(readerROOT->GetEventId()) << endl;
- cout << "EventHeader : fReaderROOT->GetLDCId() :" << readerROOT->GetLDCId() << endl;
- cout << "EventHeader : fReaderROOT->GetGDCId() :" << readerROOT->GetGDCId() << endl;
+ // Dump Event header for format
+
+ cout << "EventHeader : fReader->GetEquipmentSize() :" << reader->GetEquipmentSize() << endl;
+ cout << "EventHeader : fReader->GetType() :" << reader->GetType() << endl;
+ cout << "EventHeader : fReader->GetRunNumber() :" << reader->GetRunNumber() << endl;
+ cout << "EventHeader : fReader->GetEventId() :" << *(reader->GetEventId()) << endl;
+ cout << "EventHeader : fReader->GetLDCId() :" << reader->GetLDCId() << endl;
+ cout << "EventHeader : fReader->GetGDCId() :" << reader->GetGDCId() << endl;
}
-//__________________________________________________________________
-void AliTPCMonitor::DumpHeader(AliTPCMonitorDateFormat* dateform) const
-{
- // Dump Event header for DATE format
-
- cout << "EquipmentHeader : dateform->GetEquipmentSize() :" << dateform->GetEquipmentSize() << endl;
- cout << "EquipmentHeader : dateform->GetEquipmentType() :" << dateform->GetEquipmentType() << endl;
- cout << "EquipmentHeader : dateform->GetEquipmentID() :" << dateform->GetEquipmentID() << endl;
- cout << "EquipmentHeader : dateform->GetEquipmentTypeAttribute() :" << dateform->GetEquipmentTypeAttribute() << endl;
- cout << "EquipmentHeader : dateform->GetEquipmentBasicSize() :" << dateform->GetEquipmentBasicSize() << endl;
- cout << "EquipmentHeader : dateform->GetEquipmentHeaderSize() :" << dateform->GetEquipmentHeaderSize() << endl;
- cout << "EquipmentHeader : dateform->IsLastSubEventHeader() :" << dateform->IsLastSubEventHeader() << endl;
-}
//__________________________________________________________________
Double_t AliTPCMonitor::Gamma4(Double_t* x, Double_t* par) {
// Gamma4 function used to fit signals
- // Defined in sections: diverging branch set to 0
+ // Defined in sections: diverging branch set to 0
- Double_t val = 0.0;
+ Double_t val = 0.0;
if(x[0] > par[1])
val = par[0]*exp(4.0)* pow((x[0]-par[1])/par[2],4)*exp(-4.0*(x[0]-par[1])/par[2])+ par[3];
- else
+ else
val = 0;
return val;
}
//__________________________________________________________________
TCanvas* AliTPCMonitor::CreateCanvas(const Char_t* name)
{
- // Create Canvases
+ // Create Canvases
TCanvas* canv =0;
if( strcmp(name,"coroc" )==0) { canv = new TCanvas("coroc" ,"coroc" , -1+xoffset,(Int_t)(yspace+0.5*ysize) ,(Int_t)(1.5*xsize),(Int_t)(1.5*ysize)); return canv; }
else if(strcmp(name,"ciroc" )==0) { canv = new TCanvas("ciroc" ,"ciroc" , -1+xoffset, 0 ,(Int_t)(1.5*xsize),(Int_t)(1.5*ysize)); return canv; }
// ROC 2dim rms distribution
- else if(strcmp(name,"crmsoroc" )==0) { canv = new TCanvas("crmsoroc" ,"crmsoroc" , -1+xoffset,(Int_t)(yspace+0.5*ysize) ,(Int_t)(1.5*xsize),(Int_t)(1.5*ysize)); return canv; }
+ else if(strcmp(name,"crmsoroc" )==0) { canv = new TCanvas("crmsoroc" ,"crmsoroc" , -1+xoffset,(Int_t)(yspace+0.5*ysize) ,(Int_t)(1.5*xsize),(Int_t)(1.5*ysize)); return canv; }
else if(strcmp(name,"crmsiroc" )==0) { canv = new TCanvas("crmsiroc" ,"crmsiroc" , -1+xoffset, 0 ,(Int_t)(1.5*xsize),(Int_t)(1.5*ysize)); return canv; }
// Global ADC max Histos
else if(strcmp(name,"SIDE C all" )==0) { canv = new TCanvas("SIDE C all" ,"SIDE C all" , (Int_t)(3*xspace+ xoffset),(Int_t)(yspace+0.5*ysize) ,(Int_t)(1.5*xsize),(Int_t)(1.5*ysize)); return canv; }
// 1 dim max sum basekine distribution
else if(strcmp(name,"cmax" )==0) { canv = new TCanvas("cmax" ,"cmax" , -1+xoffset, 3*yspace ,(Int_t)(1.0*xsize),(Int_t)(1.0*ysize)); return canv; }
else if(strcmp(name,"csum" )==0) { canv = new TCanvas("csum" ,"csum" , xspace+xoffset, 3*yspace ,(Int_t)(1.0*xsize),(Int_t)(1.0*ysize)); return canv; }
- else if(strcmp(name,"cbasemean" )==0) { canv = new TCanvas("cbasemean" ,"cbasemean" , 2*xspace+xoffset, 3*yspace ,(Int_t)(1.0*xsize),(Int_t)(1.0*ysize)); return canv; }
+ else if(strcmp(name,"cbasemean" )==0) { canv = new TCanvas("cbasemean" ,"cbasemean" , 2*xspace+xoffset, 3*yspace ,(Int_t)(1.0*xsize),(Int_t)(1.0*ysize)); return canv; }
else if(strcmp(name,"cbaserms" )==0) { canv = new TCanvas("cbaserms" ,"cbaserms" , 3*xspace+xoffset, 3*yspace ,(Int_t)(1.0*xsize),(Int_t)(1.0*ysize)); return canv; }
// Projections of single channel
else if(strcmp(name,"coroc_ch" )==0) { canv = new TCanvas("coroc_ch" ,"coroc_ch" , (Int_t)(1.5*xspace+xoffset),(Int_t)(yspace+0.5*ysize),(Int_t)(1.5*xsize),(Int_t)(1.5*ysize)); return canv; }
void AliTPCMonitor::WriteHistos()
{
// Writes all available histograms to a file in current directory
- // File name will be specified by : sector, side, runid and eventnumber
+ // File name will be specified by : sector, side, runid and eventnumber
if(GetEventProcessed())
+ {
+ AliInfo("Write histos to file");
+ Char_t name[256];
+ sprintf(name,"SIDE_%i_SECTOR_%02i_RUN_%05i_EventID_%06i.root",(GetLastSector()/18),(GetLastSector()%18),fRunId,fEventNumber);
+ TFile* f = new TFile(name,"recreate");
+ for(Int_t i =0; i<fHistList->GetEntries(); i++)
{
- AliInfo("Write histos to file");
- Char_t name[256];
- sprintf(name,"SIDE_%i_SECTOR_%02i_RUN_%05i_EventID_%06i.root",(GetLastSector()/18),(GetLastSector()%18),fRunId,fEventNumber);
- TFile* f = new TFile(name,"recreate");
- for(Int_t i =0; i<fHistList->GetEntries(); i++)
- {
- if(((TH1*)fHistList->At(i))!=0)
- {
- ((TH1*)fHistList->At(i))->Write();
- }
- }
- f->ls();
- f->Close();
+ if(((TH1*)fHistList->At(i))!=0)
+ {
+ ((TH1*)fHistList->At(i))->Write();
+ }
}
+ f->ls();
+ f->Close();
+ }
else
- {
- AliError("No Event Processed : Chose Format , File and push 'Next Event' ");
- }
+ {
+ AliError("No Event Processed : Chose Format , File and push 'Next Event' ");
+ }
}
TH1* hist = 0;
if((TH1*)fHistList->FindObject(histname))
- {
- hist = (TH1*)fHistList->FindObject(histname);
- }
- else
- {
- cout << " AliTPCMonitor::GetHisto :: Can not find histo with name " << histname << endl;
- }
+ {
+ hist = (TH1*)fHistList->FindObject(histname);
+ }
+ else
+ {
+ cout << " AliTPCMonitor::GetHisto :: Can not find histo with name " << histname << endl;
+ }
return hist ;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff