- change from charge/pT to charge
authorjotwinow <jotwinow@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 17 Mar 2011 20:29:32 +0000 (20:29 +0000)
committerjotwinow <jotwinow@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 17 Mar 2011 20:29:32 +0000 (20:29 +0000)
M    TPC/AliPerformanceTPC.cxx
M    TPC/AliPerformanceTPC.h
- set proper ranges for charge selection
M    TPC/AliTPCPerformanceSummary.cxx

PWG1/TPC/AliPerformanceTPC.cxx
PWG1/TPC/AliPerformanceTPC.h
PWG1/TPC/AliTPCPerformanceSummary.cxx

index 907d70a..5794738 100644 (file)
@@ -229,9 +229,12 @@ void AliPerformanceTPC::Init()
   Float_t scaleDCA = 1.0;
   if(fAnalysisMode !=0) scaleDCA = 0.1; 
   // nTPCClust:chi2PerTPCClust:nTPCClustFindRatio:DCAr:DCAz:eta:phi:pt:charge:vertStatus
-   Int_t binsTPCTrackHisto[10]=  { 160,  20,  60,  30, 30,  30,   144,             nPtBins,   nCOverPtBins, 2 };
-   Double_t minTPCTrackHisto[10]={ 0.,   0.,  0., -3*scaleDCA, -3.*scaleDCA, -1.5, 0.,             ptMin,   coverptMin, -0.5 };
-   Double_t maxTPCTrackHisto[10]={ 160., 5., 1.2, 3*scaleDCA,  3.*scaleDCA,  1.5, 2.*TMath::Pi(), ptMax,    coverptMax,  1.5 };
+   //Int_t binsTPCTrackHisto[10]=  { 160,  20,  60,  30, 30,  30,   144,             nPtBins,   nCOverPtBins, 2 };
+   //Double_t minTPCTrackHisto[10]={ 0.,   0.,  0., -3*scaleDCA, -3.*scaleDCA, -1.5, 0.,             ptMin,   coverptMin, -0.5 };
+   //Double_t maxTPCTrackHisto[10]={ 160., 5., 1.2, 3*scaleDCA,  3.*scaleDCA,  1.5, 2.*TMath::Pi(), ptMax,    coverptMax,  1.5 };
+   Int_t binsTPCTrackHisto[10]=  { 160,  20,  60,  30, 30,  30,   144,             nPtBins,   3, 2 };
+   Double_t minTPCTrackHisto[10]={ 0.,   0.,  0., -3*scaleDCA, -3.*scaleDCA, -1.5, 0.,             ptMin,  -1.5, -0.5 };
+   Double_t maxTPCTrackHisto[10]={ 160., 5., 1.2, 3*scaleDCA,  3.*scaleDCA,  1.5, 2.*TMath::Pi(), ptMax,    1.5,  1.5 };
   
   fTPCTrackHisto = new THnSparseF("fTPCTrackHisto","nClust:chi2PerClust:nClust/nFindableClust:DCAr:DCAz:eta:phi:pt:charge/pt:vertStatus",10,binsTPCTrackHisto,minTPCTrackHisto,maxTPCTrackHisto);
   fTPCTrackHisto->SetBinEdges(7,binsPt);
@@ -245,7 +248,8 @@ void AliPerformanceTPC::Init()
   fTPCTrackHisto->GetAxis(5)->SetTitle("#eta");
   fTPCTrackHisto->GetAxis(6)->SetTitle("#phi (rad)");
   fTPCTrackHisto->GetAxis(7)->SetTitle("p_{T} (GeV/c)");
-  fTPCTrackHisto->GetAxis(8)->SetTitle("charge/pt");
+  //fTPCTrackHisto->GetAxis(8)->SetTitle("charge/pt");
+  fTPCTrackHisto->GetAxis(8)->SetTitle("charge");
   fTPCTrackHisto->GetAxis(9)->SetTitle("vertStatus");
   //fTPCTrackHisto->Sumw2();
 
@@ -654,7 +658,7 @@ void AliPerformanceTPC::Analyse()
     // Track histograms 
     // 
     // all with vertex
