@@ -276,36 +276,36 @@ struct qVectorsCorrection {
276276
277277 if (nmode == 2 ) {
278278 for (int ishift = 1 ; ishift <= 10 ; ishift++) {
279- 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 )));
280- 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 )));
279+ 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])));
280+ 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])));
281281
282- 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 )));
283- 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 )));
282+ 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])));
283+ 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])));
284284
285- 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 )));
286- 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 )));
285+ 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])));
286+ 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])));
287287 }
288288 } else if (nmode == 3 ) {
289289 for (int ishift = 1 ; ishift <= 10 ; ishift++) {
290- 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 )));
291- 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 )));
290+ 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])));
291+ 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])));
292292
293- 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 )));
294- 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 )));
293+ 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])));
294+ 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])));
295295
296- 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 )));
297- 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 )));
296+ 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])));
297+ 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])));
298298 }
299299 } else if (nmode == 4 ) {
300300 for (int ishift = 1 ; ishift <= 10 ; ishift++) {
301- 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 )));
302- 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 )));
301+ 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])));
302+ 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])));
303303
304- 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 )));
305- 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 )));
304+ 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])));
305+ 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])));
306306
307- 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 )));
308- 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 )));
307+ 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])));
308+ 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])));
309309 }
310310 }
311311 }
0 commit comments