removal of obsolete classes - cleanup of AliITSClusterFinder.cxx
authormasera <masera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 14 Mar 2007 10:27:24 +0000 (10:27 +0000)
committermasera <masera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 14 Mar 2007 10:27:24 +0000 (10:27 +0000)
ITS/AliITSClusterFinder.cxx
ITS/AliITSFindClustersV2.cxx [deleted file]
ITS/AliITSFindClustersV2.h [deleted file]
ITS/ITSrecLinkDef.h
ITS/libITSrec.pkg

index e4306ea..3ec86d3 100644 (file)
@@ -400,38 +400,6 @@ void AliITSClusterFinder::Print(ostream *os) const{
     *os << fClusterSize<<",";
     *os << fNPeaks<<endl;
 }
-/*
-//______________________________________________________________________
-void AliITSClusterFinder::RecPoints2Clusters
-(const TClonesArray *points, Int_t idx, TClonesArray *clusters) {
-  //------------------------------------------------------------
-  // Conversion AliITSRecPoint -> AliITSclusterV2 for the ITS 
-  // subdetector indexed by idx 
-  //------------------------------------------------------------
-  if(!fDetTypeRec->GetITSgeom()) {
-    Error("RecPoints2Clusters","ITS geom is null!");
-    return;
-  }
-  Int_t lastSPD1=fDetTypeRec->GetITSgeom()->GetModuleIndex(2,1,1)-1;
-  TClonesArray &cl=*clusters;
-  Int_t ncl=points->GetEntriesFast();
-  for (Int_t i=0; i<ncl; i++) {
-    AliITSRecPoint *p = (AliITSRecPoint *)points->UncheckedAt(i);
-    Float_t lp[5];
-    lp[0]=-(-p->GetX()+fYshift[idx]); if (idx<=lastSPD1) lp[0]*=-1; //SPD1
-    lp[1]=  -p->GetZ()+fZshift[idx];
-    lp[2]=p->GetSigmaX2();
-    lp[3]=p->GetSigmaZ2();
-    lp[4]=p->GetQ()*36./23333.;  //electrons -> ADC
-    Int_t lab[4]; 
-    lab[0]=p->GetLabel(0); lab[1]=p->GetLabel(1); lab[2]=p->GetLabel(2);
-    lab[3]=fNdet[idx];
-    CheckLabels(lab);
-    Int_t dummy[3]={0,0,0};
-    new (cl[i]) AliITSclusterV2(lab,lp, dummy);
-  }  
-} 
-*/
 //______________________________________________________________________
 void AliITSClusterFinder::Read(istream *is)  {
     //Standard input for this class
@@ -479,29 +447,3 @@ istream &operator>>(istream &is,AliITSClusterFinder &source){
     source.Read(&is);
     return is;
 }
-/*
-void AliITSClusterFinder::RecPoints2Clusters
-(const TClonesArray *points, Int_t idx, TClonesArray *clusters) {
-  //------------------------------------------------------------
-  // Conversion AliITSRecPoint -> AliITSclusterV2 for the ITS 
-  // subdetector indexed by idx 
-  //------------------------------------------------------------
-  TClonesArray &cl=*clusters;
-  Int_t ncl=points->GetEntriesFast();
-  for (Int_t i=0; i<ncl; i++) {
-    AliITSRecPoint *p = (AliITSRecPoint *)points->UncheckedAt(i);
-    Float_t lp[5];
-    lp[0]=-(-p->GetX()+fYshift[idx]); if (idx<=fLastSPD1) lp[0]*=-1; //SPD1
-    lp[1]=  -p->GetZ()+fZshift[idx];
-    lp[2]=p->GetSigmaX2();
-    lp[3]=p->GetSigmaZ2();
-    lp[4]=p->GetQ()*36./23333.;  //electrons -> ADC
-    Int_t lab[4]; 
-    lab[0]=p->GetLabel(0); lab[1]=p->GetLabel(1); lab[2]=p->GetLabel(2);
-    lab[3]=fNdet[idx];
-    CheckLabels(lab);
-    Int_t dummy[3]={0,0,0};
-    new (cl[i]) AliITSclusterV2(lab,lp, dummy);
-  }  
-} 
-*/
diff --git a/ITS/AliITSFindClustersV2.cxx b/ITS/AliITSFindClustersV2.cxx
deleted file mode 100644 (file)
index 7a5904c..0000000
+++ /dev/null
@@ -1,413 +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.                  *
- **************************************************************************/
-/* $Id$ */
-///////////////////////////////////////////////////////////////////
-//Class for reconstruction of clusters V2                        //
-//                                                               //
-//                                                               //
-///////////////////////////////////////////////////////////////////
-
-#include <TFile.h>
-#include <TTree.h>
-#include <TBranch.h>
-#include <TClonesArray.h>
-#include <TParticle.h>
-
-#include "AliRun.h"
-#include "AliHeader.h"
-
-#include "AliITSRecPoint.h"
-#include "AliITSFindClustersV2.h"
-#include "AliITSclusterV2.h"
-#include "AliITSgeom.h"
-#include "AliMC.h"
-
-ClassImp(AliITSFindClustersV2)
-
-//______________________________________________________________________
-AliITSFindClustersV2::AliITSFindClustersV2():
-fAr(0),
-fDeletfAr(kFALSE),
-fGeom(0),
-fInFileName(0),
-fOutFileName(0),
-fIn(0),
-fOut(0),
-fInit(kFALSE),
-fSlowFast(kFALSE){
-    // Default constructor.
-    // Inputs:
-    //  none.
-    // Outputs:
-    //   none.
-    // Return:
-    //    A zero-ed constructed AliITSFindClustersV2 class.
-
-}
-//______________________________________________________________________
-AliITSFindClustersV2::AliITSFindClustersV2(const TString infile,
-                                           const TString outfile):
-fAr(0),
-fDeletfAr(kFALSE),
-fGeom(0),
-fInFileName(0),
-fOutFileName(0),
-fIn(0),
-fOut(0),
-fInit(kFALSE),
-fSlowFast(kFALSE){
-    // Standard constructor.
-    // Inputs:
-    //  const TString infile   Input file name where the RecPoints are
-    //                         to be read from.
-    //  const TString outfile  Output file where V2 tracking clulsters
-    //                         are to be written. if =="" writen to the
-    //                         same file as input.
-    // Outputs:
-    //   none.
-    // Return:
-    //    A standardly constructed AliITSFindClustersV2 class.
-
-    fInFileName = new TString(infile);
-    if(outfile.CompareTo("")!=0){
-       fOutFileName = new TString(outfile);
-    } // end if outfile.CompareeTo("")!=0
-    
-    if(fOutFileName!=0){
-       fIn = new TFile(fInFileName->Data(),"READ");
-       fOut = new TFile(fOutFileName->Data(),"UPDATE");
-    }else{ // open fIn file for update
-       fIn = new TFile(fInFileName->Data(),"UPDATE");
-    } // end if
-
-    fAr  = (AliRun*) fIn->Get("gAlice");
-    if(!fAr){
-       Warning("AliITSFindClusterV2",
-               "Can't fine gAlice in file %s. Aborting.",fIn->GetName());
-       return;
-    } // end if !fAr
-    fDeletfAr = kTRUE; // Since gAlice was read in, delete it.
-    /*
-    AliITS *its = (AliITS*) fAr->GetModule("ITS");
-    if(!its){
-       Warning("AliITSFindClusterV2",
-               "Can't fine the ITS in gAlice. Aborting.");
-       return;
-    } // end if !its
-    fGeom = its->GetITSgeom();
-    if(!fGeom){
-       Warning("AliITSFindClusterV2",
-               "Can't fine the ITS geometry in gAlice. Aborting.");
-       return;
-    } // end if !fGeom
-    */
-    AliRunLoader* rl = AliRunLoader::Open("galice.root");
-    rl->CdGAFile();
-    fGeom = (AliITSgeom*)gDirectory->Get("AliITSgeom");
-    if(!fGeom){
-       Warning("AliITSFindClusterV2",
-               "Can't fine the ITS geometry in gAlice. Aborting.");
-       return;
-    } // end if !fGeom
-    delete rl;
-
-    if(fOut) fOut->cd();
-    fInit = kTRUE;
-}
-//______________________________________________________________________
-AliITSFindClustersV2::AliITSFindClustersV2(TFile *in,
-                                          TFile *out):
-fAr(0),
-fDeletfAr(kFALSE),
-fGeom(0),
-fInFileName(0),
-fOutFileName(0),
-fIn(0),
-fOut(0),
-fInit(kFALSE),
-fSlowFast(kFALSE){
-    // Standard constructor.
-    // Inputs:
-    //  const TFile *in   Input file where the RecPoints are
-    //                         to be read from.
-    //  const Tfile *out  Output file where V2 tracking clulsters
-    //                         are to be written. if ==0 writen to the
-    //                         same file as input.
-    // Outputs:
-    //   none.
-    // Return:
-    //    A standardly constructed AliITSFindClustersV2 class.
-
-    fIn  = in;
-    fOut = out;
-    fAr  = (AliRun*) fIn->Get("gAlice");
-    if(!fAr){
-       Warning("AliITSFindClusterV2",
-               "Can't fine gAlice in file %s. Aborting.",fIn->GetName());
-       return;
-    } // end if !fAr
-    fDeletfAr = kTRUE; // Since gAlice was read in, delete it.
-    /*
-    AliITS *its = (AliITS*) fAr->GetModule("ITS");
-    if(!its){
-       Warning("AliITSFindClusterV2",
-               "Can't fine the ITS in gAlice. Aborting.");
-       return;
-    } // end if !its
-    fGeom = its->GetITSgeom();
-    */
-    AliRunLoader* rl = AliRunLoader::Open("galice.root");
-    rl->CdGAFile();
-    fGeom = (AliITSgeom*)gDirectory->Get("AliITSgeom");
-    if(!fGeom){
-       Warning("AliITSFindClusterV2",
-               "Can't fine the ITS geometry in gAlice. Aborting.");
-       return;
-    } // end if !fGeom
-    delete rl;
-    if(fOut) fOut->cd();
-    fInit = kTRUE;
-}
-//______________________________________________________________________
-AliITSFindClustersV2::AliITSFindClustersV2(AliRun *ar,
-                                          const TString outfile):
-fAr(0),
-fDeletfAr(kFALSE),
-fGeom(0),
-fInFileName(0),
-fOutFileName(0),
-fIn(0),
-fOut(0),
-fInit(kFALSE),
-fSlowFast(kFALSE){
-    // Standard constructor.
-    // Inputs:
-    //  const TString infile   Input file name where the RecPoints are
-    //                         to be read from.
-    //  const TString outfile  Output file where V2 tracking clulsters
-    //                         are to be written. if =="" writen to the
-    //                         same file as input.
-    // Outputs:
-    //   none.
-    // Return:
-    //    A standardly constructed AliITSFindClustersV2 class.
-
-    if(outfile.CompareTo("")!=0){
-       fOutFileName = new TString(outfile);
-    } // end if outfile.CompareeTo("")!=0
-    
-    if(fOutFileName!=0){
-       fOut = new TFile(fOutFileName->Data(),"UPDATE");
-    } // end if
-
-    fAr  = ar;
-    if(!fAr){
-       Warning("AliITSFindClusterV2",
-               "ar==0. Aborting.");
-       return;
-    } // end if !fAr
-
-    /*   AliITS *its = (AliITS*) fAr->GetModule("ITS");
-    if(!its){
-       Warning("AliITSFindClusterV2",
-               "Can't fine the ITS in gAlice. Aborting.");
-       return;
-    } // end if !its
-    fGeom = its->GetITSgeom();
-    */
-    AliRunLoader* rl = AliRunLoader::Open("galice.root");
-    rl->CdGAFile();
-    fGeom = (AliITSgeom*)gDirectory->Get("AliITSgeom");
-    if(!fGeom){
-       Warning("AliITSFindClusterV2",
-               "Can't fine the ITS geometry in gAlice. Aborting.");
-       return;
-    } // end if !fGeom
-    delete rl;
-
-    if(fOut) fOut->cd();
-    fInit = kTRUE;
-}
-//______________________________________________________________________
-AliITSFindClustersV2::AliITSFindClustersV2(const AliITSFindClustersV2 &rec):TTask(rec),
-fAr(rec.fAr),
-fDeletfAr(rec.fDeletfAr),
-fGeom(rec.fGeom),
-fInFileName(rec.fInFileName),
-fOutFileName(rec.fOutFileName),
-fIn(rec.fIn),
-fOut(rec.fOut),
-fInit(rec.fInit),
-fSlowFast(rec.fSlowFast){
-    // Copy constructor. 
-
-  
-}
-//______________________________________________________________________
-AliITSFindClustersV2& AliITSFindClustersV2::operator=(const AliITSFindClustersV2& source){
-    // Assignment operator. 
-  this->~AliITSFindClustersV2();
-  new(this) AliITSFindClustersV2(source);
-  return *this;
-}
-
-//______________________________________________________________________
-AliITSFindClustersV2::~AliITSFindClustersV2(){
-    // Default constructor.
-    // Inputs:
-    //  none.
-    // Outputs:
-    //   none.
-    // Return:
-    //    A destroyed AliITSFindclustersV2 class.
-
-    fGeom = 0; // Deleted when AliRun/ITS is deleted.
-    if(fInFileName!=0){ // input file opened by AliITSFindClustersV2
-       if(fIn!=0){
-           if(fIn->IsOpen()) fIn->Close();
-           delete fIn;
-           fIn = 0;
-       } // end if
-       delete fInFileName;
-       fInFileName = 0;
-    } // end if
-
-    if(fOutFileName!=0){ // input file opened by AliITSFindClustersV2
-       if(fOut!=0){
-           if(fOut->IsOpen()) fOut->Close();
-           delete fOut;
-           fOut = 0;
-       } // end if
-       delete fOutFileName;
-       fOutFileName = 0;
-    } // end if
-    if(fDeletfAr && !fAr){
-       cout << "deleting fAr."<< endl;
-       delete fAr;
-       fAr = 0;
-       cout << "fAr deleted OK."<< endl;
-    } // end if fDeletfAr
-}
-//______________________________________________________________________ 
-void AliITSFindClustersV2::Exec(const Option_t *opt){
-    // Main FindclustersV2 function.
-    // Inputs:
-    //      Option_t * opt   list of subdetector to digitize. =0 all.
-    // Outputs:
-    //      none.
-    // Return:
-    //      none.
-    Char_t name[50];
-
-    if(!fInit){
-       Warning("Exec","(opt=%s) Initilization not succesfull. Aborting.\n",opt);
-       return;
-    } // end if !fInit
-
-    fGeom->Write();
-
-    fAr->GetEvent(0);
-    TTree *pTree = fAr->TreeR();
-    if(!pTree){
-       Warning("Exec","Error getting TreeR. TreeR=%p",pTree);
-       return;
-    } // end if !pTree
-    TBranch *branch = 0;
-    if(fSlowFast) sprintf(name,"ITSRecPointsF");
-    else sprintf(name,"ITSRecPoints");
-    branch = pTree->GetBranch(name);
-    if(!branch){
-       Warning("Exec","Can't find branch %s in TreeR fSlowFast=%d",
-               name,fSlowFast);
-       return;
-    } // end if !branch
-    TClonesArray *points = new TClonesArray("AliITSRecPoint",10000);
-    branch->SetAddress(&points);
-    Int_t nentr = (Int_t) branch->GetEntries();
-
-    if(fOut) fOut->cd();
-    TClonesArray *cluster = new TClonesArray("AliITSclusterV2",10000);
-    sprintf(name,"TreeC_ITS_%d",fAr->GetHeader()->GetEvent());
-    TTree *cTree = new TTree(name,"ITS clusters");
-    cTree->Branch("Clusters",&cluster);
-    TClonesArray &cl = *cluster;
-
-    Float_t lp[5];
-    Int_t lab[6];
-    Int_t i,j,lay,lad,det,nclusters=0,ncl;
-    Float_t kmip,x,y,zshift,yshift;
-    Double_t rot[9];
-    AliITSRecPoint *p;
-    TParticle *part;
-
-    for(i=0;i<nentr;i++){
-       points->Clear();
-       branch->GetEvent(i);
-       ncl = points->GetEntriesFast();
-       if(ncl<=0) {
-           cTree->Fill();
-           continue;
-       } // end if ncl<=0
-       nclusters += ncl;
-       fGeom->GetModuleId(i,lay,lad,det);
-       fGeom->GetTrans(i,x,y,zshift);
-       fGeom->GetRotMatrix(i,rot);
-       yshift = x*rot[0] + y*rot[1];
-       kmip = 1.0; // fGeom->GetModuletype(i)==kSPD
-       if(fGeom->GetModuleType(i)==kSDD) kmip = 280.0;
-       if(fGeom->GetModuleType(i)==kSSD) kmip = 38.0;
-       for(j=0;j<ncl;j++){
-           p = (AliITSRecPoint*)(points->UncheckedAt(j));
-           lp[0] = - (p->GetDetLocalX() + yshift);
-           if(lay==1) lp[0] = -lp[0];
-           lp[1] = p->GetZ() + zshift;
-           lp[2] = p->GetSigmaDetLocX2();
-           lp[3] = p->GetSigmaZ2();
-           lp[4] = p->GetQ()/kmip;
-           lab[0] = p->GetLabel(0);
-           lab[1] = p->GetLabel(1);
-           lab[2] = p->GetLabel(2);
-           lab[3] = i;
-           lad = lab[0];
-           if(lad>=0){
-               part = (TParticle*) fAr->GetMCApp()->Particle(lad);
-               lad = -3;
-               while(part->P() < 0.005){
-                   if(part->GetFirstMother()<0){
-                       Warning("Exec","Primary momentum: %e",part->P());
-                       break;
-                   } // end if part->GetFirstMother()<0
-                   lad = part->GetFirstMother();
-                   part = (TParticle*) fAr->GetMCApp()->Particle(lad);
-               } // end while part->P() < 0.005
-               if(lab[1]<0) lab[1] = lad;
-               else if(lab[2]<0) lab[2] = lad;
-               else Warning("Exec","No empty lables!");
-           } // end if lab>=0
-           Int_t info[3]={0,0,0};
-           new(cl[j]) AliITSclusterV2(lab,lp,info);
-       } // end for j
-       cTree->Fill();
-       cluster->Delete();
-       points->Delete();
-    } // end for i
-    cTree->Write();
-    
-    delete cTree;
-    delete cluster;
-    delete points;
-}
diff --git a/ITS/AliITSFindClustersV2.h b/ITS/AliITSFindClustersV2.h
deleted file mode 100644 (file)
index fb1c62f..0000000
+++ /dev/null
@@ -1,50 +0,0 @@
-#ifndef ALIITSFINDCLUSTERSV2_H
-#define ALIITSFINDCLUSTERSV2_H
-/* Copyright (c) 1998-2001, ALICE Experiment at CERN, All rights reserved *
- * See cxx source for full Copyright notice                               */
-/*
-  $Id$
- */
-///////////////////////////////////////////////////////////////////
-//Class for reconstruction of clusters V2                        //
-///////////////////////////////////////////////////////////////////
-#include <TTask.h>
-
-class TString;
-class TFile;
-class AliRun;
-class AliITSgeom;
-
-class AliITSFindClustersV2 : public TTask{
- public:
-    AliITSFindClustersV2(); // default constructor
-    // standard constructor files not opened, by file name
-    AliITSFindClustersV2(const TString infile,const TString outfile = "");
-    // standard constructor for files already opened.
-    AliITSFindClustersV2(TFile *in,TFile *out=0);
-    // Standard constructor for AliRun already read in.
-    AliITSFindClustersV2(AliRun *ar, const TString outfile = "");
-    AliITSFindClustersV2(const AliITSFindClustersV2& rec);
-    AliITSFindClustersV2& operator=(const AliITSFindClustersV2 &source);
-
-    virtual ~AliITSFindClustersV2();//Destructor
-    virtual Bool_t FastSimulation() const {return fSlowFast;}
-    virtual void SetSlowSimulation(){fSlowFast = kFALSE;}
-    virtual void SetFastSimulation(){fSlowFast = kTRUE;}
-    virtual void Exec(const Option_t *opt="");
- private:
-    AliRun *fAr;           //! Pointer of type AliRun
-    Bool_t fDeletfAr;      //! Logical to indecate if fAr should be deleted.
-    AliITSgeom *fGeom;     //! Pointer to ITS geometry
-    TString *fInFileName;  //! Pointer to input file name string.
-    TString *fOutFileName; //! Pointer to output file name string.
-    TFile   *fIn;          //! pointer to input file
-    TFile   *fOut;         //! pointer to output file
-    Bool_t fInit;          //! true if Init was successfull
-    Bool_t fSlowFast;      //! if true then using fast ITS simulation.
-
-    ClassDef(AliITSFindClustersV2,1) // Task to Reconstruct ITS from Digits.
-
-};
-#endif
index 99ba0a6..8f067d3 100644 (file)
@@ -54,7 +54,6 @@
 #pragma link C++ class AliITStrackerANN+;
 // Tasks
 #pragma link C++ class AliITSreconstruction+;
-#pragma link C++ class AliITSFindClustersV2+;
 //#pragma link C++ class DisplayITSv11+;
 
 #pragma link C++ class AliITSclusterTable+;
index c7aaa35..bc67496 100644 (file)
@@ -27,7 +27,6 @@ SRCS =        AliITSDetTypeRec.cxx \
                AliITSPid.cxx \
                AliITStrackV2Pid.cxx \
                AliITSreconstruction.cxx \
-               AliITSFindClustersV2.cxx \
                AliITSRiemannFit.cxx \
                AliITSclusterTable.cxx \
                AliITStrackSA.cxx \