@@ -264,26 +264,26 @@ struct SGSpectraAnalyzer {
264264 }
265265 }
266266 }
267- }
268267 }
269- if (goodtracks > 1 ) {
270- float W_gPb = TMath::Sqrt ( 2 * 2680 * sum. M () * TMath::Exp (sum. Rapidity ()));
271- if ( sum.M () < . 2 )
272- std::cout << goodtracks << " \t " << sum.M () << " \t " << sum. Pt () << std::endl;
273- if (gapSide < 2 ) {
274- registry. get <TH2>( HIST ( " E_mult_SG " ))-> Fill (goodtracks, W_gPb);
275- registry.get <TH2>(HIST (" E_M_SG " ))->Fill (sum. M () , W_gPb);
276- registry.get <TH2>(HIST (" E_Y_SG " ))->Fill (sum.Rapidity (), W_gPb);
277- registry.get <TH1 >(HIST (" all_tracks_SG " ))->Fill (alltracks );
278- registry.get <TH1>(HIST (" good_tracks_SG " ))->Fill (goodtracks );
279- } else {
280- registry. get <TH2>( HIST ( " E_mult_DG " ))-> Fill (goodtracks, W_gPb);
281- registry.get <TH2>(HIST (" E_M_DG " ))->Fill (sum. M () , W_gPb);
282- registry.get <TH2>(HIST (" E_Y_DG " ))->Fill (sum.Rapidity (), W_gPb);
283- registry.get <TH1 >(HIST (" all_tracks_DG " ))->Fill (alltracks );
284- registry.get <TH1>(HIST (" good_tracks_DG " ))->Fill (goodtracks );
285- }
268+ }
269+ if (goodtracks > 1 ) {
270+ float W_gPb = TMath::Sqrt ( 2 * 2680 * sum.M () * TMath::Exp (sum. Rapidity ()));
271+ if ( sum.M () < . 2 )
272+ std::cout << goodtracks << " \t " << sum. M () << " \t " << sum. Pt () << std::endl;
273+ if (gapSide < 2 ) {
274+ registry.get <TH2>(HIST (" E_mult_SG " ))->Fill (goodtracks , W_gPb);
275+ registry.get <TH2>(HIST (" E_M_SG " ))->Fill (sum.M (), W_gPb);
276+ registry.get <TH2 >(HIST (" E_Y_SG " ))->Fill (sum. Rapidity (), W_gPb );
277+ registry.get <TH1>(HIST (" all_tracks_SG " ))->Fill (alltracks );
278+ registry. get <TH1>( HIST ( " good_tracks_SG " ))-> Fill (goodtracks);
279+ } else {
280+ registry.get <TH2>(HIST (" E_mult_DG " ))->Fill (goodtracks , W_gPb);
281+ registry.get <TH2>(HIST (" E_M_DG " ))->Fill (sum.M (), W_gPb);
282+ registry.get <TH2 >(HIST (" E_Y_DG " ))->Fill (sum. Rapidity (), W_gPb );
283+ registry.get <TH1>(HIST (" all_tracks_DG " ))->Fill (alltracks );
284+ registry. get <TH1>( HIST ( " good_tracks_DG " ))-> Fill (goodtracks);
286285 }
286+ }
287287 }
288288
289289 void fillHistograms (const std::string& particleType, float pt, float eta, int gapSide)
0 commit comments