[u/mrichter/AliRoot.git] / PWG4 / AliNeutralMesonSelection.cxx

 /**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
/* $Id$ */

/* History of cvs commits:
 *
 * $Log$
 * Revision 1.3  2007/10/29 13:48:42  gustavo
 * Corrected coding violations
 *
 * Revision 1.2  2007/08/17 12:40:04  schutz
 * New analysis classes by Gustavo Conesa
 *
 * Revision 1.1.2.1  2007/07/26 10:32:09  schutz
 * new analysis classes in the the new analysis framework
 *
 *
 */

//_________________________________________________________________________
// Class that contains methods to select candidate pairs to neutral meson 
// 2 main selections, invariant mass around pi0 (also any other mass),
// apperture angle to distinguish from combinatorial.
//-- Author: Gustavo Conesa (INFN-LNF)

// --- ROOT system ---
#include <TLorentzVector.h>
#include <TH2.h>
#include <TList.h>
#include <TArrayD.h>
 
//---- AliRoot system ----
#include "AliNeutralMesonSelection.h" 
#include "AliLog.h"

ClassImp(AliNeutralMesonSelection)
  
  
//____________________________________________________________________________
  AliNeutralMesonSelection::AliNeutralMesonSelection() : 
    TObject(), fM(0),
    fInvMassMaxCut(0.), fInvMassMinCut(0.),
    fAngleMaxParam(),  fKeepNeutralMesonHistos(0), 
    fhAnglePairNoCut(0), fhAnglePairOpeningAngleCut(0), 
    fhAnglePairAllCut(0), 
    fhInvMassPairNoCut(0), fhInvMassPairOpeningAngleCut(0), 
    fhInvMassPairAllCut(0) 
{
  //Default Ctor
  
  //Initialize parameters
  
  // kGammaHadron and kGammaJet 
  fAngleMaxParam.Set(4) ;
  fAngleMaxParam.Reset(0.);
  
  //Initialize parameters
  InitParameters();
}

//____________________________________________________________________________
AliNeutralMesonSelection::AliNeutralMesonSelection(const AliNeutralMesonSelection & g) :   
  TObject(), fM(g.fM),
  fInvMassMaxCut(g.fInvMassMaxCut), fInvMassMinCut(g.fInvMassMinCut),
  fAngleMaxParam(g.fAngleMaxParam),
  fKeepNeutralMesonHistos(g.fKeepNeutralMesonHistos),
  fhAnglePairNoCut(g. fhAnglePairNoCut), 
  fhAnglePairOpeningAngleCut(g. fhAnglePairOpeningAngleCut), 
  fhAnglePairAllCut(g. fhAnglePairAllCut), 
  fhInvMassPairNoCut(g.fhInvMassPairNoCut),  
  fhInvMassPairOpeningAngleCut(g.fhInvMassPairOpeningAngleCut), 
  fhInvMassPairAllCut(g.fhInvMassPairAllCut)
{
  // cpy ctor
}

//_________________________________________________________________________
AliNeutralMesonSelection & AliNeutralMesonSelection::operator = (const AliNeutralMesonSelection & source)
{
  // assignment operator
  
  if(this == &source)return *this;
  ((TObject *)this)->operator=(source);

  fM = source.fM ;
  fInvMassMaxCut = source.fInvMassMaxCut ; 
  fInvMassMinCut = source.fInvMassMinCut ;
  fAngleMaxParam = source.fAngleMaxParam ;
  fKeepNeutralMesonHistos = source.fKeepNeutralMesonHistos;
 
  fhAnglePairNoCut = source. fhAnglePairNoCut ; 
  fhAnglePairOpeningAngleCut = source. fhAnglePairOpeningAngleCut ; 
  fhAnglePairAllCut = source. fhAnglePairAllCut ; 
  fhInvMassPairNoCut = source.fhInvMassPairNoCut ; 
  fhInvMassPairOpeningAngleCut = source.fhInvMassPairOpeningAngleCut ; 
  fhInvMassPairAllCut = source.fhInvMassPairAllCut ; 
  
  return *this;
  
}

