[u/mrichter/AliRoot.git] / TRD / AliTRDQADataMakerSim.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * 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.                  *
 **************************************************************************/

/* $Id$ */

////////////////////////////////////////////////////////////////////////////
//                                                                        //
//  Produces the data needed to calculate the quality assurance.          //
//  All data must be mergeable objects.                                   //
//                                                                        //
//  Author:                                                               //
//    Sylwester Radomski (radomski@physi.uni-heidelberg.de)               //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

// --- ROOT system ---
#include <TClonesArray.h>
#include <TFile.h> 
#include <TH1F.h> 
#include <TH2F.h>
#include <TH3F.h>
#include <TProfile.h>
#include <TF1.h>
#include <TCanvas.h>
#include <TTree.h>

// --- AliRoot header files ---
//#include "AliESDEvent.h"
#include "AliLog.h"
#include "AliTRDdigit.h"
#include "AliTRDhit.h"
//#include "AliTRDcluster.h"
#include "AliTRDQADataMakerSim.h"
#include "AliTRDdigitsManager.h"
#include "AliTRDgeometry.h"
#include "AliTRDarrayADC.h"
#include "AliTRDarraySignal.h"
//#include "AliTRDrawStream.h"

#include "AliQAChecker.h"

ClassImp(AliTRDQADataMakerSim)

//____________________________________________________________________________ 
  AliTRDQADataMakerSim::AliTRDQADataMakerSim() : 
  AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kTRD), "TRD Quality Assurance Data Maker"),
  fTmpHits(NULL) 
{
  //
  // Default constructor
}

//____________________________________________________________________________ 
AliTRDQADataMakerSim::AliTRDQADataMakerSim(const AliTRDQADataMakerSim& qadm) :
  AliQADataMakerSim(), 
  fTmpHits(NULL)
{
  //
  // Copy constructor 
  //

  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle()); 

}

//__________________________________________________________________
AliTRDQADataMakerSim& AliTRDQADataMakerSim::operator=(const AliTRDQADataMakerSim& qadm)
{
  //
  // Equal operator.
  //

  this->~AliTRDQADataMakerSim();
  new(this) AliTRDQADataMakerSim(qadm);
  return *this;

}

//____________________________________________________________________________ 
AliTRDQADataMakerSim::~AliTRDQADataMakerSim()
{
  if (fTmpHits) {
    fTmpHits->Clear() ; 
    delete fTmpHits ; 
  }
}

//____________________________________________________________________________ 
void AliTRDQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
{
  //
  // Detector specific actions at end of cycle
  //
  ResetEventTrigClasses();
  //AliDebug(AliQAv1::GetQADebugLevel(), Form("EndOfCycle", "Fitting RecPoints %d", task));

  // call the checker
  AliQAChecker::Instance()->Run(AliQAv1::kTRD, task, list) ;    


}

//____________________________________________________________________________ 
void AliTRDQADataMakerSim::InitHits()
{
  //
  // Create Hits histograms in Hits subdir
  //
  const Bool_t expert   = kTRUE ; 
  const Bool_t image    = kTRUE ; 
  
  const Int_t kNhist = 4;
  TH1F *hist[kNhist];

  hist[0] = new TH1F("qaTRD_hits_det", "TRD hits det;Detector Id of the hit;Counts", 540, -0.5, 539.5) ; 

  hist[1] = new TH1F("qaTRD_hist_Qdrift", "TRD hits Qdrift;Charge from tracks;Counts", 100, 0, 100);
  hist[2] = new TH1F("qaTRD_hist_Qamp", "TRD hits Qamp;Charge from TRD photon;Counts", 100, 0, 100);
  hist[3] = new TH1F("qaTRD_hist_Qphoton", "TRD hits Qphoton;Charge from TRD photon;Counts", 100, 0, 100);

  for(Int_t i=0; i<kNhist; i++) {
    //hist[i]->Sumw2();
    Add2HitsList(hist[i], i, !expert, image);
  }
  //
  ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
}

