***************************************************************************
*
* $Log$
+ * Revision 1.4 2007/05/03 09:46:10 akisiel
+ * Fixing Effective C++ warnings
+ *
* Revision 1.3 2007/04/27 07:25:59 akisiel
* Make revisions needed for compilation from the main AliRoot tree
*
***************************************************************************
*
* $Log$
+ * Revision 1.1.1.1 2007/04/25 15:38:41 panos
+ * Importing the HBT code dir
+ *
* Revision 1.4 2007/03/20 09:37:13 mchojnacki
* *** empty log message ***
*
//#include "StMaker.h"
//#endif
-#include "Base/AliFemtoTrackCut.h"
+#include "AliFemtoTrackCut.h"
class AliFemtoESDTrackCut : public AliFemtoTrackCut
{
--- /dev/null
+////////////////////////////////////////////////////////////////////////////////
+/// ///
+/// AliFemtoModelBPLCMSCorrFctn - the class for correlation function which ///
+/// uses the model framework and weight generation and calculated the 3D ///
+/// correlation function in the Bertsh-Pratt LCMS system ///
+/// Authors: Adam Kisiel, kisiel@mps.ohio-state.edu ///
+/// ///
+////////////////////////////////////////////////////////////////////////////////
+#include "AliFemtoModelBPLCMSCorrFctn.h"
+#include <cstdio>
+
+#ifdef __ROOT__
+ClassImp(AliFemtoModelBPLCMSCorrFctn)
+#endif
+
+//____________________________
+AliFemtoModelBPLCMSCorrFctn::AliFemtoModelBPLCMSCorrFctn(char* title, const int& nbins, const float& QLo, const float& QHi)
+ :
+ fNumerator(0),
+ fDenominator(0),
+ fQinvHisto(0)
+{
+
+ // set up numerator
+ char TitNum[100] = "Num";
+ strcat(TitNum,title);
+ fNumerator = new TH3D(TitNum,title,nbins,QLo,QHi,nbins,QLo,QHi,nbins,QLo,QHi);
+ // set up denominator
+ char TitDen[100] = "Den";
+ strcat(TitDen,title);
+ fDenominator = new TH3D(TitDen,title,nbins,QLo,QHi,nbins,QLo,QHi,nbins,QLo,QHi);
+ // set up ave qInv
+ char TitQinv[100] = "Qinv";
+ strcat(TitQinv,title);
+ fQinvHisto = new TH3D(TitQinv,title,nbins,QLo,QHi,nbins,QLo,QHi,nbins,QLo,QHi);
+
+ // to enable error bar calculation...
+ fNumerator->Sumw2();
+ fDenominator->Sumw2();
+}
+
+AliFemtoModelBPLCMSCorrFctn::AliFemtoModelBPLCMSCorrFctn(const AliFemtoModelBPLCMSCorrFctn& aCorrFctn) :
+ fNumerator(0),
+ fDenominator(0),
+ fQinvHisto(0)
+{
+ fNumerator = new TH3D(*aCorrFctn.fNumerator);
+ fDenominator = new TH3D (*aCorrFctn.fDenominator);
+ fQinvHisto = new TH3D(*aCorrFctn.fQinvHisto);
+}
+//____________________________
+AliFemtoModelBPLCMSCorrFctn::~AliFemtoModelBPLCMSCorrFctn(){
+ delete fNumerator;
+ delete fDenominator;
+ delete fQinvHisto;
+}
+//_________________________
+AliFemtoModelBPLCMSCorrFctn& AliFemtoModelBPLCMSCorrFctn::operator=(const AliFemtoModelBPLCMSCorrFctn& aCorrFctn)
+{
+ if (this == &aCorrFctn)
+ return *this;
+ if (fNumerator) delete fNumerator;
+ fNumerator = new TH3D(*aCorrFctn.fNumerator);
+ if (fDenominator) delete fDenominator;
+ fDenominator = new TH3D(*aCorrFctn.fDenominator);
+ if (fQinvHisto) delete fQinvHisto;
+ fQinvHisto = new TH3D(*aCorrFctn.fQinvHisto);
+
+ return *this;
+}
+
+//_________________________
+void AliFemtoModelBPLCMSCorrFctn::WriteOutHistos(){
+
+ fNumerator->Write();
+ fDenominator->Write();
+ fQinvHisto->Write();
+}
+
+//_________________________
+void AliFemtoModelBPLCMSCorrFctn::Finish(){
+ fQinvHisto->Divide(fDenominator);
+}
+
+//____________________________
+AliFemtoString AliFemtoModelBPLCMSCorrFctn::Report(){
+ string stemp = "LCMS Frame Bertsch-Pratt 3D Correlation Function Report:\n";
+ char ctemp[100];
+ sprintf(ctemp,"Number of entries in numerator:\t%E\n",fNumerator->GetEntries());
+ stemp += ctemp;
+ sprintf(ctemp,"Number of entries in denominator:\t%E\n",fDenominator->GetEntries());
+ stemp += ctemp;
+ sprintf(ctemp,"Number of entries in ratio:\t%E\n",fRatio->GetEntries());
+ stemp += ctemp;
+ sprintf(ctemp,"Normalization region in Qinv was:\t%E\t%E\n",fQinvNormLo,fQinvNormHi);
+ stemp += ctemp;
+ sprintf(ctemp,"Number of pairs in Normalization region was:\n");
+ stemp += ctemp;
+ sprintf(ctemp,"In numerator:\t%lu\t In denominator:\t%lu\n",fNumRealsNorm,fNumMixedNorm);
+ stemp += ctemp;
+ /* if (fCorrection)
+ {
+ float radius = fCorrection->GetRadius();
+ sprintf(ctemp,"Coulomb correction used radius of\t%E\n",radius);
+ }
+ else
+ {
+ sprintf(ctemp,"No Coulomb Correction applied to this CorrFctn\n");
+ }
+ stemp += ctemp;
+ */
+
+ if (fPairCut){
+ sprintf(ctemp,"Here is the PairCut specific to this CorrFctn\n");
+ stemp += ctemp;
+ stemp += fPairCut->Report();
+ }
+ else{
+ sprintf(ctemp,"No PairCut specific to this CorrFctn\n");
+ stemp += ctemp;
+ }
+
+ //
+ AliFemtoString returnThis = stemp;
+ return returnThis;
+}
+//____________________________
+void AliFemtoModelBPLCMSCorrFctn::AddRealPair( AliFemtoPair* pair){
+
+ if (fPairCut){
+ if (!(fPairCut->Pass(pair))) return;
+ }
+
+ double Qinv = fabs(pair->qInv()); // note - qInv() will be negative for identical pairs...
+ if ((Qinv < fQinvNormHi) && (Qinv > fQinvNormLo)) fNumRealsNorm++;
+ double qOut = fabs(pair->qOutCMS());
+ double qSide = fabs(pair->qSideCMS());
+ double qLong = fabs(pair->qLongCMS());
+
+ fNumerator->Fill(qOut,qSide,qLong);
+}
+//____________________________
+void AliFemtoModelBPLCMSCorrFctn::AddMixedPair( AliFemtoPair* pair){
+
+ if (fPairCut){
+ if (!(fPairCut->Pass(pair))) return;
+ }
+
+ // double CoulombWeight = (fCorrection ? fCorrection->CoulombCorrect(pair) : 1.0);
+ double CoulombWeight = 1.0;
+
+ double Qinv = fabs(pair->qInv()); // note - qInv() will be negative for identical pairs...
+ if ((Qinv < fQinvNormHi) && (Qinv > fQinvNormLo)) fNumMixedNorm++;
+ double qOut = fabs(pair->qOutCMS());
+ double qSide = fabs(pair->qSideCMS());
+ double qLong = fabs(pair->qLongCMS());
+
+ fDenominator->Fill(qOut,qSide,qLong,CoulombWeight);
+ // fUncorrectedDenominator->Fill(qOut,qSide,qLong,1.0);
+ fQinvHisto->Fill(qOut,qSide,qLong,Qinv);
+
+ /*
+ // now for the momentum resolution stuff...
+ if (fSmearPair){
+ double CorrWeight = 1.0 +
+ fLambda*exp((-qOut*qOut*fRout2 -qSide*qSide*fRside2 -qLong*qLong*fRlong2)/0.038936366329);
+ CorrWeight *= CoulombWeight; // impt.
+
+ fIDNumHisto->Fill(qOut,qSide,qLong,CorrWeight);
+ fIDDenHisto->Fill(qOut,qSide,qLong,CoulombWeight);
+
+ fSmearPair->SetUnsmearedPair(pair);
+ double qOut_prime = fabs(fSmearPair->SmearedPair().qOutCMS());
+ double qSide_prime = fabs(fSmearPair->SmearedPair().qSideCMS());
+ double qLong_prime = fabs(fSmearPair->SmearedPair().qLongCMS());
+
+ fSMNumHisto->Fill(qOut_prime,qSide_prime,qLong_prime,CorrWeight);
+
+ double SmearedCoulombWeight = ( fCorrection ?
+ fCorrection->CoulombCorrect(&(fSmearPair->SmearedPair())) :
+ 1.0);
+
+ fSMDenHisto->Fill(qOut_prime,qSide_prime,qLong_prime,SmearedCoulombWeight);
+ }
+ */
+}
+
+
--- /dev/null
+////////////////////////////////////////////////////////////////////////////////
+/// ///
+/// AliFemtoModelBPLCMSCorrFctn - the class for correlation function which ///
+/// uses the model framework and weight generation and calculated the 3D ///
+/// correlation function in the Bertsh-Pratt LCMS system ///
+/// Authors: Adam Kisiel, kisiel@mps.ohio-state.edu ///
+/// ///
+////////////////////////////////////////////////////////////////////////////////
+#ifndef AliFemtoModelBPLCMS3DCorrFctn_hh
+#define AliFemtoModelBPLCMS3DCorrFctn_hh
+
+#include "AliFemtoCorrFctn.h"
+#include "AliFemtoModelCorrFctn.h"
+#include "AliFemtoPairCut.h"
+#include "TH3D.h"
+
+class AliFemtoModelBPLCMS3DCorrFctn : public AliFemtoModelCorrFctn {
+ public:
+ AliFemtoModelBPLCMSCorrFctn();
+ AliFemtoModelBPLCMSCorrFctn(char* title, const int& nbins, const float& QLo, const float& QHi);
+ AliFemtoModelBPLCMSCorrFctn(const AliFemtoModelBPLCMSCorrFctn& aCorrFctn);
+ virtual ~AliFemtoModelBPLCMSCorrFctn();
+
+ AliFemtoModelBPLCMSCorrFctn& operator=(const AliFemtoModelBPLCMSCorrFctn& aCorrFctn);
+
+ virtual AliFemtoString Report();
+ virtual void AddRealPair(AliFemtoPair*);
+ virtual void AddMixedPair(AliFemtoPair*);
+
+ virtual void Finish();
+
+ virtual void Write();
+
+ virtual AliFemtoModelCorrFctnSource* Clone();
+
+ TH3D* Numerator();
+ TH3D* Denominator();
+ TH3D* QinvHisto();
+
+private:
+ TH3D* fNumerator;
+ TH3D* fDenominator;
+
+ TH3D* fQinvHisto;
+
+#ifdef __ROOT__
+ ClassDef(AliFemtoModelBPLCMSCorrFctn, 1)
+#endif
+};
+
+inline TH3D* AliFemtoModelBPLCMSCorrFctn::Numerator(){return fNumerator;}
+inline TH3D* AliFemtoModelBPLCMSCorrFctn::Denominator(){return fDenominator;}
+inline TH3D* AliFemtoModelBPLCMSCorrFctn::QinvHisto(){return fQinvHisto;}
+
+#endif
+
ClassImp(AliFemtoModelCorrFctnSource, 1)
#endif
-#include "Model/AliFemtoModelGausLCMSFreezeOutGenerator.h"
-#include "Model/AliFemtoModelHiddenInfo.h"
-#include "Model/AliFemtoModelCorrFctnSource.h"
+#include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
+#include "AliFemtoModelHiddenInfo.h"
+#include "AliFemtoModelCorrFctnSource.h"
//_______________________
AliFemtoModelCorrFctnSource::AliFemtoModelCorrFctnSource():
AliFemtoModelCorrFctn::Write();
}
//_______________________
-AliFemtoModelCorrFctnSource* AliFemtoModelCorrFctnSource::Clone()
+AliFemtoModelCorrFctn* AliFemtoModelCorrFctnSource::Clone()
{
AliFemtoModelCorrFctnSource *tCopy = new AliFemtoModelCorrFctnSource(*this);
#ifndef AliFemtoModelCorrFctnSource_hh
#define AliFemtoModelCorrFctnSource_hh
-#include "Base/AliFemtoCorrFctn.h"
-#include "Infrastructure/AliFemtoPair.h"
-#include "Model/AliFemtoModelManager.h"
-#include "Model/AliFemtoModelCorrFctn.h"
+#include "AliFemtoCorrFctn.h"
+#include "AliFemtoPair.h"
+#include "AliFemtoModelManager.h"
+#include "AliFemtoModelCorrFctn.h"
class AliFemtoModelCorrFctnSource: public AliFemtoModelCorrFctn {
virtual void Write();
- virtual AliFemtoModelCorrFctnSource* Clone();
+ virtual AliFemtoModelCorrFctn* Clone();
protected:
#include "math.h"
#include "AliFemtoModelGausRinvFreezeOutGenerator.h"
-#include "Model/AliFemtoModelHiddenInfo.h"
+#include "AliFemtoModelHiddenInfo.h"
#include "AliFemtoLorentzVector.h"
//_______________________
// Generate two particle emission points with respect
// to their pair momentum
// The source is the 3D Gaussian ellipsoid in the LCMS frame
- AliFemtoModelHiddenInfo *inf1 = (AliFemtoModelHiddenInfo *) aPair->track1()->HiddenInfo();
- AliFemtoModelHiddenInfo *inf2 = (AliFemtoModelHiddenInfo *) aPair->track2()->HiddenInfo();
+ AliFemtoModelHiddenInfo *inf1 = (AliFemtoModelHiddenInfo *) aPair->Track1()->HiddenInfo();
+ AliFemtoModelHiddenInfo *inf2 = (AliFemtoModelHiddenInfo *) aPair->Track2()->HiddenInfo();
if ((!inf1) || (!inf2)) { cout << "Hidden info not created! " << endl; exit(kFALSE); }
#ifndef AliFemtoModelGausRinvFreezeOutGenerator_hh
#define AliFemtoModelGausRinvFreezeOutGenerator_hh
-#include "Model/AliFemtoModelFreezeOutGenerator.h"
+#include "AliFemtoModelFreezeOutGenerator.h"
#include "TRandom.h"
***************************************************************************
*
* $Log$
+ * Revision 1.4 2007/05/03 09:46:10 akisiel
+ * Fixing Effective C++ warnings
+ *
* Revision 1.3 2007/04/27 07:25:59 akisiel
* Make revisions needed for compilation from the main AliRoot tree
*
{
//bool temp = true;
//temp ? fNPairsPassed++ : fNPairsFailed++;
- if ((fabs(pair->qLongCMS())<fQlong[0])||(fabs(pair->qLongCMS())>fQlong[1]))
+ if ((fabs(pair->QLongCMS())<fQlong[0])||(fabs(pair->QLongCMS())>fQlong[1]))
{
fNPairsFailed++;
return false;
}
- if ((fabs(pair->qOutCMS())<fQout[0])||(fabs(pair->qOutCMS())>fQout[1]))
+ if ((fabs(pair->QOutCMS())<fQout[0])||(fabs(pair->QOutCMS())>fQout[1]))
{
fNPairsFailed++;
return false;
}
- if ((fabs(pair->qSideCMS())<fQside[0])||(fabs(pair->qSideCMS())>fQside[1]))
+ if ((fabs(pair->QSideCMS())<fQside[0])||(fabs(pair->QSideCMS())>fQside[1]))
{
fNPairsFailed++;
return false;
***************************************************************************
*
* $Log$
+ * Revision 1.1.1.1 2007/04/25 15:38:41 panos
+ * Importing the HBT code dir
+ *
* Revision 1.1.1.1 2007/03/07 10:14:49 mchojnacki
* First version on CVS
*
//#include "StMaker.h"
//#endif
-#include "Base/AliFemtoPairCut.h"
+#include "AliFemtoPairCut.h"
class AliFemtoQPairCut : public AliFemtoPairCut{
public:
////////////////////////////////////////////////////////////////////////////////
#include "AliFemtoShareQualityCorrFctn.h"
-//#include "Infrastructure/AliFemtoHisto.hh"
+//#include "AliFemtoHisto.hh"
#include <cstdio>
#ifdef __ROOT__
return returnThis;
}
//____________________________
-void AliFemtoShareQualityCorrFctn::AddRealPair(const AliFemtoPair* pair){
- double Qinv = fabs(pair->qInv()); // note - qInv() will be negative for identical pairs...
+void AliFemtoShareQualityCorrFctn::AddRealPair( AliFemtoPair* pair){
+ double Qinv = fabs(pair->QInv()); // note - qInv() will be negative for identical pairs...
Int_t nh = 0;
Int_t an = 0;
Int_t ns = 0;
- for (unsigned int imap=0; imap<pair->track1()->Track()->TPCclusters().GetNbits(); imap++) {
+ for (unsigned int imap=0; imap<pair->Track1()->Track()->TPCclusters().GetNbits(); imap++) {
// If both have clusters in the same row
- if (pair->track1()->Track()->TPCclusters().TestBitNumber(imap) &&
- pair->track2()->Track()->TPCclusters().TestBitNumber(imap)) {
+ if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap) &&
+ pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) {
// Do they share it ?
- if (pair->track1()->Track()->TPCsharing().TestBitNumber(imap) ||
- pair->track2()->Track()->TPCsharing().TestBitNumber(imap))
+ if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap) &&
+ pair->Track2()->Track()->TPCsharing().TestBitNumber(imap))
{
if (Qinv < 0.01) {
cout << "Shared cluster in row " << imap << endl;
nh+=2;
}
}
- else if (pair->track1()->Track()->TPCclusters().TestBitNumber(imap) ||
- pair->track2()->Track()->TPCclusters().TestBitNumber(imap)) {
+ else if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap) ||
+ pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) {
// One track has a hit, the other does not
an++;
nh++;
if (Qinv < 0.01) {
cout << "Qinv of the pair is " << Qinv << endl;
cout << "Clusters: " << endl;
- for (unsigned int imap=0; imap<pair->track1()->Track()->TPCclusters().GetNbits(); imap++) {
+ for (unsigned int imap=0; imap<pair->Track1()->Track()->TPCclusters().GetNbits(); imap++) {
cout << imap ;
- if (pair->track1()->Track()->TPCclusters().TestBitNumber(imap)) cout << " 1 ";
+ if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap)) cout << " 1 ";
else cout << " 0 " ;
- if (pair->track2()->Track()->TPCclusters().TestBitNumber(imap)) cout << " 1 ";
+ if (pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) cout << " 1 ";
else cout << " 0 " ;
cout << " ";
- if (pair->track1()->Track()->TPCsharing().TestBitNumber(imap)) cout << " S ";
+ if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap)) cout << " S ";
else cout << " X ";
- if (pair->track2()->Track()->TPCsharing().TestBitNumber(imap)) cout << " S ";
+ if (pair->Track2()->Track()->TPCsharing().TestBitNumber(imap)) cout << " S ";
else cout << " X ";
cout << endl;
}
}
if (Qinv < 0.01) {
- cout << "Quality Sharity " << hsmval << " " << hsfval << " " << pair->track1()->Track() << " " << pair->track2()->Track() << endl;
+ cout << "Quality Sharity " << hsmval << " " << hsfval << " " << pair->Track1()->Track() << " " << pair->Track2()->Track() << endl;
}
fShareNumerator->Fill(Qinv, hsfval);
fQualityNumerator->Fill(Qinv, hsmval);
// cout << "AliFemtoShareQualityCorrFctn::AddRealPair : " << pair->qInv() << " " << Qinv <<
- //" " << pair->track1().FourMomentum() << " " << pair->track2().FourMomentum() << endl;
+ //" " << pair->Track1().FourMomentum() << " " << pair->Track2().FourMomentum() << endl;
}
//____________________________
-void AliFemtoShareQualityCorrFctn::AddMixedPair(const AliFemtoPair* pair){
+void AliFemtoShareQualityCorrFctn::AddMixedPair( AliFemtoPair* pair){
double weight = 1.0;
- double Qinv = fabs(pair->qInv()); // note - qInv() will be negative for identical pairs...
+ double Qinv = fabs(pair->QInv()); // note - qInv() will be negative for identical pairs...
Int_t nh = 0;
Int_t an = 0;
Int_t ns = 0;
- for (unsigned int imap=0; imap<pair->track1()->Track()->TPCclusters().GetNbits(); imap++) {
+ for (unsigned int imap=0; imap<pair->Track1()->Track()->TPCclusters().GetNbits(); imap++) {
// If both have clusters in the same row
- if (pair->track1()->Track()->TPCclusters().TestBitNumber(imap) &&
- pair->track2()->Track()->TPCclusters().TestBitNumber(imap)) {
+ if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap) &&
+ pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) {
// Do they share it ?
- if (pair->track1()->Track()->TPCsharing().TestBitNumber(imap) ||
- pair->track2()->Track()->TPCsharing().TestBitNumber(imap))
+ if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap) ||
+ pair->Track2()->Track()->TPCsharing().TestBitNumber(imap))
{
// cout << "A shared cluster !!!" << endl;
// cout << "imap idx1 idx2 "
nh+=2;
}
}
- else if (pair->track1()->Track()->TPCclusters().TestBitNumber(imap) ||
- pair->track2()->Track()->TPCclusters().TestBitNumber(imap)) {
+ else if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap) ||
+ pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) {
// One track has a hit, the other does not
an++;
nh++;
#include "TH1D.h"
#include "TH2D.h"
-#include "Base/AliFemtoCorrFctn.h"
+#include "AliFemtoCorrFctn.h"
class AliFemtoShareQualityCorrFctn : public AliFemtoCorrFctn {
public:
AliFemtoShareQualityCorrFctn& operator=(const AliFemtoShareQualityCorrFctn& aCorrFctn);
virtual AliFemtoString Report();
- virtual void AddRealPair(const AliFemtoPair*);
- virtual void AddMixedPair(const AliFemtoPair*);
+ virtual void AddRealPair(AliFemtoPair*);
+ virtual void AddMixedPair(AliFemtoPair*);
virtual void Finish();
Int_t an = 0;
Int_t ns = 0;
- for (unsigned int imap=0; imap<pair->track1()->Track()->TPCclusters().GetNbits(); imap++) {
+ for (unsigned int imap=0; imap<pair->Track1()->Track()->TPCclusters().GetNbits(); imap++) {
// If both have clusters in the same row
- if (pair->track1()->Track()->TPCclusters().TestBitNumber(imap) &&
- pair->track2()->Track()->TPCclusters().TestBitNumber(imap)) {
+ if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap) &&
+ pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) {
// Do they share it ?
- if (pair->track1()->Track()->TPCsharing().TestBitNumber(imap) &&
- pair->track2()->Track()->TPCsharing().TestBitNumber(imap))
+ if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap) &&
+ pair->Track2()->Track()->TPCsharing().TestBitNumber(imap))
{
// cout << "A shared cluster !!!" << endl;
// cout << "imap idx1 idx2 "
nh+=2;
}
}
- else if (pair->track1()->Track()->TPCclusters().TestBitNumber(imap) ||
- pair->track2()->Track()->TPCclusters().TestBitNumber(imap)) {
+ else if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap) ||
+ pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) {
// One track has a hit, the other does not
an++;
nh++;
}
// if (hsmval > -0.4) {
cout << "Pair quality: " << hsmval << " " << an << " " << nh << " "
- << (pair->track1()->Track()) << " "
- << (pair->track2()->Track()) << endl;
- cout << "Bits: " << pair->track1()->Track()->TPCclusters().GetNbits() << endl;
+ << (pair->Track1()->Track()) << " "
+ << (pair->Track2()->Track()) << endl;
+ cout << "Bits: " << pair->Track1()->Track()->TPCclusters().GetNbits() << endl;
// }
if (hsfval > 0.0) {
cout << "Pair sharity: " << hsfval << " " << ns << " " << nh << " " << hsmval << " " << an << " " << nh << endl;
//#include "StMaker.h"
//#endif
-#include "Base/AliFemtoPairCut.h"
+#include "AliFemtoPairCut.h"
class AliFemtoShareQualityPairCut : public AliFemtoPairCut{
public: