[u/mrichter/AliRoot.git] / PMD / AliPMDQADataMakerRec.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.
  B.K. Nandi
*/

// --- ROOT system ---
#include <TClonesArray.h>
#include <TFile.h> 
#include <TH1F.h> 
#include <TH1I.h> 
#include <TH2F.h> 

// --- Standard library ---

// --- AliRoot header files ---

#include "AliESDEvent.h"
#include "AliLog.h"
#include "AliPMDQADataMakerRec.h"
#include "AliQAChecker.h"
#include "AliPMDrecpoint1.h" 
#include "AliPMDRawStream.h"
#include "AliPMDddldata.h"
#include "AliPMDUtility.h"
#include "AliESDPmdTrack.h"
//#include "AliPMDRecoParam.h"

ClassImp(AliPMDQADataMakerRec)
           
//____________________________________________________________________________ 
  AliPMDQADataMakerRec::AliPMDQADataMakerRec() : 
  AliQADataMakerRec(AliQA::GetDetName(AliQA::kPMD), "PMD Quality Assurance Data Maker")
{
  // ctor
}

//____________________________________________________________________________ 
AliPMDQADataMakerRec::AliPMDQADataMakerRec(const AliPMDQADataMakerRec& qadm) :
  AliQADataMakerRec()
{
  //copy ctor 
  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle()); 
}

//__________________________________________________________________
AliPMDQADataMakerRec& AliPMDQADataMakerRec::operator = (const AliPMDQADataMakerRec& qadm )
{
  // Equal operator.
  this->~AliPMDQADataMakerRec();
  new(this) AliPMDQADataMakerRec(qadm);
  return *this;
}
 

