[u/mrichter/AliRoot.git] / PWGPP / TPC / macros / compareTracks.C

/*

  //
  // Macro to compare  tracks used 2 different reconstruction algorithms (Offline and HLT)
  // (In general the approach  can be used for any pair of the esds - track containers)
  //   
  // Main functions - Dump  content of HLT and OFFLINE esd files  - full esd tracks 
  //  1. Dump to the tree the pair of close tracks ( n sigma distance at the TPC entrance)
  //  2. Dump to the tree the track and corrsponding counter of tracks from other container 
  //    Compare two events - track by track comparison
  //       0. Filter "Good" tpc tracks - see function IsSelected(track)
  //       1. Dump OFFline -> HLT tracks into trees  
  //           a.) tree "offhlt"  - counters+ Ofline track + corresponding  Hlt track
  //           b.) tree "offhlt0" - offline track + hlt track counter 
  //       2. HLT and offline
  //     
  //  3. To save the CPU and disk space the input tracks to compare are pt downscaled
  //      see function IsDownscaled();
  //
  //  Visualization part is currently work in progress
  //  
  // Responsible:
  // dumping part:         marian.ivanov@cern.ch
  // visualization part:   
  // 
  Example usage:
  aliroot -b -q $ALICE_ROOT/TPC/macros/compareTracks.C+
  //
  //
  .x $HOME/rootlogon.C
  gSystem->Load("libANALYSIS");
  gSystem->Load("libANALYSISalice");
  gSystem->Load("libPWG0base");
  gSystem->Load("libPWG0dep");
  gSystem->Load("libPWG0selectors");
  .L $ALICE_ROOT/PWGPP/TPC/macros/compareTracks.C+
  //
  compareTracks("compare.list");
  //
  TFile f("dump.root");
*/
#include <exception>
#include <fstream>
#include "TFile.h"
#include "TTree.h"
#include "TChain.h"
#include "TCut.h"
#include "TStopwatch.h"
#include "TSystem.h"
#include "AliESDEvent.h"
#include "TTreeStream.h"
#include "TRandom.h"
#include "AliXRDPROOFtoolkit.h"
#include "AliSysInfo.h"

//
//
//
void CompareFile(TTreeSRedirector *pcstream,  TFile *fOff, TFile *fHLT);
void CompareEvents(AliESDEvent *evoff, AliESDEvent *evhlt, TTreeSRedirector *pcstream);
Bool_t AreTracksCloseFast(AliESDtrack *track1, AliESDtrack *track2, AliESDEvent * event, AliExternalTrackParam &inner2);
Bool_t IsSelected(AliESDtrack *track);
Bool_t IsDownscaled(AliESDtrack *track);
//
Double_t gptdownscale=100;  //pt downscale parameter
//
TTree * chainOFFHLT=0;
TTree * chainOFFHLT0=0;
TTree * chainHLTOFF=0;
TTree * chainHLTOFF0=0;

void compareTracks(const char * flistOFF="compareOFF.list", const char * flistHLT="compareHLT.list", Double_t downscale=10000){
    //
    // Compare HLT and OFFLINE esd files 
    // Input: 
    //   flist - the ascii files with the filename - ESD form the offline
    //         - hlt path constructed  from the offline ESD file repalcing the OFFLINE with HLT
    //
    gptdownscale=downscale;
    TTreeSRedirector *pcstream = new TTreeSRedirector("dump.root");
    ifstream inOFF;
    inOFF.open(flistOFF);
    ifstream inHLT;
    inHLT.open(flistHLT);


    // Read the input list of files and add them to the chain
    TString currentFileOff;
    TString currentFileHlt;
    {  while( (inOFF.good()) && (inHLT.good())) {
		        inOFF >> currentFileOff;
			inHLT >> currentFileHlt;
	                //currentFileHlt=currentFileOff;
	                //currentFileHlt.ReplaceAll("/OFFLINE/","/HLT/");
	                printf("%s\n%s\n\n",currentFileOff.Data(), currentFileHlt.Data());
	                TFile *fOff = TFile::Open(currentFileOff.Data(),"READ");
		        TFile *fHlt = TFile::Open(currentFileHlt.Data(),"READ");
	                if ( fOff && fHlt){
		           CompareFile(pcstream, fOff, fHlt);
		        }
		      delete fOff;
		      delete fHlt;
    }}
    delete pcstream;
		        /* 
			 *      TTreeSRedirector *pcstream = new TTreeSRedirector("dump.root");
			 *           TFile *fOff = TFile::Open("/lustre/alice/mknichel/cpass1/2011-12-13_0154/OFFLINE/167693/11000167693001.11/AliESDs.root","READ");
			 *                TFile *fHLT = TFile::Open("/lustre/alice/mknichel/cpass1/2011-12-13_0154/HLT/167693/11000167693001.11/AliESDs.root","READ");
			 *                     CompareFile(pcstream, fOff, fHLT);
			 *                          delete pcstream;
			 *                            */
}

void CompareFile(TTreeSRedirector *pcstream,  TFile *fOff, TFile *fHLT){
  //
  // Compare 2 esd files - fOff and fHLT
  // Pairs of tracks dumpt in the trees - stored in the filee connected to the pcstream
  // 
  TTree *tOff = (TTree*)fOff->Get("esdTree");
  TTree *tHLT = (TTree*)fHLT->Get("esdTree");
  //    
  Int_t nevents=tOff->GetEntries();
  //
  AliESDEvent *evoff = new AliESDEvent();
  AliESDEvent *evhlt = new AliESDEvent();  
  evoff->ReadFromTree(tOff);
  evhlt->ReadFromTree(tHLT);
  //
  for(Long64_t iev = 0; iev < nevents; iev++){
    tOff->GetEvent(iev);
    tHLT->GetEvent(iev);
    CompareEvents(evoff,evhlt,pcstream);    
    AliSysInfo::AddStamp(fOff->GetName(),iev);
  }
  delete evoff;
  delete evhlt;
}


