[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 "AliRawReader.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

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