//-----------------------------------------------------------------------------------------
//---            MACRO PER IL CALCOLO DELL'EFFICIENZA DEI CONTATORI DEL CSX
//---            USANDO low PT LEPTONS da J/Psi o Yupsilon
//---
//---
//---
//---
//-----------------------------------------------------------------------------------------
//
{
// 
gROOT->Reset();
//
//
//-----------------
   int laptop = 0;
   int cuts_on= 1; 
//-----------------
//
//
//
// ================================================================
// ===			    DATA SAMPLES			===
// ================================================================
//
  if (laptop == 1)
    {
//
//===  ROOT ON LAPOTOP  ===
//
TFile RF_h110("c:/ROOT/MCSamples/h110.root");
TTree *h110 = (TTree*)RF_h110 ->Get("ReadHmt/Hmt");
//
TFile RF_h120("c:/ROOT/MCSamples/h120.root");
TTree *h120 = (TTree*)RF_h120 ->Get("ReadHmt/Hmt");
//
TFile RF_h130("c:/ROOT/MCSamples/h130.root");
TTree *h130 = (TTree*)RF_h130 ->Get("ReadHmt/Hmt");
//
//  TFile RF_h140("c:/ROOT/MCSamples/h140.root");
//  TTree *h140 = (TTree*)h140 ->Get("ReadHmt/Hmt");
//
//  TFile RFh150("c:/ROOT/MCSamples/h150.root");
//  TTree *h150 = (TTree*)h150 ->Get("ReadHmt/Hmt");
//
TFile RF_z("c:/ROOT/MCSamples/ztt.root");
TTree *z = (TTree*)RF_z ->Get("ReadHmt/Hmt");
//
TFile RF_ww("c:/ROOT/MCSamples/ww.root");
TTree *ww = (TTree*)RF_ww ->Get("ReadHmt/Hmt");
}
//
  else 
{
//===  ROOT ON POOL
//
//TFile RF_h110("/home/users/pagliar/CMEX/DATA/JPSI.root");
//TTree *h110 = (TTree*)RF_h110 ->Get("csx");
//
TChain *h110 = new TChain("csx");
h110->Add("*.root");
h110->GetEntry(0);
//
}
//
//==================================================================
//
//
//===    HERE I PRINT INFOs on THE NTUPLEs    ===
//
h110->Print();
//
  if (cuts_on == 1)
{
//
//===================================================================
//===== 		HERE I DEFINE MY CUTS			=====
//===================================================================
//
// Esistenza di almeno 2 muoni
// ---------------------------
//
TCut cta0 = "nMUO>1";
//
// Z Vertex dell'Evento
// --------------------
//
TCut cta1 = "abs(Zvtx)<60.";
//
//
// Impact Parameter
// ----------------
//
TCut cta2a = "abs(MUO.trk.D0[0])<0.6";
TCut cta2b = "abs(MUO.trk.D0[1])<0.6";
TCut cta2  = cta2a && cta2b;
//
//
//
// Muon PT
// -------
//
TCut cta3a = "MUO.Pt[0]>=  1.5";
TCut cta3b = "";
//TCut cta3b = "MUO.Pt[1]>=  1.5";
TCut  cta3 = cta3a && cta3b;
//
//
//
// COT-trk vs Stub Match
// ---------------------
//
TCut cta4a = "abs(MUO.dx[0]) < 12.";
TCut cta4b = "";
TCut cta4b = "abs(MUO.dx[1]) < 12.";
TCut  cta4 = cta4a && cta4b;
//TCut  cta4 ="";
//
//
//
// Muon Isolation
// --------------
//
TCut cta5a = "MUO.R4Et[0]<2.0"; 
TCut cta5b = "MUO.R4Et[1]<2.0"; 
TCut  cta5 = cta5a && cta5b; 
 TCut  cta5 ="";
//
//
//
// MIP conditions
// --------------
//
TCut cta6a = "MUO.hdEt[0]<6.";
TCut cta6b = "MUO.hdEt[1]<6.";
TCut cta6  = cta6a && cta6b;
//
TCut cta7a = "MUO.emEt[0]<2.";
TCut cta7b = "MUO.emEt[1]<2.";
TCut cta7  = cta7a && cta7b;
//
//
//== Distance between the two Muon vertex
//
//TCut cta8 = """;
TCut cta8  = "abs(MUO.trk.Z0[0]-MUO.trk.Z0[1]) <5.";
//
//
//== Cosmic Rays Removal
TCut cta9 = "";
//
TCut CUTa = cta0&&cta1&&cta2&&cta3&&cta4&&cta5&&cta6&&cta7&&cta8;
//
//==================================
// J/Psi, Yupsilon & Z selection:
//
TCut ctb1 = "MUO.trk.Charge[0]*MUO.trk.Charge[1]<0.";
// Cerri CDFNote 5838 [2.7,3.6]
//  3096 MeV +/- 200 MeV di mass window
// 
TCut ctb2 = "mZ0mu>2.896";
TCut ctb3 = "mZ0mu<3.296";
//
/////////////////////////////////////////////MUO.nt1 == 1";

////////////////////////////////////////////////////////////////////////////////////////////////////
////
////                                         BUD RUNs
////
////////////////////////////////////////////////////////////////////////////////////////////////////
TCut  R1 = "Run!=144720.&&Run!=143891.&&Run!=143827.&&Run!=143720.&&Run!=143719.";
TCut  R2 = "Run!=143718.&&Run!=143717.&&Run!=143705.&&Run!=143703.&&Run!=143637.";
TCut  R3 = "Run!=143634.&&Run!=143633.&&Run!=143606.&&Run!=143593.&&Run!=143592.";
TCut  R4 = "Run!=143591.&&Run!=143285.&&Run!=143257.&&Run!=143254.&&Run!=143250.";
TCut  R5 = "Run!=143247.&&Run!=143245.&&Run!=143244.&&Run!=142753.&&Run!=142718.";
//
TCut  R6 = "Run!=142670.&&Run!=142667.&&Run!=142648.&&Run!=142235.&&Run!=142234.";
TCut  R7 = "Run!=142233.&&Run!=142226.&&Run!=142228.&&Run!=142105.&&Run!=141986.";
TCut  R8 = "Run!=141507.&&Run!=141434.&&Run!=140853.&&Run!=140744.&&Run!=140743.";
TCut  R9 = "Run!=140726.&&Run!=139801.&&Run!=139335.&&Run!=139296.&&Run!=139249.";
TCut R10 = "Run!=139177.&&Run!=138815.&&Run!=138809.";
////////////////////////////////////////////////////////////////////////////////////////////////////
 TCut RRR = R1&&R2&&R3&&R4&&R5&&R6&&R7&&R8&&R9&&R10;
////////////////////////////////////////////////////////////////////////////////////////////////////
// == CONDIZIONE SUL NUMERO DI RUN PER TENERE CONTO DEGLI ERRORI SU I PRIMI 12 CONTATORI
//
 //TCut ctb4 = "Run > 143500.";
 //  
TCut ctb4 = "";
TCut ctb5 = "";
TCut ctb6 = "";
TCut ctb7 = "";
TCut ctb8 = "";
TCut ctb9 = "";
//
//
//
 TCut CUTb = ctb1&&ctb2&&ctb3&&ctb4&&ctb5&&ctb6&&ctb7&&ctb8&&ctb9;
//
//
//
 //TCut CUT  = CUTa&&CUTb&&RRR;
TCut CUT  = CUTa&&CUTb;
//
}
else TCut CUT ="";
//
// OUTSIDE DETECTORS
//---------------------------
TCut T1d = "MUO.T1[0]>100.";
TCut T1u = "MUO.T1[0]<400.";
TCut T1  = T1d&&T1u;
TCut T1c = "MUO.T1[0]>400.";
//
// INSIDE DETECTORS
//---------------------------
TCut T0d = "MUO.T0[0]>100.";
TCut T0u = "MUO.T0[0]<400.";
TCut T0  = T0d&&T0u;
TCut T0c = "MUO.T0[0]>400.";
//--------------------------
// 
TCut  EST = "MUO.Eta[0] > 0.";
TCut WEST = "MUO.Eta[0] < 0.";
//
//
//TCut C012   = "MUO.C0[0] > 12.";
//TCut C112   = "MUO.C1[0] > 12.";
//
//
TCut INTEST = CUT&&EST&&T0;
TCut EXTEST = CUT&&EST&&T1;
TCut ALLEST = CUT&&EST;
//
TCut EST0C   = CUT&&EST&&T0c;
TCut EST1C   = CUT&&EST&&T1c;
//
TCut INTWEST = CUT&&WEST&&T0;
TCut EXTWEST = CUT&&WEST&&T1;
TCut ALLWEST = CUT&&WEST;
//
TCut WEST0C   = CUT&&WEST&&T0c;
TCut WEST1C   = CUT&&WEST&&T1c;
//
//===================================================================
//						     END OF CUTsDEF =
//=================================================================== 
//
//
//
//
//
//
//
//
//
//
//
//
//
//
//
char text[12];
char *symbol[];
//
//   
c1 = new TCanvas("mupt","CSX Est Counter Situation",10,40,800,600);
c1->Range(0,0,25,18);
c1->SetFillColor(40);
//
//
//TPaveLabel *pl = new TPaveLabel(1,16.1,24,17.3,"CSX Est Counter Efficiencies","br");
TPaveLabel *pl = new TPaveLabel(1,16.1,24,17.3,"CSX Est Counter Situation","br");
pl->SetFillColor(18);
pl->SetTextFont(32);
pl->SetTextColor(49);
pl->Draw();
//
TText *t = new TText();
t->SetTextFont(32);
t->SetTextColor(1);
t->SetTextSize(0.03);
t->SetTextAlign(12);
t->DrawText(4.3,15.5,"CSX Internal Est");
t->DrawText(16.9,15.5,"CSX External Est");
//
//
// Qui creo e disegno 2 PAD
//
//                                                            V- colore
   pad1 = new TPad("pad1","This is pad1",0.02,0.02,0.48,0.83,33);
   pad2 = new TPad("pad2","This is pad2",0.52,0.02,0.98,0.83,33);
   pad1->Draw();
   pad2->Draw();
//
//
// Background
//

//
//  (1) bottom/left; (2)  bottom/right; (3) top/right
//
//gStyle->SetOptDate(3);
//
// Change default style for the statistics box
//
gStyle->SetOptStat(1111110);
gStyle->SetOptStat(0000000);
//
gStyle->SetStatW(0.19);
gStyle->SetStatH(0.19);
// gStyle->SetStatFont(70)
gStyle->SetStatColor(42);
gStyle->SetTitleColor(42);
// gStyle->SetTitleTextColor(45);
// pad1->SetGrid();
// pad1->SetLogy();
c1->GetFrame()->SetFillColor(45);
//
  c1->Modified();
  c1->Update();
//
//
//================================================================
//====         CANVAS n.1    EST  COUNTERS                    ====
//================================================================
//
// hint0->Delete();
// 
// 2= Rosso  -  3= GREEN  -  4= BLUE - 5= YELLOW 
//
//
//
//EMPTY HISTO
//hpx2 = new TH2F ("hpx","",14,0.,14.,100,0.,1.0);
hpx2 = new TH2F ("hpx","",60,14.,74.,100,0.,1.0);
//
//== EST
//
TH1F *H_goodIE  = new TH1F("H_goodIE","",75,0.,75.);
TH1F *H_goodEE  = new TH1F("H_goodEE","",75,0.,75.);
//
TH1F   *H_allIE  = new TH1F("H_allIE","",75,0.,75.);
TH1F   *H_allEE  = new TH1F("H_allEE","",75,0.,75.);
//
   pad1->cd();
   pad1->Range(-0.255174,-19.25,2.29657,-6.75);
//
//== Scala X ed Y su entrambi i lati del plot
//
   pad1->SetTickx();
   pad1->SetTicky();
//
//== Scala LOGARITMICA
//
//   pad1->SetLogx();
//   pad1->SetLogy();
     pad1->SetGrid();
//
//
//  GOOD CSX INTERNAL EST
//
//H_goodIE->SetNormFactor();
H_goodIE->SetLineWidth(5);
H_goodIE->SetLineColor(46);
H_goodIE->SetFillColor(46);
h110->Draw("MUO.C0[0]>>H_goodIE",INTEST);
//
//
//  ALL CSX ITERNAL
//
//H_allIE->SetNormFactor();
H_allIE->SetLineWidth(5);
H_allIE->SetLineColor(38);
H_allIE->SetFillColor(0);
h110->Draw("MUO.C0[0]>>H_allIE",ALLEST);
//
//
/// CALCOLO L'EFFICIENZA = CONTATORI INTERNAL EST
//
//
TH1F *EFF0 = (TH1F*)H_goodIE->Clone();
EFF0->SetLineWidth(1);
EFF0->SetLineColor(4);
EFF0->SetFillColor(0);
EFF0->SetMarkerStyle(20);
EFF0->SetMarkerColor(4);
//
//
//
//
//H_allIE->Sumw2();
//H_goodIE->Sumw2();
//
//EFF0->Divide(H_goodIE,H_allIE);
//
//hpx2->Draw();
//EFF0->Draw("same");
//
//=== REDRAW AXIS
//
gPad->RedrawAxis();
//
//
//
//getchar();
//
TH1F *H_goodIEc  = new TH1F("H_goodIEc","",75,0.,75.);
h110->Draw("MUO.C0[0]>>H_goodIEc",EST0C);
TH1F *H_allIEc = (TH1F*)H_allIE->Clone();
H_allIEc->Add(H_allIE,H_goodIEc,1.0,-1.0);
//
H_allIEc->Sumw2();
H_goodIE->Sumw2();
//
EFF0->Divide(H_goodIE,H_allIEc);
//
hpx2->Draw();
EFF0->Draw("same");
//
//
//
//------------------------------------------------------------------------------
//
   pad2->cd();
   pad2->Range(-0.43642,-23.75,3.92778,-6.25);
//
//== Scala X ed Y su entrambi i lati del plot
//
   pad2->SetTickx();
   pad2->SetTicky();
//
//== Scala LOGARITMICA
//
//   pad2->SetLogx();
//   pad2->SetLogy();
     pad2->SetGrid();
//
//H_goodEE->SetNormFactor();
H_goodEE->SetLineWidth(5);
H_goodEE->SetLineColor(46);
H_goodEE->SetFillColor(46);
h110->Draw("MUO.C1[0]>>H_goodEE",EXTEST);
//
//  ALL CSX ITERNAL
//
//H_allEE->SetNormFactor();
H_allEE->SetLineWidth(5);
H_allEE->SetLineColor(38);
H_allEE->SetFillColor(0);

h110->Draw("MUO.C1[0]>>H_allEE",ALLEST);
//
//
TH1F *EFF1 = (TH1F*)H_goodEE->Clone();
EFF1->SetLineWidth(1);
EFF1->SetLineColor(46);
EFF1->SetFillColor(0);
EFF1->SetMarkerStyle(20);
EFF1->SetMarkerColor(46);
//
//H_allEE->Sumw2();
//H_goodEE->Sumw2();
//
//EFF1->Divide(H_goodEE,H_allEE);
//
TH1F *H_goodEEc  = new TH1F("H_goodEEc","",75,0.,75.);
h110->Draw("MUO.C1[0]>>H_goodEEc",EST1C);
TH1F *H_allEEc = (TH1F*)H_allEE->Clone();
H_allEEc->Add(H_allEE,H_goodEEc,1.0,-1.0);
//
H_allEEc->Sumw2();
H_goodEE->Sumw2();
//
EFF1->Divide(H_goodEE,H_allEEc);
//
hpx2->Draw();
EFF1->Draw("same");
//
gPad->RedrawAxis();
//
//
//
//
 if(laptop == 1)
     c1->Print("C:/RootAna/HMT/JPSI_eff-EST.ps");
 else
     c1->Print("Eff_EST-13-73.ps");
// 
EFF0->Delete();
EFF1->Delete();
//
//================================================================================
//===
//=== CANVAS n. 2  -  EST COUNTERS - Details
//===
//================================================================================
//
//
//-----------------------------------------------------------------
   c2 = new TCanvas("c2","Est Counters Details",200,10,700,780);
   c2->GetFrame()->SetFillColor(45);
   //
   pad1 = new TPad("pad1","This is pad1",0.02,0.52,0.48,0.98,21);
   pad2 = new TPad("pad2","This is pad2",0.52,0.52,0.98,0.98,21);
   pad3 = new TPad("pad3","This is pad3",0.02,0.02,0.48,0.48,21);
   pad4 = new TPad("pad4","This is pad4",0.52,0.02,0.98,0.48,21);
   pad1->Draw();
   pad2->Draw();
   pad3->Draw();
   pad4->Draw();
   //
   //
   gStyle->SetStatW(0.30);
   gStyle->SetStatH(0.20);
   gStyle->SetOptStat(1111110);
   gStyle->SetStatColor(42);
//
//hpx0 = new TH2F ("hpx","",200,-0.1,0.1,100,0.,2000);
//hpx2 = new TH2F ("hpx","",200,-0.1,0.1,100,0.,2000);
//
//------------------------------
//---   PAD 1
//------------------------------
   pad1->cd();
   pad1->SetTickx();
   pad1->SetTicky();
   pad1->SetGrid();
   //pad1->SetLogy();
   //pad1->SetLogx();
   pad1->GetFrame()->SetFillColor(15);
//   
//H_goodIE->SetNormFactor(15000);
H_goodIE->SetLineWidth(3);
H_goodIE->SetLineColor(4);
H_goodIE->SetFillColor(4);
H_goodIE->Draw("");
//
//------------------------------
//---   PAD 2
//------------------------------
   pad2->cd();
   pad2->SetTickx();
   pad2->SetTicky();
   pad2->SetGrid();
   //pad2->SetLogy();
   //pad2->SetLogx();
   pad2->GetFrame()->SetFillColor(15);
//   
//H_goodEE->SetNormFactor(15000);
H_goodEE->SetLineWidth(3);
H_goodEE->SetLineColor(46);
H_goodEE->SetFillColor(46);
H_goodEE->Draw("");
//
//------------------------------
//---   PAD 3
//------------------------------
   pad3->cd();
   pad3->SetTickx();
   pad3->SetTicky();
   pad3->SetGrid();
   //pad3->SetLogy();
   //pad3->SetLogx();
   pad3->GetFrame()->SetFillColor(15);
//   
//
// ALL EST
//
//H_allIE->SetNormFactor(15000);
H_allIEc->SetLineWidth(1);
H_allIEc->SetLineColor(5);
H_allIEc->SetFillColor(5);
H_allIEc->Draw("");
// 
//------------------------------
//---   PAD 4
//------------------------------
   pad4->cd();
   pad4->SetTickx();
   pad4->SetTicky();
   pad4->SetGrid();
   //pad4->SetLogy();
   //pad4->SetLogx();
   pad4->GetFrame()->SetFillColor(15);
//   
// 
//
//H_allEE->SetNormFactor(15000);
H_allEEc->SetLineWidth(1);
H_allEEc->SetLineColor(38);
H_allEEc->SetFillColor(38);
H_allEEc->Draw("");
// 
 if(laptop == 1)
     c2->Print("C:/RootAna/HMT/JPSI_EST.ps");
 else
     c2->Print("EST-13-73.ps");
//
//
getchar();
c1->Close();
c2->Close();
//
//
//===================================================================
//					         END OF CANVAS n. 2 =
//===================================================================
// 
//
//
//===================================================================
//====                          CSX  WEST                        ====
//===================================================================
//
//
//
char text[12];
char *symbol[];
//

c3 = new TCanvas("mupt","CSX West Counter Situation",10,40,800,600);
c3->Range(0,0,25,18);
c3->SetFillColor(40);
//
//
TPaveLabel *pl = new TPaveLabel(1,16.1,24,17.3,"CSX West Counter Situation","br");
pl->SetFillColor(18);
pl->SetTextFont(32);
pl->SetTextColor(49);
pl->Draw();
//
TText *t = new TText();
t->SetTextFont(32);
t->SetTextColor(1);
t->SetTextSize(0.03);
t->SetTextAlign(12);
t->DrawText(4.3,15.5,"CSX Internal West");
t->DrawText(16.9,15.5,"CSX External West");
//
//
// Qui creo e disegno 2 PAD
//
//                                                            V- colore
   pad1 = new TPad("pad1","This is pad1",0.02,0.02,0.48,0.83,33);
   pad2 = new TPad("pad2","This is pad2",0.52,0.02,0.98,0.83,33);
   pad1->Draw();
   pad2->Draw();
//
//
// Background
//

//
//  (1) bottom/left; (2)  bottom/right; (3) top/right
//
//gStyle->SetOptDate(3);
//
// Change default style for the statistics box
//
gStyle->SetOptStat(1111110);
gStyle->SetOptStat(0000000);
//
gStyle->SetStatW(0.19);
gStyle->SetStatH(0.19);
// gStyle->SetStatFont(70)
gStyle->SetStatColor(42);
gStyle->SetTitleColor(42);
// gStyle->SetTitleTextColor(45);
// pad1->SetGrid();
// pad1->SetLogy();
c3->GetFrame()->SetFillColor(45);
//
  c3->Modified();
  c3->Update();
//
//
//
//================================================================
//====         CANVAS n.3    WEST  COUNTERS                    ====
//================================================================
//
// hint0->Delete();
// 
// 2= Rosso  -  3= GREEN  -  4= BLUE - 5= YELLOW 
//
//
//
//EMPTY HISTO
//hpx2 = new TH2F ("hpx","",75,0.,75.,100,0.,1.0);
//
//== WEST
TH1F *H_goodIW  = new TH1F("H_goodIW","",75,0.,75.);
TH1F *H_goodEW  = new TH1F("H_goodEW","",75,0.,75.);
TH1F   *H_allIW  = new TH1F("H_allIW","",75,0.,75.);
TH1F   *H_allEW  = new TH1F("H_allEW","",75,0.,75.);
//
   pad1->cd();
   pad1->Range(-0.255174,-19.25,2.29657,-6.75);
//
//== Scala X ed Y su entrambi i lati del plot
//
   pad1->SetTickx();
   pad1->SetTicky();
//
//== Scala LOGARITMICA
//
//   pad1->SetLogx();
//   pad1->SetLogy();
     pad1->SetGrid();
//
//
//  GOOD CSX INTERNAL EST
//
//H_goodIW->SetNormFactor();
H_goodIW->SetLineWidth(5);
H_goodIW->SetLineColor(46);
H_goodIW->SetFillColor(46);
h110->Draw("MUO.C0[0]>>H_goodIW",INTWEST);
//
//
//  ALL CSX ITERNAL WEST
//
//H_allIW->SetNormFactor();
H_allIW->SetLineWidth(5);
H_allIW->SetLineColor(38);
H_allIW->SetFillColor(0);
h110->Draw("MUO.C0[0]>>H_allIW",ALLWEST);
//
//
/// CALCOLO L'EFFICIENZA = CONTATORI INTERNAL WEST
//
//
TH1F *EFF0 = (TH1F*)H_goodIE->Clone();
EFF0->SetLineWidth(1);
EFF0->SetLineColor(4);
EFF0->SetFillColor(0);
EFF0->SetMarkerStyle(20);
EFF0->SetMarkerColor(4);
//
//H_allIW->Sumw2();
//H_goodIW->Sumw2();
//
//EFF0->Divide(H_goodIW,H_allIW);
//
TH1F *H_goodIWc  = new TH1F("H_goodIWc","",75,0.,75.);
h110->Draw("MUO.C0[0]>>H_goodIWc",WEST0C);
TH1F *H_allIWc = (TH1F*)H_allIW->Clone();
H_allIWc->Add(H_allIW,H_goodIWc,1.0,-1.0);
//
H_allIWc->Sumw2();
H_goodIW->Sumw2();
//
EFF0->Divide(H_goodIW,H_allIWc);
//
hpx2->Draw();
EFF0->Draw("same");
//
//
//=== REDRAW AXIS
//
gPad->RedrawAxis();
//
//------------------------------------------------------------------------------
//
   pad2->cd();
   pad2->Range(-0.43642,-23.75,3.92778,-6.25);
//
//== Scala X ed Y su entrambi i lati del plot
//
   pad2->SetTickx();
   pad2->SetTicky();
//
//== Scala LOGARITMICA
//
//   pad2->SetLogx();
//   pad2->SetLogy();
     pad2->SetGrid();
//
//H_goodEW->SetNormFactor();
H_goodEW->SetLineWidth(5);
H_goodEW->SetLineColor(46);
H_goodEW->SetFillColor(46);
h110->Draw("MUO.C1[0]>>H_goodEW",EXTWEST);
//
//  ALL CSX ITERNAL
//
//H_allEW->SetNormFactor();
H_allEW->SetLineWidth(5);
H_allEW->SetLineColor(38);
H_allEW->SetFillColor(0);
h110->Draw("MUO.C1[0]>>H_allEW",ALLWEST);
//
//
TH1F *EFF1 = (TH1F*)H_goodEE->Clone();
EFF1->SetLineWidth(1);
EFF1->SetLineColor(46);
EFF1->SetFillColor(0);
EFF1->SetMarkerStyle(20);
EFF1->SetMarkerColor(46);
//
//H_allEW->Sumw2();
//H_goodEW->Sumw2();
//
//EFF1->Divide(H_goodEW,H_allEW);
//
//hpx2->Draw();
//EFF1->Draw("same");
//
//
TH1F *H_goodEWc  = new TH1F("H_goodEWc","",75,0.,75.);
h110->Draw("MUO.C1[0]>>H_goodEWc",WEST1C);
TH1F *H_allEWc = (TH1F*)H_allEW->Clone();
H_allEWc->Add(H_allEW,H_goodEWc,1.0,-1.0);
//
H_allEWc->Sumw2();
H_goodEW->Sumw2();
//
EFF1->Divide(H_goodEW,H_allEWc);
//
hpx2->Draw();
EFF1->Draw("same");
//
//
//
//
gPad->RedrawAxis();
//
//
//
//
 if(laptop == 1)
     c3->Print("C:/RootAna/HMT/JPSI_eff-WEST.ps");
 else
     c3->Print("Eff_WEST-13-73.ps");
// 
//
//================================================================================
//===
//=== CANVAS n. 4  -  WEST COUNTERS - Details
//===
//================================================================================
//
//
//-----------------------------------------------------------------
   c4 = new TCanvas("c4","The Ntuple canvas",200,10,700,780);
   c4->GetFrame()->SetFillColor(45);
   //
   pad1 = new TPad("pad1","This is pad1",0.02,0.52,0.48,0.98,21);
   pad2 = new TPad("pad2","This is pad2",0.52,0.52,0.98,0.98,21);
   pad3 = new TPad("pad3","This is pad3",0.02,0.02,0.48,0.48,21);
   pad4 = new TPad("pad4","This is pad4",0.52,0.02,0.98,0.48,21);
   pad1->Draw();
   pad2->Draw();
   pad3->Draw();
   pad4->Draw();
   //
   //
   gStyle->SetStatW(0.30);
   gStyle->SetStatH(0.20);
   gStyle->SetOptStat(1111110);
   gStyle->SetStatColor(42);
//
//hpx0 = new TH2F ("hpx","",200,-0.1,0.1,100,0.,2000);
//hpx2 = new TH2F ("hpx","",200,-0.1,0.1,100,0.,2000);
//
//------------------------------
//---   PAD 1
//------------------------------
   pad1->cd();
   pad1->SetTickx();
   pad1->SetTicky();
   pad1->SetGrid();
   //pad1->SetLogy();
   //pad1->SetLogx();
   pad1->GetFrame()->SetFillColor(15);
//   
//H_goodIW->SetNormFactor(15000);
H_goodIW->SetLineWidth(3);
H_goodIW->SetLineColor(4);
H_goodIW->SetFillColor(4);
H_goodIW->Draw("");
//
//------------------------------
//---   PAD 2
//------------------------------
   pad2->cd();
   pad2->SetTickx();
   pad2->SetTicky();
   pad2->SetGrid();
   //pad2->SetLogy();
   //pad2->SetLogx();
   pad2->GetFrame()->SetFillColor(15);
//   
//H_goodEW->SetNormFactor(15000);
H_goodEW->SetLineWidth(3);
H_goodEW->SetLineColor(46);
H_goodEW->SetFillColor(46);
H_goodEW->Draw("");
//
//------------------------------
//---   PAD 3
//------------------------------
   pad3->cd();
   pad3->SetTickx();
   pad3->SetTicky();
   pad3->SetGrid();
   //pad3->SetLogy();
   //pad3->SetLogx();
   pad3->GetFrame()->SetFillColor(15);
//   
//
// ALL EST
//
//H_allIWc->SetNormFactor(15000);
H_allIWc->SetLineWidth(1);
H_allIWc->SetLineColor(5);
H_allIWc->SetFillColor(5);
H_allIWc->Draw("");
// 
//------------------------------
//---   PAD 4
//------------------------------
   pad4->cd();
   pad4->SetTickx();
   pad4->SetTicky();
   pad4->SetGrid();
   //pad4->SetLogy();
   //pad4->SetLogx();
   pad4->GetFrame()->SetFillColor(15);
//   
// 
//
//H_allEWc->SetNormFactor(15000);
H_allEWc->SetLineWidth(1);
H_allEWc->SetLineColor(38);
H_allEWc->SetFillColor(38);
H_allEWc->Draw("");
// 
 if(laptop == 1)
     c4->Print("C:/RootAna/HMT/JPSI_WEST.ps");
 else
     c4->Print("WEST-13-73.ps");
//
//
getchar();
c3->Close();
c4->Close();
//
//===================================================================
//					         END OF CANVAS n. 4 =
//=================================================================== 
}