git://git.uio.no / u / mrichter / AliRoot.git / blame
| Commit | Line | Data |
|---|---|---|
| 00c77045 | 1 | // $Id$ |
| 2 | // | |
| 3 | // Emcal particle trigger class, which can contain either | |
| 4 | // | |
| 5 | // Author: J.Kral | |
| 6 | ||
| 7 | #include "AliEmcalTriggerPatchInfo.h" | |
| 8 | #include "AliLog.h" | |
| cf9ee6ed | 9 | #include "AliEMCALGeometry.h" |
| 10 | #include "TArrayI.h" | |
| 00c77045 | 11 | |
| 12 | //_________________________________________________________________________________________________ | |
| 13 | AliEmcalTriggerPatchInfo::AliEmcalTriggerPatchInfo() : | |
| 14 | TObject(), | |
| 15 | fCenterGeo(), | |
| 16 | fCenterMass(), | |
| 17 | fEdge1(), | |
| 18 | fEdge2(), | |
| cf9ee6ed | 19 | fADCAmp(0), |
| 7ea0d3f4 MV |
20 | fTriggerBits(0), |
| 21 | fOffSet(kTriggerTypeEnd) | |
| 00c77045 | 22 | { |
| 23 | // Default constructor. | |
| cf9ee6ed | 24 | fEdgeCell[0] = -1; |
| 25 | fEdgeCell[1] = -1; | |
| 00c77045 | 26 | } |
| 27 | ||
| 28 | ||
| 29 | //_________________________________________________________________________________________________ | |
| 30 | AliEmcalTriggerPatchInfo::AliEmcalTriggerPatchInfo(const AliEmcalTriggerPatchInfo &p) : | |
| 31 | TObject(p), | |
| 32 | fCenterGeo(p.fCenterGeo), | |
| 33 | fCenterMass(p.fCenterMass), | |
| 34 | fEdge1(p.fEdge1), | |
| 35 | fEdge2(p.fEdge2), | |
| 36 | fADCAmp(p.fADCAmp), | |
| 7ea0d3f4 MV |
37 | fTriggerBits(p.fTriggerBits), |
| 38 | fOffSet(p.fOffSet) | |
| 00c77045 | 39 | { |
| 40 | // Copy constructor. | |
| cf9ee6ed | 41 | fEdgeCell[0] = p.fEdgeCell[0]; |
| 42 | fEdgeCell[1] = p.fEdgeCell[1]; | |
| 00c77045 | 43 | } |
| 44 | ||
| 45 | //_________________________________________________________________________________________________ | |
| 46 | AliEmcalTriggerPatchInfo::~AliEmcalTriggerPatchInfo() | |
| 47 | { | |
| 48 | // Destructor. | |
| 49 | } | |
| 50 | ||
| 51 | //_________________________________________________________________________________________________ | |
| 52 | AliEmcalTriggerPatchInfo &AliEmcalTriggerPatchInfo::operator=(const AliEmcalTriggerPatchInfo &p) | |
| 53 | { | |
| 54 | // Assignment operator. | |
| 55 | ||
| 56 | if (this != &p) { | |
| 57 | fCenterGeo = p.fCenterGeo; | |
| 58 | fCenterMass = p.fCenterMass; | |
| 59 | fEdge1 = p.fEdge1; | |
| 60 | fEdge2 = p.fEdge2; | |
| 61 | fADCAmp = p.fADCAmp; | |
| 62 | fTriggerBits = p.fTriggerBits; | |
| cf9ee6ed | 63 | fEdgeCell[0] = p.fEdgeCell[0]; |
| 64 | fEdgeCell[1] = p.fEdgeCell[1]; | |
| 00c77045 | 65 | } |
| 66 | ||
| 67 | return *this; | |
| 68 | } | |
| 69 | ||
| cf9ee6ed | 70 | //_________________________________________________________________________________________________ |
| 71 | void AliEmcalTriggerPatchInfo::GetCellIndices( AliEMCALGeometry *geom, TArrayI *cells ){ | |
| 72 | ||
| 73 | // return cell indices of the given patch in hte cell array | |
| 74 | Int_t globCol, globRow, i, j, k, absId, cellAbsId[4];; | |
| 75 | ||
| 76 | cells->Set( 1024 ); | |
| 77 | ||
| 78 | // get corner, convert from cells to trigger channels | |
| 79 | globCol = GetEdgeCellX() / 2; | |
| 80 | globRow = GetEdgeCellY() / 2; | |
| 81 | ||
| 82 | // get the absolute trigger ID | |
| 83 | geom->GetAbsFastORIndexFromPositionInEMCAL( globCol, globRow, absId ); | |
| 84 | // convert to the 4 absId of the cells composing the trigger channel | |
| 85 | geom->GetCellIndexFromFastORIndex( absId, cellAbsId ); | |
| 86 | ||
| 87 | // sum the available energy in the 32/32 window of cells | |
| 88 | // step over trigger channels and get all the corresponding cells | |
| 89 | for( i = 0; i < 16; i++ ){ | |
| 90 | for( j = 0; j < 16; j++ ){ | |
| 91 | // get the 4 cells composing the trigger channel | |
| 92 | geom->GetAbsFastORIndexFromPositionInEMCAL( globCol+i, globRow+j, absId ); | |
| 93 | geom->GetCellIndexFromFastORIndex( absId, cellAbsId ); | |
| 94 | // add amplitudes and find patch edges | |
| 95 | for( k = 0; k < 4; k++ ){ | |
| 96 | cells->SetAt( cellAbsId[k], i*16*4+j*4+k ); | |
| 97 | } | |
| 98 | } | |
| 99 | } // 32x32 cell window | |
| 100 | ||
| 101 | ||
| 102 | } | |
| 103 | ||
| 104 | ||
| 00c77045 | 105 | //_________________________________________________________________________________________________ |
| 106 | void AliEmcalTriggerPatchInfo::SetLorentzVector( TLorentzVector &lv, TVector3 &v, Double_t e ){ | |
| 107 | // sets the vector | |
| 108 | Double_t r = TMath::Sqrt(v[0]*v[0]+v[1]*v[1]+v[2]*v[2] ) ; | |
| 109 | ||
| 110 | lv.SetPxPyPzE( e*v[0]/r, e*v[1]/r, e*v[2]/r, e) ; | |
| 111 | } | |
| 112 |