void CompareEvents(AliESDEvent *evoff, AliESDEvent *evhlt, TTreeSRedirector *pcstream){
  //
  // Main function to dump events:
  //
  // compare two events - track by track comparison
  // 0. Filter "Good" tpc tracks - see function IsSelected(track)
  // 1. Dump OFFline -> HLT tracks into trees  
  //    a.) tree "offhlt" - counters+ Ofline track + corresponding  Hlt track
  //    b.) tree "offhlt" - counters +Offline track
  // 2. HLT and offline
  //
  // 3. To save the CPU and disk space the input tracks to compare are pt downscaled
  //    see function Isdownscaled();

  Int_t nOffTracks = evoff->GetNumberOfTracks();
  Int_t nHLTTracks = evhlt->GetNumberOfTracks();
  Int_t iev= evoff->GetEventNumberInFile();
  if (nOffTracks==0) return;
  AliESDtrack *offlinetrack=0, *hlttrack=0;
  Double_t predchi2 = 0;  
  //for(Long64_t iev = 120; iev < 125; iev++)
  TObjArray arrayOFF(20000);
  TObjArray arrayHLT(20000);
  AliExternalTrackParam  inner2; //working track ref
  //
  //
  //
  Int_t nselOFF = 0;
  Int_t nselHLT = 0;
  // filter offline
  for(Int_t ioff = 0; ioff < nOffTracks; ioff++){
    offlinetrack = evoff->GetTrack(ioff);
    if (!IsSelected(offlinetrack)) continue;
    arrayOFF.AddAt(offlinetrack,nselOFF);	
    nselOFF++;
  }
  //filter hlt
  for(Int_t ihlt = 0; ihlt< nHLTTracks; ihlt++){
    hlttrack = evhlt->GetTrack(ihlt);
    if (!IsSelected(hlttrack)) continue;
    arrayHLT.AddAt(hlttrack,nselHLT);	
    nselHLT++;
  }
  printf("%d\t%d\t%d\t%d\n", nOffTracks, nselOFF,  nHLTTracks, nselHLT);
  //
  // - compare offline - hlt
  TStopwatch timer;
  //
  for (Int_t ioff = 0; ioff < nselOFF; ioff++){
    offlinetrack = (AliESDtrack*)arrayOFF.At(ioff);
    Float_t ncl21off= offlinetrack->GetTPCClusterInfo(2,1);
    Float_t ncl20off= offlinetrack->GetTPCClusterInfo(2,0);
    if (IsDownscaled(offlinetrack)) continue;
    Int_t counter=0;
    //
    for(Int_t ihlt = 0; ihlt < nselHLT; ihlt++){
      hlttrack = (AliESDtrack*)arrayHLT.At(ihlt);
      Bool_t close = AreTracksCloseFast(offlinetrack,hlttrack,evhlt,inner2);	   
      if (close){
	Float_t ncl21hlt= hlttrack->GetTPCClusterInfo(2,1);
	Float_t ncl20hlt= hlttrack->GetTPCClusterInfo(2,0);
	counter++;
	predchi2=offlinetrack->GetInnerParam()->GetPredictedChi2(&inner2);
	(*pcstream)<<"offhlt"<<
	  //multiplicity
	  "nHLT="<<nselHLT<<
	  "nOFF="<<nselOFF<<
	  "ncl20off="<<ncl20off<<
	  "ncl21off="<<ncl21off<<	
	  "ncl20hlt="<<ncl20hlt<<
	  "ncl21hlt="<<ncl21hlt<<	       
	  "counter="<<counter<<
	  //
	  "chi2="<<predchi2<<
	  "track1.="<<offlinetrack<<
	  "track2.="<<hlttrack<<
	  "inner2.="<<&inner2<<
	  "\n";
      }
    }
    Bool_t isDownscaled=IsDownscaled(offlinetrack);	
    (*pcstream)<<"offhlt0"<<
      "iev="<<iev<<	       
      "nHLT="<<nselHLT<<
      "nOFF="<<nselOFF<<
      //
      "ncl20="<<ncl20off<<
      "ncl21="<<ncl21off<<
      "track.="<<offlinetrack<<
      "counter="<<counter<<
      "isDownscaled="<<isDownscaled<<
      "\n";
  }
  timer.Print();
  //
  for(Int_t ihlt = 0; ihlt < nselHLT; ihlt++){
    Int_t counter=0;      
    hlttrack = (AliESDtrack*)arrayHLT.At(ihlt);	
    Float_t ncl21hlt= hlttrack->GetTPCClusterInfo(2,1);
    Float_t ncl20hlt= hlttrack->GetTPCClusterInfo(2,0);
    if (IsDownscaled(hlttrack)) continue;      
    for (Int_t ioff = 0; ioff < nselOFF; ioff++){
      offlinetrack = (AliESDtrack*)arrayOFF.At(ioff);	
      Bool_t close = AreTracksCloseFast(hlttrack,offlinetrack,evhlt,inner2);	   
      if (close){
	Float_t ncl21off= offlinetrack->GetTPCClusterInfo(2,1);
	Float_t ncl20off= offlinetrack->GetTPCClusterInfo(2,0);
	predchi2=hlttrack->GetInnerParam()->GetPredictedChi2(&inner2);
	counter++;
	(*pcstream)<<"hltoff"<<
	  //multiplicity
	  "nHLT="<<nselHLT<<
	  "nOFF="<<nselOFF<<
	  "ncl20off="<<ncl20off<<
	  "ncl21off="<<ncl21off<<
	  "ncl20hlt="<<ncl20hlt<<
	  "ncl21hlt="<<ncl21hlt<<	       
	  "counter="<<counter<<
	  //
	  "chi2="<<predchi2<<
	  "track2.="<<offlinetrack<<
	  "track1.="<<hlttrack<<
	  "inner2.="<<&inner2<<
	  "\n";
      }
    }
    Bool_t isDownscaled=IsDownscaled(hlttrack);	
    (*pcstream)<<"hltoff0"<<
      "iev="<<iev<<	       
      "nHLT="<<nselHLT<<
      "nOFF="<<nselOFF<<
      //
      "ncl20="<<ncl20hlt<<
      "ncl21="<<ncl21hlt<<
      "track.="<<hlttrack<<
      "counter="<<counter<<
      "isDownscaled="<<isDownscaled<<
      "\n";
  } 
  timer.Print();
}

Bool_t IsSelected(AliESDtrack *track){
  //
  // modified by:
  // philipp.luettig@cern.ch
  // 
  //
  if (track->IsOn(0x40)==0) return kFALSE;                // Refit
  if (TMath::Abs(track->GetTgl())>1.1)  return kFALSE;    // tangent lambda
  if (track->GetTPCClusterInfo(2,1)<50) return kFALSE;    // cluster information number of crossed rows >50
  return kTRUE;
}

