removing old files
authorjgrosseo <jgrosseo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 25 Feb 2009 10:06:11 +0000 (10:06 +0000)
committerjgrosseo <jgrosseo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 25 Feb 2009 10:06:11 +0000 (10:06 +0000)
PWG0/AliRecoVsGeneCheck.cxx [deleted file]
PWG0/AliRecoVsGeneCheck.h [deleted file]
PWG0/README [deleted file]
PWG0/esdTrackCuts/AliTestESDtrackCutsSelector.cxx [deleted file]
PWG0/esdTrackCuts/AliTestESDtrackCutsSelector.h [deleted file]
PWG0/esdTrackCuts/testESDtrackCuts.C [deleted file]

diff --git a/PWG0/AliRecoVsGeneCheck.cxx b/PWG0/AliRecoVsGeneCheck.cxx
deleted file mode 100644 (file)
index a984376..0000000
+++ /dev/null
@@ -1,143 +0,0 @@
-#include "AliRecoVsGeneCheck.h"
-
-#include <Riostream.h>
-
-//____________________________________________________________________
-ClassImp(AliRecoVsGeneCheck);
-
-//____________________________________________________________________
-AliRecoVsGeneCheck::AliRecoVsGeneCheck() {
-
-  fhVtzZRecoVsMC     = new TH2F("vtx_z_reco_vs_mc","",100,-20,20,100,-20,20);
-
-  fhVtxZRes          = new TH1F("vtx_z_res","",100,0,0.05);
-  fhVtxZResVsZ       = new TH2F("vtx_z_res_vs_z","",100,-20,20,100,0,0.05);
-  fhVtxZResVsNPart   = new TH2F("vtx_z_res_vs_npart","",100,-0.5,99.5,100,0,0.05);
-
-  fhVtxDzNorm        = new TH1F("vtx_dz_norm","",100,-10,10);
-  fhVtxDzNormVsZ     = new TH2F("vtx_dz_norm_vs_z","",100,-20,20,100,-10,10);
-  fhVtxDzNormVsNPart = new TH2F("vtx_dz_norm_vs_npart","",100,-0.5,99.5,100,-10,10);
-
-  fhVtxZMC           = new TH1F("vtx_z_mc","",100,-20,20);
-  fhVtxZReco         = new TH1F("vtx_z_reco","",100,-20,20);
-
-  fhNPart            = new TH1F("n_part","",100,-0.5,99.5);
-
-  fhDPtVsPtVsEta     = new TH3F("dpt_vs_pt_vs_eta","", 30,-1.5,1.5,30,0,3,120,-0.3,0.3);
-  fhDEtaVsPtVsEta    = new TH3F("deta_vs_pt_vs_eta","",30,-1.5,1.5,30,0,3,120,-0.3,0.3);
-  
-  fhVtxZMC  ->Sumw2();
-  fhVtxZReco->Sumw2();
-  fhNPart   ->Sumw2();
-
-  fhVtzZRecoVsMC    ->SetXTitle("vertex z_{mc}"); 
-  fhVtzZRecoVsMC    ->SetYTitle("vertex z_{reco}"); 
-
-  fhVtxZRes         ->SetXTitle("vertex #sigma_{z_{reco}}");
-  fhVtxZResVsZ      ->SetXTitle("vertex z_{mc}");
-  fhVtxZResVsZ      ->SetYTitle("vertex #sigma_{z_{reco}}");            
-  fhVtxZResVsNPart  ->SetXTitle("n part");
-  fhVtxZResVsNPart  ->SetYTitle("vertex #sigma_{z_{reco}}");
-    
-  fhVtxDzNorm       ->SetXTitle("vertex (z_{mc} - z_{reco}) / #sigma_{z_{reco}}");
-  fhVtxDzNormVsZ    ->SetXTitle("vertex z_{mc}");
-  fhVtxDzNormVsZ    ->SetYTitle("vertex (z_{mc} - z_{reco}) / #sigma_{z_{reco}}");
-  fhVtxDzNormVsNPart->SetXTitle("npart");                                   
-  fhVtxDzNormVsNPart->SetYTitle("vertex (z_{mc} - z_{reco}) / #sigma_{z_{reco}}"); 
-
-  fhNPart           ->SetXTitle("n part");             
-  fhVtxZMC          ->SetXTitle("vertex z_{mc}");
-  fhVtxZReco        ->SetXTitle("vertex z_{reco}");
-
-  fhDPtVsPtVsEta ->SetXTitle("#eta");
-  fhDPtVsPtVsEta ->SetYTitle("p_{T} [GeV/c]");
-  fhDPtVsPtVsEta ->SetZTitle("#deltap_{T} [GeV/c]");
-  
-  fhDEtaVsPtVsEta ->SetXTitle("#eta");
-  fhDEtaVsPtVsEta ->SetYTitle("p_{T} [GeV/c]");
-  fhDEtaVsPtVsEta ->SetZTitle("#delta#eta");
-  
-}
-
-//____________________________________________________________________
-void
-AliRecoVsGeneCheck::Event(Double_t* vtx, Double_t* vtx_res, Double_t* mcvtx, Int_t n_part) {
-
-  fhVtzZRecoVsMC     ->Fill(mcvtx[2],vtx[2]);     
-                    
-  fhVtxZRes          ->Fill(vtx_res[2]);
-  fhVtxZResVsZ       ->Fill(mcvtx[2], vtx_res[2]);
-  fhVtxZResVsNPart   ->Fill(n_part, vtx_res[2]);
-                    
-  if (vtx_res[2]!=0) {
-    Float_t dzNorm = (mcvtx[2] - vtx[2])/vtx_res[2];
-
-    fhVtxDzNorm       ->Fill(dzNorm);       
-    fhVtxDzNormVsZ    ->Fill(mcvtx[2], dzNorm);         
-    fhVtxDzNormVsNPart->Fill(n_part, dzNorm);         
-
-  }                 
-  fhVtxZMC      ->Fill(mcvtx[2]);
-  fhVtxZReco    ->Fill(vtx[2]);
-  fhNPart       ->Fill(n_part);     
-}
-
-
-//____________________________________________________________________
-void
-AliRecoVsGeneCheck::Track(AliESDtrack* esdTrack, TParticle* mcParticle) {
-
-  Double_t p[3];
-  esdTrack->GetConstrainedPxPyPz(p); // ### TODO or GetInnerPxPyPy / GetOuterPxPyPy
-  TVector3 vector(p);
-  
-  Float_t reco_p      = esdTrack->GetP();
-  Float_t reco_theta  = vector.Theta();
-  //  Float_t reco_phi    = 
-  Float_t reco_eta    = -TMath::Log(TMath::Tan(reco_theta/2.));
-  Float_t reco_pt     = reco_p*TMath::Sin(reco_theta);
-
-  Float_t mc_eta = mcParticle->Eta();
-  Float_t mc_pt = mcParticle->Pt();
-
-  fhDPtVsPtVsEta->Fill(mc_eta, mc_pt, mc_pt - reco_pt);
-
-  fhDEtaVsPtVsEta->Fill(mc_eta, mc_pt, mc_eta - reco_eta);
-
-}
-
-//____________________________________________________________________
-void 
-AliRecoVsGeneCheck::SaveHistograms(Char_t* dir) {
-
-  gDirectory->mkdir(dir);
-  gDirectory->cd(dir);
-
-  gDirectory->mkdir("event");
-  gDirectory->cd("event");
-
-  fhVtzZRecoVsMC      ->Write();
-  fhVtxZRes           ->Write();
-  fhVtxZResVsZ        ->Write();
-  fhVtxZResVsNPart    ->Write();                         
-  fhVtxDzNorm         ->Write();
-  fhVtxDzNormVsZ      ->Write();
-  fhVtxDzNormVsNPart  ->Write();
-  fhVtxZMC            ->Write();
-  fhVtxZReco          ->Write();
-  fhNPart             ->Write();
-
-  gDirectory->cd("../");
-
-
-  gDirectory->mkdir("track");
-  gDirectory->cd("track");
-  
-  fhDPtVsPtVsEta      ->Write();
-  fhDEtaVsPtVsEta     ->Write();
-  
-  gDirectory->cd("../");
-
-  gDirectory->cd("../");
-}
diff --git a/PWG0/AliRecoVsGeneCheck.h b/PWG0/AliRecoVsGeneCheck.h
deleted file mode 100644 (file)
index 988db56..0000000
+++ /dev/null
@@ -1,70 +0,0 @@
-//
-// Class to check the reconstruction versus the generated
-// 
-// TODO:
-// - add a lot of histograms (much more on tracks vs mc particle) 
-// - add documentation
-//
-
-
-#ifndef ROOT_TObject
-#include "TObject.h"
-#endif
-#ifndef ROOT_TTree
-#include "TTree.h"
-#endif
-#ifndef ROOT_TH3
-#include "TH3.h"
-#endif
-#ifndef ROOT_TParticle
-#include "TParticle.h"
-#endif
-
-#include "AliESDtrack.h"
-
-class AliRecoVsGeneCheck : public TObject
-{
-
-public:
-  AliRecoVsGeneCheck();
-  
-  void   Event(Double_t* vtx, Double_t* vtx_res, Double_t* mcvtx, Int_t n_part=-1);
-
-  void   Track(AliESDtrack* esdTrack, TParticle* mcParticle);
-
-  void   SetNPartAxisTitle(Char_t* t) 
-    {fhVtxZResVsNPart  ->SetXTitle(t); 
-    fhVtxDzNormVsNPart->SetXTitle(t);
-    fhNPart->SetXTitle(t);}
-
-  void   SaveHistograms(Char_t* dir="reco_vs_gene");
-
-protected:
-
-  // event specific histograms
-  TH2F* fhVtzZRecoVsMC;     // z pos of vertex reco vs MC
-  
-  TH1F* fhVtxZRes;          // estimated z res. of vertex
-  TH2F* fhVtxZResVsZ;       // estimated z res. of vertex vs z
-  TH2F* fhVtxZResVsNPart;   // estimated z res. of vertex vs n part
-  
-  TH1F* fhVtxDzNorm;        // (z_mc - z_reco)/estimated res.
-  TH2F* fhVtxDzNormVsZ;     // (z_mc - z_reco)/estimated res. vs vtx z
-  TH2F* fhVtxDzNormVsNPart; // (z_mc - z_reco)/estimated res. vs n part.
-  
-  TH1F* fhVtxZMC;           // vtz z mc 
-  TH1F* fhVtxZReco;         // vtx z reco
-  
-  TH1F* fhNPart;            // n charged primary particles 
-
-  // track specific histograms
-  TH3F* fhDPtVsPtVsEta;
-  TH3F* fhDEtaVsPtVsEta;  
-
-  private:
-    AliRecoVsGeneCheck(const AliRecoVsGeneCheck&);
-    AliRecoVsGeneCheck& operator=(const AliRecoVsGeneCheck&);
-
-  ClassDef(AliRecoVsGeneCheck, 1)
-};
-
diff --git a/PWG0/README b/PWG0/README
deleted file mode 100644 (file)
index bec685b..0000000
+++ /dev/null
@@ -1 +0,0 @@
-empty readme file...
\ No newline at end of file
diff --git a/PWG0/esdTrackCuts/AliTestESDtrackCutsSelector.cxx b/PWG0/esdTrackCuts/AliTestESDtrackCutsSelector.cxx
deleted file mode 100644 (file)
index 3a49e5d..0000000
+++ /dev/null
@@ -1,306 +0,0 @@
-/* $Id$ */
-
-#include "AliTestESDtrackCutsSelector.h"
-
-#include <TStyle.h>
-#include <TSystem.h>
-#include <TCanvas.h>
-#include <TParticle.h>
-#include <TParticlePDG.h>
-#include <TVector3.h>
-#include <TChain.h>
-#include <TFile.h>
-#include <TH1F.h>
-#include <TH3F.h>
-
-#include <TSelector.h>
-#include <TFile.h>
-
-#include <AliLog.h>
-#include <AliESD.h>
-#include <AliRunLoader.h>
-#include <AliStack.h>
-
-#include "AliESDtrackCuts.h"
-#include "AliPWG0Helper.h"
-
-
-ClassImp(AliTestESDtrackCutsSelector)
-
-AliTestESDtrackCutsSelector::AliTestESDtrackCutsSelector() :
-  AliSelectorRL(),
-  fEsdTrackCutsAll(0),
-  fEsdTrackCutsNoVtx(0),
-  fEsdTrackCutsPri(0),
-  fEsdTrackCutsSec(0),
-  fEsdTrackCutsPlusZ(0),
-  fEsdTrackCutsMinusZ(0),
-  fEsdTrackCutsPos(0),
-  fEsdTrackCutsNeg(0),
-  fPIDAfterCutNoVtx(0),
-  fPIDAfterCutAll(0),
-  fVertex(0)
-{
-  //
-  // Constructor. Initialization of pointers
-  //
-}
-
-AliTestESDtrackCutsSelector::~AliTestESDtrackCutsSelector()
-{
-  //
-  // Destructor
-  //
-
-  // histograms are in the output list and deleted when the output
-  // list is deleted by the TSelector dtor
-}
-
-void AliTestESDtrackCutsSelector::Begin(TTree* tree)
-{
-  // Begin function
-
-  AliSelectorRL::Begin(tree);
-
-  ReadUserObjects(tree);
-}
-
-void AliTestESDtrackCutsSelector::ReadUserObjects(TTree* tree)
-{
-  // read the user objects, called from slavebegin and begin
-
-  // only do it once
-  if (fEsdTrackCutsAll)
-    return;
-
-  if (!fEsdTrackCutsAll && fInput)
-    fEsdTrackCutsAll = dynamic_cast<AliESDtrackCuts*> (fInput->FindObject("esdTrackCutsAll")->Clone());
-
-  if (!fEsdTrackCutsAll && tree)
-    fEsdTrackCutsAll = dynamic_cast<AliESDtrackCuts*> (tree->GetUserInfo()->FindObject("esdTrackCutsAll"));
-
-  if (!fEsdTrackCutsAll)
-     AliDebug(AliLog::kError, "ERROR: Could not read fEsdTrackCutsAll from input list.");
-
-  fEsdTrackCutsNoVtx =    dynamic_cast<AliESDtrackCuts*> (fEsdTrackCutsAll->Clone("fEsdTrackCutsNoVtx"));
-  fEsdTrackCutsNoVtx->SetRequireSigmaToVertex(kFALSE);
-
-  fEsdTrackCutsPri =    dynamic_cast<AliESDtrackCuts*> (fEsdTrackCutsAll->Clone("fEsdTrackCutsPri"));
-  fEsdTrackCutsSec =    dynamic_cast<AliESDtrackCuts*> (fEsdTrackCutsAll->Clone("fEsdTrackCutsSec"));
-  fEsdTrackCutsPlusZ =  dynamic_cast<AliESDtrackCuts*> (fEsdTrackCutsAll->Clone("fEsdTrackCutsPlusZ"));
-  fEsdTrackCutsMinusZ = dynamic_cast<AliESDtrackCuts*> (fEsdTrackCutsAll->Clone("fEsdTrackCutsMinusZ"));
-  fEsdTrackCutsPos =    dynamic_cast<AliESDtrackCuts*> (fEsdTrackCutsAll->Clone("fEsdTrackCutsPos"));
-  fEsdTrackCutsNeg =    dynamic_cast<AliESDtrackCuts*> (fEsdTrackCutsAll->Clone("fEsdTrackCutsNeg"));
-}
-
-void AliTestESDtrackCutsSelector::SlaveBegin(TTree* tree)
-{
-  // The SlaveBegin() function is called after the Begin() function.
-  // When running with PROOF SlaveBegin() is called on each slave server.
-  // The tree argument is deprecated (on PROOF 0 is passed).
-
-  AliSelectorRL::SlaveBegin(tree);
-
-  ReadUserObjects(tree);
-
-  fPIDAfterCutNoVtx = new TH1F("fPIDAfterCutNoVtx", "fPIDAfterCutNoVtx", 5001, -2500.5, 2500.5);
-  fPIDAfterCutAll   = new TH1F("fPIDAfterCutAll", "fPIDAfterCutAll", 5001, -2500.5, 2500.5);
-
-  fVertex = new TH3F("fVertex", "fVertex", 100, -10, 10, 100, -10, 10, 100, -10, 10);
-}
-
-void AliTestESDtrackCutsSelector::Init(TTree* tree)
-{
-  // read the user objects
-
-  AliSelectorRL::Init(tree);
-
-  // Enable only the needed branches
-  if (tree)
-  {
-    tree->SetBranchStatus("*", 0);
-    tree->SetBranchStatus("fTriggerMask", 1);
-    tree->SetBranchStatus("fSPDVertex*", 1);
-    tree->SetBranchStatus("fTracks.fLabel", 1);
-
-    AliESDtrackCuts::EnableNeededBranches(tree);
-  }
-}
-
-Bool_t AliTestESDtrackCutsSelector::Process(Long64_t entry)
-{
-  // The Process() function is called for each entry in the tree (or possibly
-  // keyed object in the case of PROOF) to be processed. The entry argument
-  // specifies which entry in the currently loaded tree is to be processed.
-  // It can be passed to either TTree::GetEntry() or TBranch::GetEntry()
-  // to read either all or the required parts of the data. When processing
-  // keyed objects with PROOF, the object is already loaded and is available
-  // via the fObject pointer.
-  //
-  // This function should contain the "body" of the analysis. It can contain
-  // simple or elaborate selection criteria, run algorithms on the data
-  // of the event and typically fill histograms.
-
-  // WARNING when a selector is used with a TChain, you must use
-  //  the pointer to the current TTree to call GetEntry(entry).
-  //  The entry is always the local entry number in the current tree.
-  //  Assuming that fTree is the pointer to the TChain being processed,
-  //  use fTree->GetTree()->GetEntry(entry).
-
-  if (AliSelectorRL::Process(entry) == kFALSE)
-    return kFALSE;
-
-  // Check prerequisites
-  if (!fESD) {
-    AliDebug(AliLog::kError, "ESD branch not available");
-    return kFALSE;
-  }
-
-  if (!AliPWG0Helper::IsVertexReconstructed(fESD)) {
-    AliDebug(AliLog::kDebug+5, "Vertex is not reconstructed");
-    return kFALSE;
-  }
-
-  // check if the esd track cut objects are there
-  if (!fEsdTrackCutsAll || !fEsdTrackCutsPri || !fEsdTrackCutsSec || !fEsdTrackCutsPlusZ || !fEsdTrackCutsMinusZ || !fEsdTrackCutsPos || !fEsdTrackCutsNeg) {
-    AliDebug(AliLog::kError, "fEsdTrackCutsXXX not available");
-    return kFALSE;
-  }
-
-  // get particle stack
-  AliStack* stack = GetStack();
-  if (!stack) {
-    AliDebug(AliLog::kError, "Stack not available");
-    return kFALSE;
-  }
-  Int_t nPrim  = stack->GetNprimary();
-  
-  // ########################################################
-  // loop over esd tracks
-  Int_t nTracks = fESD->GetNumberOfTracks();
-
-  // count the number of "good" tracks as parameter for vertex reconstruction efficiency
-  for (Int_t t=0; t<nTracks; t++) {
-    AliDebug(AliLog::kDebug+1, Form("ESD Loop: Processing track %d.", t));
-
-    AliESDtrack* esdTrack = fESD->GetTrack(t);
-
-    Bool_t passed = fEsdTrackCutsAll->AcceptTrack(esdTrack);
-
-    // using the properties of the mc particle
-    Int_t label = TMath::Abs(esdTrack->GetLabel());
-    if (label == 0) {
-      AliDebug(AliLog::kWarning, Form("WARNING: cannot find corresponding mc part for track %d.", t));
-      continue;
-    }
-    TParticle* particle = stack->Particle(label);
-    if (!particle) {
-      AliDebug(AliLog::kError, Form("UNEXPECTED: part with label %d not found in stack (track loop).", label));
-      continue;
-    }
-
-    if (label < nPrim)
-      fEsdTrackCutsPri->AcceptTrack(esdTrack);
-    else
-    {
-      fEsdTrackCutsSec->AcceptTrack(esdTrack);
-      if (passed)
-      {
-        fPIDAfterCutAll->Fill(particle->GetPdgCode());
-        fVertex->Fill(particle->Vx(), particle->Vy(), particle->Vz());
-      }
-
-      if (fEsdTrackCutsNoVtx->AcceptTrack(esdTrack))
-        fPIDAfterCutNoVtx->Fill(particle->GetPdgCode());
-    }
-
-    TParticlePDG* pdgPart = particle->GetPDG();
-    if (pdgPart)
-    {
-      if (pdgPart->Charge() > 0)
-        fEsdTrackCutsPos->AcceptTrack(esdTrack);
-      else if (pdgPart->Charge() < 0)
-        fEsdTrackCutsNeg->AcceptTrack(esdTrack);
-    }
-
-    if (particle->Eta() < 0)
-      fEsdTrackCutsPlusZ->AcceptTrack(esdTrack);
-    else
-      fEsdTrackCutsMinusZ->AcceptTrack(esdTrack);
-  }
-  
-  return kTRUE;
-}
-
-void AliTestESDtrackCutsSelector::SlaveTerminate()
-{
-  // The SlaveTerminate() function is called after all entries or objects
-  // have been processed. When running with PROOF SlaveTerminate() is called
-  // on each slave server.
-
-  AliSelectorRL::SlaveTerminate();
-  
-  // Add the histograms to the output on each slave server
-  if (!fOutput)
-  {
-    AliDebug(AliLog::kError, Form("ERROR: Output list not initialized."));
-    return;
-  }
-
-  fOutput->Add(fEsdTrackCutsAll);
-  fOutput->Add(fEsdTrackCutsNoVtx);
-  fOutput->Add(fEsdTrackCutsPri);
-  fOutput->Add(fEsdTrackCutsSec);
-  fOutput->Add(fEsdTrackCutsPlusZ);
-  fOutput->Add(fEsdTrackCutsMinusZ);
-  fOutput->Add(fEsdTrackCutsPos);
-  fOutput->Add(fEsdTrackCutsNeg);
-  fOutput->Add(fPIDAfterCutNoVtx);
-  fOutput->Add(fPIDAfterCutAll);
-  fOutput->Add(fVertex);
-}
-
-void AliTestESDtrackCutsSelector::Terminate()
-{
-  // The Terminate() function is the last function to be called during
-  // a query. It always runs on the client, it can be used to present
-  // the results graphically or save the results to file.
-
-  AliSelectorRL::Terminate();
-
-  fEsdTrackCutsAll = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("esdTrackCutsAll"));
-  fEsdTrackCutsNoVtx = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("fEsdTrackCutsNoVtx"));
-  fEsdTrackCutsPri = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("fEsdTrackCutsPri"));
-  fEsdTrackCutsSec = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("fEsdTrackCutsSec"));
-  fEsdTrackCutsPlusZ = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("fEsdTrackCutsPlusZ"));
-  fEsdTrackCutsMinusZ = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("fEsdTrackCutsMinusZ"));
-  fEsdTrackCutsPos = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("fEsdTrackCutsPos"));
-  fEsdTrackCutsNeg = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("fEsdTrackCutsNeg"));
-  fPIDAfterCutNoVtx = dynamic_cast<TH1F*> (fOutput->FindObject("fPIDAfterCutNoVtx"));
-  fPIDAfterCutAll = dynamic_cast<TH1F*> (fOutput->FindObject("fPIDAfterCutAll"));
-  fVertex = dynamic_cast<TH3F*> (fOutput->FindObject("fVertex"));
-
-  // check if the esd track cut objects are there
-  if (!fEsdTrackCutsAll || !fEsdTrackCutsPri || !fEsdTrackCutsSec || !fEsdTrackCutsPlusZ || !fEsdTrackCutsMinusZ || !fEsdTrackCutsPos || !fEsdTrackCutsNeg) {
-    AliDebug(AliLog::kError, Form("fEsdTrackCutsXXX not available %p %p %p %p %p %p %p", fEsdTrackCutsAll, fEsdTrackCutsPri, fEsdTrackCutsSec, fEsdTrackCutsPlusZ, fEsdTrackCutsMinusZ, fEsdTrackCutsPos, fEsdTrackCutsNeg));
-    return;
-  }
-
-  TFile* file = TFile::Open("trackCuts.root", "RECREATE");
-
-  fEsdTrackCutsAll->SaveHistograms();
-  fEsdTrackCutsNoVtx->SaveHistograms();
-  fEsdTrackCutsPri->SaveHistograms();
-  fEsdTrackCutsSec->SaveHistograms();
-  fEsdTrackCutsPlusZ->SaveHistograms();
-  fEsdTrackCutsMinusZ->SaveHistograms();
-  fEsdTrackCutsPos->SaveHistograms();
-  fEsdTrackCutsNeg->SaveHistograms();
-  fPIDAfterCutNoVtx->Write();
-  fPIDAfterCutAll->Write();
-  fVertex->Write();
-
-  file->Close();
-
-       fEsdTrackCutsAll->DrawHistograms();
-}
diff --git a/PWG0/esdTrackCuts/AliTestESDtrackCutsSelector.h b/PWG0/esdTrackCuts/AliTestESDtrackCutsSelector.h
deleted file mode 100644 (file)
index 53b93b0..0000000
+++ /dev/null
@@ -1,50 +0,0 @@
-/* $Id$ */
-
-#ifndef ALITESTESDTRACKCUTSSELECTOR_H
-#define ALITESTESDTRACKCUTSSELECTOR_H
-
-#include "AliSelectorRL.h"
-
-class AliESDtrackCuts;
-class TH1F;
-class TH3F;
-
-class AliTestESDtrackCutsSelector : public AliSelectorRL {
-  public:
-    AliTestESDtrackCutsSelector();
-    virtual ~AliTestESDtrackCutsSelector();
-
-    virtual void    Begin(TTree* tree);
-    virtual void    SlaveBegin(TTree *tree);
-    virtual void    Init(TTree *tree);
-    virtual Bool_t  Process(Long64_t entry);
-    virtual void    SlaveTerminate();
-    virtual void    Terminate();
-
- protected:
-    void ReadUserObjects(TTree* tree);
-
-    AliESDtrackCuts*  fEsdTrackCutsAll;  // esd track cuts for all tracks   
-    AliESDtrackCuts*  fEsdTrackCutsNoVtx;  // all cuts except vtx
-
-    AliESDtrackCuts*  fEsdTrackCutsPri;  // cuts for tracks from primary particles
-    AliESDtrackCuts*  fEsdTrackCutsSec;  // cuts for tracks from secondary particles
-    AliESDtrackCuts*  fEsdTrackCutsPlusZ;  // cuts for tracks that go to z > 0
-    AliESDtrackCuts*  fEsdTrackCutsMinusZ;  // cuts for tracks that go to z < 0
-    AliESDtrackCuts*  fEsdTrackCutsPos;  // cuts for tracks from positive particles
-    AliESDtrackCuts*  fEsdTrackCutsNeg;  // cuts for tracks from negative particles
-
-    TH1F*             fPIDAfterCutNoVtx;      // true PID of tracks that passed all cuts except vtx
-    TH1F*             fPIDAfterCutAll;        // true PID of tracks that passed all cuts incl. vtx
-
-    TH3F*             fVertex;                // originating vertex of specific particles
-
-
- private:
-    AliTestESDtrackCutsSelector(const AliTestESDtrackCutsSelector&);
-    AliTestESDtrackCutsSelector& operator=(const AliTestESDtrackCutsSelector&);
-
-  ClassDef(AliTestESDtrackCutsSelector, 0);
-};
-
-#endif
diff --git a/PWG0/esdTrackCuts/testESDtrackCuts.C b/PWG0/esdTrackCuts/testESDtrackCuts.C
deleted file mode 100644 (file)
index 9ee1fa1..0000000
+++ /dev/null
@@ -1,75 +0,0 @@
-/* $Id$ */
-
-//
-// script to run the AliMultiplicityESDSelector
-//
-
-#include "../CreateESDChain.C"
-#include "../PWG0Helper.C"
-
-TChain* testESDtrackCuts(Char_t* data, Int_t nRuns=20, Int_t offset=0, Bool_t aDebug = kFALSE, Bool_t aProof = kFALSE, Char_t* proofServer = "lxb6046")
-{
-  if (aProof)
-    connectProof(proofServer);
-
-  TString libraries("libEG;libGeom;libESD;libPWG0base");
-  TString packages("PWG0base");
-
-  if (!prepareQuery(libraries, packages, 1))
-    return;
-
-  // selection of esd tracks
-  AliESDtrackCuts* esdTrackCutsAll = new AliESDtrackCuts("esdTrackCutsAll");
-
-  esdTrackCutsAll->DefineHistograms(1);
-  esdTrackCutsAll->SetMinNClustersTPC(50);
-  esdTrackCutsAll->SetMaxChi2PerClusterTPC(3.5);
-  esdTrackCutsAll->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
-  esdTrackCutsAll->SetRequireTPCRefit(kTRUE);
-  esdTrackCutsAll->SetMinNsigmaToVertex(3);
-  esdTrackCutsAll->SetRequireSigmaToVertex(kTRUE);
-  esdTrackCutsAll->SetAcceptKingDaughters(kFALSE);
-
-  TList inputList;
-  inputList.Add(esdTrackCutsAll);
-
-  TChain* chain = CreateESDChain(data, nRuns, offset);
-
-  TString selectorName = "AliTestESDtrackCutsSelector";
-  AliLog::SetClassDebugLevel(selectorName, AliLog::kInfo);
-
-  selectorName += ".cxx++";
-
-  if (aDebug != kFALSE)
-    selectorName += "g";
-
-  Int_t result = executeQuery(chain, &inputList, selectorName);
-
-  if (result != 0)
-  {
-    printf("ERROR: Executing process failed with %d.\n", result);
-    return;
-  }
-}
-
-void draw(const char* dir, const char* fileName = "trackCuts.root")
-{
-  /*
-   draw("esdTrackCutsAll")
-   draw("fEsdTrackCutsPri")
-   draw("fEsdTrackCutsSec")
-   draw("fEsdTrackCutsPlusZ")
-   draw("fEsdTrackCutsMinusZ")
-   draw("fEsdTrackCutsPos")
-   draw("fEsdTrackCutsNeg")
-  */
-
-  gSystem->Load("libPWG0base");
-
-  TFile::Open(fileName);
-
-  AliESDtrackCuts* cuts = new AliESDtrackCuts(dir, dir);
-  cuts->LoadHistograms();
-
-  cuts->DrawHistograms();
-}