//____________________________________________________________________________
AliNeutralMesonSelection::~AliNeutralMesonSelection() 
{
  //dtor

  if(!fKeepNeutralMesonHistos){
    //Histograms initialized and filled but not passed to output container
    //delete here, I am not sure this is correct
    
    if(fhAnglePairNoCut) delete fhAnglePairNoCut;
    if(fhAnglePairOpeningAngleCut) delete fhAnglePairOpeningAngleCut; 
    if(fhAnglePairAllCut) delete fhAnglePairAllCut;
    if(fhInvMassPairNoCut) delete fhInvMassPairNoCut;
    if(fhInvMassPairOpeningAngleCut) delete fhInvMassPairOpeningAngleCut;
    if(fhInvMassPairAllCut) delete fhInvMassPairAllCut; 

  }
  
}
//________________________________________________________________________
TList *  AliNeutralMesonSelection::GetCreateOutputObjects()
{  

  // Create histograms to be saved in output file and 
  // store them in outputContainer of the analysis class that calls this class.

  TList * outputContainer = new TList() ; 
  outputContainer->SetName("MesonDecayHistos") ; 
  
  fhAnglePairNoCut  = new TH2F
    ("AnglePairNoCut",
     "Angle between all #gamma pair vs E_{#pi^{0}}",200,0,50,200,0,0.2); 
  fhAnglePairNoCut->SetYTitle("Angle (rad)");
  fhAnglePairNoCut->SetXTitle("E_{ #pi^{0}} (GeV)");
  
  fhAnglePairOpeningAngleCut  = new TH2F
    ("AnglePairOpeningAngleCut",
     "Angle between all #gamma pair (opening angle + azimuth cut) vs E_{#pi^{0}}"
     ,200,0,50,200,0,0.2); 
  fhAnglePairOpeningAngleCut->SetYTitle("Angle (rad)");
  fhAnglePairOpeningAngleCut->SetXTitle("E_{ #pi^{0}} (GeV)");
  
  fhAnglePairAllCut  = new TH2F
    ("AnglePairAllCut",
     "Angle between all #gamma pair (opening angle + inv mass cut+azimuth) vs E_{#pi^{0}}"
     ,200,0,50,200,0,0.2); 
  fhAnglePairAllCut->SetYTitle("Angle (rad)");
  fhAnglePairAllCut->SetXTitle("E_{ #pi^{0}} (GeV)");    
  
  //
  fhInvMassPairNoCut  = new TH2F
    ("InvMassPairNoCut","Invariant Mass of all #gamma pair vs E_{#pi^{0}}",
     120,0,120,360,0,0.5); 
  fhInvMassPairNoCut->SetYTitle("Invariant Mass (GeV/c^{2})");
  fhInvMassPairNoCut->SetXTitle("E_{ #pi^{0}} (GeV)");
  
  fhInvMassPairOpeningAngleCut  = new TH2F
    ("InvMassPairOpeningAngleCut",
     "Invariant Mass of #gamma pair (angle cut) vs E_{#pi^{0}}",
     120,0,120,360,0,0.5); 
  fhInvMassPairOpeningAngleCut->SetYTitle("Invariant Mass (GeV/c^{2})");
  fhInvMassPairOpeningAngleCut->SetXTitle(" E_{#pi^{0}}(GeV)");
  
  fhInvMassPairAllCut  = new TH2F
    ("InvMassPairAllCut",
     "Invariant Mass of #gamma pair (opening angle+invmass cut) vs E_{#pi^{0}}",
     120,0,120,360,0,0.5); 
  fhInvMassPairAllCut->SetYTitle("Invariant Mass (GeV/c^{2})");
  fhInvMassPairAllCut->SetXTitle("E_{#pi^{0}}(GeV)");

  outputContainer->Add(fhAnglePairNoCut) ; 
  outputContainer->Add(fhAnglePairOpeningAngleCut) ;
  outputContainer->Add(fhAnglePairAllCut) ; 
  
  outputContainer->Add(fhInvMassPairNoCut) ; 
  outputContainer->Add(fhInvMassPairOpeningAngleCut) ; 
  outputContainer->Add(fhInvMassPairAllCut) ; 
  
  return outputContainer;
}

 //____________________________________________________________________________
void AliNeutralMesonSelection::InitParameters()
{
 
  //Initialize the parameters of the analysis.
  fKeepNeutralMesonHistos = kFALSE ;

  //-------------kHadron, kJetLeadCone-----------------
  fAngleMaxParam.Set(4) ;
  fAngleMaxParam.AddAt(0.4,0);//={0.4,-0.25,0.025,-2e-4};
  fAngleMaxParam.AddAt(-0.25,1) ;
  fAngleMaxParam.AddAt(0.025,2) ;
  fAngleMaxParam.AddAt(-2e-4,3) ;

  fInvMassMaxCut  = 0.16 ;
  fInvMassMinCut  = 0.11 ;

  fM = 0.1349766;//neutralMeson mass
}

