SDL::segments& segmentsInGPU = (*event->getSegments()); SDL::miniDoublets& miniDoubletsInGPU = (*event->getMiniDoublets()); SDL::quintuplets& quintupletsInGPU = (*event->getQuintuplets()); SDL::triplets& tripletsInGPU = (*event->getTriplets()); float x1 = miniDoubletsInGPU.anchorX[firstMDIndex]; float x2 = miniDoubletsInGPU.anchorX[secondMDIndex]; float x3 = miniDoubletsInGPU.anchorX[thirdMDIndex]; float x4 = miniDoubletsInGPU.anchorX[fourthMDIndex]; float x5 = miniDoubletsInGPU.anchorX[fifthMDIndex]; float y1 = miniDoubletsInGPU.anchorY[firstMDIndex]; float y2 = miniDoubletsInGPU.anchorY[secondMDIndex]; float y3 = miniDoubletsInGPU.anchorY[thirdMDIndex]; float y4 = miniDoubletsInGPU.anchorY[fourthMDIndex]; float y5 = miniDoubletsInGPU.anchorY[fifthMDIndex]; float g, f; innerRadius = computeRadiusFromThreeAnchorHits(x1, y1, x2, y2, x3, y3, g, f); outerRadius = computeRadiusFromThreeAnchorHits(x3, y3, x4, y4, x5, y5, g, f); bridgeRadius = computeRadiusFromThreeAnchorHits(x2, y2, x3, y3, x4, y4, g, f); ana.tx->pushbackToBranch("t5_lst_dnn_input_0", 2.99792458e-3f*3.8f*innerRadius); // inner t3_pt ana.tx->pushbackToBranch("t5_lst_dnn_input_1", miniDoubletsInGPU.anchorEta[firstMDIndex]); // inner (hit 1) t3_0_eta ana.tx->pushbackToBranch("t5_lst_dnn_input_2", miniDoubletsInGPU.anchorPhi[firstMDIndex]); // inner (hit 1) t3_0_phi ana.tx->pushbackToBranch("t5_lst_dnn_input_3", miniDoubletsInGPU.anchorZ[firstMDIndex]); // inner (hit 1) t3_0_z ana.tx->pushbackToBranch("t5_lst_dnn_input_4", sqrtf(x1*x1 + y1*y1)); // inner (hit 1) t3_0_r ana.tx->pushbackToBranch("t5_lst_dnn_input_5", (float)modulesInGPU.layers[lowerModuleIndex1]); // inner (hit 1) t3_0_layer ana.tx->pushbackToBranch("t5_lst_dnn_input_6", miniDoubletsInGPU.anchorEta[secondMDIndex]); // inner (hit 2) t3_2_eta ana.tx->pushbackToBranch("t5_lst_dnn_input_7", miniDoubletsInGPU.anchorPhi[secondMDIndex]); // inner (hit 2) t3_2_phi ana.tx->pushbackToBranch("t5_lst_dnn_input_8", miniDoubletsInGPU.anchorZ[secondMDIndex]); // inner (hit 2) t3_2_z ana.tx->pushbackToBranch("t5_lst_dnn_input_9", sqrtf(x2*x2 + y2*y2)); // inner (hit 2) t3_2_r ana.tx->pushbackToBranch("t5_lst_dnn_input_10", (float)modulesInGPU.layers[lowerModuleIndex2], // inner (hit 2) t3_2_layer ana.tx->pushbackToBranch("t5_lst_dnn_input_11", miniDoubletsInGPU.anchorEta[thirdMDIndex]); // inner (hit 3) t3_4_eta ana.tx->pushbackToBranch("t5_lst_dnn_input_12", miniDoubletsInGPU.anchorPhi[thirdMDIndex]); // inner (hit 3) t3_4_phi ana.tx->pushbackToBranch("t5_lst_dnn_input_13", miniDoubletsInGPU.anchorZ[thirdMDIndex]); // inner (hit 3) t3_4_z ana.tx->pushbackToBranch("t5_lst_dnn_input_15", sqrtf(x3*x3 + y3*y3)); // inner (hit 3) t3_4_r ana.tx->pushbackToBranch("t5_lst_dnn_input_16", (float)modulesInGPU.layers[lowerModuleIndex3], // inner (hit 3) t3_4_layer ana.tx->pushbackToBranch("t5_lst_dnn_input_17", 2.99792458e-3f*3.8f*innerRadius); // outer t3_pt ana.tx->pushbackToBranch("t5_lst_dnn_input_18", miniDoubletsInGPU.anchorEta[thirdMDIndex]); // outer (hit 4) t3_0_eta ana.tx->pushbackToBranch("t5_lst_dnn_input_19", miniDoubletsInGPU.anchorPhi[thirdMDIndex]); // outer (hit 4) t3_0_phi ana.tx->pushbackToBranch("t5_lst_dnn_input_20", miniDoubletsInGPU.anchorZ[thirdMDIndex]); // outer (hit 3) t3_0_z ana.tx->pushbackToBranch("t5_lst_dnn_input_21", sqrtf(x3*x3 + y3*y3)); // outer (hit 3) t3_0_r ana.tx->pushbackToBranch("t5_lst_dnn_input_22", (float)modulesInGPU.layers[lowerModuleIndex3]); // outer (hit 3) t3_0_layer ana.tx->pushbackToBranch("t5_lst_dnn_input_23", miniDoubletsInGPU.anchorEta[fourthMDIndex]); // outer (hit 4) t3_2_eta ana.tx->pushbackToBranch("t5_lst_dnn_input_24", miniDoubletsInGPU.anchorPhi[fourthMDIndex]); // outer (hit 4) t3_2_phi ana.tx->pushbackToBranch("t5_lst_dnn_input_25", miniDoubletsInGPU.anchorZ[fourthMDIndex]); // outer (hit 4) t3_2_z ana.tx->pushbackToBranch("t5_lst_dnn_input_26", sqrtf(x4*x4 + y4*y4)); // outer (hit 4) t3_2_r ana.tx->pushbackToBranch("t5_lst_dnn_input_27", (float)modulesInGPU.layers[lowerModuleIndex4]); // outer (hit 4) t3_2_layer ana.tx->pushbackToBranch("t5_lst_dnn_input_28", miniDoubletsInGPU.anchorEta[fifthMDIndex]); // outer (hit 5) t3_4_eta ana.tx->pushbackToBranch("t5_lst_dnn_input_29", miniDoubletsInGPU.anchorPhi[fifthMDIndex]); // outer (hit 5) t3_4_phi ana.tx->pushbackToBranch("t5_lst_dnn_input_30", miniDoubletsInGPU.anchorZ[fifthMDIndex]); // outer (hit 5) t3_4_z ana.tx->pushbackToBranch("t5_lst_dnn_input_31", sqrtf(x5*x5 + y5*y5)); // outer (hit 5) t3_4_r ana.tx->pushbackToBranch("t5_lst_dnn_input_32", (float)modulesInGPU.layers[lowerModuleIndex5]); // outer (hit 5) t3_4_layer ana.tx->pushbackToBranch("t5_lst_dnn_input_33", (innerRadius+outerRadius)*3.8f*1.602f/(2*100*5.39f)); // t5_pt ana.tx->pushbackToBranch("t5_lst_dnn_input_34", miniDoubletsInGPU.anchorEta[segmentsInGPU.mdIndices[2*tripletsInGPU.segmentIndices[2*innerTripletIndex+layer2_adjustment]]]); // t5_eta ana.tx->pushbackToBranch("t5_lst_dnn_input_35", miniDoubletsInGPU.anchorPhi[segmentsInGPU.mdIndices[2*tripletsInGPU.segmentIndices[2*innerTripletIndex+layer2_adjustment]]]); // t5_phi ana.tx->pushbackToBranch("t5_lst_dnn_input_36", innerRadius); // t5_innerRadius ana.tx->pushbackToBranch("t5_lst_dnn_input_37", bridgeRadius); // t5_bridgeRadius ana.tx->pushbackToBranch("t5_lst_dnn_input_38", outerRadius); // t5_outerRadius