@@ -59,7 +59,6 @@ struct PCMQCMC {
5959 Configurable<float > maxRgen{" maxRgen" , 90 .f , " maximum radius for generated particles" };
6060 Configurable<float > margin_z_mc{" margin_z_mc" , 7.0 , " margin for z cut in cm for MC" };
6161 Configurable<bool > cfgRequireTrueAssociation{" cfgRequireTrueAssociation" , false , " flag to require true mc collision association" };
62- Configurable<bool > cfg_fill_resolution{" cfg_fill_resoltion" , false , " flag to fill resolution histogram" };
6362
6463 EMPhotonEventCut fEMEventCut ;
6564 struct : ConfigurableGroup {
@@ -204,22 +203,9 @@ struct PCMQCMC {
204203 fRegistry .add (" V0/primary/hPtGen_DeltaEta" , " photon #eta resolution;p_{T}^{gen} (GeV/c);#eta^{rec} - #eta^{gen}" , kTH2F , {{1000 , 0 , 10 }, {100 , -0 .5f , 0 .5f }}, true );
205204 fRegistry .add (" V0/primary/hPtGen_DeltaPhi" , " photon #varphi resolution;p_{T}^{gen} (GeV/c);#varphi^{rec} - #varphi^{gen} (rad.)" , kTH2F , {{1000 , 0 , 10 }, {100 , -0 .5f , 0 .5f }}, true );
206205 fRegistry .add (" V0/primary/hRxyGen_DeltaPtOverPtGen" , " photon p_{T} resolution; R_{xy}^{gen} (cm);(p_{T}^{rec} - p_{T}^{gen})/p_{T}^{gen}" , kTH2F , {{100 , 0 , 100 }, {200 , -1 .0f , 1 .0f }}, true );
207- fRegistry .add (" V0/primary/hsPhotonResolution" ,
208- " Photon resolution;p_{T};#eta;R_{xy};Z_{conv};Z_{vtx};#Deltap_{T}/p_{T};#Delta#eta;#Delta#phi" ,
209- kTHnSparseF ,
210- {{100 , 0 , 10 },
211- {80 , -1.6 , 1.6 },
212- {100 , 0 , 100 },
213- {100 , -50 , 50 },
214- {100 , -50 , 50 },
215- {200 , -1 , 1 },
216- {100 , -0.5 , 0.5 },
217- {100 , -0.5 , 0.5 }},
218- false );
219206 fRegistry .add (" V0/primary/hRxyGen_DeltaEta" , " photon #eta resolution;R_{xy}^{gen} (cm);#eta^{rec} - #eta^{gen}" , kTH2F , {{100 , 0 , 100 }, {100 , -0 .5f , 0 .5f }}, true );
220207 fRegistry .add (" V0/primary/hRxyGen_DeltaPhi" , " photon #varphi resolution;R_{xy}^{gen} (cm);#varphi^{rec} - #varphi^{gen} (rad.)" , kTH2F , {{100 , 0 , 100 }, {100 , -0 .5f , 0 .5f }}, true );
221208 fRegistry .add (" V0/primary/hRxyGen_DeltaR" , " photon #varphi resolution;R_{xy}^{gen} (cm);#varphi^{rec} - #varphi^{gen} (rad.)" , kTH2F , {{100 , 0 , 100 }, {100 , 0 , 100 }}, true );
222- fRegistry .add (" V0/primary/hsConvVtxZPtR" , " z_{vtx} vs p_{T} vs R_{xy};z_{vtx} (cm);p_{T} (GeV/c);R_{xy} (cm)" , kTHnSparseF , {{100 , -20 .0f , +20 .0f }, {100 , 0 .0f , 10 .0f }, {100 , 0 , 100 }}, false );
223209 fRegistry .add (" V0/primary/hXY_MC" , " X vs. Y of true photon conversion point.;X (cm);Y (cm)" , kTH2F , {{400 , -100 .0f , +100 }, {400 , -100 , +100 }}, true );
224210 fRegistry .add (" V0/primary/hRZ_MC" , " R vs. Z of true photon conversion point;Z (cm);R_{xy} (cm)" , kTH2F , {{200 , -100 .0f , +100 }, {200 , 0 , 100 }}, true );
225211 fRegistry .add (" V0/primary/hsConvPoint" , " photon conversion point;r_{xy} (cm);#varphi (rad.);#eta;" , kTHnSparseF , {{100 , 0 .0f , 100 }, {90 , 0 , 2 * M_PI}, {80 , -2 , +2 }}, false );
@@ -253,21 +239,6 @@ struct PCMQCMC {
253239 fRegistry .add (" V0Leg/primary/hZY" , " Z vs. Y;Z (cm);Y (cm)" , kTH2F , {{200 , -100 , 100 }, {40 , -20 , 20 }}, false );
254240 fRegistry .add (" V0Leg/primary/hPtGen_DeltaPtOverPtGen" , " electron p_{T} resolution;p_{T}^{gen} (GeV/c);(p_{T}^{rec} - p_{T}^{gen})/p_{T}^{gen}" , kTH2F , {{1000 , 0 , 10 }, {200 , -1 .0f , 1 .0f }}, true );
255241 fRegistry .add (" V0Leg/primary/hPtGen_DeltaEta" , " electron #eta resolution;p_{T}^{gen} (GeV/c);#eta^{rec} - #eta^{gen}" , kTH2F , {{1000 , 0 , 10 }, {100 , -0 .5f , 0 .5f }}, true );
256- if (cfg_fill_resolution) {
257- fRegistry .add (" V0Leg/primary/hsPhotonResolution" ,
258- " Photon resolution;p_{T};#eta;R_{xy};Z_{conv};Z_{vtx};#Deltap_{T}/p_{T};#Delta#eta;#Delta#phi" ,
259- kTHnSparseF ,
260- {{100 , 0 , 10 },
261- {80 , -1.6 , 1.6 },
262- {100 , 0 , 100 },
263- {100 , -50 , 50 },
264- {100 , -50 , 50 },
265- {200 , -1 , 1 },
266- {100 , -0.5 , 0.5 },
267- {100 , -0.5 , 0.5 }},
268- false );
269- }
270-
271242 fRegistry .add (" V0Leg/primary/hPtGen_DeltaPhi" , " electron #varphi resolution;p_{T}^{gen} (GeV/c);#varphi^{rec} - #varphi^{gen} (rad.)" , kTH2F , {{1000 , 0 , 10 }, {100 , -0 .5f , 0 .5f }}, true );
272243 fRegistry .add (" V0Leg/primary/hRxyGen_DeltaPtOverPtGen" , " photon p_{T} resolution; R_{xy}^{gen} (cm);(p_{T}^{rec} - p_{T}^{gen})/p_{T}^{gen}" , kTH2F , {{100 , 0 , 100 }, {200 , -1 .0f , 1 .0f }}, true );
273244 fRegistry .add (" V0Leg/primary/hRxyGen_DeltaEta" , " photon #eta resolution;R_{xy}^{gen} (cm);#eta^{rec} - #eta^{gen}" , kTH2F , {{100 , 0 , 100 }, {100 , -0 .5f , 0 .5f }}, true );
@@ -391,9 +362,6 @@ struct PCMQCMC {
391362 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hMassGamma" ), v0.v0radius (), v0.mGamma ());
392363 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hKFChi2vsM" ), v0.mGamma (), v0.chiSquareNDF ());
393364 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hKFChi2vsR" ), v0.v0radius (), v0.chiSquareNDF ());
394- fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hsConvVtxZPtR" ),
395- v0.vz (), v0.pt (), v0.v0radius ());
396-
397365 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hKFChi2vsX" ), v0.vx (), v0.chiSquareNDF ());
398366 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hKFChi2vsY" ), v0.vy (), v0.chiSquareNDF ());
399367 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hKFChi2vsZ" ), v0.vz (), v0.chiSquareNDF ());
@@ -402,17 +370,6 @@ struct PCMQCMC {
402370 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hPtGen_DeltaPhi" ), mcphoton.pt (), v0.phi () - mcphoton.phi ());
403371 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hRxyGen_DeltaPtOverPtGen" ), std::sqrt (std::pow (mcleg.vx (), 2 ) + std::pow (mcleg.vy (), 2 )), (v0.pt () - mcphoton.pt ()) / mcphoton.pt ());
404372 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hRxyGen_DeltaEta" ), std::sqrt (std::pow (mcleg.vx (), 2 ) + std::pow (mcleg.vy (), 2 )), v0.eta () - mcphoton.eta ());
405- if (cfg_fill_resolution) {
406- fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hsPhotonResolution" ),
407- mcphoton.pt (),
408- mcphoton.eta (),
409- std::sqrt (mcleg.vx () * mcleg.vx () + mcleg.vy () * mcleg.vy ()),
410- mcleg.vz (),
411- v0.vz (),
412- (v0.pt () - mcphoton.pt ()) / mcphoton.pt (),
413- v0.eta () - mcphoton.eta (),
414- TVector2::Phi_mpi_pi (v0.phi () - mcphoton.phi ()));
415- }
416373 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hRxyGen_DeltaPhi" ), std::sqrt (std::pow (mcleg.vx (), 2 ) + std::pow (mcleg.vy (), 2 )), v0.phi () - mcphoton.phi ());
417374 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hRxyGen_DeltaR" ), std::sqrt (std::pow (mcleg.vx (), 2 ) + std::pow (mcleg.vy (), 2 )), v0.v0radius () - std::sqrt (std::pow (mcleg.vx (), 2 ) + std::pow (mcleg.vy (), 2 )));
418375 fRegistry .fill (HIST (" V0/" ) + HIST (mcphoton_types[mctype]) + HIST (" hConvPoint_diffX" ), mcleg.vx (), v0.vx () - mcleg.vx ());
@@ -457,17 +414,6 @@ struct PCMQCMC {
457414 fRegistry .fill (HIST (" V0Leg/" ) + HIST (mcphoton_types[mctype]) + HIST (" hPtGen_DeltaEta" ), mcleg.pt (), leg.eta () - mcleg.eta ());
458415 fRegistry .fill (HIST (" V0Leg/" ) + HIST (mcphoton_types[mctype]) + HIST (" hPtGen_DeltaPhi" ), mcleg.pt (), leg.phi () - mcleg.phi ());
459416 fRegistry .fill (HIST (" V0Leg/" ) + HIST (mcphoton_types[mctype]) + HIST (" hRxyGen_DeltaPtOverPtGen" ), std::sqrt (std::pow (mcleg.vx (), 2 ) + std::pow (mcleg.vy (), 2 )), (leg.pt () - mcleg.pt ()) / mcleg.pt ());
460- if (cfg_fill_resolution) {
461- fRegistry .fill (HIST (" V0Leg/" ) + HIST (mcphoton_types[mctype]) + HIST (" hsPhotonResolution" ),
462- mcleg.pt (),
463- mcleg.eta (),
464- std::sqrt (mcleg.vx () * mcleg.vx () + mcleg.vy () * mcleg.vy ()),
465- mcleg.vz (),
466- leg.z (),
467- (leg.pt () - mcleg.pt ()) / mcleg.pt (),
468- leg.eta () - mcleg.eta (),
469- TVector2::Phi_mpi_pi (leg.phi () - mcleg.phi ()));
470- }
471417 fRegistry .fill (HIST (" V0Leg/" ) + HIST (mcphoton_types[mctype]) + HIST (" hRxyGen_DeltaEta" ), std::sqrt (std::pow (mcleg.vx (), 2 ) + std::pow (mcleg.vy (), 2 )), leg.eta () - mcleg.eta ());
472418 fRegistry .fill (HIST (" V0Leg/" ) + HIST (mcphoton_types[mctype]) + HIST (" hRxyGen_DeltaPhi" ), std::sqrt (std::pow (mcleg.vx (), 2 ) + std::pow (mcleg.vy (), 2 )), leg.phi () - mcleg.phi ());
473419 }
0 commit comments