From f9374725164d8dbe9adea975f1c6133ab9f8f0ce Mon Sep 17 00:00:00 2001 From: morsch Date: Thu, 30 Jul 2009 09:12:51 +0000 Subject: [PATCH] Obsolete code removed. --- JETAN/AliJetESDmcReader.cxx | 175 ------------------------------------ JETAN/AliJetESDmcReader.h | 36 -------- JETAN/JETANMCLinkDef.h | 10 --- JETAN/libJETANMC.pkg | 12 --- 4 files changed, 233 deletions(-) delete mode 100644 JETAN/AliJetESDmcReader.cxx delete mode 100644 JETAN/AliJetESDmcReader.h delete mode 100644 JETAN/JETANMCLinkDef.h delete mode 100644 JETAN/libJETANMC.pkg diff --git a/JETAN/AliJetESDmcReader.cxx b/JETAN/AliJetESDmcReader.cxx deleted file mode 100644 index 258088c47d7..00000000000 --- a/JETAN/AliJetESDmcReader.cxx +++ /dev/null @@ -1,175 +0,0 @@ -/************************************************************************** - * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * - * * - * Author: The ALICE Off-line Project. * - * Contributors are mentioned in the code where appropriate. * - * * - * Permission to use, copy, modify and distribute this software and its * - * documentation strictly for non-commercial purposes is hereby granted * - * without fee, provided that the above copyright notice appears in all * - * copies and that both the copyright notice and this permission notice * - * appear in the supporting documentation. The authors make no claims * - * about the suitability of this software for any purpose. It is * - * provided "as is" without express or implied warranty. * - **************************************************************************/ - -// Jet ESD Reader -// ESD reader for jet analysis (it reads the esd and the MC trees) -// Author: Mercedes Lopez Noriega (mercedes.lopez.noriega@cern.ch) - -#include -#include -#include -#include -#include -#include -#include -#include -#include "AliJetESDmcReader.h" -#include "AliJetESDReaderHeader.h" -#include "AliESDEvent.h" -#include "AliESDtrack.h" -#include "AliHeader.h" - -ClassImp(AliJetESDmcReader) - -AliJetESDmcReader::AliJetESDmcReader(): - AliJetESDReader(), - fChainMC(0x0), - fAliHeader(0) -{ - // Constructor -} - -//____________________________________________________________________________ - -AliJetESDmcReader::~AliJetESDmcReader() -{ - // Destructor - delete fChainMC; -} - -//____________________________________________________________________________ - -void AliJetESDmcReader::OpenInputFiles() - -{ - // chain for the ESDs - fChain = new TChain("esdTree"); - fChainMC = new TChain("mcStackTree"); - - // get directory and pattern name from the header - const char* dirName=fReaderHeader->GetDirectory(); - const char* pattern=fReaderHeader->GetPattern(); - - // Add files matching patters to the chain - void *dir = gSystem->OpenDirectory(dirName); - const char *name = 0x0; - int nesd = ((AliJetESDReaderHeader*)fReaderHeader)->GetNesd(); - int a = 0; - while ((name = gSystem->GetDirEntry(dir))){ - if (a>=nesd) continue; - if (strstr(name,pattern)){ - printf("Adding %s\n",name); - char path[256]; - sprintf(path,"%s/%s",dirName,name); - fChain->AddFile(path,-1); - fChainMC->AddFile(path,-1); - a++; - } - } - printf("%d ESDs added\n",a); - - gSystem->FreeDirectory(dir); - fChain->SetBranchAddress("ESD", &fESD); - fChainMC->SetBranchAddress("Header", &fAliHeader); - fChainMC->SetBranchAddress("Stack", &fArrayMC); - - int nMax = fChain->GetEntries(); - printf("\nTotal number of events in chain= %d\n",nMax); - - // set number of events in header - if (fReaderHeader->GetLastEvent() == -1) - fReaderHeader->SetLastEvent(nMax); - else { - Int_t nUsr = fReaderHeader->GetLastEvent(); - fReaderHeader->SetLastEvent(TMath::Min(nMax,nUsr)); - } -} - -//____________________________________________________________________________ - -Bool_t AliJetESDmcReader::FillMomentumArray(Int_t event) -{ - // Fill momentum array - - Int_t goodTrack = 0; - Int_t nt = 0; - Int_t pdgCode = 0; - Float_t pt, eta; - Float_t energyMC, pxMC, pyMC, pzMC, ptMC; // Monte Carlo - TVector3 p3; - - // clear momentum array - ClearArray(); - - // get event from chain - fChain->GetEntry(event); - fChainMC->GetEntry(event); - - // get number of tracks in event (for the loop) - nt = fESD->GetNumberOfTracks(); - - // temporary storage of signal and pt cut flag - Int_t* sflag = new Int_t[nt]; - Int_t* cflag = new Int_t[nt]; - - // get cuts set by user - Float_t ptMin = fReaderHeader->GetPtCut(); - Float_t etaMin = fReaderHeader->GetFiducialEtaMin(); - Float_t etaMax = fReaderHeader->GetFiducialEtaMax(); - - //loop over tracks - for (Int_t it = 0; it < nt; it++) { - AliESDtrack *track = fESD->GetTrack(it); - UInt_t status = track->GetStatus(); - - Double_t mom[3]; - track->GetPxPyPz(mom); - p3.SetXYZ(mom[0],mom[1],mom[2]); - pt = p3.Pt(); - if (((status & AliESDtrack::kITSrefit) == 0) || - ((status & AliESDtrack::kTPCrefit) == 0)) continue; // quality check - if (((AliJetESDReaderHeader*) fReaderHeader)->ReadSignalOnly() - && TMath::Abs(track->GetLabel()) > 10000) continue; // quality check - if (((AliJetESDReaderHeader*) fReaderHeader)->ReadBkgdOnly() - && TMath::Abs(track->GetLabel()) < 10000) continue; // quality check - eta = p3.Eta(); - if ( (eta > etaMax) || (eta < etaMin)) continue; // checking eta cut - - new ((*fMomentumArray)[goodTrack]) TLorentzVector(p3,p3.Mag()); - - sflag[goodTrack]=0; - if (TMath::Abs(track->GetLabel()) < 10000) sflag[goodTrack]=1; - cflag[goodTrack]=0; - if (pt > ptMin) cflag[goodTrack]=1; // pt cut - - goodTrack++; - - // Monte Carlo information - Int_t label = TMath::Abs(track->GetLabel()); - TClonesArray &arrayMC = *fArrayMC; - TParticle *part = (TParticle*)arrayMC[label]; //particle - ptMC = part->Pt(); - pdgCode = part->GetPdgCode(); - energyMC = part->Energy(); - pxMC = part->Px(); pyMC = part->Py(); pzMC = part->Pz(); - } - // set the signal flags - fSignalFlag.Set(goodTrack,sflag); - fCutFlag.Set(goodTrack,cflag); - return kTRUE; - -} - - diff --git a/JETAN/AliJetESDmcReader.h b/JETAN/AliJetESDmcReader.h deleted file mode 100644 index 31969f8e2c6..00000000000 --- a/JETAN/AliJetESDmcReader.h +++ /dev/null @@ -1,36 +0,0 @@ -#ifndef ALIJETESDMCREADER_H -#define ALIJETESDMCREADER_H - -/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * - * See cxx source for full Copyright notice */ - -// Jet ESD Reader -// ESD reader for jet analysis (it reads the esd and the MC trees) -// Author: Mercedes Lopez Noriega (mercedes.lopez.noriega@cern.ch) - -#include "AliJetESDReader.h" -class AliJetESDReaderHeader; -class AliHeader; - - - -class AliJetESDmcReader : public AliJetESDReader -{ - public: - AliJetESDmcReader(); - virtual ~AliJetESDmcReader(); - // Setters - Bool_t FillMomentumArray(Int_t event); - Bool_t FillMomentumArray() {return kTRUE;} - void OpenInputFiles(); - private: - AliJetESDmcReader(const AliJetESDmcReader &det); - AliJetESDmcReader &operator=(const AliJetESDmcReader &det); - - protected: - TChain* fChainMC; // chain for mc information - AliHeader* fAliHeader; //! Event header - ClassDef(AliJetESDmcReader,1) -}; - -#endif diff --git a/JETAN/JETANMCLinkDef.h b/JETAN/JETANMCLinkDef.h deleted file mode 100644 index 3d641c8a32d..00000000000 --- a/JETAN/JETANMCLinkDef.h +++ /dev/null @@ -1,10 +0,0 @@ -#ifdef __CINT__ - -#pragma link off all globals; -#pragma link off all classes; -#pragma link off all functions; - - -#pragma link C++ class AliJetESDmcReader; - -#endif diff --git a/JETAN/libJETANMC.pkg b/JETAN/libJETANMC.pkg deleted file mode 100644 index 763dc178f74..00000000000 --- a/JETAN/libJETANMC.pkg +++ /dev/null @@ -1,12 +0,0 @@ -#-*-Mode: Makefile-*- - -SRCS = AliJetESDmcReader.cxx - -HDRS:= $(SRCS:.cxx=.h) -EINCLUDE:= PYTHIA6 -DHDR= JETANMCLinkDef.h - -ifeq (win32gcc,$(ALICE_TARGET)) -PACKSOFLAGS:= $(SOFLAGS) -L$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET) \ - -lJETAN -lESD -endif -- 2.43.0