Bool_t IsDownscaled(AliESDtrack *track){
  //
  // Downscale randomly low pt tracks
  //
  //return kFALSE;
  Double_t scalempt= TMath::Min(1./TMath::Abs(track->GetParameter()[4]),10.);
  if (TMath::Exp(2*scalempt)<gptdownscale*gRandom->Rndm()) return kTRUE;
  return kFALSE;
}


Bool_t AreTracksCloseFast(AliESDtrack *track1, AliESDtrack *track2, AliESDEvent * event, AliExternalTrackParam &inner2){
  //
  // 
  // Fast comparison uning track param close to the prim vertex and at the inner wall of the TPC
  //
  // 1. Fast cut on the invariant (under rotation) variable P3 and P4 (means pz/pt and 1/pt)
  // 2. Slower cuts - parameters at the entrance of the TPC (tracks to be propagated)
  //
  // In case the tracks are relativelaly close -the inner2 parameters are created
  //                                           -track 2 propagated and rotated to the same position as the track1 
  // 
  const Double_t absCut[5] ={5,10, 0.02,0.02,1};   // abs cut values  
  const Double_t pullCut[5]={6,100,6,   100,6};    // pull cut values
  //   *              "External" track parametrisation class                       *
  //   *                                                                           *
  //   *      external param0:   local Y-coordinate of a track (cm)                *
  //   *      external param1:   local Z-coordinate of a track (cm)                *
  //   *      external param2:   local sine of the track momentum azimuthal angle  *
  //   *      external param3:   tangent of the track momentum dip angle           *
  //   *      external param4:   1/pt (1/(GeV/c))                                  *

  //
  // 
  const Double_t kTglCut=0.1;    
  const Double_t k1PtCut=0.5;
  //  const Double_t kAlphaCut=0.2;
  const Double_t kTglCutSigma=10;
  const Double_t k1PtCutSigma=10;
  //
  //
  if(!track1) return kFALSE;
  if(!track2) return kFALSE;  
  const Double_t *param1 = track1->GetParameter();
  const Double_t *param2 = track2->GetParameter();
  const Double_t *param1I = track1->GetInnerParam()->GetParameter();
  const Double_t *param2I = track2->GetInnerParam()->GetParameter();
  const Double_t *covar1 = track1->GetCovariance();
  const Double_t *covar2 = track2->GetCovariance();
  const Double_t *covar1I = track1->GetInnerParam()->GetCovariance();
  const Double_t *covar2I = track2->GetInnerParam()->GetCovariance();
  //
  if (TMath::Abs(param1[3]-param2[3])>kTglCut) return kFALSE;	
  //if (TMath::Abs(param1[4]-param2[4])>k1PtCut) return kFALSE;	
  if (TMath::Abs(param1I[3]-param2I[3])>kTglCut) return kFALSE;	
  if (TMath::Abs(param1I[4]-param2I[4])>k1PtCut) return kFALSE;	
  Double_t dalpha = TMath::Abs(track1->GetAlpha()-track2->GetAlpha());
  if (dalpha>TMath::Pi()) dalpha-=TMath::Abs(dalpha-TMath::TwoPi());
  //if (dalpha > kAlphaCut) return kFALSE;

  //
  Int_t index22=track1->GetIndex(2,2);	       
  Int_t index33=track1->GetIndex(3,3);	       
  Int_t index44=track1->GetIndex(4,4);	       
  if (TMath::Abs(param1[3]-param2[3])/TMath::Sqrt(TMath::Max(covar1[index33],covar2[index33]))>kTglCutSigma) return kFALSE;
  //if (TMath::Abs(param1[4]-param2[4])/TMath::Sqrt(TMath::Max(covar1[index44],covar2[index44]))>k1PtCutSigma) return kFALSE;  
  if (TMath::Abs(param1I[3]-param2I[3])/TMath::Sqrt(TMath::Max(covar1I[index33],covar2I[index33]))>kTglCutSigma) return kFALSE;
  if (TMath::Abs(param1I[4]-param2I[4])/TMath::Sqrt(TMath::Max(covar1I[index44],covar2I[index44]))>k1PtCutSigma) return kFALSE;  
  if (TMath::Abs(dalpha)/TMath::Sqrt(TMath::Max(covar1[index22],covar2[index22]))>k1PtCutSigma) return kFALSE;  
  //
  // 2. Slow cuts on the paramters at the entrance of the TPC
  //
  inner2=*(track2->GetInnerParam());
  inner2.Rotate(track1->GetInnerParam()->GetAlpha());
  inner2.PropagateTo(track1->GetInnerParam()->GetX(),  event->GetMagneticField());
  const Double_t *pinner2  = inner2.GetParameter();
  const Double_t *pcovar2  = inner2.GetCovariance();
  //
  Bool_t isOK = kTRUE;
  for (Int_t ipar=0; ipar<5; ipar++){
    // 
    Int_t index=track1->GetIndex(ipar,ipar);	      
    if (TMath::Abs(pinner2[ipar]-param1I[ipar])>absCut[ipar]) isOK=kFALSE;
    if (TMath::Abs(pinner2[ipar]-param1I[ipar])>pullCut[ipar]*TMath::Sqrt(TMath::Max(covar1I[index],pcovar2[index]))) isOK=kFALSE;
  }
  if (!isOK) return kFALSE; 
  return kTRUE;
}


void MakeChain(){
  //
  //
  AliXRDPROOFtoolkit toolkit;
  chainOFFHLT= toolkit.MakeChainRandom("dumpHLTOFFLINE.list","offhlt",0,100);
  chainOFFHLT0=toolkit.MakeChainRandom("dumpHLTOFFLINE.list","offhlt0",0,100);
  chainHLTOFF=toolkit.MakeChainRandom("dumpHLTOFFLINE.list","hltoff",0,100);
  chainHLTOFF0=toolkit.MakeChainRandom("dumpHLTOFFLINE.list","hltoff0",0,100);

 
}

