[u/mrichter/AliRoot.git] / ACORDE / AliACORDEQADataMakerSim.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.                  *
 **************************************************************************/


//---
//  Produces the data needed to calculate the quality assurance. 
//  All data must be mergeable objects.

//  Authors:
//
//  Luciano Diaz Gonzalez <luciano.diaz@nucleares.unam.mx> (ICN-UNAM)
//  Mario Rodriguez Cahuantzi <mrodrigu@mail.cern.ch> (FCFM-BUAP)
//  Arturo Fernandez Tellez <afernan@mail.cern.ch (FCFM-BUAP)
//
//  Created: June 13th 2008
//---
// Last Update: Aug. 27th 2008 ---> Implementation to declare QA expert histogram 


// --- ROOT system ---
#include <TClonesArray.h>
#include <TFile.h> 
#include <TH1F.h> 
#include <TH2F.h>
#include <TDirectory.h>
// --- Standard library ---

// --- AliRoot header files ---
#include "AliESDEvent.h"
#include "AliLog.h"
#include "AliACORDEdigit.h" 
#include "AliACORDEhit.h"
#include "AliACORDEQADataMakerSim.h"
#include "AliQAChecker.h"
#include "AliACORDERawReader.h"
ClassImp(AliACORDEQADataMakerSim)
           
//____________________________________________________________________________ 
AliACORDEQADataMakerSim::AliACORDEQADataMakerSim():AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kACORDE), "ACORDE Quality Assurance Data Maker")
{
}
//____________________________________________________________________________ 
AliACORDEQADataMakerSim::AliACORDEQADataMakerSim(const AliACORDEQADataMakerSim& qadm) :
  AliQADataMakerSim() 
{
  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle()); 
}
//__________________________________________________________________
AliACORDEQADataMakerSim& AliACORDEQADataMakerSim::operator = (const AliACORDEQADataMakerSim& qadm )
{
  // Equal operator.
  this->~AliACORDEQADataMakerSim();
  new(this) AliACORDEQADataMakerSim(qadm);
  return *this;
}
//____________________________________________________________________________
void AliACORDEQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
{
  //Detector specific actions at end of cycle
  // do the QA checking
  ResetEventTrigClasses(); // reset triggers list to select all histos
  AliDebug(AliQAv1::GetQADebugLevel(), "ACORDE---->Detector specific actions at END of cycle\n................\n");

  AliQAChecker::Instance()->Run(AliQAv1::kACORDE, task, list) ;
}
//____________________________________________________________________________
void AliACORDEQADataMakerSim::StartOfDetectorCycle()
{
  //Detector specific actions at start of cycle
  AliDebug(AliQAv1::GetQADebugLevel(), "ACORDE---->Detector specific actions at START of cycle\n................\n");
}
//____________________________________________________________________________ 
void AliACORDEQADataMakerSim::InitHits()
{
  // create Hits histograms in Hits subdir
  
  const Bool_t expert   = kTRUE ; 
  const Bool_t image    = kTRUE ; 
  
  TH1F * fHitsACORDE;  
  
  fHitsACORDE = new TH1F("ACORDEBitPatternfromHits","Distribution of ACORDE fired modules from HITS",60,1,60);
  Add2HitsList(fHitsACORDE,0,!expert,image);
  
  const char *acoModule[60]={"0_0","0_1","0_2","0_3","0_4","0_5","0_6","0_7","0_8","0_9",
			     "1_0","1_1","1_2","1_3","1_4","1_5","1_6","1_7","1_8","1_9",
			     "2_0","2_1","2_2","2_3","2_4","2_5","2_6","2_7","2_8","2_9",
			     "3_0","3_1","3_2","3_3","3_4","3_5","3_6","3_7","3_8","3_9",
			     "4_0","4_1","4_2","4_3","4_4","4_5","4_6","4_7","4_8","4_9",
			     "5_0","5_1","5_2","5_3","5_4","5_5","5_6","5_7","5_8","5_9"};
  
  
  fHitsACORDE->SetXTitle("Modules");
  fHitsACORDE->SetYTitle("Counts");
  for (Int_t i=0;i<60;i++)
    {
      fHitsACORDE->GetXaxis()->SetBinLabel(i+1,acoModule[i]);
    }
  //
  //
  ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
}