//__________________________________________________________________________-
Bool_t AliNeutralMesonSelection::IsAngleInWindow(const Float_t angle,const Float_t e) const {
  //Check if the opening angle of the candidate pairs is inside 
  //our selection windowd

  Bool_t result = kFALSE;
  Double_t max =  fAngleMaxParam.At(0)*TMath::Exp(fAngleMaxParam.At(1)*e)
    +fAngleMaxParam.At(2)+fAngleMaxParam.At(3)*e;
  Double_t arg = (e*e-2*fM*fM)/(e*e);
  Double_t min = 100. ;
  if(arg>0.)
    min = TMath::ACos(arg);

  if((angle<max)&&(angle>=min))
    result = kTRUE;
 
  return result;
}

//____________________________________________________________________________
Bool_t  AliNeutralMesonSelection::SelectPair(TLorentzVector gammai, TLorentzVector gammaj)  
{  
  
  //Search for the neutral pion within selection cuts
  Bool_t goodpair = kFALSE ;
  
  Double_t pt  = (gammai+gammaj).Pt();
  Double_t phi = (gammai+gammaj).Phi();
  if(phi < 0)
    phi+=TMath::TwoPi();
  Double_t invmass = (gammai+gammaj).M();
  Double_t angle   = gammaj.Angle(gammai.Vect());
  Double_t e       = (gammai+gammaj).E();
  
  //Fill histograms with no cuts applied.
  fhAnglePairNoCut->Fill(e,angle);
  fhInvMassPairNoCut->Fill(e,invmass);
    
  //Cut on the aperture of the pair
  if(IsAngleInWindow(angle,e)){
    fhAnglePairOpeningAngleCut     ->Fill(e,angle);
    fhInvMassPairOpeningAngleCut->Fill(e,invmass);
    AliDebug(2,Form("Angle cut: pt %f, phi %f",pt,phi));
    
    //Cut on the invariant mass of the pair
    if((invmass>fInvMassMinCut) && (invmass<fInvMassMaxCut)){ 
      fhInvMassPairAllCut  ->Fill(e,invmass);
      fhAnglePairAllCut       ->Fill(e,angle);
      goodpair = kTRUE;
      AliDebug(2,Form("IM cut: pt %f, phi %f",pt,phi));
    }//(invmass>0.125) && (invmass<0.145)
  }//Opening angle cut
  
  return goodpair; 
  
}

//__________________________________________________________________
void AliNeutralMesonSelection::Print(const Option_t * opt) const
{

  //Print some relevant parameters set for the analysis
  if(! opt)
    return;
  
  Info("Print", "%s %s", GetName(), GetTitle() ) ;

  printf("mass : %f  \n", fM );
  printf("Invariant mass limits : %f < m < %f \n", fInvMassMinCut , fInvMassMinCut );
  printf("Angle selection param: \n");
  printf("p0 :     %f", fAngleMaxParam.At(0));
  printf("p1 :     %f", fAngleMaxParam.At(1));
  printf("p2 :     %f", fAngleMaxParam.At(2));
  printf("p3 :     %f", fAngleMaxParam.At(3));

  printf("Keep Neutral Meson Histos = %d\n",fKeepNeutralMesonHistos);

}
Commit	Line	Data
bdcfac30	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15	/* $Id$ */
	16
	17	/* History of cvs commits:
	18	*
	19	* $Log$
3788def4	20	* Revision 1.3 2007/10/29 13:48:42 gustavo
	21	* Corrected coding violations
	22	*
3bb2c538	23	* Revision 1.2 2007/08/17 12:40:04 schutz
	24	* New analysis classes by Gustavo Conesa
	25	*
bdcfac30	26	* Revision 1.1.2.1 2007/07/26 10:32:09 schutz
	27	* new analysis classes in the the new analysis framework
	28	*
	29	*
	30	*/
	31
	32	//_________________________________________________________________________
	33	// Class that contains methods to select candidate pairs to neutral meson
3bb2c538	34	// 2 main selections, invariant mass around pi0 (also any other mass),
3bb2c538	35	// apperture angle to distinguish from combinatorial.
bdcfac30	36	//-- Author: Gustavo Conesa (INFN-LNF)
	37
	38	// --- ROOT system ---