void DrawDiffPt(){
  //
  // Draw difference between the HLT and offline tracks
  //
  TCut cut="sqrt(chi2)<10&&ncl21off>120";
  TCut cutNoiseEvent = "abs(nHLT/nOFF-1)<0.2";   //mainly laser events
  //
  // 1. check the edge effect 1/pt resolution TPC only pull
  // ...
  chainOFFHLT->Draw("(track1.fIp.fP[4]-track2.fIp.fP[4])/sqrt(max(track1.fIp.fC[14],track2.fIp.fC[14])):sign(inner2.fP[4])*inner2.fP[0]/inner2.fX>>hisTPCEdge(50,-0.18,0.18,100,-6,6)",cut+"abs(track1.fP[4])<0.25","colz",200000);
  /*
    hisTPCEdge->FitSlicesY();
    hisTPCEdge_2->GetXaxis()->SetTitle("q*ly/lx");
    hisTPCEdge_2->GetYaxis()->SetTitle("#Delta_{1/pt}/#sigma_{1/pt}");
    hisTPCEdge_2->Draw();
  */
  // 2. check the edge effect 1/pt resolution combined 
  chainOFFHLT->Draw("(track1.fP[4]-track2.fP[4])/sqrt(max(track1.fC[14],track2.fC[14])):sign(inner2.fP[4])*inner2.fP[0]/inner2.fX>>hisCombEdge(50,-0.18,0.18,100,-6,6)",cut+"abs(track1.fP[4])<0.25","colz",200000);
  /*
    hisCombEdge->FitSlicesY();
    hisCombEdge_2->Draw();
  */
  // 3. Combined momentum resolution as function of the inverse moment
  chainOFFHLT->Draw("(track1.fIp.fP[4]-track2.fIp.fP[4])/sqrt(max(track1.fIp.fC[14],track2.fIp.fC[14])):abs(track1.fP[4])>>hisTPCP4(20,-0.0,1,100,-6,6)",cut+"abs(track1.fP[4])<1","colz",200000);
  /*
    hisTPCP4->FitSlicesY();  
    hisTPCP4_2->GetXaxis()->SetTitle("1/p_{t} (1/GeV))");
    hisTPCP4_2->GetYaxis()->SetTitle("#Delta_{1/pt}/#sigma_{1/pt}");
    hisTPCP4_2->Draw();
  */
  // 4. Combined momentum resolution as function of the inverse moment
  chainOFFHLT->Draw("(track1.fP[4]-track2.fP[4])/sqrt(max(track1.fC[14],track2.fC[14])):abs(track1.fP[4])>>hisCombP4(20,-0.0,1,100,-6,6)",cut+"abs(track1.fP[4])<1","colz",200000);
  /*
    hisCombP4->FitSlicesY();
    hisCombP4_2->Draw();
  */

}
//
void DrawDiffEff(){
  //
  //
  //
  TCut cutEff = "abs(track.fIp.fP[4])<5&&abs(track.fIp.fP[1])<90";
  TCut cutNoiseEvent = "abs(nHLT/nOFF-1)<0.2"; // HLT cluster finder more sensitive to the noise
  //
  // 
  //
  chainOFFHLT0->Draw("counter==0:(nOFF+nHLT)/2.>>effOccuOFFHLT(20,0,8000)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1&&ncl21>120","prof",50000);
  chainHLTOFF0->Draw("counter==0:(nOFF+nHLT)/2.>>effOccuHLTOFF(20,0,8000)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1&&ncl21>120","prof",50000);

 chainOFFHLT0->Draw("counter==0:track.fTPCncls>>effNCLOFFHLT(40,0,160)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1","prof",50000);
 chainHLTOFF0->Draw("counter==0:track.fTPCncls>>effNCLHLTOFF(40,0,160)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1","prof",50000);
 /*
   effNCLOFFHLT->SetMarkerStyle(25);
   effNCLHLTOFF->SetMarkerStyle(25);
   effNCLOFFHLT->SetMarkerColor(2);
   effNCLHLTOFF->SetMarkerColor(4);
   effNCLOFFHLT->Draw();
   effNCLHLTOFF->Draw("same");
  */


  chainHLTOFF0->Draw("counter==0:sign(track.fIp.fP[4])*track.fIp.fP[0]/track.fIp.fX>>profTPCEdge(50,-0.18,0.18)",cutNoiseEvent+cutEff+"abs(track.fP[4])<0.25","prof",50000);

  chainOFFHLT0->Draw("counter==0:sign(track.fIp.fP[4])*track.fIp.fP[0]/track.fIp.fX>>profTPCEdge(50,-0.18,0.18)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1","prof",50000);

  chainOFFHLT0->Draw("track.fTPCncls:sign(track.fIp.fP[4])*track.fIp.fP[0]/track.fIp.fX>>profTPCEdge(50,-0.18,0.18)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1","prof",50000);


}