//____________________________________________________________________________ 
void AliACORDEQADataMakerSim::InitDigits()
{
  // create Digits histograms in Digits subdir

  const Bool_t expert   = kTRUE ; 
  const Bool_t image    = kTRUE ; 
  
  TH1F *    fhDigitsModule;
  fhDigitsModule = new TH1F("ACORDEBitPatternfromDigits","Distribution of ACORDE fired modules from DIGITS",60,1,60);
  Add2DigitsList(fhDigitsModule,0,!expert,image);
  const char *acoModule[60]={"0_0","0_1","0_2","0_3","0_4","0_5","0_6","0_7","0_8","0_9",
                        "1_0","1_1","1_2","1_3","1_4","1_5","1_6","1_7","1_8","1_9",
                        "2_0","2_1","2_2","2_3","2_4","2_5","2_6","2_7","2_8","2_9",
                        "3_0","3_1","3_2","3_3","3_4","3_5","3_6","3_7","3_8","3_9",
                        "4_0","4_1","4_2","4_3","4_4","4_5","4_6","4_7","4_8","4_9",
                        "5_0","5_1","5_2","5_3","5_4","5_5","5_6","5_7","5_8","5_9"};
  
  
  fhDigitsModule->SetXTitle("Modules");
  fhDigitsModule->SetYTitle("Counts");
  for (Int_t i=0;i<60;i++)
    {
      fhDigitsModule->GetXaxis()->SetBinLabel(i+1,acoModule[i]);
    }
    //
  ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
}

//____________________________________________________________________________
void AliACORDEQADataMakerSim::MakeHits(TTree *hitTree)
{
  // Here we fill the QA histos for Hits declared above
  IncEvCountCycleHits();
  IncEvCountTotalHits();

	if (fHitsArray) 
    	fHitsArray->Clear() ; 
  	else
    	fHitsArray = new TClonesArray("AliACORDEhit",1000);
	TBranch * branch = hitTree->GetBranch("ACORDE");
	if (!branch) 
	{
		AliWarning("ACORDE branch in Hit Tree not found");
	} else 
	{
    		branch->SetAddress(&fHitsArray);
		for(Int_t track = 0 ; track < branch->GetEntries() ; track++) 
		{
			branch->GetEntry(track);
			for(Int_t ihit=0 ; ihit < fHitsArray->GetEntriesFast() ; ihit++) 
			{
				AliACORDEhit *AcoHit = (AliACORDEhit*) fHitsArray->UncheckedAt(ihit);
				if(!AcoHit) 
				{
					AliError("The unchecked hit doesn't exist");
					continue ;
				}
				FillHitsData(0,AcoHit->GetModule());
      			}			
    		}
  	}

}
//____________________________________________________________________________
void AliACORDEQADataMakerSim::MakeDigits( TTree *digitsTree)
{
  //fills QA histos for Digits
  IncEvCountCycleDigits();
  IncEvCountTotalDigits();

  if (fDigitsArray) 
    fDigitsArray->Clear() ; 
  else
    fDigitsArray = new TClonesArray("AliACORDEdigit",1000);
  
  TBranch * branch = digitsTree->GetBranch("ACORDEdigit");
  if (!branch) {
    AliWarning("ACORDE branch in Digits Tree not found");
  } else {
   branch->SetAddress(&fDigitsArray);
    for(Int_t track = 0 ; track < branch->GetEntries() ; track++) {
      branch->GetEntry(track);
      for(Int_t idigit = 0 ; idigit < fDigitsArray->GetEntriesFast() ; idigit++) {
        AliACORDEdigit *AcoDigit = (AliACORDEdigit*) fDigitsArray->UncheckedAt(idigit);
        if (!AcoDigit) {
          AliError("The unchecked digit doesn't exist");
          continue ;
        }
        FillDigitsData(0,AcoDigit->GetModule());
      }
    }
  }
}
Commit	Line	Data
51504028	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
	18	//---
	19	// Produces the data needed to calculate the quality assurance.
	20	// All data must be mergeable objects.
	21
	22	// Authors:
	23	//
	24	// Luciano Diaz Gonzalez <luciano.diaz@nucleares.unam.mx> (ICN-UNAM)
	25	// Mario Rodriguez Cahuantzi <mrodrigu@mail.cern.ch> (FCFM-BUAP)
	26	// Arturo Fernandez Tellez <afernan@mail.cern.ch (FCFM-BUAP)
	27	//
	28	// Created: June 13th 2008
	29	//---
