AliCDBId's in the list of retrieved OCDB parameters in the ESD's user info

[u/mrichter/AliRoot.git] / HLT / PHOS / AliHLTPHOSChannelCounter.cxx

 /**************************************************************************
 * This file is property of and copyright by the ALICE HLT Project        * 
 * All rights reserved.                                                   *
 *                                                                        *
 * Primary Authors: Oystein Djuvsland                                     *
 *                                                                        *
 * 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   *
 * 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          * 
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

#include "AliHLTPHOSChannelCounter.h"
#include "AliHLTPHOSRcuCellEnergyDataStruct.h"
#include "AliHLTPHOSValidCellDataStruct.h"
#include "AliHLTPHOSConstants.h"
#include "TFile.h"
#include "TH2F.h"
#include "TH2D.h"

//ClassImp(AliHLTPHOSChannelCounter);

using namespace PhosHLTConst;

AliHLTPHOSChannelCounter::AliHLTPHOSChannelCounter() :
  AliHLTPHOSBase(),
  // fChannelArrayPtr(0),
  fHistHighGainPtr(0),
  fHistLowGainPtr(0),
  fHistHighRatioPtr(0),
  fHistLowRatioPtr(0)
{
  //comment
  //  fChannelArrayPtr = new UInt_t[N_XCOLUMNS_MOD][N_ZROWS_MOD][N_GAINS];
  
  fHistHighGainPtr = new TH2I("highchannelcount", "High gain channel count", 
			      64, 0, 63, 56, 0, 56); 
  fHistLowGainPtr = new TH2I("lowchannelcount", "Low gain channel count", 
			      64, 0, 63, 56, 0, 56); 
  fHistHighRatioPtr = new TH2F("highoutofsync", "High gain channel count divided by number of events", 
			      64, 0, 63, 56, 0, 56); 
  fHistLowRatioPtr = new TH2F("lowoutofsync", "Low gain channel count divided by number of events", 
			     64, 0, 63, 56, 0, 56); 
  
   for(UInt_t x = 0; x < N_XCOLUMNS_MOD; x++)
     {
       for(UInt_t z = 0; z < N_ZROWS_MOD; z++)
	 {
	   for(UInt_t gain = 0; gain < 2; gain++)
	     {
	       fChannelArrayPtr[x][z][gain] = 0;
	     }
	 }
     }
}

AliHLTPHOSChannelCounter::~AliHLTPHOSChannelCounter()
{
//comment
}

void
AliHLTPHOSChannelCounter::CountChannels(AliHLTPHOSRcuCellEnergyDataStruct* channelDataPtr)
{
  //comment
  Int_t tmp[N_XCOLUMNS_MOD][N_ZROWS_MOD][N_GAINS];
  for(UInt_t x = 0; x < N_XCOLUMNS_MOD; x++)
  {
    for(UInt_t z = 0; z < N_ZROWS_MOD; z++)
    {
      for(UInt_t gain = 0; gain < 2; gain++)
      {
	tmp[x][z][gain] = 0;
      }
    }
  }
  for(Int_t i = 0; i < channelDataPtr->fCnt; i++)
    {
      AliHLTPHOSValidCellDataStruct *validDataPtr = &(channelDataPtr->fValidData[i]);
      
      tmp[validDataPtr->fX + channelDataPtr->fRcuX*N_XCOLUMNS_RCU]
	  [validDataPtr->fZ + channelDataPtr->fRcuZ*N_ZROWS_RCU]
	  [validDataPtr->fGain]++;
    }
  for(UInt_t x = 0; x < N_XCOLUMNS_MOD; x++)
    {
       for(UInt_t z = 0; z < N_ZROWS_MOD; z++)
	{
	  for(UInt_t gain = 0; gain < N_ZROWS_MOD; gain++)
	  {
	    if(tmp[x][z][gain] > 1)
	      fChannelArrayPtr[x][z][gain] = fChannelArrayPtr[x][z][gain] + 1;
	  }
	}
    }
}

void
AliHLTPHOSChannelCounter::PrintOutOfSyncChannels(Int_t nEvents)
{
  //comment
  printf("After %d events:\n", nEvents);
  for(UInt_t x = 0; x < N_XCOLUMNS_MOD; x++)
  //for(UInt_t x = 0; x < 63; x++)
    {
      for(UInt_t z = 0; z < N_ZROWS_MOD; z++)
      //for(UInt_t z = 0; z < 55; z++)
	{
	  for(UInt_t gain = 0; gain < N_GAINS; gain++)
	  //for(UInt_t gain = 0; gain < 2; gain++)
	    {
	      printf("x = %d -- z = %d -- gain = %d:   %d\n", x, z, gain, (fChannelArrayPtr[x][z][gain]));
	    }
	}
    }
  printf("\n");
}

void
AliHLTPHOSChannelCounter::FillHistograms(Int_t nEvents)
{
  //comment
  printf("Filling histograms...");
  //for(Int_t x = 0; x < N_XCOLUMNS_MOD; x++)
  for(Int_t x = 0; x < 63; x++)
    {
      //  for(Int_t z = 0; z < N_ZROWS_MOD; z++)
      for(Int_t z = 0; z < 55; z++)
	{
	  for(Int_t gain = 0; gain < 2; gain++)
	    {
	      
	      if(gain == 0)  //TODO: this line was "if(gain = 0)" I presume that was a bug. Please check.
		{
		  fHistHighGainPtr->SetBinContent(x, z, fChannelArrayPtr[x][z][gain]);
		  fHistHighRatioPtr->SetBinContent(x, z, fChannelArrayPtr[x][z][gain]/(float)nEvents);
		  continue;
		}
	      
	      fHistLowGainPtr->SetBinContent(x, z, fChannelArrayPtr[x][z][gain]);
	      fHistLowRatioPtr->SetBinContent(x, z, fChannelArrayPtr[x][z][gain]/(float)nEvents);
	    }
	}
    }
  printf("Done!\n");
}

void 
AliHLTPHOSChannelCounter::WriteHistograms(const char* filename)
{
  //comment
  printf("Writing histograms to file...");

  TFile *outfile = new TFile(filename,"recreate");

  fHistHighGainPtr->Write();
  fHistLowGainPtr->Write();
  fHistHighRatioPtr->Write();
  fHistLowRatioPtr->Write();

  outfile->Close();

  delete outfile;
  outfile = 0;
  printf("Done!\n");
}