[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemtoUser / AliFemtoCorrFctn3DSphericalEMCIC.cxx

/***************************************************************************
 *
 * $Id: AliFemtoCorrFctn3DSphercicalEMCIC.h  $
 *
 * Author: Nicolas Bock, Ohio State University, bock@mps.ohio-state.edu
 ***************************************************************************
 *
 * Description: Calculates of the 3D Correlation Function in Spherical
 *              coordinates, and also produces histograms to calculate 
 *              Energy Momentum Conservation Induced Correlations  (EMCICs)
 *
 * This Class produces the following histograms as function of Q, theta, phi
 * (for both real and mixed pairs):
 *        1)   E1 + E2
 *        2)   E1 * E2
 *        3)   Pt1*Pt2
 *        4)   Pz1*Pz2
 *  
 * This class is similar to AliFemtoCorrFctn3DSpherical, but it uses Q
 * instead of K to do the binning. 
 * 
 * NOTE: The EMCIC histograms are not averaged in this class, to obtain 
 * the average, the user needs to divide the real pair histograms by 
 * the numerator, and the mixed pair histograms by the denominator.
 *
 ***************************************************************************
 *
 **************************************************************************/

#include "AliFemtoCorrFctn3DSphericalEMCIC.h"
#include "AliFemtoKTPairCut.h"
#include <TMath.h>
#include <TVector2.h>
#include <cstdio>

#ifdef __ROOT__ 
ClassImp(AliFemtoCorrFctn3DSphericalEMCIC)
#endif


//____________________________
AliFemtoCorrFctn3DSphericalEMCIC::AliFemtoCorrFctn3DSphericalEMCIC(char* title, const int& nqbins, const float& QLo, const float& QHi, const int& nphibins, const int& ncthetabins):
  AliFemtoCorrFctn(),
  fNumerator(0),
  fDenominator(0),
/*fEnergyTotalReal(0),  
  fEnergyMultReal(0),        
  fPzMultReal(0),      
  fPtMultReal(0),*/            
  fEnergyTotalMix (0),      
  fEnergyMultMix (0),      
  fPzMultMix(0),            
  fPtMultMix(0),
  fPairCut(0x0)
{
  
  // To better sample on phi shift bin low edge by binsize/2 = Pi/numBins.
  Double_t shiftPhi=TMath::Pi()/nphibins;
  // set up numerator
  char tTitNum[101] = "Num";
  strncat(tTitNum,title, 100);
  fNumerator = new TH3D(tTitNum,title,nqbins,QLo,QHi,nphibins,
			-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
  // set up denominator
  char tTitDen[101] = "Den";
  strncat(tTitDen,title, 100);
  fDenominator = new TH3D(tTitDen,title,nqbins,QLo,QHi,nphibins,
			  -TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);

  //Added histograms to calculate EMCICs , Nicolas Bock 19.01.2010
  //Setup EnergyTotalReal
  /*char tTitNum1[101] = "ESumReal";
  strncat(tTitNum1,title, 100);
  fEnergyTotalReal = new TH3D(tTitNum1,title,nqbins,QLo,QHi,nphibins,
			      -TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
 
  //Setup EnergyMultReal
  char tTitNum2[101] = "EMultReal";
  strncat(tTitNum2,title, 100);
  fEnergyMultReal = new TH3D(tTitNum2,title,nqbins,QLo,QHi,nphibins,
			     -TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
  
  //Setup Pz MultReal
  char tTitNum3[101] = "PzMultReal";
  strncat(tTitNum3,title, 100);
  fPzMultReal = new TH3D(tTitNum3,title,nqbins,QLo,QHi,nphibins,
			 -TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);

  //Setup Pt MultReal
  char tTitNum4[101] = "PtMultReal";
  strncat(tTitNum4,title, 100);
  fPtMultReal = new TH3D(tTitNum4,title,nqbins,QLo,QHi,nphibins,
			 -TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
  */ 


  //Setup EnergyTotalMix
  char tTitNum5[101] = "ESumMix";
  strncat(tTitNum5,title, 100);
  fEnergyTotalMix = new TH3D(tTitNum5,title,nqbins,QLo,QHi,nphibins,
			     -TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
  
  //Setup EnergyMultMix
  char tTitNum6[101] = "EMultMix";
  strncat(tTitNum6,title, 100);
  fEnergyMultMix = new TH3D(tTitNum6,title,nqbins,QLo,QHi,nphibins,
			    -TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
  
  //Setup Pz MultMix
  char tTitNum7[101] = "PzMultMix";
  strncat(tTitNum7,title, 100);
  fPzMultMix = new TH3D(tTitNum7,title,nqbins,QLo,QHi,nphibins,
			-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);

  //Setup Pt MultMix
  char tTitNum8[101] = "PtMultMix";
  strncat(tTitNum8,title, 100);
  fPtMultMix = new TH3D(tTitNum8,title,nqbins,QLo,QHi,nphibins,
			-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,
			ncthetabins,-1.0,1.0);

  // to enable error bar calculation...
  fNumerator->Sumw2();
  fDenominator->Sumw2();
  /*fEnergyTotalReal->Sumw2();
  fEnergyMultReal->Sumw2();
  fPzMultReal->Sumw2();
  fPtMultReal->Sumw2();  */
  fEnergyTotalMix->Sumw2();
  fEnergyMultMix->Sumw2();
  fPzMultMix->Sumw2();
  fPtMultMix->Sumw2();


}

AliFemtoCorrFctn3DSphericalEMCIC::AliFemtoCorrFctn3DSphericalEMCIC(const AliFemtoCorrFctn3DSphericalEMCIC& aCorrFctn) :
  AliFemtoCorrFctn(),
  fNumerator(0),
  fDenominator(0),
  /*fEnergyTotalReal(0),
  fEnergyMultReal(0),        
  fPzMultReal(0),      
  fPtMultReal(0),            */
  fEnergyTotalMix (0),      
  fEnergyMultMix (0),      
  fPzMultMix(0),            
  fPtMultMix(0),
  fPairCut(0x0)
{
  // Copy constructor
  fNumerator = new TH3D(*aCorrFctn.fNumerator);
  fDenominator = new TH3D(*aCorrFctn.fDenominator);
  /*fEnergyTotalReal = new TH3D(*aCorrFctn.fEnergyTotalReal);
  fEnergyMultReal = new TH3D(*aCorrFctn.fEnergyMultReal);
  fPzMultReal = new TH3D(*aCorrFctn.fPzMultReal);
  fPtMultReal = new TH3D(*aCorrFctn.fPtMultReal);*/
  fEnergyTotalMix = new TH3D(*aCorrFctn.fEnergyTotalMix);
  fEnergyMultMix = new TH3D(*aCorrFctn.fEnergyMultMix);
  fPzMultMix = new TH3D(*aCorrFctn.fPzMultMix);
  fPtMultMix = new TH3D(*aCorrFctn.fPtMultMix);
  fPairCut = aCorrFctn.fPairCut;
}
//____________________________
AliFemtoCorrFctn3DSphericalEMCIC::~AliFemtoCorrFctn3DSphericalEMCIC(){
  // Destructor
  delete fNumerator;
  delete fDenominator;
  /*delete fEnergyTotalReal;
  delete fEnergyMultReal;        
  delete fPzMultReal;     
  delete fPtMultReal;            */
  delete fEnergyTotalMix;      
  delete fEnergyMultMix; 
  delete fPzMultMix;   
  delete fPtMultMix;
}
//_________________________
AliFemtoCorrFctn3DSphericalEMCIC& AliFemtoCorrFctn3DSphericalEMCIC::operator=(const AliFemtoCorrFctn3DSphericalEMCIC& aCorrFctn)
{
  // assignment operator
  if (this == &aCorrFctn)
    return *this;

  if (fNumerator) delete fNumerator;
  fNumerator = new TH3D(*aCorrFctn.fNumerator);
  if (fDenominator) delete fDenominator;
  fDenominator = new TH3D(*aCorrFctn.fDenominator);
  /*  if (fEnergyTotalReal) delete fEnergyTotalReal;
  fEnergyTotalReal = new TH3D(*aCorrFctn.fEnergyTotalReal);
  if (fEnergyMultReal) delete fEnergyMultReal;
  fEnergyMultReal = new TH3D(*aCorrFctn.fEnergyMultReal);
  if (fPzMultReal) delete fPzMultReal;
  fPzMultReal = new TH3D(*aCorrFctn.fPzMultReal);
  if (fPtMultReal) delete fPtMultReal;
  fPtMultReal = new TH3D(*aCorrFctn.fPtMultReal);  */
  if (fEnergyTotalMix) delete fEnergyTotalMix;
  fEnergyTotalMix = new TH3D(*aCorrFctn.fEnergyTotalMix);
  if (fEnergyMultMix) delete fEnergyMultMix;
  fEnergyMultMix = new TH3D(*aCorrFctn.fEnergyMultMix);
  if (fPzMultMix) delete fPzMultMix;
  fPzMultMix = new TH3D(*aCorrFctn.fPzMultMix);
  if (fPtMultMix) delete fPtMultMix;
  fPtMultMix = new TH3D(*aCorrFctn.fPtMultMix);
  fPairCut = aCorrFctn.fPairCut;
  
  return *this;
}

//_________________________
void AliFemtoCorrFctn3DSphericalEMCIC::WriteOutHistos(){
  // Write out all histograms to file
  fNumerator->Write();
  fDenominator->Write();
  /*fEnergyTotalReal->Write();
  fEnergyMultReal->Write();        
  fPzMultReal->Write();      
  fPtMultReal->Write();            */
  fEnergyTotalMix->Write();      
  fEnergyMultMix->Write();      
  fPzMultMix->Write();            
  fPtMultMix->Write();
}
//______________________________
TList* AliFemtoCorrFctn3DSphericalEMCIC::GetOutputList()
{
  // Prepare the list of objects to be written to the output
  TList *tOutputList = new TList();

  tOutputList->Add(fNumerator); 
  tOutputList->Add(fDenominator);  
  /*  tOutputList->Add(fEnergyTotalReal);
  tOutputList->Add(fEnergyMultReal);        
  tOutputList->Add(fPzMultReal);      
  tOutputList->Add(fPtMultReal);            */
  tOutputList->Add(fEnergyTotalMix );      
  tOutputList->Add(fEnergyMultMix );      
  tOutputList->Add(fPzMultMix);            
  tOutputList->Add(fPtMultMix);
  return tOutputList;
}

//_________________________
void AliFemtoCorrFctn3DSphericalEMCIC::Finish(){
  // here is where we should normalize, fit, etc...
}

//____________________________
AliFemtoString AliFemtoCorrFctn3DSphericalEMCIC::Report(){
  // Construct the report
  string stemp = "PRF Frame SphericalEMCIC 3D Correlation Function Report:\n";
  char ctemp[100];
  snprintf(ctemp , 100, "Number of entries in numerator:\t%E\n",fNumerator->GetEntries());
  stemp += ctemp;
  snprintf(ctemp , 100, "Number of entries in denominator:\t%E\n",fDenominator->GetEntries());
  stemp += ctemp;

  if (fPairCut){
    snprintf(ctemp , 100, "Here is the PairCut specific to this CorrFctn\n");
    stemp += ctemp;
    stemp += fPairCut->Report();
  }
  else{
    snprintf(ctemp , 100, "No PairCut specific to this CorrFctn\n");
    stemp += ctemp;
  }

  //  
  AliFemtoString returnThis = stemp;
  return returnThis;
}
//____________________________
void AliFemtoCorrFctn3DSphericalEMCIC::AddRealPair( AliFemtoPair* pair){
  // perform operations on real pairs
  if (fPairCut){
    AliFemtoKTPairCut *ktc = dynamic_cast<AliFemtoKTPairCut *>(fPairCut);
    if (!ktc){
      if (!(fPairCut->Pass(pair))) return;
    }
    else
      if (!(ktc->Pass(pair))) return;
  }

  //                          
  double tQO = pair->QOutCMS();  
  double tQS = pair->QSideCMS();  
  double tQL = pair->QLongCMS();  

  double tQR = sqrt(tQO*tQO + tQS*tQS + tQL*tQL);
  double tQC = 0;
  if ( fabs(tQR) < 1e-10 ) tQC = 0.0;
  else tQC = tQL/tQR;
  double tQP = atan2(tQS,tQO);

  fNumerator->Fill(tQR,tQP,tQC);
  
  // EMCICs  
  /*AliFemtoLorentzVector tMom1 = pair->Track1()->FourMomentum();
  AliFemtoLorentzVector tMom2 = pair->Track2()->FourMomentum();
  double tE1 = tMom1.e();
  double tE2 = tMom2.e();
  double tPz1 = tMom1.pz();
  double tPz2 = tMom2.pz();
  
  TVector2 tPt1;  
  TVector2 tPt2; 
  tPt1.Set(tMom1.px(),tMom1.py());
  tPt2.Set(tMom2.px(),tMom2.py());
  double tPt1DotPt2 = tPt1*tPt2;
  
  fEnergyTotalReal->Fill(tQR,tQP,tQC,tE1+tE2);
  fEnergyMultReal->Fill(tQR,tQP,tQC,tE1*tE2);
  fPzMultReal->Fill(tQR,tQP,tQC,tPz1*tPz2);
  fPtMultReal->Fill(tQR,tQP,tQC,tPt1DotPt2);*/
   

}
//____________________________
void AliFemtoCorrFctn3DSphericalEMCIC::AddMixedPair( AliFemtoPair* pair){
  // perform operations on mixed pairs
  if (fPairCut){
    AliFemtoKTPairCut *ktc = dynamic_cast<AliFemtoKTPairCut *>(fPairCut);
    if (!ktc){
      if (!(fPairCut->Pass(pair))) return;
    }
    else
      if (!(ktc->Pass(pair))) return;
  }
  

 //                          //Changed K to Q to be in LCMS, N. Bock
  double tQO = pair->QOutCMS();  
  double tQS = pair->QSideCMS();   
  double tQL = pair->QLongCMS();  

  double tQR = sqrt(tQO*tQO + tQS*tQS + tQL*tQL);
  double tQC;
  if ( fabs(tQR) < 1e-10 ) tQC = 0.0;
  else tQC=tQL/tQR;
  double tQP=atan2(tQS,tQO);

  fDenominator->Fill(tQR,tQP,tQC);

  // EMCICs   
  AliFemtoLorentzVector tMom1 = pair->Track1()->FourMomentum();
  AliFemtoLorentzVector tMom2 = pair->Track2()->FourMomentum();
  double tE1 = tMom1.e();
  double tE2 = tMom2.e();
  double tPz1 = tMom1.pz();
  double tPz2 = tMom2.pz();
  
  TVector2 tPt1;  
  TVector2 tPt2; 
  tPt1.Set(tMom1.px(),tMom1.py());
  tPt2.Set(tMom2.px(),tMom2.py());
  double tPt1DotPt2 = tPt1*tPt2;
  
  fEnergyTotalMix->Fill(tQR,tQP,tQC,tE1+tE2);
  fEnergyMultMix->Fill(tQR,tQP,tQC,tE1*tE2);
  fPzMultMix->Fill(tQR,tQP,tQC,tPz1*tPz2);
  fPtMultMix->Fill(tQR,tQP,tQC,tPt1DotPt2);
  
}
Commit	Line	Data
fa35b516	1	/***************************************************************************
	2	*
	3	* $Id: AliFemtoCorrFctn3DSphercicalEMCIC.h $
	4	*
	5	* Author: Nicolas Bock, Ohio State University, bock@mps.ohio-state.edu
	6	***************************************************************************
	7	*
	8	* Description: Calculates of the 3D Correlation Function in Spherical
	9	* coordinates, and also produces histograms to calculate
	10	* Energy Momentum Conservation Induced Correlations (EMCICs)
	11	*
	12	* This Class produces the following histograms as function of Q, theta, phi
	13	* (for both real and mixed pairs):
	14	* 1) E1 + E2
	15	* 2) E1 * E2
	16	* 3) Pt1*Pt2
	17	* 4) Pz1*Pz2
	18	*
	19	* This class is similar to AliFemtoCorrFctn3DSpherical, but it uses Q
	20	* instead of K to do the binning.
	21	*
	22	* NOTE: The EMCIC histograms are not averaged in this class, to obtain
	23	* the average, the user needs to divide the real pair histograms by
	24	* the numerator, and the mixed pair histograms by the denominator.
	25	*
	26	***************************************************************************
	27	*
	28	**************************************************************************/
	29
	30	#include "AliFemtoCorrFctn3DSphericalEMCIC.h"
746bf34d	31	#include "AliFemtoKTPairCut.h"
fa35b516	32	#include <TMath.h>
	33	#include <TVector2.h>
	34	#include <cstdio>
	35
	36	#ifdef __ROOT__
	37	ClassImp(AliFemtoCorrFctn3DSphericalEMCIC)
	38	#endif
	39
	40
	41	//____________________________
	42	AliFemtoCorrFctn3DSphericalEMCIC::AliFemtoCorrFctn3DSphericalEMCIC(char* title, const int& nqbins, const float& QLo, const float& QHi, const int& nphibins, const int& ncthetabins):
	43	AliFemtoCorrFctn(),
	44	fNumerator(0),
	45	fDenominator(0),
943d9a85	46	/*fEnergyTotalReal(0),
fa35b516	47	fEnergyMultReal(0),
fa35b516	48	fPzMultReal(0),
943d9a85	49	fPtMultReal(0),*/
fa35b516	50	fEnergyTotalMix (0),
	51	fEnergyMultMix (0),
	52	fPzMultMix(0),
	53	fPtMultMix(0),
	54	fPairCut(0x0)
	55	{
	56
	57	// To better sample on phi shift bin low edge by binsize/2 = Pi/numBins.
	58	Double_t shiftPhi=TMath::Pi()/nphibins;
	59	// set up numerator
9f36b0a8	60	char tTitNum[101] = "Num";
adecdc37	61	strncat(tTitNum,title, 100);
fa35b516	62	fNumerator = new TH3D(tTitNum,title,nqbins,QLo,QHi,nphibins,
	63	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
	64	// set up denominator
9f36b0a8	65	char tTitDen[101] = "Den";
adecdc37	66	strncat(tTitDen,title, 100);
fa35b516	67	fDenominator = new TH3D(tTitDen,title,nqbins,QLo,QHi,nphibins,
	68	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
	69
	70	//Added histograms to calculate EMCICs , Nicolas Bock 19.01.2010
	71	//Setup EnergyTotalReal
9f36b0a8	72	/*char tTitNum1[101] = "ESumReal";
adecdc37	73	strncat(tTitNum1,title, 100);
fa35b516	74	fEnergyTotalReal = new TH3D(tTitNum1,title,nqbins,QLo,QHi,nphibins,
	75	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
	76
	77	//Setup EnergyMultReal
9f36b0a8	78	char tTitNum2[101] = "EMultReal";
adecdc37	79	strncat(tTitNum2,title, 100);
fa35b516	80	fEnergyMultReal = new TH3D(tTitNum2,title,nqbins,QLo,QHi,nphibins,
	81	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
	82
	83	//Setup Pz MultReal
9f36b0a8	84	char tTitNum3[101] = "PzMultReal";
adecdc37	85	strncat(tTitNum3,title, 100);
fa35b516	86	fPzMultReal = new TH3D(tTitNum3,title,nqbins,QLo,QHi,nphibins,
	87	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
	88
	89	//Setup Pt MultReal
9f36b0a8	90	char tTitNum4[101] = "PtMultReal";
adecdc37	91	strncat(tTitNum4,title, 100);
fa35b516	92	fPtMultReal = new TH3D(tTitNum4,title,nqbins,QLo,QHi,nphibins,
fa35b516	93	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
943d9a85	94	*/
fa35b516	95
	96
	97	//Setup EnergyTotalMix
9f36b0a8	98	char tTitNum5[101] = "ESumMix";
adecdc37	99	strncat(tTitNum5,title, 100);
fa35b516	100	fEnergyTotalMix = new TH3D(tTitNum5,title,nqbins,QLo,QHi,nphibins,
	101	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
	102
	103	//Setup EnergyMultMix
9f36b0a8	104	char tTitNum6[101] = "EMultMix";
adecdc37	105	strncat(tTitNum6,title, 100);
fa35b516	106	fEnergyMultMix = new TH3D(tTitNum6,title,nqbins,QLo,QHi,nphibins,
	107	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
	108
	109	//Setup Pz MultMix
9f36b0a8	110	char tTitNum7[101] = "PzMultMix";
adecdc37	111	strncat(tTitNum7,title, 100);
fa35b516	112	fPzMultMix = new TH3D(tTitNum7,title,nqbins,QLo,QHi,nphibins,
	113	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,ncthetabins,-1.0,1.0);
	114
	115	//Setup Pt MultMix
9f36b0a8	116	char tTitNum8[101] = "PtMultMix";
adecdc37	117	strncat(tTitNum8,title, 100);
fa35b516	118	fPtMultMix = new TH3D(tTitNum8,title,nqbins,QLo,QHi,nphibins,
	119	-TMath::Pi()-shiftPhi,TMath::Pi()-shiftPhi,
	120	ncthetabins,-1.0,1.0);
	121
	122	// to enable error bar calculation...
	123	fNumerator->Sumw2();
	124	fDenominator->Sumw2();
943d9a85	125	/*fEnergyTotalReal->Sumw2();
fa35b516	126	fEnergyMultReal->Sumw2();
fa35b516	127	fPzMultReal->Sumw2();
943d9a85	128	fPtMultReal->Sumw2(); */
fa35b516	129	fEnergyTotalMix->Sumw2();
	130	fEnergyMultMix->Sumw2();
	131	fPzMultMix->Sumw2();
	132	fPtMultMix->Sumw2();
	133
	134
	135	}
	136
	137	AliFemtoCorrFctn3DSphericalEMCIC::AliFemtoCorrFctn3DSphericalEMCIC(const AliFemtoCorrFctn3DSphericalEMCIC& aCorrFctn) :
	138	AliFemtoCorrFctn(),
	139	fNumerator(0),
	140	fDenominator(0),
943d9a85	141	/*fEnergyTotalReal(0),
fa35b516	142	fEnergyMultReal(0),
fa35b516	143	fPzMultReal(0),
943d9a85	144	fPtMultReal(0), */
fa35b516	145	fEnergyTotalMix (0),
	146	fEnergyMultMix (0),
	147	fPzMultMix(0),
	148	fPtMultMix(0),
	149	fPairCut(0x0)
	150	{
	151	// Copy constructor
	152	fNumerator = new TH3D(*aCorrFctn.fNumerator);
	153	fDenominator = new TH3D(*aCorrFctn.fDenominator);
943d9a85	154	/fEnergyTotalReal = new TH3D(aCorrFctn.fEnergyTotalReal);
fa35b516	155	fEnergyMultReal = new TH3D(*aCorrFctn.fEnergyMultReal);
fa35b516	156	fPzMultReal = new TH3D(*aCorrFctn.fPzMultReal);
943d9a85	157	fPtMultReal = new TH3D(aCorrFctn.fPtMultReal);/
fa35b516	158	fEnergyTotalMix = new TH3D(*aCorrFctn.fEnergyTotalMix);
	159	fEnergyMultMix = new TH3D(*aCorrFctn.fEnergyMultMix);
	160	fPzMultMix = new TH3D(*aCorrFctn.fPzMultMix);
	161	fPtMultMix = new TH3D(*aCorrFctn.fPtMultMix);
	162	fPairCut = aCorrFctn.fPairCut;
	163	}
	164	//____________________________
	165	AliFemtoCorrFctn3DSphericalEMCIC::~AliFemtoCorrFctn3DSphericalEMCIC(){
	166	// Destructor
	167	delete fNumerator;
	168	delete fDenominator;
943d9a85	169	/*delete fEnergyTotalReal;
fa35b516	170	delete fEnergyMultReal;
fa35b516	171	delete fPzMultReal;
943d9a85	172	delete fPtMultReal; */
fa35b516	173	delete fEnergyTotalMix;
	174	delete fEnergyMultMix;
	175	delete fPzMultMix;
	176	delete fPtMultMix;
	177	}
	178	//_________________________
	179	AliFemtoCorrFctn3DSphericalEMCIC& AliFemtoCorrFctn3DSphericalEMCIC::operator=(const AliFemtoCorrFctn3DSphericalEMCIC& aCorrFctn)
	180	{
	181	// assignment operator
	182	if (this == &aCorrFctn)
	183	return *this;
	184
	185	if (fNumerator) delete fNumerator;
	186	fNumerator = new TH3D(*aCorrFctn.fNumerator);
	187	if (fDenominator) delete fDenominator;
	188	fDenominator = new TH3D(*aCorrFctn.fDenominator);
943d9a85	189	/* if (fEnergyTotalReal) delete fEnergyTotalReal;
fa35b516	190	fEnergyTotalReal = new TH3D(*aCorrFctn.fEnergyTotalReal);
	191	if (fEnergyMultReal) delete fEnergyMultReal;
	192	fEnergyMultReal = new TH3D(*aCorrFctn.fEnergyMultReal);
	193	if (fPzMultReal) delete fPzMultReal;
	194	fPzMultReal = new TH3D(*aCorrFctn.fPzMultReal);
	195	if (fPtMultReal) delete fPtMultReal;
943d9a85	196	fPtMultReal = new TH3D(aCorrFctn.fPtMultReal); /
fa35b516	197	if (fEnergyTotalMix) delete fEnergyTotalMix;
	198	fEnergyTotalMix = new TH3D(*aCorrFctn.fEnergyTotalMix);
	199	if (fEnergyMultMix) delete fEnergyMultMix;
	200	fEnergyMultMix = new TH3D(*aCorrFctn.fEnergyMultMix);
	201	if (fPzMultMix) delete fPzMultMix;
	202	fPzMultMix = new TH3D(*aCorrFctn.fPzMultMix);
	203	if (fPtMultMix) delete fPtMultMix;
	204	fPtMultMix = new TH3D(*aCorrFctn.fPtMultMix);
	205	fPairCut = aCorrFctn.fPairCut;
	206
	207	return *this;
	208	}
	209
	210	//_________________________
	211	void AliFemtoCorrFctn3DSphericalEMCIC::WriteOutHistos(){
	212	// Write out all histograms to file
	213	fNumerator->Write();
	214	fDenominator->Write();
943d9a85	215	/*fEnergyTotalReal->Write();
fa35b516	216	fEnergyMultReal->Write();
fa35b516	217	fPzMultReal->Write();
943d9a85	218	fPtMultReal->Write(); */
fa35b516	219	fEnergyTotalMix->Write();
	220	fEnergyMultMix->Write();
	221	fPzMultMix->Write();
	222	fPtMultMix->Write();
	223	}
	224	//______________________________
	225	TList* AliFemtoCorrFctn3DSphericalEMCIC::GetOutputList()
	226	{
	227	// Prepare the list of objects to be written to the output
	228	TList *tOutputList = new TList();
	229
	230	tOutputList->Add(fNumerator);
	231	tOutputList->Add(fDenominator);
943d9a85	232	/* tOutputList->Add(fEnergyTotalReal);
fa35b516	233	tOutputList->Add(fEnergyMultReal);
fa35b516	234	tOutputList->Add(fPzMultReal);
943d9a85	235	tOutputList->Add(fPtMultReal); */
fa35b516	236	tOutputList->Add(fEnergyTotalMix );
	237	tOutputList->Add(fEnergyMultMix );
	238	tOutputList->Add(fPzMultMix);
	239	tOutputList->Add(fPtMultMix);
	240	return tOutputList;
	241	}
	242
	243	//_________________________
	244	void AliFemtoCorrFctn3DSphericalEMCIC::Finish(){
	245	// here is where we should normalize, fit, etc...
	246	}
	247
	248	//____________________________
	249	AliFemtoString AliFemtoCorrFctn3DSphericalEMCIC::Report(){
	250	// Construct the report
	251	string stemp = "PRF Frame SphericalEMCIC 3D Correlation Function Report:\n";
	252	char ctemp[100];
adecdc37	253	snprintf(ctemp , 100, "Number of entries in numerator:\t%E\n",fNumerator->GetEntries());
fa35b516	254	stemp += ctemp;
adecdc37	255	snprintf(ctemp , 100, "Number of entries in denominator:\t%E\n",fDenominator->GetEntries());
fa35b516	256	stemp += ctemp;
	257
	258	if (fPairCut){
adecdc37	259	snprintf(ctemp , 100, "Here is the PairCut specific to this CorrFctn\n");
fa35b516	260	stemp += ctemp;
	261	stemp += fPairCut->Report();
	262	}
	263	else{
adecdc37	264	snprintf(ctemp , 100, "No PairCut specific to this CorrFctn\n");
fa35b516	265	stemp += ctemp;
	266	}
	267
	268	//
	269	AliFemtoString returnThis = stemp;
	270	return returnThis;
	271	}
	272	//____________________________
	273	void AliFemtoCorrFctn3DSphericalEMCIC::AddRealPair( AliFemtoPair* pair){
	274	// perform operations on real pairs
	275	if (fPairCut){
746bf34d	276	AliFemtoKTPairCut ktc = dynamic_cast<AliFemtoKTPairCut >(fPairCut);
	277	if (!ktc){
	278	if (!(fPairCut->Pass(pair))) return;
	279	}
	280	else
	281	if (!(ktc->Pass(pair))) return;
fa35b516	282	}
fa35b516	283
943d9a85	284	//
	285	double tQO = pair->QOutCMS();
	286	double tQS = pair->QSideCMS();
	287	double tQL = pair->QLongCMS();
fa35b516	288
943d9a85	289	double tQR = sqrt(tQOtQO + tQStQS + tQL*tQL);
	290	double tQC = 0;
	291	if ( fabs(tQR) < 1e-10 ) tQC = 0.0;
	292	else tQC = tQL/tQR;
	293	double tQP = atan2(tQS,tQO);
fa35b516	294
943d9a85	295	fNumerator->Fill(tQR,tQP,tQC);
fa35b516	296
fa35b516	297	// EMCICs
943d9a85	298	/*AliFemtoLorentzVector tMom1 = pair->Track1()->FourMomentum();
fa35b516	299	AliFemtoLorentzVector tMom2 = pair->Track2()->FourMomentum();
	300	double tE1 = tMom1.e();
	301	double tE2 = tMom2.e();
	302	double tPz1 = tMom1.pz();
	303	double tPz2 = tMom2.pz();
943d9a85	304
	305	TVector2 tPt1;
	306	TVector2 tPt2;
fa35b516	307	tPt1.Set(tMom1.px(),tMom1.py());
fa35b516	308	tPt2.Set(tMom2.px(),tMom2.py());
943d9a85	309	double tPt1DotPt2 = tPt1*tPt2;
fa35b516	310
943d9a85	311	fEnergyTotalReal->Fill(tQR,tQP,tQC,tE1+tE2);
	312	fEnergyMultReal->Fill(tQR,tQP,tQC,tE1*tE2);
	313	fPzMultReal->Fill(tQR,tQP,tQC,tPz1*tPz2);
	314	fPtMultReal->Fill(tQR,tQP,tQC,tPt1DotPt2);*/
fa35b516	315
	316
	317	}
	318	//____________________________
	319	void AliFemtoCorrFctn3DSphericalEMCIC::AddMixedPair( AliFemtoPair* pair){
	320	// perform operations on mixed pairs
	321	if (fPairCut){
746bf34d	322	AliFemtoKTPairCut ktc = dynamic_cast<AliFemtoKTPairCut >(fPairCut);
	323	if (!ktc){
	324	if (!(fPairCut->Pass(pair))) return;
	325	}
	326	else
	327	if (!(ktc->Pass(pair))) return;
fa35b516	328	}
746bf34d	329
	330
	331
fa35b516	332	// //Changed K to Q to be in LCMS, N. Bock
943d9a85	333	double tQO = pair->QOutCMS();
	334	double tQS = pair->QSideCMS();
	335	double tQL = pair->QLongCMS();
fa35b516	336
943d9a85	337	double tQR = sqrt(tQOtQO + tQStQS + tQL*tQL);
	338	double tQC;
	339	if ( fabs(tQR) < 1e-10 ) tQC = 0.0;
	340	else tQC=tQL/tQR;
	341	double tQP=atan2(tQS,tQO);
fa35b516	342
943d9a85	343	fDenominator->Fill(tQR,tQP,tQC);
fa35b516	344
	345	// EMCICs
	346	AliFemtoLorentzVector tMom1 = pair->Track1()->FourMomentum();
	347	AliFemtoLorentzVector tMom2 = pair->Track2()->FourMomentum();
	348	double tE1 = tMom1.e();
	349	double tE2 = tMom2.e();
	350	double tPz1 = tMom1.pz();
	351	double tPz2 = tMom2.pz();
943d9a85	352
	353	TVector2 tPt1;
	354	TVector2 tPt2;
fa35b516	355	tPt1.Set(tMom1.px(),tMom1.py());
fa35b516	356	tPt2.Set(tMom2.px(),tMom2.py());
943d9a85	357	double tPt1DotPt2 = tPt1*tPt2;
fa35b516	358
943d9a85	359	fEnergyTotalMix->Fill(tQR,tQP,tQC,tE1+tE2);
	360	fEnergyMultMix->Fill(tQR,tQP,tQC,tE1*tE2);
	361	fPzMultMix->Fill(tQR,tQP,tQC,tPz1*tPz2);
	362	fPtMultMix->Fill(tQR,tQP,tQC,tPt1DotPt2);
fa35b516	363
fa35b516	364	}
746bf34d	365