StRoot: StSpinPool/StJets/StJet.cxx Source File

00001 
00002 //
00003 // $Id: StJet.cxx,v 1.3 2009/09/05 22:15:16 pibero Exp $
00004 // $Log: StJet.cxx,v $
00005 // Revision 1.3  2009/09/05 22:15:16  pibero
00006 // Add tracks and towers references to jet
00007 //
00008 // Revision 1.2  2009/09/05 18:18:05  pibero
00009 // Add utility functions for trigger patches
00010 //
00011 // Revision 1.1  2008/06/01 03:41:44  tai
00012 // moved StJet, StJets, and TrackToJetIndex to StSpinPool/StJets
00013 //
00014 // Revision 1.4  2008/03/27 02:25:04  tai
00015 // moved the definitions of the construcors from .h to .cxx
00016 // set jetEt and 3 other variables in a constructor
00017 //
00018 // Revision 1.3  2007/06/07 01:36:06  kocolosk
00019 // fix minor AutoBuild warnings
00020 //
00021 // Revision 1.2  2007/06/05 21:56:13  kocolosk
00022 // added data members for zVertex and geometric trigger associations, plus methods for detEta (barrel only)
00023 //
00024 // Revision 1.1  2004/07/08 15:41:03  mmiller
00025 // First release of StJetMaker.  Mike's commit of full clean of StJetFinder, StJetMaker, and StSpinMaker.  builds and runs in dev.
00026 //
00027 // Revision 1.5  2003/09/11 18:14:18  thenry
00028 // *** empty log message ***
00029 //
00030 // Revision 1.4  2003/09/02 17:59:01  perev
00031 // gcc 3.2 updates + WarnOff
00032 //
00033 // Revision 1.3  2002/12/04 20:28:07  thenry
00034 // StppuDstMaker was modified to allow multiple jet analysis modules to be
00035 // run simultaneosly with various parameters while the Maker loads the events
00036 // and analyses them.  Four different jet analyzers exist:
00037 //
00038 // Konstanin's Analyzers:
00039 //     Kt type: StppKonstKtJetAnalyzer
00040 //     Cone type: StppKonstConeJetAnalyzer
00041 //
00042 // Mike's Analyzers:
00043 //     Kt type: StppMikeKtJetAnalyzer
00044 //     Cone type: StppMikeConeJetAnalyzer
00045 //
00046 // These modules all require the StJetFinder modules.
00047 //
00048 // Revision 1.2  2002/06/24 13:22:59  akio
00049 // numerous bug fix & updates
00050 //
00051 // Revision 1.1  2002/02/11 20:30:48  akio
00052 // Many updates, including very first version of jet finder.
00053 //
00054 //
00055 // Revision 1.0  2001/06/14 Akio Ogawa
00056 //
00058 
00059 #include <iostream>
00060 #include <stdio.h>
00061 
00062 #include "TrackToJetIndex.h"
00063 #include "StJet.h"
00064 
00065 ClassImp(StJet);
00066 
00067 StJet::StJet()
00068   : TLorentzVector(0,0,0,0)
00069   , nCell(0)
00070   , charge(0)
00071   , nTracks(0)
00072   , nBtowers(0)
00073   , nEtowers(0)
00074   , tpcEtSum(0.0)
00075   , btowEtSum(0.0)
00076   , etowEtSum(0.0)
00077   , jetEt(0.0)
00078   , jetPt(0.0)
00079   , jetEta(0.0)
00080   , jetPhi(0.0)
00081   , zVertex(-999)
00082 {
00083 }
00084 
00085 StJet::StJet(double lE, double lpx, double lpy, double lpz, Int_t size, int c)
00086   : TLorentzVector(lpx, lpy, lpz, lE)
00087   , nCell(size)
00088   , charge(c)
00089   , nTracks(0)
00090   , nBtowers(0)
00091   , nEtowers(0)
00092   , tpcEtSum(0.0)
00093   , btowEtSum(0.0)
00094   , etowEtSum(0.0)
00095   , zVertex(-999)
00096 {
00097   jetEt = Et();
00098   jetPt = Pt();
00099   jetEta = Eta();
00100   jetPhi = Phi();
00101 }
00102 
00103 TrackToJetIndex* StJet::leadingChargedParticle() const
00104 {
00105   TrackToJetIndex* lcp = 0;
00106   for (int iTrack = 0; iTrack < numberOfTracks(); ++iTrack) {
00107     TrackToJetIndex* t = track(iTrack);
00108     if (!lcp || t->Pt() > lcp->Pt()) lcp = t;
00109   }
00110   return lcp;
00111 }
00112 
00113 Float_t StJet::detEta() const {
00114     if (zVertex > -998) return detEta(zVertex);
00115     return -999;
00116 }
00117 
00118 Float_t StJet::detEta(float vz, float r) const {
00119     float hold(0.),denom(0.);
00120     if (Theta()==TMath::PiOver2()) { // if Jet Theta = 90 then tan is undefined
00121         if (vz==0) {hold = TMath::PiOver2();}
00122         else {hold = atan2(r,vz);}
00123     }
00124     else
00125     {
00126         denom = (r/tan(Theta()))+vz;
00127         if (denom==0.) {hold = TMath::PiOver2();}
00128         if (denom!=0.) {hold = atan2(r,denom);}
00129     }
00130     return -TMath::Log(TMath::Tan(hold/2));
00131 }
00132 
00133 void StJet::addGeomTrigger(int trigId) {
00134     //check if the trigId is already in the vector -- not really necessary
00135     for(unsigned int i=0; i<mGeomTriggers.size(); i++) {
00136         if (mGeomTriggers[i] == trigId) return;
00137     }
00138     mGeomTriggers.push_back(trigId);
00139 }
00140 
00141 bool StJet::geomTrigger(int trigId) const {
00142     for(unsigned int i=0; i<mGeomTriggers.size(); i++) {
00143         if (mGeomTriggers[i] == trigId) return true;
00144     }
00145     return false;
00146 }
00147 
00148 bool StJet::getJetPatchEtaPhi(int jetPatch, float& eta, float& phi)
00149 {
00150   //
00151   // Pibero Djawotho <pibero@tamu.edu>
00152   // Texas A&M University
00153   // 5 September 2009
00154   //
00155 
00156   // Sanity check
00157   if (jetPatch < 0 || jetPatch >= 12) return false;
00158 
00159   //
00160   // The jet patches are numbered starting with JP0 centered at 150 degrees
00161   // looking from the West into the IR (intersection region) and increasing
00162   // clockwise, i.e. JP1 at 90 degrees, JP2 at 30 degrees, etc. On the East
00163   // side the numbering picks up at JP6 centered again at 150 degrees and
00164   // increasing clockwise (again as seen from the *West* into the IR). Thus
00165   // JP0 and JP6 are in the same phi location in the STAR coordinate system.
00166   // So are JP1 and JP7, etc.
00167   //
00168   // Jet Patch#  Eta  Phi  Quadrant
00169   //
00170   //          0  0.5  150       10'
00171   //          1  0.5   90       12'
00172   //          2  0.5   30        2'
00173   //          3  0.5  -30        4'
00174   //          4  0.5  -90        6'
00175   //          5  0.5 -150        8'
00176   //          6 -0.5  150       10'
00177   //          7 -0.5   90       12'
00178   //          8 -0.5   30        2'
00179   //          9 -0.5  -30        4'
00180   //         10 -0.5  -90        6'
00181   //         11 -0.5 -150        8'
00182   //
00183   // http://www.nikhef.nl/~ogrebeny/emc/files/Towers%20Layout.pdf
00184   // http://www.nikhef.nl/~ogrebeny/emc/files/BEMC.pdf
00185   // http://drupal.star.bnl.gov/STAR/system/files/BEMC_y2004.pdf
00186   //
00187   eta = (jetPatch < 6) ? 0.5 : -0.5;
00188   phi = 150 - (jetPatch % 6) * 60; // Degrees
00189 
00190   // Degrees to radians
00191   phi *= TMath::DegToRad();
00192 
00193   // Map phi into [-pi,pi]
00194   phi = TVector2::Phi_mpi_pi(phi);
00195 
00196   return true;
00197 }
00198 
00199 bool StJet::getJetPatchId(float eta, float phi, int& id)
00200 {
00201   //
00202   // Pibero Djawotho <pibero@tamu.edu>
00203   // Texas A&M University
00204   // 5 September 2009
00205   //
00206 
00207   // Jet patch id is left at -1 on failure
00208   id = -1;
00209 
00210   // Check range of eta
00211   if (eta < -1 || eta > 1) return false;
00212 
00213   // Map phi into [-pi,pi] if necessary
00214   if (phi < -M_PI || phi > M_PI) phi = TVector2::Phi_mpi_pi(phi);
00215 
00216   // Get jet patch id
00217   static const double PI_OVER_3 = M_PI/3;
00218 
00219   if (0 <= eta && eta <= 1) {
00220     if ( 2*PI_OVER_3 <= phi && phi <         M_PI) id = 0;
00221     if (   PI_OVER_3 <= phi && phi <  2*PI_OVER_3) id = 1;
00222     if (           0 <= phi && phi <    PI_OVER_3) id = 2;
00223     if (  -PI_OVER_3 <= phi && phi <            0) id = 3;
00224     if (-2*PI_OVER_3 <= phi && phi <   -PI_OVER_3) id = 4;
00225     if (       -M_PI <= phi && phi < -2*PI_OVER_3) id = 5;
00226   }
00227 
00228   if (-1 <= eta && eta < 0) {
00229     if ( 2*PI_OVER_3 <= phi && phi <         M_PI) id = 6;
00230     if (   PI_OVER_3 <= phi && phi <  2*PI_OVER_3) id = 7;
00231     if (           0 <= phi && phi <    PI_OVER_3) id = 8;
00232     if (  -PI_OVER_3 <= phi && phi <            0) id = 9;
00233     if (-2*PI_OVER_3 <= phi && phi <   -PI_OVER_3) id = 10;
00234     if (       -M_PI <= phi && phi < -2*PI_OVER_3) id = 11;
00235   }
00236 
00237   return (0 <= id && id < 12);
00238 }