//____________________________________________________________________________ 
void AliTRDQADataMakerSim::InitDigits()
{
  //
  // Create Digits histograms in Digits subdir
  //
  const Bool_t expert   = kTRUE ; 
  const Bool_t image    = kTRUE ; 
  
  const Int_t kNhist = 3;
  TH1F *hist[kNhist];

  hist[0] = new TH1F("qaTRD_digits_det", "TRD digits det;Detector Id of the digit;Counts", 540, -0.5, 539.5);
  hist[1] = new TH1F("qaTRD_digits_time", "TRDdigits time;Time bin;Counts", 40, -0.5, 39.5);
  hist[2] = new TH1F("qaTRD_digits_amp", "TRD digits amp;Amplitude;Counts", 100, -5.5, 94.5);

  for(Int_t i=0; i<kNhist; i++) {
    hist[i]->Sumw2();
    Add2DigitsList(hist[i], i, !expert, image);
  }
  //
  ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
}

//____________________________________________________________________________ 
void AliTRDQADataMakerSim::InitSDigits()
{
  //
  // Create SDigits histograms in SDigits subdir
  //
  const Bool_t expert   = kTRUE ; 
  const Bool_t image    = kTRUE ; 
  
  const Int_t kNhist = 3;
  TH1F *hist[kNhist];

  hist[0] = new TH1F("qaTRD_sdigits_det", "TRD sdigits det;Detector Id of the digit;Counts", 540, -0.5, 539.5);
  hist[1] = new TH1F("qaTRD_sdigits_time", "TRD sdigits time;Time bin;Counts", 40, -0.5, 39.5);
  hist[2] = new TH1F("qaTRD_sdigits_amp", "TRD sdigits amp;Amplitude;Counts", 100, 0, 1e7);

  for(Int_t i=0; i<kNhist; i++) {
    hist[i]->Sumw2();
    Add2SDigitsList(hist[i], i, !expert, image);
  }
  //
  ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line
}

//____________________________________________________________________________
void AliTRDQADataMakerSim::MakeHits()
{
  //
  // Make QA data from Hits
  //
  
  TIter next(fHitsArray); 
  AliTRDhit * hit; 

  while ( (hit = dynamic_cast<AliTRDhit *>(next())) ) {
    FillHitsData(0,hit->GetDetector());
    Double_t q = TMath::Abs(hit->GetCharge());

    if (hit->FromDrift()) FillHitsData(1,q);
    if (hit->FromAmplification()) FillHitsData(2,q);
    if (hit->FromTRphoton()) FillHitsData(3,q);
  }
}

//____________________________________________________________________________
void AliTRDQADataMakerSim::MakeHits(TTree * hitTree)
{
  //
  // Make QA data from Hits
  //

  if (!CheckPointer(hitTree, "TRD hits tree")) return;

  TBranch *branch = hitTree->GetBranch("TRD");
  if (!CheckPointer(branch, "TRD hits branch")) return;
  
  Int_t nhits = (Int_t)(hitTree->GetTotBytes()/sizeof(AliTRDhit));
  if (fHitsArray)
    fHitsArray->Clear() ; 
  else
   fHitsArray = new TClonesArray("AliTRDhit", nhits+1000);
  
  //Int_t index = 0;
  Int_t nEntries = (Int_t)branch->GetEntries();
  for(Int_t i = 0; i < nEntries; i++) {
    branch->GetEntry(i);
    MakeHits();
    fHitsArray->Clear() ; 
  }
  //
  IncEvCountCycleHits();
  IncEvCountTotalHits();
  //
}

