StRoot: StPmdClusterMaker/StPmdClustering.h Source File

00001 
00005 /******************************************************
00006  *
00007  * $Id: StPmdClustering.h,v 1.11 2010/05/29 00:47:10 rashmi Exp $
00008  *
00009  * Author: Dr. S.C. Phatak
00010  *         Dipak Mishra
00011  ******************************************************
00012  *
00013  * Description: Base class for PMD clusters
00014  *
00015  * $Log: StPmdClustering.h,v $
00016  * Revision 1.11  2010/05/29 00:47:10  rashmi
00017  * Added a call to new clustering routines in StPmdClustering
00018  *
00019  * Revision 1.10  2010/04/15 06:52:32  rashmi
00020  * Clustering with option to turn calibration refineclustering on/off
00021  *
00022  * Revision 1.9  2007/11/02 11:00:14  rashmi
00023  * Applying hitcalibration; eta,phi wrt primary vertex
00024  *
00025  * Revision 1.8  2007/08/31 10:52:00  rashmi
00026  * Defined cutoff, setting it using inline SetAdcCutOff()
00027  *
00028  * Revision 1.7  2004/09/07 23:02:43  fisyak
00029  * Add missing default ctors
00030  *
00031  * Revision 1.6  2004/07/19 13:23:27  subhasis
00032  * checks applied on clust_cell dimension
00033  *
00034  * Revision 1.5  2004/06/24 13:43:02  subhasis
00035  * several changes in clustering code
00036  *
00037  *
00038  ******************************************************/
00039 #ifndef STAR_StPmdClustering
00040 #define STAR_StPmdClustering
00041 #include <TH2.h>
00042 #include<math.h>
00043 #include <TH1.h>
00044 #include <TCanvas.h>
00045 #include <TArrayF.h>
00046 #include <TArrayI.h>
00047 #include "StPmdAbsClustering.h"
00048 #include "StThreeVectorF.hh"
00049 class StPmdHit;
00050 class StPmdCluster;
00051 class StPmdModule;
00052 class StPmdDetector;
00053 class StPmdClusterCollection;
00054 class StPmdClustering:public StPmdAbsClustering
00055 {
00056 
00057   private:
00058  
00059   // These options are passed to this maker by StPmdClusterMaker
00060   Bool_t mOptCalibrate;
00061   Bool_t mOptSimulate;
00062   Bool_t mOptRefineCluster;
00063   Double_t cutoff;
00064   StThreeVectorF mVertexPos;
00065 
00066   StPmdClusterCollection * pmdclus;
00067 
00068  //Made static so that I to declare this only once and not for each module
00069   static const Int_t NColumnMax = 96;; //Maximum number of cell along x axis
00070   static const Int_t NRowMax = 72; //Maximum number of cells along y axis
00071   // static StPmdHit* hitmap[NRowMax][NColumnMax];//stores addresses of hits  
00072   static const Int_t CutOff = 7;  
00073   static const Int_t MaxHits = 15000; // max hits in an event
00074   static const Int_t MaxSuperSize = 2000;
00075   static const Int_t MaxLocalPeaks = 700;
00076   static const Int_t MaxClusterSize = 2000;
00077   static const Int_t MaxNeighLP = 35;
00078   //  static const Float_t logweight = 4.5;
00079   //  Float_t logweight;
00080   static const Int_t debug = 0; // 0 for no cout; 1 for cout;
00081 
00082   //  TClonesArray * mPmdSuperClusters;
00083   void CellXY(Int_t, Int_t, Float_t&, Float_t&);
00084 
00085   //checks whether the hit position in  question is within the SM boundary
00086   Int_t CheckBoundary(Int_t, Int_t, Int_t); 
00087 
00088 
00089  protected:
00090   
00091   
00092  public:
00093   // functions for clustering
00095   StPmdClustering(StPmdDetector *pmd_det = 0, StPmdDetector* cpv_det = 0);
00097   virtual ~StPmdClustering();
00098   
00099   // These are the new cluster maker functions
00100 
00102   //  void findPmdClusters();  // new method
00103   Bool_t findPmdClusters2(StPmdDetector *);
00104   // Make clusters our of superclusters => refclust
00105   Int_t MakeClusters(TClonesArray *);
00106   //Int_t MakeClusters(TClonesArray *, StPmdClusterCollection*);
00107   // Make SuperClusters => crclust
00108   Int_t MakeSuperClusters(Int_t, StPmdModule*,TClonesArray*);
00109   // Helping BreakSuperClusters
00110   Int_t GetLocalPeaks(Int_t*,TList *);
00111   // Breaks the supercluster by associating cells with each localPeak
00112   Bool_t BreakSuperCluster(Int_t, Int_t*, TList * );
00113   //Bool_t BreakSuperCluster(Int_t, Int_t*, TList * , StPmdClusterCollection *);
00114   Float_t BuildCluster(StPmdCluster*, Float_t*);
00115   // GetCluster Properties from hits =>CentroidCal
00116   void Cell_eta_phi(Float_t, Float_t, Float_t &, Float_t &);
00117 
00118   // These are the old cluster maker functions
00119 
00121   //  void findPmdClusters(); // old method
00122   void findPmdClusters(StPmdDetector *);
00123   
00125   Int_t crclust(Double_t , Double_t, Int_t, Int_t); 
00127   void refclust(StPmdDetector*,Int_t, Int_t, Int_t,StPmdClusterCollection*);
00128   
00130   void order(Int_t);  
00132   void arrange(Int_t);
00133   Double_t Dist(Double_t, Double_t, Double_t, Double_t);  
00134   
00136   //  void gaussfit(Int_t, Int_t, Double_t &, Double_t &, Double_t &, Double_t &, Double_t &,Double_t &,Double_t & ); 
00138   void printclust(Int_t,Int_t, StPmdCluster*); 
00140   Int_t CentroidCal(Int_t,Int_t,Double_t&,Double_t&,Double_t&,Double_t&,Double_t&,Double_t&,Double_t&,Double_t&,Double_t&);
00141   void SetAdcCutOff(Double_t adccutoff);
00142   void SetVertexPos(const StThreeVectorF&);
00143   void Cluster_Eta_Phi(Float_t, Float_t, Float_t,Float_t&, Float_t&);
00144   
00145 
00146   void SetOptCalibrate(Bool_t a=kTRUE){mOptCalibrate = a;}  // Default is on; YES Calibrate
00147   void SetOptSimulate(Bool_t a=kFALSE){mOptSimulate = a;}    // Default is off; No Simulation
00148   void SetOptRefineCluster(Bool_t a=kFALSE){mOptRefineCluster = a;} // Default is off; NO Refine Clustering
00149   
00150   
00152   //  StPmdHit* GetHit(StPmdDetector*, Int_t, Int_t, Int_t);
00153   StPmdHit* GetHit(StPmdDetector*, Int_t, Double_t, Double_t);
00154   
00155   ClassDef(StPmdClustering, 1) 
00156     };
00157     
00158     inline void StPmdClustering::SetAdcCutOff(Double_t adccutoff){
00159       cutoff = adccutoff;
00160     }
00161 
00162 inline void StPmdClustering::SetVertexPos(const StThreeVectorF& vertexPos){
00163   mVertexPos = vertexPos;
00164 }
00165 
00166 
00167 
00168 
00169 #endif
00170 
00171 
00172 
00173 
00174