/*
  This is  shell script real example  to submit jobs for the track comparison: 
  //
  //
  rm -rf list*
  rm -rf dirlist*
  wdir=`pwd`
  split offline.list --lines=50 list -d
  for a in `ls list*`; do
    mkdir dir$a
    cd dir$a
    mv ../$a compare.list
    bsub -q proof -oo outcompare.log aliroot -b -q $ALICE_ROOT/PWGPP/TPC/macroscompareTracks.C+
    cd $wdir
  done;

  find  `pwd`/  | grep .root > dumpHLTOFFLINE.list

*/
Commit	Line	Data
78a54db2	1	/*
	2
	3	//
	4	// Macro to compare tracks used 2 different reconstruction algorithms (Offline and HLT)
	5	// (In general the approach can be used for any pair of the esds - track containers)
	6	//
	7	// Main functions - Dump content of HLT and OFFLINE esd files - full esd tracks
	8	// 1. Dump to the tree the pair of close tracks ( n sigma distance at the TPC entrance)
	9	// 2. Dump to the tree the track and corrsponding counter of tracks from other container
	10	// Compare two events - track by track comparison
	11	// 0. Filter "Good" tpc tracks - see function IsSelected(track)
	12	// 1. Dump OFFline -> HLT tracks into trees
	13	// a.) tree "offhlt" - counters+ Ofline track + corresponding Hlt track
	14	// b.) tree "offhlt0" - offline track + hlt track counter
	15	// 2. HLT and offline
	16	//
	17	// 3. To save the CPU and disk space the input tracks to compare are pt downscaled
	18	// see function IsDownscaled();
	19	//
	20	// Visualization part is currently work in progress
cdd9da81	21	//
	22	// Responsible:
	23	// dumping part: marian.ivanov@cern.ch
	24	// visualization part:
78a54db2	25	//
	26	Example usage:
	27	aliroot -b -q $ALICE_ROOT/TPC/macros/compareTracks.C+
	28	//
	29	//
	30	.x $HOME/rootlogon.C
	31	gSystem->Load("libANALYSIS");
	32	gSystem->Load("libANALYSISalice");
	33	gSystem->Load("libPWG0base");
	34	gSystem->Load("libPWG0dep");
	35	gSystem->Load("libPWG0selectors");
2bfe5463	36	.L $ALICE_ROOT/PWGPP/TPC/macros/compareTracks.C+
78a54db2	37	//
	38	compareTracks("compare.list");
	39	//
	40	TFile f("dump.root");
	41	*/
	42	#include <exception>
	43	#include <fstream>
	44	#include "TFile.h"
	45	#include "TTree.h"
	46	#include "TChain.h"
	47	#include "TCut.h"
	48	#include "TStopwatch.h"
	49	#include "TSystem.h"
	50	#include "AliESDEvent.h"
	51	#include "TTreeStream.h"
	52	#include "TRandom.h"
	53	#include "AliXRDPROOFtoolkit.h"
	54	#include "AliSysInfo.h"
	55
	56	//
	57	//
	58	//
	59	void CompareFile(TTreeSRedirector pcstream, TFile fOff, TFile *fHLT);
	60	void CompareEvents(AliESDEvent evoff, AliESDEvent evhlt, TTreeSRedirector *pcstream);
	61	Bool_t AreTracksCloseFast(AliESDtrack track1, AliESDtrack track2, AliESDEvent * event, AliExternalTrackParam &inner2);
	62	Bool_t IsSelected(AliESDtrack *track);
	63	Bool_t IsDownscaled(AliESDtrack *track);
	64	//
	65	Double_t gptdownscale=100; //pt downscale parameter
	66	//
	67	TTree * chainOFFHLT=0;
	68	TTree * chainOFFHLT0=0;
	69	TTree * chainHLTOFF=0;
	70	TTree * chainHLTOFF0=0;
	71
