[u/mrichter/AliRoot.git] / PMD / AliPMDQualitycheck.C

//********************************************************************
//     Macro that checks PMD data Quality 
//                                         using the ESD classes
//     Modified by P.V.K.S.BABA (J.U.)/Ajay(IOP) for doing Quality check.
//********************************************************************

#if !defined( __CINT__) || defined(__MAKECINT__)
#include <Riostream.h>
#include "TTree.h"
#include "TFile.h"
#include "TH1F.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TStopwatch.h"

#include "AliESD.h"
#endif
double Xcon[96] =
  { 75.133, 54.204, 53.254, 32.326, 31.376,10.447,
    75.133, 54.204, 53.254, 32.326, 31.376,10.447,
    75.133, 54.204, 53.254, 32.326, 31.376,10.447,
    75.133, 54.204, 53.254, 32.326, 31.376,10.447,
    -75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
    -75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
    -75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
    -75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
    9.167, -32.543, -33.493, -75.133,
    9.167, -32.543, -33.493, -75.133,
    9.167, -32.543, -33.493, -75.133,
    9.167, -32.543, -33.493, -75.133,
    9.167, -32.543, -33.493, -75.133,
    9.167, -32.543, -33.493, -75.133,
    -9.167, 32.543, 33.493, 75.133,
    -9.167, 32.543, 33.493, 75.133,
    -9.167, 32.543, 33.493, 75.133,
    -9.167, 32.543, 33.493, 75.133,
    -9.167, 32.543, 33.493, 75.133,
    -9.167, 32.543, 33.493, 75.133,};

double Ycon[96] =
  {86.475,  86.475,  86.475, 86.475,  86.475,  86.475,
   38.225,  38.225,  38.225, 38.225,  38.225,  38.225,
   37.325,  37.325,  37.325, 37.325,  37.325,  37.325,
   -10.925, -10.925, -10.925, -10.925, -10.925, -10.925,
   -86.475, -86.475, -86.475, -86.475, -86.475, -86.475,
   -38.225,  -38.225,  -38.225, -38.225, -38.225, -38.225,
   -37.325,  -37.325,  -37.325, -37.325,  -37.325,  -37.325
   10.925, 10.925, 10.925, 10.925, 10.925, 10.925,

   86.475,  86.475, 86.475,  86.475,
   62.225,  62.225, 62.225,  62.225,
   61.325,  61.325, 61.325,  61.325,
   37.075, 37.075, 37.075, 37.075
   36.175,  36.175, 36.175,  36.175,
   11.925, 11.925, 11.925 , 11.925,
   -86.475,  -86.475, -86.475,  -86.475,
   -62.225,  -62.225, -62.225,  -62.225,
   -61.325,  -61.325, -61.325,  -61.325,
   -37.075,  -37.075, -37.075,  -37.075
   -36.175,  -36.175, -36.175,  -36.175,
   -11.925, -11.925,  -11.925 , -11.925 };


extern TStyle *gStyle;

