[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

  const Bool_t saveForCorr = kTRUE;
    // 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,saveForCorr);//Ajay
    //Add2RawsList(h48, 48);//Ajay
    TH2F * h49 = new TH2F("hCpvXY","CPV plane",200,-100.,100.,200,-100.,100.);
    Add2RawsList(h49, 49,saveForCorr);//Ajay
    //Add2RawsList(h49, 49);//Ajay

}
//____________________________________________________________________________ 
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"
78328afd	43	//#include "AliPMDRecoParam.h"
a9523c3e	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
5f7ce82b	78	const Bool_t saveForCorr = kTRUE;
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.);
5f7ce82b	279	Add2RawsList(h48, 48,saveForCorr);//Ajay
5f7ce82b	280	//Add2RawsList(h48, 48);//Ajay
2c1131dd	281	TH2F * h49 = new TH2F("hCpvXY","CPV plane",200,-100.,100.,200,-100.,100.);
5f7ce82b	282	Add2RawsList(h49, 49,saveForCorr);//Ajay
5f7ce82b	283	//Add2RawsList(h49, 49);//Ajay
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
5f7ce82b	378	TObjArray *pmdddlcont = 0x0;
2c1131dd	379	pmdddlcont = new TObjArray();
a9523c3e	380	AliPMDRawStream stream(rawReader);
a9523c3e	381
2c1131dd	382	AliPMDddldata *pmdddl = 0x0;
a9523c3e	383
2c1131dd	384	Int_t iddl = -1;
	385	Int_t xpad = -1;
	386	Int_t ypad = -1;
	387	Float_t xx, yy;
a9523c3e	388
2c1131dd	389	AliPMDUtility cc;
	390
	391	while ((iddl = stream.DdlData(pmdddlcont)) >=0)
a9523c3e	392	{
2c1131dd	393	Int_t ientries = pmdddlcont->GetEntries();
a9523c3e	394	//printf(" ======= DDLNO = %d ientries = %d \n", iddl, ientries);
2c1131dd	395
a9523c3e	396	for (Int_t ient = 0; ient < ientries; ient++)
2c1131dd	397	{
	398	//AliPMDddldata pmdddl = (AliPMDddldata)pmdddlcont->UncheckedAt(ient);
	399	pmdddl = (AliPMDddldata*)pmdddlcont->UncheckedAt(ient);
	400
	401	Int_t det = pmdddl->GetDetector();
	402	Int_t smn = pmdddl->GetSMN();
	403	Int_t mcm = pmdddl->GetMCM();
a9523c3e	404	//Int_t chno = pmdddl->GetChannel();
2c1131dd	405	Int_t row = pmdddl->GetRow();
2c1131dd	406	Int_t col = pmdddl->GetColumn();
a9523c3e	407	Int_t sig = pmdddl->GetSignal();
a9523c3e	408
2c1131dd	409	if (mcm == 0 \|\| smn == -1) continue;
	410
	411	if (det == 0)
	412	{
	413	GetRawsData(smn)->Fill(sig);
	414	if(smn < 12)
	415	{
	416	xpad = col;
	417	ypad = row;
	418	}
	419	else if(smn >= 12 && smn < 24)
	420	{
	421	xpad = row;
	422	ypad = col;
	423	}
	424	cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
	425	GetRawsData(48)->Fill(xx,yy);
	426	}
	427	if (det == 1)
	428	{
	429	GetRawsData(24+smn)->Fill(sig);
	430	if(smn < 12)
	431	{
	432	xpad = col;
	433	ypad = row;
	434	}
	435	else if(smn >= 12 && smn < 24)
	436	{
	437	xpad = row;
	438	ypad = col;
	439	}
	440
	441	cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
	442	GetRawsData(49)->Fill(xx,yy);
	443
	444	}
	445
	446	}
	447
	448	pmdddlcont->Delete();
a9523c3e	449	}
2c1131dd	450
	451	delete pmdddlcont;
	452	pmdddlcont = 0x0;
	453
