ACORDE/AliACORDEDataDCS.cxx

   1 #include "AliACORDEDataDCS.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(AliACORDEDataDCS)
  17
  18 //---------------------------------------------------------------
  19 AliACORDEDataDCS::AliACORDEDataDCS():
  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 AliACORDEDataDCS::AliACORDEDataDCS(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 AliACORDEDataDCS::~AliACORDEDataDCS() {
  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 AliACORDEDataDCS::ProcessData(TMap& aliasMap){
  59
  60         if(!(fHv[0])) Init();
  61
  62         TObjArray *aliasArr;
  63         AliDCSValue* aValue;
  64
  65         for(int j=0; j<kNAliases; j++)
  66         {
  67                 aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
  68                 if(!aliasArr)
  69                 {
  70                         AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
  71                         continue;
  72                 }
  73                 Introduce(j, aliasArr);
  74
  75                 if(aliasArr->GetEntries()<2)
  76                 {
  77                         AliError(Form("Alias %s has just %d entries!",
  78                                         fAliasNames[j].Data(),aliasArr->GetEntries()));
  79                         continue;
  80                 }
  81
  82                 TIter iterarray(aliasArr);
  83
  84                 Double_t *time = new Double_t[aliasArr->GetEntries()];
  85                 Double_t *val  = new Double_t[aliasArr->GetEntries()];
  86
  87                 UInt_t ne=0;
  88
  89                 while ((aValue = (AliDCSValue*) iterarray.Next()))
  90                 {
  91                 val[ne] = aValue->GetFloat();
  92                 time[ne] = (Double_t) (aValue->GetTimeStamp());
  93                 fHv[j]->Fill(val[ne]);
  94                 ne++;
  95                 }
  96
  97
  98
  99                 CreateGraph(j, aliasArr->GetEntries(), time, val);
 100                 delete[] val;
 101                 delete[] time;
 102         }
 103
 104
 105         // calculate mean and rms of the first two histos
 106         for(int i=0;i<kNHistos;i++)
 107         {
 108                 fMean[i] = fHv[i]->GetMean();
 109                 fWidth[i] = fHv[i]->GetRMS();
 110         }
 111
 112
 113         fIsProcessed=kTRUE;
 114
 115
 116 }
 117
 118 //---------------------------------------------------------------
 119 void AliACORDEDataDCS::Init(){
 120
 121         TH1::AddDirectory(kFALSE);
 122
 123         fGraphs.SetOwner(1);
 124
 125         TString aliasName;
 126
 127         for(int i=0;i<kNAliases;i++){
 128
 129                 aliasName.Form("ACO_HV_MODULE%02d_VMON",i);
 130                 fAliasNames[i] = aliasName;
 131         }
 132
 133         for(int i=0;i<kNHistos;i++)
 134         {
 135                 fHv[i] = new TH1F(fAliasNames[i].Data(),fAliasNames[i].Data(), 20, kHvMin, kHvMax);
 136                 fHv[i]->GetXaxis()->SetTitle("Hv");
 137         }
 138 }
 139
 140 //---------------------------------------------------------------
 141 void AliACORDEDataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)
 142 {
 143
 144         int entries=aliasArr->GetEntries();
 145         AliInfo(Form("************ Alias: %s **********",fAliasNames[numAlias].Data()));
 146         AliInfo(Form("          %d DP values collected",entries));
 147
 148 }
 149
 150 //---------------------------------------------------------------
 151 void AliACORDEDataDCS::CreateGraph(int i, int dim, const Double_t *x, const Double_t *y)
 152 {
 153
 154         TGraph *gr = new(fGraphs[fGraphs.GetEntriesFast()]) TGraph(dim, x, y);
 155
 156         gr->GetXaxis()->SetTimeDisplay(1);
 157         gr->SetTitle(fAliasNames[i].Data());
 158
 159         AliInfo(Form("Array entries: %d",fGraphs.GetEntriesFast()));
 160
 161
 162 }
 163
 164 //---------------------------------------------------------------
 165 void AliACORDEDataDCS::Draw(const Option_t* /*option*/)
 166 {
 167 // Draw all histos and graphs
 168
 169   if(!fIsProcessed) return;
 170
 171   TString canvasHistoName;
 172   TCanvas *ch[10];
 173
 174   for (int i=0;i<10;i++)
 175   {
 176   canvasHistoName.Form("ACO_HV_MODULE");
 177   ch[i]=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);
 178   ch[i]->Divide(2,3);
 179
 180     for(int j=0;j<6;j++)
 181     {
 182     ch[i]->cd(j+1);
 183     ((TGraph*) fGraphs.UncheckedAt(i*6+j))->SetMarkerStyle(20);
 184     ((TGraph*) fGraphs.UncheckedAt(i*6+j))->Draw("alp");
 185     }
 186
 187   }
 188
 189
 190 }
 191