StEmcPosition.h

//******************************************************************************
//                                                                            
// StEmcPosition.h
//
// Authors: Marcia Maria de Moura
//
// Initial version: 2001/12/21
//
//******************************************************************************

#ifndef StEmcPosition_H
#define StEmcPosition_H

#include "TObject.h"

#include "StarClassLibrary/StThreeVectorD.hh"
#include "StarClassLibrary/StThreeVectorF.hh"
#include "StPhysicalHelixD.hh"

class StTrack;
class StMuTrack;
class StMcTrack;
class StEmcGeom;
class StVertex;
class StMcVertex;

class StEmcPosition : public TObject
{
  public:            

             StEmcPosition();
    virtual  ~StEmcPosition();

    Bool_t            projTrack(StThreeVectorD* atFinal,StThreeVectorD* momentumAtFinal,const StTrack* const track,Double_t magField,Double_t radius=225.405,Int_t option=1) const;    
    Bool_t            projTrack(StThreeVectorD* atFinal,StThreeVectorD* momentumAtFinal,const StMuTrack* const track,Double_t magField,Double_t radius=225.405,Int_t option=1) const;    
    Bool_t            projTrack(StThreeVectorD* atFinal,StThreeVectorD* momentumAtFinal,const StMcTrack* const mcTrack,Double_t magField,Double_t radius=225.405,Int_t option=1) const;  
    Bool_t            projTrack(StThreeVectorD* atFinal,StThreeVectorD* momentumAtFinal,const StPhysicalHelixD* const helix,Double_t magField,Double_t radius=225.405,Int_t option=1) const;  
    
    Bool_t            trackOnEmc(StThreeVectorD* position,StThreeVectorD* momentum,const StTrack* const track,Double_t magField,Double_t emcRadius=225.405) const;   
    Bool_t            trackOnEmc(StThreeVectorD* position,StThreeVectorD* momentum,const StMuTrack* const track,double magField,double emcRadius=225.405) const;   
    Bool_t            trackOnEmc(StThreeVectorD* position,StThreeVectorD* momentum,const StMcTrack* const mcTrack,Double_t magField,Double_t emcRadius=225.405) const; 
    
    Int_t             getTowerEtaPhi(Double_t eta, Double_t phi, Float_t* towerEta, Float_t* towerPhi) const; 

    Int_t             getNextTowerId(Float_t eta, Float_t phi, Int_t nTowersdEta, Int_t nTowersdPhi) const; 
    Int_t             getNextTowerId(Int_t softId, Int_t nTowersdEta, Int_t nTowersdPhi) const; 
    Int_t             getNextTowerId(Int_t m, Int_t e, Int_t s, Int_t nTowersdEta, Int_t nTowersdPhi) const; 
        Int_t             getNextId(Int_t det, Int_t m, Int_t e, Int_t s, Int_t nEta, Int_t nPhi) const; 
        Int_t getNextId(Int_t det, Int_t softId, Int_t nEta, Int_t nPhi)const;

        Float_t           getDistTowerToTrack(Double_t trackEta, Double_t trackPhi, Int_t nTowersdEta, Int_t nTowersdPhi) const; 

    StThreeVectorF    getPosFromVertex(const StVertex* const vertex, Int_t TowerId) const; 
    StThreeVectorF    getPosFromVertex(const StThreeVectorF& position, Int_t TowerId) const; 
    StThreeVectorF    getPosFromVertex(const StMcVertex* const vertex, Int_t TowerId) const; 
    Float_t           getThetaFromVertex(const StVertex* const vertex, Int_t TowerId) const; 
    Float_t           getThetaFromVertex(const StThreeVectorF& vertex, Int_t TowerId) const; 
    Float_t           getThetaFromVertex(const StMcVertex* const vertex, Int_t TowerId) const; 
    Float_t           getEtaFromVertex(const StVertex* const vertex, Int_t TowerId) const; 
    Float_t           getEtaFromVertex(const StThreeVectorF& vertex, Int_t TowerId) const; 
    Float_t           getEtaFromVertex(const StMcVertex* const vertex, Int_t TowerId) const; 
    Float_t           getPhiFromVertex(const StVertex* const vertex, Int_t TowerId) const; 
    Float_t           getPhiFromVertex(const StThreeVectorF& vertex, Int_t TowerId) const; 
    Float_t           getPhiFromVertex(const StMcVertex* const vertex, Int_t TowerId) const; 

  protected:     

    StEmcGeom* mGeom[4];   

  ClassDef(StEmcPosition,2)

};
#endif