[u/mrichter/AliRoot.git] / ACORDE / AliACORDEQADataMakerRec.cxx

/**************************************************************************\r
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
 *                                                                        *\r
 * Author: The ALICE Off-line Project.                                    *\r
 * Contributors are mentioned in the code where appropriate.              *\r
 *                                                                        *\r
 * Permission to use, copy, modify and distribute this software and its   *\r
 * documentation strictly for non-commercial purposes is hereby granted   *\r
 * without fee, provided that the above copyright notice appears in all   *\r
 * copies and that both the copyright notice and this permission notice   *\r
 * appear in the supporting documentation. The authors make no claims     *\r
 * about the suitability of this software for any purpose. It is          *\r
 * provided "as is" without express or implied warranty.                  *\r
 **************************************************************************/\r
//---\r
//  Produces the data needed to calculate the quality assurance. \r
//  All data must be mergeable objects.\r
\r
\r
//  Authors:\r
//\r
//  Luciano Diaz Gonzalez <luciano.diaz@nucleares.unam.mx> (ICN-UNAM)\r
//  Mario Rodriguez Cahuantzi <mrodrigu@mail.cern.ch> (FCFM-BUAP)\r
//  Arturo Fernandez Tellez <afernan@mail.cern.ch (FCFM-BUAP)\r
//\r
//  Created: June 13th 2008\r
//---\r
// Last update: Sept. 29th. 2011 (by Mario RC: mrodrigu@mail.cern.ch) \r
//	--> ACOMultiSL0_DQM_Shifter filling histogram fixed\r
//	--> Expert histogram updated: 2 histograms (Checks the hits for SL0 and AMU mode)\r
//	--> To be include in the next update: threshold settings from AliACORDEQAThreshold class (not yet)\r
// Last update: May 5th. 2011 (by Mario RC: mrodrigu@mail.cern.ch) -->Creates QA expert histograms \r
// and QA-shifter histograms also with threshold lines and visual alarm\r
// Last Update: Aug. 27th 2008 --> Implementation to declare QA expert histogram\r
// Last update: Nov. 14t 2009 --> MRC <mrodrigu@mail.cern.ch> (FCFM-BUAP) \r
//			|--> Change in Multiplicity histogram for AMORE (to detect empty triggers events of ACORDE)\r
\r
\r
\r
// --- ROOT system ---\r
#include <TClonesArray.h>\r
#include <TFile.h> \r
#include <TH1F.h> \r
#include <TDirectory.h>\r
#include <TPaveText.h>\r
// --- Standard library ---\r
\r
// --- AliRoot header files ---\r
#include "AliESDEvent.h"\r
#include "AliLog.h"\r
#include "AliACORDEdigit.h" \r
#include "AliACORDEhit.h"\r
#include "AliACORDEQADataMakerRec.h"\r
#include "AliQAChecker.h"\r
#include "AliACORDERawReader.h"\r
#include "AliACORDERawStream.h"\r
ClassImp(AliACORDEQADataMakerRec)\r
           \r
//____________________________________________________________________________ \r
 AliACORDEQADataMakerRec::AliACORDEQADataMakerRec():AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kACORDE), "ACORDE Quality Assurance Data Maker")\r
{\r
\r
}\r
//____________________________________________________________________________ \r
AliACORDEQADataMakerRec::AliACORDEQADataMakerRec(const AliACORDEQADataMakerRec& qadm):\r
  AliQADataMakerRec()\r
{\r
  SetName((const char*)qadm.GetName()) ; \r
  SetTitle((const char*)qadm.GetTitle()); \r
}\r
\r
//__________________________________________________________________\r
AliACORDEQADataMakerRec::~AliACORDEQADataMakerRec()\r
{\r
}\r
\r
//__________________________________________________________________\r
AliACORDEQADataMakerRec& AliACORDEQADataMakerRec::operator = (const AliACORDEQADataMakerRec& qadm )\r
{\r
  // Equal operator.\r
  this->~AliACORDEQADataMakerRec();\r
  new(this) AliACORDEQADataMakerRec(qadm);\r
  return *this;\r
}\r
//____________________________________________________________________________\r
void AliACORDEQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)\r
{\r
  //Detector specific actions at end of cycle\r
  // do the QA checking\r
  ResetEventTrigClasses(); // reset triggers list to select all histos\r
  // Update for DQM GUI\r
  //\r
\r
  for (Int_t specie = 0; specie < AliRecoParam::kNSpecies ; specie++) {\r
    if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) ) continue ;\r
    // \r
    // RS Set event specie\r
    SetEventSpecie(AliRecoParam::ConvertIndex(specie));\r
    //\r
    for (int itc=-1;itc<GetNTrigClasses();itc++) { // RS Loop over the trigger classes\r
      //\r
      	TObjArray * parr = GetRawsDataOfTrigClass(itc);\r
      	if (!parr) continue;\r
      	TObjArray &harr = *parr;\r
      	TH1* h0 = (TH1*)harr[0];\r
      	TH1* h2 = (TH1*)harr[2];\r
      	if (!h0 || !h2) continue;\r
\r
      	Double_t integralSL0 = 0;\r
      	Double_t integralAMU = 0;\r
\r
	if ((itc==-1) && ( (integralAMU=h2->Integral()==0) || (integralSL0=h0->Integral()==0) ) ) continue;\r
\r
	Float_t maxSL0 = 1.1*h0->GetMaximum();\r
	Float_t scaleSL0 = 0.;\r
	if (maxSL0!=0) scaleSL0 = 1./maxSL0;\r
	else scaleSL0 = 1.;\r
	h0->Scale(scaleSL0);\r
\r
\r
	Float_t maxAMU = 1.1*h2->GetMaximum();\r
	Float_t scaleAMU = 0.;\r
	if (maxAMU!=0) scaleAMU = 1./maxAMU;\r
	else scaleAMU = 1.;\r
	h2->Scale(scaleAMU);\r
      \r
	Int_t indexActiveModuleSL0 = 0;\r
	Int_t indexActiveModuleAMU = 0;\r
\r
	for(Int_t iModule = 0; iModule < 60; iModule++){\r
		if (h0->GetBinContent(iModule) > 0) indexActiveModuleSL0++;\r
		if (h2->GetBinContent(iModule) > 0) indexActiveModuleAMU++;\r
	}\r
	\r
	Float_t meanHitsSL0 = 0.;\r
	Float_t meanHitsAMU = 0.;\r
	if ((indexActiveModuleSL0==0) || (indexActiveModuleAMU == 0)) continue;\r
\r
	meanHitsSL0 = h0->Integral()/indexActiveModuleSL0;\r
	meanHitsAMU = h2->Integral()/indexActiveModuleAMU;\r
\r
        TH1* h4 = (TH1*)harr[4];\r
        TH1* h5 = (TH1*)harr[5];\r
\r
	if (h4 && h5 && meanHitsAMU!=0 && meanHitsSL0!=0){\r
		for (Int_t iModule = 0; iModule < 60; iModule++){\r
			h4->Fill(h0->GetBinContent(iModule)/meanHitsSL0-1);\r
			h5->Fill(h2->GetBinContent(iModule)/meanHitsAMU-1);\r
		}\r
	}\r
\r
    } // end of trigger classes loop\r
  } // end specie loop\r
  // QA Checker standar (to be updated)\r
  //\r
  AliQAChecker::Instance()->Run(AliQAv1::kACORDE, task, list) ;\r
}\r
\r
//____________________________________________________________________________\r
void AliACORDEQADataMakerRec::StartOfDetectorCycle()\r
{\r
  //Detector specific actions at start of cycle\r
\r
}\r
 \r
