@@ -25,17 +25,6 @@ class UserHooks_promptgamma : public Pythia8::UserHooks
2525 bool doVetoPartonLevel (const Pythia8 ::Event & event ) override {
2626
2727// printf("Event, size %d\n", event.size());
28- //
29- // for (Int_t ida = 0; ida < 10; ida++) {
30- // printf("parton %d, PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n", ida,
31- // event[ida].id(),
32- // event[ida].status(),
33- // event[ida].mother1(),
34- // event[ida].e(),
35- // event[ida].pT(),
36- // event[ida].eta(),
37- // event[ida].phi()*TMath::RadToDeg());
38- // }
3928
4029 // Get the outgoing 2->2 partons.
4130 // The photon and the associated outgoing parton are in position 5 or 6.
@@ -51,6 +40,18 @@ class UserHooks_promptgamma : public Pythia8::UserHooks
5140 if ( event [idGam ].id () != 22 )
5241 {
5342 printf ("No direct photon found in the parton list!\n" );
43+
44+ for (Int_t ida = 0 ; ida < 10 ; ida ++ ) {
45+ printf ("parton %d, PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n" , ida ,
46+ event [ida ].id (),
47+ event [ida ].status (),
48+ event [ida ].mother1 (),
49+ event [ida ].e (),
50+ event [ida ].pT (),
51+ event [ida ].eta (),
52+ event [ida ].phi ()* TMath ::RadToDeg ());
53+ }
54+
5455 return true;
5556 }
5657
@@ -69,18 +70,18 @@ class UserHooks_promptgamma : public Pythia8::UserHooks
6970 //
7071 if ( detector_acceptance (mAcceptance , event [idGam ].phi (), event [idGam ].eta ()) )
7172 {
72- printf ("+++ Accepted event +++ \n" );
73- printf ("gamma, PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n" ,
73+ // printf("+++ Accepted event +++ \n");
74+ printf ("Selected gamma, id %d, PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n" , idGam ,
7475 event [idGam ].id () , event [idGam ].status (), event [idGam ].mother1 (),
7576 event [idGam ].e () , event [idGam ].pT (),
7677 event [idGam ].eta (), event [idGam ].phi ()* TMath ::RadToDeg ());
7778
78- // printf("parton, PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n",
79+ // printf("Back-to-back parton, id %d, PDG %d, status %d, mother %d, E %2.2f, pT %2.2f, eta %2.2f, phi %2.2f\n", idPar,
7980// event[idPar].id() , event[idPar].status(), event[idPar].mother1(),
8081// event[idPar].e() , event[idPar].pT(),
8182// event[idPar].eta(), event[idPar].phi()*TMath::RadToDeg());
82-
83- // Check difference in pT and azimuthal angle, it should be 0 and +-180 degrees, respectively.
83+ //
84+ // // Check difference in pT and azimuthal angle, it should be 0 and +-180 degrees, respectively.
8485// printf("parton-photon, Delta E %2.2f, Delta pT %2.2f, Delta eta %2.2f, Delta phi %2.2f\n",
8586// event[idPar].e() - event[idGam].e(),
8687// event[idPar].pT() - event[idGam].pT(),
0 commit comments