VZERO/AliVZERODataDCS.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16
  17 #include "AliVZERODataDCS.h"
  18
  19 #include "AliDCSValue.h"
  20 #include "AliLog.h"
  21
  22 #include "TGraph.h"
  23 #include "TAxis.h"
  24 #include "TCanvas.h"
  25 #include "TTimeStamp.h"
  26 #include "TMap.h"
  27
  28 class TH2;
  29 class AliCDBMetaData;
  30 class TDatime;
  31
  32 // AliVZERODataDCS class
  33 // main aim to introduce the aliases for the VZERO DCS
  34 // data points to be then
  35 // stored in the OCDB, and to process them.
  36 // ProcessData() method called by VZEROPreprocessor
  37
  38 ClassImp(AliVZERODataDCS)
  39
  40 //_____________________________________________________________________________
  41 AliVZERODataDCS::AliVZERODataDCS():
  42         TObject(),
  43         fRun(0),
  44         fStartTime(0),
  45         fEndTime(0),
  46         fGraphs("TGraph",kNGraphs),
  47         fIsProcessed(kFALSE)
  48 {
  49   // Default constructor
  50 }
  51
  52 //_____________________________________________________________________________
  53 AliVZERODataDCS::AliVZERODataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime):
  54         TObject(),
  55         fRun(nRun),
  56         fStartTime(startTime),
  57         fEndTime(endTime),
  58         fGraphs("TGraph",kNGraphs),
  59         fIsProcessed(kFALSE)
  60 {
  61
  62   // constructor with arguments
  63
  64         AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", nRun,
  65         TTimeStamp(startTime).AsString(),
  66         TTimeStamp(endTime).AsString()));
  67
  68         Init();
  69
  70 }
  71
  72 //_____________________________________________________________________________
  73 AliVZERODataDCS::~AliVZERODataDCS() {
  74
  75   // destructor
  76   fGraphs.Clear("C");
  77
  78 }
  79
  80 //_____________________________________________________________________________
  81 void AliVZERODataDCS::ProcessData(TMap& aliasMap){
  82
  83   // method to process the data
  84
  85   if(!(fAliasNames[0])) Init();
  86
  87   TObjArray *aliasArr;
  88   AliDCSValue* aValue;
  89
  90   // starting loop on aliases
  91   for(int iAlias=0; iAlias<kNAliases; iAlias++){
  92
  93     aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[iAlias].Data());
  94     if(!aliasArr){
  95       AliError(Form("Alias %s not found!", fAliasNames[iAlias].Data()));
  96       return;
  97     }
  98
  99     //Introduce(iAlias, aliasArr);
 100
 101     if(aliasArr->GetEntries()<2){
 102       AliError(Form("Alias %s has just %d entries!",
 103                     fAliasNames[iAlias].Data(),aliasArr->GetEntries()));
 104       continue;
 105     }
 106
 107     TIter iterarray(aliasArr);
 108
 109     Int_t nentries = aliasArr->GetEntries();
 110
 111     Double_t *Times = new Double_t[nentries];
 112     Double_t *Values = new Double_t[nentries];
 113
 114     UInt_t iValue=0;
 115     while((aValue = (AliDCSValue*) iterarray.Next())) {
 116                 Values[iValue] = aValue->GetFloat();
 117                 Times[iValue] = (Double_t) (aValue->GetTimeStamp());
 118                 fHv[iAlias]->Fill(Values[iValue]);
 119                 printf("%s %f\n",fAliasNames[iAlias].Data(),Values[iValue]);
 120                 iValue++;
 121     }
 122     CreateGraph(iAlias, aliasArr->GetEntries(), Times, Values); // fill graphs
 123
 124     delete[] Values;
 125     delete[] Times;
 126   }
 127
 128         // calculate mean and rms of the first two histos
 129         // and convert index to aliroot channel
 130         for(int i=0;i<kNAliases;i++){
 131                 Int_t iChannel     = GetOfflineChannel(i);
 132                 fMeanHV[iChannel]  = fHv[i]->GetMean();
 133                 fWidthHV[iChannel] = fHv[i]->GetRMS();
 134         }
 135
 136   fIsProcessed=kTRUE;
 137 }
 138
 139 //_____________________________________________________________________________
 140 void AliVZERODataDCS::Init(){
 141
 142   // initialization of aliases and DCS data
 143
 144   TString sindex;
 145   int iAlias = 0;
 146
 147   for(int iSide = 0; iSide<2 ; iSide++){
 148         for(int iRing = 0; iRing<4 ; iRing++){
 149                 for(int iSector = 0; iSector<8 ; iSector++){
 150                         if(iSide == 0) fAliasNames[iAlias] = "V00/HV/V0A/SECTOR";
 151                         else fAliasNames[iAlias] = "V00/HV/V0C/SECTOR";
 152                         sindex.Form("%d/RING%d",iSector,iRing);
 153                         fAliasNames[iAlias] += sindex;
 154
 155                         fHv[iAlias] = new TH1F(fAliasNames[iAlias].Data(),fAliasNames[iAlias].Data(), 2000, kHvMin, kHvMax);
 156                         fHv[iAlias]->GetXaxis()->SetTitle("Hv");
 157                         iAlias++;
 158                 }
 159         }
 160   }
 161   if(iAlias!=kNAliases)
 162               AliError(Form("Number of DCS Aliases defined not correct"));
 163
 164 }
 165
 166 //_____________________________________________________________________________
 167 void AliVZERODataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)const
 168 {
 169
 170   // method to introduce new aliases
 171
 172   int entries=aliasArr->GetEntries();
 173   AliInfo(Form("************ Alias: %s **********",fAliasNames[numAlias].Data()));
 174   AliInfo(Form("        %d DP values collected",entries));
 175
 176 }
 177
 178 //_____________________________________________________________________________
 179 void AliVZERODataDCS::CreateGraph(int i, int dim, const Double_t *x, const Double_t *y)
 180 {
 181
 182    // Create graphics
 183
 184    TGraph *gr = new(fGraphs[fGraphs.GetEntriesFast()]) TGraph(dim, x, y);
 185
 186    gr->GetXaxis()->SetTimeDisplay(1);
 187    gr->SetTitle(fAliasNames[i].Data());
 188
 189    AliInfo(Form("Array entries: %d",fGraphs.GetEntriesFast()));
 190
 191 }
 192
 193
 194 //_____________________________________________________________________________
 195 void AliVZERODataDCS::Draw(const Option_t* /*option*/)
 196 {
 197 // Draw all histos and graphs
 198
 199   if(!fIsProcessed) return;
 200
 201   if(fGraphs.GetEntries()==0)  return;
 202
 203   TString CanvasName;
 204   TCanvas *cHV[8];
 205
 206   for(int iSide = 0 ;iSide<2;iSide++){
 207         for(int iRing=0;iRing<4;iRing++){
 208                 if(iSide == 0)  CanvasName = "V0A_Ring";
 209                 else  CanvasName = "V0C_Ring";
 210                 CanvasName += iRing;
 211                 int iCanvas = iSide*4 + iRing;
 212                 cHV[iCanvas] = new TCanvas(CanvasName,CanvasName);
 213                 cHV[iCanvas]->Divide(3,3);
 214                 for(int iSector=0;iSector<8;iSector++){
 215                         cHV[iCanvas]->cd(iSector+1);
 216                         int iChannel = iSide*32 + iRing*8 + iSector;
 217                         ((TGraph*) fGraphs.UncheckedAt(iChannel))->SetMarkerStyle(20);
 218                         ((TGraph*) fGraphs.UncheckedAt(iChannel))->Draw("ALP");
 219
 220                 }
 221
 222         }
 223   }
 224
 225 }
 226