add option to use / neglect clusters which are not in chamber

update macro run.C to load also libTENDER.so

diff --git a/PWG1/TRD/AliTRDcheckDET.cxx b/PWG1/TRD/AliTRDcheckDET.cxx
index 3b34d5dd9607924b430a384d71f2ba87b5356225..d3a006622354a55462877c41fe5907c8d5773915 100644 (file)
--- a/PWG1/TRD/AliTRDcheckDET.cxx
+++ b/PWG1/TRD/AliTRDcheckDET.cxx
@@ -76,6 +76,7 @@ AliTRDcheckDET::AliTRDcheckDET():
   ,fTriggerNames(0x0)
   ,fReconstructor(0x0)
   ,fGeo(0x0)
   ,fTriggerNames(0x0)
   ,fReconstructor(0x0)
   ,fGeo(0x0)

+  ,fFlags(0)

 {
   //
   // Default constructor
 {
   //
   // Default constructor


@@ -873,7 +874,7 @@ TH1 *AliTRDcheckDET::PlotPHt(const AliTRDtrackV1 *track){
     Int_t detector = tracklet->GetDetector();
     tracklet->ResetClusterIter();
     while((c = tracklet->NextCluster())){
     Int_t detector = tracklet->GetDetector();
     tracklet->ResetClusterIter();
     while((c = tracklet->NextCluster())){

-      if(!c->IsInChamber()) continue;
+      if(!IsUsingClustersOutsideChamber() && !c->IsInChamber()) continue;

       Int_t localtime        = c->GetLocalTimeBin();
       Double_t absoluteCharge = TMath::Abs(c->GetQ());
       h->Fill(localtime, absoluteCharge);
       Int_t localtime        = c->GetLocalTimeBin();
       Double_t absoluteCharge = TMath::Abs(c->GetQ());
       h->Fill(localtime, absoluteCharge);


@@ -936,7 +937,7 @@ TH1 *AliTRDcheckDET::PlotPHx(const AliTRDtrackV1 *track){
     if(!(tracklet = fkTrack->GetTracklet(itl)) || !(tracklet->IsOK())) continue;
     tracklet->ResetClusterIter();
     while((c = tracklet->NextCluster())){
     if(!(tracklet = fkTrack->GetTracklet(itl)) || !(tracklet->IsOK())) continue;
     tracklet->ResetClusterIter();
     while((c = tracklet->NextCluster())){

-      if(!c->IsInChamber()) continue;
+      if(!IsUsingClustersOutsideChamber() && !c->IsInChamber()) continue;

       x = c->GetX()-AliTRDcluster::GetXcorr(c->GetLocalTimeBin());
       y = c->GetY()-AliTRDcluster::GetYcorr(AliTRDgeometry::GetLayer(c->GetDetector()), c->GetCenter());
 
       x = c->GetX()-AliTRDcluster::GetXcorr(c->GetLocalTimeBin());
       y = c->GetY()-AliTRDcluster::GetYcorr(AliTRDgeometry::GetLayer(c->GetDetector()), c->GetCenter());
 


@@ -1207,4 +1208,4 @@ Bool_t AliTRDcheckDET::MakeBarPlot(TH1 *histo, Int_t color){
   histo->SetBarWidth(.6);
   histo->Draw("bar1");
   return kTRUE;
   histo->SetBarWidth(.6);
   histo->Draw("bar1");
   return kTRUE;

-}
+}
\ No newline at end of file

diff --git a/PWG1/TRD/AliTRDcheckDET.h b/PWG1/TRD/AliTRDcheckDET.h
index 04fef04dd431338b81c5d637ad1d532f54d01175..b5b55e1ece99d497aaea2347a9df867af86c30f2 100644 (file)
--- a/PWG1/TRD/AliTRDcheckDET.h
+++ b/PWG1/TRD/AliTRDcheckDET.h
@@ -75,9 +75,14 @@ public:
   virtual Bool_t PostProcess();
   virtual Bool_t GetRefFigure(Int_t ifig);
   
   virtual Bool_t PostProcess();
   virtual Bool_t GetRefFigure(Int_t ifig);
   

+  Bool_t IsUsingClustersOutsideChamber() const {return TESTBIT(fFlags, kUseClustersOutsideChamber);}
+  void UseClustersOutsideChamber(Bool_t b = kTRUE) {if(b) SETBIT(fFlags, kUseClustersOutsideChamber); else CLRBIT(fFlags, kUseClustersOutsideChamber);}

   void SetRecoParam(AliTRDrecoParam *r);
 
 private:
   void SetRecoParam(AliTRDrecoParam *r);
 
 private:

+  enum{
+    kUseClustersOutsideChamber
+  };

   AliTRDcheckDET(const AliTRDcheckDET &);
   AliTRDcheckDET& operator=(const AliTRDcheckDET &);
   void GetDistanceToTracklet(Double_t *dist, AliTRDseedV1 * const tracklet, AliTRDcluster * const c);
   AliTRDcheckDET(const AliTRDcheckDET &);
   AliTRDcheckDET& operator=(const AliTRDcheckDET &);
   void GetDistanceToTracklet(Double_t *dist, AliTRDseedV1 * const tracklet, AliTRDcluster * const c);


@@ -90,8 +95,8 @@ private:
   TMap *fTriggerNames;                 //! Containing trigger class names
   AliTRDReconstructor *fReconstructor; // TRD Reconstructor
   AliTRDgeometry *fGeo;                // TRD Geometry object
   TMap *fTriggerNames;                 //! Containing trigger class names
   AliTRDReconstructor *fReconstructor; // TRD Reconstructor
   AliTRDgeometry *fGeo;                // TRD Geometry object

+  UChar_t fFlags;                      // Flags for setting

     
   ClassDef(AliTRDcheckDET, 1)
 };
 #endif
     
   ClassDef(AliTRDcheckDET, 1)
 };
 #endif

-

diff --git a/PWG1/TRD/run.C b/PWG1/TRD/run.C
index c29d9834cb3cf8d6615a04c7d924fe62b002928d..8838dff9cd237389a60801bfa04b17cdeac3ab95 100644 (file)
--- a/PWG1/TRD/run.C
+++ b/PWG1/TRD/run.C
@@ -99,6 +99,7 @@ void run(Char_t *trd="ALL", const Char_t *files=0x0, Long64_t nev=1234567890, Lo
   AliLog::SetGlobalLogLevel(AliLog::kError);
   if(gSystem->Load("libANALYSIS.so")<0) return;
   if(gSystem->Load("libANALYSISalice.so")<0) return;
   AliLog::SetGlobalLogLevel(AliLog::kError);
   if(gSystem->Load("libANALYSIS.so")<0) return;
   if(gSystem->Load("libANALYSISalice.so")<0) return;

+  if(gSystem->Load("libTENDER.so")<0) return;

   if(gSystem->Load("libPWG1.so")<0) return;
 
   Bool_t fHasMCdata =  HasReadMCData(trd);
   if(gSystem->Load("libPWG1.so")<0) return;
 
   Bool_t fHasMCdata =  HasReadMCData(trd);