//____________________________________________________________________________ 
void AliPMDQADataMakerRec::InitRaws()
{
  // create Raws histograms in Raws subdir


    // Preshower plane

    TH1I * h0 = new TH1I("hPreEdepM0","ADC Distribution PRE - Module 0", 100, 0, 2000);
    h0->Sumw2();
    Add2RawsList(h0, 0); 

    TH1I * h1 = new TH1I("hPreEdepM1","ADC Distribution PRE - Module 1", 100, 0, 2000);
    h1->Sumw2();
    Add2RawsList(h1, 1); 

    TH1I * h2 = new TH1I("hPreEdepM2","ADC Distribution PRE - Module 2", 100, 0, 2000);
    h2->Sumw2();
    Add2RawsList(h2, 2); 
    
    TH1I * h3 = new TH1I("hPreEdepM3","ADC Distribution PRE - Module 3", 100, 0, 2000);
    h3->Sumw2(); 
    Add2RawsList(h3, 3); 

    TH1I * h4 = new TH1I("hPreEdepM4","ADC Distribution PRE - Module 4", 100, 0, 2000);
    h4->Sumw2(); 
    Add2RawsList(h4, 4); 

    TH1I * h5 = new TH1I("hPreEdepM5","ADC Distribution PRE - Module 5", 100, 0, 2000);
    h5->Sumw2(); 
    Add2RawsList(h5, 5); 

    TH1I * h6 = new TH1I("hPreEdepM6","ADC Distribution PRE - Module 6", 100, 0, 2000);
    h6->Sumw2(); 
    Add2RawsList(h6, 6); 

    TH1I * h7 = new TH1I("hPreEdepM7","ADC Distribution PRE - Module 7", 100, 0, 2000);
    h7->Sumw2(); 
    Add2RawsList(h7, 7); 

    TH1I * h8 = new TH1I("hPreEdepM8","ADC Distribution PRE - Module 8", 100, 0, 2000);
    h8->Sumw2(); 
    Add2RawsList(h8, 8); 

    TH1I * h9 = new TH1I("hPreEdepM9","ADC Distribution PRE - Module 9", 100, 0, 2000);
    h9->Sumw2(); 
    Add2RawsList(h9, 9); 

    TH1I * h10 = new TH1I("hPreEdepM10","ADC Distribution PRE - Module 10", 100, 0, 2000);
    h10->Sumw2(); 
    Add2RawsList(h10, 10); 

    TH1I * h11 = new TH1I("hPreEdepM11","ADC Distribution PRE - Module 11", 100, 0, 2000);
    h11->Sumw2(); 
    Add2RawsList(h11, 11); 

    TH1I * h12 = new TH1I("hPreEdepM12","ADC Distribution PRE - Module 12", 100, 0, 2000);
    h12->Sumw2(); 
    Add2RawsList(h12, 12); 

    TH1I * h13 = new TH1I("hPreEdepM13","ADC Distribution PRE - Module 13", 100, 0, 2000);
    h13->Sumw2(); 
    Add2RawsList(h13, 13); 

    TH1I * h14 = new TH1I("hPreEdepM14","ADC Distribution PRE - Module 14", 100, 0, 2000);
    h14->Sumw2(); 
    Add2RawsList(h14, 14); 

    TH1I * h15 = new TH1I("hPreEdepM15","ADC Distribution PRE - Module 15", 100, 0, 2000);
    h15->Sumw2(); 
    Add2RawsList(h15, 15); 

    TH1I * h16 = new TH1I("hPreEdepM16","ADC Distribution PRE - Module 16", 100, 0, 2000);
    h16->Sumw2(); 
    Add2RawsList(h16, 16); 

    TH1I * h17 = new TH1I("hPreEdepM17","ADC Distribution PRE - Module 17", 100, 0, 2000);
    h17->Sumw2(); 
    Add2RawsList(h17, 17); 

    TH1I * h18 = new TH1I("hPreEdepM18","ADC Distribution PRE - Module 18", 100, 0, 2000);
    h18->Sumw2(); 
    Add2RawsList(h18, 18); 

    TH1I * h19 = new TH1I("hPreEdepM19","ADC Distribution PRE - Module 19", 100, 0, 2000);
    h19->Sumw2(); 
    Add2RawsList(h19, 19); 

    TH1I * h20 = new TH1I("hPreEdepM20","ADC Distribution PRE - Module 20", 100, 0, 2000);
    h20->Sumw2(); 
    Add2RawsList(h20, 20);
 
    TH1I * h21 = new TH1I("hPreEdepM21","ADC Distribution PRE - Module 21", 100, 0, 2000);
    h21->Sumw2(); 
    Add2RawsList(h21, 21); 

    TH1I * h22 = new TH1I("hPreEdepM22","ADC Distribution PRE - Module 22", 100, 0, 2000);
    h22->Sumw2(); 
    Add2RawsList(h22, 22); 

    TH1I * h23 = new TH1I("hPreEdepM23","ADC Distribution PRE - Module 23", 100, 0, 2000);
    h23->Sumw2(); 
    Add2RawsList(h23, 23);

    // CPV histos

    TH1I * h24 = new TH1I("hCpvEdepM24","ADC Distribution CPV - Module 24", 100, 0, 2000);
    h24->Sumw2();
    Add2RawsList(h24, 24); 

    TH1I * h25 = new TH1I("hCpvEdepM25","ADC Distribution CPV - Module 25", 100, 0, 2000);
    h25->Sumw2(); 
    Add2RawsList(h25, 25); 

    TH1I * h26 = new TH1I("hCpvEdepM26","ADC Distribution CPV - Module 26", 100, 0, 2000);
    h26->Sumw2(); 
    Add2RawsList(h26, 26); 

    TH1I * h27 = new TH1I("hCpvEdepM27","ADC Distribution CPV - Module 27", 100, 0, 2000);
    h27->Sumw2(); 
    Add2RawsList(h27, 27); 

    TH1I * h28 = new TH1I("hCpvEdepM28","ADC Distribution CPV - Module 28", 100, 0, 2000);
    h28->Sumw2(); 
    Add2RawsList(h28, 28); 

    TH1I * h29 = new TH1I("hCpvEdepM29","ADC Distribution CPV - Module 29", 100, 0, 2000);
    h29->Sumw2(); 
    Add2RawsList(h29, 29); 

    TH1I * h30 = new TH1I("hCpvEdepM30","ADC Distribution CPV - Module 30", 100, 0, 2000);
    h30->Sumw2(); 
    Add2RawsList(h30, 30); 

    TH1I * h31 = new TH1I("hCpvEdepM31","ADC Distribution CPV - Module 31", 100, 0, 2000);
    h31->Sumw2(); 
    Add2RawsList(h31, 31); 

    TH1I * h32 = new TH1I("hCpvEdepM32","ADC Distribution CPV - Module 32", 100, 0, 2000);
    h32->Sumw2();
    Add2RawsList(h32, 32); 

    TH1I * h33 = new TH1I("hCpvEdepM33","ADC Distribution CPV - Module 33", 100, 0, 2000);
    h33->Sumw2();
    Add2RawsList(h33, 33); 

    TH1I * h34 = new TH1I("hCpvEdepM34","ADC Distribution CPV - Module 34", 100, 0, 2000);
    h34->Sumw2(); 
    Add2RawsList(h34, 34); 

    TH1I * h35 = new TH1I("hCpvEdepM35","ADC Distribution CPV - Module 35", 100, 0, 2000);
    h35->Sumw2(); 
    Add2RawsList(h35, 35); 

    TH1I * h36 = new TH1I("hCpvEdepM36","ADC Distribution CPV - Module 36", 100, 0, 2000);
    h36->Sumw2(); 
    Add2RawsList(h36, 36); 

    TH1I * h37 = new TH1I("hCpvEdepM37","ADC Distribution CPV - Module 37", 100, 0, 2000);
    h37->Sumw2(); 
    Add2RawsList(h37, 37); 

    TH1I * h38 = new TH1I("hCpvEdepM38","ADC Distribution CPV - Module 38", 100, 0, 2000);
    h38->Sumw2(); 
    Add2RawsList(h38, 38); 

    TH1I * h39 = new TH1I("hCpvEdepM39","ADC Distribution CPV - Module 39", 100, 0, 2000);
    h39->Sumw2(); 
    Add2RawsList(h39, 39); 

    TH1I * h40 = new TH1I("hCpvEdepM40","ADC Distribution CPV - Module 40", 100, 0, 2000);
    h40->Sumw2(); 
    Add2RawsList(h40, 40); 

    TH1I * h41 = new TH1I("hCpvEdepM41","ADC Distribution CPV - Module 41", 100, 0, 2000);
    h41->Sumw2();
    Add2RawsList(h41, 41); 

    TH1I * h42 = new TH1I("hCpvEdepM42","ADC Distribution CPV - Module 42", 100, 0, 2000);
    h42->Sumw2(); 
    Add2RawsList(h42, 42); 

    TH1I * h43 = new TH1I("hCpvEdepM43","ADC Distribution CPV - Module 43", 100, 0, 2000);
    h43->Sumw2(); 
    Add2RawsList(h43, 43); 

    TH1I * h44 = new TH1I("hCpvEdepM44","ADC Distribution CPV - Module 44", 100, 0, 2000);
    h44->Sumw2(); 
    Add2RawsList(h44, 44); 

    TH1I * h45 = new TH1I("hCpvEdepM45","ADC Distribution CPV - Module 45", 100, 0, 2000);
    h45->Sumw2();
    Add2RawsList(h45, 45); 

    TH1I * h46 = new TH1I("hCpvEdepM46","ADC Distribution CPV - Module 46", 100, 0, 2000);
    h46->Sumw2(); 
    Add2RawsList(h46, 46); 

    TH1I * h47 = new TH1I("hCpvEdepM47","ADC Distribution CPV - Module 47", 100, 0, 2000);
    h47->Sumw2(); 
    Add2RawsList(h47, 47); 


    // Y vs. X for PRE and CPV planes

    TH2F * h48 = new TH2F("hPreXY","PRE plane",200,-100.,100.,200,-100.,100.);
    Add2RawsList(h48, 48);
    TH2F * h49 = new TH2F("hCpvXY","CPV plane",200,-100.,100.,200,-100.,100.);
    Add2RawsList(h49, 49);

}
//____________________________________________________________________________ 
void AliPMDQADataMakerRec::InitRecPoints()
{
  // create Reconstructed Points histograms in RecPoints subdir
  
  /*
    TH2F * h0 = new TH2F("hPreXY","RecPoints Y vs X PRE(PMD)", 100,-100.,100.,100,-100.,100.);
    Add2RecPointsList(h0,0);
    
    TH2F * h1 = new TH2F("hCpvXY","RecPoints Y vs X CPV(PMD)", 100,-100.,100.,100,-100.,100.);
    Add2RecPointsList(h1,1);
  */

    //  Ncell distribution in a cluster
  
    // PRE plane

  TH1F * h0 = new TH1F("hPreDdl0Ncell","PRE: Ddl0 Ncell in a cluster",50,0.,50.);
  h0->Sumw2();
  Add2RecPointsList(h0, 0);


  TH1F * h1 = new TH1F("hPreDdl1Ncell","PRE: Ddl1 Ncell in a cluster",50,0.,50.);
  h1->Sumw2();
  Add2RecPointsList(h1, 1);


  TH1F * h2 = new TH1F("hPreDdl2Ncell","PRE: Ddl2 Ncell in a cluster",50,0.,50.);
  h2->Sumw2();
  Add2RecPointsList(h2, 2);


  TH1F * h3 = new TH1F("hPreDdl3Ncell","PRE: Ddl3 Ncell in a cluster",50,0.,50.);
  h3->Sumw2();
  Add2RecPointsList(h3, 3);

  // CPV plane

  TH1F * h4 = new TH1F("hCpvDdl4Ncell","CPV: Ddl4 Ncell in a cluster",50,0.,50.);
  h4->Sumw2();
  Add2RecPointsList(h4, 4);

  TH1F * h5 = new TH1F("hCpvDdl5Ncell","CPV: Ddl5 Ncell in a cluster",50,0.,50.);
  h5->Sumw2();
  Add2RecPointsList(h5, 5);

  // Correlation plot

  TH2I *h6 = new TH2I("hPre10","Cluster - DDL1 vs DDL0", 100,0,200,100,0,200);
  Add2RecPointsList(h6,6);

  TH2I *h7 = new TH2I("hPre32","Cluster - DDL3 vs DDL2", 100,0,200,100,0,200);
  Add2RecPointsList(h7,7);

  TH2I *h8 = new TH2I("hCpv54","Cluster - DDL5 vs DDL4", 100,0,200,100,0,200);
  Add2RecPointsList(h8,8);


}

