From: alla <alla@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Mon, 18 Apr 2011 07:22:49 +0000 (+0000)
Subject: efficiency histograms added
X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=750c312d30256c4aa6b34c3fd0fdd8a7fbd4e0bb;hp=17ab437aef33eb3fc9b8c771ddc62c0aa1e0b868

efficiency histograms added
---

diff --git a/PWG1/T0/AliT0AnalysisTaskQA.cxx b/PWG1/T0/AliT0AnalysisTaskQA.cxx
index f65560b5b77..8140f5f3654 100644
--- a/PWG1/T0/AliT0AnalysisTaskQA.cxx
+++ b/PWG1/T0/AliT0AnalysisTaskQA.cxx
@@ -73,18 +73,22 @@ void AliT0AnalysisTaskQA::UserCreateOutputObjects()
  fTimeVSAmplitude = new TH2F*[NPMT0];
 
  for (Int_t i=0; i<NPMT0; i++) {
-    fTimeVSAmplitude[i]= new TH2F (Form("fTimeVSAmplitude%d",i+1),"fTimeVsAmplitude",60, -10, 50,500,6000,7000);
+    fTimeVSAmplitude[i]= new TH2F (Form("fTimeVSAmplitude%d",i+1),"fTimeVsAmplitude",60, -10, 50,1500,1000,7000);
   }
 
   fTzeroORAplusORC = new TH1F("fTzeroORAplusORC","ORA+ORC /2",100,-2000,2000);   //or A plus or C 
   fResolution      = new TH1F("fResolution","fResolution",100,-2000,2000);// or A minus or C spectrum
   fTzeroORA        = new TH1F("fTzeroORA","fTzeroORA",100,-2000,2000);// or A spectrum
   fTzeroORC        = new TH1F("fTzeroORC","fTzeroORC",100,-2000,2000);// or C spectrum
-  fCFDVSPmtId      = new TH2F("fCFDVSPmtId","fCFDVSPmtId",24,0,24,250,6000,7000);  // 
+  fCFDVSPmtId      = new TH2F("fCFDVSPmtId","fCFDVSPmtId",24,0,24,1500,1000,7000);  // 
   fSPDVertexVST0Vertex = new TH2F("fSPDVertexVST0Vertex","fSPDVertexVST0Vertex",30,-30,30,30,-30,30);
   fOrAvsNtracks    = new TH2F("fAvstracks", "Avstracks",200, 0, 1000, 500, -1000, 1000);
   fOrCvsNtracks    = new TH2F("fCvstracks", "Cvstracks",200, 0, 1000, 500, -1000, 1000);
   fT0vsNtracks     = new TH2F("fT0ACvstrackles", "T0ACvstracks",200, 0, 1000, 500, -1000, 1000); 
+   fEffAC = new TH1F("EffAC","T0Eff ", 200, 0, 1000);
+   fEffA = new TH1F("EffA","T0AEff ", 200, 0, 1000);
+   fEffC = new TH1F("EffC","T0CEff ", 200, 0, 1000);
+   ftracksEffSPD= new TH1F("ftracksEffSPD","SPDeff", 200, 0, 1000);
 
   fTzeroObject     = new TObjArray();
   fTzeroObject->SetOwner(kTRUE);
@@ -101,6 +105,10 @@ void AliT0AnalysisTaskQA::UserCreateOutputObjects()
   fTzeroObject->AddAtAndExpand(fT0vsNtracks, 30);
   fTzeroObject->AddAtAndExpand(fOrAvsNtracks,31);
   fTzeroObject->AddAtAndExpand(fOrCvsNtracks, 32);
+  fTzeroObject->AddAtAndExpand(fEffC, 33);
+  fTzeroObject->AddAtAndExpand(fEffA, 34);
+  fTzeroObject->AddAtAndExpand(fEffAC, 35);
+  fTzeroObject->AddAtAndExpand(ftracksEffSPD, 36);
 
   PostData(1, fTzeroObject);
   // Called once
@@ -139,36 +147,37 @@ void AliT0AnalysisTaskQA::UserExec(Option_t *)
   Double32_t orC = mean[2];
   Int_t ntracks = fESD->GetNumberOfTracks(); 
 
-  if(orA<9999){
+  if(orA<99999){
     fTzeroORA->Fill(orA);
     fOrAvsNtracks->Fill(ntracks, orA);
+    fEffA->Fill(ntracks);
   }
-  if(orC<9999) {
+  if(orC<99999) {
     fTzeroORC->Fill(orC);
     fOrCvsNtracks->Fill(ntracks, orC);
-  }
-  if(orA<9999 && orC<9999) {
+    fEffC->Fill(ntracks);
+   }
+  if(orA<99999 && orC<99999) {
     fResolution->Fill((orA-orC)/2.);
     fTzeroORAplusORC->Fill(mean[0]);
+    fEffAC->Fill(ntracks);
     fT0vsNtracks->Fill(ntracks, mean[0]);
   }
   
   Double32_t t0vertex = fESD->GetT0zVertex();
-  //  cout << "t0 vertex "<<t0vertex<<endl;
   Double32_t esdzvertex;
   const AliESDVertex * esdvertex = fESD->GetPrimaryVertex();
   Int_t nofcontrib=-1;
-  if(esdvertex && t0vertex<999)
+  nofcontrib=esdvertex->GetNContributors();
+  if(nofcontrib>0)    ftracksEffSPD->Fill(ntracks);
+   if(esdvertex && t0vertex<999)
     {
-      nofcontrib=esdvertex->GetNContributors();
-      if(nofcontrib>1)
+      if(nofcontrib>0)
 	{
 	  esdzvertex=esdvertex->GetZv();
-	  //	  cout << "esd vertex "<<esdzvertex<<endl;
 	  fSPDVertexVST0Vertex->Fill(t0vertex,esdzvertex);
 	}
     }
-  // printf("%f   %f  %f\n",orA,orC,time);
   PostData(1, fTzeroObject);
 }      
  //________________________________________________________________________
diff --git a/PWG1/T0/AliT0AnalysisTaskQA.h b/PWG1/T0/AliT0AnalysisTaskQA.h
index 6d6330ed7d0..5a1c97cb48d 100644
--- a/PWG1/T0/AliT0AnalysisTaskQA.h
+++ b/PWG1/T0/AliT0AnalysisTaskQA.h
@@ -39,7 +39,11 @@ class AliT0AnalysisTaskQA : public AliAnalysisTaskSE {
   TH2F        *fSPDVertexVST0Vertex; //! SPD vertex vs T0 vertex   
   TH2F        *fOrAvsNtracks; //! T0A vs Ntracks
   TH2F        *fOrCvsNtracks; //! T0C vs Ntracks
-  TH2F        *fT0vsNtracks; //! T0A vs Ntracks
+  TH2F        *fT0vsNtracks;  //! T0A vs Ntracks
+  TH1F        *fEffAC;        //T0AC eff
+  TH1F        *fEffA;         //T0A eff
+  TH1F        *fEffC;         //T0C eff
+  TH1F        *ftracksEffSPD;  // ntracks in SPD 
   
  
   AliT0AnalysisTaskQA(const AliT0AnalysisTaskQA&); // not implemented