Added plots fpr PID qa from PIDqa task + addedd mean and sigma for Nsigma_TOF + cleanup

[u/mrichter/AliRoot.git] / TPC / Calib / AliTPCcalibLaser.cxx

diff --git a/TPC/Calib/AliTPCcalibLaser.cxx b/TPC/Calib/AliTPCcalibLaser.cxx
index 7e1aa7f5fc387c393df92ab716298f35c1852ffe..cac5c4a668fec2f877ae658f9f0bddca4b9374af 100644 (file)
--- a/TPC/Calib/AliTPCcalibLaser.cxx
+++ b/TPC/Calib/AliTPCcalibLaser.cxx
@@ -90,15 +90,9 @@
 #include "TLinearFitter.h"
 #include "AliTPCcalibLaser.h"
 #include "AliExternalTrackParam.h"
-//#include "AliESDEvent.h"
-//#include "AliESDfriend.h"
-//#include "AliESDtrack.h"
-
-#include "AliVEvent.h"
-#include "AliVTrack.h"
-#include "AliVfriendEvent.h"
-#include "AliVfriendTrack.h"
-
+#include "AliESDEvent.h"
+#include "AliESDfriend.h"
+#include "AliESDtrack.h"
 #include "AliTPCTracklet.h"
 #include "TH1D.h"
 #include "TH1F.h"
@@ -142,11 +136,11 @@ ClassImp(AliTPCcalibLaser)
 
 AliTPCcalibLaser::AliTPCcalibLaser():
   AliTPCcalibBase(),
-  fEvent(0),
-  fEventFriend(0),
+  fESD(0),
+  fESDfriend(0),
   fNtracks(0),
   fTracksMirror(336),
-  fTracks(336),
+  fTracksEsd(336),
   fTracksEsdParam(336),
   fTracksTPC(336),
   fFullCalib(kTRUE),
@@ -251,11 +245,11 @@ AliTPCcalibLaser::AliTPCcalibLaser():
 
 AliTPCcalibLaser::AliTPCcalibLaser(const Text_t *name, const Text_t *title, Bool_t full):
   AliTPCcalibBase(),
-  fEvent(0),
-  fEventFriend(0),
+  fESD(0),
+  fESDfriend(0),
   fNtracks(0),
   fTracksMirror(336),
-  fTracks(336),
+  fTracksEsd(336),
   fTracksEsdParam(336),
   fTracksTPC(336),
   fFullCalib(full),
