T0/AliT0DataDCS.cxx

   1 #include "AliT0DataDCS.h"
   2
   3 #include "AliCDBMetaData.h"
   4 #include "AliDCSValue.h"
   5 #include "AliLog.h"
   6
   7 #include <TTimeStamp.h>
   8 #include <TObjString.h>
   9 #include <TH2F.h>
  10 #include <TProfile.h>
  11 #include <TGraph.h>
  12 #include <TDatime.h>
  13 #include <TStyle.h>
  14 #include <TCanvas.h>
  15
  16 ClassImp(AliT0DataDCS)
  17
  18 //---------------------------------------------------------------
  19 AliT0DataDCS::AliT0DataDCS():
  20         TObject(),
  21         fRun(0),
  22         fStartTime(0),
  23         fEndTime(0),
  24         fGraphs("TGraph",kNGraphs),
  25         fIsProcessed(kFALSE)
  26 {
  27         for(int i=0;i<kNHistos;i++) fHv[i]=0x0;
  28         fFunc = 0;
  29 }
  30
  31 //---------------------------------------------------------------
  32 AliT0DataDCS::AliT0DataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime):
  33         TObject(),
  34         fRun(nRun),
  35         fStartTime(startTime),
  36         fEndTime(endTime),
  37         fGraphs("TGraph",kNGraphs),
  38         fIsProcessed(kFALSE)
  39 {
  40         AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", nRun,
  41         TTimeStamp(startTime).AsString(),
  42         TTimeStamp(endTime).AsString()));
  43
  44         fFunc = 0;
  45         Init();
  46
  47 }
  48
  49 //---------------------------------------------------------------
  50 AliT0DataDCS::~AliT0DataDCS() {
  51
  52         for(int i=0;i<kNHistos;i++) {delete fHv[i]; fHv[i]=0;}
  53         fGraphs.Clear("C");
  54         fFunc=0;
  55 }
  56
  57 //---------------------------------------------------------------
  58 void AliT0DataDCS::ProcessData(TMap& aliasMap){
  59
  60         if(!(fHv[0])) Init();
  61
  62         TObjArray *aliasArr;
  63         AliDCSValue* aValue;
  64         for(int j=0; j<kNAliases; j++){
  65                 aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
  66                 if(!aliasArr){
  67                         AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
  68                         continue;
  69                 }
  70                 Introduce(j, aliasArr);
  71
  72                 if(aliasArr->GetEntries()<2){
  73                         AliError(Form("Alias %s has just %d entries!",
  74                                         fAliasNames[j].Data(),aliasArr->GetEntries()));
  75                         continue;
  76                 }
  77
  78                 TIter iterarray(aliasArr);
  79
  80                 Double_t *time = new Double_t[aliasArr->GetEntries()];
  81                 Double_t *val = new Double_t[aliasArr->GetEntries()];
  82
  83                 UInt_t ne=0;
  84                 while ((aValue = (AliDCSValue*) iterarray.Next())) {
  85
  86                 val[ne] = aValue->GetFloat();
  87                 time[ne] = (Double_t) (aValue->GetTimeStamp());
  88                 // fill histos (alias 0-2)
  89                 if(j < 3) fHv[j]->Fill(val[ne]);
  90                 ne++;
  91                 }
  92                 // fill graphs (alias 3-5)
  93                 if(j >= 3) CreateGraph(j, aliasArr->GetEntries(), time, val);
  94                 delete[] val;
  95                 delete[] time;
  96         }
  97
  98         // calculate mean and rms of the first two histos
  99         for(int i=0;i<kNHistos;i++){
 100                 fMean[i] = fHv[i]->GetMean();
 101                 fWidth[i] = fHv[i]->GetRMS();
 102         }
 103
 104         // pol1 fit of the first graph
 105         if(fGraphs.GetEntries() > 0){
 106                 ((TGraph*) fGraphs.UncheckedAt(0))->Fit("pol1");
 107                 fFunc = ((TGraph*) fGraphs.UncheckedAt(0))->GetFunction("pol1");
 108         }
 109
 110         fIsProcessed=kTRUE;
 111
 112
 113 }
 114
 115 //---------------------------------------------------------------
 116 void AliT0DataDCS::Init(){
 117
 118         TH1::AddDirectory(kFALSE);
 119
 120         fGraphs.SetOwner(1);
 121
 122         for(int i=0;i<kNAliases;i++){
 123                 fAliasNames[i] = "DCSAlias";
 124                 fAliasNames[i] += i;
 125         }
 126
 127         for(int i=0;i<kNHistos;i++){
 128                 fHv[i] = new TH1F(fAliasNames[i].Data(),fAliasNames[i].Data(), 20, kHvMin, kHvMax);
 129                 fHv[i]->GetXaxis()->SetTitle("Hv");
 130         }
 131 }
 132
 133 //---------------------------------------------------------------
 134 void AliT0DataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr){
 135
 136         int entries=aliasArr->GetEntries();
 137         AliInfo(Form("************ Alias: %s **********",fAliasNames[numAlias].Data()));
 138         AliInfo(Form("          %d DP values collected",entries));
 139
 140 }
 141
 142 //---------------------------------------------------------------
 143 void AliT0DataDCS::CreateGraph(int i, int dim, const Double_t *x, const Double_t *y)
 144 {
 145
 146         TGraph *gr = new(fGraphs[fGraphs.GetEntriesFast()]) TGraph(dim, x, y);
 147
 148         gr->GetXaxis()->SetTimeDisplay(1);
 149         gr->SetTitle(fAliasNames[i].Data());
 150
 151         AliInfo(Form("Array entries: %d",fGraphs.GetEntriesFast()));
 152
 153
 154 }
 155
 156 //---------------------------------------------------------------
 157 void AliT0DataDCS::Draw(const Option_t* /*option*/)
 158 {
 159 // Draw all histos and graphs
 160
 161   if(!fIsProcessed) return;
 162
 163   TCanvas *ch;
 164   TString canvasHistoName="Histos";
 165   ch=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);
 166   ch->Divide(2,2);
 167   ch->cd(1);
 168   fHv[0]->Draw();
 169   ch->cd(2);
 170   fHv[1]->Draw();
 171   ch->cd(3);
 172   fHv[2]->Draw();
 173
 174
 175   if(fGraphs.GetEntries() == 0) return;
 176
 177   TCanvas *cg;
 178   TString canvasGraphName="Graphs";
 179   cg=new TCanvas(canvasGraphName,canvasGraphName,40,40,600,600);
 180   cg->Divide(2,2);
 181   cg->cd(1);
 182   ((TGraph*) fGraphs.UncheckedAt(0))->Draw("alp");
 183
 184   cg->cd(2);
 185   ((TGraph*) fGraphs.UncheckedAt(1))->Draw("alp");
 186   cg->cd(3);
 187   ((TGraph*) fGraphs.UncheckedAt(2))->Draw("alp");
 188
 189   if(fFunc){
 190         cg->cd(4);
 191         fFunc->Draw("l");
 192   }
 193
 194 }
 195