-    fTPCTrackHisto->GetAxis(8)->SetRangeUser(-10,10);
+    fTPCTrackHisto->GetAxis(8)->SetRangeUser(-1.5,1.5);
     fTPCTrackHisto->GetAxis(9)->SetRangeUser(0.5,1.5);
     selString = "all_recVertex";
     for(Int_t i=0; i <= 9; i++) {
@@ -670,7 +674,7 @@ void AliPerformanceTPC::Analyse()
 
 
     // Track histograms (pos with vertex)
-    fTPCTrackHisto->GetAxis(8)->SetRangeUser(0,10);
+    fTPCTrackHisto->GetAxis(8)->SetRangeUser(0,1.5);
     selString = "pos_recVertex";
     for(Int_t i=0; i <= 9; i++) {
         AddProjection(aFolderObj, "track", fTPCTrackHisto, i, &selString);
@@ -686,7 +690,7 @@ void AliPerformanceTPC::Analyse()
     AddProjection(aFolderObj, "track", fTPCTrackHisto, 5, 6, 7, &selString);
   
     // Track histograms (neg with vertex)
-    fTPCTrackHisto->GetAxis(8)->SetRangeUser(-10,0);
+    fTPCTrackHisto->GetAxis(8)->SetRangeUser(-1.5,0);
     selString = "neg_recVertex";
     for(Int_t i=0; i <= 9; i++) {
         AddProjection(aFolderObj, "track", fTPCTrackHisto, i, &selString);
@@ -702,7 +706,7 @@ void AliPerformanceTPC::Analyse()
     AddProjection(aFolderObj, "track", fTPCTrackHisto, 5, 6, 7, &selString);
 
     //restore cuts
-    fTPCTrackHisto->GetAxis(8)->SetRangeUser(-10,10);
+    fTPCTrackHisto->GetAxis(8)->SetRangeUser(-1.5,1.5);
     fTPCTrackHisto->GetAxis(9)->SetRangeUser(-0.5,1.5);
   
   
index 1845c0b..60490d2 100644 (file)
@@ -106,7 +106,7 @@ private:
   AliPerformanceTPC(const AliPerformanceTPC&); // not implemented
   AliPerformanceTPC& operator=(const AliPerformanceTPC&); // not implemented
 
-  ClassDef(AliPerformanceTPC,9);
+  ClassDef(AliPerformanceTPC,10);
 };
 
 #endif
index fd3a5ea..c26bc80 100644 (file)
@@ -531,9 +531,9 @@ Int_t AliTPCPerformanceSummary::AnalyzeDCARPhiPos(const AliPerformanceTPC* pTPC,
     if (his3D && !fgForceTHnSparse) { 
         his2D = dynamic_cast<TH2*>(his3D->Project3D("xy")); 
     } else {    
-        pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(0,10);        
+        pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(0,1.5);        
         his2D = pTPC->GetTPCTrackHisto()->Projection(3,5);
-        pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-10,10);
+        pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-1.5,1.5);
     }            
     
     his2D->FitSlicesY(0,0,-1,10,"QNR",&arrayFit);
@@ -619,9 +619,9 @@ Int_t AliTPCPerformanceSummary::AnalyzeDCARPhiNeg(const AliPerformanceTPC* pTPC,
     if (his3D && !fgForceTHnSparse) {
         his2D = dynamic_cast<TH2*>(his3D->Project3D("xy")); 
     } else {    
-        pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-10,0);        
+        pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-1.5,0);        
         his2D = pTPC->GetTPCTrackHisto()->Projection(3,5);
-        pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-10,10);
+        pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-1.5,1.5);
     }            
     his2D->FitSlicesY(0,0,-1,10,"QNR",&arrayFit);
     delete his2D;    
@@ -937,8 +937,6 @@ Int_t AliTPCPerformanceSummary::AnalyzeDriftPos(const AliPerformanceTPC* pTPC, T
     TH2* his2D=0;
     TH3* his3D=0;
     
-    //pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(0,10);    
-    
    if (pTPC->GetHistos()->FindObject("h_tpc_track_pos_recvertex_4_5_7")) {    
         his3D = dynamic_cast<TH3*>(pTPC->GetHistos()->FindObject("h_tpc_track_pos_recvertex_4_5_7"));
         his3D->GetYaxis()->SetRangeUser(-1,1);
@@ -954,7 +952,6 @@ Int_t AliTPCPerformanceSummary::AnalyzeDriftPos(const AliPerformanceTPC* pTPC, T
     TObjArray arrayFit;
     his2D->FitSlicesY(0,0,-1,10,"QNR",&arrayFit);
     delete his2D;
-    //pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-10,10);
     
     his1D = (TH1*) arrayFit.At(1);
     his1D->Fit(fpol1,"QNRROB=0.8","QNR",-0.8,-0.1);
@@ -1027,7 +1024,6 @@ Int_t AliTPCPerformanceSummary::AnalyzeDriftNeg(const AliPerformanceTPC* pTPC, T
     TH2* his2D=0;
     TH3* his3D=0;
     
-    //pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-10,0);    
     
    if (pTPC->GetHistos()->FindObject("h_tpc_track_neg_recvertex_4_5_7")) {    
         his3D = dynamic_cast<TH3*>(pTPC->GetHistos()->FindObject("h_tpc_track_neg_recvertex_4_5_7"));
@@ -1045,7 +1041,6 @@ Int_t AliTPCPerformanceSummary::AnalyzeDriftNeg(const AliPerformanceTPC* pTPC, T
     his2D->FitSlicesY(0,0,-1,10,"QNR",&arrayFit);
     delete his2D;
     
-    // pTPC->GetTPCTrackHisto()->GetAxis(8)->SetRangeUser(-10,10);
     his1D = (TH1*) arrayFit.At(1);
     his1D->Fit(fpol1,"QNRROB=0.8","QNR",-0.8,-0.1);
     offsetdZCNeg=fpol1->GetParameter(0);