fAODObjArrayName(""), fAddToHistogramsName(""),
fAODCaloCells(0x0),//fAODCaloClusters(0x0),
fCaloPID(0x0), fFidCut(0x0), fIC(0x0),fMCUtils(0x0), fNMS(0x0),
- fCaloUtils(0x0),
- //fAnaOutContainer(0x0),
+ fCaloUtils(0x0),
+ //fAnaOutContainer(0x0),
fHistoPtBins(0), fHistoPtMax(0.), fHistoPtMin(0.),
fHistoPhiBins(0), fHistoPhiMax(0.), fHistoPhiMin(0.),
fHistoEtaBins(0), fHistoEtaMax(0.), fHistoEtaMin(0.),
fHistoMassBins(0), fHistoMassMax(0.), fHistoMassMin(0.),
- fHistoAsymBins(0), fHistoAsymMax(0.), fHistoAsymMin(0.)
+ fHistoAsymBins(0), fHistoAsymMax(0.), fHistoAsymMin(0.)
{
//Default Ctor
//Initialize parameters
InitParameters();
}
-
+/*
//___________________________________________________________
AliAnaPartCorrBaseClass::AliAnaPartCorrBaseClass(const AliAnaPartCorrBaseClass & abc) :
TObject(), fDataMC(abc.fDataMC), fDebug(abc.fDebug),
// cpy ctor
}
-
+*/
+ /*
//_________________________________________________________________________
AliAnaPartCorrBaseClass & AliAnaPartCorrBaseClass::operator = (const AliAnaPartCorrBaseClass & abc)
{
return *this;
}
-
+ */
//____________________________________________________________________________
AliAnaPartCorrBaseClass::~AliAnaPartCorrBaseClass()
{
fMinPt = 0.1 ; //Min pt in particle analysis
fMaxPt = 300. ; //Max pt in particle analysis
- fReader = new AliCaloTrackReader();
- fCaloUtils = new AliCalorimeterUtils();
+ //fReader = new AliCaloTrackReader(); //Initialized in maker
+ //fCaloUtils = new AliCalorimeterUtils();//Initialized in maker
fCaloPID = new AliCaloPID();
fFidCut = new AliFiducialCut();
fIC = new AliIsolationCut();
class AliAnaPartCorrBaseClass : public TObject {
-public:
-
+ public:
AliAnaPartCorrBaseClass() ; // default ctor
+ virtual ~AliAnaPartCorrBaseClass() ; //virtual dtor
+
+ private:
AliAnaPartCorrBaseClass(const AliAnaPartCorrBaseClass & g) ; // cpy ctor
AliAnaPartCorrBaseClass & operator = (const AliAnaPartCorrBaseClass & g) ;//cpy assignment
- virtual ~AliAnaPartCorrBaseClass() ; //virtual dtor
-// virtual void AddAODCaloCluster(AliAODCaloCluster calo) ;
+ public:
+ // virtual void AddAODCaloCluster(AliAODCaloCluster calo) ;
virtual void AddAODParticle(AliAODPWG4Particle part) ;
// virtual void ConnectAODCaloClusters();
TObject(),
fOutputContainer(new TList ), fAnalysisContainer(new TList ),
fMakeHisto(kFALSE), fMakeAOD(kFALSE), fMakeMixing(kFALSE), fAnaDebug(0),
-fReader(0x0), fCaloUtils(0x0), fAODBranchList(new TList ),
-fhNEvents(0x0)
+fReader(new AliCaloTrackReader()), fCaloUtils(new AliCalorimeterUtils()),
+fAODBranchList(new TList ), fhNEvents(0x0)
{
//Default Ctor
if(fAnaDebug > 1 ) printf("*** Analysis Maker Constructor *** \n");
//Initialize parameters, pointers and histograms
- if(!fReader)
- fReader = new AliCaloTrackReader();
- if(!fCaloUtils)
- fCaloUtils = new AliCalorimeterUtils();
-
InitParameters();
}
fOutputContainer(new TList()), fAnalysisContainer(new TList()),
fMakeHisto(maker.fMakeHisto), fMakeAOD(maker.fMakeAOD), fMakeMixing(maker.fMakeMixing),
fAnaDebug(maker.fAnaDebug),
-fReader(new AliCaloTrackReader(*maker.fReader)),
-fCaloUtils(new AliCalorimeterUtils(*maker.fCaloUtils)),
+fReader(),//new AliCaloTrackReader(*maker.fReader)),
+fCaloUtils(),//(new AliCalorimeterUtils(*maker.fCaloUtils)),
fAODBranchList(new TList()),
fhNEvents(maker.fhNEvents)
{
class AliAnaPartCorrMaker : public TObject {
public:
-
AliAnaPartCorrMaker() ; // default ctor
- AliAnaPartCorrMaker(const AliAnaPartCorrMaker & maker) ; // cpy ctor
virtual ~AliAnaPartCorrMaker() ; //virtual dtor
+ AliAnaPartCorrMaker(const AliAnaPartCorrMaker & maker) ; // cpy ctor
private:
AliAnaPartCorrMaker & operator = (const AliAnaPartCorrMaker & ) ;//cpy assignment
fTrackStatus|=AliESDtrack::kITSrefit;
}
-
+/*
//____________________________________________________________________________
AliCaloTrackAODReader::AliCaloTrackAODReader(const AliCaloTrackAODReader & aodr) :
AliCaloTrackReader(aodr), fWriteOutputStdAOD(aodr.fWriteOutputStdAOD)
{
// cpy ctor
}
-
+*/
//_________________________________________________________________________
//AliCaloTrackAODReader & AliCaloTrackAODReader::operator = (const AliCaloTrackAODReader & source)
//{
public:
AliCaloTrackAODReader() ; // ctor
- AliCaloTrackAODReader(const AliCaloTrackAODReader & g) ; // cpy ctor
+ //AliCaloTrackAODReader(const AliCaloTrackAODReader & g) ; // cpy ctor
//AliCaloTrackAODReader & operator = (const AliCaloTrackAODReader & g) ;//cpy assignment
virtual ~AliCaloTrackAODReader() {;} //virtual dtor
AliCaloTrackESDReader::AliCaloTrackESDReader() :
AliCaloTrackReader()
{
- //Default Ctor
-
- //Initialize parameters
- fDataType=kESD;
- fReadStack = kTRUE;
- fReadAODMCParticles = kFALSE;
- //We want tracks fitted in the detectors:
- fTrackStatus=AliESDtrack::kTPCrefit;
- fTrackStatus|=AliESDtrack::kITSrefit;
-
+ //Default Ctor
+
+ //Initialize parameters
+ fDataType=kESD;
+ fReadStack = kTRUE;
+ fReadAODMCParticles = kFALSE;
+ //We want tracks fitted in the detectors:
+ fTrackStatus=AliESDtrack::kTPCrefit;
+ fTrackStatus|=AliESDtrack::kITSrefit;
+ //This reader creates the AOD objects (Tracks-clusters), so it owns them
+ //fAODCTS ->SetOwner(kTRUE);
+ //fAODEMCAL->SetOwner(kTRUE);
+ //fAODPHOS ->SetOwner(kTRUE);
}
-
+/*
//____________________________________________________________________________
AliCaloTrackESDReader::AliCaloTrackESDReader(const AliCaloTrackESDReader & g) :
AliCaloTrackReader(g)
{
// cpy ctor
}
-
+*/
//_________________________________________________________________________
//AliCaloTrackESDReader & AliCaloTrackESDReader::operator = (const AliCaloTrackESDReader & source)
//{
Double_t bfield = ((AliESDEvent*)fInputEvent)->GetMagneticField();
Double_t timezero = 0; //TO BE FIXED
+
+ //List of output tracks
+ TClonesArray &tracks = *(fOutputEvent->GetTracks());
//To be replaced by call to AliEMCALGeoUtils when the class becomes available
Double_t radius = 441.0; //[cm] EMCAL radius +13cm
+
if(fDebug > 1) printf("AliCaloTrackESDReader::FillInputCTS() - org entries %d\n", nTracks);
for (Int_t itrack = 0; itrack < nTracks; itrack++) {////////////// track loop
AliESDtrack * track = (AliESDtrack*) ((AliESDEvent*)fInputEvent)->GetTrack(itrack) ; // retrieve track from esd
if (impactXY<3) {
// track inside the beam pipe
//Put new aod object in file in AOD tracks array
- AliAODTrack *aodTrack = new((*(fOutputEvent->GetTracks()))[naod++])
+ AliAODTrack *aodTrack = new(tracks[naod++])
+ //AliAODTrack *aodTrack = new((*(fOutputEvent->GetTracks()))[naod++])
AliAODTrack(track->GetID(), track->GetLabel(), p, kTRUE, pos, kFALSE,covTr, (Short_t)track->GetSign(), track->GetITSClusterMap(),
pid,
0x0,//primary,
public:
AliCaloTrackESDReader() ; // ctor
- AliCaloTrackESDReader(const AliCaloTrackESDReader & g) ; // cpy ctor
+ //AliCaloTrackESDReader(const AliCaloTrackESDReader & g) ; // cpy ctor
//AliCaloTrackESDReader & operator = (const AliCaloTrackESDReader & g) ;//cpy assignment
virtual ~AliCaloTrackESDReader() {;} //virtual dtor
fReadAODMCParticles = kFALSE;
}
-
+/*
//____________________________________________________________________________
AliCaloTrackMCReader::AliCaloTrackMCReader(const AliCaloTrackMCReader & g) :
AliCaloTrackReader(g), fDecayPi0(g.fDecayPi0),
{
// cpy ctor
}
-
+*/
//_________________________________________________________________________
//AliCaloTrackMCReader & AliCaloTrackMCReader::operator = (const AliCaloTrackMCReader & source)
//{
class AliCaloTrackMCReader : public AliCaloTrackReader {
public:
-
AliCaloTrackMCReader() ; // ctor
- AliCaloTrackMCReader(const AliCaloTrackMCReader & g) ; // cpy ctor
virtual ~AliCaloTrackMCReader() ;//virtual dtor
-
- private :
+ private :
+ AliCaloTrackMCReader(const AliCaloTrackMCReader & g) ; // cpy ctor
AliCaloTrackMCReader & operator = (const AliCaloTrackMCReader & g) ;//cpy assignment
public:
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliCaloTrackReader::AliCaloTrackReader(const AliCaloTrackReader & reader) :
TObject(reader), fEventNumber(reader.fEventNumber), fCurrentFileName(reader.fCurrentFileName),
{
// cpy ctor
}
-
+*/
//_________________________________________________________________________
//AliCaloTrackReader & AliCaloTrackReader::operator = (const AliCaloTrackReader & source)
//{
class AliCaloTrackReader : public TObject {
public:
-
AliCaloTrackReader() ; // ctor
- AliCaloTrackReader(const AliCaloTrackReader & g) ; // cpy ctor
virtual ~AliCaloTrackReader() ;//virtual dtor
-
private:
+ AliCaloTrackReader(const AliCaloTrackReader & g) ; // cpy ctor
AliCaloTrackReader & operator = (const AliCaloTrackReader & g) ;//cpy assignment
public:
fEMCALGeo(0x0), fPHOSGeo(0x0),
fEMCALGeoMatrixSet(kFALSE), fPHOSGeoMatrixSet(kFALSE),
fRemoveBadChannels(kFALSE),
- fEMCALBadChannelMap(new TObjArray()),fPHOSBadChannelMap(new TObjArray()),
+ fEMCALBadChannelMap(0x0),fPHOSBadChannelMap(0x0),
fNCellsFromEMCALBorder(0), fNCellsFromPHOSBorder(0),
fNoEMCALBorderAtEta0(kFALSE),fRecalibration(kFALSE),
- fEMCALRecalibrationFactors(new TObjArray()), fPHOSRecalibrationFactors(new TObjArray())
+ fEMCALRecalibrationFactors(), fPHOSRecalibrationFactors()
{
//Ctor
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliCalorimeterUtils::AliCalorimeterUtils(const AliCalorimeterUtils & calo) :
TObject(calo), fDebug(calo.fDebug),
{
// cpy ctor
}
-
+*/
//_________________________________
AliCalorimeterUtils::~AliCalorimeterUtils() {
}
if(fEMCALRecalibrationFactors) {
- fEMCALRecalibrationFactors->Clear();
- delete fEMCALRecalibrationFactors;
+ fEMCALRecalibrationFactors->Clear();
+ delete fEMCALRecalibrationFactors;
}
if(fPHOSRecalibrationFactors) {
- fPHOSRecalibrationFactors->Clear();
- delete fPHOSRecalibrationFactors;
+ fPHOSRecalibrationFactors->Clear();
+ delete fPHOSRecalibrationFactors;
}
}
if (!fRemoveBadChannels) return kFALSE;
- if(calorimeter == "EMCAL" && !fEMCALBadChannelMap->GetEntries()) return kFALSE;
- if(calorimeter == "PHOS" && !fPHOSBadChannelMap ->GetEntries()) return kFALSE;
+ if(calorimeter == "EMCAL" && !fEMCALBadChannelMap) return kFALSE;
+ if(calorimeter == "PHOS" && !fPHOSBadChannelMap) return kFALSE;
Int_t icol = -1;
Int_t irow = -1;
//In order to avoid rewriting the same histograms
Bool_t oldStatus = TH1::AddDirectoryStatus();
TH1::AddDirectory(kFALSE);
-
- //for (int i = 0; i < 12; i++) fEMCALBadChannelMap->Add(new TH2I(Form("EMCALBadChannelMap_SM%d_%s",i,fTaskName.Data()),Form("EMCALBadChannelMap_SM%d",i), 48, 0, 48, 24, 0, 24));
- for (int i = 0; i < 12; i++) fEMCALBadChannelMap->Add(new TH2I(Form("EMCALBadChannelMap_SM%d",i),Form("EMCALBadChannelMap_SM%d",i), 48, 0, 48, 24, 0, 24));
+
+ fEMCALBadChannelMap = new TObjArray(12);
+ //TH2F * hTemp = new TH2I("EMCALBadChannelMap","EMCAL SuperModule bad channel map", 48, 0, 48, 24, 0, 24);
+ for (int i = 0; i < 12; i++) {
+ fEMCALBadChannelMap->Add(new TH2I(Form("EMCALBadChannelMap_Mod%d",i),Form("EMCALBadChannelMap_Mod%d",i), 48, 0, 48, 24, 0, 24));
+ //fEMCALBadChannelMap->Add((TH2I*)hTemp->Clone(Form("EMCALBadChannelMap_Mod%d",i)));
+ }
+
+ //delete hTemp;
fEMCALBadChannelMap->SetOwner(kTRUE);
fEMCALBadChannelMap->Compress();
-
+
//In order to avoid rewriting the same histograms
TH1::AddDirectory(oldStatus);
}
//In order to avoid rewriting the same histograms
Bool_t oldStatus = TH1::AddDirectoryStatus();
TH1::AddDirectory(kFALSE);
-
- //for (int i = 0; i < 5; i++)fPHOSBadChannelMap->Add(new TH2I(Form("PHOSBadChannelMap_Mod%d_%s",i,fTaskName.Data()),Form("PHOSBadChannelMap_Mod%d",i), 56, 0, 56, 64, 0, 64));
+
+ fPHOSBadChannelMap = new TObjArray(5);
for (int i = 0; i < 5; i++)fPHOSBadChannelMap->Add(new TH2I(Form("PHOSBadChannelMap_Mod%d",i),Form("PHOSBadChannelMap_Mod%d",i), 56, 0, 56, 64, 0, 64));
fPHOSBadChannelMap->SetOwner(kTRUE);
//In order to avoid rewriting the same histograms
Bool_t oldStatus = TH1::AddDirectoryStatus();
TH1::AddDirectory(kFALSE);
-
+
+ fEMCALRecalibrationFactors = new TObjArray(12);
for (int i = 0; i < 12; i++) fEMCALRecalibrationFactors->Add(new TH2F(Form("EMCALRecalFactors_SM%d",i),Form("EMCALRecalFactors_SM%d",i), 48, 0, 48, 24, 0, 24));
//Init the histograms with 1
for (Int_t sm = 0; sm < 12; sm++) {
//In order to avoid rewriting the same histograms
Bool_t oldStatus = TH1::AddDirectoryStatus();
TH1::AddDirectory(kFALSE);
+
+ fPHOSRecalibrationFactors = new TObjArray(5);
for (int i = 0; i < 5; i++)fPHOSRecalibrationFactors->Add(new TH2F(Form("PHOSRecalFactors_Mod%d",i),Form("PHOSRecalFactors_Mod%d",i), 56, 0, 56, 64, 0, 64));
//Init the histograms with 1
for (Int_t m = 0; m < 5; m++) {
class AliCalorimeterUtils : public TObject {
-public:
-
+ public:
AliCalorimeterUtils() ; // ctor
- AliCalorimeterUtils(const AliCalorimeterUtils & g) ; // cpy ctor
virtual ~AliCalorimeterUtils() ;//virtual dtor
-
-private:
+ private:
+ AliCalorimeterUtils(const AliCalorimeterUtils & g) ; // cpy ctor
AliCalorimeterUtils & operator = (const AliCalorimeterUtils & g) ;//cpy assignment
-public:
+ public:
virtual void InitParameters();
virtual void Print(const Option_t * opt) const;
void InitPHOSRecalibrationFactors () ;
Float_t GetEMCALChannelRecalibrationFactor(Int_t iSM , Int_t iCol, Int_t iRow) const {
- if(fEMCALRecalibrationFactors->GetEntriesFast()>0) return (Float_t) ((TH2F*)fEMCALRecalibrationFactors->At(iSM))->GetBinContent(iCol,iRow);
- else return 1;}
-
+ if(fEMCALRecalibrationFactors) return (Float_t) ((TH2F*)fEMCALRecalibrationFactors->At(iSM))->GetBinContent(iCol,iRow);
+ else return 1;}
+
Float_t GetPHOSChannelRecalibrationFactor (Int_t imod, Int_t iCol, Int_t iRow) const {
- if(fPHOSRecalibrationFactors->GetEntriesFast()>0)return (Float_t) ((TH2F*)fPHOSRecalibrationFactors->At(imod))->GetBinContent(iCol,iRow);
- else return 1;}
-
+ if(fPHOSRecalibrationFactors)return (Float_t) ((TH2F*)fPHOSRecalibrationFactors->At(imod))->GetBinContent(iCol,iRow);
+ else return 1;}
+
void SetEMCALChannelRecalibrationFactor(Int_t iSM , Int_t iCol, Int_t iRow, Double_t c = 1) {
- ((TH2F*)fEMCALRecalibrationFactors->At(iSM))->SetBinContent(iCol,iRow,c);}
+ if(!fEMCALRecalibrationFactors) InitEMCALRecalibrationFactors();
+ ((TH2F*)fEMCALRecalibrationFactors->At(iSM))->SetBinContent(iCol,iRow,c);}
void SetPHOSChannelRecalibrationFactor (Int_t imod, Int_t iCol, Int_t iRow, Double_t c = 1) {
- ((TH2F*)fPHOSRecalibrationFactors->At(imod))->SetBinContent(iCol,iRow,c);}
+ if(!fPHOSRecalibrationFactors) InitPHOSRecalibrationFactors();
+ ((TH2F*)fPHOSRecalibrationFactors->At(imod))->SetBinContent(iCol,iRow,c);}
void SetEMCALChannelRecalibrationFactors(Int_t iSM , TH2F* h) {fEMCALRecalibrationFactors->AddAt(h,iSM);}
void SetPHOSChannelRecalibrationFactors(Int_t imod , TH2F* h) {fPHOSRecalibrationFactors ->AddAt(h,imod);}
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliFiducialCut::AliFiducialCut(const AliFiducialCut & g) :
TObject(g),
return *this;
}
-
+*/
//_________________________________
AliFiducialCut::~AliFiducialCut() {
//Dtor
class AliFiducialCut : public TObject {
-public:
-
+ public:
AliFiducialCut() ; // ctor
+ virtual ~AliFiducialCut() ;//virtual dtor
+
+ private:
AliFiducialCut(const AliFiducialCut & g) ; // cpy ctor
AliFiducialCut & operator = (const AliFiducialCut & g) ;//cpy assignment
- virtual ~AliFiducialCut() ;//virtual dtor
+
+ public:
void InitParameters();
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliIsolationCut::AliIsolationCut(const AliIsolationCut & g) :
TObject(g),
return *this;
}
-
+*/
//____________________________________________________________________________
TString AliIsolationCut::GetICParametersList()
{
class AliIsolationCut : public TObject {
public:
-
AliIsolationCut() ; // default ctor
+ virtual ~AliIsolationCut() {;} //virtual dtalr
+ private:
AliIsolationCut(const AliIsolationCut & g) ; // cpy ctor
AliIsolationCut & operator = (const AliIsolationCut & g) ;//cpy assignment
- virtual ~AliIsolationCut() {;} //virtual dtalr
-
+
+ public:
+
enum type {kPtThresIC, kSumPtIC, kPtFracIC, kSumPtFracIC};
Float_t GetConeSize() const {return fConeSize ; }
{
//Ctor
}
-
+/*
//____________________________________________________________________________
AliMCAnalysisUtils::AliMCAnalysisUtils(const AliMCAnalysisUtils & mcutils) :
TObject(mcutils), fCurrentEvent(mcutils.fCurrentEvent), fDebug(mcutils.fDebug),
// cpy ctor
}
+*/
//_________________________________________________________________________
//AliMCAnalysisUtils & AliMCAnalysisUtils::operator = (const AliMCAnalysisUtils & mcutils)
class AliMCAnalysisUtils : public TObject {
-public:
-
- AliMCAnalysisUtils() ; // ctor
- AliMCAnalysisUtils(const AliMCAnalysisUtils & mcu) ; // cpy ctor
- virtual ~AliMCAnalysisUtils() ;//virtual dtor
-
-private:
- AliMCAnalysisUtils & operator = (const AliMCAnalysisUtils & ) ;//cpy assignment
+ public:
+ AliMCAnalysisUtils() ; // ctor
+ virtual ~AliMCAnalysisUtils() ;//virtual dtor
+ private:
+ AliMCAnalysisUtils & operator = (const AliMCAnalysisUtils & ) ;//cpy assignment
+ AliMCAnalysisUtils(const AliMCAnalysisUtils & mcu) ; // cpy ctor
-public:
+ public:
//"Mostly" photon parent types on line 1,
//
//then "mostly" electron parent types on line 2, (e.g. electrons can
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliNeutralMesonSelection::AliNeutralMesonSelection(const AliNeutralMesonSelection & g) :
TObject(), fM(g.fM),
return *this;
}
-
+*/
//____________________________________________________________________________
AliNeutralMesonSelection::~AliNeutralMesonSelection()
{
class AliNeutralMesonSelection : public TObject {
public:
-
AliNeutralMesonSelection() ; // default ctor
+ virtual ~AliNeutralMesonSelection() ; //virtual dtor
+ private:
AliNeutralMesonSelection(const AliNeutralMesonSelection & g) ; // cpy ctor
AliNeutralMesonSelection & operator = (const AliNeutralMesonSelection & g) ;//cpy assignment
- virtual ~AliNeutralMesonSelection() ; //virtual dtor
+ public:
+
TList * GetCreateOutputObjects();
Double_t GetAngleMaxParam(Int_t i) const {return fAngleMaxParam.At(i) ; }
fHistoXBins(100), fHistoXMax(1000), fHistoXMin(-1000),
fHistoYBins(100), fHistoYMax(1000), fHistoYMin(-1000),
fHistoZBins(100), fHistoZMax(1000), fHistoZMin(-1000),
- fHistoSSBins(40), fHistoSSMax(10), fHistoSSMin(0),
+ fHistoSSBins(40), fHistoSSMax(10), fHistoSSMin(0),
fhE(0),fhPt(0),fhPhi(0),fhEta(0), fhEtaPhiE(0),
fhECharged(0),fhPtCharged(0),fhPhiCharged(0),fhEtaCharged(0), fhEtaPhiECharged(0),
//fhEChargedNoOut(0),fhPtChargedNoOut(0),fhPhiChargedNoOut(0),fhEtaChargedNoOut(0), fhEtaPhiChargedNoOut(0),
fhDeltaE(0), fhDeltaPt(0),fhDeltaPhi(0),fhDeltaEta(0), fhRatioE(0), fhRatioPt(0),fhRatioPhi(0),fhRatioEta(0),
fh2E(0),fh2Pt(0),fh2Phi(0),fh2Eta(0),
- fhLambda(0), fhDispersion(0),
- fhIM(0), fhIMCellCut(0),fhAsym(0),
+ fhLambda(0), fhDispersion(0),
+ fhIM(0), fhIMCellCut(0),fhAsym(0),
fhNCellsPerCluster(0),fhNCellsPerClusterMIP(0), fhNCellsPerClusterMIPCharged(0), fhNClusters(0),
fhClusterTimeEnergy(0),fhCellTimeSpreadRespectToCellMax(0),fhCellIdCellLargeTimeSpread(0),
fhRNCells(0),fhXNCells(0),fhYNCells(0),fhZNCells(0),
fhDeltaCellClusterRNCells(0),fhDeltaCellClusterXNCells(0),fhDeltaCellClusterYNCells(0),fhDeltaCellClusterZNCells(0),
fhDeltaCellClusterRE(0), fhDeltaCellClusterXE(0), fhDeltaCellClusterYE(0), fhDeltaCellClusterZE(0),
fhNCells(0), fhAmplitude(0), fhAmpId(0), fhEtaPhiAmp(0),
- fhTime(0), fhTimeId(0), fhTimeAmp(0), //fhT0Time(0), fhT0TimeId(0), fhT0TimeAmp(0),
+ fhTime(0), fhTimeId(0), fhTimeAmp(0), //fhT0Time(0), fhT0TimeId(0), fhT0TimeAmp(0),
fhCaloCorrNClusters(0), fhCaloCorrEClusters(0), fhCaloCorrNCells(0), fhCaloCorrECells(0),
fhEMod(0), fhNClustersMod(0), fhNCellsPerClusterMod(0), fhNCellsMod(0),
fhGridCellsMod(0), fhGridCellsEMod(0), fhGridCellsTimeMod(0),
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliAnaCalorimeterQA::AliAnaCalorimeterQA(const AliAnaCalorimeterQA & qa) :
AliAnaPartCorrBaseClass(qa), fCalorimeter(qa.fCalorimeter), fStyleMacro(qa.fStyleMacro),
// cpy ctor
}
-
+*/
//________________________________________________________________________________________________________________________________________________
//AliAnaCalorimeterQA::~AliAnaCalorimeterQA() {
// dtor
outputContainer->Add(fhEtaPhiECharged);
-// fhEChargedNoOut = new TH1F ("hEChargedNoOut","E reconstructed clusters, matched with track, no output track params", nptbins,ptmin,ptmax); \
-// fhEChargedNoOut->SetXTitle("E (GeV)"); \
-// outputContainer->Add(fhEChargedNoOut); \
-// \
-// fhPtChargedNoOut = new TH1F ("hPtChargedNoOut","p_{T} reconstructed clusters, matched with track, no output track params", nptbins,ptmin,ptmax); \
-// fhPtChargedNoOut->SetXTitle("p_{T} (GeV/c)"); \
-// outputContainer->Add(fhPtChargedNoOut); \
-// \
-// fhPhiChargedNoOut = new TH1F ("hPhiChargedNoOut","#phi reconstructed clusters, matched with track, no output track params",nphibins,phimin,phimax); \
-// fhPhiChargedNoOut->SetXTitle("#phi (rad)"); \
-// outputContainer->Add(fhPhiChargedNoOut); \
-// \
-// fhEtaChargedNoOut = new TH1F ("hEtaChargedNoOut","#eta reconstructed clusters, matched with track, no output track params",netabins,etamin,etamax); \
-// fhEtaChargedNoOut->SetXTitle("#eta "); \
-// outputContainer->Add(fhEtaChargedNoOut); \
-// \
-// fhEtaPhiChargedNoOut = new TH2F ("hEtaPhiChargedNoOut","#eta vs #phi, reconstructed clusters, matched with track, no output track params",netabins,etamin,etamax,nphibins,phimin,phimax); \
-// fhEtaPhiChargedNoOut->SetXTitle("#eta "); \
-// fhEtaPhiChargedNoOut->SetYTitle("#phi "); \
+// fhEChargedNoOut = new TH1F ("hEChargedNoOut","E reconstructed clusters, matched with track, no output track params", nptbins,ptmin,ptmax);
+// fhEChargedNoOut->SetXTitle("E (GeV)");
+// outputContainer->Add(fhEChargedNoOut);
+//
+// fhPtChargedNoOut = new TH1F ("hPtChargedNoOut","p_{T} reconstructed clusters, matched with track, no output track params", nptbins,ptmin,ptmax);
+// fhPtChargedNoOut->SetXTitle("p_{T} (GeV/c)");
+// outputContainer->Add(fhPtChargedNoOut);
+//
+// fhPhiChargedNoOut = new TH1F ("hPhiChargedNoOut","#phi reconstructed clusters, matched with track, no output track params",nphibins,phimin,phimax);
+// fhPhiChargedNoOut->SetXTitle("#phi (rad)");
+// outputContainer->Add(fhPhiChargedNoOut);
+//
+// fhEtaChargedNoOut = new TH1F ("hEtaChargedNoOut","#eta reconstructed clusters, matched with track, no output track params",netabins,etamin,etamax);
+// fhEtaChargedNoOut->SetXTitle("#eta ");
+// outputContainer->Add(fhEtaChargedNoOut);
+//
+// fhEtaPhiChargedNoOut = new TH2F ("hEtaPhiChargedNoOut","#eta vs #phi, reconstructed clusters, matched with track, no output track params",netabins,etamin,etamax,nphibins,phimin,phimax);
+// fhEtaPhiChargedNoOut->SetXTitle("#eta ");
+// fhEtaPhiChargedNoOut->SetYTitle("#phi ");
// outputContainer->Add(fhEtaPhiChargedNoOut);
fh1pOverE = new TH2F("h1pOverE","TRACK matches p/E",nptbins,ptmin,ptmax, nPoverEbins,pOverEmin,pOverEmax);
class AliAnaCalorimeterQA : public AliAnaPartCorrBaseClass {
public:
-
- AliAnaCalorimeterQA() ; // default ctor
- AliAnaCalorimeterQA(const AliAnaCalorimeterQA & g) ; // cpy ctor
-
-private:
+ AliAnaCalorimeterQA() ; // default ctor
+ virtual ~AliAnaCalorimeterQA() {;} //virtual dtor
+ private:
AliAnaCalorimeterQA & operator = (const AliAnaCalorimeterQA & g) ;//cpy assignment
-
+ AliAnaCalorimeterQA(const AliAnaCalorimeterQA & g) ; // cpy ctor
+
public:
- virtual ~AliAnaCalorimeterQA() {;} //virtual dtor
void ClusterHistograms(const TLorentzVector mom, const Double_t tof, Float_t *pos, Float_t * showerShape,
const Int_t nCaloCellsPerCluster, const Int_t nModule,
InitParameters();
}
+/*
//____________________________________________________________________________
AliAnaChargedParticles::AliAnaChargedParticles(const AliAnaChargedParticles & ch) :
AliAnaPartCorrBaseClass(ch), fPdg(ch.fPdg), fhPt(ch.fhPt), fhPhi(ch.fhPhi),fhEta(ch.fhEta),
return *this;
}
+*/
//________________________________________________________________________
TList * AliAnaChargedParticles::GetCreateOutputObjects()
class AliAnaChargedParticles : public AliAnaPartCorrBaseClass {
public:
-
AliAnaChargedParticles() ; // default ctor
+ virtual ~AliAnaChargedParticles() {;} //virtual dtor
+ private:
AliAnaChargedParticles(const AliAnaChargedParticles & g) ; // cpy ctor
AliAnaChargedParticles & operator = (const AliAnaChargedParticles & g) ;//cpy assignment
- virtual ~AliAnaChargedParticles() {;} //virtual dtor
-
+
+ public:
+
TList * GetCreateOutputObjects();
void Init();
InitParameters();\r
\r
}\r
-\r
+/*\r
//____________________________________________________________________________\r
AliAnaElectron::AliAnaElectron(const AliAnaElectron & g) \r
: AliAnaPartCorrBaseClass(g),fCalorimeter(g.fCalorimeter),\r
return *this;\r
\r
}\r
+*/\r
\r
//____________________________________________________________________________\r
AliAnaElectron::~AliAnaElectron() \r
\r
class AliAnaElectron : public AliAnaPartCorrBaseClass {\r
\r
-public: \r
-\r
+ public: \r
AliAnaElectron() ; // default ctor\r
+ virtual ~AliAnaElectron() ; //virtual dtor\r
+ private:\r
AliAnaElectron(const AliAnaElectron & g) ; // cpy ctor\r
AliAnaElectron & operator = (const AliAnaElectron & g) ;//cpy assignment\r
- virtual ~AliAnaElectron() ; //virtual dtor\r
+ \r
+ public:\r
\r
TList * GetCreateOutputObjects();\r
\r
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliAnaExample::AliAnaExample(const AliAnaExample & ex) :
AliAnaPartCorrBaseClass(ex), fPdg(ex.fPdg), fDetector(ex.fDetector), fhPt(ex.fhPt), fhPhi(ex.fhPhi),fhEta(ex.fhEta),
return *this;
}
-
+*/
// //____________________________________________________________________________
// AliAnaExample::~AliAnaExample()
// {
class AliAnaExample : public AliAnaPartCorrBaseClass {
- public:
-
+ public:
AliAnaExample() ; // default ctor
+ virtual ~AliAnaExample() {;} //virtual dtor
+
+ private :
AliAnaExample(const AliAnaExample & g) ; // cpy ctor
AliAnaExample & operator = (const AliAnaExample & g) ;//cpy assignment
- virtual ~AliAnaExample() {;} //virtual dtor
+
+ public:
TList * GetCreateOutputObjects();
//Default Ctor\r
InitParameters();\r
}\r
-\r
+/*\r
//______________________________________________________________________________\r
AliAnaOmegaToPi0Gamma::AliAnaOmegaToPi0Gamma(const AliAnaOmegaToPi0Gamma & ex) : AliAnaPartCorrBaseClass(ex),\r
fInputAODPi0(new TClonesArray (*ex.fInputAODPi0)),\r
// cpy ctor\r
//Do not need it\r
}\r
+*/\r
/*\r
//______________________________________________________________________________\r
AliAnaOmegaToPi0Gamma & AliAnaOmegaToPi0Gamma::operator = (const AliAnaOmegaToPi0Gamma & ex)\r
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliAnaParticleHadronCorrelation::AliAnaParticleHadronCorrelation(const AliAnaParticleHadronCorrelation & g):
AliAnaPartCorrBaseClass(g),
return *this;
}
-
+*/
//________________________________________________________________________
TList * AliAnaParticleHadronCorrelation::GetCreateOutputObjects()
{
class AliAnaParticleHadronCorrelation : public AliAnaPartCorrBaseClass {
public:
-
AliAnaParticleHadronCorrelation() ; // default ctor
+ virtual ~AliAnaParticleHadronCorrelation() {;} //virtual dtor
+ private:
AliAnaParticleHadronCorrelation(const AliAnaParticleHadronCorrelation & ph) ; // cpy ctor
AliAnaParticleHadronCorrelation & operator = (const AliAnaParticleHadronCorrelation & ph) ;//cpy assignment
- virtual ~AliAnaParticleHadronCorrelation() {;} //virtual dtor
-
+
+ public:
+
TList * GetCreateOutputObjects();
Double_t GetDeltaPhiMaxCut() const {return fDeltaPhiMaxCut ; }
}
-
+/*
//____________________________________________________________________________
AliAnaParticleIsolation::AliAnaParticleIsolation(const AliAnaParticleIsolation & g) :
AliAnaPartCorrBaseClass(g), fCalorimeter(g.fCalorimeter),
return *this;
}
+*/
//____________________________________________________________________________
AliAnaParticleIsolation::~AliAnaParticleIsolation()
class AliAnaParticleIsolation : public AliAnaPartCorrBaseClass {
-public:
-
+ public:
AliAnaParticleIsolation() ; // default ctor
+ virtual ~AliAnaParticleIsolation() ; //virtual dtor
+
+ private:
AliAnaParticleIsolation(const AliAnaParticleIsolation & g) ; // cpy ctor
AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & g) ;//cpy assignment
- virtual ~AliAnaParticleIsolation() ; //virtual dtor
-
+
+ public:
+
Bool_t CheckInvMass(const Int_t icalo,const AliAODPWG4Particle * ph) const ;
TList * GetCreateOutputObjects();
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliAnaParticleJetFinderCorrelation::AliAnaParticleJetFinderCorrelation(const AliAnaParticleJetFinderCorrelation & pjf) :
AliAnaPartCorrBaseClass(pjf),
return *this;
}
-
+*/
//____________________________________________________________________________
//AliAnaParticleJetFinderCorrelation::~AliAnaParticleJetFinderCorrelation()
//{
class AliAnaParticleJetFinderCorrelation : public AliAnaPartCorrBaseClass {
- public:
-
+ public:
AliAnaParticleJetFinderCorrelation() ; // default ctor
+ virtual ~AliAnaParticleJetFinderCorrelation() {;} //virtual dtor
+ private:
AliAnaParticleJetFinderCorrelation(const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor
AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ;//cpy assignment
- virtual ~AliAnaParticleJetFinderCorrelation() {;} //virtual dtor
-
+
+ public:
+
TList * GetCreateOutputObjects();
Bool_t OnlyIsolated() const {return fSelectIsolated ; }
//_________________________________________________________________________
// Class that contains the algorithm for the reconstruction of jet, cone around leading particle
// The seed is a backward particle (direct photon)
-// 1)Take the trigger particle stored in AliAODPWG4ParticleCorrelation,
+// 1) Take the trigger particle stored in AliAODPWG4ParticleCorrelation,
// 2) Search for the highest pt leading particle opposite to the photon within a phi, pt window
// 3) Take all particles around leading in a cone R with pt larger than threshold and construct the jet
//
fhNeutralLeadingPt(0),fhNeutralLeadingPhi(0),fhNeutralLeadingEta(0),
fhNeutralLeadingDeltaPt(0),fhNeutralLeadingDeltaPhi(0),fhNeutralLeadingDeltaEta(0),
fhNeutralLeadingRatioPt(0),fhChargedLeadingXi(0), fhNeutralLeadingXi(0),
- fhChargedLeadingDeltaPhiRatioPt30(0), fhNeutralLeadingDeltaPhiRatioPt30(0),
- fhChargedLeadingDeltaPhiRatioPt50(0), fhNeutralLeadingDeltaPhiRatioPt50(0),
+ fhChargedLeadingDeltaPhiRatioPt30(0), fhNeutralLeadingDeltaPhiRatioPt30(0),
+ fhChargedLeadingDeltaPhiRatioPt50(0), fhNeutralLeadingDeltaPhiRatioPt50(0),
//Jet
fhJetPt(0),fhJetRatioPt(0),fhJetDeltaPhi(0), fhJetDeltaEta(0),
fhJetLeadingRatioPt(0),fhJetLeadingDeltaPhi(0),fhJetLeadingDeltaEta(0),
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliAnaParticleJetLeadingConeCorrelation::AliAnaParticleJetLeadingConeCorrelation(const AliAnaParticleJetLeadingConeCorrelation & jetlc) :
AliAnaPartCorrBaseClass(jetlc), fJetsOnlyInCTS(jetlc.fJetsOnlyInCTS), fPbPb(jetlc.fPbPb),
return *this;
}
-
+*/
//____________________________________________________________________________
AliAnaParticleJetLeadingConeCorrelation::~AliAnaParticleJetLeadingConeCorrelation()
{
//_________________________________________________________________________
// Class that contains the algorithm for the reconstruction of jet, cone around leading particle
// The seed is a backward particle (direct photon)
-// 1)Take the a trigger particle found stored in AliAODPWG4ParticleCorrelation,
+// 1) Take the a trigger particle found stored in AliAODPWG4ParticleCorrelation,
// 2) Search for the highest pt leading particle opposite to the trigger within a phi, pt window
// 3) Take all particles around leading in a cone R with pt larger than threshold and construct the jet
//
class AliAnaParticleJetLeadingConeCorrelation : public AliAnaPartCorrBaseClass {
public:
-
AliAnaParticleJetLeadingConeCorrelation() ; // default ctor
+ virtual ~AliAnaParticleJetLeadingConeCorrelation() ; //virtual dtor
+ private:
AliAnaParticleJetLeadingConeCorrelation(const AliAnaParticleJetLeadingConeCorrelation & g) ; // cpy ctor
AliAnaParticleJetLeadingConeCorrelation & operator = (const AliAnaParticleJetLeadingConeCorrelation & g) ;//cpy assignment
- virtual ~AliAnaParticleJetLeadingConeCorrelation() ; //virtual dtor
-
+
+ public:
+
TList * GetCreateOutputObjects();
void InitParameters();
//Initialize parameters
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliAnaParticlePartonCorrelation::AliAnaParticlePartonCorrelation(const AliAnaParticlePartonCorrelation & g) :
AliAnaPartCorrBaseClass(g),
return *this;
}
-
+*/
//________________________________________________________________________
TList * AliAnaParticlePartonCorrelation::GetCreateOutputObjects()
class AliAnaParticlePartonCorrelation : public AliAnaPartCorrBaseClass {
- public:
-
- AliAnaParticlePartonCorrelation() ; // default ctor
- AliAnaParticlePartonCorrelation(const AliAnaParticlePartonCorrelation & g) ; // cpy ctor
- AliAnaParticlePartonCorrelation & operator = (const AliAnaParticlePartonCorrelation & g) ;//cpy assignment
- virtual ~AliAnaParticlePartonCorrelation() {;} //virtual dtor
-
+ public:
+ AliAnaParticlePartonCorrelation() ; // default ctor
+ virtual ~AliAnaParticlePartonCorrelation() {;} //virtual dtor
+ private:
+ AliAnaParticlePartonCorrelation(const AliAnaParticlePartonCorrelation & g) ; // cpy ctor
+ AliAnaParticlePartonCorrelation & operator = (const AliAnaParticlePartonCorrelation & g) ;//cpy assignment
+
+ public:
+
TList * GetCreateOutputObjects();
void InitParameters();
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliAnaPhoton::AliAnaPhoton(const AliAnaPhoton & g) :
AliAnaPartCorrBaseClass(g), fCalorimeter(g.fCalorimeter),
return *this;
}
-
+*/
//____________________________________________________________________________
AliAnaPhoton::~AliAnaPhoton()
{
//Check origin of the candidates
if(IsDataMC()){
- aodph.SetTag(GetMCAnalysisUtils()->CheckOrigin(calo->GetLabel(0),GetReader(), aodph.GetInputFileIndex()));
+ aodph.SetTag(GetMCAnalysisUtils()->CheckOrigin(calo->GetLabels(),calo->GetNLabel(),GetReader(), aodph.GetInputFileIndex()));
if(GetDebug() > 0) printf("AliAnaPhoton::MakeAnalysisFillAOD() - Origin of candidate, bit map %d\n",aodph.GetTag());
}//Work with stack also
class AliAnaPhoton : public AliAnaPartCorrBaseClass {
-public:
-
+ public:
AliAnaPhoton() ; // default ctor
+ virtual ~AliAnaPhoton() ; //virtual dtor
+ private:
AliAnaPhoton(const AliAnaPhoton & g) ; // cpy ctor
AliAnaPhoton & operator = (const AliAnaPhoton & g) ;//cpy assignment
- virtual ~AliAnaPhoton() ; //virtual dtor
-
+
+ public:
+
TList * GetCreateOutputObjects();
void Init();
InitParameters();
}
-
+/*
//________________________________________________________________________________________________________________________________________________
AliAnaPi0::AliAnaPi0(const AliAnaPi0 & ex) : AliAnaPartCorrBaseClass(ex),
fNCentrBin(ex.fNCentrBin),fNZvertBin(ex.fNZvertBin),fNrpBin(ex.fNrpBin),
// cpy ctor
//Do not need it
}
+*/
//
////________________________________________________________________________________________________________________________________________________
//AliAnaPi0 & AliAnaPi0::operator = (const AliAnaPi0 & ex)
class AliAnaPi0 : public AliAnaPartCorrBaseClass {
-public:
+ public:
AliAnaPi0() ; // default ctor
- AliAnaPi0(const char *name) ; // default ctor
- AliAnaPi0(const AliAnaPi0 & g) ; // cpy ctor
virtual ~AliAnaPi0() ;//virtual dtor
-
-private:
+ private:
+ AliAnaPi0(const AliAnaPi0 & g) ; // cpy ctor
AliAnaPi0 & operator = (const AliAnaPi0 & api0) ;//cpy assignment
-public:
+ public:
+
TList * GetCreateOutputObjects();
void Print(const Option_t * opt) const;
InitParameters();
}
-
+/*
//____________________________________________________________________________
AliAnaPi0EbE::AliAnaPi0EbE(const AliAnaPi0EbE & p) :
AliAnaPartCorrBaseClass(p), fAnaType(p.fAnaType), fCalorimeter(p.fCalorimeter),
return *this;
}
-
+*/
//____________________________________________________________________________
AliAnaPi0EbE::~AliAnaPi0EbE()
{
class AliAnaPi0EbE : public AliAnaPartCorrBaseClass {
public:
-
AliAnaPi0EbE() ; // default ctor
+ virtual ~AliAnaPi0EbE() ; //virtual dtor
+ private:
AliAnaPi0EbE(const AliAnaPi0EbE & g) ; // cpy ctor
AliAnaPi0EbE & operator = (const AliAnaPi0EbE & g) ;//cpy assignment
- virtual ~AliAnaPi0EbE() ; //virtual dtor
-
+
+ public:
+
enum anaTypes {kIMCalo, kSSCalo, kIMCaloTracks};
TList * GetCreateOutputObjects();
maker->SetCaloUtils(cu); //pointer to calorimeter utils
maker->AddAnalysis(emcalQA,0);
maker->AddAnalysis(phosQA,1);
- maker->SetAnaDebug(0) ; // 0 to at least print the event number
+ maker->SetAnaDebug(-1) ; // 0 to at least print the event number
maker->SwitchOnHistogramsMaker() ;
maker->SwitchOffAODsMaker() ;
if(kPrintSettings) maker->Print("");