bdcfac30	39	#include <TLorentzVector.h>
	40	#include <TH2.h>
	41	#include <TList.h>
3bb2c538	42	#include <TArrayD.h>
3bb2c538	43
bdcfac30	44	//---- AliRoot system ----
bdcfac30	45	#include "AliNeutralMesonSelection.h"
bdcfac30	46	#include "AliLog.h"
	47
	48	ClassImp(AliNeutralMesonSelection)
	49
	50
	51	//____________________________________________________________________________
	52	AliNeutralMesonSelection::AliNeutralMesonSelection() :
2c8efea8	53	TObject(), fM(0),
bdcfac30	54	fInvMassMaxCut(0.), fInvMassMinCut(0.),
2cc4014d	55	fAngleMaxParam(), fKeepNeutralMesonHistos(0),
	56	fhAnglePairNoCut(0), fhAnglePairOpeningAngleCut(0),
	57	fhAnglePairAllCut(0),
2c8efea8	58	fhInvMassPairNoCut(0), fhInvMassPairOpeningAngleCut(0),
2c8efea8	59	fhInvMassPairAllCut(0)
bdcfac30	60	{
	61	//Default Ctor
	62
	63	//Initialize parameters
	64
	65	// kGammaHadron and kGammaJet
	66	fAngleMaxParam.Set(4) ;
	67	fAngleMaxParam.Reset(0.);
	68
	69	//Initialize parameters
	70	InitParameters();
	71	}
	72
	73	//____________________________________________________________________________
	74	AliNeutralMesonSelection::AliNeutralMesonSelection(const AliNeutralMesonSelection & g) :
2c8efea8	75	TObject(), fM(g.fM),
bdcfac30	76	fInvMassMaxCut(g.fInvMassMaxCut), fInvMassMinCut(g.fInvMassMinCut),
2c8efea8	77	fAngleMaxParam(g.fAngleMaxParam),
bdcfac30	78	fKeepNeutralMesonHistos(g.fKeepNeutralMesonHistos),
bdcfac30	79	fhAnglePairNoCut(g. fhAnglePairNoCut),
bdcfac30	80	fhAnglePairOpeningAngleCut(g. fhAnglePairOpeningAngleCut),
	81	fhAnglePairAllCut(g. fhAnglePairAllCut),
	82	fhInvMassPairNoCut(g.fhInvMassPairNoCut),
bdcfac30	83	fhInvMassPairOpeningAngleCut(g.fhInvMassPairOpeningAngleCut),
	84	fhInvMassPairAllCut(g.fhInvMassPairAllCut)
	85	{
	86	// cpy ctor
	87	}
	88
	89	//_________________________________________________________________________
	90	AliNeutralMesonSelection & AliNeutralMesonSelection::operator = (const AliNeutralMesonSelection & source)
	91	{
	92	// assignment operator
	93
	94	if(this == &source)return *this;
	95	((TObject *)this)->operator=(source);
	96
bdcfac30	97	fM = source.fM ;
	98	fInvMassMaxCut = source.fInvMassMaxCut ;
	99	fInvMassMinCut = source.fInvMassMinCut ;
	100	fAngleMaxParam = source.fAngleMaxParam ;
bdcfac30	101	fKeepNeutralMesonHistos = source.fKeepNeutralMesonHistos;
	102
	103	fhAnglePairNoCut = source. fhAnglePairNoCut ;
bdcfac30	104	fhAnglePairOpeningAngleCut = source. fhAnglePairOpeningAngleCut ;
	105	fhAnglePairAllCut = source. fhAnglePairAllCut ;
	106	fhInvMassPairNoCut = source.fhInvMassPairNoCut ;
bdcfac30	107	fhInvMassPairOpeningAngleCut = source.fhInvMassPairOpeningAngleCut ;
	108	fhInvMassPairAllCut = source.fhInvMassPairAllCut ;
	109
	110	return *this;
	111
	112	}
	113
	114	//____________________________________________________________________________
	115	AliNeutralMesonSelection::~AliNeutralMesonSelection()
	116	{
2c8efea8	117	//dtor
bdcfac30	118
2c8efea8	119	if(!fKeepNeutralMesonHistos){
	120	//Histograms initialized and filled but not passed to output container
	121	//delete here, I am not sure this is correct
	122
	123	if(fhAnglePairNoCut) delete fhAnglePairNoCut;
	124	if(fhAnglePairOpeningAngleCut) delete fhAnglePairOpeningAngleCut;
	125	if(fhAnglePairAllCut) delete fhAnglePairAllCut;
	126	if(fhInvMassPairNoCut) delete fhInvMassPairNoCut;
	127	if(fhInvMassPairOpeningAngleCut) delete fhInvMassPairOpeningAngleCut;
	128	if(fhInvMassPairAllCut) delete fhInvMassPairAllCut;
bdcfac30	129
2c8efea8	130	}
	131
	132	}