//____________________________________________________________________________ \r
void AliACORDEQADataMakerRec::InitRaws()\r
{\r
  // create Raw histograms in Raw subdir\r
\r
  	const Bool_t expert   = kTRUE ; \r
  	const Bool_t saveCorr = kTRUE ; \r
  	const Bool_t image    = kTRUE ; \r
	TH1F * fhACORDEBitPatternDQM = new TH1F("ACOHitsSL0_DQM_Shifter","Distribution of ACORDE fired modules for DQM shifter; No. of module; Counts",60,-0.5,59.5);// Hits histogram for QA-shifter ACO-SL0 trigger mode\r
        TH1F * fhACORDEMultiplicitySL0DQM = new TH1F("ACOMultiSL0_DQM_Shifter","Multiplicity of ACORDE fired modules for DQM shifter; No. of fired modules; No. of events",62,-1,60); // Multiplicity histo. for QA-shifter ACO-SL0 trigger mode\r
        TH1F * fhACORDEBitPatternAMUDQM = new TH1F("ACOHitsAMU_DQM_Shifter","Distribution of ACORDE fired modules for DQM shifter; No. of module; Counts",60,-0.5,59.5);// Hits histogram for QA-shifter ACO-SL0 trigger mode\r
        TH1F * fhACORDEMultiplicityAMUDQM = new TH1F("ACOMultiAMU_DQM_Shifter","Multiplicity of ACORDE fired modules for DQM shifter; No. of fired modules; No. of events",62,-1,60); // Multiplicity histo. for QA-shifter ACO-SL0 trigger mode\r
        TH1F * fhACORDEBitPatternCheckDQMSL0 = new TH1F("ACOHitsTriggerCheck_DQMExpertSL0","Check the activity of ACORDE's modules; Hits per module/mean of Hits; Counts",100,-3,5); // Check the trigger status of ACORDE (SL0 vs AMU)\r
        TH1F * fhACORDEBitPatternCheckDQMAMU = new TH1F("ACOHitsTriggerCheck_DQMExpertAMU","Check the activity of ACORDE's modules; Hits per module/mean of Hits; Counts",100,-3,5); // Check the trigger status of ACORDE (SL0 vs AMU)\r
         // Expert histograms\r
         // Check the hits multiplicity from trigger configuration\r
         Add2RawsList(fhACORDEBitPatternCheckDQMSL0,4,expert,image,!saveCorr);\r
         fhACORDEBitPatternCheckDQMSL0->SetFillColor(kOrange);\r
	 Add2RawsList(fhACORDEBitPatternCheckDQMAMU,5,expert,image,!saveCorr);\r
         fhACORDEBitPatternCheckDQMAMU->SetFillColor(kRed+2);\r
\r
	\r
	// AMORE diplay settings for shifter on GUI\r
 \r
        // For SL0 ACO trigger mode\r
 \r
         Add2RawsList(fhACORDEBitPatternDQM,0,!expert,image,!saveCorr);\r
	 ForbidCloning(fhACORDEBitPatternDQM);\r
         Add2RawsList(fhACORDEMultiplicitySL0DQM,1,!expert,image,!saveCorr);\r
 	 ForbidCloning(fhACORDEMultiplicitySL0DQM);\r
\r
         // For Hits distribution on ACORDE\r
 \r
         fhACORDEBitPatternDQM->SetFillColor(kMagenta+2);\r
 \r
         // For ACORDE Multiplicity distribution of fired modules\r
 \r
         fhACORDEMultiplicitySL0DQM->SetFillColor(kMagenta);\r
 \r
         // For AMU ACO trigger mode\r
 \r
         Add2RawsList(fhACORDEBitPatternAMUDQM,2,!expert,image,!saveCorr);\r
         Add2RawsList(fhACORDEMultiplicityAMUDQM,3,!expert,image,!saveCorr);\r
 \r
         // For Hits distribution on ACORDE\r
 \r
         fhACORDEBitPatternAMUDQM->SetFillColor(kViolet+7);\r
 \r
         // For ACORDE Multiplicity distribution of fired modules\r
 \r
         fhACORDEMultiplicityAMUDQM->SetFillColor(kViolet+6);\r
 \r
  //\r
  ClonePerTrigClass(AliQAv1::kRAWS); // this should be the last line\r
}\r
\r
//____________________________________________________________________________ \r
void AliACORDEQADataMakerRec::InitDigits()\r
{\r
  // create Digits histograms in Digits subdir\r
  \r
  const Bool_t expert   = kTRUE ; \r
  const Bool_t image    = kTRUE ; \r
  TH1F *    fhDigitsModule;\r
  const char *acoModule[60]={"0_0","0_1","0_2","0_3","0_4","0_5","0_6","0_7","0_8","0_9",\r
                        "1_0","1_1","1_2","1_3","1_4","1_5","1_6","1_7","1_8","1_9",\r
                        "2_0","2_1","2_2","2_3","2_4","2_5","2_6","2_7","2_8","2_9",\r
                        "3_0","3_1","3_2","3_3","3_4","3_5","3_6","3_7","3_8","3_9",\r
                        "4_0","4_1","4_2","4_3","4_4","4_5","4_6","4_7","4_8","4_9",\r
                        "5_0","5_1","5_2","5_3","5_4","5_5","5_6","5_7","5_8","5_9"};\r
\r
\r
  fhDigitsModule = new TH1F("ACORDEBitPatternfromDigits","Distribution of ACORDE from DIGITS;Modules;Counts",60,1,60);\r
  Add2DigitsList(fhDigitsModule,0,!expert,image);\r
  for (Int_t i=0;i<60;i++) fhDigitsModule->GetXaxis()->SetBinLabel(i+1,acoModule[i]); \r
  //\r
  ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line\r
}\r
\r
//____________________________________________________________________________ \r
\r
void AliACORDEQADataMakerRec::InitRecPoints()\r
{\r
  // create cluster histograms in RecPoint subdir\r
  // Not needed for ACORDE by now !!!\r
  //\r
  ClonePerTrigClass(AliQAv1::kRECPOINTS); // this should be the last line\r
}\r
\r
//____________________________________________________________________________\r
void AliACORDEQADataMakerRec::InitESDs()\r
{\r
  //create ESDs histograms in ESDs subdir\r
\r
  const Bool_t expert   = kTRUE ; \r
  const Bool_t image    = kTRUE ; \r
  \r
  TH1F *    fhESDsSingle;\r
  TH1F *    fhESDsMulti;\r
  TH1F *	fhESDsMultiplicity;\r
  const char *acoModule[60]={"0_0","0_1","0_2","0_3","0_4","0_5","0_6","0_7","0_8","0_9",\r
                        "1_0","1_1","1_2","1_3","1_4","1_5","1_6","1_7","1_8","1_9",\r
                        "2_0","2_1","2_2","2_3","2_4","2_5","2_6","2_7","2_8","2_9",\r
                        "3_0","3_1","3_2","3_3","3_4","3_5","3_6","3_7","3_8","3_9",\r
                        "4_0","4_1","4_2","4_3","4_4","4_5","4_6","4_7","4_8","4_9",\r
                        "5_0","5_1","5_2","5_3","5_4","5_5","5_6","5_7","5_8","5_9"};\r
\r
\r
   fhESDsSingle = new TH1F("ACORDEBitPatternfromESDsSingle","Distribution of ACORDE fired modules from ESDs-Single;Modules;Counts",60,1,60);\r
   Add2ESDsList(fhESDsSingle,0,!expert,image);\r
\r
   fhESDsMulti = new TH1F("ACORDEBitPatternfromESDsMulti","Distribution of ACORDE fired modules from ESDs-Multi;Modules;Counts",60,1,60);\r
   Add2ESDsList(fhESDsMulti,1,!expert,image);\r
   \r
   fhESDsMultiplicity = new TH1F("ACORDEMultiplicityfromESD","Number of fired ACORDE modules; No. of fired ACORDE modules;No. of events in ACORDE",60,-0.5,60);\r
   Add2ESDsList(fhESDsMultiplicity,2,!expert,image);	\r
   for (Int_t i=0;i<60;i++)\r
   {\r
	fhESDsSingle->GetXaxis()->SetBinLabel(i+1,acoModule[i]);\r
	fhESDsMulti->GetXaxis()->SetBinLabel(i+1,acoModule[i]);\r
   }\r
   //\r
   ClonePerTrigClass(AliQAv1::kESDS); // this should be the last line\r
}\r
//____________________________________________________________________________\r
void AliACORDEQADataMakerRec::MakeRaws(AliRawReader* rawReader)\r
{\r
  //fills QA histos for RAW\r
  IncEvCountCycleRaws();\r
  IncEvCountTotalRaws();\r
  rawReader->Reset();\r
  AliACORDERawStream rawStream(rawReader);\r
  size_t contSingle=0;\r
  size_t contMulti=0;\r
  UInt_t dy[4];\r
\r
  bool kroSingle[60],kroMulti[60];\r
  UInt_t tmpDy;\r
\r
  for(Int_t m=0;m<60;m++) {kroSingle[m]=0;kroMulti[m]=0;}\r
\r
if(rawStream.Next())\r
{\r
        dy[0]=rawStream.GetWord(0);\r
        dy[1]=rawStream.GetWord(1);\r
        dy[2]=rawStream.GetWord(2);\r
        dy[3]=rawStream.GetWord(3);\r
        tmpDy=dy[0];\r
        for(Int_t r=0;r<30;r++)\r
        {\r
                kroSingle[r] = tmpDy & 1;\r
                tmpDy>>=1;\r
        }\r
        tmpDy=dy[1];\r
        for(Int_t r=30;r<60;r++)\r
        {\r
                kroSingle[r] = tmpDy & 1;\r
                tmpDy>>=1;\r
        }\r
        tmpDy=dy[2];\r
        for(Int_t r=0;r<30;r++)\r
        {\r
                kroMulti[r] = tmpDy & 1;\r
                tmpDy>>=1;\r
        }\r
        tmpDy=dy[3];\r
        for(Int_t r=30;r<60;r++)\r
        {\r
                kroMulti[r] = tmpDy & 1;\r
                tmpDy>>=1;\r
        }\r
        contSingle=0;\r
	contMulti=0;\r
        for(Int_t r=0;r<60;r++)\r
        {\r
			if(kroSingle[r]==1)\r
			{\r
			  FillRawsData(0,r);\r
			  //FillRawsData(4,r);\r
			  contSingle++;\r
			}\r
			if(kroMulti[r]==1)\r
			{\r
			  FillRawsData(2,r);\r
			  //FillRawsData(6,r);\r
			  contMulti++;\r
			}\r
			\r
        } \r
	FillRawsData(1,contSingle); \r
	FillRawsData(3,contMulti); \r
	//	FillRawsData(7,contMulti);\r
}\r
}\r
//____________________________________________________________________________\r
void AliACORDEQADataMakerRec::MakeDigits( TTree *digitsTree)\r
{\r
  //fills QA histos for Digits\r
  IncEvCountCycleDigits();\r
  IncEvCountTotalDigits();\r
\r
  if (fDigitsArray) \r
    fDigitsArray->Clear() ; \r
  else \r
    fDigitsArray = new TClonesArray("AliACORDEdigit",1000);\r
  TBranch * branch = digitsTree->GetBranch("ACORDEdigit");\r
  if (!branch) {\r
    AliWarning("ACORDE branch in Digits Tree not found");\r
  } else {\r
    branch->SetAddress(&fDigitsArray);\r
    for(Int_t track = 0 ; track < branch->GetEntries() ; track++) {\r
      branch->GetEntry(track);\r
      for(Int_t idigit = 0 ; idigit < fDigitsArray->GetEntriesFast() ; idigit++) {\r
        AliACORDEdigit *AcoDigit = (AliACORDEdigit*) fDigitsArray->UncheckedAt(idigit);\r
        if (!AcoDigit) {\r
          AliError("The unchecked digit doesn't exist");\r
          continue ;\r
        }\r
        FillDigitsData(0,AcoDigit->GetModule()-1);\r
      }\r
    }\r
  }\r
}\r
\r
//____________________________________________________________________________\r
void AliACORDEQADataMakerRec::MakeESDs(AliESDEvent * esd)\r
{\r
  //fills QA histos for ESD\r
  IncEvCountCycleESDs();\r
  IncEvCountTotalESDs();\r
  //\r
  AliESDACORDE * fESDACORDE= esd->GetACORDEData();\r
  Int_t acoMulti=0;\r
  for(int i=0;i<60;i++)\r
  {\r
    if(fESDACORDE->GetHitChannel(i)) \r
	  {\r
	  FillESDsData(0,i+1);\r
	  FillESDsData(1,i+1);\r
	  acoMulti++;\r
	}\r
  } FillESDsData(2,acoMulti);\r
\r
}\r
Commit	Line	Data
dc7f1e9e	1	/**************************************************************************\r
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
	3	* *\r
	4	* Author: The ALICE Off-line Project. *\r
	5	* Contributors are mentioned in the code where appropriate. *\r
	6	* *\r
	7	* Permission to use, copy, modify and distribute this software and its *\r
	8	* documentation strictly for non-commercial purposes is hereby granted *\r
	9	* without fee, provided that the above copyright notice appears in all *\r
	10	* copies and that both the copyright notice and this permission notice *\r
	11	* appear in the supporting documentation. The authors make no claims *\r
	12	* about the suitability of this software for any purpose. It is *\r
	13	* provided "as is" without express or implied warranty. *\r
	14	**************************************************************************/\r
	15	//---\r
	16	// Produces the data needed to calculate the quality assurance. \r
	17	// All data must be mergeable objects.\r
	18	\r
	19	\r
	20	// Authors:\r
	21	//\r
	22	// Luciano Diaz Gonzalez <luciano.diaz@nucleares.unam.mx> (ICN-UNAM)\r
	23	// Mario Rodriguez Cahuantzi <mrodrigu@mail.cern.ch> (FCFM-BUAP)\r
	24	// Arturo Fernandez Tellez <afernan@mail.cern.ch (FCFM-BUAP)\r
	25	//\r
	26	// Created: June 13th 2008\r
	27	//---\r