//____________________________________________________________________________ 

void AliPMDQADataMakerRec::InitESDs()
{
  //Create histograms to controll ESD

  TH1F *h0 = new TH1F("hPreClADC","Cluster ADC of PRE plane",500,0.,10000.);
  h0->Sumw2();
  Add2ESDsList(h0, 0)  ;

  TH1F *h1 = new TH1F("hCpvClADC","Cluster ADC of CPV plane",500,0.,10000.);
  h1->Sumw2();
  Add2ESDsList(h1, 1)  ;

  TH2I *h2 = new TH2I("hPmdClMult","Cluster Multiplicity: PRE vs. CPVplane",100,0,1000,100,0,1000);
  h2->Sumw2();
  Add2ESDsList(h2, 2)  ;

/*
  TH1I * h3 = new TH1I("hCpvClMult","Cluster Multiplicity of CPV plane",100,0.,1000.);
  h3->Sumw2();
  Add2ESDsList(h3, 3)  ;
*/

}

//____________________________________________________________________________
void AliPMDQADataMakerRec::MakeRaws(AliRawReader* rawReader)
{
    //Fill prepared histograms with Raw digit properties

    TObjArray *pmdddlcont = 0x0;
    pmdddlcont = new TObjArray();
    AliPMDRawStream stream(rawReader);
    
    AliPMDddldata *pmdddl = 0x0;

    Int_t   iddl = -1;
    Int_t   xpad = -1;
    Int_t   ypad = -1;
    Float_t xx, yy;

    AliPMDUtility cc;

    while ((iddl = stream.DdlData(pmdddlcont)) >=0)
    {
	Int_t ientries = pmdddlcont->GetEntries();
	//printf(" ======= DDLNO = %d ientries = %d \n", iddl, ientries);
	
	for (Int_t ient = 0; ient < ientries; ient++)
	  {
	      //AliPMDddldata *pmdddl = (AliPMDddldata*)pmdddlcont->UncheckedAt(ient);
	    pmdddl = (AliPMDddldata*)pmdddlcont->UncheckedAt(ient);

	    Int_t det = pmdddl->GetDetector();
	    Int_t smn = pmdddl->GetSMN();
	    Int_t mcm = pmdddl->GetMCM();
	    //Int_t chno = pmdddl->GetChannel();
	    Int_t row = pmdddl->GetRow();
	    Int_t col = pmdddl->GetColumn();
	    Int_t sig = pmdddl->GetSignal();
	    
	    if (mcm == 0 || smn == -1) continue;

	    if (det == 0)
	    {
		GetRawsData(smn)->Fill(sig);
		if(smn < 12)
		{
		    xpad = col;
		    ypad = row;
		}
		else if(smn >= 12 && smn < 24)
		{
		    xpad = row;
		    ypad = col;
		}
		cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
		GetRawsData(48)->Fill(xx,yy);
	    }
	    if (det == 1)
	    {
		GetRawsData(24+smn)->Fill(sig);
		if(smn < 12)
		{
		    xpad = col;
		    ypad = row;
		}
		else if(smn >= 12 && smn < 24)
		{
		    xpad = row;
		    ypad = col;
		}

		cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
		GetRawsData(49)->Fill(xx,yy);

	    }

	  }

	pmdddlcont->Delete();
    }

    delete pmdddlcont;
    pmdddlcont = 0x0;

}
//____________________________________________________________________________
void AliPMDQADataMakerRec::MakeRecPoints(TTree * clustersTree)
{
    // makes data from RecPoints

    Int_t multDdl0 = 0, multDdl1 = 0, multDdl2 = 0;
    Int_t multDdl3 = 0, multDdl4 = 0, multDdl5 = 0;

    AliPMDrecpoint1 * recpoint; 

    TClonesArray * recpoints = 0x0;
    recpoints = new TClonesArray("AliPMDrecpoint1", 1000) ; 
    
    TBranch * branch = clustersTree->GetBranch("PMDRecpoint") ;
    branch->SetAddress(&recpoints) ;

    if ( ! branch )
    {
	AliWarning("PMD branch in Recpoints Tree not found") ; 
    }
    else
    {
	for (Int_t imod = 0; imod < branch->GetEntries(); imod++)
	{
	    branch->GetEntry(imod) ;

	    TIter next(recpoints) ; 

	    while ( (recpoint = dynamic_cast<AliPMDrecpoint1 *>(next())) )
	      {
		//Float_t xpos = recpoint->GetClusX();
		//Float_t ypos = recpoint->GetClusY();
		//Int_t smn = recpoint->GetSMNumber();
		
		  if(recpoint->GetDetector() == 0)
		  {
		    if(recpoint->GetSMNumber() >= 0 && recpoint->GetSMNumber() < 6)
		      {
			GetRecPointsData(0)->Fill(recpoint->GetClusCells());
			multDdl0++;
		      }
		    if(recpoint->GetSMNumber() >= 6 && recpoint->GetSMNumber() < 12)
		      {
			GetRecPointsData(1)->Fill(recpoint->GetClusCells());
			multDdl1++;
		      }
		    if(recpoint->GetSMNumber() >= 12 && recpoint->GetSMNumber() < 18)
		      {
			GetRecPointsData(2)->Fill(recpoint->GetClusCells());
			multDdl2++;
		      }
		    if(recpoint->GetSMNumber() >= 18 && recpoint->GetSMNumber() < 24)
		      {
			GetRecPointsData(3)->Fill(recpoint->GetClusCells());
			multDdl3++;
		      }
		  }

		if(recpoint->GetDetector() == 1)
		  {
		    if(recpoint->GetSMNumber() >= 0 && recpoint->GetSMNumber() < 6 || 
		       recpoint->GetSMNumber() >= 18 && recpoint->GetSMNumber() < 24)
		      {
			GetRecPointsData(4)->Fill(recpoint->GetClusCells());
			multDdl4++;
		      }
		    if(recpoint->GetSMNumber() >= 6 && recpoint->GetSMNumber() < 18 )
		      {
			GetRecPointsData(5)->Fill(recpoint->GetClusCells());
			multDdl5++;
		      }
		  }
	      } 
	}
    }
    
    GetRecPointsData(6)->Fill(multDdl0,multDdl1);
    GetRecPointsData(7)->Fill(multDdl2,multDdl3);
    GetRecPointsData(8)->Fill(multDdl4,multDdl5);

    delete recpoints;

}