bdcfac30	133	//________________________________________________________________________
	134	TList * AliNeutralMesonSelection::GetCreateOutputObjects()
	135	{
	136
	137	// Create histograms to be saved in output file and
2c8efea8	138	// store them in outputContainer of the analysis class that calls this class.
2c8efea8	139
bdcfac30	140	TList * outputContainer = new TList() ;
	141	outputContainer->SetName("MesonDecayHistos") ;
	142
	143	fhAnglePairNoCut = new TH2F
	144	("AnglePairNoCut",
	145	"Angle between all #gamma pair vs E_{#pi^{0}}",200,0,50,200,0,0.2);
	146	fhAnglePairNoCut->SetYTitle("Angle (rad)");
	147	fhAnglePairNoCut->SetXTitle("E_{ #pi^{0}} (GeV)");
	148
	149	fhAnglePairOpeningAngleCut = new TH2F
	150	("AnglePairOpeningAngleCut",
	151	"Angle between all #gamma pair (opening angle + azimuth cut) vs E_{#pi^{0}}"
	152	,200,0,50,200,0,0.2);
	153	fhAnglePairOpeningAngleCut->SetYTitle("Angle (rad)");
	154	fhAnglePairOpeningAngleCut->SetXTitle("E_{ #pi^{0}} (GeV)");
	155
	156	fhAnglePairAllCut = new TH2F
	157	("AnglePairAllCut",
	158	"Angle between all #gamma pair (opening angle + inv mass cut+azimuth) vs E_{#pi^{0}}"
	159	,200,0,50,200,0,0.2);
	160	fhAnglePairAllCut->SetYTitle("Angle (rad)");
	161	fhAnglePairAllCut->SetXTitle("E_{ #pi^{0}} (GeV)");
	162
	163	//
	164	fhInvMassPairNoCut = new TH2F
	165	("InvMassPairNoCut","Invariant Mass of all #gamma pair vs E_{#pi^{0}}",
	166	120,0,120,360,0,0.5);
	167	fhInvMassPairNoCut->SetYTitle("Invariant Mass (GeV/c^{2})");
	168	fhInvMassPairNoCut->SetXTitle("E_{ #pi^{0}} (GeV)");
	169
	170	fhInvMassPairOpeningAngleCut = new TH2F
	171	("InvMassPairOpeningAngleCut",
	172	"Invariant Mass of #gamma pair (angle cut) vs E_{#pi^{0}}",
	173	120,0,120,360,0,0.5);
	174	fhInvMassPairOpeningAngleCut->SetYTitle("Invariant Mass (GeV/c^{2})");
	175	fhInvMassPairOpeningAngleCut->SetXTitle(" E_{#pi^{0}}(GeV)");
	176
	177	fhInvMassPairAllCut = new TH2F
	178	("InvMassPairAllCut",
	179	"Invariant Mass of #gamma pair (opening angle+invmass cut) vs E_{#pi^{0}}",
	180	120,0,120,360,0,0.5);
	181	fhInvMassPairAllCut->SetYTitle("Invariant Mass (GeV/c^{2})");
	182	fhInvMassPairAllCut->SetXTitle("E_{#pi^{0}}(GeV)");
	183
bdcfac30	184	outputContainer->Add(fhAnglePairNoCut) ;
	185	outputContainer->Add(fhAnglePairOpeningAngleCut) ;
	186	outputContainer->Add(fhAnglePairAllCut) ;
	187
	188	outputContainer->Add(fhInvMassPairNoCut) ;
	189	outputContainer->Add(fhInvMassPairOpeningAngleCut) ;
	190	outputContainer->Add(fhInvMassPairAllCut) ;
bdcfac30	191
	192	return outputContainer;
	193	}
	194
	195	//____________________________________________________________________________
	196	void AliNeutralMesonSelection::InitParameters()
	197	{
	198
	199	//Initialize the parameters of the analysis.
2c8efea8	200	fKeepNeutralMesonHistos = kFALSE ;
bdcfac30	201
	202	//-------------kHadron, kJetLeadCone-----------------
	203	fAngleMaxParam.Set(4) ;
	204	fAngleMaxParam.AddAt(0.4,0);//={0.4,-0.25,0.025,-2e-4};
	205	fAngleMaxParam.AddAt(-0.25,1) ;
	206	fAngleMaxParam.AddAt(0.025,2) ;
	207	fAngleMaxParam.AddAt(-2e-4,3) ;
	208
	209	fInvMassMaxCut = 0.16 ;
	210	fInvMassMinCut = 0.11 ;
	211
	212	fM = 0.1349766;//neutralMeson mass
bdcfac30	213	}
	214
	215	//__________________________________________________________________________-
