Adding chi2 of the fit.

(Marian)

diff --git a/PWG1/TPC/macros/MakeReportTPC.C b/PWG1/TPC/macros/MakeReportTPC.C
index f24da2a278c1adc1f99a22108fc3ecc8a25319e5..492cd4ba53892cc9c6184a3961cdc1c5b9d2fe86 100644 (file)
--- a/PWG1/TPC/macros/MakeReportTPC.C
+++ b/PWG1/TPC/macros/MakeReportTPC.C
@@ -8,11 +8,15 @@
    gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+")
    AliXRDPROOFtoolkit tool;
    TChain * chain = tool.MakeChain("summaryTPCQA.txt","tpcQA",0,200000);
+
+   TTree *tree = chain->CopyTree("1");
 */
 
 
 #if !defined(__CINT__) || defined(__MAKECINT__)
 #include "TSystem.h"
+#include "TMath.h"
+#include "TVectorD.h"
 #include "TFile.h"
 #include "TGrid.h"
 #include "TF1.h"
@@ -22,15 +26,18 @@
 #include "THnSparse.h"
 #include "TTreeStream.h"
 #include "AliPerformanceTPC.h"
+#include "AliPerformanceDEdx.h"
 #endif
 
 TObject *pTPCObject=0;
+TObject *pTPCObjectGain=0;
 TTreeSRedirector  *pcstream=0;
 void Init();
 void ReadObjects(const char * path=0);
 void MakeReportTPC(Int_t run, const char *path=0 );
 void AnalyzeNCL();
 void AnalyzeDrift();
+void AnalyzeGain();
 
 void Init(){
   //
@@ -45,9 +52,9 @@ void Init(){
   gSystem->Load("libPWG0selectors.so");
   gSystem->Load("libPWG1.so");
   gSystem->Load("libTPCcalib"); 
-  gSystem->Setenv("alien_CLOSE_SE","ALICE::GSI::SE");
-  TGrid::Connect("alien://",0,0,"t"); 
-  gSystem->Setenv("alien_CLOSE_SE","ALICE::GSI::SE");    
+//   gSystem->Setenv("alien_CLOSE_SE","ALICE::GSI::SE");
+//   TGrid::Connect("alien://",0,0,"t"); 
+//   gSystem->Setenv("alien_CLOSE_SE","ALICE::GSI::SE");    
 }
 
 
@@ -68,6 +75,7 @@ void ReadObjects(const char * path){
     return;
   }
   pTPCObject= (AliPerformanceTPC*)list->FindObject("AliPerformanceTPC");
+  pTPCObjectGain = (AliPerformanceDEdx*)list->FindObject("AliPerformanceDEdxTPCInner");
 }
 
 
@@ -83,6 +91,7 @@ void MakeReportTPC(Int_t run, const char *path ){
   (*pcstream)<<"tpcQA"<<"run="<<run;
   AnalyzeNCL();
   AnalyzeDrift();
+  AnalyzeGain();
   (*pcstream)<<"tpcQA"<<"\n";
   delete pcstream;
 }