9a19fee3	72	void compareTracks(const char * flistOFF="compareOFF.list", const char * flistHLT="compareHLT.list", Double_t downscale=10000){
	73	//
	74	// Compare HLT and OFFLINE esd files
	75	// Input:
	76	// flist - the ascii files with the filename - ESD form the offline
	77	// - hlt path constructed from the offline ESD file repalcing the OFFLINE with HLT
	78	//
	79	gptdownscale=downscale;
	80	TTreeSRedirector *pcstream = new TTreeSRedirector("dump.root");
	81	ifstream inOFF;
	82	inOFF.open(flistOFF);
	83	ifstream inHLT;
	84	inHLT.open(flistHLT);
	85
78a54db2	86
9a19fee3	87	// Read the input list of files and add them to the chain
	88	TString currentFileOff;
	89	TString currentFileHlt;
	90	{ while( (inOFF.good()) && (inHLT.good())) {
	91	inOFF >> currentFileOff;
	92	inHLT >> currentFileHlt;
	93	//currentFileHlt=currentFileOff;
	94	//currentFileHlt.ReplaceAll("/OFFLINE/","/HLT/");
	95	printf("%s\n%s\n\n",currentFileOff.Data(), currentFileHlt.Data());
	96	TFile *fOff = TFile::Open(currentFileOff.Data(),"READ");
	97	TFile *fHlt = TFile::Open(currentFileHlt.Data(),"READ");
	98	if ( fOff && fHlt){
	99	CompareFile(pcstream, fOff, fHlt);
	100	}
	101	delete fOff;
	102	delete fHlt;
78a54db2	103	}}
9a19fee3	104	delete pcstream;
	105	/*
	106	* TTreeSRedirector *pcstream = new TTreeSRedirector("dump.root");
	107	* TFile *fOff = TFile::Open("/lustre/alice/mknichel/cpass1/2011-12-13_0154/OFFLINE/167693/11000167693001.11/AliESDs.root","READ");
	108	* TFile *fHLT = TFile::Open("/lustre/alice/mknichel/cpass1/2011-12-13_0154/HLT/167693/11000167693001.11/AliESDs.root","READ");
	109	* CompareFile(pcstream, fOff, fHLT);
	110	* delete pcstream;
	111	* */
78a54db2	112	}
	113
	114	void CompareFile(TTreeSRedirector pcstream, TFile fOff, TFile *fHLT){
	115	//
	116	// Compare 2 esd files - fOff and fHLT
	117	// Pairs of tracks dumpt in the trees - stored in the filee connected to the pcstream
	118	//
	119	TTree tOff = (TTree)fOff->Get("esdTree");
	120	TTree tHLT = (TTree)fHLT->Get("esdTree");
	121	//
	122	Int_t nevents=tOff->GetEntries();
	123	//
	124	AliESDEvent *evoff = new AliESDEvent();
	125	AliESDEvent *evhlt = new AliESDEvent();
	126	evoff->ReadFromTree(tOff);
	127	evhlt->ReadFromTree(tHLT);
	128	//
	129	for(Long64_t iev = 0; iev < nevents; iev++){
	130	tOff->GetEvent(iev);
	131	tHLT->GetEvent(iev);
	132	CompareEvents(evoff,evhlt,pcstream);
	133	AliSysInfo::AddStamp(fOff->GetName(),iev);
	134	}
	135	delete evoff;
	136	delete evhlt;
	137	}
	138
	139
	140
	141	void CompareEvents(AliESDEvent evoff, AliESDEvent evhlt, TTreeSRedirector *pcstream){
	142	//
	143	// Main function to dump events:
	144	//
	145	// compare two events - track by track comparison
	146	// 0. Filter "Good" tpc tracks - see function IsSelected(track)
	147	// 1. Dump OFFline -> HLT tracks into trees
	148	// a.) tree "offhlt" - counters+ Ofline track + corresponding Hlt track
	149	// b.) tree "offhlt" - counters +Offline track
	150	// 2. HLT and offline
	151	//
	152	// 3. To save the CPU and disk space the input tracks to compare are pt downscaled
	153	// see function Isdownscaled();
	154
	155	Int_t nOffTracks = evoff->GetNumberOfTracks();
	156	Int_t nHLTTracks = evhlt->GetNumberOfTracks();
	157	Int_t iev= evoff->GetEventNumberInFile();
	158	if (nOffTracks==0) return;
	159	AliESDtrack offlinetrack=0, hlttrack=0;
	160	Double_t predchi2 = 0;
	161	//for(Long64_t iev = 120; iev < 125; iev++)
	162	TObjArray arrayOFF(20000);
	163	TObjArray arrayHLT(20000);
	164	AliExternalTrackParam inner2; //working track ref
	165	//
	166	//
	167	//
	168	Int_t nselOFF = 0;
	169	Int_t nselHLT = 0;
	170	// filter offline
	171	for(Int_t ioff = 0; ioff < nOffTracks; ioff++){
	172	offlinetrack = evoff->GetTrack(ioff);
	173	if (!IsSelected(offlinetrack)) continue;
	174	arrayOFF.AddAt(offlinetrack,nselOFF);
	175	nselOFF++;
176	}
177	//filter hlt
178	for(Int_t ihlt = 0; ihlt< nHLTTracks; ihlt++){
179	hlttrack = evhlt->GetTrack(ihlt);
180	if (!IsSelected(hlttrack)) continue;
181	arrayHLT.AddAt(hlttrack,nselHLT);
182	nselHLT++;
183	}
184	printf("%d\t%d\t%d\t%d\n", nOffTracks, nselOFF, nHLTTracks, nselHLT);
185	//
186	// - compare offline - hlt
187	TStopwatch timer;
188	//
189	for (Int_t ioff = 0; ioff < nselOFF; ioff++){
190	offlinetrack = (AliESDtrack*)arrayOFF.At(ioff);
191	Float_t ncl21off= offlinetrack->GetTPCClusterInfo(2,1);
192	Float_t ncl20off= offlinetrack->GetTPCClusterInfo(2,0);
193	if (IsDownscaled(offlinetrack)) continue;
194	Int_t counter=0;
195	//
196	for(Int_t ihlt = 0; ihlt < nselHLT; ihlt++){
197	hlttrack = (AliESDtrack*)arrayHLT.At(ihlt);
198	Bool_t close = AreTracksCloseFast(offlinetrack,hlttrack,evhlt,inner2);
199	if (close){
200	Float_t ncl21hlt= hlttrack->GetTPCClusterInfo(2,1);
201	Float_t ncl20hlt= hlttrack->GetTPCClusterInfo(2,0);
202	counter++;
203	predchi2=offlinetrack->GetInnerParam()->GetPredictedChi2(&inner2);
204	(*pcstream)<<"offhlt"<<
205	//multiplicity
206	"nHLT="<<nselHLT<<
207	"nOFF="<<nselOFF<<
208	"ncl20off="<<ncl20off<<
209	"ncl21off="<<ncl21off<<
210	"ncl20hlt="<<ncl20hlt<<
211	"ncl21hlt="<<ncl21hlt<<
212	"counter="<<counter<<
213	//
214	"chi2="<<predchi2<<
215	"track1.="<<offlinetrack<<
216	"track2.="<<hlttrack<<
217	"inner2.="<<&inner2<<
218	"\n";
219	}
220	}
221	Bool_t isDownscaled=IsDownscaled(offlinetrack);
222	(*pcstream)<<"offhlt0"<<
223	"iev="<<iev<<
224	"nHLT="<<nselHLT<<
225	"nOFF="<<nselOFF<<
226	//
227	"ncl20="<<ncl20off<<
228	"ncl21="<<ncl21off<<
229	"track.="<<offlinetrack<<
230	"counter="<<counter<<
231	"isDownscaled="<<isDownscaled<<
232	"\n";
233	}
234	timer.Print();
235	//
236	for(Int_t ihlt = 0; ihlt < nselHLT; ihlt++){
237	Int_t counter=0;
238	hlttrack = (AliESDtrack*)arrayHLT.At(ihlt);
239	Float_t ncl21hlt= hlttrack->GetTPCClusterInfo(2,1);
240	Float_t ncl20hlt= hlttrack->GetTPCClusterInfo(2,0);
241	if (IsDownscaled(hlttrack)) continue;
242	for (Int_t ioff = 0; ioff < nselOFF; ioff++){
243	offlinetrack = (AliESDtrack*)arrayOFF.At(ioff);
244	Bool_t close = AreTracksCloseFast(hlttrack,offlinetrack,evhlt,inner2);
245	if (close){
246	Float_t ncl21off= offlinetrack->GetTPCClusterInfo(2,1);
247	Float_t ncl20off= offlinetrack->GetTPCClusterInfo(2,0);
248	predchi2=hlttrack->GetInnerParam()->GetPredictedChi2(&inner2);
249	counter++;
250	(*pcstream)<<"hltoff"<<
251	//multiplicity
252	"nHLT="<<nselHLT<<
253	"nOFF="<<nselOFF<<
254	"ncl20off="<<ncl20off<<
255	"ncl21off="<<ncl21off<<
256	"ncl20hlt="<<ncl20hlt<<
257	"ncl21hlt="<<ncl21hlt<<
258	"counter="<<counter<<
259	//
260	"chi2="<<predchi2<<
261	"track2.="<<offlinetrack<<
262	"track1.="<<hlttrack<<
263	"inner2.="<<&inner2<<
264	"\n";
265	}
266	}
267	Bool_t isDownscaled=IsDownscaled(hlttrack);
268	(*pcstream)<<"hltoff0"<<
269	"iev="<<iev<<
270	"nHLT="<<nselHLT<<
271	"nOFF="<<nselOFF<<
272	//
273	"ncl20="<<ncl20hlt<<
274	"ncl21="<<ncl21hlt<<
275	"track.="<<hlttrack<<
276	"counter="<<counter<<
277	"isDownscaled="<<isDownscaled<<
278	"\n";
279	}
280	timer.Print();
281	}
282
283	Bool_t IsSelected(AliESDtrack *track){
284	//
cdd9da81	285	// modified by:
cdd9da81	286	// philipp.luettig@cern.ch
78a54db2	287	//
	288	//
	289	if (track->IsOn(0x40)==0) return kFALSE; // Refit
	290	if (TMath::Abs(track->GetTgl())>1.1) return kFALSE; // tangent lambda
	291	if (track->GetTPCClusterInfo(2,1)<50) return kFALSE; // cluster information number of crossed rows >50
	292	return kTRUE;
	293	}
	294
	295	Bool_t IsDownscaled(AliESDtrack *track){
	296	//
	297	// Downscale randomly low pt tracks
	298	//
	299	//return kFALSE;
	300	Double_t scalempt= TMath::Min(1./TMath::Abs(track->GetParameter()[4]),10.);
	301	if (TMath::Exp(2scalempt)<gptdownscalegRandom->Rndm()) return kTRUE;
	302	return kFALSE;
	303	}
	304
	305
	306	Bool_t AreTracksCloseFast(AliESDtrack track1, AliESDtrack track2, AliESDEvent * event, AliExternalTrackParam &inner2){
	307	//
	308	//
	309	// Fast comparison uning track param close to the prim vertex and at the inner wall of the TPC
	310	//
cdd9da81	311	// 1. Fast cut on the invariant (under rotation) variable P3 and P4 (means pz/pt and 1/pt)
78a54db2	312	// 2. Slower cuts - parameters at the entrance of the TPC (tracks to be propagated)
	313	//
	314	// In case the tracks are relativelaly close -the inner2 parameters are created
	315	// -track 2 propagated and rotated to the same position as the track1
cdd9da81	316	//
78a54db2	317	const Double_t absCut[5] ={5,10, 0.02,0.02,1}; // abs cut values
78a54db2	318	const Double_t pullCut[5]={6,100,6, 100,6}; // pull cut values
cdd9da81	319	// * "External" track parametrisation class *
	320	// * *
	321	// * external param0: local Y-coordinate of a track (cm) *
	322	// * external param1: local Z-coordinate of a track (cm) *
	323	// * external param2: local sine of the track momentum azimuthal angle *
	324	// * external param3: tangent of the track momentum dip angle *
	325	// * external param4: 1/pt (1/(GeV/c)) *
	326
78a54db2	327	//
	328	//
	329	const Double_t kTglCut=0.1;
	330	const Double_t k1PtCut=0.5;
	331	// const Double_t kAlphaCut=0.2;
	332	const Double_t kTglCutSigma=10;
	333	const Double_t k1PtCutSigma=10;
	334	//
	335	//
	336	if(!track1) return kFALSE;
	337	if(!track2) return kFALSE;
	338	const Double_t *param1 = track1->GetParameter();
	339	const Double_t *param2 = track2->GetParameter();
	340	const Double_t *param1I = track1->GetInnerParam()->GetParameter();
	341	const Double_t *param2I = track2->GetInnerParam()->GetParameter();
	342	const Double_t *covar1 = track1->GetCovariance();
	343	const Double_t *covar2 = track2->GetCovariance();
	344	const Double_t *covar1I = track1->GetInnerParam()->GetCovariance();
	345	const Double_t *covar2I = track2->GetInnerParam()->GetCovariance();
	346	//
	347	if (TMath::Abs(param1[3]-param2[3])>kTglCut) return kFALSE;
	348	//if (TMath::Abs(param1[4]-param2[4])>k1PtCut) return kFALSE;
	349	if (TMath::Abs(param1I[3]-param2I[3])>kTglCut) return kFALSE;
	350	if (TMath::Abs(param1I[4]-param2I[4])>k1PtCut) return kFALSE;
	351	Double_t dalpha = TMath::Abs(track1->GetAlpha()-track2->GetAlpha());
	352	if (dalpha>TMath::Pi()) dalpha-=TMath::Abs(dalpha-TMath::TwoPi());
	353	//if (dalpha > kAlphaCut) return kFALSE;
	354
	355	//
	356	Int_t index22=track1->GetIndex(2,2);
	357	Int_t index33=track1->GetIndex(3,3);
	358	Int_t index44=track1->GetIndex(4,4);
	359	if (TMath::Abs(param1[3]-param2[3])/TMath::Sqrt(TMath::Max(covar1[index33],covar2[index33]))>kTglCutSigma) return kFALSE;
	360	//if (TMath::Abs(param1[4]-param2[4])/TMath::Sqrt(TMath::Max(covar1[index44],covar2[index44]))>k1PtCutSigma) return kFALSE;
	361	if (TMath::Abs(param1I[3]-param2I[3])/TMath::Sqrt(TMath::Max(covar1I[index33],covar2I[index33]))>kTglCutSigma) return kFALSE;
	362	if (TMath::Abs(param1I[4]-param2I[4])/TMath::Sqrt(TMath::Max(covar1I[index44],covar2I[index44]))>k1PtCutSigma) return kFALSE;
	363	if (TMath::Abs(dalpha)/TMath::Sqrt(TMath::Max(covar1[index22],covar2[index22]))>k1PtCutSigma) return kFALSE;
	364	//
	365	// 2. Slow cuts on the paramters at the entrance of the TPC
	366	//
	367	inner2=*(track2->GetInnerParam());
	368	inner2.Rotate(track1->GetInnerParam()->GetAlpha());
	369	inner2.PropagateTo(track1->GetInnerParam()->GetX(), event->GetMagneticField());
	370	const Double_t *pinner2 = inner2.GetParameter();
	371	const Double_t *pcovar2 = inner2.GetCovariance();
	372	//
	373	Bool_t isOK = kTRUE;
	374	for (Int_t ipar=0; ipar<5; ipar++){
	375	//
	376	Int_t index=track1->GetIndex(ipar,ipar);
	377	if (TMath::Abs(pinner2[ipar]-param1I[ipar])>absCut[ipar]) isOK=kFALSE;
	378	if (TMath::Abs(pinner2[ipar]-param1I[ipar])>pullCut[ipar]*TMath::Sqrt(TMath::Max(covar1I[index],pcovar2[index]))) isOK=kFALSE;
	379	}
	380	if (!isOK) return kFALSE;
	381	return kTRUE;
	382	}
	383
	384
	385	void MakeChain(){
	386	//
	387	//
	388	AliXRDPROOFtoolkit toolkit;
	389	chainOFFHLT= toolkit.MakeChainRandom("dumpHLTOFFLINE.list","offhlt",0,100);
	390	chainOFFHLT0=toolkit.MakeChainRandom("dumpHLTOFFLINE.list","offhlt0",0,100);
391	chainHLTOFF=toolkit.MakeChainRandom("dumpHLTOFFLINE.list","hltoff",0,100);
392	chainHLTOFF0=toolkit.MakeChainRandom("dumpHLTOFFLINE.list","hltoff0",0,100);
cdd9da81	393
cdd9da81	394
78a54db2	395	}
	396
	397	void DrawDiffPt(){
	398	//
	399	// Draw difference between the HLT and offline tracks
	400	//
	401	TCut cut="sqrt(chi2)<10&&ncl21off>120";
cdd9da81	402	TCut cutNoiseEvent = "abs(nHLT/nOFF-1)<0.2"; //mainly laser events
78a54db2	403	//
78a54db2	404	// 1. check the edge effect 1/pt resolution TPC only pull
cdd9da81	405	// ...
cdd9da81	406	chainOFFHLT->Draw("(track1.fIp.fP[4]-track2.fIp.fP[4])/sqrt(max(track1.fIp.fC[14],track2.fIp.fC[14])):sign(inner2.fP[4])*inner2.fP[0]/inner2.fX>>hisTPCEdge(50,-0.18,0.18,100,-6,6)",cut+"abs(track1.fP[4])<0.25","colz",200000);
78a54db2	407	/*
cdd9da81	408	hisTPCEdge->FitSlicesY();
	409	hisTPCEdge_2->GetXaxis()->SetTitle("q*ly/lx");
	410	hisTPCEdge_2->GetYaxis()->SetTitle("#Delta_{1/pt}/#sigma_{1/pt}");
	411	hisTPCEdge_2->Draw();
78a54db2	412	*/
	413	// 2. check the edge effect 1/pt resolution combined
	414	chainOFFHLT->Draw("(track1.fP[4]-track2.fP[4])/sqrt(max(track1.fC[14],track2.fC[14])):sign(inner2.fP[4])*inner2.fP[0]/inner2.fX>>hisCombEdge(50,-0.18,0.18,100,-6,6)",cut+"abs(track1.fP[4])<0.25","colz",200000);
	415	/*
	416	hisCombEdge->FitSlicesY();
	417	hisCombEdge_2->Draw();
	418	*/
	419	// 3. Combined momentum resolution as function of the inverse moment
	420	chainOFFHLT->Draw("(track1.fIp.fP[4]-track2.fIp.fP[4])/sqrt(max(track1.fIp.fC[14],track2.fIp.fC[14])):abs(track1.fP[4])>>hisTPCP4(20,-0.0,1,100,-6,6)",cut+"abs(track1.fP[4])<1","colz",200000);
	421	/*
cdd9da81	422	hisTPCP4->FitSlicesY();
	423	hisTPCP4_2->GetXaxis()->SetTitle("1/p_{t} (1/GeV))");
	424	hisTPCP4_2->GetYaxis()->SetTitle("#Delta_{1/pt}/#sigma_{1/pt}");
78a54db2	425	hisTPCP4_2->Draw();
	426	*/
	427	// 4. Combined momentum resolution as function of the inverse moment
	428	chainOFFHLT->Draw("(track1.fP[4]-track2.fP[4])/sqrt(max(track1.fC[14],track2.fC[14])):abs(track1.fP[4])>>hisCombP4(20,-0.0,1,100,-6,6)",cut+"abs(track1.fP[4])<1","colz",200000);
	429	/*
	430	hisCombP4->FitSlicesY();
	431	hisCombP4_2->Draw();
	432	*/
	433
	434	}