//____________________________________________________________________________
void AliTRDQADataMakerSim::MakeDigits()
{
  //
  // Makes data from Digits
  //

  if (!fDigitsArray)
    return ; 
  
  TIter next(fDigitsArray) ; 
  AliTRDdigit * digit ; 
  
  // Info("Make digits", "From the arrya");
  
  while ( (digit = dynamic_cast<AliTRDdigit *>(next())) ) {
    if (digit->GetAmp() < 1) continue;
    FillDigitsData(0,digit->GetDetector());
    FillDigitsData(1,digit->GetTime());
    FillDigitsData(2,digit->GetAmp());
  }

}

//____________________________________________________________________________
void AliTRDQADataMakerSim::MakeDigits(TTree * digits)
{
  //
  // Makes data from digits tree
  //
  // Info("Make digits", "From a tree");

  AliTRDdigitsManager *digitsManager = new AliTRDdigitsManager();
  digitsManager->CreateArrays();
  digitsManager->ReadDigits(digits);
  static TObjArray histDet,histTime,histSignal;
  //
  GetMatchingDigitsData(0,&histDet);
  GetMatchingDigitsData(1,&histTime);
  GetMatchingDigitsData(2,&histSignal);

  for (Int_t i = 0; i < AliTRDgeometry::kNdet; i++) 
    {
      AliTRDarrayADC *digitsIn = (AliTRDarrayADC *) digitsManager->GetDigits(i);      

      // This is to take care of switched off super modules
      if (digitsIn->GetNtime() == 0) continue;
      
      digitsIn->Expand();
      
      //AliTRDSignalIndex* indexes = digitsManager->GetIndexes(i);
      //if (indexes->IsAllocated() == kFALSE) digitsManager->BuildIndexes(i);
      
      Int_t nRows = digitsIn->GetNrow();
      Int_t nCols = digitsIn->GetNcol();
      Int_t nTbins = digitsIn->GetNtime();
      
      for(Int_t row = 0; row < nRows; row++) 
        for(Int_t col = 0; col < nCols; col++) 
          for(Int_t time = 0; time < nTbins; time++) 
            {   
              Float_t signal = digitsIn->GetData(row,col,time);
              if (signal < 1) continue;
              for (int ih=histDet.GetEntriesFast();ih--;) ((TH1*)histDet.UncheckedAt(ih))->Fill(i);
              for (int ih=histTime.GetEntriesFast();ih--;) ((TH1*)histTime.UncheckedAt(ih))->Fill(time);
              for (int ih=histSignal.GetEntriesFast();ih--;) ((TH1*)histSignal.UncheckedAt(ih))->Fill(signal);
            }
      
      //delete digitsIn;
    }
  delete digitsManager;
  //
  IncEvCountCycleDigits();
  IncEvCountTotalDigits();
  //
}

//____________________________________________________________________________
void AliTRDQADataMakerSim::MakeSDigits()
{
  //
  // Makes data from Digits
  //

  if (!fSDigitsArray)
    return ; 

  TIter next(fSDigitsArray) ; 
  AliTRDdigit * digit ; 
  while ( (digit = dynamic_cast<AliTRDdigit *>(next())) ) {
    FillDigitsData(0,digit->GetDetector());
    FillDigitsData(1,digit->GetTime());
    FillDigitsData(2,digit->GetAmp());
  }

}

