@@ -275,36 +275,36 @@ struct qVectorsCorrection {
275275
276276 if (nmode == 2 ) {
277277 for (int ishift = 1 ; ishift <= 10 ; ishift++) {
278- histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * DetId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd]) / static_cast < float >(nmode )));
279- histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * DetId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd]) / static_cast < float >(nmode )));
278+ histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * DetId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd])));
279+ histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * DetId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd])));
280280
281- histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * RefAId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd]) / static_cast < float >(nmode )));
282- histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * RefAId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd]) / static_cast < float >(nmode )));
281+ histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * RefAId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd])));
282+ histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * RefAId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd])));
283283
284- histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * RefBId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd]) / static_cast < float >(nmode )));
285- histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * RefBId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd]) / static_cast < float >(nmode )));
284+ histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * RefBId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd])));
285+ histosQA.fill (HIST (" histShiftV2" ), vec.cent (), 2.0 * RefBId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd])));
286286 }
287287 } else if (nmode == 3 ) {
288288 for (int ishift = 1 ; ishift <= 10 ; ishift++) {
289- histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * DetId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd]) / static_cast < float >(nmode )));
290- histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * DetId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd]) / static_cast < float >(nmode )));
289+ histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * DetId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd])));
290+ histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * DetId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd])));
291291
292- histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * RefAId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd]) / static_cast < float >(nmode )));
293- histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * RefAId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd]) / static_cast < float >(nmode )));
292+ histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * RefAId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd])));
293+ histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * RefAId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd])));
294294
295- histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * RefBId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd]) / static_cast < float >(nmode )));
296- histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * RefBId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd]) / static_cast < float >(nmode )));
295+ histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * RefBId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd])));
296+ histosQA.fill (HIST (" histShiftV3" ), vec.cent (), 2.0 * RefBId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd])));
297297 }
298298 } else if (nmode == 4 ) {
299299 for (int ishift = 1 ; ishift <= 10 ; ishift++) {
300- histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * DetId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd]) / static_cast < float >(nmode )));
301- histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * DetId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd]) / static_cast < float >(nmode )));
300+ histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * DetId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd])));
301+ histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * DetId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[DetInd], vec.qvecRe ()[DetInd])));
302302
303- histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * RefAId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd]) / static_cast < float >(nmode )));
304- histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * RefAId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd]) / static_cast < float >(nmode )));
303+ histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * RefAId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd])));
304+ histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * RefAId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[RefAInd], vec.qvecRe ()[RefAInd])));
305305
306- histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * RefBId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd]) / static_cast < float >(nmode )));
307- histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * RefBId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * static_cast < float >(nmode) * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd]) / static_cast < float >(nmode )));
306+ histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * RefBId + 0.5 , ishift - 0.5 , TMath::Sin (ishift * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd])));
307+ histosQA.fill (HIST (" histShiftV4" ), vec.cent (), 2.0 * RefBId + 1.5 , ishift - 0.5 , TMath::Cos (ishift * TMath::ATan2 (vec.qvecIm ()[RefBInd], vec.qvecRe ()[RefBInd])));
308308 }
309309 }
310310 }
0 commit comments