a6bcaa19	28	// Last update: Sept. 29th. 2011 (by Mario RC: mrodrigu@mail.cern.ch) \r
	29	// --> ACOMultiSL0_DQM_Shifter filling histogram fixed\r
	30	// --> Expert histogram updated: 2 histograms (Checks the hits for SL0 and AMU mode)\r
	31	// --> To be include in the next update: threshold settings from AliACORDEQAThreshold class (not yet)\r
92664bc8	32	// Last update: May 5th. 2011 (by Mario RC: mrodrigu@mail.cern.ch) -->Creates QA expert histograms \r
92664bc8	33	// and QA-shifter histograms also with threshold lines and visual alarm\r
dc7f1e9e	34	// Last Update: Aug. 27th 2008 --> Implementation to declare QA expert histogram\r
92664bc8	35	// Last update: Nov. 14t 2009 --> MRC <mrodrigu@mail.cern.ch> (FCFM-BUAP) \r
1b368017	36	// \|--> Change in Multiplicity histogram for AMORE (to detect empty triggers events of ACORDE)\r
1b368017	37	\r
dc7f1e9e	38	\r
	39	\r
	40	// --- ROOT system ---\r
	41	#include <TClonesArray.h>\r
	42	#include <TFile.h> \r
	43	#include <TH1F.h> \r
	44	#include <TDirectory.h>\r
