36 auto trackindex = std::make_unique<std::vector<int>>();
37 auto rawid = std::make_unique<std::vector<unsigned int>>();
38 auto localdirx = std::make_unique<std::vector<double>>();
39 auto localdiry = std::make_unique<std::vector<double>>();
40 auto localdirz = std::make_unique<std::vector<double>>();
41 auto firststrip = std::make_unique<std::vector<unsigned short>>();
42 auto nstrips = std::make_unique<std::vector<unsigned short>>();
43 auto saturation = std::make_unique<std::vector<bool>>();
44 auto overlapping = std::make_unique<std::vector<bool>>();
45 auto farfromedge = std::make_unique<std::vector<bool>>();
46 auto charge = std::make_unique<std::vector<unsigned int>>();
47 auto path = std::make_unique<std::vector<double>>();
48 #ifdef ExtendedCALIBTree
49 auto chargeoverpath = std::make_unique<std::vector<double>>();
51 auto amplitude = std::make_unique<std::vector<unsigned char>>();
52 auto gainused = std::make_unique<std::vector<double>>();
53 auto gainusedTick = std::make_unique<std::vector<double>>();
68 vector<TrajectoryMeasurement> measurements = traj->
measurements();
69 for (vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin();
70 measurement_it != measurements.end();
86 for (
unsigned int h = 0;
h < 2;
h++) {
87 if (!sistripmatchedhit &&
h == 1) {
89 }
else if (sistripmatchedhit &&
h == 0) {
92 }
else if (sistripmatchedhit &&
h == 1) {
96 }
else if (sistripsimplehit) {
97 StripCluster = (sistripsimplehit->
cluster()).
get();
99 }
else if (sistripsimple1dhit) {
100 StripCluster = (sistripsimple1dhit->
cluster()).
get();
102 }
else if (sipixelhit) {
110 double cosine = trackDirection.
z() / trackDirection.
mag();
111 bool Saturation =
false;
115 double PrevGain = -1;
116 double PrevGainTick = -1;
121 const auto& Ampls = StripCluster->
amplitudes();
123 NStrips = Ampls.size();
124 int APVId = FirstStrip / 128;
131 for (
unsigned int a = 0;
a < Ampls.size();
a++) {
140 if (FirstStrip == 128)
142 if (FirstStrip == 256)
144 if (FirstStrip == 384)
146 if (FirstStrip == 512)
148 if (FirstStrip == 640)
151 if (FirstStrip <= 127 && FirstStrip + Ampls.size() > 127)
153 if (FirstStrip <= 255 && FirstStrip + Ampls.size() > 255)
155 if (FirstStrip <= 383 && FirstStrip + Ampls.size() > 383)
157 if (FirstStrip <= 511 && FirstStrip + Ampls.size() > 511)
159 if (FirstStrip <= 639 && FirstStrip + Ampls.size() > 639)
162 if (FirstStrip + Ampls.size() == 127)
164 if (FirstStrip + Ampls.size() == 255)
166 if (FirstStrip + Ampls.size() == 383)
168 if (FirstStrip + Ampls.size() == 511)
170 if (FirstStrip + Ampls.size() == 639)
172 if (FirstStrip + Ampls.size() == 767)
174 }
else if (PixelCluster) {
175 const auto& Ampls = PixelCluster->
pixelADC();
178 FirstStrip = ((FirstRow / 80) << 3 | (FirstCol / 52)) * 128;
183 for (
unsigned int a = 0;
a < Ampls.size();
a++) {
189 #ifdef ExtendedCALIBTree
190 double ChargeOverPath = (double)
Charge /
Path;
194 rawid->push_back(
DetId);
195 localdirx->push_back(trackDirection.
x());
196 localdiry->push_back(trackDirection.
y());
197 localdirz->push_back(trackDirection.
z());
198 firststrip->push_back(FirstStrip);
199 nstrips->push_back(NStrips);
205 #ifdef ExtendedCALIBTree
206 chargeoverpath->push_back(ChargeOverPath);
208 gainused->push_back(PrevGain);
209 gainusedTick->push_back(PrevGainTick);
226 #ifdef ExtendedCALIBTree