@@ -97,10 +106,14 @@ void AnalyzeNCL(){
   static Double_t rmsTPCnclF=0;
   static Double_t slopeATPCnclF=0;
   static Double_t slopeCTPCnclF=0;
+  static Double_t slopeATPCnclFErr=0;
+  static Double_t slopeCTPCnclFErr=0;
   static Double_t meanTPCncl=0;
   static Double_t rmsTPCncl=0;
   static Double_t slopeATPCncl=0;
   static Double_t slopeCTPCncl=0;
+  static Double_t slopeATPCnclErr=0;
+  static Double_t slopeCTPCnclErr=0;
   AliPerformanceTPC * pTPC=  (AliPerformanceTPC *)pTPCObject;
 
   TH1* his1D=0;
@@ -119,8 +132,10 @@ void AnalyzeNCL(){
   hprof = pTPC->GetTPCTrackHisto()->Projection(0,5)->ProfileX();
   hprof->Fit(fpol1,"QNR","QNR",0,0.8);
   slopeATPCncl= fpol1->GetParameter(1);
+  slopeATPCnclErr= fpol1->GetParError(1);
   hprof->Fit(fpol1,"QNR","QNR",-0.8,0.0);
   slopeCTPCncl= fpol1->GetParameter(1);
+  slopeCTPCnclErr= fpol1->GetParameter(1);
   delete hprof;
   //
   // findable clusters
@@ -129,12 +144,14 @@ void AnalyzeNCL(){
   his1D = pTPC->GetTPCTrackHisto()->Projection(2);
   meanTPCnclF= his1D->GetMean();
   rmsTPCnclF= his1D->GetRMS();
-   delete his1D;
+  delete his1D;
   his1D = pTPC->GetTPCTrackHisto()->Projection(2,5)->ProfileX();
   his1D->Fit(fpol1,"QNR","QNR",0,0.8);
   slopeATPCnclF= fpol1->GetParameter(1);
+  slopeATPCnclFErr= fpol1->GetParError(1);
   his1D->Fit(fpol1,"QNR","QNR",-0.8,0.0);
   slopeCTPCnclF= fpol1->GetParameter(1);
+  slopeCTPCnclFErr= fpol1->GetParameter(1);
   delete his1D;
   printf("Cluster QA report\n");
   printf("meanTPCnclF=\t%f\n",meanTPCnclF);
@@ -149,14 +166,18 @@ void AnalyzeNCL(){
   // dump results to the tree
   //
   (*pcstream)<<"tpcQA"<<
-    "meanTPCnclF="<<meanTPCnclF <<
+    "meanTPCnclF="<<meanTPCnclF <<   
     "rmsTPCnclF="<<rmsTPCnclF <<
     "slopeATPCnclF="<< slopeATPCnclF<<
     "slopeCTPCnclF="<< slopeCTPCnclF<<
+    "slopeATPCnclFErr="<< slopeATPCnclFErr<<
+    "slopeCTPCnclFErr="<< slopeCTPCnclFErr<<
     "meanTPCncl="<<meanTPCncl <<
     "rmsTPCncl="<< rmsTPCncl<<
     "slopeATPCncl="<< slopeATPCncl<<
-    "slopeCTPCncl="<< slopeCTPCncl;
+    "slopeCTPCncl="<< slopeCTPCncl<<
+    "slopeATPCnclErr="<< slopeATPCnclErr<<
+    "slopeCTPCnclErr="<< slopeCTPCnclErr;
 
     
 }
@@ -170,6 +191,14 @@ void AnalyzeDrift(){
   static Double_t slopedZA=0;
   static Double_t offsetdZC=0;
   static Double_t slopedZC=0;
+  static Double_t offsetdZAErr=0;
+  static Double_t slopedZAErr=0;
+  static Double_t offsetdZCErr=0;
+  static Double_t slopedZCErr=0;
+  static Double_t offsetdZAchi2=0;
+  static Double_t slopedZAchi2=0;
+  static Double_t offsetdZchi2=0;
+  static Double_t slopedZCchi2=0;
   AliPerformanceTPC * pTPC=  (AliPerformanceTPC *)pTPCObject;
 
   TH1* his1D=0;
@@ -180,12 +209,19 @@ void AnalyzeDrift(){
   his2D->FitSlicesY(0,0,-1,10,"QNR",&arrayFit);
   delete his2D;
   his1D = (TH1*) arrayFit.At(1);
-  his1D->Fit(fpol1,"QNR","QNR",-0.8,-0.1);
+  his1D->Fit(fpol1,"QNRROB=0.8","QNR",-0.8,-0.1);
   offsetdZC=fpol1->GetParameter(0);
   slopedZC=fpol1->GetParameter(1);
-  his1D->Fit(fpol1,"QNR","QNR",0.1,0.8);
+  offsetdZCchi2=fpol1->GetChisquare();
+  slopedZCchi2=fpol1->GetChisquare();
+  //
+  his1D->Fit(fpol1,"QNRROB=0.8","QNR",0.1,0.8);
   offsetdZA=fpol1->GetParameter(0);
   slopedZA=fpol1->GetParameter(1);
+  offsetdZAErr=fpol1->GetParError(0);
+  slopedZAErr=fpol1->GetParError(1);
+  offsetdZAchi2=fpol1->GetChisquare();
+  slopedZAchi2=fpol1->GetChisquare();
   //
   printf("Drift velocity QA report\n");
   printf("offsetdZA\t%f\n",offsetdZA);
@@ -199,6 +235,116 @@ void AnalyzeDrift(){
     "offsetdZA="<< offsetdZA<<
     "slopedZA="<< slopedZA<<
     "offsetdZC="<< offsetdZC<<
-    "slopedZC="<<slopedZC;
+    "slopedZC="<<slopedZC<<
+    //
+    "offsetdZAErr="<< offsetdZAErr<<
+    "slopedZAErr="<< slopedZAErr<<
+    "offsetdZCErr="<< offsetdZCErr<<
+    "slopedZCErr="<<slopedZCErr<<
+    //
+    "offsetdZAchi2="<< offsetdZAchi2<<
+    "slopedZAchi2="<< slopedZAchi2<<
+    "offsetdZCchi2="<< offsetdZCchi2<<
+    "slopedZCchi2="<<slopedZCchi2;
   
 }
+
+
+void AnalyzeGain(){
+  //
+  // gain statistic
+  //
+  AliPerformanceDEdx * pTPCgain =  (AliPerformanceDEdx *)pTPCObjectGain;
+  static TVectorD MIPvsSector(36);
+  static TVectorD sector(36);
+  static Float_t MIPmean = 0;
+  static Float_t MIPresolution = 0;
+  static Float_t attachSlopeC = 0;
+  static Float_t attachSlopeA = 0;
+
+  MIPvsSector.Zero();
+  //
+  // select MIP particles
+  //
+  pTPCgain->GetDeDxHisto()->GetAxis(7)->SetRangeUser(0.4,0.55);
+  pTPCgain->GetDeDxHisto()->GetAxis(0)->SetRangeUser(35,60);
+  pTPCgain->GetDeDxHisto()->GetAxis(6)->SetRangeUser(80,160);
+  pTPCgain->GetDeDxHisto()->GetAxis(5)->SetRangeUser(-1,1);
+  //
+  // MIP position and resolution
+  //
+  TF1 gausFit("gausFit","gaus");
+  TH1 * hisProj1D = pTPCgain->GetDeDxHisto()->Projection(0);
+  hisProj1D->Fit(&gausFit,"QN","QN");
+  delete hisProj1D;
+  MIPmean = gausFit.GetParameter(1);
+  MIPresolution = 0;
+  if (MIPmean!=0) MIPresolution = gausFit.GetParameter(2)/MIPmean;
+  //
+  // MIP position vs. dip angle (attachment)
+  //
+  pTPCgain->GetDeDxHisto()->GetAxis(5)->SetRangeUser(-3,0);
+  TH2* his2D = pTPCgain->GetDeDxHisto()->Projection(0,5);
+  TF1 * fpol = new TF1("fpol","pol1");
+  TObjArray arrayFit;
+  his2D->FitSlicesY(0,0,-1,10,"QN",&arrayFit);
+  delete his2D;
+  TH1 * his1D = (TH1*) arrayFit.At(1);
+  his1D->Fit(fpol,"QNROB=0.8","QNR",-1,0);
+  attachSlopeC = fpol->GetParameter(1);
+  //
+  pTPCgain->GetDeDxHisto()->GetAxis(5)->SetRangeUser(0,3);
+  TH2* his2DA = pTPCgain->GetDeDxHisto()->Projection(0,5);
+  TF1 * fpolA = new TF1("fpolA","pol1");
+  TObjArray arrayFitA;
+  his2DA->FitSlicesY(0,0,-1,10,"QN",&arrayFit);
+  delete his2DA;
+  TH1 * his1DA = (TH1*) arrayFit.At(1);
+  his1DA->Fit(fpolA,"QNROB=0.8","QN",0,1);
+  attachSlopeA = fpolA->GetParameter(1);
+  //
+  // MIP position vs. sector
+  //
+  pTPCgain->GetDeDxHisto()->GetAxis(5)->SetRangeUser(-3,0); // C side
+  for(Int_t i = 0; i < 18; i++) { // loop over sectors; correct mapping to be checked!
+    TH1* his1D=0;
+    Float_t phiLow = -TMath::Pi() + i*(20./360.)*(2*TMath::Pi());
+    Float_t phiUp  = -TMath::Pi() + (i+1)*(20./360.)*(2*TMath::Pi());
+    pTPCgain->GetDeDxHisto()->GetAxis(1)->SetRangeUser(phiLow,phiUp);
+    his1D = pTPCgain->GetDeDxHisto()->Projection(0);
+    TF1 gausFunc("gausFunc","gaus");
+    his1D->Fit(&gausFunc, "QN");
+    MIPvsSector(i) = gausFunc.GetParameter(1);
+    sector(i)=i;
+    delete his1D;
+  }
+  //
+  pTPCgain->GetDeDxHisto()->GetAxis(5)->SetRangeUser(0,3); // A side
+  for(Int_t i = 0; i < 18; i++) { // loop over sectors; correct mapping to be checked!
+    TH1* his1D=0;
+    Float_t phiLow = -TMath::Pi() + i*(20./360.)*(2*TMath::Pi());
+    Float_t phiUp  = -TMath::Pi() + (i+1)*(20./360.)*(2*TMath::Pi());
+    pTPCgain->GetDeDxHisto()->GetAxis(1)->SetRangeUser(phiLow,phiUp);
+    his1D = pTPCgain->GetDeDxHisto()->Projection(0);
+    TF1 gausFunc("gausFunc","gaus");
+    his1D->Fit(&gausFunc, "QN");
+    MIPvsSector(i+18) = gausFunc.GetParameter(1);
+    sector(i+18)=i+18;
+    delete his1D;
+  }
+  //
+  printf("Gain QA report\n");
+  printf("MIP mean\t%f\n",MIPmean);
+  printf("MIP resolution\t%f\n",MIPresolution);
+  printf("MIPslopeA\t%f\n",attachSlopeA);
+  printf("MIPslopeC\t%f\n",attachSlopeC);
+  //
+  (*pcstream)<<"tpcQA"<<
+    "MIPattachSlopeC="<<attachSlopeC<<
+    "MIPattachSlopeA="<<attachSlopeA<<
+    "MIPresolution="<<MIPresolution<<
+    "MIPvsSector.="<<&MIPvsSector<<
+    "sector.="<<&sector<<
+    "MIPmean="<<MIPmean;
+
+}