cdd9da81	435	//
	436	void DrawDiffEff(){
	437	//
	438	//
	439	//
	440	TCut cutEff = "abs(track.fIp.fP[4])<5&&abs(track.fIp.fP[1])<90";
	441	TCut cutNoiseEvent = "abs(nHLT/nOFF-1)<0.2"; // HLT cluster finder more sensitive to the noise
	442	//
	443	//
	444	//
	445	chainOFFHLT0->Draw("counter==0:(nOFF+nHLT)/2.>>effOccuOFFHLT(20,0,8000)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1&&ncl21>120","prof",50000);
	446	chainHLTOFF0->Draw("counter==0:(nOFF+nHLT)/2.>>effOccuHLTOFF(20,0,8000)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1&&ncl21>120","prof",50000);
	447
	448	chainOFFHLT0->Draw("counter==0:track.fTPCncls>>effNCLOFFHLT(40,0,160)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1","prof",50000);
	449	chainHLTOFF0->Draw("counter==0:track.fTPCncls>>effNCLHLTOFF(40,0,160)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1","prof",50000);
	450	/*
	451	effNCLOFFHLT->SetMarkerStyle(25);
	452	effNCLHLTOFF->SetMarkerStyle(25);
	453	effNCLOFFHLT->SetMarkerColor(2);
	454	effNCLHLTOFF->SetMarkerColor(4);
	455	effNCLOFFHLT->Draw();
	456	effNCLHLTOFF->Draw("same");
	457	*/
	458
78a54db2	459
cdd9da81	460	chainHLTOFF0->Draw("counter==0:sign(track.fIp.fP[4])*track.fIp.fP[0]/track.fIp.fX>>profTPCEdge(50,-0.18,0.18)",cutNoiseEvent+cutEff+"abs(track.fP[4])<0.25","prof",50000);
	461
	462	chainOFFHLT0->Draw("counter==0:sign(track.fIp.fP[4])*track.fIp.fP[0]/track.fIp.fX>>profTPCEdge(50,-0.18,0.18)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1","prof",50000);
	463
	464	chainOFFHLT0->Draw("track.fTPCncls:sign(track.fIp.fP[4])*track.fIp.fP[0]/track.fIp.fX>>profTPCEdge(50,-0.18,0.18)",cutNoiseEvent+cutEff+"abs(track.fP[4])<1","prof",50000);
	465
	466
	467	}
78a54db2	468
	469	/*
	470	This is shell script real example to submit jobs for the track comparison:
	471	//
	472	//
	473	rm -rf list*
	474	rm -rf dirlist*
	475	wdir=`pwd`
	476	split offline.list --lines=50 list -d
	477	for a in `ls list*`; do
	478	mkdir dir$a
	479	cd dir$a
	480	mv ../$a compare.list
2bfe5463	481	bsub -q proof -oo outcompare.log aliroot -b -q $ALICE_ROOT/PWGPP/TPC/macroscompareTracks.C+
78a54db2	482	cd $wdir
	483	done;
	484
	485	find `pwd`/ \| grep .root > dumpHLTOFFLINE.list
	486
	487	*/