exclusive bins for high pT dihadron study

[u/mrichter/AliRoot.git] / ITS / MakeITSPlaneEfficiencySPD.C

diff --git a/ITS/MakeITSPlaneEfficiencySPD.C b/ITS/MakeITSPlaneEfficiencySPD.C
index 9508d097bf8bed8db58318ba286951b700934caf..b78054dbe680a35fc1010835adf3c5f315382bff 100644 (file)
--- a/ITS/MakeITSPlaneEfficiencySPD.C
+++ b/ITS/MakeITSPlaneEfficiencySPD.C
@@ -1,34 +1,58 @@
-void MakeITSPlaneEfficiencySPD(Int_t firstRun=0,Int_t lastRun=AliCDBRunRange::Infinity()){
+#if !defined(__CINT__) || defined(__MAKECINT__)
+#include <TRandom3.h>
+#include "AliCDBManager.h"
+#include "AliCDBMetaData.h"
+#include "AliCDBStorage.h"
+#include "AliCDBRunRange.h"
+#include "AliCDBId.h"
+#include "AliITSPlaneEffSPD.h"
+#endif
+
+void MakeITSPlaneEfficiencySPD(Int_t firstRun=0, Int_t lastRun=AliCDBRunRange::Infinity(), 
+ Double_t eff=0.99, Int_t nTried=1000){
   
+  if(eff<0 || eff > 1) {
+   printf("Efficiency must be in the range [0,1]: nothing done");
+  }
   if(!AliCDBManager::Instance()->IsDefaultStorageSet()) {
-    AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
+    AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
   }
   
   AliCDBMetaData *md1= new AliCDBMetaData(); // metaData describing the object
   md1->SetObjectClassName("AliITSPlaneEff");
   md1->SetResponsible("Giuseppe Bruno");
   md1->SetBeamPeriod(0);
-  md1->SetAliRootVersion("head 19/11/07"); //root version
+  md1->SetAliRootVersion("head 16/02/08"); //root version
 
   AliCDBId idplaneeffSPD("ITS/PlaneEff/PlaneEffSPD",firstRun, lastRun);
   
   AliITSPlaneEffSPD* planeeffSPD = new AliITSPlaneEffSPD();
+  TRandom3 *gran = new TRandom3();
 
 //  planeeffSPD->SetOwner(kFALSE);
   
-
+//  Int_t nTried=1000;
+  Double_t limit=nTried;
+  limit*=(1-eff);
+  printf("limit = %f",limit);
   // loop over SPD chip
   Bool_t BFound=kFALSE;
-  //for(Int_t key=0;key<1200;key++){
-  for(UInt_t mod=0;mod<240;mod++){
-  for(UInt_t chip=0;chip<5;chip++){
+  for(UInt_t key=0;key<planeeffSPD->Nblock();key++){
+  //for(UInt_t mod=0;mod<240;mod++){
+  //for(UInt_t chip=0;chip<5;chip++){
   // suppose to have 1000 tracks in each chip and an average efficiency of 99%
-    for(Int_t j=0; j<1000; j++) {
+    for(Int_t j=0; j<nTried; j++) {
       BFound=kFALSE;
-      if (gRandom->Uniform(0,1000)>10) BFound=kTRUE;
-      planeeffSPD->UpDatePlaneEff(BFound,mod,chip);
+      //if (gRandom->Uniform(0,1000)>10) BFound=kTRUE;
+      if (nTried*gran->Uniform()>limit) BFound=kTRUE;
+      //planeeffSPD->UpDatePlaneEff(BFound,mod,chip);
+      planeeffSPD->UpDatePlaneEff(BFound,key);
     }
-  }}
+  //}}
+  }
   if(AliCDBManager::Instance()->GetDefaultStorage()->Put(planeeffSPD, idplaneeffSPD, md1))
-  {cout << "Riuscito " << endl;}
-}
\ No newline at end of file
+  printf("Local CDB file with random SPD plane efficiencies written \n");
+ delete gran;
+ delete planeeffSPD;
+ delete md1;
+}