//____________________________________________________________________________

void AliPMDQADataMakerRec::MakeESDs(AliESDEvent * esd)
{
  // make QA data from ESDs

    Int_t premul = 0, cpvmul = 0;

  for (Int_t icl = 0; icl < esd->GetNumberOfPmdTracks(); icl++)
    {
      AliESDPmdTrack *pmdtr = esd->GetPmdTrack(icl);
      
      //Int_t   det   = pmdtr->GetDetector(); 
      //Float_t clsX  = pmdtr->GetClusterX();
      //Float_t clsY  = pmdtr->GetClusterY();
      //Float_t clsZ  = pmdtr->GetClusterZ();
      //Float_t ncell = pmdtr->GetClusterCells();
      Float_t adc   = pmdtr->GetClusterADC();
      //Float_t pid   = pmdtr->GetClusterPID();

      if (pmdtr->GetDetector() == 0)
	{
	  GetESDsData(0)->Fill(adc);
	  premul++;
	}
      if (pmdtr->GetDetector() == 1)
	{
	  GetESDsData(1)->Fill(adc) ;
	  cpvmul++;
	}
    }
  
  GetESDsData(2)->Fill(cpvmul,premul) ;
  //GetESDsData(3)->Fill(cpvmul) ;  
}

//____________________________________________________________________________ 

void AliPMDQADataMakerRec::StartOfDetectorCycle()
{
  //Detector specific actions at start of cycle
  
}
//____________________________________________________________________________ 
void AliPMDQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
{
  //Detector specific actions at end of cycle
  // do the QA checking
  AliQAChecker::Instance()->Run(AliQA::kPMD, task, list) ;  
}
Commit	Line	Data
a9523c3e	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	Produces the data needed to calculate the quality assurance.
	19	All data must be mergeable objects.
	20	B.K. Nandi
	21	*/
	22
	23	// --- ROOT system ---
	24	#include <TClonesArray.h>
	25	#include <TFile.h>
	26	#include <TH1F.h>
	27	#include <TH1I.h>
	28	#include <TH2F.h>
	29
	30	// --- Standard library ---
	31
	32	// --- AliRoot header files ---
	33
	34	#include "AliESDEvent.h"
	35	#include "AliLog.h"
	36	#include "AliPMDQADataMakerRec.h"
	37	#include "AliQAChecker.h"
	38	#include "AliPMDrecpoint1.h"
	39	#include "AliPMDRawStream.h"
	40	#include "AliPMDddldata.h"
	41	#include "AliPMDUtility.h"
	42	#include "AliESDPmdTrack.h"
	43	//#include "AliPMDRecoParam.h"
	44
	45	ClassImp(AliPMDQADataMakerRec)
	46
	47	//____________________________________________________________________________
	48	AliPMDQADataMakerRec::AliPMDQADataMakerRec() :
	49	AliQADataMakerRec(AliQA::GetDetName(AliQA::kPMD), "PMD Quality Assurance Data Maker")
	50	{
	51	// ctor
	52	}
	53
	54	//____________________________________________________________________________
	55	AliPMDQADataMakerRec::AliPMDQADataMakerRec(const AliPMDQADataMakerRec& qadm) :
	56	AliQADataMakerRec()
	57	{
	58	//copy ctor
	59	SetName((const char*)qadm.GetName()) ;
	60	SetTitle((const char*)qadm.GetTitle());
	61	}
	62
	63	//__________________________________________________________________
	64	AliPMDQADataMakerRec& AliPMDQADataMakerRec::operator = (const AliPMDQADataMakerRec& qadm )