// Default is that it will take 1000 Events sample. Want to process less give (Number) .
Int_t AliPMDQualitycheck(Int_t nevt=1000) { 
AliCDBManager::Instance()->SetRun(0);
  TStopwatch timer;
  char PRINTSMN ;

//  printf("Whether You want plotof All SMN , Give y or n ? \n " ) ;
//  scanf("%1s" , &PRINTSMN );
 
  TH2F *hP1 = new TH2F("hP1","XY of Clusters",100,-100.,100.,100,-100.,100.);
  TH1F *hC2 = new TH1F("hC2","CPV  PHI",200,-1,9);
  TH1F *hP2 = new TH1F("hP2","PRE  PHI",200,-1,9);
  TH1F *hC3 = new TH1F("hC3","CPV  Clus",30,0.,500.);
  TH1F *hP3 = new TH1F("hP3","PRE  N-gammalike",20,0.,500.);
  TH1F *hP4 = new TH1F("hP4","PRE  EDEP",30,0.,1000.);
  TH1F *hC5 = new TH1F("hC5","CPV  n-cell",20,0.,100.);
  TH1F *hP5 = new TH1F("hP5","PMD  n-cell",20,0.,100.);
  TH2F *hCP0 = new TH2F("hCP0","PRE CLUS Quad.1 vs 2",150,0.,300.,150,0.,300.);
  TH2F *hCP1 = new TH2F("hCP1","PRE CLUS Quad.3 vs 4",150,0.,300.,150,0.,300.);
  TH2F *hCP2 = new TH2F("hCP2","PRE EDEP Quad.3 vs 4",50,0.,300.,50,0.,300.);
  TH2F *hCP3 = new TH2F("hCP3","PRE EDEP vs Tot Clus ",10,0.,1000.,10,0.,300.);
  TH2F *hCP4 = new TH2F("hCP4","PRE Clus vs CPV Clus ",150,0.,200.,150,0.,200.);

  TH2F *hSM1 = new TH2F("hSM1","PRE Cluster XY",200,-100,100,200,-100,100);
  TH2F *hSM2 = new TH2F("hSM2","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM3 = new TH2F("hSM3","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM4 = new TH2F("hSM4","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM5 = new TH2F("hSM5","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM6 = new TH2F("hSM6","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM7 = new TH2F("hSM7","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM8 = new TH2F("hSM8","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM9 = new TH2F("hSM9","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM10 = new TH2F("hSM10","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM11 = new TH2F("hSM11","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM12 = new TH2F("hSM12","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM13 = new TH2F("hSM13","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM14 = new TH2F("hSM14","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM15 = new TH2F("hSM15","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM16 = new TH2F("hSM16","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM17 = new TH2F("hSM17","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM18 = new TH2F("hSM18","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM19 = new TH2F("hSM19","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM20 = new TH2F("hSM20","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM21 = new TH2F("hSM21","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM22 = new TH2F("hSM22","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM23 = new TH2F("hSM23","",999,-100.0,100.0,999,-100.0,100.0);
  TH2F *hSM24 = new TH2F("hSM24","",999,-100.0,100.0,999,-100.0,100.0);
  TH1F *hSM = new TH1F("hSM","Plot of all 24 Super Modules",24,0,24);

// Star processing.
  
  AliPMDUtility *cc = new AliPMDUtility(); 
 
  TFile *ef=TFile::Open("AliESDtree.root");
  if (!ef || !ef->IsOpen()) {cerr<<"Cant open AliESDtree.root !\n"; return 1;}
  AliESD* event = new AliESD;
  TTree* tree = (TTree*) ef->Get("esdTree");
  if (!tree) {cerr<<"no ESD tree found\n"; return 1;};
  tree->SetBranchAddress("ESD", &event);
  Int_t SMN;
  Int_t n=0;
  Float_t TotCPVClus     ;
  Float_t TotPREClus     ;
  Float_t TotPREEdep  ;
  Float_t TotCPVCell     ;
  Float_t TotPRECell     ;
  Float_t PREcluQUAD[4]  ;
  Float_t CPVcluQUAD[4]  ;
  Float_t PREadcQUAD[4]  ;
  Float_t CPVADCQUAD[4]  ;
  Float_t PREcelQUAD[4]  ;
  Float_t CPVcelQUAD[4]  ;
 
 
  //******* The loop over events
  while (tree->GetEvent(n)) {
    cout<<endl<<"Processing event number : "<<n++<<endl;
    
    
    Int_t CPVhits = 0 ;
    Int_t PMDclus = 0 ;
    Int_t npmdcl=event->GetNumberOfPmdTracks();
    //cout<<"Number of PMD tracks : "<<npmdcl<<endl; 
    
    //****** The loop over PMD clusters

    for (Int_t kk=0; kk<4;kk++) {
        CPVcluQUAD[kk] = 0.0 ;
        PREcluQUAD[kk] = 0.0 ;
        CPVcelQUAD[kk] = 0.0 ;
        PREcelQUAD[kk] = 0.0 ;
        PREadcQUAD[kk] = 0.0 ;
     }
    while (npmdcl--) {
      AliESDPmdTrack *pmdtr = event->GetPmdTrack(npmdcl);
      
      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();

//    cout << " "  << det << " " << ncell << " " << adc << " " << pid << endl ;      
      cc->SetXYZ(clsX,clsY,clsZ);
      //cc->SetPxPyPz(clsX,clsY,clsZ);
      cc->CalculateEta();
      cc->CalculatePhi();
      Float_t eta = cc->GetEta();
      Float_t phi = cc->GetPhi();
 
// Calculating S.Module Number from Cluster .
            CalculateSMN(clsX, clsY,SMN);
      if( det == 1)
	{
          if(SMN >= 0 && SMN <= 5) {
                 ++CPVcluQUAD[0] ;
           CPVADCQUAD[0] =+ adc  ;
           CPVcelQUAD[0] =+ ncell ;
          }
          if(SMN >= 6 && SMN <=11) {
                 ++CPVcluQUAD[1] ;
           CPVADCQUAD[1] =+ adc  ;
           CPVcelQUAD[1] =+ ncell ;
          }
          if(SMN >=12 && SMN <=17) {
                 ++CPVcluQUAD[2] ;
           CPVADCQUAD[2] =+ adc  ;
           CPVcelQUAD[2] =+ ncell ;
          }
          if(SMN >=18 && SMN <=23) {
                 ++CPVcluQUAD[3] ;
           CPVADCQUAD[3] =+ adc  ;
           CPVcelQUAD[3] =+ ncell ;
          }

	  if(eta >= 2.3 && eta <= 3.5)
	  {
	    hC2->Fill(phi);
	  }
	}
      if( det == 0)
	{
          if(SMN >= 0 && SMN <= 5) { 
                 ++PREcluQUAD[0] ;
           PREadcQUAD[0] =+ adc  ;    
           PREcelQUAD[0] =+ ncell ;    
          }    
          if(SMN >= 6 && SMN <=11) { 
                 ++PREcluQUAD[1] ;
           PREadcQUAD[1] =+ adc  ;    
           PREcelQUAD[1] =+ ncell ;    
          }    
          if(SMN >=12 && SMN <=17) { 
                 ++PREcluQUAD[2] ;
           PREadcQUAD[2] =+ adc  ;    
           PREcelQUAD[2] =+ ncell ;    
          }    
          if(SMN >=18 && SMN <=23) { 
                 ++PREcluQUAD[3] ;
           PREadcQUAD[3] =+ adc  ;    
           PREcelQUAD[3] =+ ncell ;    
          }    
          if ( n <= 100 ) { 
                         hSM->Fill(SMN);
            if(SMN == 0)hSM1->Fill(-clsX,clsY);
            if(SMN == 0)hSM1->Fill(-clsX,clsY);
            if(SMN == 1)hSM2->Fill(-clsX,clsY);
            if(SMN == 2)hSM3->Fill(-clsX,clsY);
            if(SMN == 3)hSM4->Fill(-clsX,clsY);
            if(SMN == 4)hSM5->Fill(-clsX,clsY);
            if(SMN == 5)hSM6->Fill(-clsX,clsY);
            if(SMN == 6)hSM7->Fill(-clsX,clsY);
            if(SMN == 7)hSM8->Fill(-clsX,clsY);
            if(SMN == 8)hSM9->Fill(-clsX,clsY);
            if(SMN == 9)hSM10->Fill(-clsX,clsY);
            if(SMN ==10)hSM11->Fill(-clsX,clsY);
            if(SMN ==11)hSM12->Fill(-clsX,clsY);
            if(SMN ==12)hSM13->Fill(-clsX,clsY);
            if(SMN ==13)hSM14->Fill(-clsX,clsY);
            if(SMN ==14)hSM15->Fill(-clsX,clsY);
            if(SMN ==15)hSM16->Fill(-clsX,clsY);
            if(SMN ==16)hSM17->Fill(-clsX,clsY);
            if(SMN ==17)hSM18->Fill(-clsX,clsY);
            if(SMN ==18)hSM19->Fill(-clsX,clsY);
            if(SMN ==19)hSM20->Fill(-clsX,clsY);
            if(SMN ==20)hSM21->Fill(-clsX,clsY);
            if(SMN ==21)hSM22->Fill(-clsX,clsY);
            if(SMN ==22)hSM23->Fill(-clsX,clsY);
            if(SMN ==23)hSM24->Fill(-clsX,clsY);
          }     
	  if(eta >= 2.3 && eta <= 3.5)
	  {
	    hP2->Fill(phi);
	  }
	    hP1->Fill(clsX,clsY);
	}
    }
    for (Int_t k=0; k<4;k++) {
        TotCPVClus =+ CPVcluQUAD [k] ;
        TotPREClus =+ PREcluQUAD [k] ;
        TotCPVCell =+ CPVcelQUAD [k] ;
        TotPRECell =+ PREcelQUAD [k] ;
        TotPREEdep =+ PREadcQUAD [k] ;     
    }
        Float_t TotCPVPREClus = TotPREClus + TotCPVClus ;
        Float_t TotCPVPRECell = TotPRECell + TotCPVCell ;
	  if(eta >= 2.3 && eta <= 3.5)
	  {
	    hC3->Fill(TotCPVClus);
	    hP3->Fill(TotPREClus);
	    hP4->Fill(TotPREEdep);
	    hP5->Fill(TotPRECell);
           hCP0->Fill(PREcluQUAD[0],PREcluQUAD[1]);
           hCP1->Fill(PREcluQUAD[2],PREcluQUAD[3]);
           hCP2->Fill(PREadcQUAD[2],PREadcQUAD[3]);
           hCP3->Fill(TotPREEdep,TotCPVPREClus);
           hCP4->Fill(TotPREClus,TotCPVClus);
	  }
           TotCPVClus == 0.0; 
           TotPREClus == 0.0; 
           TotCPVCell == 0.0; 
           TotPRECell == 0.0; 
           TotPREEdep == 0.0; 
    if(n >= nevt)break;    
  }

  gStyle->SetOptStat(110000);
  gStyle->SetOptFit(1);

//  if ( PRINTSMN == 'y' ) {
  TCanvas *cP1 = new TCanvas("cP1","Cluster XY", 10,10, 600, 600);
  cP1->Range(-100, -100,100 ,100 );
  hSM1->SetMarkerColor(2);
  hSM1->Draw();
  hSM1->GetXaxis()->SetTitle("Cluster X");
  hSM1->GetYaxis()->SetTitle("Cluster Y");
  hSM2->SetMarkerColor(2);
  hSM2->Draw("same");
  hSM3->SetMarkerColor(2);
  hSM3->Draw("same");
  hSM4->SetMarkerColor(2);
  hSM4->Draw("same");
  hSM5->SetMarkerColor(2);
  hSM5->Draw("same");
  hSM6->SetMarkerColor(2);
  hSM6->Draw("same");
  hSM7->SetMarkerColor(4);
  hSM7->Draw("same");
  hSM8->SetMarkerColor(4);
  hSM8->Draw("same");
  hSM9->SetMarkerColor(4);
  hSM9->Draw("same");
  hSM10->SetMarkerColor(4);
  hSM10->Draw("same");
  hSM11->SetMarkerColor(4);
  hSM11->Draw("same");
  hSM12->SetMarkerColor(4);
  hSM12->Draw("same");
  hSM13->SetMarkerColor(6);
  hSM13->Draw("same");
  hSM14->SetMarkerColor(6);
  hSM14->Draw("same");
  hSM15->SetMarkerColor(6);
  hSM15->Draw("same");
  hSM16->SetMarkerColor(6);
  hSM16->Draw("same");
  hSM17->SetMarkerColor(6);
  hSM17->Draw("same");
  hSM18->SetMarkerColor(6);
  hSM18->Draw("same");
  hSM19->SetMarkerColor(8);
  hSM19->Draw("same");
  hSM20->SetMarkerColor(8);
  hSM20->Draw("same");
  hSM21->SetMarkerColor(8);
  hSM21->Draw("same");
  hSM22->SetMarkerColor(8);
  hSM22->Draw("same");
  hSM23->SetMarkerColor(8);
  hSM23->Draw("same");
  hSM24->SetMarkerColor(8);
  hSM24->Draw("same");

  Int_t linsav = gStyle->GetLineWidth();
  DrawPMDBoundary();
  DrawPMDBoundarySM1();
  DrawPMDBoundarySM2();
  DrawPMDBoundarySM3();
  DrawPMDBoundarySM4();
  cP1->Print("ClusterXY.gif");
  

  TCanvas *cP2 = new TCanvas("PHI CPV / PMD","",10,10, 600,600);
  cP2->Divide(1,2);
  cP2->cd(1);
  cP2->SetFillColor(0);
  hC2->SetLineColor(4);
  hC2->Draw();
  cP2->cd(2);
  hP2->SetLineColor(2);
  hP2->Draw();
  cP2->Print("CPVPREphi.gif");

// }

  TCanvas *cP2 = new TCanvas("cP2","",10,10,600,600);
  cP2->cd();
  hSM->SetFillColor(2);
  hSM->Draw();
  cP2->Print("AllSMN.gif");


  TCanvas *cP3 = new TCanvas("CPV Clus  PRE Clus Correlations", " ",10,10, 600,600);
  cP3->Divide(2,2);
  cP3->cd(1);
  hCP0->SetMarkerColor(9);
  hCP0->Draw();
  cP3->cd(2);
  hCP1->SetMarkerColor(6);
  hCP1->Draw();
  cP3->cd(3);
  hP3->SetLineColor(2);
  hP3->Draw();
  cP3->cd(4);
  hCP4->SetMarkerColor(3);
  hCP4->Draw();
  cP3->Print("CPVPREClus.gif");

  TCanvas *cP6 = new TCanvas("CPV Clus PRE Adc"," ",10,10, 600,600);
  cP6->Divide(1,2);
  cP6->cd(1);
  hC3->SetLineColor(4);
  hC3->Draw();
  cP6->cd(2);
  hP4->SetLineColor(2);
  hP4->Draw();
  cP6->Print("CPVPREAdc.gif");

  timer.Stop();
  timer.Print();
  return 0;
  }

void CalculateSMN( Float_t clsX, Float_t clsY, Int_t & smn)
{
  //smn = 0;

  //---------------------------------------------------------------------
  if((clsX <= Xcon[0]) && (clsX >= Xcon[1]) &&
    (clsY <= Ycon[0]) && (clsY >= Ycon[6])) smn = 0 ;

  else if((clsX <=Xcon[2]) && (clsX >= Xcon[3]) &&
  (clsY <= Ycon[1]) && (clsY >= Ycon[7]))smn = 1 ;

  else if((clsX <=Xcon[4]) && (clsX >= Xcon[5]) &&
  (clsY <= Ycon[3]) && (clsY >= Ycon[8]))smn = 2 ;

  else if((clsX <= Xcon[0]) && (clsX >= Xcon[1]) &&
  (clsY <= Ycon[12]) && (clsY >= Ycon[18])) smn = 3 ;

  else if((clsX <=Xcon[2]) && (clsX >= Xcon[3]) &&
  (clsY <= Ycon[12]) && (clsY >= Ycon[18]))smn = 4 ;

  else if((clsX <=Xcon[4]) && (clsX >= Xcon[5]) &&
  (clsY <= Ycon[12]) && (clsY >= Ycon[18]))smn = 5 ;
  //------------------------------------------------------------------
  else if((clsX >= Xcon[24]) && (clsX <= Xcon[25]) &&
  (clsY >= Ycon[24]) && (clsY <= Ycon[30])) smn = 6 ;

  else if((clsX >=Xcon[26]) && (clsX <= Xcon[27]) &&
  (clsY >= Ycon[25]) && (clsY <= Ycon[31]))smn = 7 ;

  else if((clsX >=Xcon[28]) && (clsX <= Xcon[29]) &&
  (clsY >= Ycon[26]) && (clsY <= Ycon[32]))smn = 8 ;

  else if((clsX >= Xcon[24]) && (clsX <= Xcon[25]) &&
  (clsY >= Ycon[36]) && (clsY <= Ycon[42])) smn = 9 ;

  else if((clsX >=Xcon[26]) && (clsX <= Xcon[27]) &&
  (clsY >= Ycon[36]) && (clsY <= Ycon[42]))smn = 10;

  else if((clsX >=Xcon[28]) && (clsX <= Xcon[29]) &&
  (clsY >= Ycon[36]) && (clsY <= Ycon[42]))smn = 11;
  //------------------------------------------------------------------
  else if((clsX <= Xcon[48]) && (clsX >= Xcon[49]) &&
  (clsY <= Ycon[48]) && (clsY >= Ycon[52])) smn = 12 ;

  else if((clsX <=Xcon[50]) && (clsX >= Xcon[51]) &&
  (clsY <= Ycon[48]) && (clsY >= Ycon[52]))smn = 13 ;

  else if((clsX <=Xcon[48]) && (clsX >= Xcon[49]) &&
  (clsY <= Ycon[56]) && (clsY >= Ycon[60]))smn = 14 ;

  else if((clsX <=Xcon[50]) && (clsX >= Xcon[51]) &&
  (clsY <= Ycon[56]) && (clsY >= Ycon[60]))smn = 15 ;

  else if((clsX <=Xcon[48]) && (clsX >= Xcon[49]) &&
  (clsY <= Ycon[64]) && (clsY >= Ycon[68]))smn = 16 ;

  else if((clsX <=Xcon[50]) && (clsX >= Xcon[51]) &&
  (clsY <= Ycon[64]) && (clsY >= Ycon[68]))smn = 17 ;
  //--------------------------------------------------------------
  else if((clsX >= Xcon[72]) && (clsX <= Xcon[73]) &&
  (clsY >= Ycon[72]) && (clsY <= Ycon[76])) smn = 18 ;

  else if((clsX >=Xcon[74]) && (clsX <= Xcon[75]) &&
  (clsY >= Ycon[72]) && (clsY <= Ycon[76]))smn = 19 ;

  else if((clsX >=Xcon[72]) && (clsX <= Xcon[73]) &&
  (clsY >= Ycon[80]) && (clsY <= Ycon[84]))smn = 20 ;

  else if((clsX >=Xcon[74]) && (clsX <= Xcon[75]) &&
  (clsY >= Ycon[80]) && (clsY <= Ycon[84]))smn = 21;

  else if((clsX >= Xcon[72]) && (clsX <= Xcon[73]) &&
  (clsY >= Ycon[88]) && (clsY <= Ycon[92])) smn = 22 ;

  else if((clsX >=Xcon[74]) && (clsX <= Xcon[75]) &&
  (clsY >= Ycon[88]) && (clsY <= Ycon[92]))smn = 23 ;
  else smn = 111;

 }

void DrawPMDBoundary()
{
  TH2F *h = new TH2F("h","",200,-100,100,200,-100,100);
  gStyle->SetLineWidth(2);
  gStyle->SetLineColor(2);
  TLine * l;
  l = new TLine(75.1333, 86.475, -75.1333, 86.475); l->Draw("same");
  l = new TLine(-75.1333, 86.470,-75.1333, -86.475); l->Draw("same");
  l = new TLine(-75.1333, -86.475,75.1333, -86.475); l->Draw("same");
  l = new TLine(75.1333, -86.475,75.1333, 86.475); l->Draw("same");
}
void DrawPMDBoundarySM1()
{
  TH2F *hsm1 = new TH2F("hsm1","",200,-100,100,200,-100,100);
  gStyle->SetLineWidth(1);
  gStyle->SetLineColor(4);
  TLine * l;
  l = new TLine(-75.1333, 86.475, -10.447,  86.475); l->Draw("same");
  l = new TLine(-10.447, 86.475, -10.446, -10.925); l->Draw("same");
  l = new TLine(-10.446, -10.925, -75.1333,-10.925); l->Draw("same");
  l = new TLine(-75.1333,-10.925, -75.1333, 86.475); l->Draw("same");
}
void DrawPMDBoundarySM2()
{
  TH2F *hsm2 = new TH2F("hsm2","",200,-100,100,200,-100,100);
  gStyle->SetLineWidth(1);
  gStyle->SetLineColor(4);
  TLine * l;
  l = new TLine(75.1333, -86.475, 10.446,  -86.475); l->Draw("same");
  l = new TLine(10.446,  -86.475, 10.446,  10.925); l->Draw("same");
  l = new TLine(10.446,   10.925, 75.1333, 10.925); l->Draw("same");
  l = new TLine(75.1333,  10.925, 75.1333, -86.475); l->Draw("same");
}

void DrawPMDBoundarySM3()
{
  TH2F *hsm3 = new TH2F("hsm3","",200,-100,100,200,-100,100);
  gStyle->SetLineWidth(1);
  gStyle->SetLineColor(1);
  TLine * l;
  l = new TLine(  -9.167, 86.475, 75.1333, 86.475); l->Draw("same");
  l = new TLine(75.1333,86.475, 75.1333, 11.925); l->Draw("same");
  l = new TLine(75.1333,11.925,   -9.167,  11.925); l->Draw("same");
  l = new TLine(  -9.167, 11.925,   -9.167,  86.475); l->Draw("same");
}
void DrawPMDBoundarySM4()
{
  TH2F *hsm4 = new TH2F("hsm4","",200,-100,100,200,-100,100);
  gStyle->SetLineWidth(1);
  gStyle->SetLineColor(1);
  TLine * l;
  l = new TLine(9.167, -86.475, -75.1333,-86.475); l->Draw("same");
  l = new TLine(-75.1333,-86.475, -75.1333,-11.925); l->Draw("same");
  l = new TLine(-75.1333,-11.925, 9.167, -11.925); l->Draw("same");
  l = new TLine(9.167, -11.925, 9.167, -86.475); l->Draw("same");
}
Commit	Line	Data
92534737	1	//********************************************************************
	2	// Macro that checks PMD data Quality
	3	// using the ESD classes
	4	// Modified by P.V.K.S.BABA (J.U.)/Ajay(IOP) for doing Quality check.
	5	//********************************************************************
	6
	7	#if !defined( __CINT__) \|\| defined(__MAKECINT__)
	8	#include <Riostream.h>
	9	#include "TTree.h"
	10	#include "TFile.h"
	11	#include "TH1F.h"
	12	#include "TCanvas.h"
	13	#include "TStyle.h"
	14	#include "TStopwatch.h"
	15
	16	#include "AliESD.h"
	17	#endif
	18	double Xcon[96] =
	19	{ 75.133, 54.204, 53.254, 32.326, 31.376,10.447,
	20	75.133, 54.204, 53.254, 32.326, 31.376,10.447,
	21	75.133, 54.204, 53.254, 32.326, 31.376,10.447,
	22	75.133, 54.204, 53.254, 32.326, 31.376,10.447,
	23	-75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
	24	-75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
	25	-75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
	26	-75.133, -54.204, -53.254, -32.326, -31.376,-10.447,
	27	9.167, -32.543, -33.493, -75.133,
	28	9.167, -32.543, -33.493, -75.133,
	29	9.167, -32.543, -33.493, -75.133,
	30	9.167, -32.543, -33.493, -75.133,
	31	9.167, -32.543, -33.493, -75.133,
	32	9.167, -32.543, -33.493, -75.133,
	33	-9.167, 32.543, 33.493, 75.133,
	34	-9.167, 32.543, 33.493, 75.133,
	35	-9.167, 32.543, 33.493, 75.133,
	36	-9.167, 32.543, 33.493, 75.133,
	37	-9.167, 32.543, 33.493, 75.133,
	38	-9.167, 32.543, 33.493, 75.133,};
	39
	40	double Ycon[96] =
	41	{86.475, 86.475, 86.475, 86.475, 86.475, 86.475,
	42	38.225, 38.225, 38.225, 38.225, 38.225, 38.225,
	43	37.325, 37.325, 37.325, 37.325, 37.325, 37.325,
	44	-10.925, -10.925, -10.925, -10.925, -10.925, -10.925,
	45	-86.475, -86.475, -86.475, -86.475, -86.475, -86.475,
	46	-38.225, -38.225, -38.225, -38.225, -38.225, -38.225,
	47	-37.325, -37.325, -37.325, -37.325, -37.325, -37.325
	48	10.925, 10.925, 10.925, 10.925, 10.925, 10.925,
	49
	50	86.475, 86.475, 86.475, 86.475,
	51	62.225, 62.225, 62.225, 62.225,
	52	61.325, 61.325, 61.325, 61.325,
	53	37.075, 37.075, 37.075, 37.075
	54	36.175, 36.175, 36.175, 36.175,
	55	11.925, 11.925, 11.925 , 11.925,
	56	-86.475, -86.475, -86.475, -86.475,
	57	-62.225, -62.225, -62.225, -62.225,
	58	-61.325, -61.325, -61.325, -61.325,
	59	-37.075, -37.075, -37.075, -37.075
	60	-36.175, -36.175, -36.175, -36.175,
	61	-11.925, -11.925, -11.925 , -11.925 };
	62
	63
	64	extern TStyle *gStyle;
65
66	// Default is that it will take 1000 Events sample. Want to process less give (Number) .
67	Int_t AliPMDQualitycheck(Int_t nevt=1000) {
68	AliCDBManager::Instance()->SetRun(0);
69	TStopwatch timer;
70	char PRINTSMN ;
71
72	// printf("Whether You want plotof All SMN , Give y or n ? \n " ) ;
73	// scanf("%1s" , &PRINTSMN );
74
75	TH2F *hP1 = new TH2F("hP1","XY of Clusters",100,-100.,100.,100,-100.,100.);
76	TH1F *hC2 = new TH1F("hC2","CPV PHI",200,-1,9);
77	TH1F *hP2 = new TH1F("hP2","PRE PHI",200,-1,9);
78	TH1F *hC3 = new TH1F("hC3","CPV Clus",30,0.,500.);
79	TH1F *hP3 = new TH1F("hP3","PRE N-gammalike",20,0.,500.);
80	TH1F *hP4 = new TH1F("hP4","PRE EDEP",30,0.,1000.);
81	TH1F *hC5 = new TH1F("hC5","CPV n-cell",20,0.,100.);
82	TH1F *hP5 = new TH1F("hP5","PMD n-cell",20,0.,100.);
83	TH2F *hCP0 = new TH2F("hCP0","PRE CLUS Quad.1 vs 2",150,0.,300.,150,0.,300.);
84	TH2F *hCP1 = new TH2F("hCP1","PRE CLUS Quad.3 vs 4",150,0.,300.,150,0.,300.);
85	TH2F *hCP2 = new TH2F("hCP2","PRE EDEP Quad.3 vs 4",50,0.,300.,50,0.,300.);
86	TH2F *hCP3 = new TH2F("hCP3","PRE EDEP vs Tot Clus ",10,0.,1000.,10,0.,300.);
87	TH2F *hCP4 = new TH2F("hCP4","PRE Clus vs CPV Clus ",150,0.,200.,150,0.,200.);
88
89	TH2F *hSM1 = new TH2F("hSM1","PRE Cluster XY",200,-100,100,200,-100,100);
90	TH2F *hSM2 = new TH2F("hSM2","",999,-100.0,100.0,999,-100.0,100.0);
91	TH2F *hSM3 = new TH2F("hSM3","",999,-100.0,100.0,999,-100.0,100.0);
92	TH2F *hSM4 = new TH2F("hSM4","",999,-100.0,100.0,999,-100.0,100.0);
93	TH2F *hSM5 = new TH2F("hSM5","",999,-100.0,100.0,999,-100.0,100.0);
94	TH2F *hSM6 = new TH2F("hSM6","",999,-100.0,100.0,999,-100.0,100.0);
95	TH2F *hSM7 = new TH2F("hSM7","",999,-100.0,100.0,999,-100.0,100.0);
96	TH2F *hSM8 = new TH2F("hSM8","",999,-100.0,100.0,999,-100.0,100.0);
97	TH2F *hSM9 = new TH2F("hSM9","",999,-100.0,100.0,999,-100.0,100.0);
98	TH2F *hSM10 = new TH2F("hSM10","",999,-100.0,100.0,999,-100.0,100.0);
99	TH2F *hSM11 = new TH2F("hSM11","",999,-100.0,100.0,999,-100.0,100.0);
100	TH2F *hSM12 = new TH2F("hSM12","",999,-100.0,100.0,999,-100.0,100.0);
101	TH2F *hSM13 = new TH2F("hSM13","",999,-100.0,100.0,999,-100.0,100.0);
102	TH2F *hSM14 = new TH2F("hSM14","",999,-100.0,100.0,999,-100.0,100.0);
103	TH2F *hSM15 = new TH2F("hSM15","",999,-100.0,100.0,999,-100.0,100.0);
104	TH2F *hSM16 = new TH2F("hSM16","",999,-100.0,100.0,999,-100.0,100.0);
105	TH2F *hSM17 = new TH2F("hSM17","",999,-100.0,100.0,999,-100.0,100.0);
106	TH2F *hSM18 = new TH2F("hSM18","",999,-100.0,100.0,999,-100.0,100.0);
107	TH2F *hSM19 = new TH2F("hSM19","",999,-100.0,100.0,999,-100.0,100.0);
108	TH2F *hSM20 = new TH2F("hSM20","",999,-100.0,100.0,999,-100.0,100.0);
109	TH2F *hSM21 = new TH2F("hSM21","",999,-100.0,100.0,999,-100.0,100.0);
110	TH2F *hSM22 = new TH2F("hSM22","",999,-100.0,100.0,999,-100.0,100.0);
111	TH2F *hSM23 = new TH2F("hSM23","",999,-100.0,100.0,999,-100.0,100.0);
112	TH2F *hSM24 = new TH2F("hSM24","",999,-100.0,100.0,999,-100.0,100.0);
113	TH1F *hSM = new TH1F("hSM","Plot of all 24 Super Modules",24,0,24);
114
115	// Star processing.
116
117	AliPMDUtility *cc = new AliPMDUtility();
118
119	TFile *ef=TFile::Open("AliESDtree.root");
120	if (!ef \|\| !ef->IsOpen()) {cerr<<"Cant open AliESDtree.root !\n"; return 1;}
121	AliESD* event = new AliESD;
122	TTree* tree = (TTree*) ef->Get("esdTree");
123	if (!tree) {cerr<<"no ESD tree found\n"; return 1;};
124	tree->SetBranchAddress("ESD", &event);
125	Int_t SMN;
126	Int_t n=0;
127	Float_t TotCPVClus ;
128	Float_t TotPREClus ;
129	Float_t TotPREEdep ;
130	Float_t TotCPVCell ;
131	Float_t TotPRECell ;
132	Float_t PREcluQUAD[4] ;
133	Float_t CPVcluQUAD[4] ;
134	Float_t PREadcQUAD[4] ;
135	Float_t CPVADCQUAD[4] ;
136	Float_t PREcelQUAD[4] ;
137	Float_t CPVcelQUAD[4] ;
138
139
140	//******* The loop over events
141	while (tree->GetEvent(n)) {
142	cout<<endl<<"Processing event number : "<<n++<<endl;
143
144
145	Int_t CPVhits = 0 ;
146	Int_t PMDclus = 0 ;
147	Int_t npmdcl=event->GetNumberOfPmdTracks();
148	//cout<<"Number of PMD tracks : "<<npmdcl<<endl;
149
150	//****** The loop over PMD clusters
151
152	for (Int_t kk=0; kk<4;kk++) {
153	CPVcluQUAD[kk] = 0.0 ;
154	PREcluQUAD[kk] = 0.0 ;
155	CPVcelQUAD[kk] = 0.0 ;
156	PREcelQUAD[kk] = 0.0 ;
157	PREadcQUAD[kk] = 0.0 ;
158	}
159	while (npmdcl--) {
160	AliESDPmdTrack *pmdtr = event->GetPmdTrack(npmdcl);
161
162	Int_t det = pmdtr->GetDetector();
163	Float_t clsX = pmdtr->GetClusterX();
164	Float_t clsY = pmdtr->GetClusterY();
165	Float_t clsZ = pmdtr->GetClusterZ();
166	Float_t ncell = pmdtr->GetClusterCells();
167	Float_t adc = pmdtr->GetClusterADC();
168	Float_t pid = pmdtr->GetClusterPID();
169
170	// cout << " " << det << " " << ncell << " " << adc << " " << pid << endl ;
171	cc->SetXYZ(clsX,clsY,clsZ);
172	//cc->SetPxPyPz(clsX,clsY,clsZ);
173	cc->CalculateEta();
174	cc->CalculatePhi();
175	Float_t eta = cc->GetEta();
176	Float_t phi = cc->GetPhi();
177
178	// Calculating S.Module Number from Cluster .
179	CalculateSMN(clsX, clsY,SMN);
180	if( det == 1)
181	{
182	if(SMN >= 0 && SMN <= 5) {
183	++CPVcluQUAD[0] ;
184	CPVADCQUAD[0] =+ adc ;
185	CPVcelQUAD[0] =+ ncell ;
186	}
187	if(SMN >= 6 && SMN <=11) {
188	++CPVcluQUAD[1] ;
189	CPVADCQUAD[1] =+ adc ;
190	CPVcelQUAD[1] =+ ncell ;
191	}
192	if(SMN >=12 && SMN <=17) {
193	++CPVcluQUAD[2] ;
194	CPVADCQUAD[2] =+ adc ;
195	CPVcelQUAD[2] =+ ncell ;
196	}
197	if(SMN >=18 && SMN <=23) {
198	++CPVcluQUAD[3] ;
199	CPVADCQUAD[3] =+ adc ;
200	CPVcelQUAD[3] =+ ncell ;
201	}
202
203	if(eta >= 2.3 && eta <= 3.5)
204	{
205	hC2->Fill(phi);
206	}
207	}
208	if( det == 0)
209	{
210	if(SMN >= 0 && SMN <= 5) {
211	++PREcluQUAD[0] ;
212	PREadcQUAD[0] =+ adc ;
213	PREcelQUAD[0] =+ ncell ;
214	}
215	if(SMN >= 6 && SMN <=11) {
216	++PREcluQUAD[1] ;
217	PREadcQUAD[1] =+ adc ;
218	PREcelQUAD[1] =+ ncell ;
219	}
220	if(SMN >=12 && SMN <=17) {
221	++PREcluQUAD[2] ;
222	PREadcQUAD[2] =+ adc ;
223	PREcelQUAD[2] =+ ncell ;
224	}
225	if(SMN >=18 && SMN <=23) {
226	++PREcluQUAD[3] ;
227	PREadcQUAD[3] =+ adc ;
228	PREcelQUAD[3] =+ ncell ;
229	}
230	if ( n <= 100 ) {
231	hSM->Fill(SMN);
232	if(SMN == 0)hSM1->Fill(-clsX,clsY);
233	if(SMN == 0)hSM1->Fill(-clsX,clsY);
234	if(SMN == 1)hSM2->Fill(-clsX,clsY);
235	if(SMN == 2)hSM3->Fill(-clsX,clsY);
236	if(SMN == 3)hSM4->Fill(-clsX,clsY);
237	if(SMN == 4)hSM5->Fill(-clsX,clsY);
238	if(SMN == 5)hSM6->Fill(-clsX,clsY);
239	if(SMN == 6)hSM7->Fill(-clsX,clsY);
240	if(SMN == 7)hSM8->Fill(-clsX,clsY);
241	if(SMN == 8)hSM9->Fill(-clsX,clsY);
242	if(SMN == 9)hSM10->Fill(-clsX,clsY);
243	if(SMN ==10)hSM11->Fill(-clsX,clsY);
244	if(SMN ==11)hSM12->Fill(-clsX,clsY);
245	if(SMN ==12)hSM13->Fill(-clsX,clsY);
246	if(SMN ==13)hSM14->Fill(-clsX,clsY);
247	if(SMN ==14)hSM15->Fill(-clsX,clsY);
248	if(SMN ==15)hSM16->Fill(-clsX,clsY);
249	if(SMN ==16)hSM17->Fill(-clsX,clsY);
250	if(SMN ==17)hSM18->Fill(-clsX,clsY);
251	if(SMN ==18)hSM19->Fill(-clsX,clsY);
252	if(SMN ==19)hSM20->Fill(-clsX,clsY);
253	if(SMN ==20)hSM21->Fill(-clsX,clsY);
254	if(SMN ==21)hSM22->Fill(-clsX,clsY);
255	if(SMN ==22)hSM23->Fill(-clsX,clsY);
256	if(SMN ==23)hSM24->Fill(-clsX,clsY);
257	}
258	if(eta >= 2.3 && eta <= 3.5)
259	{
260	hP2->Fill(phi);
261	}
262	hP1->Fill(clsX,clsY);
263	}
264	}
265	for (Int_t k=0; k<4;k++) {
266	TotCPVClus =+ CPVcluQUAD [k] ;
267	TotPREClus =+ PREcluQUAD [k] ;
268	TotCPVCell =+ CPVcelQUAD [k] ;
269	TotPRECell =+ PREcelQUAD [k] ;
270	TotPREEdep =+ PREadcQUAD [k] ;
271	}
272	Float_t TotCPVPREClus = TotPREClus + TotCPVClus ;
273	Float_t TotCPVPRECell = TotPRECell + TotCPVCell ;
274	if(eta >= 2.3 && eta <= 3.5)
275	{
276	hC3->Fill(TotCPVClus);
277	hP3->Fill(TotPREClus);
278	hP4->Fill(TotPREEdep);
279	hP5->Fill(TotPRECell);
280	hCP0->Fill(PREcluQUAD[0],PREcluQUAD[1]);
281	hCP1->Fill(PREcluQUAD[2],PREcluQUAD[3]);
282	hCP2->Fill(PREadcQUAD[2],PREadcQUAD[3]);
283	hCP3->Fill(TotPREEdep,TotCPVPREClus);
284	hCP4->Fill(TotPREClus,TotCPVClus);
285	}
286	TotCPVClus == 0.0;
287	TotPREClus == 0.0;
288	TotCPVCell == 0.0;
289	TotPRECell == 0.0;
290	TotPREEdep == 0.0;
291	if(n >= nevt)break;
292	}
293
294	gStyle->SetOptStat(110000);
295	gStyle->SetOptFit(1);
296
297	// if ( PRINTSMN == 'y' ) {
298	TCanvas *cP1 = new TCanvas("cP1","Cluster XY", 10,10, 600, 600);
299	cP1->Range(-100, -100,100 ,100 );
300	hSM1->SetMarkerColor(2);
301	hSM1->Draw();
302	hSM1->GetXaxis()->SetTitle("Cluster X");
303	hSM1->GetYaxis()->SetTitle("Cluster Y");
304	hSM2->SetMarkerColor(2);
305	hSM2->Draw("same");
306	hSM3->SetMarkerColor(2);
307	hSM3->Draw("same");
308	hSM4->SetMarkerColor(2);
309	hSM4->Draw("same");
310	hSM5->SetMarkerColor(2);
311	hSM5->Draw("same");
312	hSM6->SetMarkerColor(2);
313	hSM6->Draw("same");
314	hSM7->SetMarkerColor(4);
315	hSM7->Draw("same");
316	hSM8->SetMarkerColor(4);
317	hSM8->Draw("same");
318	hSM9->SetMarkerColor(4);
319	hSM9->Draw("same");
320	hSM10->SetMarkerColor(4);
321	hSM10->Draw("same");
322	hSM11->SetMarkerColor(4);
323	hSM11->Draw("same");
324	hSM12->SetMarkerColor(4);
325	hSM12->Draw("same");
326	hSM13->SetMarkerColor(6);
327	hSM13->Draw("same");
328	hSM14->SetMarkerColor(6);
329	hSM14->Draw("same");
330	hSM15->SetMarkerColor(6);
331	hSM15->Draw("same");
332	hSM16->SetMarkerColor(6);
333	hSM16->Draw("same");
334	hSM17->SetMarkerColor(6);
335	hSM17->Draw("same");
336	hSM18->SetMarkerColor(6);
337	hSM18->Draw("same");
338	hSM19->SetMarkerColor(8);
339	hSM19->Draw("same");
340	hSM20->SetMarkerColor(8);
341	hSM20->Draw("same");
342	hSM21->SetMarkerColor(8);
343	hSM21->Draw("same");
344	hSM22->SetMarkerColor(8);
345	hSM22->Draw("same");
346	hSM23->SetMarkerColor(8);
347	hSM23->Draw("same");
348	hSM24->SetMarkerColor(8);
349	hSM24->Draw("same");
350
351	Int_t linsav = gStyle->GetLineWidth();
352	DrawPMDBoundary();
353	DrawPMDBoundarySM1();
354	DrawPMDBoundarySM2();
355	DrawPMDBoundarySM3();
356	DrawPMDBoundarySM4();
357	cP1->Print("ClusterXY.gif");
358
359
360	TCanvas *cP2 = new TCanvas("PHI CPV / PMD","",10,10, 600,600);
361	cP2->Divide(1,2);
362	cP2->cd(1);
363	cP2->SetFillColor(0);
364	hC2->SetLineColor(4);
365	hC2->Draw();
366	cP2->cd(2);
367	hP2->SetLineColor(2);
368	hP2->Draw();
369	cP2->Print("CPVPREphi.gif");
370
371	// }
372
373	TCanvas *cP2 = new TCanvas("cP2","",10,10,600,600);
374	cP2->cd();
375	hSM->SetFillColor(2);
376	hSM->Draw();
377	cP2->Print("AllSMN.gif");
378
379
380
381	TCanvas *cP3 = new TCanvas("CPV Clus PRE Clus Correlations", " ",10,10, 600,600);
382	cP3->Divide(2,2);
383	cP3->cd(1);
384	hCP0->SetMarkerColor(9);
385	hCP0->Draw();
386	cP3->cd(2);
387	hCP1->SetMarkerColor(6);
388	hCP1->Draw();
389	cP3->cd(3);
390	hP3->SetLineColor(2);
391	hP3->Draw();
392	cP3->cd(4);
393	hCP4->SetMarkerColor(3);
394	hCP4->Draw();
395	cP3->Print("CPVPREClus.gif");
396
397	TCanvas *cP6 = new TCanvas("CPV Clus PRE Adc"," ",10,10, 600,600);
398	cP6->Divide(1,2);
399	cP6->cd(1);
400	hC3->SetLineColor(4);
401	hC3->Draw();
402	cP6->cd(2);
403	hP4->SetLineColor(2);
404	hP4->Draw();
405	cP6->Print("CPVPREAdc.gif");
406
407	timer.Stop();
408	timer.Print();
409	return 0;
410	}
411
412	void CalculateSMN( Float_t clsX, Float_t clsY, Int_t & smn)
413	{
414	//smn = 0;
415
416	//---------------------------------------------------------------------
417	if((clsX <= Xcon[0]) && (clsX >= Xcon[1]) &&
418	(clsY <= Ycon[0]) && (clsY >= Ycon[6])) smn = 0 ;
419
420	else if((clsX <=Xcon[2]) && (clsX >= Xcon[3]) &&
421	(clsY <= Ycon[1]) && (clsY >= Ycon[7]))smn = 1 ;
422
423	else if((clsX <=Xcon[4]) && (clsX >= Xcon[5]) &&
424	(clsY <= Ycon[3]) && (clsY >= Ycon[8]))smn = 2 ;
425
426	else if((clsX <= Xcon[0]) && (clsX >= Xcon[1]) &&
427	(clsY <= Ycon[12]) && (clsY >= Ycon[18])) smn = 3 ;
428
429	else if((clsX <=Xcon[2]) && (clsX >= Xcon[3]) &&
430	(clsY <= Ycon[12]) && (clsY >= Ycon[18]))smn = 4 ;
431
432	else if((clsX <=Xcon[4]) && (clsX >= Xcon[5]) &&
433	(clsY <= Ycon[12]) && (clsY >= Ycon[18]))smn = 5 ;
434	//------------------------------------------------------------------
435	else if((clsX >= Xcon[24]) && (clsX <= Xcon[25]) &&
436	(clsY >= Ycon[24]) && (clsY <= Ycon[30])) smn = 6 ;
437
438	else if((clsX >=Xcon[26]) && (clsX <= Xcon[27]) &&
439	(clsY >= Ycon[25]) && (clsY <= Ycon[31]))smn = 7 ;
440
441	else if((clsX >=Xcon[28]) && (clsX <= Xcon[29]) &&
442	(clsY >= Ycon[26]) && (clsY <= Ycon[32]))smn = 8 ;
443
444	else if((clsX >= Xcon[24]) && (clsX <= Xcon[25]) &&
445	(clsY >= Ycon[36]) && (clsY <= Ycon[42])) smn = 9 ;
446
447	else if((clsX >=Xcon[26]) && (clsX <= Xcon[27]) &&
448	(clsY >= Ycon[36]) && (clsY <= Ycon[42]))smn = 10;
449
450	else if((clsX >=Xcon[28]) && (clsX <= Xcon[29]) &&
451	(clsY >= Ycon[36]) && (clsY <= Ycon[42]))smn = 11;
452	//------------------------------------------------------------------
453	else if((clsX <= Xcon[48]) && (clsX >= Xcon[49]) &&
454	(clsY <= Ycon[48]) && (clsY >= Ycon[52])) smn = 12 ;
455
456	else if((clsX <=Xcon[50]) && (clsX >= Xcon[51]) &&
457	(clsY <= Ycon[48]) && (clsY >= Ycon[52]))smn = 13 ;
458
459	else if((clsX <=Xcon[48]) && (clsX >= Xcon[49]) &&
460	(clsY <= Ycon[56]) && (clsY >= Ycon[60]))smn = 14 ;
461
462	else if((clsX <=Xcon[50]) && (clsX >= Xcon[51]) &&
463	(clsY <= Ycon[56]) && (clsY >= Ycon[60]))smn = 15 ;
464
465	else if((clsX <=Xcon[48]) && (clsX >= Xcon[49]) &&
466	(clsY <= Ycon[64]) && (clsY >= Ycon[68]))smn = 16 ;
467
468	else if((clsX <=Xcon[50]) && (clsX >= Xcon[51]) &&
469	(clsY <= Ycon[64]) && (clsY >= Ycon[68]))smn = 17 ;
470	//--------------------------------------------------------------
471	else if((clsX >= Xcon[72]) && (clsX <= Xcon[73]) &&
472	(clsY >= Ycon[72]) && (clsY <= Ycon[76])) smn = 18 ;
473
474	else if((clsX >=Xcon[74]) && (clsX <= Xcon[75]) &&
475	(clsY >= Ycon[72]) && (clsY <= Ycon[76]))smn = 19 ;
476
477	else if((clsX >=Xcon[72]) && (clsX <= Xcon[73]) &&
478	(clsY >= Ycon[80]) && (clsY <= Ycon[84]))smn = 20 ;
479
480	else if((clsX >=Xcon[74]) && (clsX <= Xcon[75]) &&
481	(clsY >= Ycon[80]) && (clsY <= Ycon[84]))smn = 21;
482
483	else if((clsX >= Xcon[72]) && (clsX <= Xcon[73]) &&
484	(clsY >= Ycon[88]) && (clsY <= Ycon[92])) smn = 22 ;
485
486	else if((clsX >=Xcon[74]) && (clsX <= Xcon[75]) &&
487	(clsY >= Ycon[88]) && (clsY <= Ycon[92]))smn = 23 ;
488	else smn = 111;
489
490	}
491
492	void DrawPMDBoundary()
493	{
494	TH2F *h = new TH2F("h","",200,-100,100,200,-100,100);
495	gStyle->SetLineWidth(2);
496	gStyle->SetLineColor(2);
497	TLine * l;
498	l = new TLine(75.1333, 86.475, -75.1333, 86.475); l->Draw("same");
499	l = new TLine(-75.1333, 86.470,-75.1333, -86.475); l->Draw("same");
500	l = new TLine(-75.1333, -86.475,75.1333, -86.475); l->Draw("same");
501	l = new TLine(75.1333, -86.475,75.1333, 86.475); l->Draw("same");
502	}
503	void DrawPMDBoundarySM1()
504	{
505	TH2F *hsm1 = new TH2F("hsm1","",200,-100,100,200,-100,100);
506	gStyle->SetLineWidth(1);
507	gStyle->SetLineColor(4);
508	TLine * l;
509	l = new TLine(-75.1333, 86.475, -10.447, 86.475); l->Draw("same");
510	l = new TLine(-10.447, 86.475, -10.446, -10.925); l->Draw("same");
511	l = new TLine(-10.446, -10.925, -75.1333,-10.925); l->Draw("same");
512	l = new TLine(-75.1333,-10.925, -75.1333, 86.475); l->Draw("same");
513	}
514	void DrawPMDBoundarySM2()
515	{
516	TH2F *hsm2 = new TH2F("hsm2","",200,-100,100,200,-100,100);
517	gStyle->SetLineWidth(1);
518	gStyle->SetLineColor(4);
519	TLine * l;
520	l = new TLine(75.1333, -86.475, 10.446, -86.475); l->Draw("same");
521	l = new TLine(10.446, -86.475, 10.446, 10.925); l->Draw("same");
522	l = new TLine(10.446, 10.925, 75.1333, 10.925); l->Draw("same");
523	l = new TLine(75.1333, 10.925, 75.1333, -86.475); l->Draw("same");
524	}
525
526	void DrawPMDBoundarySM3()
527	{
528	TH2F *hsm3 = new TH2F("hsm3","",200,-100,100,200,-100,100);
529	gStyle->SetLineWidth(1);
530	gStyle->SetLineColor(1);
531	TLine * l;
532	l = new TLine( -9.167, 86.475, 75.1333, 86.475); l->Draw("same");
533	l = new TLine(75.1333,86.475, 75.1333, 11.925); l->Draw("same");
534	l = new TLine(75.1333,11.925, -9.167, 11.925); l->Draw("same");
535	l = new TLine( -9.167, 11.925, -9.167, 86.475); l->Draw("same");
536	}
537	void DrawPMDBoundarySM4()
538	{
539	TH2F *hsm4 = new TH2F("hsm4","",200,-100,100,200,-100,100);
540	gStyle->SetLineWidth(1);
541	gStyle->SetLineColor(1);
542	TLine * l;
543	l = new TLine(9.167, -86.475, -75.1333,-86.475); l->Draw("same");
544	l = new TLine(-75.1333,-86.475, -75.1333,-11.925); l->Draw("same");
545	l = new TLine(-75.1333,-11.925, 9.167, -11.925); l->Draw("same");
546	l = new TLine(9.167, -11.925, 9.167, -86.475); l->Draw("same");
547	}
548