//____________________________________________________________________________
void AliTRDQADataMakerSim::MakeSDigits(TTree * digits)
{
  //
  // Makes data from SDigits
  //
  // Check id histograms already created for this Event Specie
  if ( ! GetSDigitsData(0) )
    InitSDigits() ;
  
  AliTRDdigitsManager *digitsManager = new AliTRDdigitsManager();
  digitsManager->SetSDigits(kTRUE);
  digitsManager->CreateArrays();
  digitsManager->ReadDigits(digits);

  static TObjArray histDet,histTime,histSignal;
  //
  GetMatchingSDigitsData(0,&histDet);
  GetMatchingSDigitsData(1,&histTime);
  GetMatchingSDigitsData(2,&histSignal);

  for (Int_t i = 0; i < AliTRDgeometry::kNdet; i++) 
    {
      AliTRDarraySignal *digitsIn = (AliTRDarraySignal *) digitsManager->GetSDigits(i);
  
    // This is to take care of switched off super modules
      if (digitsIn->GetNtime() == 0) continue;
      
    digitsIn->Expand(); 
      
      //AliTRDSignalIndex* indexes = digitsManager->GetIndexes(i);
      //if (indexes->IsAllocated() == kFALSE) digitsManager->BuildIndexes(i);
      Int_t nRows = digitsIn->GetNrow();
      Int_t nCols = digitsIn->GetNcol();
      Int_t nTbins = digitsIn->GetNtime();

      for(Int_t row = 0; row < nRows; row++) 
	{
	  for(Int_t col = 0; col < nCols; col++) 
	    {
	      for(Int_t time = 0; time < nTbins; time++) 
		{	    
		  Float_t signal = digitsIn->GetData(row,col,time);
		  if (signal < 1) continue;
		  for (int ih=histDet.GetEntriesFast();ih--;) ((TH1*)histDet.UncheckedAt(ih))->Fill(i);
		  for (int ih=histTime.GetEntriesFast();ih--;) ((TH1*)histTime.UncheckedAt(ih))->Fill(time);
		  for (int ih=histSignal.GetEntriesFast();ih--;) ((TH1*)histSignal.UncheckedAt(ih))->Fill(signal);
		}
	    }
	}
      // delete digitsIn;
    }
  delete digitsManager;
  //
  IncEvCountCycleSDigits();
  IncEvCountTotalSDigits();  
  //
}

//____________________________________________________________________________ 
void AliTRDQADataMakerSim::StartOfDetectorCycle()
{
  //
  // Detector specific actions at start of cycle
  //

}

//__________________________________________________________________________
Int_t AliTRDQADataMakerSim::CheckPointer(TObject *obj, const char *name) 
{
  //
  // Checks initialization of pointers
  //

  if (!obj) AliWarning(Form("null pointer: %s", name));
  return !!obj;

}
Commit	Line	Data
04236e67	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	/* $Id$ */
	17
	18	////////////////////////////////////////////////////////////////////////////
	19	// //
	20	// Produces the data needed to calculate the quality assurance. //
	21	// All data must be mergeable objects. //
	22	// //
	23	// Author: //
	24	// Sylwester Radomski (radomski@physi.uni-heidelberg.de) //
	25	// //
	26	////////////////////////////////////////////////////////////////////////////
	27
	28	// --- ROOT system ---
	29	#include <TClonesArray.h>
	30	#include <TFile.h>
92664bc8	31	#include <TH1F.h>
	32	#include <TH2F.h>
	33	#include <TH3F.h>
04236e67	34	#include <TProfile.h>
	35	#include <TF1.h>
	36	#include <TCanvas.h>
	37	#include <TTree.h>
	38
	39	// --- AliRoot header files ---
	40	//#include "AliESDEvent.h"
	41	#include "AliLog.h"
	42	#include "AliTRDdigit.h"
	43	#include "AliTRDhit.h"
	44	//#include "AliTRDcluster.h"
	45	#include "AliTRDQADataMakerSim.h"
	46	#include "AliTRDdigitsManager.h"
	47	#include "AliTRDgeometry.h"
b65e5048	48	#include "AliTRDarrayADC.h"
b65e5048	49	#include "AliTRDarraySignal.h"
987ba9a3	50	//#include "AliTRDrawStream.h"
04236e67	51
	52	#include "AliQAChecker.h"
	53
	54	ClassImp(AliTRDQADataMakerSim)
	55
	56	//____________________________________________________________________________
	57	AliTRDQADataMakerSim::AliTRDQADataMakerSim() :
