Take from Andrews git
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 4 Oct 2010 20:34:12 +0000 (20:34 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 4 Oct 2010 20:34:12 +0000 (20:34 +0000)
PWG4/TwoPartCorr/AutoCorr.C
PWG4/TwoPartCorr/AutoCorr.h

index 0d10e1a..bbb0237 100644 (file)
@@ -6,7 +6,8 @@ ClassImp(AutoCorr)
 
 Int_t AutoCorr::InitEventPools(Int_t depth, 
                               Int_t nMultBins, Double_t multbin[], 
-                              Int_t nZvtxBins, Double_t zvtxbin[])
+                              Int_t nZvtxBins, Double_t zvtxbin[],
+                              Double_t ptMin, Double_t ptMax)
 {
   // First assign AutoCorr members
   fNMultBins = nMultBins;
@@ -17,7 +18,8 @@ Int_t AutoCorr::InitEventPools(Int_t depth,
     for (Int_t iZ=0; iZ<nZvtxBins; iZ++) {
       evp.push_back(new EventPool(depth, 
                                  multbin[iM], multbin[iM+1], 
-                                 zvtxbin[iZ], zvtxbin[iZ+1] ));
+                                 zvtxbin[iZ], zvtxbin[iZ+1],
+                                 ptMin, ptMax ));
     }
     fEvPool.push_back(evp);
   }
@@ -136,7 +138,7 @@ Bool_t AutoCorr::IsPairOk(const MyPart &t1, const MyPart &t2) const
   Double_t deta = DeltaEta(t1, t2);
   Double_t dphi = DeltaPhi(t1, t2);
   Double_t dpmax = 0.03;
-  Double_t demax = 0.01;
+  Double_t demax = 0.03;
   Double_t dr = dphi*dphi/(dpmax*dpmax) + deta*deta/(demax*demax);
   return (dr > 1);
 }
@@ -145,8 +147,8 @@ Bool_t AutoCorr::IsMixedPairOk(const MyPart &t1, const MyPart &t2) const
 {
   Double_t deta = DeltaEta(t1, t2);
   Double_t dphi = DeltaPhi(t1, t2);
-  Double_t dpmax = 0.04;
-  Double_t demax = 0.04;
+  Double_t dpmax = 0.03;
+  Double_t demax = 0.03;
   Double_t dr = dphi*dphi/(dpmax*dpmax) + deta*deta/(demax*demax);
   return (dr > 1);
 }
index b09f696..1b31096 100644 (file)
@@ -38,8 +38,9 @@ public:
   EventPool* GetEventPool(Int_t iMult, Int_t iZvtx)             const;
   Int_t      InitEventPools(Int_t depth, 
                             Int_t nmultbins, Double_t multbins[], 
-                            Int_t nzvtxbins, Double_t zvtxbins[]);
-  Int_t      UpdatePools(Int_t iEvent, const MyHeader *ev, TClonesArray *trk);
+                            Int_t nzvtxbins, Double_t zvtxbins[],
+                           Double_t ptMin, Double_t ptMax);
+  Int_t      UpdatePools(Int_t iEvent, const MyHeader* ev, TClonesArray* trk);
 
 protected:
   Int_t      fNMultBins;                         // mult bins