@@ -367,11 +361,11 @@ AliTPCcalibLaser::AliTPCcalibLaser(const Text_t *name, const Text_t *title, Bool
 
 AliTPCcalibLaser::AliTPCcalibLaser(const AliTPCcalibLaser& calibLaser):
   AliTPCcalibBase(calibLaser), 
-  fEvent(0),
-  fEventFriend(0),
+  fESD(0),
+  fESDfriend(0),
   fNtracks(0),
   fTracksMirror(336),
-  fTracks(336),
+  fTracksEsd(336),
   fTracksEsdParam(336),
   fTracksTPC(336),
   fFullCalib(calibLaser.fFullCalib),
@@ -581,7 +575,7 @@ AliTPCcalibLaser::~AliTPCcalibLaser() {
 
 
 
-void AliTPCcalibLaser::Process(AliVEvent *event) {
+void AliTPCcalibLaser::Process(AliESDEvent * event) {
   //
   //
   // Loop over tracks and call  Process function
@@ -589,25 +583,24 @@ void AliTPCcalibLaser::Process(AliVEvent *event) {
   const Int_t  kMinTracks=20;
   const Int_t  kMinClusters=40;
 
-  fEvent = event;
-  if (!fEvent) {
-      //Printf("****ERROR AliTPCcalibLaser::Process(): Event not available!!********");
+  fESD = event;
+  if (!fESD) {
     return;
   }
-  fEventFriend=fEvent->FindFriend();
-  if (!fEventFriend) {
+  fESDfriend=static_cast<AliESDfriend*>(fESD->FindListObject("AliESDfriend"));
+  if (!fESDfriend) {
     return;
   }
-  if (fEventFriend->TestSkipBit()) return;
-  if (fEvent->GetNumberOfTracks()<kMinTracks) return; //not enough tracks
+  if (fESDfriend->TestSkipBit()) return;
+  if (fESD->GetNumberOfTracks()<kMinTracks) return; //not enough tracks
   AliDebug(4,Form("Event number in current file: %d",event->GetEventNumberInFile()));
   //
   // find CE background if present
   //
   if (AliTPCLaserTrack::GetTracks()==0) AliTPCLaserTrack::LoadTracks();
   TH1D hisCE("hhisCE","hhisCE",100,-100,100);  
-  for (Int_t i=0;i<fEvent->GetNumberOfTracks();++i) {
-    AliVTrack *track=fEvent->GetVTrack(i);
+  for (Int_t i=0;i<fESD->GetNumberOfTracks();++i) {
+    AliESDtrack *track=fESD->GetTrack(i);
     if (!track) continue;
     hisCE.Fill(track->GetZ());
     hisCE.Fill(track->GetZ()+2);
@@ -618,7 +611,7 @@ void AliTPCcalibLaser::Process(AliVEvent *event) {
 
 
   fTracksTPC.Clear();
-  fTracks.Clear();
+  fTracksEsd.Clear();
   fTracksEsdParam.Delete();
   for (Int_t id=0; id<336;id++) {
     fCounter[id]=0;
@@ -626,12 +619,12 @@ void AliTPCcalibLaser::Process(AliVEvent *event) {
     fClusterSatur[id]=0;
   }
   //
-  Int_t n=fEvent->GetNumberOfTracks();
+  Int_t n=fESD->GetNumberOfTracks();
   Int_t counter=0;
   for (Int_t i=0;i<n;++i) {
-    const AliVfriendTrack *friendTrack=fEventFriend->GetTrack(i);
+    AliESDfriendTrack *friendTrack=fESDfriend->GetTrack(i);
     if (!friendTrack) continue;
-    AliVTrack *track= fEvent->GetVTrack(i);
+    AliESDtrack *track=fESD->GetTrack(i);
     if (!track) continue;
     Double_t binC = hisCE.GetBinContent(hisCE.FindBin(track->GetZ()));
     if (binC>336) continue; //remove CE background
@@ -696,7 +689,7 @@ void AliTPCcalibLaser::MakeDistHisto(Int_t id){
 
     AliExternalTrackParam *param=(AliExternalTrackParam*)fTracksEsdParam.At(id);
     AliTPCLaserTrack *ltrp = ( AliTPCLaserTrack*)fTracksMirror.At(id);
-    AliVTrack   *track    = (AliVTrack*)fTracks.At(id);
+    AliESDtrack   *track    = (AliESDtrack*)fTracksEsd.At(id);
     if (!param) return;
     if (!ltrp) return;
     if (!track) return;
@@ -803,7 +796,7 @@ void AliTPCcalibLaser::FitDriftV(){
       if (!AcceptLaser(id)) continue;
       if ( fClusterSatur[id]>kSaturCut)  continue;
       if ( fClusterCounter[id]<kMinClusters)  continue;
-      AliVTrack   *track    = (AliVTrack*)fTracks.At(id);
+      AliESDtrack   *track    = (AliESDtrack*)fTracksEsd.At(id);
       if (track->GetTPCsignal()<kMinSignal) continue;
       AliExternalTrackParam *param=(AliExternalTrackParam*)fTracksEsdParam.At(id);
       AliTPCLaserTrack *ltrp = ( AliTPCLaserTrack*)fTracksMirror.At(id);
@@ -1106,7 +1099,7 @@ Bool_t  AliTPCcalibLaser::FitDriftV(Float_t minFraction){
       if (!AcceptLaser(id)) continue;
       if ( fClusterSatur[id]>kSaturCut)  continue;
       if ( fClusterCounter[id]<kMinClusters)  continue;
-      AliVTrack   *track    = (AliVTrack*)fTracks.At(id);
+      AliESDtrack   *track    = (AliESDtrack*)fTracksEsd.At(id);
       if (track->GetTPCsignal()<kMinSignal) continue;
       AliExternalTrackParam *param=(AliExternalTrackParam*)fTracksEsdParam.At(id);
       AliTPCLaserTrack *ltrp = ( AliTPCLaserTrack*)fTracksMirror.At(id);
@@ -1363,7 +1356,7 @@ Bool_t  AliTPCcalibLaser::AcceptLaser(Int_t id){
   return kTRUE;
 }
 
-Int_t  AliTPCcalibLaser::FindMirror(AliVTrack *track, AliTPCseed *seed){
+Int_t  AliTPCcalibLaser::FindMirror(AliESDtrack *track, AliTPCseed *seed){
   //
   // Find corresponding mirror
   // add the corresponding tracks
@@ -1443,7 +1436,7 @@ Int_t  AliTPCcalibLaser::FindMirror(AliVTrack *track, AliTPCseed *seed){
   if (accept){
     fClusterCounter[id]=countercl;
     fTracksEsdParam.AddAt(param.Clone(),id);
-    fTracks.AddAt(track,id);
+    fTracksEsd.AddAt(track,id);
     fTracksTPC.AddAt(seed,id);
   }
   return id;
@@ -1455,7 +1448,7 @@ void AliTPCcalibLaser::DumpLaser(Int_t id) {
   //
   //  Dump Laser info to the tree
   //
-  AliVTrack   *track    = (AliVTrack*)fTracks.At(id);
+  AliESDtrack   *track    = (AliESDtrack*)fTracksEsd.At(id);
   AliExternalTrackParam *param=(AliExternalTrackParam*)fTracksEsdParam.At(id);
   AliTPCLaserTrack *ltrp = ( AliTPCLaserTrack*)fTracksMirror.At(id);
   //
@@ -1476,7 +1469,7 @@ void AliTPCcalibLaser::DumpLaser(Int_t id) {
   
   if (fStreamLevel>0){
     TTreeSRedirector *cstream = GetDebugStreamer();
-    Int_t time = fEvent->GetTimeStamp();
+    Int_t time = fESD->GetTimeStamp();
     Bool_t accept = AcceptLaser(id);
     if (cstream){
       (*cstream)<<"Track"<<