3bb2c538	216	Bool_t AliNeutralMesonSelection::IsAngleInWindow(const Float_t angle,const Float_t e) const {
bdcfac30	217	//Check if the opening angle of the candidate pairs is inside
	218	//our selection windowd
	219
	220	Bool_t result = kFALSE;
	221	Double_t max = fAngleMaxParam.At(0)TMath::Exp(fAngleMaxParam.At(1)e)
	222	+fAngleMaxParam.At(2)+fAngleMaxParam.At(3)*e;
	223	Double_t arg = (ee-2fMfM)/(ee);
	224	Double_t min = 100. ;
	225	if(arg>0.)
	226	min = TMath::ACos(arg);
	227
	228	if((angle<max)&&(angle>=min))
	229	result = kTRUE;
	230
	231	return result;
	232	}
	233
	234	//____________________________________________________________________________
2c8efea8	235	Bool_t AliNeutralMesonSelection::SelectPair(TLorentzVector gammai, TLorentzVector gammaj)
bdcfac30	236	{
	237
	238	//Search for the neutral pion within selection cuts
bdcfac30	239	Bool_t goodpair = kFALSE ;
	240
	241	Double_t pt = (gammai+gammaj).Pt();
	242	Double_t phi = (gammai+gammaj).Phi();
	243	if(phi < 0)
	244	phi+=TMath::TwoPi();
	245	Double_t invmass = (gammai+gammaj).M();
	246	Double_t angle = gammaj.Angle(gammai.Vect());
	247	Double_t e = (gammai+gammaj).E();
	248
	249	//Fill histograms with no cuts applied.
	250	fhAnglePairNoCut->Fill(e,angle);
	251	fhInvMassPairNoCut->Fill(e,invmass);
bdcfac30	252
2c8efea8	253	//Cut on the aperture of the pair
	254	if(IsAngleInWindow(angle,e)){
	255	fhAnglePairOpeningAngleCut ->Fill(e,angle);
	256	fhInvMassPairOpeningAngleCut->Fill(e,invmass);
	257	AliDebug(2,Form("Angle cut: pt %f, phi %f",pt,phi));
bdcfac30	258
2c8efea8	259	//Cut on the invariant mass of the pair
	260	if((invmass>fInvMassMinCut) && (invmass<fInvMassMaxCut)){
	261	fhInvMassPairAllCut ->Fill(e,invmass);
	262	fhAnglePairAllCut ->Fill(e,angle);
	263	goodpair = kTRUE;
	264	AliDebug(2,Form("IM cut: pt %f, phi %f",pt,phi));
	265	}//(invmass>0.125) && (invmass<0.145)
	266	}//Opening angle cut
bdcfac30	267
	268	return goodpair;
	269
	270	}
	271
	272	//__________________________________________________________________
	273	void AliNeutralMesonSelection::Print(const Option_t * opt) const
	274	{
	275
	276	//Print some relevant parameters set for the analysis
	277	if(! opt)
	278	return;
	279
	280	Info("Print", "%s %s", GetName(), GetTitle() ) ;
	281
	282	printf("mass : %f \n", fM );
	283	printf("Invariant mass limits : %f < m < %f \n", fInvMassMinCut , fInvMassMinCut );
	284	printf("Angle selection param: \n");
	285	printf("p0 : %f", fAngleMaxParam.At(0));
	286	printf("p1 : %f", fAngleMaxParam.At(1));
	287	printf("p2 : %f", fAngleMaxParam.At(2));
	288	printf("p3 : %f", fAngleMaxParam.At(3));
	289
bdcfac30	290	printf("Keep Neutral Meson Histos = %d\n",fKeepNeutralMesonHistos);
	291
	292	}