7626fda4	45	#include <TPaveText.h>\r
dc7f1e9e	46	// --- Standard library ---\r
	47	\r
	48	// --- AliRoot header files ---\r
	49	#include "AliESDEvent.h"\r
	50	#include "AliLog.h"\r
	51	#include "AliACORDEdigit.h" \r
	52	#include "AliACORDEhit.h"\r
	53	#include "AliACORDEQADataMakerRec.h"\r
	54	#include "AliQAChecker.h"\r
	55	#include "AliACORDERawReader.h"\r
	56	#include "AliACORDERawStream.h"\r
	57	ClassImp(AliACORDEQADataMakerRec)\r
	58	\r
	59	//____________________________________________________________________________ \r
dd1da44f	60	AliACORDEQADataMakerRec::AliACORDEQADataMakerRec():AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kACORDE), "ACORDE Quality Assurance Data Maker")\r
dc7f1e9e	61	{\r
	62	\r
	63	}\r
	64	//____________________________________________________________________________ \r
92664bc8	65	AliACORDEQADataMakerRec::AliACORDEQADataMakerRec(const AliACORDEQADataMakerRec& qadm):\r
dd1da44f	66	AliQADataMakerRec()\r
dc7f1e9e	67	{\r
	68	SetName((const char*)qadm.GetName()) ; \r
	69	SetTitle((const char*)qadm.GetTitle()); \r
	70	}\r