491ceb3e	30	// Last Update: Aug. 27th 2008 ---> Implementation to declare QA expert histogram
491ceb3e	31
51504028	32
	33	// --- ROOT system ---
	34	#include <TClonesArray.h>
	35	#include <TFile.h>
	36	#include <TH1F.h>
	37	#include <TH2F.h>
	38	#include <TDirectory.h>
	39	// --- Standard library ---
	40
	41	// --- AliRoot header files ---
	42	#include "AliESDEvent.h"
	43	#include "AliLog.h"
	44	#include "AliACORDEdigit.h"
	45	#include "AliACORDEhit.h"
	46	#include "AliACORDEQADataMakerSim.h"
	47	#include "AliQAChecker.h"
	48	#include "AliACORDERawReader.h"
	49	ClassImp(AliACORDEQADataMakerSim)
	50
	51	//____________________________________________________________________________
4e25ac79	52	AliACORDEQADataMakerSim::AliACORDEQADataMakerSim():AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kACORDE), "ACORDE Quality Assurance Data Maker")
51504028	53	{
	54	}
	55	//____________________________________________________________________________
	56	AliACORDEQADataMakerSim::AliACORDEQADataMakerSim(const AliACORDEQADataMakerSim& qadm) :
	57	AliQADataMakerSim()
	58	{
	59	SetName((const char*)qadm.GetName()) ;
	60	SetTitle((const char*)qadm.GetTitle());
	61	}
	62	//__________________________________________________________________
	63	AliACORDEQADataMakerSim& AliACORDEQADataMakerSim::operator = (const AliACORDEQADataMakerSim& qadm )
	64	{
	65	// Equal operator.
	66	this->~AliACORDEQADataMakerSim();
	67	new(this) AliACORDEQADataMakerSim(qadm);
	68	return *this;
	69	}
	70	//____________________________________________________________________________