6252ceeb	58	AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kTRD), "TRD Quality Assurance Data Maker"),
6252ceeb	59	fTmpHits(NULL)
04236e67	60	{
	61	//
	62	// Default constructor
	63	}
	64
	65	//____________________________________________________________________________
	66	AliTRDQADataMakerSim::AliTRDQADataMakerSim(const AliTRDQADataMakerSim& qadm) :
6252ceeb	67	AliQADataMakerSim(),
6252ceeb	68	fTmpHits(NULL)
04236e67	69	{
	70	//
	71	// Copy constructor
	72	//
	73
	74	SetName((const char*)qadm.GetName()) ;
	75	SetTitle((const char*)qadm.GetTitle());
	76
	77	}
	78
	79	//__________________________________________________________________
	80	AliTRDQADataMakerSim& AliTRDQADataMakerSim::operator=(const AliTRDQADataMakerSim& qadm)
	81	{
	82	//
	83	// Equal operator.
	84	//
	85
	86	this->~AliTRDQADataMakerSim();
	87	new(this) AliTRDQADataMakerSim(qadm);
	88	return *this;
	89
	90	}
	91
6252ceeb	92	//____________________________________________________________________________
	93	AliTRDQADataMakerSim::~AliTRDQADataMakerSim()
	94	{
	95	if (fTmpHits) {
	96	fTmpHits->Clear() ;
	97	delete fTmpHits ;
	98	}
	99	}
	100
04236e67	101	//____________________________________________________________________________
4e25ac79	102	void AliTRDQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
04236e67	103	{
	104	//
	105	// Detector specific actions at end of cycle
	106	//
92664bc8	107	ResetEventTrigClasses();
5379c4a3	108	//AliDebug(AliQAv1::GetQADebugLevel(), Form("EndOfCycle", "Fitting RecPoints %d", task));
04236e67	109
04236e67	110	// call the checker
4e25ac79	111	AliQAChecker::Instance()->Run(AliQAv1::kTRD, task, list) ;
04236e67	112
	113
	114	}
	115
	116	//____________________________________________________________________________
	117	void AliTRDQADataMakerSim::InitHits()
	118	{
	119	//
	120	// Create Hits histograms in Hits subdir
	121	//
7d297381	122	const Bool_t expert = kTRUE ;
	123	const Bool_t image = kTRUE ;
	124
04236e67	125	const Int_t kNhist = 4;
92664bc8	126	TH1F *hist[kNhist];
04236e67	127
92664bc8	128	hist[0] = new TH1F("qaTRD_hits_det", "TRD hits det;Detector Id of the hit;Counts", 540, -0.5, 539.5) ;
04236e67	129
92664bc8	130	hist[1] = new TH1F("qaTRD_hist_Qdrift", "TRD hits Qdrift;Charge from tracks;Counts", 100, 0, 100);
	131	hist[2] = new TH1F("qaTRD_hist_Qamp", "TRD hits Qamp;Charge from TRD photon;Counts", 100, 0, 100);
	132	hist[3] = new TH1F("qaTRD_hist_Qphoton", "TRD hits Qphoton;Charge from TRD photon;Counts", 100, 0, 100);
04236e67	133
	134	for(Int_t i=0; i<kNhist; i++) {
	135	//hist[i]->Sumw2();
7d297381	136	Add2HitsList(hist[i], i, !expert, image);
04236e67	137	}
92664bc8	138	//
92664bc8	139	ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
04236e67	140	}
	141
	142	//____________________________________________________________________________
	143	void AliTRDQADataMakerSim::InitDigits()
	144	{
	145	//
	146	// Create Digits histograms in Digits subdir
	147	//
7d297381	148	const Bool_t expert = kTRUE ;
	149	const Bool_t image = kTRUE ;
	150
04236e67	151	const Int_t kNhist = 3;
92664bc8	152	TH1F *hist[kNhist];
04236e67	153
92664bc8	154	hist[0] = new TH1F("qaTRD_digits_det", "TRD digits det;Detector Id of the digit;Counts", 540, -0.5, 539.5);
	155	hist[1] = new TH1F("qaTRD_digits_time", "TRDdigits time;Time bin;Counts", 40, -0.5, 39.5);
	156	hist[2] = new TH1F("qaTRD_digits_amp", "TRD digits amp;Amplitude;Counts", 100, -5.5, 94.5);
04236e67	157
	158	for(Int_t i=0; i<kNhist; i++) {
	159	hist[i]->Sumw2();
7d297381	160	Add2DigitsList(hist[i], i, !expert, image);
04236e67	161	}
92664bc8	162	//
92664bc8	163	ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
04236e67	164	}
	165
	166	//____________________________________________________________________________
	167	void AliTRDQADataMakerSim::InitSDigits()
	168	{
	169	//
	170	// Create SDigits histograms in SDigits subdir
	171	//
7d297381	172	const Bool_t expert = kTRUE ;
	173	const Bool_t image = kTRUE ;
	174
04236e67	175	const Int_t kNhist = 3;
92664bc8	176	TH1F *hist[kNhist];
04236e67	177
92664bc8	178	hist[0] = new TH1F("qaTRD_sdigits_det", "TRD sdigits det;Detector Id of the digit;Counts", 540, -0.5, 539.5);
	179	hist[1] = new TH1F("qaTRD_sdigits_time", "TRD sdigits time;Time bin;Counts", 40, -0.5, 39.5);
	180	hist[2] = new TH1F("qaTRD_sdigits_amp", "TRD sdigits amp;Amplitude;Counts", 100, 0, 1e7);
04236e67	181
	182	for(Int_t i=0; i<kNhist; i++) {
	183	hist[i]->Sumw2();
7d297381	184	Add2SDigitsList(hist[i], i, !expert, image);
04236e67	185	}
92664bc8	186	//
92664bc8	187	ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line
04236e67	188	}
	189
	190	//____________________________________________________________________________
