#pragma link C++ class AliAnaScale+;
#pragma link C++ class AliNeutralMesonSelection+;
-#pragma link C++ class AliFidutialCut+;
+#pragma link C++ class AliFiducialCut+;
#pragma link C++ class AliCaloPID+;
#pragma link C++ class AliMCAnalysisUtils+;
#pragma link C++ class AliIsolationCut+;
#include "AliAnaPartCorrBaseClass.h"
#include "AliCaloTrackReader.h"
#include "AliCaloPID.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliIsolationCut.h"
#include "AliMCAnalysisUtils.h"
#include "AliNeutralMesonSelection.h"
fReader = new AliCaloTrackReader();
fCaloPID = new AliCaloPID();
- fFidCut = new AliFidutialCut();
+ fFidCut = new AliFiducialCut();
fIC = new AliIsolationCut();
fMCUtils = new AliMCAnalysisUtils();
parList+=onePar ;
sprintf(onePar,"fDataMC =%d (Check MC information, on/off) \n",fDataMC) ;
parList+=onePar ;
- sprintf(onePar,"fCheckFidCut=%d (Check Fidutial cut selection on/off) \n",fCheckFidCut) ;
+ sprintf(onePar,"fCheckFidCut=%d (Check Fiducial cut selection on/off) \n",fCheckFidCut) ;
parList+=onePar ;
sprintf(onePar,"fCheckCaloPID =%d (Use Bayesian PID in calorimetes, on/off) \n",fCheckCaloPID) ;
parList+=onePar ;
fMaxPt = 300. ; //Max pt in particle analysis
fCaloPID = new AliCaloPID ;
- fFidCut = new AliFidutialCut;
+ fFidCut = new AliFiducialCut;
fIC = new AliIsolationCut;
fNMS = new AliNeutralMesonSelection;
fNewAOD = kFALSE ;
printf("Max Photon pT = %3.2f\n", fMaxPt) ;
printf("Check PID = %d\n", fCheckCaloPID) ;
printf("Recalculate PID = %d\n", fRecalculateCaloPID) ;
- printf("Check Fidutial cut = %d\n", fCheckFidCut) ;
+ printf("Check Fiducial cut = %d\n", fCheckFidCut) ;
printf("Check MC labels = %d\n", fDataMC);
printf("Debug Level = %d\n", fDebug);
printf("Histograms: %3.1f < pT < %3.1f, Nbin = %d\n", fHistoPtMin, fHistoPtMax, fHistoNPtBins);
class AliAODCaloCells;
class AliCaloTrackReader ;
class AliCaloPID ;
-class AliFidutialCut ;
+class AliFiducialCut ;
class AliIsolationCut ;
class AliMCAnalysisUtils ;
class AliNeutralMesonSelection ;
virtual AliCaloPID * GetCaloPID() const {return fCaloPID ;}
virtual void SetCaloPID(AliCaloPID * const pid) { fCaloPID = pid ;}
- virtual AliFidutialCut * GetFidutialCut() const {return fFidCut ;}
- virtual void SetFidutialCut(AliFidutialCut * const fc) { fFidCut = fc ;}
+ virtual AliFiducialCut * GetFiducialCut() const {return fFidCut ;}
+ virtual void SetFiducialCut(AliFiducialCut * const fc) { fFidCut = fc ;}
virtual AliIsolationCut * GetIsolationCut() const {return fIC ;}
virtual void SetIsolationCut(AliIsolationCut * const ic) { fIC = ic ;}
virtual void SwitchOnDataMC() {fDataMC = kTRUE ; }
virtual void SwitchOffDataMC() {fDataMC = kFALSE ; }
- virtual Bool_t IsFidutialCutOn() const {return fCheckFidCut ; }
- virtual void SwitchOnFidutialCut() { fCheckFidCut = kTRUE;}
- virtual void SwitchOffFidutialCut() { fCheckFidCut = kFALSE;}
+ virtual Bool_t IsFiducialCutOn() const {return fCheckFidCut ; }
+ virtual void SwitchOnFiducialCut() { fCheckFidCut = kTRUE;}
+ virtual void SwitchOffFiducialCut() { fCheckFidCut = kFALSE;}
virtual Bool_t IsCaloPIDOn() const {return fCheckCaloPID ; }
virtual void SwitchOnCaloPID() { fCheckCaloPID = kTRUE;}
//Analysis helper classes access pointers
AliCaloPID * fCaloPID; // PID calculation
- AliFidutialCut * fFidCut; // Acceptance cuts
+ AliFiducialCut * fFidCut; // Acceptance cuts
AliIsolationCut * fIC; // Isolation cut
AliMCAnalysisUtils * fMCUtils; // MonteCarlo Analysis utils
AliNeutralMesonSelection * fNMS; // Neutral Meson Selection
#include "AliAODCaloCluster.h"
#include "AliAODTrack.h"
#include "AliESDtrack.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliAODInputHandler.h"
#include "AliAnalysisManager.h"
track->GetPxPyPz(p) ;
TLorentzVector momentum(p[0],p[1],p[2],0);
- if(fCTSPtMin < momentum.Pt() && fFidutialCut->IsInFidutialCut(momentum,"CTS")){
+ if(fCTSPtMin < momentum.Pt() && fFiducialCut->IsInFiducialCut(momentum,"CTS")){
if(fDebug > 2 && momentum.Pt() > 0.1) printf("AliCaloTrackAODReader::FillInputCTS() - Selected tracks E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
fAODCTS->Add(newtrack); //Use AOD stored in output for references.
}
else fAODCTS->Add(track);
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}// track loop
fAODCTSNormalInputEntries = fAODCTS->GetEntriesFast();
track->GetPxPyPz(p) ;
TLorentzVector momentum(p[0],p[1],p[2],0);
- if(fCTSPtMin < momentum.Pt() && fFidutialCut->IsInFidutialCut(momentum,"CTS")){
+ if(fCTSPtMin < momentum.Pt() && fFiducialCut->IsInFiducialCut(momentum,"CTS")){
if(fDebug > 2 && momentum.Pt() > 0.1) printf("AliCaloTrackAODReader::FillInputCTS() - Selected tracks E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
}
else fAODCTS->Add(track);
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}// track loop
if(fDebug > 1) printf("AliCaloTrackAODReader::FillInputCTS() - aod normal entries %d, after second input %d\n", fAODCTSNormalInputEntries, fAODCTS->GetEntriesFast());
GetVertex(v);
Int_t naod = (fOutputEvent->GetCaloClusters())->GetEntriesFast();
- //Loop to select clusters in fidutial cut and fill container with aodClusters
+ //Loop to select clusters in fiducial cut and fill container with aodClusters
Int_t nclusters = ((AliAODEvent*)fInputEvent)->GetNCaloClusters();
for (Int_t iclus = 0; iclus < nclusters; iclus++) {
AliAODCaloCluster * clus = 0;
TLorentzVector momentum ;
clus->GetMomentum(momentum, v);
- if(fEMCALPtMin < momentum.Pt() && fFidutialCut->IsInFidutialCut(momentum,"EMCAL")){
+ if(fEMCALPtMin < momentum.Pt() && fFiducialCut->IsInFiducialCut(momentum,"EMCAL")){
if(fDebug > 2 && momentum.E() > 0.1) printf("AliCaloTrackAODReader::FillInputEMCAL() - Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
fAODEMCAL->Add(newclus);
}
else fAODEMCAL->Add(clus);
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}//EMCAL cluster
}// cluster exists
}// cluster loop
TLorentzVector momentum ;
clus->GetMomentum(momentum, v);
- if(fEMCALPtMin < momentum.Pt() && fFidutialCut->IsInFidutialCut(momentum,"EMCAL")){
+ if(fEMCALPtMin < momentum.Pt() && fFiducialCut->IsInFiducialCut(momentum,"EMCAL")){
if(fDebug > 2 && momentum.E() > 0.1) printf("AliCaloTrackAODReader::FillInputEMCAL() - Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
fAODEMCAL->Add(newclus);
}
else fAODEMCAL->Add(clus);
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}//EMCAL cluster
}// cluster exists
}// cluster loop
GetVertex(v);
Int_t naod = (fOutputEvent->GetCaloClusters())->GetEntriesFast();
- //Loop to select clusters in fidutial cut and fill container with aodClusters
+ //Loop to select clusters in fiducial cut and fill container with aodClusters
Int_t nclusters = ((AliAODEvent*)fInputEvent)->GetNCaloClusters();
for (Int_t iclus = 0; iclus < nclusters; iclus++) {
AliAODCaloCluster * clus = 0;
TLorentzVector momentum ;
clus->GetMomentum(momentum, v);
- if(fPHOSPtMin < momentum.Pt() && fFidutialCut->IsInFidutialCut(momentum,"PHOS")){
+ if(fPHOSPtMin < momentum.Pt() && fFiducialCut->IsInFiducialCut(momentum,"PHOS")){
if(fDebug > 2 && momentum.E() > 0.1) printf("AliCaloTrackAODReader::FillInputPHOS() - Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
fAODPHOS->Add(newclus);
}
else fAODPHOS->Add(clus);
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}//PHOS cluster
}//cluster exists
}//esd cluster loop
TLorentzVector momentum ;
clus->GetMomentum(momentum, v);
- if(fPHOSPtMin < momentum.Pt() && fFidutialCut->IsInFidutialCut(momentum,"PHOS")){
+ if(fPHOSPtMin < momentum.Pt() && fFiducialCut->IsInFiducialCut(momentum,"PHOS")){
if(fDebug > 2 && momentum.E() > 0.1) printf("AliCaloTrackAODReader::FillInputPHOS() - Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
fAODPHOS->Add(newclus);
}
else fAODPHOS->Add(clus);
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}//PHOS cluster
}// cluster exists
}// cluster loop
#include "AliAODTrack.h"
#include "AliAODPid.h"
#include "AliAODEvent.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
ClassImp(AliCaloTrackESDReader)
track->GetPxPyPz(p) ;
TLorentzVector momentum(p[0],p[1],p[2],0);
- if(fCTSPtMin < momentum.Pt() &&fFidutialCut->IsInFidutialCut(momentum,"CTS")){
+ if(fCTSPtMin < momentum.Pt() &&fFiducialCut->IsInFiducialCut(momentum,"CTS")){
if(fDebug > 3 && momentum.Pt() > 0.2) printf("AliCaloTrackESDReader::FillInputCTS() - Selected tracks E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
aodTrack->SetDetPID(aodpid);
}
else continue; // outside the beam pipe: orphan track
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}// track loop
//Put references to selected tracks in array
Int_t naod = (fOutputEvent->GetCaloClusters())->GetEntriesFast();
Int_t nTracks = fInputEvent->GetNumberOfTracks() ;
- //Loop to select clusters in fidutial cut and fill container with aodClusters
+ //Loop to select clusters in fiducial cut and fill container with aodClusters
for (Int_t iclus = 0; iclus < ((AliESDEvent*)fInputEvent)->GetNumberOfCaloClusters(); iclus++) {
AliESDCaloCluster * clus = 0;
if ( (clus = ((AliESDEvent*)fInputEvent)->GetCaloCluster(iclus)) ) {
clus->GetMomentum(momentum, v);
if(fDebug > 3 && momentum.E() > 0.2) printf("AliCaloTrackESDReader::FillInputEMCAL() - all clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
- if(fEMCALPtMin < momentum.Pt() && fFidutialCut->IsInFidutialCut(momentum,"EMCAL")){
+ if(fEMCALPtMin < momentum.Pt() && fFiducialCut->IsInFiducialCut(momentum,"EMCAL")){
if(fDebug > 2 && momentum.E() > 0.2) printf("AliCaloTrackESDReader::FillInputEMCAL() - Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
}
}
//Fill reference array
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}//EMCAL cluster
}//cluster exists
}//esd cluster loop
Int_t naod = (fOutputEvent->GetCaloClusters())->GetEntriesFast();
Int_t nTracks = fInputEvent->GetNumberOfTracks() ;
- //Loop to select clusters in fidutial cut and fill container with aodClusters
+ //Loop to select clusters in fiducial cut and fill container with aodClusters
for (Int_t iclus = 0; iclus < ((AliESDEvent*)fInputEvent)->GetNumberOfCaloClusters(); iclus++) {
AliESDCaloCluster * clus = 0;
if ( (clus = ((AliESDEvent*)fInputEvent)->GetCaloCluster(iclus)) ) {
clus->GetMomentum(momentum, v);
if(fDebug > 3 && momentum.E() > 0.1)printf("AliCaloTrackESDReader::FillInputPHOS() - all clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
- if(fPHOSPtMin < momentum.Pt() && fFidutialCut->IsInFidutialCut(momentum,"PHOS")){
+ if(fPHOSPtMin < momentum.Pt() && fFiducialCut->IsInFiducialCut(momentum,"PHOS")){
if(fDebug > 2 && momentum.E() > 0.1)printf("AliCaloTrackESDReader::FillInputPHOS() - Selected clusters E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
momentum.E(),momentum.Pt(),momentum.Phi()*TMath::RadToDeg(),momentum.Eta());
caloCluster->AddTrackMatched((fInputEvent->GetTrack(iESDtrack)));
}
}
- }//Pt and Fidutial cut passed.
+ }//Pt and Fiducial cut passed.
}//PHOS cluster
}//cluster exists
}//esd cluster loop
#include "AliAODCaloCluster.h"
#include "AliAODTrack.h"
#include "AliAODEvent.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliMCAnalysisUtils.h"
ClassImp(AliCaloTrackMCReader)
//In PHOS
if(fFillPHOS && momentum.Pt() > fPHOSPtMin){
- if(!fFidutialCut->IsInFidutialCut(momentum,"PHOS")) return;
+ if(!fFiducialCut->IsInFiducialCut(momentum,"PHOS")) return;
Int_t index = iParticle ;
Int_t pdg = TMath::Abs(particle->GetPdgCode());
//In EMCAL
if(fFillEMCAL && momentum.Pt() > fEMCALPtMin){
- if(!fFidutialCut->IsInFidutialCut(momentum,"EMCAL")) return;
+ if(!fFiducialCut->IsInFiducialCut(momentum,"EMCAL")) return;
Int_t index = iParticle ;
Int_t pdg = TMath::Abs(particle->GetPdgCode());
if(fFillCTS && (momentum.Pt() > fCTSPtMin)){
//Particles in CTS acceptance
- if(!fFidutialCut->IsInFidutialCut(momentum,"CTS")) continue;
+ if(!fFiducialCut->IsInFiducialCut(momentum,"CTS")) continue;
if(fDebug > 3 && momentum.Pt() > 0.2)
printf("AliCaloTrackMCReader::FillInputEvent() - CTS : Selected tracks E %3.2f, pt %3.2f, phi %3.2f, eta %3.2f\n",
//_________________________________________________________________________
// Base class for reading data: MonteCarlo, ESD or AOD, of PHOS EMCAL and
// Central Barrel Tracking detectors (CTS).
-// Not all MC particles/tracks/clusters are kept, some kinematical/fidutial restrictions are done.
+// Not all MC particles/tracks/clusters are kept, some kinematical/fiducial restrictions are done.
// Mother class of : AliCaloTrackESDReader: Fills ESD data in 3 TObjArrays (PHOS, EMCAL, CTS)
// : AliCaloTrackMCReader: Fills Kinematics data in 3 TObjArrays (PHOS, EMCAL, CTS)
// : AliCaloTrackAODReader: Fills AOD data in 3 TObjArrays (PHOS, EMCAL, CTS)
//---- ANALYSIS system ----
#include "AliCaloTrackReader.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliMCEvent.h"
#include "AliAODMCHeader.h"
#include "AliGenPythiaEventHeader.h"
//____________________________________________________________________________
AliCaloTrackReader::AliCaloTrackReader() :
TObject(), fEventNumber(-1), fCurrentFileName(""),fDataType(0), fDebug(0),
- fFidutialCut(0x0), fComparePtHardAndJetPt(kFALSE), fPtHardAndJetPtFactor(7),
+ fFiducialCut(0x0), fComparePtHardAndJetPt(kFALSE), fPtHardAndJetPtFactor(7),
fCTSPtMin(0), fEMCALPtMin(0),fPHOSPtMin(0),
fAODCTS(new TObjArray()), fAODEMCAL(new TObjArray()), fAODPHOS(new TObjArray()),
fEMCALCells(0x0), fPHOSCells(0x0),
AliCaloTrackReader::AliCaloTrackReader(const AliCaloTrackReader & g) :
TObject(g), fEventNumber(g.fEventNumber), fCurrentFileName(g.fCurrentFileName),
fDataType(g.fDataType), fDebug(g.fDebug),
- fFidutialCut(g.fFidutialCut),
+ fFiducialCut(g.fFiducialCut),
fComparePtHardAndJetPt(g.fComparePtHardAndJetPt),
fPtHardAndJetPtFactor(g.fPtHardAndJetPtFactor),
fCTSPtMin(g.fCTSPtMin), fEMCALPtMin(g.fEMCALPtMin),fPHOSPtMin(g.fPHOSPtMin),
fDebug = source.fDebug ;
fEventNumber = source.fEventNumber ;
fCurrentFileName = source.fCurrentFileName ;
- fFidutialCut = source.fFidutialCut;
+ fFiducialCut = source.fFiducialCut;
fComparePtHardAndJetPt = source.fComparePtHardAndJetPt;
fPtHardAndJetPtFactor = source.fPtHardAndJetPtFactor;
AliCaloTrackReader::~AliCaloTrackReader() {
//Dtor
- if(fFidutialCut) delete fFidutialCut ;
+ if(fFiducialCut) delete fFiducialCut ;
if(fAODCTS){
fAODCTS->Clear() ;
fFillEMCALCells = kFALSE;
fFillPHOSCells = kFALSE;
- fFidutialCut = new AliFidutialCut();
+ fFiducialCut = new AliFiducialCut();
fSecondInputFileName = "" ;
fSecondInputFirstEvent = 0 ;
fReadStack = kFALSE; // Check in the constructor of the other readers if it was set or in the configuration file
class AliVEvent;
class AliAODEvent;
class AliMCEvent;
-class AliFidutialCut;
+class AliFiducialCut;
class AliAODMCHeader;
class AliCaloTrackReader : public TObject {
virtual void ResetLists();
- virtual AliFidutialCut * GetFidutialCut() const {return fFidutialCut ;}
- virtual void SetFidutialCut(AliFidutialCut * const fc) { fFidutialCut = fc ;}
+ virtual AliFiducialCut * GetFiducialCut() const {return fFiducialCut ;}
+ virtual void SetFiducialCut(AliFiducialCut * const fc) { fFiducialCut = fc ;}
virtual void SetInputOutputMCEvent(AliVEvent* /*esd*/, AliAODEvent* /*aod*/, AliMCEvent* /*mc*/) {;}
TString fCurrentFileName; // Current file name under analysis
Int_t fDataType ; // Select MC:Kinematics, Data:ESD/AOD, MCData:Both
Int_t fDebug; // Debugging level
- AliFidutialCut * fFidutialCut; // Acceptance cuts
+ AliFiducialCut * fFiducialCut; // Acceptance cuts
Bool_t fComparePtHardAndJetPt; // In MonteCarlo, jet events, reject fake events with wrong jet energy.
Float_t fPtHardAndJetPtFactor; // Factor between ptHard and jet pT to reject/accept event.
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-/* $Id: AliFidutialCut.cxx 21839 2007-10-29 13:49:42Z gustavo $ */
+/* $Id: AliFiducialCut.cxx 21839 2007-10-29 13:49:42Z gustavo $ */
//_________________________________________________________________________
// Class for track/cluster acceptance selection
//#include <TArrayF.h>
//---- ANALYSIS system ----
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
-ClassImp(AliFidutialCut)
+ClassImp(AliFiducialCut)
//____________________________________________________________________________
-AliFidutialCut::AliFidutialCut() :
+AliFiducialCut::AliFiducialCut() :
TObject(),
- fEMCALFidutialCut(0), fPHOSFidutialCut(0), fCTSFidutialCut(0),
+ fEMCALFiducialCut(0), fPHOSFiducialCut(0), fCTSFiducialCut(0),
fCTSFidCutMinEta(0x0),fCTSFidCutMinPhi(0x0),fCTSFidCutMaxEta(0x0), fCTSFidCutMaxPhi(0x0),
fEMCALFidCutMinEta(0x0),fEMCALFidCutMinPhi(0x0),fEMCALFidCutMaxEta(0x0), fEMCALFidCutMaxPhi(0x0),
fPHOSFidCutMinEta(0x0),fPHOSFidCutMinPhi(0x0),fPHOSFidCutMaxEta(0x0), fPHOSFidCutMaxPhi(0x0)
}
//____________________________________________________________________________
-AliFidutialCut::AliFidutialCut(const AliFidutialCut & g) :
+AliFiducialCut::AliFiducialCut(const AliFiducialCut & g) :
TObject(g),
- fEMCALFidutialCut(g.fEMCALFidutialCut), fPHOSFidutialCut(g.fPHOSFidutialCut), fCTSFidutialCut(g. fCTSFidutialCut),
+ fEMCALFiducialCut(g.fEMCALFiducialCut), fPHOSFiducialCut(g.fPHOSFiducialCut), fCTSFiducialCut(g. fCTSFiducialCut),
fCTSFidCutMinEta(g.fCTSFidCutMinEta?new TArrayF(*g.fCTSFidCutMinEta):0x0),
fCTSFidCutMinPhi(g.fCTSFidCutMinPhi?new TArrayF(*g.fCTSFidCutMinPhi):0x0),
fCTSFidCutMaxEta(g.fCTSFidCutMaxEta?new TArrayF(*g.fCTSFidCutMaxEta):0x0),
}
//_________________________________________________________________________
-AliFidutialCut & AliFidutialCut::operator = (const AliFidutialCut & source)
+AliFiducialCut & AliFiducialCut::operator = (const AliFiducialCut & source)
{
// assignment operator
if(&source == this) return *this;
- fEMCALFidutialCut = source.fEMCALFidutialCut;
- fPHOSFidutialCut = source.fPHOSFidutialCut;
- fCTSFidutialCut = source.fCTSFidutialCut;
+ fEMCALFiducialCut = source.fEMCALFiducialCut;
+ fPHOSFiducialCut = source.fPHOSFiducialCut;
+ fCTSFiducialCut = source.fCTSFiducialCut;
fCTSFidCutMinEta = source.fCTSFidCutMinEta?new TArrayF(*source.fCTSFidCutMinEta):0x0;
fCTSFidCutMinPhi = source.fCTSFidCutMinPhi?new TArrayF(*source.fCTSFidCutMinPhi):0x0;
}
//_________________________________
-AliFidutialCut::~AliFidutialCut() {
+AliFiducialCut::~AliFiducialCut() {
//Dtor
delete fCTSFidCutMinEta ;
//_______________________________________________________________
-Bool_t AliFidutialCut::IsInFidutialCut(const TLorentzVector momentum, const TString det) const
+Bool_t AliFiducialCut::IsInFiducialCut(const TLorentzVector momentum, const TString det) const
{
//Selects EMCAL or PHOS cluster or CTS track if it is inside eta-phi defined regions
if(det == "CTS"){
- if(!fCTSFidutialCut) return kTRUE; //Fidutial cut not requested, accept all tracks
- else return CheckFidutialRegion(momentum, fCTSFidCutMinPhi , fCTSFidCutMaxPhi , fCTSFidCutMinEta , fCTSFidCutMaxEta );
+ if(!fCTSFiducialCut) return kTRUE; //Fiducial cut not requested, accept all tracks
+ else return CheckFiducialRegion(momentum, fCTSFidCutMinPhi , fCTSFidCutMaxPhi , fCTSFidCutMinEta , fCTSFidCutMaxEta );
}
else if(det == "EMCAL") {
- if(!fEMCALFidutialCut) return kTRUE; //Fidutial cut not requested, accept all clusters
- else return CheckFidutialRegion(momentum, fEMCALFidCutMinPhi, fEMCALFidCutMaxPhi, fEMCALFidCutMinEta, fEMCALFidCutMaxEta);
+ if(!fEMCALFiducialCut) return kTRUE; //Fiducial cut not requested, accept all clusters
+ else return CheckFiducialRegion(momentum, fEMCALFidCutMinPhi, fEMCALFidCutMaxPhi, fEMCALFidCutMinEta, fEMCALFidCutMaxEta);
}
else if(det == "PHOS") {
- if(!fPHOSFidutialCut) return kTRUE; //Fidutial cut not requested, accept all clusters
- else return CheckFidutialRegion(momentum, fPHOSFidCutMinPhi , fPHOSFidCutMaxPhi , fPHOSFidCutMinEta , fPHOSFidCutMaxEta );
+ if(!fPHOSFiducialCut) return kTRUE; //Fiducial cut not requested, accept all clusters
+ else return CheckFiducialRegion(momentum, fPHOSFidCutMinPhi , fPHOSFidCutMaxPhi , fPHOSFidCutMinEta , fPHOSFidCutMaxEta );
}
else{
- printf("AliFidutialCut::IsInFidutialCut() - Wrong detector name = %s\n", det.Data());
+ printf("AliFiducialCut::IsInFiducialCut() - Wrong detector name = %s\n", det.Data());
return kFALSE;
}
}
//_______________________________________________________________
-Bool_t AliFidutialCut::CheckFidutialRegion(const TLorentzVector momentum, const TArrayF* minphi, const TArrayF* maxphi, const TArrayF* mineta, const TArrayF* maxeta) const {
+Bool_t AliFiducialCut::CheckFiducialRegion(const TLorentzVector momentum, const TArrayF* minphi, const TArrayF* maxphi, const TArrayF* mineta, const TArrayF* maxeta) const {
//Given the selection regions in Eta and Phi, check if particle is in this region.
Double_t phi = momentum.Phi();
if(phi < 0) phi+=TMath::TwoPi() ;
Double_t eta = momentum.Eta();
- //printf("IsInFidutialCut::Det: %s, phi = %f, eta = %f\n", det.Data(),phi*TMath::RadToDeg(), eta);
+ //printf("IsInFiducialCut::Det: %s, phi = %f, eta = %f\n", det.Data(),phi*TMath::RadToDeg(), eta);
Int_t netaregions = fCTSFidCutMaxEta->GetSize();
Int_t nphiregions = fCTSFidCutMaxPhi->GetSize();
if(netaregions != mineta->GetSize() || nphiregions != minphi->GetSize())
- printf("AliFidutialCut::IsInFidutialCut() - Wrong number of fidutial cut regions: nmaxeta %d != nmineta %d; nmaxphi %d != nminphi %d\n",
+ printf("AliFiducialCut::IsInFiducialCut() - Wrong number of fiducial cut regions: nmaxeta %d != nmineta %d; nmaxphi %d != nminphi %d\n",
netaregions, mineta->GetSize(), nphiregions, minphi->GetSize());
- //Eta fidutial cut
+ //Eta fiducial cut
Bool_t bInEtaFidCut = kFALSE;
for(Int_t ieta = 0; ieta < netaregions; ieta++)
if(eta > mineta->GetAt(ieta) && eta < maxeta->GetAt(ieta)) bInEtaFidCut = kTRUE;
if(bInEtaFidCut){
//printf("Eta cut passed\n");
- //Phi fidutial cut
+ //Phi fiducial cut
Bool_t bInPhiFidCut = kFALSE;
for(Int_t iphi = 0; iphi < nphiregions; iphi++)
if(phi > minphi->GetAt(iphi) *TMath::DegToRad()&& phi < maxphi->GetAt(iphi)*TMath::DegToRad()) bInPhiFidCut = kTRUE ;
if(bInPhiFidCut) {
- //printf("IsInFidutialCut:: %s cluster/track accepted\n",det.Data());
+ //printf("IsInFiducialCut:: %s cluster/track accepted\n",det.Data());
return kTRUE;
}
else return kFALSE;
//_______________________________________________________________
-void AliFidutialCut::InitParameters()
+void AliFiducialCut::InitParameters()
{
//Initialize the parameters of the analysis.
- fEMCALFidutialCut = kTRUE ;
- fPHOSFidutialCut = kTRUE ;
- fCTSFidutialCut = kTRUE ;
+ fEMCALFiducialCut = kTRUE ;
+ fPHOSFiducialCut = kTRUE ;
+ fCTSFiducialCut = kTRUE ;
fCTSFidCutMinEta = new TArrayF(1);
fCTSFidCutMinEta->SetAt(-0.9,0);
//________________________________________________________________
-void AliFidutialCut::Print(const Option_t * opt) const
+void AliFiducialCut::Print(const Option_t * opt) const
{
//Print some relevant parameters set for the analysis
printf("***** Print: %s %s ******\n", GetName(), GetTitle() ) ;
- if(fCTSFidutialCut){
+ if(fCTSFiducialCut){
Int_t netaregions = fCTSFidCutMaxEta->GetSize();
Int_t nphiregions = fCTSFidCutMaxPhi->GetSize();
- printf(">> CTS Fidutial regions : phi %d eta %d\n", netaregions, nphiregions) ;
+ printf(">> CTS Fiducial regions : phi %d eta %d\n", netaregions, nphiregions) ;
for(Int_t ieta = 0; ieta < netaregions; ieta++)
printf(" region %d : %3.2f < eta < %3.2f\n", ieta, fCTSFidCutMinEta->GetAt(ieta), fCTSFidCutMaxEta->GetAt(ieta)) ;
for(Int_t iphi = 0; iphi < nphiregions; iphi++)
printf(" region %d : %3.1f < phi < %3.1f\n", iphi, fCTSFidCutMinPhi->GetAt(iphi), fCTSFidCutMaxPhi->GetAt(iphi)) ;
}
- else printf(">>No fidutial cuts in CTS\n");
+ else printf(">>No fiducial cuts in CTS\n");
- if(fEMCALFidutialCut){
+ if(fEMCALFiducialCut){
Int_t netaregions = fEMCALFidCutMaxEta->GetSize();
Int_t nphiregions = fEMCALFidCutMaxPhi->GetSize();
- printf(">>EMCAL Fidutial regions : phi %d eta %d\n", netaregions, nphiregions) ;
+ printf(">>EMCAL Fiducial regions : phi %d eta %d\n", netaregions, nphiregions) ;
for(Int_t ieta = 0; ieta < netaregions; ieta++)
printf(" region %d : %3.2f < eta < %3.2f\n", ieta, fEMCALFidCutMinEta->GetAt(ieta), fEMCALFidCutMaxEta->GetAt(ieta)) ;
for(Int_t iphi = 0; iphi < nphiregions; iphi++)
printf(" region %d : %3.1f < phi < %3.1f\n", iphi, fEMCALFidCutMinPhi->GetAt(iphi), fEMCALFidCutMaxPhi->GetAt(iphi)) ;
}
- else printf(">>No fidutial cuts in EMCAL\n");
+ else printf(">>No fiducial cuts in EMCAL\n");
- if(fPHOSFidutialCut){
+ if(fPHOSFiducialCut){
Int_t netaregions = fPHOSFidCutMaxEta->GetSize();
Int_t nphiregions = fPHOSFidCutMaxPhi->GetSize();
- printf(">>PHOS Fidutial regions : phi %d eta %d\n", netaregions, nphiregions) ;
+ printf(">>PHOS Fiducial regions : phi %d eta %d\n", netaregions, nphiregions) ;
for(Int_t ieta = 0; ieta < netaregions; ieta++)
printf(" region %d : %3.2f < eta < %3.2f\n", ieta, fPHOSFidCutMinEta->GetAt(ieta), fPHOSFidCutMaxEta->GetAt(ieta)) ;
for(Int_t iphi = 0; iphi < nphiregions; iphi++)
printf(" region %d : %3.1f < phi < %3.1f\n", iphi, fPHOSFidCutMinPhi->GetAt(iphi), fPHOSFidCutMaxPhi->GetAt(iphi)) ;
}
- else printf(">>No fidutial cuts in PHOS\n");
+ else printf(">>No fiducial cuts in PHOS\n");
printf(" \n") ;
}
//_______________________________________________________________
-void AliFidutialCut::SetSimpleCTSFidutialCut(const Float_t eta, const Float_t minphi, const Float_t maxphi){
+void AliFiducialCut::SetSimpleCTSFiducialCut(const Float_t eta, const Float_t minphi, const Float_t maxphi){
//Method to set simple acceptance cut to CTS
fCTSFidCutMinEta->Set(1);
//_______________________________________________________________
-void AliFidutialCut::SetSimpleEMCALFidutialCut(const Float_t eta, const Float_t minphi, const Float_t maxphi){
+void AliFiducialCut::SetSimpleEMCALFiducialCut(const Float_t eta, const Float_t minphi, const Float_t maxphi){
//Method to set simple acceptance cut to EMCAL
fEMCALFidCutMinEta->Set(1);
}
//_______________________________________________________________
-void AliFidutialCut::SetSimplePHOSFidutialCut(const Float_t eta, const Float_t minphi, const Float_t maxphi){
+void AliFiducialCut::SetSimplePHOSFiducialCut(const Float_t eta, const Float_t minphi, const Float_t maxphi){
//Method to set simple acceptance cut to PHOS
fPHOSFidCutMinEta->Set(1);
-#ifndef ALIFIDUTIALCUT_H
-#define ALIFIDUTIALCUT_H
+#ifndef ALIFIDUCIALCUT_H
+#define ALIFIDUCIALCUT_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
/* $Id: $ */
//--- AliRoot system ---
-class AliFidutialCut : public TObject {
+class AliFiducialCut : public TObject {
public:
- AliFidutialCut() ; // ctor
- AliFidutialCut(const AliFidutialCut & g) ; // cpy ctor
- AliFidutialCut & operator = (const AliFidutialCut & g) ;//cpy assignment
- virtual ~AliFidutialCut() ;//virtual dtor
+ AliFiducialCut() ; // ctor
+ AliFiducialCut(const AliFiducialCut & g) ; // cpy ctor
+ AliFiducialCut & operator = (const AliFiducialCut & g) ;//cpy assignment
+ virtual ~AliFiducialCut() ;//virtual dtor
void InitParameters();
void Print(const Option_t * opt)const;
- Bool_t CheckFidutialRegion(const TLorentzVector lv, const TArrayF* minphi, const TArrayF* maxphi, const TArrayF* mineta, const TArrayF* maxeta) const ;
- Bool_t IsInFidutialCut (const TLorentzVector lv, const TString det) const ;
+ Bool_t CheckFiducialRegion(const TLorentzVector lv, const TArrayF* minphi, const TArrayF* maxphi, const TArrayF* mineta, const TArrayF* maxeta) const ;
+ Bool_t IsInFiducialCut (const TLorentzVector lv, const TString det) const ;
- void DoCTSFidutialCut(Bool_t b) {fCTSFidutialCut = b; }
- void DoEMCALFidutialCut(Bool_t b) {fEMCALFidutialCut = b; }
- void DoPHOSFidutialCut(Bool_t b) {fPHOSFidutialCut = b; }
+ void DoCTSFiducialCut(Bool_t b) {fCTSFiducialCut = b; }
+ void DoEMCALFiducialCut(Bool_t b) {fEMCALFiducialCut = b; }
+ void DoPHOSFiducialCut(Bool_t b) {fPHOSFiducialCut = b; }
- Bool_t GetCTSFidutialCutStatus() const {return fCTSFidutialCut ; }
- Bool_t GetEMCALFidutialCut() const {return fEMCALFidutialCut ; }
- Bool_t GetPHOSFidutialCutStatus() const {return fPHOSFidutialCut ; }
+ Bool_t GetCTSFiducialCutStatus() const {return fCTSFiducialCut ; }
+ Bool_t GetEMCALFiducialCut() const {return fEMCALFiducialCut ; }
+ Bool_t GetPHOSFiducialCutStatus() const {return fPHOSFiducialCut ; }
- void SetSimpleCTSFidutialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
- void SetSimpleEMCALFidutialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
- void SetSimplePHOSFidutialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
+ void SetSimpleCTSFiducialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
+ void SetSimpleEMCALFiducialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
+ void SetSimplePHOSFiducialCut(const Float_t abseta, const Float_t phimin, const Float_t phimax) ;
void AddCTSFidCutMaxEtaArray(TArrayF & array)
{ fCTSFidCutMaxEta = new TArrayF(array) ; }
protected:
//Detector acceptance cuts
- Bool_t fEMCALFidutialCut ; // Apply fidutial cuts to EMCAL clusters
- Bool_t fPHOSFidutialCut ;// Apply fidutial cuts to PHOS clusters
- Bool_t fCTSFidutialCut ;//Apply fidutial cuts to CTS tracks
+ Bool_t fEMCALFiducialCut ; // Apply fiducial cuts to EMCAL clusters
+ Bool_t fPHOSFiducialCut ;// Apply fiducial cuts to PHOS clusters
+ Bool_t fCTSFiducialCut ;//Apply fiducial cuts to CTS tracks
TArrayF * fCTSFidCutMinEta ; //Take particles in CTS with eta > fCTSFidCutMinEta
TArrayF * fCTSFidCutMinPhi ; //Take particles in CTS with phi > fCTSFidCutMinPhi
TArrayF * fPHOSFidCutMaxEta ; //Take particles in PHOS with eta < fPHOSFidCutMaxEta
TArrayF * fPHOSFidCutMaxPhi ; //Take particles in PHOS with phi > fPHOSFidCutMaxPhi
- ClassDef(AliFidutialCut,1)
+ ClassDef(AliFiducialCut,1)
} ;
-#endif //ALIFIDUTIALCUT_H
+#endif //ALIFIDUCIALCUT_H
#include "AliAODCaloCells.h"
#include "AliESDCaloCells.h"
#include "AliAODCaloCluster.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliESDtrack.h"
#include "AliAODTrack.h"
#include "AliESDCaloCluster.h"
//Get cluster kinematics
clus2->GetMomentum(mom2,v);
}
-
+ if(GetEventNumber() > 139 && GetEventNumber() < 145)
+ printf("Event %d: epair %f mass %f\n", GetEventNumber(), (mom+mom2).E(),(mom+mom2).M());
fhIM ->Fill((mom+mom2).E(),(mom+mom2).M());
fhAsym->Fill((mom+mom2).E(),TMath::Abs((mom.E()-mom2.E())/(mom.E()+mom2.E())));
if (TMath::Abs(etaMC) > 1) return;
Bool_t in = kTRUE;
- if(IsFidutialCutOn()) in = GetFidutialCut()->IsInFidutialCut(mom,fCalorimeter) ;
+ if(IsFiducialCutOn()) in = GetFiducialCut()->IsInFiducialCut(mom,fCalorimeter) ;
if (pdg==22) {
fhGenGamPt ->Fill(ptMC);
c5->cd(1) ;
fhIM->SetLineColor(4);
- fhIM->Draw();
+ fhIM->Draw("lego");
c5->cd(2) ;
binmin = 0;
- binmax = (Int_t) (5-emin)*nebins/emax;
+ binmax = (Int_t) (1-emin)*nebins/emax;
TH1D *pyim5 = fhIM->ProjectionY("pyim5",binmin,binmax);
- pyim5->SetTitle("E_{pair} < 5 GeV");
+ pyim5->SetTitle("0 < E_{pair} < 1 GeV");
+ if(pyim5->GetEntries() > 1) gPad->SetLogy();
+ pyim5->SetTitle("All Pairs");
pyim5->SetLineColor(4);
pyim5->Draw();
c5->cd(3) ;
- binmin = (Int_t) (5-emin)*nebins/emax;
- binmax = (Int_t) (10-emin)*nebins/emax;
- TH1D *pyim510 = fhIM->ProjectionY("pyim5_10",binmin,binmax);
- pyim510->SetTitle("5 < E_{pair} < 10 GeV");
+ //binmin = (Int_t) (5-emin)*nebins/emax;
+ //binmax = (Int_t) (10-emin)*nebins/emax;
+ //TH1D *pyim510 = fhIM->ProjectionY("pyim5_10",binmin,binmax);
+ //pyim510->SetTitle("5 < E_{pair} < 10 GeV");
+ binmin = (Int_t) (1-emin)*nebins/emax;
+ binmax = (Int_t) (2-emin)*nebins/emax;
+ TH1D *pyim510 = fhIM->ProjectionY("pyim0_2",binmin,binmax);
+ if(pyim510->GetEntries() > 1) gPad->SetLogy();
+ pyim510->SetTitle("1 < E_{pair} < 2 GeV");
+
pyim510->SetLineColor(4);
pyim510->Draw();
c5->cd(4) ;
- binmin = (Int_t) (10-emin)*nebins/emax;
+ binmin = (Int_t) (2-emin)*nebins/emax;
binmax = -1;
TH1D *pyim10 = fhIM->ProjectionY("pyim10",binmin,binmax);
- pyim10->SetTitle("E_{pair} > 10 GeV");
+ pyim10->SetTitle("E_{pair} > 2 GeV");
+ if(pyim10->GetEntries() > 1) gPad->SetLogy();
pyim10->SetLineColor(4);
pyim10->Draw();
TString GetStyleMacro() const {return fStyleMacro ;}
void SetStyleMacro( TString macro ) {fStyleMacro = macro; }
- void SwithOnPlotsMaking() {fMakePlots = kTRUE;}
- void SwithOffPlotsMaking() {fMakePlots = kFALSE;}
+ void SwitchOnPlotsMaking() {fMakePlots = kTRUE;}
+ void SwitchOffPlotsMaking() {fMakePlots = kFALSE;}
void Terminate(TList * outputList);
void ReadHistograms(TList * outputList); //Fill histograms with histograms in ouput list, needed in Terminate.
#include "AliCaloTrackReader.h"
#include "AliAODPWG4Particle.h"
#include "AliStack.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliAODTrack.h"
#include "AliAODMCParticle.h"
p3.SetXYZ(mom[0],mom[1],mom[2]);
//Acceptance selection
- Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,"CTS") ;
- if(GetDebug() > 1) printf("AliAnaChargedParticles::MakeAnalysisFillAOD() - Track pt %2.2f, phi %2.2f, in fidutial cut %d\n",p3.Pt(), p3.Phi(), in);
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,"CTS") ;
+ if(GetDebug() > 1) printf("AliAnaChargedParticles::MakeAnalysisFillAOD() - Track pt %2.2f, phi %2.2f, in fiducial cut %d\n",p3.Pt(), p3.Phi(), in);
if(p3.Pt() > GetMinPt() && in) {
//Keep only particles identified with fPdg
//Selection not done for the moment
#include "AliCaloTrackReader.h"\r
#include "AliMCAnalysisUtils.h"\r
#include "AliAODCaloCluster.h"\r
-#include "AliFidutialCut.h"\r
+#include "AliFiducialCut.h"\r
#include "AliAODTrack.h"\r
#include "AliAODPid.h"\r
#include "AliCaloPID.h"\r
//Get parameters set in base class.\r
parList += GetBaseParametersList() ;\r
\r
- //Get parameters set in FidutialCut class (not available yet)\r
+ //Get parameters set in FiducialCut class (not available yet)\r
//parlist += GetFidCut()->GetFidCutParametersList() \r
\r
TObjString *oString= new TObjString(parList) ;\r
track->Eta() > -0.7 && track->Eta() < 0.7) in = kTRUE;\r
////////////////////////////\r
//THIS HAS A MEM LEAK JLK 24-Oct-09\r
- //Bool_t in = GetFidutialCut()->IsInFidutialCut(mom2,fCalorimeter) ;\r
+ //Bool_t in = GetFiducialCut()->IsInFiducialCut(mom2,fCalorimeter) ;\r
///////////////////////////\r
- if(GetDebug() > 1) printf("AliAnaElectron::MakeAnalysisFillAOD() - Track(Extrap) pt %2.2f(%2.2f), phi %2.2f(%2.2f), eta %2.2f(%2.2f) in fidutial cut %d\n",track->Pt(), mom.Pt(), track->Phi(), mom.Phi(), track->Eta(),mom.Eta(), in);\r
+ if(GetDebug() > 1) printf("AliAnaElectron::MakeAnalysisFillAOD() - Track(Extrap) pt %2.2f(%2.2f), phi %2.2f(%2.2f), eta %2.2f(%2.2f) in fiducial cut %d\n",track->Pt(), mom.Pt(), track->Phi(), mom.Phi(), track->Eta(),mom.Eta(), in);\r
\r
if(mom.Pt() > GetMinPt() && in) {\r
\r
mom.Eta() > -0.7 && mom.Eta() < 0.7) in = kTRUE;\r
/////////////////////////////////\r
//THIS HAS A MEM LEAK JLK 24-Oct-09\r
- //Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,fCalorimeter);\r
+ //Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,fCalorimeter);\r
////////////////////////////////\r
if(mom.Pt() < GetMinPt()) continue;\r
if(!in) continue;\r
#include "AliESDCaloCells.h"
#include "AliStack.h"
#include "AliCaloPID.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliAODCaloCells.h"
#include "AliAODCaloCluster.h"
#include "AliAODTrack.h"
//Acceptance selection
Bool_t in = kTRUE;
- if(IsFidutialCutOn())
- in = GetFidutialCut()->IsInFidutialCut(mom,fDetector) ;
+ if(IsFiducialCutOn())
+ in = GetFiducialCut()->IsInFiducialCut(mom,fDetector) ;
- if(GetDebug() > 1) printf("AliAnaExample::MakeAnalysisFillAOD() - Cluster pt %2.2f, phi %2.2f, pdg %d, in fidutial cut %d \n",mom.Pt(), mom.Phi(), pdg, in);
+ if(GetDebug() > 1) printf("AliAnaExample::MakeAnalysisFillAOD() - Cluster pt %2.2f, phi %2.2f, pdg %d, in fiducial cut %d \n",mom.Pt(), mom.Phi(), pdg, in);
//Select cluster if momentum, pdg and acceptance are good
if(mom.Pt() > GetMinPt() && pdg ==fPdg && in) {
p3.SetXYZ(mom[0],mom[1],mom[2]);
//Acceptance selection
- Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,"CTS") ;
- if(GetDebug() > 1) printf("AliAnaExample::MakeAnalysisFillAOD() - Track pt %2.2f, phi %2.2f, in fidutial cut %d\n",p3.Pt(), p3.Phi(), in);
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,"CTS") ;
+ if(GetDebug() > 1) printf("AliAnaExample::MakeAnalysisFillAOD() - Track pt %2.2f, phi %2.2f, in fiducial cut %d\n",p3.Pt(), p3.Phi(), in);
if(p3.Pt() > GetMinPt() && in) {
AliAODPWG4Particle tr = AliAODPWG4Particle(mom[0],mom[1],mom[2],0);
tr.SetDetector("CTS");
#include "AliCaloTrackReader.h"
#include "AliCaloPID.h"
#include "AliAODPWG4ParticleCorrelation.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliAODTrack.h"
#include "AliAODCaloCluster.h"
#include "AliMCAnalysisUtils.h"
// printf("rat = %f, xE = %f, cosi =%f \n", rat, xE, cosi);
// printf("phi = %f \n", phi);
- if(IsFidutialCutOn()){
- Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,"CTS") ;
+ if(IsFiducialCutOn()){
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,"CTS") ;
if(! in ) continue ;
}
if (aodParticle->GetDetector() == "EMCAL" && GetReader()->GetAODEMCALNormalInputEntries() <= iclus) inputi = 1 ;
else if(aodParticle->GetDetector() == "PHOS" && GetReader()->GetAODPHOSNormalInputEntries() <= iclus) inputi = 1;
- //Cluster selection, not charged, with photon or pi0 id and in fidutial cut
+ //Cluster selection, not charged, with photon or pi0 id and in fiducial cut
Int_t pdg=0;
if (inputi == 0 && !SelectCluster(calo, vertex, gammai, pdg)) continue ;
else if(inputi == 1 && !SelectCluster(calo, vertex2, gammai, pdg)) continue ;
if (aodParticle->GetDetector() == "EMCAL" && GetReader()->GetAODEMCALNormalInputEntries() <= jclus) inputj = 1;
else if(aodParticle->GetDetector() == "PHOS" && GetReader()->GetAODPHOSNormalInputEntries() <= jclus) inputj = 1;
- //Cluster selection, not charged with photon or pi0 id and in fidutial cut
+ //Cluster selection, not charged with photon or pi0 id and in fiducial cut
Int_t pdgj=0;
if (inputj == 0 && !SelectCluster(calo2, vertex, gammaj, pdgj)) continue ;
else if(inputj == 1 && !SelectCluster(calo2, vertex2, gammaj, pdgj)) continue ;
}//PID on
//Check acceptance selection
- if(IsFidutialCutOn()){
- Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,detector) ;
+ if(IsFiducialCutOn()){
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,detector) ;
if(! in ) return kFALSE ;
}
#include "AliAnaParticleJetLeadingConeCorrelation.h"
#include "AliCaloPID.h"
#include "AliAODPWG4ParticleCorrelation.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliAODTrack.h"
#include "AliAODCaloCluster.h"
if (particle->GetDetector() == "EMCAL" && GetReader()->GetAODEMCALNormalInputEntries() <= iclus) inputi = 1 ;
else if(particle->GetDetector() == "PHOS" && GetReader()->GetAODPHOSNormalInputEntries() <= iclus) inputi = 1;
- //Cluster selection, not charged, with photon or pi0 id and in fidutial cut
+ //Cluster selection, not charged, with photon or pi0 id and in fiducial cut
Int_t pdgi=0;
if (inputi == 0 && !SelectCluster(calo, vertex, gammai, pdgi)) continue ;
else if(inputi == 1 && !SelectCluster(calo, vertex2, gammai, pdgi)) continue ;
if (particle->GetDetector() == "EMCAL" && GetReader()->GetAODEMCALNormalInputEntries() <= jclus) inputj = 1;
else if(particle->GetDetector() == "PHOS" && GetReader()->GetAODPHOSNormalInputEntries() <= jclus) inputj = 1;
- //Cluster selection, not charged with photon or pi0 id and in fidutial cut
+ //Cluster selection, not charged with photon or pi0 id and in fiducial cut
Int_t pdgj=0;
if (inputj == 0 && !SelectCluster(calo2, vertex, gammaj, pdgj)) continue ;
else if(inputj == 1 && !SelectCluster(calo2, vertex2, gammaj, pdgj)) continue ;
}//CaloPID
//Check acceptance selection
- if(IsFidutialCutOn()){
- Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,"EMCAL") ;
+ if(IsFiducialCutOn()){
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,"EMCAL") ;
if(! in ) return kFALSE ;
}
#include "AliStack.h"
#include "AliCaloPID.h"
#include "AliMCAnalysisUtils.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "AliAODCaloCluster.h"
#include "AliAODMCParticle.h"
//Get parameters set in PID class.
parList += GetCaloPID()->GetPIDParametersList() ;
- //Get parameters set in FidutialCut class (not available yet)
+ //Get parameters set in FiducialCut class (not available yet)
//parlist += GetFidCut()->GetFidCutParametersList()
TObjString *oString= new TObjString(parList) ;
for(Int_t icalo = 0; icalo < pl->GetEntriesFast(); icalo++){
AliAODCaloCluster * calo = (AliAODCaloCluster*) (pl->At(icalo));
- //Cluster selection, not charged, with photon id and in fidutial cut
+ //Cluster selection, not charged, with photon id and in fiducial cut
//Input from second AOD?
Int_t input = 0;
//printf("AliAnaPhoton::Current Event %d; Current File Name : %s, E %f, pT %f, Ecl %f\n",GetReader()->GetEventNumber(),(GetReader()->GetCurrentFileName()).Data(), mom.E(), mom.Pt(),calo->E());
//Check acceptance selection
- if(IsFidutialCutOn()){
- Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,fCalorimeter) ;
+ if(IsFiducialCutOn()){
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,fCalorimeter) ;
if(! in ) continue ;
}
aodph.SetCaloLabel(calo->GetID(),-1);
aodph.SetDetector(fCalorimeter);
if(GetDebug() > 1)
- printf("AliAnaPhoton::MakeAnalysisFillAOD() - Min pt cut and fidutial cut passed: pt %3.2f, phi %2.2f, eta %1.2f\n",aodph.Pt(),aodph.Phi(),aodph.Eta());
+ printf("AliAnaPhoton::MakeAnalysisFillAOD() - Min pt cut and fiducial cut passed: pt %3.2f, phi %2.2f, eta %1.2f\n",aodph.Pt(),aodph.Phi(),aodph.Eta());
//Check Distance to Bad channel, set bit.
Double_t distBad=calo->GetDistToBadChannel() ; //Distance to bad channel
#include "AliCaloTrackReader.h"
#include "AliCaloPID.h"
#include "AliStack.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "TParticle.h"
#include "AliAODCaloCluster.h"
#include "AliVEvent.h"
TLorentzVector lv1, lv2;
phot1->Momentum(lv1);
phot2->Momentum(lv2);
- if(GetFidutialCut()->IsInFidutialCut(lv1,fCalorimeter) && GetFidutialCut()->IsInFidutialCut(lv2,fCalorimeter))
+ if(GetFiducialCut()->IsInFiducialCut(lv1,fCalorimeter) && GetFiducialCut()->IsInFiducialCut(lv2,fCalorimeter))
inacceptance = kTRUE ;
if(GetDebug() > 2) printf("In %s fiducial cut acceptance? %d\n",fCalorimeter.Data(),inacceptance);
}
TLorentzVector lv1, lv2;
phot1->Momentum(lv1);
phot2->Momentum(lv2);
- if(GetFidutialCut()->IsInFidutialCut(lv1,fCalorimeter) && GetFidutialCut()->IsInFidutialCut(lv2,fCalorimeter))
+ if(GetFiducialCut()->IsInFiducialCut(lv1,fCalorimeter) && GetFiducialCut()->IsInFiducialCut(lv2,fCalorimeter))
inacceptance = kTRUE ;
if(GetDebug() > 2) printf("In %s fiducial cut acceptance? %d\n",fCalorimeter.Data(),inacceptance);
}
#include "AliCaloPID.h"
#include "AliMCAnalysisUtils.h"
#include "AliStack.h"
-#include "AliFidutialCut.h"
+#include "AliFiducialCut.h"
#include "TParticle.h"
#include "AliAODCaloCluster.h"
#include "AliAODEvent.h"
for(Int_t icalo = 0; icalo < pl->GetEntriesFast(); icalo++){
AliAODCaloCluster * calo = (AliAODCaloCluster*) (pl->At(icalo));
- //Cluster selection, not charged, with pi0 id and in fidutial cut
+ //Cluster selection, not charged, with pi0 id and in fiducial cut
//Input from second AOD?
Int_t input = 0;
//If too small or big pt, skip it
if(mom.Pt() < GetMinPt() || mom.Pt() > GetMaxPt() ) continue ;
//Check acceptance selection
- if(IsFidutialCutOn()){
- Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,fCalorimeter) ;
+ if(IsFiducialCutOn()){
+ Bool_t in = GetFiducialCut()->IsInFiducialCut(mom,fCalorimeter) ;
if(! in ) continue ;
}
aodpi0.SetCaloLabel(calo->GetID(),-1);
aodpi0.SetDetector(fCalorimeter);
if(GetDebug() > 1)
- printf("AliAnaPi0EbE::MakeShowerShapeIdentification() - FillAOD: Min pt cut and fidutial cut passed: pt %3.2f, phi %2.2f, eta %1.2f\n",aodpi0.Pt(),aodpi0.Phi(),aodpi0.Eta());
+ printf("AliAnaPi0EbE::MakeShowerShapeIdentification() - FillAOD: Min pt cut and fiducial cut passed: pt %3.2f, phi %2.2f, eta %1.2f\n",aodpi0.Pt(),aodpi0.Phi(),aodpi0.Eta());
//Check Distance to Bad channel, set bit.
Double_t distBad=calo->GetDistToBadChannel() ; //Distance to bad channel
#-*- Mode: Makefile -*-
-SRCS = PartCorrBase/AliNeutralMesonSelection.cxx PartCorrBase/AliFidutialCut.cxx \
+SRCS = PartCorrBase/AliNeutralMesonSelection.cxx PartCorrBase/AliFiducialCut.cxx \
PartCorrBase/AliCaloPID.cxx PartCorrBase/AliMCAnalysisUtils.cxx PartCorrBase/AliIsolationCut.cxx\
PartCorrBase/AliAnaScale.cxx PartCorrBase/AliAnaPartCorrMaker.cxx \
PartCorrBase/AliAnaPartCorrBaseClass.cxx PartCorrBase/AliAnalysisTaskParticleCorrelation.cxx \
// ##### Analysis algorithm settings ####
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
AliAnaCalorimeterQA *emcalQA = new AliAnaCalorimeterQA();
//emcalQA->SetDebug(2); //10 for lots of messages
if(kUseKinematics) emcalQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.
else emcalQA->SwitchOffDataMC() ;
emcalQA->AddToHistogramsName("EMCAL_"); //Begining of histograms name
- emcalQA->SetFidutialCut(fidCut);
- emcalQA->SwitchOnFidutialCut();
+ emcalQA->SetFiducialCut(fidCut);
+ emcalQA->SwitchOnFiducialCut();
if(kPrintSettings) emcalQA->Print("");
//emcalQA->GetMCAnalysisUtils()->SetDebug(10);
if(kUseKinematics) phosQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.
else phosQA->SwitchOffDataMC() ;
phosQA->AddToHistogramsName("PHOS_");//Begining of histograms name
- phosQA->SetFidutialCut(fidCut);
- phosQA->SwitchOnFidutialCut();
+ phosQA->SetFiducialCut(fidCut);
+ phosQA->SwitchOnFiducialCut();
if(kPrintSettings)phosQA->Print("");
//phosQA->GetMCAnalysisUtils()->SetDebug(10);
pid->SetDebug(-1);
if(kPrintSettings) pid->Print("");
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
AliAnaCalorimeterQA *qa = new AliAnaCalorimeterQA();
qa->SetDebug(-1); //10 for lots of messages
else qa->SwitchOffDataMC() ;
//qa->AddToHistogramsName(Form("AnaCaloQA_%s",calorimeter.Data()));
qa->AddToHistogramsName("AnaCaloQA_");
- qa->SetFidutialCut(fidCut);
- qa->SwitchOnFidutialCut();
+ qa->SetFiducialCut(fidCut);
+ qa->SwitchOnFiducialCut();
if(kPrintSettings) qa->Print("");
- AliFidutialCut * fidCut1stYear = new AliFidutialCut();
- fidCut1stYear->DoCTSFidutialCut(kFALSE) ;
- fidCut1stYear->DoEMCALFidutialCut(kTRUE) ;
- fidCut1stYear->DoPHOSFidutialCut(kTRUE) ;
- fidCut1stYear->SetSimpleEMCALFidutialCut(0.7,80.,120.);
- fidCut1stYear->SetSimplePHOSFidutialCut(0.12,260.,320.);
+ AliFiducialCut * fidCut1stYear = new AliFiducialCut();
+ fidCut1stYear->DoCTSFiducialCut(kFALSE) ;
+ fidCut1stYear->DoEMCALFiducialCut(kTRUE) ;
+ fidCut1stYear->DoPHOSFiducialCut(kTRUE) ;
+ fidCut1stYear->SetSimpleEMCALFiducialCut(0.7,80.,120.);
+ fidCut1stYear->SetSimplePHOSFiducialCut(0.12,260.,320.);
AliAnaPhoton *anaphoton1 = new AliAnaPhoton();
anaphoton1->SetDebug(-1); //10 for lots of messages
else anaphoton1->SwitchOffDataMC() ;
anaphoton1->SwitchOffCaloPID();
anaphoton1->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton1->SwitchOnFidutialCut();
- anaphoton1->SetFidutialCut(fidCut1stYear);
+ anaphoton1->SwitchOnFiducialCut();
+ anaphoton1->SetFiducialCut(fidCut1stYear);
anaphoton1->SetOutputAODName(Form("PhotonsForIM%s",calorimeter.Data()));
if(kPrintSettings) anaphoton1->Print("");
anapi0->SetInputAODName(Form("PhotonsForIM%s",calorimeter.Data()));
anapi0->SetCaloPID(pid);
anapi0->SetCalorimeter(calorimeter);
- anapi0->SwitchOnFidutialCut();
+ anapi0->SwitchOnFiducialCut();
anapi0->SwitchOffDataMC() ;//Access MC stack and fill more histograms
if(kPrintSettings) anapi0->Print("");
else anaphoton2->SwitchOffDataMC() ;
anaphoton2->SwitchOnCaloPID();
if(calorimeter == "EMCAL") anaphoton2->SwitchOnCaloPIDRecalculation();
- anaphoton2->SwitchOffFidutialCut();
+ anaphoton2->SwitchOffFiducialCut();
anaphoton2->SetOutputAODName(Form("DirectPhotons%s",calorimeter.Data()));
anaphoton2->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
anaphoton2->AddToHistogramsName("AnaPhotonCorr_");
// ### Correlation with Jet Finder AOD output
AliAnaParticleJetFinderCorrelation *anacorrjet = new AliAnaParticleJetFinderCorrelation();
anacorrjet->SetInputAODName(Form("DirectPhotons%s",calorimeter.Data()));
- anacorrjet->SwitchOffFidutialCut();
+ anacorrjet->SwitchOffFiducialCut();
anacorrjet->SetDebug(-1);
anacorrjet->SetConeSize(1);
anacorrjet->SelectIsolated(kTRUE); // do correlation with isolated photons
anacorrhadron->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
anacorrhadron->SetDebug(-1);
anacorrhadron->SwitchOffCaloPID();
- anacorrhadron->SwitchOffFidutialCut();
+ anacorrhadron->SwitchOffFiducialCut();
anacorrhadron->SetPtCutRange(1,100);
anacorrhadron->SetDeltaPhiCutRange(1.5,4.5);
anacorrhadron->SelectIsolated(kTRUE); // do correlation with isolated photons
// Analysis algorithm
//---------------------------------------------------------------------
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
AliAnaCalorimeterQA *anaEMCAL = new AliAnaCalorimeterQA();
anaEMCAL->SwitchOnDataMC() ;//Access MC stack and fill more histograms
//anaEMCAL->SetStyleMacro("$MACROS/style.C") ;
anaEMCAL->AddToHistogramsName("AnaCaloQA_EMCAL_");
- anaEMCAL->SetFidutialCut(fidCut);
- anaEMCAL->SwitchOnFidutialCut();
+ anaEMCAL->SetFiducialCut(fidCut);
+ anaEMCAL->SwitchOnFiducialCut();
anaEMCAL->Print("");
AliAnaCalorimeterQA *anaPHOS = new AliAnaCalorimeterQA();
anaPHOS->SwitchOnDataMC() ;//Access MC stack and fill more histograms
//anaPHOS->SetStyleMacro("$MACROS/style.C") ;
anaPHOS->AddToHistogramsName("AnaCaloQA_PHOS_");
- anaPHOS->SetFidutialCut(fidCut);
- anaPHOS->SwitchOnFidutialCut();
+ anaPHOS->SetFiducialCut(fidCut);
+ anaPHOS->SwitchOnFiducialCut();
anaPHOS->Print("");
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kTRUE) ;
- fidCut->DoEMCALFidutialCut(kFALSE) ;
- fidCut->DoPHOSFidutialCut(kFALSE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kTRUE) ;
+ fidCut->DoEMCALFiducialCut(kFALSE) ;
+ fidCut->DoPHOSFiducialCut(kFALSE) ;
- fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
+ fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
fidCut->Print("");
// reader->AddStatusArray(statusArray) ;
// reader->SwitchOnStatusSelection() ;
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// //We want tracks fitted in the detectors:
// ULong_t status=AliAODTrack::kTPCrefit;
anatrack->SetDebug(-1);//10 for lots of messages
anatrack->SetMinPt(5.);
anatrack->SetOutputAODName("Charged");
- anatrack->SwitchOffFidutialCut();
+ anatrack->SwitchOffFiducialCut();
anatrack->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
anatrack->SwitchOnDataMC() ;//Access MC stack and fill more histograms
//Set Histrograms bins and ranges
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.12,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.12,220.,320.);
- // //Fidutial cut EMCAL, 5 regions
+ // //Fiducial cut EMCAL, 5 regions
// Float_t etamax[]={0.67,0.51,0.16,-0.21,-0.61};
// TArrayF etamaxarr(5,etamax);
// fidCut->AddEMCALFidCutMaxEtaArray(etamaxarr);
//Remove the temporal AODs we create.
reader->SwitchOnCleanStdAOD();
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// //We want tracks fitted in the detectors:
// ULong_t status=AliAODTrack::kTPCrefit;
// Analysis algorithm
//---------------------------------------------------------------------
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kTRUE) ;
- fidCut2->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kTRUE) ;
+ fidCut2->DoPHOSFiducialCut(kTRUE) ;
- //fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut2->SetSimpleEMCALFidutialCut(0.5,100.,150.);
- fidCut2->SetSimplePHOSFidutialCut(0.1,240.,280.);
+ //fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut2->SetSimpleEMCALFiducialCut(0.5,100.,150.);
+ fidCut2->SetSimplePHOSFiducialCut(0.1,240.,280.);
AliCaloPID * pid = new AliCaloPID();
// use selection with simple weights
AliAnaExample *ana = new AliAnaExample();
ana->SetDebug(-1);
ana->SetCaloPID(pid);
- ana->SetFidutialCut(fidCut2);
+ ana->SetFiducialCut(fidCut2);
ana->SetDetector("PHOS");
//ana->SwitchOnDataMC();
ana->SetMinPt(0.);
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kFALSE) ;
- fidCut->DoPHOSFidutialCut(kFALSE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kFALSE) ;
+ fidCut->DoPHOSFiducialCut(kFALSE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- //fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ //fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut->Print("");
reader->SetPHOSPtMin(0.2);
reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// //We want tracks fitted in the detectors:
// ULong_t status=AliAODTrack::kTPCrefit;
// Analysis algorithm
//---------------------------------------------------------------------
//>>>> First Analysis <<<<
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kFALSE) ;
- fidCut2->DoPHOSFidutialCut(kFALSE) ;
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kFALSE) ;
+ fidCut2->DoPHOSFiducialCut(kFALSE) ;
- //fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut2->SetSimplePHOSFidutialCut(0.12,220.,320.);
+ //fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut2->SetSimplePHOSFiducialCut(0.12,220.,320.);
// //Select particles in N regions of the detectors
// Float_t etamax[]={0.67,0.51,0.16,-0.21,-0.61};
anaphoton->SetMinPt(0.2);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
- anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
+ anaphoton->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("PHOS");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOnCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOnFidutialCut();
+ anaphoton->SwitchOnFiducialCut();
anaphoton->SetOutputAODName("Photons");
anaphoton->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
//Set Histrograms bins and ranges
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kTRUE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kTRUE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut->Print("");
reader->SetPHOSPtMin(0.5);
reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// //We want tracks fitted in the detectors:
// ULong_t status=AliAODTrack::kTPCrefit;
// Analysis algorithm
//---------------------------------------------------------------------
//<<<< first analysis >>> select the photons
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kTRUE) ;
- fidCut2->DoPHOSFidutialCut(kFALSE) ;
-
- fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kTRUE) ;
+ fidCut2->DoPHOSFiducialCut(kFALSE) ;
+
+ fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut2->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut2->Print("");
AliCaloPID * pid = new AliCaloPID();
anaphoton->SetMinPt(0.2);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
- anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
+ anaphoton->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("PHOS");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOnCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOnFidutialCut();
+ anaphoton->SwitchOnFiducialCut();
anaphoton->SetOutputAODName("Photons");
anaphoton->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
//Set Histrograms bins and ranges
anacorr->SetCaloPID(pid);
anacorr->SwitchOnCaloPID();
anacorr->SwitchOnCaloPIDRecalculation(); //recommended for EMCAL
- anacorr->SwitchOffFidutialCut();
+ anacorr->SwitchOffFiducialCut();
anacorr->SetPtCutRange(1,100);
anacorr->SetDeltaPhiCutRange(1.5,4.5);
//Set Histrograms bins and ranges
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kTRUE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kTRUE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut->Print("");
reader->SetPHOSPtMin(0.5);
reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// //We want tracks fitted in the detectors:
// ULong_t status=AliAODTrack::kTPCrefit;
// Analysis algorithm
//---------------------------------------------------------------------
//<<<< first analysis >>> select the photons
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kTRUE) ;
- fidCut2->DoPHOSFidutialCut(kFALSE) ;
-
- fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kTRUE) ;
+ fidCut2->DoPHOSFiducialCut(kFALSE) ;
+
+ fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut2->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut2->Print("");
AliCaloPID * pid = new AliCaloPID();
anaphoton->SetMinPt(0.2);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
- anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
+ anaphoton->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("PHOS");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOnCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOnFidutialCut();
+ anaphoton->SwitchOnFiducialCut();
anaphoton->SetOutputAODName("Photons");
anaphoton->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
//Set Histrograms bins and ranges
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kTRUE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kTRUE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut->Print("");
reader->SetPHOSPtMin(0.5);
reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// //We want tracks fitted in the detectors:
// ULong_t status=AliAODTrack::kTPCrefit;
// Analysis algorithm
//---------------------------------------------------------------------
//<<<< first analysis >>> select the photons
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kTRUE) ;
- fidCut2->DoPHOSFidutialCut(kFALSE) ;
-
- fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kTRUE) ;
+ fidCut2->DoPHOSFiducialCut(kFALSE) ;
+
+ fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut2->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut2->Print("");
AliCaloPID * pid = new AliCaloPID();
anaphoton->SetMinPt(0.2);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
- anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
+ anaphoton->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("PHOS");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOnCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOnFidutialCut();
+ anaphoton->SwitchOnFiducialCut();
anaphoton->SetOutputAODName("Photons");
anaphoton->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
//Set Histrograms bins and ranges
anacorr->SetInputAODName("Photons");
anacorr->SwitchOnCaloPID();
anacorr->SwitchOnCaloPIDRecalculation(); //recommended for EMCAL
- anacorr->SwitchOffFidutialCut();
+ anacorr->SwitchOffFiducialCut();
anacorr->SwitchOffJetsOnlyInCTS();
anacorr->SwitchOffJetsRecalculation();
//Analysis cuts for leading particle selection
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kTRUE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kTRUE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut->Print("");
reader->SetPHOSPtMin(0.5);
reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
reader->Print("");
// Analysis algorithm
//---------------------------------------------------------------------
//<<<< first analysis >>> select the photons
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kTRUE) ;
- fidCut2->DoPHOSFidutialCut(kFALSE) ;
-
- fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kTRUE) ;
+ fidCut2->DoPHOSFiducialCut(kFALSE) ;
+
+ fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut2->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut2->Print("");
AliCaloPID * pid = new AliCaloPID();
anaphoton->SetMinPt(0.2);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
- anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
+ anaphoton->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("PHOS");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOnCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOnFidutialCut();
+ anaphoton->SwitchOnFiducialCut();
anaphoton->SetOutputAODName("Photons");
anaphoton->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
//Set Histrograms bins and ranges
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kFALSE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kFALSE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
- // //Fidutial cut EMCAL, 5 regions
+ // //Fiducial cut EMCAL, 5 regions
// Float_t etamax[]={0.67,0.51,0.16,-0.21,-0.61};
// TArrayF etamaxarr(5,etamax);
// fidCut->AddEMCALFidCutMaxEtaArray(etamaxarr);
reader->SetPHOSPtMin(0.);
reader->SetCTSPtMin(.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
//Remove the temporal AODs we create.
reader->SwitchOnCleanStdAOD();
// Analysis algorithm
//---------------------------------------------------------------------
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kFALSE) ;
- fidCut2->DoPHOSFidutialCut(kTRUE) ;
- fidCut2->SetSimplePHOSFidutialCut(0.1,240.,280.);
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kFALSE) ;
+ fidCut2->DoPHOSFiducialCut(kTRUE) ;
+ fidCut2->SetSimplePHOSFiducialCut(0.1,240.,280.);
AliCaloPID * pid = new AliCaloPID();
// use selection with simple weights
AliAnaExample *ana = new AliAnaExample();
ana->SetDebug(-1);
ana->SetCaloPID(pid);
- ana->SetFidutialCut(fidCut2);
+ ana->SetFiducialCut(fidCut2);
ana->SetDetector("PHOS");
//ana->SwitchOnDataMC();
ana->SetMinPt(0.);
printf("ConfigAnalysis() \n");
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.12,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.12,220.,320.);
fidCut->Print("");
reader->SetPHOSPtMin(0.2);
//reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
//Remove the temporal AODs we create.
reader->SwitchOnCleanStdAOD();
anaphoton->SetMinPt(0.2);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
-// anaphoton->SetFidutialCut(fidCut); //More acceptance selections if needed at this level
+// anaphoton->SetFiducialCut(fidCut); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("PHOS");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOffCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOffFidutialCut();
+ anaphoton->SwitchOffFiducialCut();
anaphoton->SetOutputAODName("PhotonsPHOS");
anaphoton->SetOutputAODClassName("AliAODPWG4Particle");
anaphoton->AddToHistogramsName("AnaPhotonPHOS_");
anaomega->SetNPtBinsMinMax(200,0,20.);
anaomega->SetNMassBinsMinMas(200,0,1.);
anaomega->SetPi0MassPeakWidthCut(0.015);
- anaomega->SwitchOnFidutialCut();
+ anaomega->SwitchOnFiducialCut();
anaomega->SwitchOnDataMC() ;//Access MC stack and fill more histograms
anaomega->AddToHistogramsName("AnaNeuPHOS_");
anaomega->Print("");
anaphoton1->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton1->SwitchOffCaloPID();
anaphoton1->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton1->SwitchOffFidutialCut();
+ anaphoton1->SwitchOffFiducialCut();
anaphoton1->SetOutputAODName("PhotonsEMCAL");
anaphoton1->SetOutputAODClassName("AliAODPWG4Particle");
anaphoton1->AddToHistogramsName("AnaPhotonEMCAL_");
anaomega1->SetNPtBinsMinMax(200,0,20.);
anaomega1->SetNMassBinsMinMas(200,0,1.);
anaomega1->SetPi0MassPeakWidthCut(0.015);
- anaomega1->SwitchOnFidutialCut();
+ anaomega1->SwitchOnFiducialCut();
anaomega1->SwitchOnDataMC() ;//Access MC stack and fill more histograms
anaomega1->AddToHistogramsName("AnaNeuEMCAL_");
anaomega1->Print("");
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kFALSE) ;
- fidCut->DoPHOSFidutialCut(kFALSE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kFALSE) ;
+ fidCut->DoPHOSFiducialCut(kFALSE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- //fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ //fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut->Print("");
//reject events with large difference between ptHard and triggered jet
//reader->SetPtHardAndJetPtComparison(kTRUE);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
//Embedding/mixing/merging analysis with other events in another AOD file
reader->SetSecondInputFileName("./aod.root");
// Analysis algorithm
//---------------------------------------------------------------------
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kFALSE) ;
- fidCut2->DoPHOSFidutialCut(kFALSE) ;
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kFALSE) ;
+ fidCut2->DoPHOSFiducialCut(kFALSE) ;
- //fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- //fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ //fidCut2->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut2->Print("");
ana->SetMinPt(5.);
ana->SetMinDistanceToBadChannel(2, 4, 5);
ana->SetCaloPID(pid);
- ana->SetFidutialCut(fidCut2);
+ ana->SetFiducialCut(fidCut2);
ana->SetCalorimeter("EMCAL");
ana->SwitchOnDataMC() ;//Access MC stack and fill more histograms
ana->SwitchOffCaloPID();
ana->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
ana->SwitchOnTrackMatchRejection(); //Only in use when OnCaloPID
- ana->SwitchOffFidutialCut();
+ ana->SwitchOffFiducialCut();
ana->SetOutputAODName("Photons");
ana->SetOutputAODClassName("AliAODPWG4Particle");
//Set Histrograms bins and ranges
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kFALSE) ;
- fidCut->DoPHOSFidutialCut(kFALSE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kFALSE) ;
+ fidCut->DoPHOSFiducialCut(kFALSE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- //fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ //fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut->Print("");
//reject events with large difference between ptHard and triggered jet
//reader->SetPtHardAndJetPtComparison(kTRUE);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// Analysis with tracks, select only tracks with
// following bits
// Analysis algorithm
//---------------------------------------------------------------------
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kFALSE) ;
- fidCut2->DoPHOSFidutialCut(kFALSE) ;
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kFALSE) ;
+ fidCut2->DoPHOSFiducialCut(kFALSE) ;
- //fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- //fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ //fidCut2->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut2->Print("");
ana->SetMinPt(5.);
ana->SetMinDistanceToBadChannel(2, 4, 5);
ana->SetCaloPID(pid);
- ana->SetFidutialCut(fidCut2);
+ ana->SetFiducialCut(fidCut2);
ana->SetCalorimeter("PHOS");
ana->SwitchOnDataMC() ;//Access MC stack and fill more histograms
ana->SwitchOffCaloPID();
ana->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
ana->SwitchOnTrackMatchRejection(); //Only in use when OnCaloPID
- ana->SwitchOffFidutialCut();
+ ana->SwitchOffFiducialCut();
ana->SetOutputAODName("Photons");
ana->SetOutputAODClassName("AliAODPWG4Particle");
//Set Histrograms bins and ranges
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kTRUE) ;
- fidCut->DoEMCALFidutialCut(kTRUE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kTRUE) ;
+ fidCut->DoEMCALFiducialCut(kTRUE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- //fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ //fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut->Print("");
//reject events with large difference between ptHard and triggered jet
//reader->SetPtHardAndJetPtComparison(kTRUE);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
//Anaysis of final particles, not pi0/eta etc.
TArrayI statusArray(1) ;
// Analysis algorithm
//---------------------------------------------------------------------
- //Detector Fidutial Cuts for analysis part
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoCTSFidutialCut(kFALSE) ;
- fidCut2->DoEMCALFidutialCut(kFALSE) ;
- fidCut2->DoPHOSFidutialCut(kFALSE) ;
+ //Detector Fiducial Cuts for analysis part
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoCTSFiducialCut(kFALSE) ;
+ fidCut2->DoEMCALFiducialCut(kFALSE) ;
+ fidCut2->DoPHOSFiducialCut(kFALSE) ;
- //fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- //fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
+ //fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ //fidCut2->SetSimplePHOSFiducialCut(0.13,220.,320.);
fidCut2->Print("");
ana->SetDebug(-1);
ana->SetMinPt(5.);
ana->SetMinDistanceToBadChannel(2, 4, 5);
- ana->SetFidutialCut(fidCut2);
+ ana->SetFiducialCut(fidCut2);
ana->SetCalorimeter("PHOS");
ana->SwitchOnDataMC() ;//Access MC stack and fill more histograms
ana->SwitchOffCaloPID(); //No need with MC reader
ana->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL, no need with MC reader
ana->SwitchOffTrackMatchRejection(); //Only in use when OnCaloPID
- ana->SwitchOffFidutialCut();
+ ana->SwitchOffFiducialCut();
ana->SetOutputAODName("Photons");
ana->SetOutputAODClassName("AliAODPWG4Particle");
//Set Histrograms bins and ranges
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kFALSE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kFALSE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.12,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.12,220.,320.);
fidCut->Print("");
//reader->SetPHOSPtMin(0.5);
//reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
//Remove the temporal AODs we create.
reader->SwitchOnCleanStdAOD();
anaphoton->SetMinPt(1);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
- //anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
+ //anaphoton->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("EMCAL");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOffCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOffFidutialCut();
+ anaphoton->SwitchOffFiducialCut();
anaphoton->SetOutputAODName("Photons");
anaphoton->SetOutputAODClassName("AliAODPWG4Particle");
//Set Histrograms bins and ranges
// anaphoton->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
// anaphoton->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ; anaphoton->Print("");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoPHOSFidutialCut(kTRUE) ;
- fidCut2->SetSimplePHOSFidutialCut(0.12,220.,320.);
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoPHOSFiducialCut(kTRUE) ;
+ fidCut2->SetSimplePHOSFiducialCut(0.12,220.,320.);
fidCut->Print("");
// anapi0->SetNRPBin(6); //number of bins in reaction plain
anapi0->SetNMaxEvMix(2);//Maximal number of events for mixing
anapi0->SetCalorimeter("EMCAL");
- anapi0->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level. Not used if real geometry
- anapi0->SwitchOnFidutialCut();
+ anapi0->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level. Not used if real geometry
+ anapi0->SwitchOnFiducialCut();
anapi0->SwitchOnDataMC() ;//Access MC stack and fill more histograms
anapi0->Print("");
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kFALSE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kFALSE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- fidCut->SetSimplePHOSFidutialCut(0.12,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ fidCut->SetSimplePHOSFiducialCut(0.12,220.,320.);
fidCut->Print("");
reader->SetPHOSPtMin(0.5);
//reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// //We want tracks fitted in the detectors:
// ULong_t status=AliAODTrack::kTPCrefit;
anaphoton->SetMinPt(1);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
- //anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
+ //anaphoton->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("EMCAL");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOffCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOffFidutialCut();
+ anaphoton->SwitchOffFiducialCut();
anaphoton->SetOutputAODName("Photons");
anaphoton->SetOutputAODClassName("AliAODPWG4Particle");
//Set Histrograms bins and ranges
// anaphoton->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
// anaphoton->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ; anaphoton->Print("");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut2 = new AliFidutialCut();
- fidCut2->DoPHOSFidutialCut(kTRUE) ;
- fidCut2->SetSimplePHOSFidutialCut(0.12,220.,320.);
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut2 = new AliFiducialCut();
+ fidCut2->DoPHOSFiducialCut(kTRUE) ;
+ fidCut2->SetSimplePHOSFiducialCut(0.12,220.,320.);
fidCut->Print("");
printf("======================== \n");
- //Detector Fidutial Cuts
- AliFidutialCut * fidCut = new AliFidutialCut();
- fidCut->DoCTSFidutialCut(kFALSE) ;
- fidCut->DoEMCALFidutialCut(kFALSE) ;
- fidCut->DoPHOSFidutialCut(kTRUE) ;
+ //Detector Fiducial Cuts
+ AliFiducialCut * fidCut = new AliFiducialCut();
+ fidCut->DoCTSFiducialCut(kFALSE) ;
+ fidCut->DoEMCALFiducialCut(kFALSE) ;
+ fidCut->DoPHOSFiducialCut(kTRUE) ;
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
- //fidCut->SetSimplePHOSFidutialCut(0.12,220.,320.);
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);
+ //fidCut->SetSimplePHOSFiducialCut(0.12,220.,320.);
fidCut->Print("");
//reader->SetPHOSPtMin(0.5);
reader->SetCTSPtMin(0.2);
- reader->SetFidutialCut(fidCut);
+ reader->SetFiducialCut(fidCut);
// //We want tracks fitted in the detectors:
// ULong_t status=AliAODTrack::kTPCrefit;
anaphoton->SetMinPt(1);
anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
anaphoton->SetCaloPID(pid);
- //anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
+ //anaphoton->SetFiducialCut(fidCut2); //More acceptance selections if needed at this level
anaphoton->SetCalorimeter("EMCAL");
anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
anaphoton->SwitchOffCaloPID();
anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
- anaphoton->SwitchOffFidutialCut();
+ anaphoton->SwitchOffFiducialCut();
anaphoton->SetOutputAODName("Photons");
anaphoton->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
//Set Histrograms bins and ranges
}\r
}\r
\r
- //Detector Fidutial Cuts\r
- AliFidutialCut * fidCut = new AliFidutialCut();\r
- fidCut->DoCTSFidutialCut(kFALSE) ;\r
- fidCut->DoEMCALFidutialCut(kFALSE) ;\r
- fidCut->DoPHOSFidutialCut(kFALSE) ;\r
+ //Detector Fiducial Cuts\r
+ AliFiducialCut * fidCut = new AliFiducialCut();\r
+ fidCut->DoCTSFiducialCut(kFALSE) ;\r
+ fidCut->DoEMCALFiducialCut(kFALSE) ;\r
+ fidCut->DoPHOSFiducialCut(kFALSE) ;\r
\r
- //fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);\r
- //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);\r
- //fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);\r
+ //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);\r
+ //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);\r
+ //fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);\r
\r
fidCut->Print("");\r
\r
//reject events with large difference between ptHard and triggered jet\r
//reader->SetPtHardAndJetPtComparison(kTRUE);\r
\r
- reader->SetFidutialCut(fidCut);\r
+ reader->SetFiducialCut(fidCut);\r
\r
if(!kInputIsESD){\r
// Analysis with tracks, select only tracks with\r
reader->Print("");\r
\r
\r
- //Detector Fidutial Cuts\r
- AliFidutialCut * fidCut2 = new AliFidutialCut();\r
- fidCut2->DoEMCALFidutialCut(kTRUE) ;\r
- fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);\r
+ //Detector Fiducial Cuts\r
+ AliFiducialCut * fidCut2 = new AliFiducialCut();\r
+ fidCut2->DoEMCALFiducialCut(kTRUE) ;\r
+ fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);\r
\r
- fidCut2->DoCTSFidutialCut(kTRUE) ;\r
- fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.); \r
+ fidCut2->DoCTSFiducialCut(kTRUE) ;\r
+ fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.); \r
\r
fidCut2->Print("");\r
\r