Adding analysis task

[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisEtSelectorPhos.cxx

diff --git a/PWGLF/totEt/AliAnalysisEtSelectorPhos.cxx b/PWGLF/totEt/AliAnalysisEtSelectorPhos.cxx
index 3c6fe9998e9ba3620cfc811b82c75177eb11b936..71b33f986a6e3437466c8c523b114515e9d3380b 100644 (file)
--- a/PWGLF/totEt/AliAnalysisEtSelectorPhos.cxx
+++ b/PWGLF/totEt/AliAnalysisEtSelectorPhos.cxx
@@ -1,3 +1,11 @@
+//_________________________________________________________________________
+//  Utility Class for transverse energy studies
+//  Selector Base class for PHOS
+//  -  
+// implementation file
+//
+//*-- Authors: Oystein Djuvsland (Bergen)
+//_________________________________________________________________________
 #include "AliAnalysisEtSelectorPhos.h"
 #include "AliAnalysisEtCuts.h"
 #include "AliESDCaloCluster.h"
@@ -18,7 +26,6 @@ AliAnalysisEtSelectorPhos::AliAnalysisEtSelectorPhos(AliAnalysisEtCuts* cuts): A
 ,fBadMapM2(0)
 ,fBadMapM3(0)
 ,fBadMapM4(0)
-,fInitialized(kFALSE)
 ,fMatrixInitialized(kFALSE)
 {
   
@@ -30,7 +37,7 @@ AliAnalysisEtSelectorPhos::~AliAnalysisEtSelectorPhos()
 }
 
 TRefArray* AliAnalysisEtSelectorPhos::GetClusters()
-{
+{ // Get clusters
   if(!fClusterArray) fClusterArray = new TRefArray;
   
   if(fClusterArray)
@@ -46,7 +53,7 @@ TRefArray* AliAnalysisEtSelectorPhos::GetClusters()
 }
 
 Int_t AliAnalysisEtSelectorPhos::Init(const AliESDEvent* event)
-{
+{ // Init
   
   AliAnalysisEtSelector::Init(event);
   Printf("Initializing selector for run: %d", event->GetRunNumber());
@@ -72,19 +79,22 @@ Int_t AliAnalysisEtSelectorPhos::Init(const AliESDEvent* event)
   return 0;
 }
 
-Bool_t AliAnalysisEtSelectorPhos::CutMinEnergy(const AliESDCaloCluster& cluster) const
+Bool_t AliAnalysisEtSelectorPhos::PassMinEnergyCut(const AliESDCaloCluster& cluster) const
 {
+  
+//    std::cout << fCuts->GetReconstructedPhosClusterEnergyCut();
   return cluster.E() > fCuts->GetReconstructedPhosClusterEnergyCut();
 }
 
-Bool_t AliAnalysisEtSelectorPhos::CutMinEnergy(const TParticle& part) const
+Bool_t AliAnalysisEtSelectorPhos::PassMinEnergyCut(const TParticle& part) const
 {
+//    std::cout << fCuts->GetReconstructedPhosClusterEnergyCut();
     return part.Energy() > fCuts->GetReconstructedPhosClusterEnergyCut();
 }
 
 
-Bool_t AliAnalysisEtSelectorPhos::CutDistanceToBadChannel(const AliESDCaloCluster& cluster) const
-{
+Bool_t AliAnalysisEtSelectorPhos::PassDistanceToBadChannelCut(const AliESDCaloCluster& cluster) const
+{ // cut distance to bad channel
   if(!fMatrixInitialized)
   {
     Printf("Misalignment matrices are not initialized");
@@ -96,6 +106,7 @@ Bool_t AliAnalysisEtSelectorPhos::CutDistanceToBadChannel(const AliESDCaloCluste
     TVector3 glVec(gPos);
     fGeoUtils->GlobalPos2RelId(glVec, relId);
 
+    //std::cout << "In phos distance to bad channel cut!" << std::endl;
     TVector3 locVec;
     fGeoUtils->Global2Local(locVec, glVec, relId[0]);
 //    std::cout << fGeoUtils << std::endl;
@@ -187,8 +198,8 @@ Bool_t AliAnalysisEtSelectorPhos::CutDistanceToBadChannel(const AliESDCaloCluste
 
 }
 
-Bool_t AliAnalysisEtSelectorPhos::CutTrackMatching(const AliESDCaloCluster& cluster) const
-{
+Bool_t AliAnalysisEtSelectorPhos::PassTrackMatchingCut(const AliESDCaloCluster& cluster) const
+{ // cut track matching
 
   if(!fMatrixInitialized)
   {
@@ -199,7 +210,7 @@ Bool_t AliAnalysisEtSelectorPhos::CutTrackMatching(const AliESDCaloCluster& clus
   // cluster->GetTrackDx(), cluster->GetTrackDz(), event->GetTrack(trackMatchedIndex)->Pt(), event->GetTrack(trackMatchedIndex)->Charge(), ev
   
   Int_t nTracksMatched = cluster.GetNTracksMatched();
-  if(nTracksMatched == 0) return kFALSE;
+  if(nTracksMatched == 0) return kTRUE;
   
   Int_t trackMatchedIndex = cluster.GetTrackMatchedIndex();
   if(trackMatchedIndex < 0) return kTRUE;
@@ -243,7 +254,7 @@ Bool_t AliAnalysisEtSelectorPhos::CutTrackMatching(const AliESDCaloCluster& clus
 }
 
 int AliAnalysisEtSelectorPhos::LoadGeometry()
-{
+{ // load geometry
 
   fGeoUtils = AliPHOSGeometry::GetInstance("IHEP");
     // ifstream f("badchannels.txt", ios::in);
@@ -251,7 +262,7 @@ int AliAnalysisEtSelectorPhos::LoadGeometry()
 }
 
 int AliAnalysisEtSelectorPhos::LoadBadMaps()
-{
+{ // load bad maps
 TFile *f = TFile::Open("badchannels.root", "READ");
 
     if(!f)
@@ -282,12 +293,6 @@ TFile *f = TFile::Open("badchannels.root", "READ");
     
 }
 
-void AliAnalysisEtSelectorPhos::SetEvent(const AliESDEvent* event)
-{
-    //AliAnalysisEtSelector::SetEvent(event);
-    fEvent = event;
-    if(!fInitialized) Init(event);
-}
 
 Bool_t AliAnalysisEtSelectorPhos::CutGeometricalAcceptance(const TParticle& part) const
 {
@@ -302,3 +307,4 @@ Bool_t AliAnalysisEtSelectorPhos::CutGeometricalAcceptance(const AliVTrack& trac
            track.Phi() > fCuts->GetGeometryPhosPhiAccMaxCut()*TMath::Pi()/180. &&
            track.Phi() < fCuts->GetGeometryPhosPhiAccMinCut()*TMath::Pi()/180.;
 }
+