65	{
66	// Equal operator.
67	this->~AliPMDQADataMakerRec();
68	new(this) AliPMDQADataMakerRec(qadm);
69	return *this;
70	}
71
72
73	//____________________________________________________________________________
74	void AliPMDQADataMakerRec::InitRaws()
75	{
76	// create Raws histograms in Raws subdir
a9523c3e	77
a9523c3e	78
2c1131dd	79	// Preshower plane
	80
	81	TH1I * h0 = new TH1I("hPreEdepM0","ADC Distribution PRE - Module 0", 100, 0, 2000);
	82	h0->Sumw2();
	83	Add2RawsList(h0, 0);
	84
	85	TH1I * h1 = new TH1I("hPreEdepM1","ADC Distribution PRE - Module 1", 100, 0, 2000);
	86	h1->Sumw2();
	87	Add2RawsList(h1, 1);
a9523c3e	88
2c1131dd	89	TH1I * h2 = new TH1I("hPreEdepM2","ADC Distribution PRE - Module 2", 100, 0, 2000);
	90	h2->Sumw2();
	91	Add2RawsList(h2, 2);
	92
	93	TH1I * h3 = new TH1I("hPreEdepM3","ADC Distribution PRE - Module 3", 100, 0, 2000);
	94	h3->Sumw2();
	95	Add2RawsList(h3, 3);
a9523c3e	96
2c1131dd	97	TH1I * h4 = new TH1I("hPreEdepM4","ADC Distribution PRE - Module 4", 100, 0, 2000);
	98	h4->Sumw2();
	99	Add2RawsList(h4, 4);
a9523c3e	100
2c1131dd	101	TH1I * h5 = new TH1I("hPreEdepM5","ADC Distribution PRE - Module 5", 100, 0, 2000);
	102	h5->Sumw2();
	103	Add2RawsList(h5, 5);
	104
	105	TH1I * h6 = new TH1I("hPreEdepM6","ADC Distribution PRE - Module 6", 100, 0, 2000);
	106	h6->Sumw2();
	107	Add2RawsList(h6, 6);
	108
	109	TH1I * h7 = new TH1I("hPreEdepM7","ADC Distribution PRE - Module 7", 100, 0, 2000);
	110	h7->Sumw2();
	111	Add2RawsList(h7, 7);
	112
	113	TH1I * h8 = new TH1I("hPreEdepM8","ADC Distribution PRE - Module 8", 100, 0, 2000);
	114	h8->Sumw2();
	115	Add2RawsList(h8, 8);
	116
	117	TH1I * h9 = new TH1I("hPreEdepM9","ADC Distribution PRE - Module 9", 100, 0, 2000);
	118	h9->Sumw2();
	119	Add2RawsList(h9, 9);
	120
	121	TH1I * h10 = new TH1I("hPreEdepM10","ADC Distribution PRE - Module 10", 100, 0, 2000);
	122	h10->Sumw2();
	123	Add2RawsList(h10, 10);
	124
	125	TH1I * h11 = new TH1I("hPreEdepM11","ADC Distribution PRE - Module 11", 100, 0, 2000);
	126	h11->Sumw2();
	127	Add2RawsList(h11, 11);
	128
	129	TH1I * h12 = new TH1I("hPreEdepM12","ADC Distribution PRE - Module 12", 100, 0, 2000);
	130	h12->Sumw2();
	131	Add2RawsList(h12, 12);
	132
	133	TH1I * h13 = new TH1I("hPreEdepM13","ADC Distribution PRE - Module 13", 100, 0, 2000);
	134	h13->Sumw2();
	135	Add2RawsList(h13, 13);
	136
	137	TH1I * h14 = new TH1I("hPreEdepM14","ADC Distribution PRE - Module 14", 100, 0, 2000);
	138	h14->Sumw2();
	139	Add2RawsList(h14, 14);
	140
	141	TH1I * h15 = new TH1I("hPreEdepM15","ADC Distribution PRE - Module 15", 100, 0, 2000);
	142	h15->Sumw2();
	143	Add2RawsList(h15, 15);
	144
	145	TH1I * h16 = new TH1I("hPreEdepM16","ADC Distribution PRE - Module 16", 100, 0, 2000);
	146	h16->Sumw2();
	147	Add2RawsList(h16, 16);
	148
	149	TH1I * h17 = new TH1I("hPreEdepM17","ADC Distribution PRE - Module 17", 100, 0, 2000);
	150	h17->Sumw2();
	151	Add2RawsList(h17, 17);
	152
	153	TH1I * h18 = new TH1I("hPreEdepM18","ADC Distribution PRE - Module 18", 100, 0, 2000);
	154	h18->Sumw2();
	155	Add2RawsList(h18, 18);
	156
	157	TH1I * h19 = new TH1I("hPreEdepM19","ADC Distribution PRE - Module 19", 100, 0, 2000);
	158	h19->Sumw2();
	159	Add2RawsList(h19, 19);
	160
	161	TH1I * h20 = new TH1I("hPreEdepM20","ADC Distribution PRE - Module 20", 100, 0, 2000);
	162	h20->Sumw2();
	163	Add2RawsList(h20, 20);
	164
165	TH1I * h21 = new TH1I("hPreEdepM21","ADC Distribution PRE - Module 21", 100, 0, 2000);
166	h21->Sumw2();
167	Add2RawsList(h21, 21);
168
169	TH1I * h22 = new TH1I("hPreEdepM22","ADC Distribution PRE - Module 22", 100, 0, 2000);
170	h22->Sumw2();
171	Add2RawsList(h22, 22);
172
173	TH1I * h23 = new TH1I("hPreEdepM23","ADC Distribution PRE - Module 23", 100, 0, 2000);
174	h23->Sumw2();
175	Add2RawsList(h23, 23);
176
177	// CPV histos
178
179	TH1I * h24 = new TH1I("hCpvEdepM24","ADC Distribution CPV - Module 24", 100, 0, 2000);
180	h24->Sumw2();
181	Add2RawsList(h24, 24);
182
183	TH1I * h25 = new TH1I("hCpvEdepM25","ADC Distribution CPV - Module 25", 100, 0, 2000);
184	h25->Sumw2();
185	Add2RawsList(h25, 25);
186
187	TH1I * h26 = new TH1I("hCpvEdepM26","ADC Distribution CPV - Module 26", 100, 0, 2000);
188	h26->Sumw2();
189	Add2RawsList(h26, 26);
190
191	TH1I * h27 = new TH1I("hCpvEdepM27","ADC Distribution CPV - Module 27", 100, 0, 2000);
192	h27->Sumw2();
193	Add2RawsList(h27, 27);
194
195	TH1I * h28 = new TH1I("hCpvEdepM28","ADC Distribution CPV - Module 28", 100, 0, 2000);
196	h28->Sumw2();
197	Add2RawsList(h28, 28);
198
199	TH1I * h29 = new TH1I("hCpvEdepM29","ADC Distribution CPV - Module 29", 100, 0, 2000);
200	h29->Sumw2();
201	Add2RawsList(h29, 29);
202
203	TH1I * h30 = new TH1I("hCpvEdepM30","ADC Distribution CPV - Module 30", 100, 0, 2000);
204	h30->Sumw2();
205	Add2RawsList(h30, 30);
206
207	TH1I * h31 = new TH1I("hCpvEdepM31","ADC Distribution CPV - Module 31", 100, 0, 2000);
208	h31->Sumw2();
209	Add2RawsList(h31, 31);
210
211	TH1I * h32 = new TH1I("hCpvEdepM32","ADC Distribution CPV - Module 32", 100, 0, 2000);
212	h32->Sumw2();
213	Add2RawsList(h32, 32);
214
215	TH1I * h33 = new TH1I("hCpvEdepM33","ADC Distribution CPV - Module 33", 100, 0, 2000);
216	h33->Sumw2();
217	Add2RawsList(h33, 33);
218
219	TH1I * h34 = new TH1I("hCpvEdepM34","ADC Distribution CPV - Module 34", 100, 0, 2000);
220	h34->Sumw2();
221	Add2RawsList(h34, 34);
222
223	TH1I * h35 = new TH1I("hCpvEdepM35","ADC Distribution CPV - Module 35", 100, 0, 2000);
224	h35->Sumw2();
225	Add2RawsList(h35, 35);
226
227	TH1I * h36 = new TH1I("hCpvEdepM36","ADC Distribution CPV - Module 36", 100, 0, 2000);
228	h36->Sumw2();
229	Add2RawsList(h36, 36);
230
231	TH1I * h37 = new TH1I("hCpvEdepM37","ADC Distribution CPV - Module 37", 100, 0, 2000);
232	h37->Sumw2();
233	Add2RawsList(h37, 37);
234
235	TH1I * h38 = new TH1I("hCpvEdepM38","ADC Distribution CPV - Module 38", 100, 0, 2000);
236	h38->Sumw2();
237	Add2RawsList(h38, 38);
238
239	TH1I * h39 = new TH1I("hCpvEdepM39","ADC Distribution CPV - Module 39", 100, 0, 2000);
240	h39->Sumw2();
241	Add2RawsList(h39, 39);
242
243	TH1I * h40 = new TH1I("hCpvEdepM40","ADC Distribution CPV - Module 40", 100, 0, 2000);
244	h40->Sumw2();
245	Add2RawsList(h40, 40);
246
247	TH1I * h41 = new TH1I("hCpvEdepM41","ADC Distribution CPV - Module 41", 100, 0, 2000);
248	h41->Sumw2();
249	Add2RawsList(h41, 41);
250
251	TH1I * h42 = new TH1I("hCpvEdepM42","ADC Distribution CPV - Module 42", 100, 0, 2000);
252	h42->Sumw2();
253	Add2RawsList(h42, 42);
254
255	TH1I * h43 = new TH1I("hCpvEdepM43","ADC Distribution CPV - Module 43", 100, 0, 2000);
256	h43->Sumw2();
257	Add2RawsList(h43, 43);
258
259	TH1I * h44 = new TH1I("hCpvEdepM44","ADC Distribution CPV - Module 44", 100, 0, 2000);
260	h44->Sumw2();
261	Add2RawsList(h44, 44);
262
263	TH1I * h45 = new TH1I("hCpvEdepM45","ADC Distribution CPV - Module 45", 100, 0, 2000);
264	h45->Sumw2();
265	Add2RawsList(h45, 45);
266
267	TH1I * h46 = new TH1I("hCpvEdepM46","ADC Distribution CPV - Module 46", 100, 0, 2000);
268	h46->Sumw2();
269	Add2RawsList(h46, 46);
270
271	TH1I * h47 = new TH1I("hCpvEdepM47","ADC Distribution CPV - Module 47", 100, 0, 2000);
272	h47->Sumw2();
273	Add2RawsList(h47, 47);
274
275
276	// Y vs. X for PRE and CPV planes
277
278	TH2F * h48 = new TH2F("hPreXY","PRE plane",200,-100.,100.,200,-100.,100.);
279	Add2RawsList(h48, 48);
280	TH2F * h49 = new TH2F("hCpvXY","CPV plane",200,-100.,100.,200,-100.,100.);
281	Add2RawsList(h49, 49);
a9523c3e	282
a9523c3e	283	}
	284	//____________________________________________________________________________
	285	void AliPMDQADataMakerRec::InitRecPoints()
	286	{
	287	// create Reconstructed Points histograms in RecPoints subdir
2c1131dd	288
2c1131dd	289	/*
a9523c3e	290	TH2F * h0 = new TH2F("hPreXY","RecPoints Y vs X PRE(PMD)", 100,-100.,100.,100,-100.,100.);
2c1131dd	291	Add2RecPointsList(h0,0);
2c1131dd	292
a9523c3e	293	TH2F * h1 = new TH2F("hCpvXY","RecPoints Y vs X CPV(PMD)", 100,-100.,100.,100,-100.,100.);
2c1131dd	294	Add2RecPointsList(h1,1);
2c1131dd	295	*/
a9523c3e	296
2c1131dd	297	// Ncell distribution in a cluster
2c1131dd	298
b5144d1f	299	// PRE plane
b5144d1f	300
2c1131dd	301	TH1F * h0 = new TH1F("hPreDdl0Ncell","PRE: Ddl0 Ncell in a cluster",50,0.,50.);
	302	h0->Sumw2();
	303	Add2RecPointsList(h0, 0);
a9523c3e	304
a9523c3e	305
2c1131dd	306	TH1F * h1 = new TH1F("hPreDdl1Ncell","PRE: Ddl1 Ncell in a cluster",50,0.,50.);
	307	h1->Sumw2();
	308	Add2RecPointsList(h1, 1);
	309
	310
	311	TH1F * h2 = new TH1F("hPreDdl2Ncell","PRE: Ddl2 Ncell in a cluster",50,0.,50.);
	312	h2->Sumw2();
	313	Add2RecPointsList(h2, 2);
	314
	315
	316	TH1F * h3 = new TH1F("hPreDdl3Ncell","PRE: Ddl3 Ncell in a cluster",50,0.,50.);
	317	h3->Sumw2();
	318	Add2RecPointsList(h3, 3);
	319
b5144d1f	320	// CPV plane
	321
	322	TH1F * h4 = new TH1F("hCpvDdl4Ncell","CPV: Ddl4 Ncell in a cluster",50,0.,50.);
2c1131dd	323	h4->Sumw2();
	324	Add2RecPointsList(h4, 4);
	325
b5144d1f	326	TH1F * h5 = new TH1F("hCpvDdl5Ncell","CPV: Ddl5 Ncell in a cluster",50,0.,50.);
2c1131dd	327	h5->Sumw2();
	328	Add2RecPointsList(h5, 5);
	329
	330	// Correlation plot
	331
	332	TH2I *h6 = new TH2I("hPre10","Cluster - DDL1 vs DDL0", 100,0,200,100,0,200);
	333	Add2RecPointsList(h6,6);
	334
	335	TH2I *h7 = new TH2I("hPre32","Cluster - DDL3 vs DDL2", 100,0,200,100,0,200);
	336	Add2RecPointsList(h7,7);
	337
	338	TH2I *h8 = new TH2I("hCpv54","Cluster - DDL5 vs DDL4", 100,0,200,100,0,200);
	339	Add2RecPointsList(h8,8);
a9523c3e	340
a9523c3e	341
	342
	343	}
	344
	345	//____________________________________________________________________________
	346
	347	void AliPMDQADataMakerRec::InitESDs()
	348	{
	349	//Create histograms to controll ESD
2c1131dd	350
2c1131dd	351	TH1F *h0 = new TH1F("hPreClADC","Cluster ADC of PRE plane",500,0.,10000.);
a9523c3e	352	h0->Sumw2();
	353	Add2ESDsList(h0, 0) ;
	354
2c1131dd	355	TH1F *h1 = new TH1F("hCpvClADC","Cluster ADC of CPV plane",500,0.,10000.);
a9523c3e	356	h1->Sumw2();
	357	Add2ESDsList(h1, 1) ;
	358
2c1131dd	359	TH2I *h2 = new TH2I("hPmdClMult","Cluster Multiplicity: PRE vs. CPVplane",100,0,1000,100,0,1000);
a9523c3e	360	h2->Sumw2();
	361	Add2ESDsList(h2, 2) ;
	362
2c1131dd	363	/*
2c1131dd	364	TH1I * h3 = new TH1I("hCpvClMult","Cluster Multiplicity of CPV plane",100,0.,1000.);
a9523c3e	365	h3->Sumw2();
a9523c3e	366	Add2ESDsList(h3, 3) ;
2c1131dd	367	*/
a9523c3e	368
	369	}
	370
	371	//____________________________________________________________________________
	372	void AliPMDQADataMakerRec::MakeRaws(AliRawReader* rawReader)
	373	{
	374	//Fill prepared histograms with Raw digit properties
	375
2c1131dd	376	TObjArray *pmdddlcont = 0x0;
2c1131dd	377	pmdddlcont = new TObjArray();
a9523c3e	378	AliPMDRawStream stream(rawReader);
a9523c3e	379
2c1131dd	380	AliPMDddldata *pmdddl = 0x0;
a9523c3e	381
2c1131dd	382	Int_t iddl = -1;
	383	Int_t xpad = -1;
	384	Int_t ypad = -1;
	385	Float_t xx, yy;
a9523c3e	386
2c1131dd	387	AliPMDUtility cc;
	388
	389	while ((iddl = stream.DdlData(pmdddlcont)) >=0)
a9523c3e	390	{
2c1131dd	391	Int_t ientries = pmdddlcont->GetEntries();
a9523c3e	392	//printf(" ======= DDLNO = %d ientries = %d \n", iddl, ientries);
2c1131dd	393
a9523c3e	394	for (Int_t ient = 0; ient < ientries; ient++)
2c1131dd	395	{
	396	//AliPMDddldata pmdddl = (AliPMDddldata)pmdddlcont->UncheckedAt(ient);
	397	pmdddl = (AliPMDddldata*)pmdddlcont->UncheckedAt(ient);
	398
	399	Int_t det = pmdddl->GetDetector();
	400	Int_t smn = pmdddl->GetSMN();
	401	Int_t mcm = pmdddl->GetMCM();
a9523c3e	402	//Int_t chno = pmdddl->GetChannel();
2c1131dd	403	Int_t row = pmdddl->GetRow();
2c1131dd	404	Int_t col = pmdddl->GetColumn();
a9523c3e	405	Int_t sig = pmdddl->GetSignal();
a9523c3e	406
2c1131dd	407	if (mcm == 0 \|\| smn == -1) continue;
	408
	409	if (det == 0)
	410	{
	411	GetRawsData(smn)->Fill(sig);
	412	if(smn < 12)
	413	{
	414	xpad = col;
	415	ypad = row;
	416	}
	417	else if(smn >= 12 && smn < 24)
	418	{
	419	xpad = row;
	420	ypad = col;
	421	}
	422	cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
	423	GetRawsData(48)->Fill(xx,yy);
	424	}
	425	if (det == 1)
	426	{
	427	GetRawsData(24+smn)->Fill(sig);
	428	if(smn < 12)
	429	{
	430	xpad = col;
	431	ypad = row;
	432	}
	433	else if(smn >= 12 && smn < 24)
	434	{
	435	xpad = row;
	436	ypad = col;
	437	}
	438
	439	cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
	440	GetRawsData(49)->Fill(xx,yy);
	441
	442	}
	443
	444	}
	445
	446	pmdddlcont->Delete();
a9523c3e	447	}
2c1131dd	448
	449	delete pmdddlcont;
	450	pmdddlcont = 0x0;
	451
a9523c3e	452	}
	453	//____________________________________________________________________________
	454	void AliPMDQADataMakerRec::MakeRecPoints(TTree * clustersTree)
	455	{
	456	// makes data from RecPoints
	457
b5144d1f	458	Int_t multDdl0 = 0, multDdl1 = 0, multDdl2 = 0;
b5144d1f	459	Int_t multDdl3 = 0, multDdl4 = 0, multDdl5 = 0;
2c1131dd	460
a9523c3e	461	AliPMDrecpoint1 * recpoint;
a9523c3e	462
2c1131dd	463	TClonesArray * recpoints = 0x0;
2c1131dd	464	recpoints = new TClonesArray("AliPMDrecpoint1", 1000) ;
a9523c3e	465
	466	TBranch * branch = clustersTree->GetBranch("PMDRecpoint") ;
	467	branch->SetAddress(&recpoints) ;
	468
	469	if ( ! branch )
	470	{
b5144d1f	471	AliWarning("PMD branch in Recpoints Tree not found") ;
a9523c3e	472	}
	473	else
	474	{
	475	for (Int_t imod = 0; imod < branch->GetEntries(); imod++)
	476	{
	477	branch->GetEntry(imod) ;
	478
	479	TIter next(recpoints) ;
	480
	481	while ( (recpoint = dynamic_cast<AliPMDrecpoint1 *>(next())) )
2c1131dd	482	{
a9523c3e	483	//Float_t xpos = recpoint->GetClusX();
a9523c3e	484	//Float_t ypos = recpoint->GetClusY();
2c1131dd	485	//Int_t smn = recpoint->GetSMNumber();
	486
	487	if(recpoint->GetDetector() == 0)
	488	{
	489	if(recpoint->GetSMNumber() >= 0 && recpoint->GetSMNumber() < 6)
	490	{
	491	GetRecPointsData(0)->Fill(recpoint->GetClusCells());
b5144d1f	492	multDdl0++;
2c1131dd	493	}
	494	if(recpoint->GetSMNumber() >= 6 && recpoint->GetSMNumber() < 12)
	495	{
	496	GetRecPointsData(1)->Fill(recpoint->GetClusCells());
b5144d1f	497	multDdl1++;
2c1131dd	498	}
	499	if(recpoint->GetSMNumber() >= 12 && recpoint->GetSMNumber() < 18)
	500	{
	501	GetRecPointsData(2)->Fill(recpoint->GetClusCells());
b5144d1f	502	multDdl2++;
2c1131dd	503	}
	504	if(recpoint->GetSMNumber() >= 18 && recpoint->GetSMNumber() < 24)
	505	{
b5144d1f	506	GetRecPointsData(3)->Fill(recpoint->GetClusCells());
b5144d1f	507	multDdl3++;
2c1131dd	508	}
2c1131dd	509	}
a9523c3e	510
2c1131dd	511	if(recpoint->GetDetector() == 1)
	512	{
	513	if(recpoint->GetSMNumber() >= 0 && recpoint->GetSMNumber() < 6 \|\|
	514	recpoint->GetSMNumber() >= 18 && recpoint->GetSMNumber() < 24)
	515	{
	516	GetRecPointsData(4)->Fill(recpoint->GetClusCells());
b5144d1f	517	multDdl4++;
2c1131dd	518	}
	519	if(recpoint->GetSMNumber() >= 6 && recpoint->GetSMNumber() < 18 )
	520	{
	521	GetRecPointsData(5)->Fill(recpoint->GetClusCells());
b5144d1f	522	multDdl5++;
2c1131dd	523	}
	524	}
	525	}
a9523c3e	526	}
a9523c3e	527	}
2c1131dd	528
b5144d1f	529	GetRecPointsData(6)->Fill(multDdl0,multDdl1);
	530	GetRecPointsData(7)->Fill(multDdl2,multDdl3);
	531	GetRecPointsData(8)->Fill(multDdl4,multDdl5);
2c1131dd	532
2c1131dd	533	delete recpoints;
a9523c3e	534
a9523c3e	535	}
	536
	537	//____________________________________________________________________________
	538
	539	void AliPMDQADataMakerRec::MakeESDs(AliESDEvent * esd)
	540	{
	541	// make QA data from ESDs
	542
	543	Int_t premul = 0, cpvmul = 0;
2c1131dd	544
2c1131dd	545	for (Int_t icl = 0; icl < esd->GetNumberOfPmdTracks(); icl++)
a9523c3e	546	{
2c1131dd	547	AliESDPmdTrack *pmdtr = esd->GetPmdTrack(icl);
	548
	549	//Int_t det = pmdtr->GetDetector();
	550	//Float_t clsX = pmdtr->GetClusterX();
	551	//Float_t clsY = pmdtr->GetClusterY();
	552	//Float_t clsZ = pmdtr->GetClusterZ();
	553	//Float_t ncell = pmdtr->GetClusterCells();
	554	Float_t adc = pmdtr->GetClusterADC();
	555	//Float_t pid = pmdtr->GetClusterPID();
	556
	557	if (pmdtr->GetDetector() == 0)
a9523c3e	558	{
2c1131dd	559	GetESDsData(0)->Fill(adc);
2c1131dd	560	premul++;
a9523c3e	561	}
2c1131dd	562	if (pmdtr->GetDetector() == 1)
a9523c3e	563	{
2c1131dd	564	GetESDsData(1)->Fill(adc) ;
2c1131dd	565	cpvmul++;
a9523c3e	566	}
a9523c3e	567	}
2c1131dd	568
	569	GetESDsData(2)->Fill(cpvmul,premul) ;
	570	//GetESDsData(3)->Fill(cpvmul) ;
a9523c3e	571	}
	572
	573	//____________________________________________________________________________
	574
	575	void AliPMDQADataMakerRec::StartOfDetectorCycle()
	576	{
	577	//Detector specific actions at start of cycle
	578
	579	}
	580	//____________________________________________________________________________
	581	void AliPMDQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
	582	{
	583	//Detector specific actions at end of cycle
	584	// do the QA checking
	585	AliQAChecker::Instance()->Run(AliQA::kPMD, task, list) ;
	586	}