28 #include "StHiMicroEvent.h"
33 TClonesArray* StHiMicroEvent::mSTracks = 0;
34 TClonesArray* StHiMicroEvent::mSHits = 0;
38 StHiMicroEvent::StHiMicroEvent()
46 Int_t hitSZ = trackSZ*50;
49 cout <<
"StHiMicroEvent::StHiMicroEvent" << endl
50 <<
"\tCreating the clonesarrays..." << endl;
53 mTracks =
new TClonesArray(
"StHiMicroTrack",trackSZ);
55 mHits =
new TClonesArray(
"StHiMicroHit",hitSZ);
57 if(DEBUG) cout <<
"\t...done" <<endl;
69 StHiMicroEvent::~StHiMicroEvent()
71 if(DEBUG) cout <<
"StHiMicroEvent::~StHiMicroEvent" << endl;
83 TClonesArray &tracks = *mTracks;
95 TClonesArray &hits = *mHits;
108 StHiMicroEvent::SetCentrality(Int_t N)
113 if (this->CenterOfMassEnergy() >= 199.) {
114 if (fabs(this->MagneticField()) >= 4.) {
115 cout <<
"StHiMicroEvent->SetCentralityYear2AuAuFull: " << N << endl;
119 int cent[] = {14,30,56,94,146,217,312,431,510};
120 if (N < cent[0]) { mCentrality = 0; }
121 else if (N < cent[1]) { mCentrality = 1; }
122 else if (N < cent[2]) { mCentrality = 2; }
123 else if (N < cent[3]) { mCentrality = 3; }
124 else if (N < cent[4]) { mCentrality = 4; }
125 else if (N < cent[5]) { mCentrality = 5; }
126 else if (N < cent[6]) { mCentrality = 6; }
127 else if (N < cent[7]) { mCentrality = 7; }
128 else if (N < cent[8]) { mCentrality = 8; }
129 else { mCentrality = 9; }
132 cout <<
"StHiMicroEvent->SetCentralityYear2AuAuHalf: " << N << endl;
136 int cent[] = {14,32,59,98,149,216,302,409,474};
137 if (N < cent[0]) { mCentrality = 0; }
138 else if (N < cent[1]) { mCentrality = 1; }
139 else if (N < cent[2]) { mCentrality = 2; }
140 else if (N < cent[3]) { mCentrality = 3; }
141 else if (N < cent[4]) { mCentrality = 4; }
142 else if (N < cent[5]) { mCentrality = 5; }
143 else if (N < cent[6]) { mCentrality = 6; }
144 else if (N < cent[7]) { mCentrality = 7; }
145 else if (N < cent[8]) { mCentrality = 8; }
146 else { mCentrality = 9; }
150 cout <<
"StHiMicroEvent->SetCentrality: " << N << endl;
152 cout <<
"WARNING: MAKE SURE THE RIGHT ETA RANGE IS DEFINED BEFORE USING THESE DEFINITIONS!!!" << endl;
157 int cent[] = {20,100,180,270,360,460,560,660,870};
158 if (N < cent[0]) { mCentrality = 0; }
159 else if (N < cent[1]) { mCentrality = 1; }
160 else if (N < cent[2]) { mCentrality = 2; }
161 else if (N < cent[3]) { mCentrality = 3; }
162 else if (N < cent[4]) { mCentrality = 4; }
163 else if (N < cent[5]) { mCentrality = 5; }
164 else if (N < cent[6]) { mCentrality = 6; }
165 else if (N < cent[7]) { mCentrality = 7; }
166 else if (N < cent[8]) { mCentrality = 8; }
167 else { mCentrality = 9; }
175 StHiMicroEvent::Clear(Option_t* option)
177 mTracks->Clear(option);
178 mHits->Clear(option);