@@ -173,17 +173,17 @@ void PlotPbLumi(int runNumber = 572073, bool sum = 1, double cut = 0, bool qc =
173173 Double_t * tcetoznc = tcetozncvec .data ();
174174 Double_t * vchtoznc = vchtozncvec .data ();
175175 for (int i = i0 ; i < ilast ; i ++ ) {
176- // for (int i = 30; i < 40; i++) {
176+ // for (int i = 30; i < 40; i++) {
177177
178178 int iv = i - i0 ;
179179 x [iv ] = (double_t )(recs [i + 1 ].intRecord .orbit + recs [i ].intRecord .orbit ) / 2. - orbit0 ;
180180 x [iv ] *= 88e-6 ;
181181 // x[i] = (double_t)(recs[i+1].epochTime + recs[i].epochTime)/2.;
182182 double_t tt = (double_t )(recs [i + 1 ].intRecord .orbit - recs [i ].intRecord .orbit );
183183 tt = tt * 88e-6 ;
184- //std::cout << i << " " << iv << " " << tt << std::endl;
184+ // std::cout << i << " " << iv << " " << tt << std::endl;
185185 //
186- // std::cout << recs[i+1].scalersInps[25] << std::endl;
186+ // std::cout << recs[i+1].scalersInps[25] << std::endl;
187187 double_t znci = (double_t )(recs [i + 1 ].scalersInps [25 ] - recs [i ].scalersInps [25 ]);
188188 double_t mu = - TMath ::Log (1. - znci / tt / nbc / frev );
189189 double_t zncipp = mu * nbc * frev ;
@@ -202,14 +202,14 @@ void PlotPbLumi(int runNumber = 572073, bool sum = 1, double cut = 0, bool qc =
202202 // std::cout << recs[i+1].scalers[tce].lmBefore << std::endl;
203203 had += recs [i + 1 ].scalers [tsc ].lmBefore - recs [i ].scalers [tsc ].lmBefore ;
204204 // rat = (double_t)(had)/double_t(recs[i+1].scalersInps[25] - recs[i].scalersInps[25])*28;
205- if (zncpp [iv ] > cut ) {
205+ if (zncpp [iv ] > cut ) {
206206 tcetsctoznc [iv ] = (double_t )(had ) / zncpp [iv ] / tt ;
207207 } else {
208208 tcetsctoznc [iv ] = 0. ;
209209 }
210210 had = recs [i + 1 ].scalers [tce ].lmBefore - recs [i ].scalers [tce ].lmBefore ;
211211 // rat = (double_t)(had)/double_t(recs[i+1].scalersInps[25] - recs[i].scalersInps[25])*28;
212- if (zncpp [iv ] > cut ) {
212+ if (zncpp [iv ] > cut ) {
213213 tcetoznc [iv ] = (double_t )(had ) / zncpp [iv ] / tt ;
214214 } else {
215215 tcetoznc [iv ] = 0. ;
@@ -218,7 +218,7 @@ void PlotPbLumi(int runNumber = 572073, bool sum = 1, double cut = 0, bool qc =
218218 double_t muvch = - TMath ::Log (1. - had / tt / nbc / frev );
219219
220220 // rat = (double_t)(had)/double_t(recs[i+1].scalersInps[25] - recs[i].scalersInps[25])*28;
221- if (zncpp [iv ] > cut ) {
221+ if (zncpp [iv ] > cut ) {
222222 vchtoznc [iv ] = (double_t )(had ) / zncpp [iv ] / tt ;
223223 } else {
224224 vchtoznc [iv ] = 0. ;
0 commit comments