4e25ac79	71	void AliACORDEQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
51504028	72	{
	73	//Detector specific actions at end of cycle
	74	// do the QA checking
92664bc8	75	ResetEventTrigClasses(); // reset triggers list to select all histos
5379c4a3	76	AliDebug(AliQAv1::GetQADebugLevel(), "ACORDE---->Detector specific actions at END of cycle\n................\n");
51504028	77
4e25ac79	78	AliQAChecker::Instance()->Run(AliQAv1::kACORDE, task, list) ;
51504028	79	}
	80	//____________________________________________________________________________
	81	void AliACORDEQADataMakerSim::StartOfDetectorCycle()
	82	{
	83	//Detector specific actions at start of cycle
5379c4a3	84	AliDebug(AliQAv1::GetQADebugLevel(), "ACORDE---->Detector specific actions at START of cycle\n................\n");
51504028	85	}
	86	//____________________________________________________________________________
	87	void AliACORDEQADataMakerSim::InitHits()
	88	{
	89	// create Hits histograms in Hits subdir
92664bc8	90
7d297381	91	const Bool_t expert = kTRUE ;
7d297381	92	const Bool_t image = kTRUE ;
1093dbd7	93
1093dbd7	94	TH1F * fHitsACORDE;
92664bc8	95
	96	fHitsACORDE = new TH1F("ACORDEBitPatternfromHits","Distribution of ACORDE fired modules from HITS",60,1,60);
	97	Add2HitsList(fHitsACORDE,0,!expert,image);
	98
	99	const char *acoModule[60]={"0_0","0_1","0_2","0_3","0_4","0_5","0_6","0_7","0_8","0_9",
	100	"1_0","1_1","1_2","1_3","1_4","1_5","1_6","1_7","1_8","1_9",
	101	"2_0","2_1","2_2","2_3","2_4","2_5","2_6","2_7","2_8","2_9",
	102	"3_0","3_1","3_2","3_3","3_4","3_5","3_6","3_7","3_8","3_9",
	103	"4_0","4_1","4_2","4_3","4_4","4_5","4_6","4_7","4_8","4_9",
	104	"5_0","5_1","5_2","5_3","5_4","5_5","5_6","5_7","5_8","5_9"};
	105
	106
	107	fHitsACORDE->SetXTitle("Modules");
	108	fHitsACORDE->SetYTitle("Counts");
	109	for (Int_t i=0;i<60;i++)
	110	{
	111	fHitsACORDE->GetXaxis()->SetBinLabel(i+1,acoModule[i]);
	112	}
	113	//
	114	//
	115	ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
51504028	116	}
92664bc8	117
51504028	118	//____________________________________________________________________________
	119	void AliACORDEQADataMakerSim::InitDigits()
	120	{
	121	// create Digits histograms in Digits subdir
	122
7d297381	123	const Bool_t expert = kTRUE ;
	124	const Bool_t image = kTRUE ;
	125
	126	TH1F * fhDigitsModule;
1093dbd7	127	fhDigitsModule = new TH1F("ACORDEBitPatternfromDigits","Distribution of ACORDE fired modules from DIGITS",60,1,60);
7d297381	128	Add2DigitsList(fhDigitsModule,0,!expert,image);
41d0241f	129	const char *acoModule[60]={"0_0","0_1","0_2","0_3","0_4","0_5","0_6","0_7","0_8","0_9",
1093dbd7	130	"1_0","1_1","1_2","1_3","1_4","1_5","1_6","1_7","1_8","1_9",
	131	"2_0","2_1","2_2","2_3","2_4","2_5","2_6","2_7","2_8","2_9",
	132	"3_0","3_1","3_2","3_3","3_4","3_5","3_6","3_7","3_8","3_9",
	133	"4_0","4_1","4_2","4_3","4_4","4_5","4_6","4_7","4_8","4_9",
	134	"5_0","5_1","5_2","5_3","5_4","5_5","5_6","5_7","5_8","5_9"};
92664bc8	135
	136
	137	fhDigitsModule->SetXTitle("Modules");
	138	fhDigitsModule->SetYTitle("Counts");
	139	for (Int_t i=0;i<60;i++)
	140	{
	141	fhDigitsModule->GetXaxis()->SetBinLabel(i+1,acoModule[i]);
	142	}
	143	//
	144	ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
51504028	145	}
51504028	146
92664bc8	147	//____________________________________________________________________________
51504028	148	void AliACORDEQADataMakerSim::MakeHits(TTree *hitTree)
	149	{
	150	// Here we fill the QA histos for Hits declared above
92664bc8	151	IncEvCountCycleHits();
92664bc8	152	IncEvCountTotalHits();
51504028	153
6252ceeb	154	if (fHitsArray)
1093dbd7	155	fHitsArray->Clear() ;
	156	else
	157	fHitsArray = new TClonesArray("AliACORDEhit",1000);
d461c877	158	TBranch * branch = hitTree->GetBranch("ACORDE");
1093dbd7	159	if (!branch)
1093dbd7	160	{
d461c877	161	AliWarning("ACORDE branch in Hit Tree not found");
1093dbd7	162	} else
	163	{
	164	branch->SetAddress(&fHitsArray);
	165	for(Int_t track = 0 ; track < branch->GetEntries() ; track++)
	166	{
d461c877	167	branch->GetEntry(track);
1093dbd7	168	for(Int_t ihit=0 ; ihit < fHitsArray->GetEntriesFast() ; ihit++)
1093dbd7	169	{
6252ceeb	170	AliACORDEhit AcoHit = (AliACORDEhit) fHitsArray->UncheckedAt(ihit);
1093dbd7	171	if(!AcoHit)
1093dbd7	172	{
d461c877	173	AliError("The unchecked hit doesn't exist");
fcdfb1ff	174	continue ;
d461c877	175	}
92664bc8	176	FillHitsData(0,AcoHit->GetModule());
1093dbd7	177	}
	178	}
	179	}
	180
51504028	181	}
	182	//____________________________________________________________________________
	183	void AliACORDEQADataMakerSim::MakeDigits( TTree *digitsTree)
	184	{
	185	//fills QA histos for Digits
92664bc8	186	IncEvCountCycleDigits();
	187	IncEvCountTotalDigits();
	188
6252ceeb	189	if (fDigitsArray)
	190	fDigitsArray->Clear() ;
	191	else
	192	fDigitsArray = new TClonesArray("AliACORDEdigit",1000);
	193
fcdfb1ff	194	TBranch * branch = digitsTree->GetBranch("ACORDEdigit");
	195	if (!branch) {
	196	AliWarning("ACORDE branch in Digits Tree not found");
	197	} else {
6252ceeb	198	branch->SetAddress(&fDigitsArray);
fcdfb1ff	199	for(Int_t track = 0 ; track < branch->GetEntries() ; track++) {
fcdfb1ff	200	branch->GetEntry(track);
6252ceeb	201	for(Int_t idigit = 0 ; idigit < fDigitsArray->GetEntriesFast() ; idigit++) {
6252ceeb	202	AliACORDEdigit AcoDigit = (AliACORDEdigit) fDigitsArray->UncheckedAt(idigit);
fcdfb1ff	203	if (!AcoDigit) {
	204	AliError("The unchecked digit doesn't exist");
	205	continue ;
51504028	206	}
92664bc8	207	FillDigitsData(0,AcoDigit->GetModule());
fcdfb1ff	208	}
	209	}
	210	}
51504028	211	}