a9523c3e	454	}
	455	//____________________________________________________________________________
	456	void AliPMDQADataMakerRec::MakeRecPoints(TTree * clustersTree)
	457	{
	458	// makes data from RecPoints
	459
b5144d1f	460	Int_t multDdl0 = 0, multDdl1 = 0, multDdl2 = 0;
b5144d1f	461	Int_t multDdl3 = 0, multDdl4 = 0, multDdl5 = 0;
2c1131dd	462
a9523c3e	463	AliPMDrecpoint1 * recpoint;
a9523c3e	464
2c1131dd	465	TClonesArray * recpoints = 0x0;
2c1131dd	466	recpoints = new TClonesArray("AliPMDrecpoint1", 1000) ;
a9523c3e	467
	468	TBranch * branch = clustersTree->GetBranch("PMDRecpoint") ;
	469	branch->SetAddress(&recpoints) ;
	470
	471	if ( ! branch )
	472	{
b5144d1f	473	AliWarning("PMD branch in Recpoints Tree not found") ;
a9523c3e	474	}
	475	else
	476	{
	477	for (Int_t imod = 0; imod < branch->GetEntries(); imod++)
	478	{
	479	branch->GetEntry(imod) ;
	480
	481	TIter next(recpoints) ;
	482
	483	while ( (recpoint = dynamic_cast<AliPMDrecpoint1 *>(next())) )
2c1131dd	484	{
a9523c3e	485	//Float_t xpos = recpoint->GetClusX();
a9523c3e	486	//Float_t ypos = recpoint->GetClusY();
2c1131dd	487	//Int_t smn = recpoint->GetSMNumber();
	488
	489	if(recpoint->GetDetector() == 0)
	490	{
	491	if(recpoint->GetSMNumber() >= 0 && recpoint->GetSMNumber() < 6)
	492	{
	493	GetRecPointsData(0)->Fill(recpoint->GetClusCells());
b5144d1f	494	multDdl0++;
2c1131dd	495	}
	496	if(recpoint->GetSMNumber() >= 6 && recpoint->GetSMNumber() < 12)
	497	{
	498	GetRecPointsData(1)->Fill(recpoint->GetClusCells());
b5144d1f	499	multDdl1++;
2c1131dd	500	}
	501	if(recpoint->GetSMNumber() >= 12 && recpoint->GetSMNumber() < 18)
	502	{
	503	GetRecPointsData(2)->Fill(recpoint->GetClusCells());
b5144d1f	504	multDdl2++;
2c1131dd	505	}
	506	if(recpoint->GetSMNumber() >= 18 && recpoint->GetSMNumber() < 24)
	507	{
b5144d1f	508	GetRecPointsData(3)->Fill(recpoint->GetClusCells());
b5144d1f	509	multDdl3++;
2c1131dd	510	}
2c1131dd	511	}
a9523c3e	512
2c1131dd	513	if(recpoint->GetDetector() == 1)
2c1131dd	514	{
d6cc2277	515	if((recpoint->GetSMNumber() >= 0 && recpoint->GetSMNumber() < 6) \|\|
d6cc2277	516	(recpoint->GetSMNumber() >= 18 && recpoint->GetSMNumber() < 24))
2c1131dd	517	{
2c1131dd	518	GetRecPointsData(4)->Fill(recpoint->GetClusCells());
b5144d1f	519	multDdl4++;
2c1131dd	520	}
	521	if(recpoint->GetSMNumber() >= 6 && recpoint->GetSMNumber() < 18 )
	522	{
	523	GetRecPointsData(5)->Fill(recpoint->GetClusCells());
b5144d1f	524	multDdl5++;
2c1131dd	525	}
	526	}
	527	}
a9523c3e	528	}
a9523c3e	529	}
2c1131dd	530
b5144d1f	531	GetRecPointsData(6)->Fill(multDdl0,multDdl1);
	532	GetRecPointsData(7)->Fill(multDdl2,multDdl3);
	533	GetRecPointsData(8)->Fill(multDdl4,multDdl5);
2c1131dd	534
2c1131dd	535	delete recpoints;
a9523c3e	536
a9523c3e	537	}
	538
	539	//____________________________________________________________________________
	540
	541	void AliPMDQADataMakerRec::MakeESDs(AliESDEvent * esd)
	542	{
	543	// make QA data from ESDs
	544
	545	Int_t premul = 0, cpvmul = 0;
2c1131dd	546
2c1131dd	547	for (Int_t icl = 0; icl < esd->GetNumberOfPmdTracks(); icl++)
a9523c3e	548	{
2c1131dd	549	AliESDPmdTrack *pmdtr = esd->GetPmdTrack(icl);
	550
	551	//Int_t det = pmdtr->GetDetector();
	552	//Float_t clsX = pmdtr->GetClusterX();
	553	//Float_t clsY = pmdtr->GetClusterY();
	554	//Float_t clsZ = pmdtr->GetClusterZ();
	555	//Float_t ncell = pmdtr->GetClusterCells();
	556	Float_t adc = pmdtr->GetClusterADC();
	557	//Float_t pid = pmdtr->GetClusterPID();
	558
	559	if (pmdtr->GetDetector() == 0)
a9523c3e	560	{
2c1131dd	561	GetESDsData(0)->Fill(adc);
2c1131dd	562	premul++;
a9523c3e	563	}
2c1131dd	564	if (pmdtr->GetDetector() == 1)
a9523c3e	565	{
2c1131dd	566	GetESDsData(1)->Fill(adc) ;
2c1131dd	567	cpvmul++;
a9523c3e	568	}
a9523c3e	569	}
2c1131dd	570
	571	GetESDsData(2)->Fill(cpvmul,premul) ;
	572	//GetESDsData(3)->Fill(cpvmul) ;
a9523c3e	573	}
	574
	575	//____________________________________________________________________________
	576
	577	void AliPMDQADataMakerRec::StartOfDetectorCycle()
	578	{
	579	//Detector specific actions at start of cycle
	580
	581	}
	582	//____________________________________________________________________________
57acd2d2	583	void AliPMDQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray ** list)
a9523c3e	584	{
	585	//Detector specific actions at end of cycle
	586	// do the QA checking
	587	AliQAChecker::Instance()->Run(AliQA::kPMD, task, list) ;
	588	}