Adding option for fixed cone (M. Zimmermann)
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 18 Jun 2012 14:48:21 +0000 (14:48 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 18 Jun 2012 14:48:21 +0000 (14:48 +0000)
PWGJE/UserTasks/AliAnalysisTaskJetHBOM.cxx
PWGJE/UserTasks/AliAnalysisTaskJetHBOM.h

index 9e04069..2ea4c70 100644 (file)
@@ -90,6 +90,9 @@ AliAnalysisTaskJetHBOM::AliAnalysisTaskJetHBOM():
   fNSkipLeadingRan(0),
   fNSkipLeadingCone(0),
   fNRandomCones(0),
+  randCone_pos(0),
+  randCone_Eta(0),
+  randCone_Phi(0),
   fNHBOM(0),
   fTrackEtaWindow(0.9),    
   //fRecEtaWindow(0.5),
@@ -164,6 +167,9 @@ AliAnalysisTaskJetHBOM::AliAnalysisTaskJetHBOM(const char* name):
   fNSkipLeadingRan(0),
   fNSkipLeadingCone(0),
   fNRandomCones(0),
+  randCone_pos(0),
+  randCone_Eta(0),
+  randCone_Phi(0),
   fNHBOM(0),
   fTrackEtaWindow(0.9),    
   //fRecEtaWindow(0.5),
@@ -614,7 +620,7 @@ void AliAnalysisTaskJetHBOM::UserExec(Option_t */*option*/)
 
 
     // generate random cones
-    if(fTCARandomConesOut){       
+    if(fTCARandomConesOut){
       // create a random jet within the acceptance
       Double_t etaMax = fTrackEtaWindow - fRparam;//0.9 - 0.4
       Int_t nCone = 0;
@@ -622,6 +628,11 @@ void AliAnalysisTaskJetHBOM::UserExec(Option_t */*option*/)
       Double_t pTC = 1; // small number
       Double_t etaC = etaMax*2.*(fRandom->Rndm()-0.5); // +- etamax
       Double_t phiC = fRandom->Rndm()*2.*TMath::Pi(); // 0 - 2pi
+      // use fixed position for random Cones
+      if(randCone_pos){
+       etaC = randCone_Eta;
+       phiC = randCone_Phi;
+      }
       // massless jet
       Double_t thetaC = 2.*TMath::ATan(TMath::Exp(-etaC));  
       Double_t pZC = pTC/TMath::Tan(thetaC);
index f2eddb5..f1bd18d 100644 (file)
@@ -70,6 +70,7 @@ class AliAnalysisTaskJetHBOM : public AliAnalysisTaskSE
     virtual void SetNSkipLeadingRan(Int_t x){fNSkipLeadingRan = x;}
     virtual void SetNSkipLeadingCone(Int_t x){fNSkipLeadingCone = x;}
     virtual void SetNRandomCones(Int_t x){fNRandomCones = x;}
+    virtual void SetRandConePos(Double_t eta, Double_t phi){randCone_pos=1;randCone_Eta=eta;randCone_Phi=phi;}
 
     virtual void SetfNHBOM(Int_t x){fNHBOM = x;};
     virtual void SetEfficiencyPt(TH1F *h){fh1efficiencyPt = h;}
@@ -151,6 +152,9 @@ class AliAnalysisTaskJetHBOM : public AliAnalysisTaskSE
     Int_t         fNSkipLeadingRan;       // number of leading tracks to be skipped in the randomized event
     Int_t         fNSkipLeadingCone;      // number of leading jets to be for the random cones
     Int_t         fNRandomCones;          // number of generated random cones
+    Bool_t        randCone_pos;           // use fixed position for random cones
+    Double_t      randCone_Eta;           // eta for random Cone at fixed position
+    Double_t      randCone_Phi;           // phi for random Cone at fixed position
     Int_t         fNHBOM;                 // number of detector runs
     Float_t       fTrackEtaWindow;        // eta window used for corraltion plots between rec and gen 
     //    Float_t       fRecEtaWindow;          // eta window used for corraltion plots between rec and gen