6252ceeb	191	void AliTRDQADataMakerSim::MakeHits()
04236e67	192	{
	193	//
	194	// Make QA data from Hits
	195	//
6252ceeb	196
6252ceeb	197	TIter next(fHitsArray);
04236e67	198	AliTRDhit * hit;
	199
	200	while ( (hit = dynamic_cast<AliTRDhit *>(next())) ) {
92664bc8	201	FillHitsData(0,hit->GetDetector());
04236e67	202	Double_t q = TMath::Abs(hit->GetCharge());
04236e67	203
92664bc8	204	if (hit->FromDrift()) FillHitsData(1,q);
	205	if (hit->FromAmplification()) FillHitsData(2,q);
	206	if (hit->FromTRphoton()) FillHitsData(3,q);
04236e67	207	}
04236e67	208	}
	209
	210	//____________________________________________________________________________
	211	void AliTRDQADataMakerSim::MakeHits(TTree * hitTree)
	212	{
	213	//
	214	// Make QA data from Hits
	215	//
	216
	217	if (!CheckPointer(hitTree, "TRD hits tree")) return;
	218
	219	TBranch *branch = hitTree->GetBranch("TRD");
	220	if (!CheckPointer(branch, "TRD hits branch")) return;
eca4fa66	221
04236e67	222	Int_t nhits = (Int_t)(hitTree->GetTotBytes()/sizeof(AliTRDhit));
6252ceeb	223	if (fHitsArray)
	224	fHitsArray->Clear() ;
	225	else
	226	fHitsArray = new TClonesArray("AliTRDhit", nhits+1000);
	227
bb5265c9	228	//Int_t index = 0;
04236e67	229	Int_t nEntries = (Int_t)branch->GetEntries();
	230	for(Int_t i = 0; i < nEntries; i++) {
	231	branch->GetEntry(i);
6252ceeb	232	MakeHits();
6252ceeb	233	fHitsArray->Clear() ;
04236e67	234	}
92664bc8	235	//
	236	IncEvCountCycleHits();
	237	IncEvCountTotalHits();
	238	//
04236e67	239	}
	240
	241	//____________________________________________________________________________