92664bc8	71	\r
7626fda4	72	//__________________________________________________________________\r
	73	AliACORDEQADataMakerRec::~AliACORDEQADataMakerRec()\r
	74	{\r
7626fda4	75	}\r
7626fda4	76	\r
dc7f1e9e	77	//__________________________________________________________________\r
	78	AliACORDEQADataMakerRec& AliACORDEQADataMakerRec::operator = (const AliACORDEQADataMakerRec& qadm )\r
	79	{\r
	80	// Equal operator.\r
	81	this->~AliACORDEQADataMakerRec();\r
	82	new(this) AliACORDEQADataMakerRec(qadm);\r
	83	return *this;\r
	84	}\r
	85	//____________________________________________________________________________\r
	86	void AliACORDEQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)\r
	87	{\r
	88	//Detector specific actions at end of cycle\r
	89	// do the QA checking\r
92664bc8	90	ResetEventTrigClasses(); // reset triggers list to select all histos\r
	91	// Update for DQM GUI\r
	92	//\r
a6bcaa19	93	\r
92664bc8	94	for (Int_t specie = 0; specie < AliRecoParam::kNSpecies ; specie++) {\r
	95	if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) ) continue ;\r
	96	// \r
	97	// RS Set event specie\r
	98	SetEventSpecie(AliRecoParam::ConvertIndex(specie));\r
	99	//\r
	100	for (int itc=-1;itc<GetNTrigClasses();itc++) { // RS Loop over the trigger classes\r
	101	//\r
a6bcaa19	102	TObjArray * parr = GetRawsDataOfTrigClass(itc);\r
	103	if (!parr) continue;\r
	104	TObjArray &harr = *parr;\r
	105	TH1* h0 = (TH1*)harr[0];\r
a6bcaa19	106	TH1* h2 = (TH1*)harr[2];\r
dd1da44f	107	if (!h0 \|\| !h2) continue;\r
a6bcaa19	108	\r
	109	Double_t integralSL0 = 0;\r
	110	Double_t integralAMU = 0;\r
	111	\r
2024a260	112	if ((itc==-1) && ( (integralAMU=h2->Integral()==0) \|\| (integralSL0=h0->Integral()==0) ) ) continue;\r
a6bcaa19	113	\r
dd1da44f	114	Float_t maxSL0 = 1.1*h0->GetMaximum();\r
	115	Float_t scaleSL0 = 0.;\r
	116	if (maxSL0!=0) scaleSL0 = 1./maxSL0;\r
	117	else scaleSL0 = 1.;\r
	118	h0->Scale(scaleSL0);\r
a6bcaa19	119	\r
a6bcaa19	120	\r
dd1da44f	121	Float_t maxAMU = 1.1*h2->GetMaximum();\r
	122	Float_t scaleAMU = 0.;\r
	123	if (maxAMU!=0) scaleAMU = 1./maxAMU;\r
	124	else scaleAMU = 1.;\r
	125	h2->Scale(scaleAMU);\r
92664bc8	126	\r
a6bcaa19	127	Int_t indexActiveModuleSL0 = 0;\r
dd1da44f	128	Int_t indexActiveModuleAMU = 0;\r
	129	\r
	130	for(Int_t iModule = 0; iModule < 60; iModule++){\r
	131	if (h0->GetBinContent(iModule) > 0) indexActiveModuleSL0++;\r
	132	if (h2->GetBinContent(iModule) > 0) indexActiveModuleAMU++;\r
a6bcaa19	133	}\r
92664bc8	134	\r
a6bcaa19	135	Float_t meanHitsSL0 = 0.;\r
a6bcaa19	136	Float_t meanHitsAMU = 0.;\r
dd1da44f	137	if ((indexActiveModuleSL0==0) \|\| (indexActiveModuleAMU == 0)) continue;\r
a6bcaa19	138	\r
dd1da44f	139	meanHitsSL0 = h0->Integral()/indexActiveModuleSL0;\r
dd1da44f	140	meanHitsAMU = h2->Integral()/indexActiveModuleAMU;\r
a6bcaa19	141	\r
dd1da44f	142	TH1* h4 = (TH1*)harr[4];\r
a6bcaa19	143	TH1* h5 = (TH1*)harr[5];\r
a6bcaa19	144	\r
dd1da44f	145	if (h4 && h5 && meanHitsAMU!=0 && meanHitsSL0!=0){\r
	146	for (Int_t iModule = 0; iModule < 60; iModule++){\r
	147	h4->Fill(h0->GetBinContent(iModule)/meanHitsSL0-1);\r
	148	h5->Fill(h2->GetBinContent(iModule)/meanHitsAMU-1);\r
a6bcaa19	149	}\r
	150	}\r
	151	\r
92664bc8	152	} // end of trigger classes loop\r
	153	} // end specie loop\r
	154	// QA Checker standar (to be updated)\r
	155	//\r
