- // reconstructed vertex status > 0
- fTPCEventHisto->GetAxis(6)->SetRange(2,2);
- for(Int_t i=0; i<6; i++)
- {
- h = (TH1F*)fTPCEventHisto->Projection(i);
- sprintf(name,"h_tpc_event_recVertex%d",i);
- h->SetName(name);
- h->GetXaxis()->SetTitle(fTPCEventHisto->GetAxis(i)->GetTitle());
- h->GetYaxis()->SetTitle("events");
- sprintf(title,"%s rec. vertex",fTPCEventHisto->GetAxis(i)->GetTitle());
- h->SetTitle(title);
-
- aFolderObj->Add(h);
- }
-
- //
- // Track histograms
- //
- for(Int_t i=0; i<9; i++)
- {
- h = (TH1F*)fTPCTrackHisto->Projection(i);
- sprintf(name,"h_tpc_track_%d",i);
- h->SetName(name);
- h->GetXaxis()->SetTitle(fTPCTrackHisto->GetAxis(i)->GetTitle());
- h->GetYaxis()->SetTitle("tracks");
- sprintf(title,"%s",fTPCTrackHisto->GetAxis(i)->GetTitle());
- h->SetTitle(title);
-
- if(i==7) h->Scale(1,"width");
- aFolderObj->Add(h);
- }
+//_____________________________________________________________________________
+void AliPerformanceTPC::Analyse()
+{
+ //
+ // Analyse comparison information and store output histograms
+ // in the folder "folderTPC"
+ //
+ TH1::AddDirectory(kFALSE);
+ TH1::SetDefaultSumw2(kFALSE);
+ TObjArray *aFolderObj = new TObjArray;
+ //aFolderObj->SetOwner(); // objects are owned by fanalysisFolder
+ TString selString;
+
+ //
+ // Cluster histograms
+ //
+ AddProjection(aFolderObj, fTPCClustHisto, 0, 1, 2);
+
+ selString = "all";
+ for(Int_t i=0; i <= 2; i++) {
+ AddProjection(aFolderObj, fTPCClustHisto, i, &selString);
+ }
+
+ //fTPCClustHisto->GetAxis(2)->SetRange(1,1); // A-side
+ //selString = "A_side";
+ //AddProjection(aFolderObj, fTPCClustHisto, 0, 1, &selString);
+
+ //fTPCClustHisto->GetAxis(2)->SetRange(2,2); // C-side
+ //selString = "C_side";
+ //AddProjection(aFolderObj, fTPCClustHisto, 0, 1, &selString);
+
+ //reset range
+ fTPCClustHisto->GetAxis(2)->SetRange(1,2);
+
+ //
+ // event histograms
+ //
+ for(Int_t i=0; i<=6; i++) {
+ AddProjection(aFolderObj, fTPCEventHisto, i);
+ }
+ AddProjection(aFolderObj, fTPCEventHisto, 4, 5);
+ AddProjection(aFolderObj, fTPCEventHisto, 0, 1);
+ AddProjection(aFolderObj, fTPCEventHisto, 0, 3);
+ AddProjection(aFolderObj, fTPCEventHisto, 1, 3);
+ AddProjection(aFolderObj, fTPCEventHisto, 2, 3);
+
+ // reconstructed vertex status > 0
+ fTPCEventHisto->GetAxis(6)->SetRange(2,2);
+ selString = "recVertex";
+ for(Int_t i=0; i<=5; i++) {
+ AddProjection(aFolderObj, fTPCEventHisto, i, &selString);
+ }
+ AddProjection(aFolderObj, fTPCEventHisto, 4, 5, &selString);
+ AddProjection(aFolderObj, fTPCEventHisto, 0, 1, &selString);
+ AddProjection(aFolderObj, fTPCEventHisto, 0, 3, &selString);
+ AddProjection(aFolderObj, fTPCEventHisto, 1, 3, &selString);
+ AddProjection(aFolderObj, fTPCEventHisto, 2, 3, &selString);
+
+ // reset cuts
+ fTPCEventHisto->GetAxis(6)->SetRange(1,2);
+
+ //
+ // Track histograms
+ //
+ // all with vertex
+ 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++) {
+ AddProjection(aFolderObj, fTPCTrackHisto, i, &selString);
+ }
+ for(Int_t i=0; i <= 4; i++) {
+ AddProjection(aFolderObj, fTPCTrackHisto, i, 5, 7, &selString);
+ }