6252ceeb	242	void AliTRDQADataMakerSim::MakeDigits()
04236e67	243	{
	244	//
	245	// Makes data from Digits
	246	//
	247
6252ceeb	248	if (!fDigitsArray)
6252ceeb	249	return ;
eca4fa66	250
6252ceeb	251	TIter next(fDigitsArray) ;
04236e67	252	AliTRDdigit * digit ;
3590bda8	253
	254	// Info("Make digits", "From the arrya");
	255
04236e67	256	while ( (digit = dynamic_cast<AliTRDdigit *>(next())) ) {
3590bda8	257	if (digit->GetAmp() < 1) continue;
92664bc8	258	FillDigitsData(0,digit->GetDetector());
	259	FillDigitsData(1,digit->GetTime());
	260	FillDigitsData(2,digit->GetAmp());
04236e67	261	}
	262
	263	}
	264
	265	//____________________________________________________________________________
	266	void AliTRDQADataMakerSim::MakeDigits(TTree * digits)
	267	{
	268	//
	269	// Makes data from digits tree
	270	//
3590bda8	271	// Info("Make digits", "From a tree");
3590bda8	272
04236e67	273	AliTRDdigitsManager *digitsManager = new AliTRDdigitsManager();
	274	digitsManager->CreateArrays();
	275	digitsManager->ReadDigits(digits);
92664bc8	276	static TObjArray histDet,histTime,histSignal;
	277	//
	278	GetMatchingDigitsData(0,&histDet);
	279	GetMatchingDigitsData(1,&histTime);
	280	GetMatchingDigitsData(2,&histSignal);
3590bda8	281
b65e5048	282	for (Int_t i = 0; i < AliTRDgeometry::kNdet; i++)
	283	{
	284	AliTRDarrayADC digitsIn = (AliTRDarrayADC ) digitsManager->GetDigits(i);
	285
	286	// This is to take care of switched off super modules
	287	if (digitsIn->GetNtime() == 0) continue;
	288
	289	digitsIn->Expand();
	290
	291	//AliTRDSignalIndex* indexes = digitsManager->GetIndexes(i);
	292	//if (indexes->IsAllocated() == kFALSE) digitsManager->BuildIndexes(i);
	293
	294	Int_t nRows = digitsIn->GetNrow();
	295	Int_t nCols = digitsIn->GetNcol();
	296	Int_t nTbins = digitsIn->GetNtime();
	297
	298	for(Int_t row = 0; row < nRows; row++)
eca4fa66	299	for(Int_t col = 0; col < nCols; col++)
	300	for(Int_t time = 0; time < nTbins; time++)
	301	{
	302	Float_t signal = digitsIn->GetData(row,col,time);
	303	if (signal < 1) continue;
92664bc8	304	for (int ih=histDet.GetEntriesFast();ih--;) ((TH1*)histDet.UncheckedAt(ih))->Fill(i);
	305	for (int ih=histTime.GetEntriesFast();ih--;) ((TH1*)histTime.UncheckedAt(ih))->Fill(time);
	306	for (int ih=histSignal.GetEntriesFast();ih--;) ((TH1*)histSignal.UncheckedAt(ih))->Fill(signal);
eca4fa66	307	}
	308
	309	//delete digitsIn;
	310	}
04236e67	311	delete digitsManager;
92664bc8	312	//
	313	IncEvCountCycleDigits();
	314	IncEvCountTotalDigits();
	315	//
04236e67	316	}
	317
	318	//____________________________________________________________________________
