add flag for the number of leading jets excluded for random cones
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 17 Mar 2011 15:03:06 +0000 (15:03 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 17 Mar 2011 15:03:06 +0000 (15:03 +0000)
JETAN/AliAnalysisTaskJetCluster.cxx
JETAN/AliAnalysisTaskJetCluster.h

index 2265781..1c36bfa 100644 (file)
@@ -96,6 +96,7 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(): AliAnalysisTaskSE(),
   fTrackTypeRec(kTrackUndef),
   fTrackTypeGen(kTrackUndef),  
   fNSkipLeadingRan(0),
+  fNSkipLeadingCone(0),
   fNRandomCones(0),
   fAvgTrials(1),
   fExternalWeight(1),
@@ -204,6 +205,7 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(const char* name):
   fTrackTypeRec(kTrackUndef),
   fTrackTypeGen(kTrackUndef),
   fNSkipLeadingRan(0),
+  fNSkipLeadingCone(0),
   fNRandomCones(0),
   fAvgTrials(1),
   fExternalWeight(1),    
@@ -353,8 +355,8 @@ void AliAnalysisTaskJetCluster::UserCreateOutputObjects()
        }
       }
       // create the branch for the random cones with the same R 
-      TString cName = Form("%sRandomCone",fNonStdBranch.Data());
-
+      TString cName = Form("%sRandomConeSkip%02d",fNonStdBranch.Data(),fNSkipLeadingCone);
+  
       if(fNRandomCones>0){
        if(!AODEvent()->FindListObject(cName.Data())){
          fTCARandomConesOut = new TClonesArray("AliAODJet", 0);
@@ -983,7 +985,7 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/)
        Double_t pC  = TMath::Sqrt(pTC*pTC+pZC*pZC); 
        AliAODJet tmpRecC (pXC,pYC,pZC, pC); 
        bool skip = false;
-       for(int jj = 0; jj < TMath::Min(nRec,2);jj++){// test for overlap with leading jets
+       for(int jj = 0; jj < TMath::Min(nRec,fNSkipLeadingCone);jj++){// test for overlap with leading jets
         AliAODJet jet (sortedJets[jj].px(), sortedJets[jj].py(), sortedJets[jj].pz(), sortedJets[jj].E());
         if(jet.DeltaR(& tmpRecC)<2.*fRparam+0.2){
           skip = true;
index 0e2bb20..c9b49a2 100644 (file)
@@ -65,6 +65,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     virtual void SetBackgroundBranch(const char* c){fBackgroundBranch = c;}
     virtual const char* GetBackgroundBranch(){return fBackgroundBranch.Data();}    
     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 SetJetOutputBranch(const char *c){fNonStdBranch = c;}
@@ -122,7 +123,8 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     UInt_t        fFilterMask;            // filter bit for slecected tracks
     Int_t         fTrackTypeRec;          // type of tracks used for FF 
     Int_t         fTrackTypeGen;          // type of tracks used for FF 
-    Int_t         fNSkipLeadingRan;        // number of leading tracks to be skipped in the randomized event
+    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
     Float_t       fAvgTrials;             // Average nimber of trials
     Float_t       fExternalWeight;        // external weight
@@ -228,7 +230,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     TList *fHistList; //!leading tracks to be skipped in the randomized event Output list
    
 
-    ClassDef(AliAnalysisTaskJetCluster, 14) 
+    ClassDef(AliAnalysisTaskJetCluster, 16) 
 };
  
 #endif