dc7f1e9e	156	AliQAChecker::Instance()->Run(AliQAv1::kACORDE, task, list) ;\r
	157	}\r
	158	\r
	159	//____________________________________________________________________________\r
	160	void AliACORDEQADataMakerRec::StartOfDetectorCycle()\r
	161	{\r
	162	//Detector specific actions at start of cycle\r
	163	\r
	164	}\r
	165	\r
	166	//____________________________________________________________________________ \r
	167	void AliACORDEQADataMakerRec::InitRaws()\r
	168	{\r
	169	// create Raw histograms in Raw subdir\r
	170	\r
a6bcaa19	171	const Bool_t expert = kTRUE ; \r
	172	const Bool_t saveCorr = kTRUE ; \r
	173	const Bool_t image = kTRUE ; \r
	174	TH1F * fhACORDEBitPatternDQM = new TH1F("ACOHitsSL0_DQM_Shifter","Distribution of ACORDE fired modules for DQM shifter; No. of module; Counts",60,-0.5,59.5);// Hits histogram for QA-shifter ACO-SL0 trigger mode\r
	175	TH1F * fhACORDEMultiplicitySL0DQM = new TH1F("ACOMultiSL0_DQM_Shifter","Multiplicity of ACORDE fired modules for DQM shifter; No. of fired modules; No. of events",62,-1,60); // Multiplicity histo. for QA-shifter ACO-SL0 trigger mode\r
	176	TH1F * fhACORDEBitPatternAMUDQM = new TH1F("ACOHitsAMU_DQM_Shifter","Distribution of ACORDE fired modules for DQM shifter; No. of module; Counts",60,-0.5,59.5);// Hits histogram for QA-shifter ACO-SL0 trigger mode\r
	177	TH1F * fhACORDEMultiplicityAMUDQM = new TH1F("ACOMultiAMU_DQM_Shifter","Multiplicity of ACORDE fired modules for DQM shifter; No. of fired modules; No. of events",62,-1,60); // Multiplicity histo. for QA-shifter ACO-SL0 trigger mode\r
dd1da44f	178	TH1F * fhACORDEBitPatternCheckDQMSL0 = new TH1F("ACOHitsTriggerCheck_DQMExpertSL0","Check the activity of ACORDE's modules; Hits per module/mean of Hits; Counts",100,-3,5); // Check the trigger status of ACORDE (SL0 vs AMU)\r
dd1da44f	179	TH1F * fhACORDEBitPatternCheckDQMAMU = new TH1F("ACOHitsTriggerCheck_DQMExpertAMU","Check the activity of ACORDE's modules; Hits per module/mean of Hits; Counts",100,-3,5); // Check the trigger status of ACORDE (SL0 vs AMU)\r
92664bc8	180	// Expert histograms\r
92664bc8	181	// Check the hits multiplicity from trigger configuration\r
a6bcaa19	182	Add2RawsList(fhACORDEBitPatternCheckDQMSL0,4,expert,image,!saveCorr);\r
	183	fhACORDEBitPatternCheckDQMSL0->SetFillColor(kOrange);\r
	184	Add2RawsList(fhACORDEBitPatternCheckDQMAMU,5,expert,image,!saveCorr);\r
	185	fhACORDEBitPatternCheckDQMAMU->SetFillColor(kRed+2);\r
	186	\r
	187	\r
	188	// AMORE diplay settings for shifter on GUI\r
92664bc8	189	\r
a6bcaa19	190	// For SL0 ACO trigger mode\r
92664bc8	191	\r
	192	Add2RawsList(fhACORDEBitPatternDQM,0,!expert,image,!saveCorr);\r
	193	ForbidCloning(fhACORDEBitPatternDQM);\r
	194	Add2RawsList(fhACORDEMultiplicitySL0DQM,1,!expert,image,!saveCorr);\r
	195	ForbidCloning(fhACORDEMultiplicitySL0DQM);\r
a6bcaa19	196	\r
92664bc8	197	// For Hits distribution on ACORDE\r
92664bc8	198	\r
a6bcaa19	199	fhACORDEBitPatternDQM->SetFillColor(kMagenta+2);\r
92664bc8	200	\r
	201	// For ACORDE Multiplicity distribution of fired modules\r
	202	\r
a6bcaa19	203	fhACORDEMultiplicitySL0DQM->SetFillColor(kMagenta);\r
92664bc8	204	\r
	205	// For AMU ACO trigger mode\r
	206	\r
	207	Add2RawsList(fhACORDEBitPatternAMUDQM,2,!expert,image,!saveCorr);\r
	208	Add2RawsList(fhACORDEMultiplicityAMUDQM,3,!expert,image,!saveCorr);\r
	209	\r
	210	// For Hits distribution on ACORDE\r
	211	\r
a6bcaa19	212	fhACORDEBitPatternAMUDQM->SetFillColor(kViolet+7);\r
92664bc8	213	\r
	214	// For ACORDE Multiplicity distribution of fired modules\r
	215	\r
a6bcaa19	216	fhACORDEMultiplicityAMUDQM->SetFillColor(kViolet+6);\r
92664bc8	217	\r
92664bc8	218	//\r
92664bc8	219	ClonePerTrigClass(AliQAv1::kRAWS); // this should be the last line\r
dc7f1e9e	220	}\r
92664bc8	221	\r
dc7f1e9e	222	//____________________________________________________________________________ \r
	223	void AliACORDEQADataMakerRec::InitDigits()\r
	224	{\r
	225	// create Digits histograms in Digits subdir\r
	226	\r
	227	const Bool_t expert = kTRUE ; \r
	228	const Bool_t image = kTRUE ; \r
	229	TH1F * fhDigitsModule;\r
41d0241f	230	const char *acoModule[60]={"0_0","0_1","0_2","0_3","0_4","0_5","0_6","0_7","0_8","0_9",\r
dc7f1e9e	231	"1_0","1_1","1_2","1_3","1_4","1_5","1_6","1_7","1_8","1_9",\r
	232	"2_0","2_1","2_2","2_3","2_4","2_5","2_6","2_7","2_8","2_9",\r
	233	"3_0","3_1","3_2","3_3","3_4","3_5","3_6","3_7","3_8","3_9",\r
	234	"4_0","4_1","4_2","4_3","4_4","4_5","4_6","4_7","4_8","4_9",\r
	235	"5_0","5_1","5_2","5_3","5_4","5_5","5_6","5_7","5_8","5_9"};\r
	236	\r
	237	\r
	238	fhDigitsModule = new TH1F("ACORDEBitPatternfromDigits","Distribution of ACORDE from DIGITS;Modules;Counts",60,1,60);\r
	239	Add2DigitsList(fhDigitsModule,0,!expert,image);\r
	240	for (Int_t i=0;i<60;i++) fhDigitsModule->GetXaxis()->SetBinLabel(i+1,acoModule[i]); \r
92664bc8	241	//\r
3542faa1	242	ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line\r
dc7f1e9e	243	}\r
	244	\r
	245	//____________________________________________________________________________ \r
	246	\r
	247	void AliACORDEQADataMakerRec::InitRecPoints()\r
	248	{\r
	249	// create cluster histograms in RecPoint subdir\r
	250	// Not needed for ACORDE by now !!!\r
92664bc8	251	//\r
92664bc8	252	ClonePerTrigClass(AliQAv1::kRECPOINTS); // this should be the last line\r
dc7f1e9e	253	}\r
	254	\r
	255	//____________________________________________________________________________\r
	256	void AliACORDEQADataMakerRec::InitESDs()\r
	257	{\r
	258	//create ESDs histograms in ESDs subdir\r
	259	\r
	260	const Bool_t expert = kTRUE ; \r
	261	const Bool_t image = kTRUE ; \r
	262	\r
	263	TH1F * fhESDsSingle;\r
	264	TH1F * fhESDsMulti;\r
1b368017	265	TH1F * fhESDsMultiplicity;\r
41d0241f	266	const char *acoModule[60]={"0_0","0_1","0_2","0_3","0_4","0_5","0_6","0_7","0_8","0_9",\r
dc7f1e9e	267	"1_0","1_1","1_2","1_3","1_4","1_5","1_6","1_7","1_8","1_9",\r
	268	"2_0","2_1","2_2","2_3","2_4","2_5","2_6","2_7","2_8","2_9",\r
	269	"3_0","3_1","3_2","3_3","3_4","3_5","3_6","3_7","3_8","3_9",\r
	270	"4_0","4_1","4_2","4_3","4_4","4_5","4_6","4_7","4_8","4_9",\r
	271	"5_0","5_1","5_2","5_3","5_4","5_5","5_6","5_7","5_8","5_9"};\r
	272	\r
	273	\r
	274	fhESDsSingle = new TH1F("ACORDEBitPatternfromESDsSingle","Distribution of ACORDE fired modules from ESDs-Single;Modules;Counts",60,1,60);\r
	275	Add2ESDsList(fhESDsSingle,0,!expert,image);\r
	276	\r
	277	fhESDsMulti = new TH1F("ACORDEBitPatternfromESDsMulti","Distribution of ACORDE fired modules from ESDs-Multi;Modules;Counts",60,1,60);\r
	278	Add2ESDsList(fhESDsMulti,1,!expert,image);\r
1b368017	279	\r
	280	fhESDsMultiplicity = new TH1F("ACORDEMultiplicityfromESD","Number of fired ACORDE modules; No. of fired ACORDE modules;No. of events in ACORDE",60,-0.5,60);\r
	281	Add2ESDsList(fhESDsMultiplicity,2,!expert,image); \r
dc7f1e9e	282	for (Int_t i=0;i<60;i++)\r
	283	{\r
	284	fhESDsSingle->GetXaxis()->SetBinLabel(i+1,acoModule[i]);\r
	285	fhESDsMulti->GetXaxis()->SetBinLabel(i+1,acoModule[i]);\r
	286	}\r
92664bc8	287	//\r
92664bc8	288	ClonePerTrigClass(AliQAv1::kESDS); // this should be the last line\r
dc7f1e9e	289	}\r
	290	//____________________________________________________________________________\r
	291	void AliACORDEQADataMakerRec::MakeRaws(AliRawReader* rawReader)\r
	292	{\r
	293	//fills QA histos for RAW\r
92664bc8	294	IncEvCountCycleRaws();\r
92664bc8	295	IncEvCountTotalRaws();\r
dc7f1e9e	296	rawReader->Reset();\r
	297	AliACORDERawStream rawStream(rawReader);\r
	298	size_t contSingle=0;\r
	299	size_t contMulti=0;\r
	300	UInt_t dy[4];\r
	301	\r
	302	bool kroSingle[60],kroMulti[60];\r
	303	UInt_t tmpDy;\r
	304	\r
	305	for(Int_t m=0;m<60;m++) {kroSingle[m]=0;kroMulti[m]=0;}\r
	306	\r
	307	if(rawStream.Next())\r
	308	{\r
	309	dy[0]=rawStream.GetWord(0);\r
	310	dy[1]=rawStream.GetWord(1);\r
	311	dy[2]=rawStream.GetWord(2);\r
	312	dy[3]=rawStream.GetWord(3);\r
	313	tmpDy=dy[0];\r
10f9e41a	314	for(Int_t r=0;r<30;r++)\r
dc7f1e9e	315	{\r
	316	kroSingle[r] = tmpDy & 1;\r
	317	tmpDy>>=1;\r
	318	}\r
	319	tmpDy=dy[1];\r
10f9e41a	320	for(Int_t r=30;r<60;r++)\r
dc7f1e9e	321	{\r
	322	kroSingle[r] = tmpDy & 1;\r
	323	tmpDy>>=1;\r
	324	}\r
	325	tmpDy=dy[2];\r
10f9e41a	326	for(Int_t r=0;r<30;r++)\r
dc7f1e9e	327	{\r
	328	kroMulti[r] = tmpDy & 1;\r
	329	tmpDy>>=1;\r
	330	}\r
	331	tmpDy=dy[3];\r
10f9e41a	332	for(Int_t r=30;r<60;r++)\r
dc7f1e9e	333	{\r
	334	kroMulti[r] = tmpDy & 1;\r
	335	tmpDy>>=1;\r
	336	}\r
	337	contSingle=0;\r
92664bc8	338	contMulti=0;\r
10f9e41a	339	for(Int_t r=0;r<60;r++)\r
dc7f1e9e	340	{\r
92664bc8	341	if(kroSingle[r]==1)\r
	342	{\r
	343	FillRawsData(0,r);\r
	344	//FillRawsData(4,r);\r
	345	contSingle++;\r
	346	}\r
1b368017	347	if(kroMulti[r]==1)\r
1b368017	348	{\r
92664bc8	349	FillRawsData(2,r);\r
	350	//FillRawsData(6,r);\r
	351	contMulti++;\r
1b368017	352	}\r
92664bc8	353	\r
92664bc8	354	} \r
47f5f498	355	FillRawsData(1,contSingle); \r
47f5f498	356	FillRawsData(3,contMulti); \r
9f9f68fb	357	// FillRawsData(7,contMulti);\r
dc7f1e9e	358	}\r
	359	}\r
	360	//____________________________________________________________________________\r
	361	void AliACORDEQADataMakerRec::MakeDigits( TTree *digitsTree)\r
	362	{\r
	363	//fills QA histos for Digits\r
3542faa1	364	IncEvCountCycleDigits();\r
3542faa1	365	IncEvCountTotalDigits();\r
92664bc8	366	\r
dc7f1e9e	367	if (fDigitsArray) \r
	368	fDigitsArray->Clear() ; \r
	369	else \r
	370	fDigitsArray = new TClonesArray("AliACORDEdigit",1000);\r
	371	TBranch * branch = digitsTree->GetBranch("ACORDEdigit");\r
	372	if (!branch) {\r
	373	AliWarning("ACORDE branch in Digits Tree not found");\r
	374	} else {\r
	375	branch->SetAddress(&fDigitsArray);\r
	376	for(Int_t track = 0 ; track < branch->GetEntries() ; track++) {\r
	377	branch->GetEntry(track);\r
	378	for(Int_t idigit = 0 ; idigit < fDigitsArray->GetEntriesFast() ; idigit++) {\r
	379	AliACORDEdigit AcoDigit = (AliACORDEdigit) fDigitsArray->UncheckedAt(idigit);\r
	380	if (!AcoDigit) {\r
	381	AliError("The unchecked digit doesn't exist");\r
	382	continue ;\r
	383	}\r
92664bc8	384	FillDigitsData(0,AcoDigit->GetModule()-1);\r
dc7f1e9e	385	}\r
	386	}\r
	387	}\r
	388	}\r
	389	\r
	390	//____________________________________________________________________________\r
	391	void AliACORDEQADataMakerRec::MakeESDs(AliESDEvent * esd)\r
	392	{\r
	393	//fills QA histos for ESD\r
92664bc8	394	IncEvCountCycleESDs();\r
	395	IncEvCountTotalESDs();\r
	396	//\r
dc7f1e9e	397	AliESDACORDE * fESDACORDE= esd->GetACORDEData();\r
1b368017	398	Int_t acoMulti=0;\r
dc7f1e9e	399	for(int i=0;i<60;i++)\r
dc7f1e9e	400	{\r
92664bc8	401	if(fESDACORDE->GetHitChannel(i)) \r
	402	{\r
	403	FillESDsData(0,i+1);\r
	404	FillESDsData(1,i+1);\r
	405	acoMulti++;\r
dc7f1e9e	406	}\r
92664bc8	407	} FillESDsData(2,acoMulti);\r
dc7f1e9e	408	\r
dc7f1e9e	409	}\r