6252ceeb	319	void AliTRDQADataMakerSim::MakeSDigits()
04236e67	320	{
	321	//
	322	// Makes data from Digits
	323	//
	324
6252ceeb	325	if (!fSDigitsArray)
	326	return ;
	327
	328	TIter next(fSDigitsArray) ;
04236e67	329	AliTRDdigit * digit ;
04236e67	330	while ( (digit = dynamic_cast<AliTRDdigit *>(next())) ) {
92664bc8	331	FillDigitsData(0,digit->GetDetector());
	332	FillDigitsData(1,digit->GetTime());
	333	FillDigitsData(2,digit->GetAmp());
04236e67	334	}
	335
	336	}
	337
	338	//____________________________________________________________________________
	339	void AliTRDQADataMakerSim::MakeSDigits(TTree * digits)
	340	{
	341	//
	342	// Makes data from SDigits
	343	//
eca4fa66	344	// Check id histograms already created for this Event Specie
	345	if ( ! GetSDigitsData(0) )
	346	InitSDigits() ;
	347
04236e67	348	AliTRDdigitsManager *digitsManager = new AliTRDdigitsManager();
c10bf383	349	digitsManager->SetSDigits(kTRUE);
04236e67	350	digitsManager->CreateArrays();
	351	digitsManager->ReadDigits(digits);
	352
92664bc8	353	static TObjArray histDet,histTime,histSignal;
	354	//
	355	GetMatchingSDigitsData(0,&histDet);
	356	GetMatchingSDigitsData(1,&histTime);
	357	GetMatchingSDigitsData(2,&histSignal);
3590bda8	358
b65e5048	359	for (Int_t i = 0; i < AliTRDgeometry::kNdet; i++)
	360	{
	361	AliTRDarraySignal digitsIn = (AliTRDarraySignal ) digitsManager->GetSDigits(i);
	362
04236e67	363	// This is to take care of switched off super modules
b65e5048	364	if (digitsIn->GetNtime() == 0) continue;
	365
	366	digitsIn->Expand();
	367
	368	//AliTRDSignalIndex* indexes = digitsManager->GetIndexes(i);
	369	//if (indexes->IsAllocated() == kFALSE) digitsManager->BuildIndexes(i);
	370	Int_t nRows = digitsIn->GetNrow();
	371	Int_t nCols = digitsIn->GetNcol();
	372	Int_t nTbins = digitsIn->GetNtime();
	373
	374	for(Int_t row = 0; row < nRows; row++)
	375	{
	376	for(Int_t col = 0; col < nCols; col++)
	377	{
	378	for(Int_t time = 0; time < nTbins; time++)
	379	{
	380	Float_t signal = digitsIn->GetData(row,col,time);
	381	if (signal < 1) continue;
92664bc8	382	for (int ih=histDet.GetEntriesFast();ih--;) ((TH1*)histDet.UncheckedAt(ih))->Fill(i);
	383	for (int ih=histTime.GetEntriesFast();ih--;) ((TH1*)histTime.UncheckedAt(ih))->Fill(time);
	384	for (int ih=histSignal.GetEntriesFast();ih--;) ((TH1*)histSignal.UncheckedAt(ih))->Fill(signal);
b65e5048	385	}
	386	}
	387	}
	388	// delete digitsIn;
	389	}
04236e67	390	delete digitsManager;
92664bc8	391	//
	392	IncEvCountCycleSDigits();
	393	IncEvCountTotalSDigits();
	394	//
04236e67	395	}
	396
	397	//____________________________________________________________________________
	398	void AliTRDQADataMakerSim::StartOfDetectorCycle()
	399	{
	400	//
	401	// Detector specific actions at start of cycle
	402	//
	403
	404	}
	405
	406	//__________________________________________________________________________
	407	Int_t AliTRDQADataMakerSim::CheckPointer(TObject obj, const char name)
	408	{
	409	//
	410	// Checks initialization of pointers
	411	//
	412
	413	if (!obj) AliWarning(Form("null pointer: %s", name));
	414	return !!obj;
	415
	416	}