Improvements on the selection of fragmentation photons. (G. Conesa)
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 8 Jan 2007 16:13:28 +0000 (16:13 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 8 Jan 2007 16:13:28 +0000 (16:13 +0000)
PYTHIA6/AliGenPythia.cxx
PYTHIA6/AliGenPythia.h

index 3a3fd63..3dedee2 100644 (file)
@@ -111,17 +111,17 @@ AliGenPythia::AliGenPythia():
     fHeader(0),  
     fRL(0),      
     fFileName(0),
-    fBremssInCalo(kFALSE),
+    fFragPhotonInCalo(kFALSE),
     fPi0InCalo(kFALSE) ,
-    fBremssPi0InEMCAL(kFALSE),
-    fBremssPi0InPHOS(kFALSE),
-    fBremssPi0MinPt(0), 
-    fPHOSMinPhi(220.),
-    fPHOSMaxPhi(320.),
-    fPHOSEta(0.12),
-    fEMCALMinPhi(80.),
-    fEMCALMaxPhi(190.),
-    fEMCALEta(0.7)
+    fCheckEMCAL(kFALSE),
+    fCheckPHOS(kFALSE),
+    fFragPhotonOrPi0MinPt(0), 
+    fPHOSMinPhi(219.),
+    fPHOSMaxPhi(321.),
+    fPHOSEta(0.13),
+    fEMCALMinPhi(79.),
+    fEMCALMaxPhi(191.),
+    fEMCALEta(0.71)
 
 {
 // Default Constructor
@@ -195,17 +195,17 @@ AliGenPythia::AliGenPythia(Int_t npart)
      fHeader(0),  
      fRL(0),      
      fFileName(0),
-     fBremssInCalo(kFALSE),
+     fFragPhotonInCalo(kFALSE),
      fPi0InCalo(kFALSE) ,
-     fBremssPi0InEMCAL(kFALSE),
-     fBremssPi0InPHOS(kFALSE),
-     fBremssPi0MinPt(0),
-     fPHOSMinPhi(220.),
-     fPHOSMaxPhi(320.),
-     fPHOSEta(0.12),
-     fEMCALMinPhi(80.),
-     fEMCALMaxPhi(190.),
-     fEMCALEta(0.7)
+     fCheckEMCAL(kFALSE),
+     fCheckPHOS(kFALSE),
+     fFragPhotonOrPi0MinPt(0),
+     fPHOSMinPhi(219.),
+     fPHOSMaxPhi(321.),
+     fPHOSEta(0.13),
+     fEMCALMinPhi(79.),
+     fEMCALMaxPhi(191.),
+     fEMCALEta(0.71)
 {
 // default charm production at 5. 5 TeV
 // semimuonic decay
@@ -286,17 +286,17 @@ AliGenPythia::AliGenPythia(const AliGenPythia & Pythia)
      fHeader(0),  
      fRL(0),      
      fFileName(0),
-     fBremssInCalo(kFALSE),
+     fFragPhotonInCalo(kFALSE),
      fPi0InCalo(kFALSE) ,
-     fBremssPi0InEMCAL(kFALSE),
-     fBremssPi0InPHOS(kFALSE),
-     fBremssPi0MinPt(0),
-     fPHOSMinPhi(220.),
-     fPHOSMaxPhi(320.),
-     fPHOSEta(0.12),
-     fEMCALMinPhi(80.),
-     fEMCALMaxPhi(190.),
-     fEMCALEta(0.7)   
+     fCheckEMCAL(kFALSE),
+     fCheckPHOS(kFALSE),
+     fFragPhotonOrPi0MinPt(0),
+     fPHOSMinPhi(219.),
+     fPHOSMaxPhi(321.),
+     fPHOSEta(0.13),
+     fEMCALMinPhi(79.),
+     fEMCALMaxPhi(191.),
+     fEMCALEta(0.71)   
 {
 // copy constructor
     Pythia.Copy(*this);
@@ -908,24 +908,28 @@ Int_t  AliGenPythia::GenerateMB()
        }
     }
 
-    // Select jets with bremsstrahlung or pi0 going to PHOS or EMCAL
-    if (fProcess == kPyJets && (fBremssInCalo || fPi0InCalo) ) {
+    // Select jets with fragmentation photon or pi0 going to PHOS or EMCAL
+    if (fProcess == kPyJets && (fFragPhotonInCalo || fPi0InCalo) ) {
 
       Bool_t ok = kFALSE;
 
       Int_t pdg  = 0; 
-      if (fBremssInCalo) pdg = 22   ; // Photon
+      if (fFragPhotonInCalo) pdg = 22   ; // Photon
       else if (fPi0InCalo) pdg = 111 ; // Pi0
 
       for (i=0; i< np; i++) {
        TParticle* iparticle = (TParticle *) fParticles->At(i);
        if(iparticle->GetStatusCode()==1 && iparticle->GetPdgCode()==pdg && 
-          iparticle->Pt() > fBremssPi0MinPt){
-         Float_t phi = iparticle->Phi()*180./TMath::Pi(); //Convert to degrees
-         Float_t eta =TMath::Abs(iparticle->Eta());//in calos etamin=-etamax     
-         if((fBremssPi0InEMCAL && IsInEMCAL(phi,eta)) ||
-             fBremssPi0InPHOS    && IsInPHOS(phi,eta)) 
-           ok =kTRUE;  
+          iparticle->Pt() > fFragPhotonOrPi0MinPt){
+         if(pdg == 111 || 
+            (pdg == 22 && iparticle->GetFirstMother() != 11))//No photon from hadron decay
+           {
+             Float_t phi = iparticle->Phi()*180./TMath::Pi(); //Convert to degrees
+             Float_t eta =TMath::Abs(iparticle->Eta());//in calos etamin=-etamax         
+             if((fCheckEMCAL && IsInEMCAL(phi,eta)) ||
+                (fCheckPHOS    && IsInPHOS(phi,eta)) )
+               ok =kTRUE;
+           }
        }
       }
       if(!ok)
@@ -1421,10 +1425,9 @@ void AliGenPythia::GetSubEventTime()
 Bool_t AliGenPythia::IsInEMCAL(Float_t phi, Float_t eta)
 {
   // Is particle in EMCAL acceptance? 
-  // Acceptance slightly larger considered.
   // phi in degrees, etamin=-etamax
-  if(phi > fEMCALMinPhi-0.1  && phi < fEMCALMaxPhi+0.1 && 
-     eta < fEMCALEta+0.01   ) 
+  if(phi > fEMCALMinPhi  && phi < fEMCALMaxPhi && 
+     eta < fEMCALEta  ) 
     return kTRUE;
   else 
     return kFALSE;
@@ -1435,8 +1438,8 @@ Bool_t AliGenPythia::IsInPHOS(Float_t phi, Float_t eta)
   // Is particle in PHOS acceptance? 
   // Acceptance slightly larger considered.
   // phi in degrees, etamin=-etamax
-  if(phi > fPHOSMinPhi-0.1  && phi < fPHOSMaxPhi+0.1 && 
-     eta < fPHOSEta+0.01   ) 
+  if(phi > fPHOSMinPhi  && phi < fPHOSMaxPhi && 
+     eta < fPHOSEta  ) 
     return kTRUE;
   else 
     return kFALSE;
index 48eb0a2..d650974 100644 (file)
@@ -77,12 +77,23 @@ class AliGenPythia : public AliGenMC
     virtual void    SetGammaPhiRange(Float_t phimin = 0., Float_t phimax = 360.)
        {fPhiMinGamma = TMath::Pi()*phimin/180.; fPhiMaxGamma = TMath::Pi()*phimax/180.;}
 
-   // Select jets with bremsstrahlung or pi0 going to PHOS or EMCAL
-    virtual void  SetBremssInCalo(Bool_t b)           {fBremssInCalo = b;}
-    virtual void  SetPi0InCalo       (Bool_t b)           {fPi0InCalo = b;}
-    virtual void  SetBremssPi0InPHOS(Bool_t b)    {fBremssPi0InPHOS = b;}
-    virtual void  SetBremssPi0InEMCAL(Bool_t b)  {fBremssPi0InEMCAL = b;}
-    virtual void  SetBremsPi0MinPt(Float_t pt)      {fBremssPi0MinPt=pt;}
+   // Select jets with fragmentation photon or pi0 going to PHOS or EMCAL
+
+    virtual void  SetFragPhotonInCalo(Bool_t b) {fFragPhotonInCalo = b;}
+    virtual void  SetPi0InCalo       (Bool_t b)      {fPi0InCalo = b;}
+    virtual void  SetCheckPHOS (Bool_t b)        {fCheckPHOS = b;}
+    virtual void  SetCheckEMCAL(Bool_t b)       {fCheckEMCAL = b;}
+
+    virtual void  SetFragPhotonInEMCAL(Bool_t b) 
+    {fCheckEMCAL = b;  fFragPhotonInCalo = b ;}
+    virtual void  SetFragPhotonInPHOS(Bool_t b) 
+    {fCheckPHOS = b ;   fFragPhotonInCalo = b ;}
+    virtual void  SetPi0InEMCAL(Bool_t b) 
+    {fCheckEMCAL = b ;  fPi0InCalo = b ;}
+    virtual void  SetPi0InPHOS(Bool_t b) 
+    {fCheckPHOS = b ;   fPi0InCalo = b ;}
+    
+    virtual void  SetFragPhotonOrPi0MinPt(Float_t pt)      {fFragPhotonOrPi0MinPt=pt;}
   
     // Trigger on a particle
     virtual void    SetTriggerParticle(Int_t particle = 0, Float_t etamax = 0.9) 
@@ -235,11 +246,11 @@ class AliGenPythia : public AliGenMC
     //                                     are counted
     //
 
-    Bool_t fBremssInCalo ; // Option ask for Bremsstrahlung in calorimeters acceptance
+    Bool_t fFragPhotonInCalo ; // Option ask for Fragmentation Photon in calorimeters acceptance
     Bool_t fPi0InCalo ; // Option ask for Pi0 in calorimeters acceptance
-    Bool_t fBremssPi0InEMCAL ; //  Option ask for Bremss or Pi0 in calorimeters EMCAL acceptance
-    Bool_t fBremssPi0InPHOS ; //  Option ask for Bremss or Pi0 in calorimeters PHOS acceptance
-    Float_t fBremssPi0MinPt ; // Minimum momentum of Bremsstrahlung or Pi0 
+    Bool_t fCheckEMCAL ; //  Option ask for FragPhoton or Pi0 in calorimeters EMCAL acceptance
+    Bool_t fCheckPHOS ; //  Option ask for FragPhoton or Pi0 in calorimeters PHOS acceptance
+    Float_t fFragPhotonOrPi0MinPt ; // Minimum momentum of Fragmentation Photon or Pi0 
     //Calorimeters eta-phi acceptance 
     Float_t fPHOSMinPhi ;
     Float_t fPHOSMaxPhi ;