Fixes for centrality < 0, adding jet finder setting with higher pT cut
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 30 Mar 2011 13:14:09 +0000 (13:14 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 30 Mar 2011 13:14:09 +0000 (13:14 +0000)
PWG4/JetTasks/AliAnalysisTaskJetServices.cxx
PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx
PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.h
PWG4/macros/AddTaskFragmentationFunction.C
PWG4/macros/AnalysisTrainPWG4Jets.C
PWG4/macros/jetspectrum/PlotNote.C
PWG4/macros/jetspectrum/rootlogon.C

index da221e0..f7d67c6 100644 (file)
@@ -455,7 +455,7 @@ void AliAnalysisTaskJetServices::UserExec(Option_t */*option*/)
       }
     }
 
-    if(cand&&esdVtxIn){
+    if(cand&&esdVtxIn&&iCl<5){
       fh2ESDTriggerCount->Fill(0.,kSelectedALICEVertexIn);
       fh2ESDTriggerCount->Fill(iCl,kSelectedALICEVertexIn);
       fh2ESDTriggerVtx->Fill(kSelectedALICEVertexIn*(iCl+1),zvtx);
@@ -505,7 +505,7 @@ void AliAnalysisTaskJetServices::UserExec(Option_t */*option*/)
        fh2TriggerVtx->Fill(kSelectedALICEVertexValid*(iCl+1),zvtx);
       }
     }
-    if(cand&&aodVtxIn){
+    if(cand&&aodVtxIn&&iCl<5){
       fh2TriggerCount->Fill(0.,kSelectedALICEVertexIn);
       fh2TriggerCount->Fill(iCl,kSelectedALICEVertexIn);
       fh2TriggerVtx->Fill(kSelectedALICEVertexIn*(iCl+1),zvtx);
index addc41f..1f2ee0c 100644 (file)
@@ -125,6 +125,7 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2():
   fh2RPCosDeltaRP(0x0),
   fh2PtFGen(0x0),
   fh2RelPtFGen(0x0),
+  fh3PhiWeights(0x0),
   fh3RPPhiTracks(0x0),
   fHistList(0x0)  
 {
@@ -226,6 +227,7 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(const char* name):
   fh2RPCosDeltaRP(0x0),
   fh2PtFGen(0x0),
   fh2RelPtFGen(0x0),
+  fh3PhiWeights(0x0),
   fh3RPPhiTracks(0x0),
   fHistList(0x0)
 {
@@ -1502,8 +1504,8 @@ Bool_t AliAnalysisTaskJetSpectrum2::CalculateReactionPlaneAngle(const TList *tra
     fh3RPPhiTracks->Fill(fCentrality,momentum,track->Phi());
     count[0]++;
 
-    //    Double_t phiweight=GetPhiWeight(track->Phi(),momentum);
-    Double_t phiweight=1; 
+    Double_t phiweight=GetPhiWeight(track->Phi(),momentum);
+    //    Double_t phiweight=1; 
     Double_t weight=2;
     if(track->Pt()<2){weight=track->Pt();}
     
@@ -1595,4 +1597,11 @@ Int_t AliAnalysisTaskJetSpectrum2::GetPhiBin(Double_t phi)
     return phibin;
 }
 
+
+
+Double_t AliAnalysisTaskJetSpectrum2::GetPhiWeight(Double_t phi,Double_t signedpt){
+  if(!fh3PhiWeights)return 1;
+  else return fh3PhiWeights->GetBinContent(fh3PhiWeights->GetXaxis()->FindBin(fCentrality),fh3PhiWeights->GetYaxis()->FindBin(signedpt),fh3PhiWeights->GetZaxis()->FindBin(phi));
+}
+
  //________________________________________________________________________
index f916e9f..3b0643f 100644 (file)
@@ -72,6 +72,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     virtual void SetTrackTypeRec(Int_t i){fTrackTypeRec = i;}
     virtual void SetFilterMask(UInt_t i){fFilterMask = i;}
     virtual void SetEventSelectionMask(UInt_t i){fEventSelectionMask = i;}
+    virtual void SetPhiWeights(TH3F *phiw){fh3PhiWeights = phiw;}
     virtual void SetNonStdFile(char* c){fNonStdFile = c;} 
 
 
@@ -108,6 +109,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     Float_t GetCentrality();
     Bool_t  CalculateReactionPlaneAngle(const TList *trackList);
     Int_t   GetPhiBin(Double_t phi);
+    Double_t GetPhiWeight(Double_t phi,Double_t signedpt);
     Int_t   GetListOfJets(TList *list,TClonesArray* jarray,Int_t type);
     void    FillJetHistos(TList &jetsList,TList &particlesList,Int_t iType);
 
@@ -187,6 +189,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TH2F*         fh2PtFGen;                //! found vs generated 
     TH2F*         fh2RelPtFGen;             //! relative difference between generated and found 
 
+    TH3F*         fh3PhiWeights;  // RP phi weights, need to be set externally
     TH3F*         fh3RPPhiTracks; //! RP angle
     
 
index bf62ba4..252415f 100644 (file)
@@ -233,8 +233,6 @@ AliAnalysisTaskFragmentationFunction *AddTaskFragmentationFunction(
    task->SetKindSlices();      // default: kindSlice = 1 (inv mass)
    task->SetFFRadius();        // default: R = 0.4
    task->SetFFBckgRadius();    // default: R = 0.7
-   task->SetBckgMode();        // default: bgMode = 1
-   task->SetBckgType();        // default: 0,1,2
    task->SetBckgSubMethod();   // default: subMethod = O, 1 = leading jet removed for rho extraction, 2 = 2 leading jets removed
    task->SetIJMode(0);          // default: ijMode = 1
    task->SetQAMode();          // default: qaMode = 3
@@ -245,6 +243,15 @@ AliAnalysisTaskFragmentationFunction *AddTaskFragmentationFunction(
    task->SetHighPtThreshold(); // default: pt > 5 Gev
    task->UseRecEffRecJetPtBins(); // efficiency in bins of rec/gen jet pt - default: kTRUE  
 
+   task->SetBckgMode(1);        // default: bgMode = 1 
+   task->SetBckgType(AliAnalysisTaskFragmentationFunction::kBckgPerp,
+                    AliAnalysisTaskFragmentationFunction::kBckgOutLJStat,
+                    AliAnalysisTaskFragmentationFunction::kBckgClustersOutLeading,
+                    AliAnalysisTaskFragmentationFunction::kBckgClusters);  
+   task->SetBranchRecBackClusters("clustersAOD_KT04_B0_Filter00256_Cut00150_Skip00");
+
+
+
    // Define histo bins
    task->SetFFHistoBins();
    task->SetQAJetHistoBins();
index 71928b8..9f22916 100644 (file)
@@ -84,8 +84,10 @@ Int_t       iJETAN             = 1;      // Jet analysis (PWG4) // 1 write stand
 Int_t       iJETSUBTRACT        = 1;      // Jet background subtration
 TString     kDefaultJetBranch     = "";      // is currently set when filled (iJETAN or clusters) or from config macro 
 TString     kDefaultJetBackgroundBranch            = "";      // is currently set when filled (jet clsuters  
+TString     kDefaultJetBackgroundBranchCut1        = "";      // is currently set when filled (jet clsuters  
 TString     kDefaultJetBackgroundBranch_extra     = "";      // is currently set when filled (jet clsuters) 
 TString     kJetSubtractBranches     = "";      // is currently set when filled (jet clsuters  
+TString     kJetSubtractBranchesCut1     = "";      // is currently set when filled (jet clsuters  
 TString     kJetSubtractBranches_extra     = "";      // is currently set when filled (jet clsuters  
 
 TString     kDefaultJetBranchMC     = "";      // is currently set when filled (iJETAN or clusters) or from config macro 
@@ -432,6 +434,12 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
       gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/AddTaskESDFilterPWG4Train.C");
       AliAnalysisTaskESDfilter *taskesdfilter = AddTaskESDFilter(kUseKinefilter,kUseMuonfilter);
       taskesdfilter->SetEnableFillAOD(!kFilterAOD);
+      taskesdfilter->DisableV0s();
+      taskesdfilter->DisableCascades();
+      taskesdfilter->DisableKinks();
+      taskesdfilter->DisablePmdClusters();
+      taskesdfilter->DisableCaloClusters();
+      taskesdfilter->DisableCells();
 
       if(kIsMC){
        mgr->RegisterExtraFile("pyxsec_hists.root");
@@ -540,10 +548,19 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
         taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelectionFlag,"KT",0.4,0,1, kDeltaAODJetName.Data(),0.15,fTrackEtaWindow,0); // this one is for the background and random jets, random cones with no skip
         taskCl->SetBackgroundCalc(kTRUE);
         taskCl->SetNRandomCones(10);
+        taskCl->SetDebugLevel(11);
         taskCl->SetCentralityCut(fCenLo,fCenUp);
         taskCl->SetGhostEtamax(fTrackEtaWindow);
         kDefaultJetBackgroundBranch = Form("%s_%s",AliAODJetEventBackground::StdBranchName(),taskCl->GetJetOutputBranch());
-        
+
+
+        taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelectionFlag,"KT",0.4,0,1, kDeltaAODJetName.Data(),1.0,fTrackEtaWindow,0); // this one is for the background and random jets, random cones with no skip
+        taskCl->SetNRandomCones(10);
+        taskCl->SetBackgroundCalc(kTRUE);
+        taskCl->SetCentralityCut(fCenLo,fCenUp);
+        taskCl->SetGhostEtamax(fTrackEtaWindow);
+        kDefaultJetBackgroundBranchCut1 = Form("%s_%s",AliAODJetEventBackground::StdBranchName(),taskCl->GetJetOutputBranch());
+
          if (iPWG4FastEmbedding) {
            AliAnalysisTaskJetCluster *taskClEmb = 0;
            taskClEmb = AddTaskJetCluster("AODextra","",kHighPtFilterMask,iPhysicsSelectionFlag,"KT",0.4,0,1, kDeltaAODJetName.Data(),0.15,fTrackEtaWindow); // this one is for the background and random jets
@@ -589,6 +606,10 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
        kDefaultJetBranch = taskCl->GetJetOutputBranch();
        kJetSubtractBranches += Form("%s ",taskCl->GetJetOutputBranch());
 
+       taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelectionFlag,"ANTIKT",0.4,2,1,kDeltaAODJetName.Data(),1.0);
+       taskCl->SetCentralityCut(fCenLo,fCenUp);
+       if(kIsPbPb)taskCl->SetBackgroundBranch(kDefaultJetBackgroundBranchCut1.Data());
+       kJetSubtractBranchesCut1 += Form("%s ",taskCl->GetJetOutputBranch());
 
        taskCl = AddTaskJetCluster("AOD","",kHighPtFilterMask,iPhysicsSelectionFlag,"ANTIKT",0.2,0,1,kDeltaAODJetName.Data(),0.15);
        taskCl->SetCentralityCut(fCenLo,fCenUp);
@@ -666,6 +687,12 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
        if(kJetSubtractBranches.Contains(kDefaultJetBranch.Data())&&kIsPbPb){
         kDefaultJetBranch.ReplaceAll(taskSubtract->GetToReplace(),Form(taskSubtract->GetReplacementMask(),taskSubtract->GetSubtractionMethod()));
        }
+
+       // cut1
+       taskSubtract = AddTaskJetBackgroundSubtract(kJetSubtractBranchesCut1,1,kJetSubtractMask1.Data(),kJetSubtractMask2.Data(),"Cut1000");
+       taskSubtract->SetBackgroundBranch(kDefaultJetBackgroundBranchCut1.Data());              
+       taskSubtract->SelectCollisionCandidates(iPhysicsSelectionFlag);
+       if(kDeltaAODJetName.Length()>0)taskSubtract->SetNonStdOutputFile(kDeltaAODJetName.Data());
      }
      if(kJetSubtractBranches_extra.Length()){
        taskSubtract = AddTaskJetBackgroundSubtract(kJetSubtractBranches_extra,1,kJetSubtractMask1.Data(),kJetSubtractMask2.Data(),"extra");
@@ -733,7 +760,7 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
      }
    }
 
-if(iPWG4JetSpectrum){
+   if(iPWG4JetSpectrum){
      gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/AddTaskJetSpectrum2.C");
      AliAnalysisTaskJetSpectrum2 *taskjetSpectrum = 0;
      TString bkgClusters = kDefaultJetBackgroundBranch.Data(); 
@@ -755,10 +782,54 @@ if(iPWG4JetSpectrum){
           //      taskjetSpectrum->SetMinJetPt(10);
           taskjetSpectrum->SetTrackEtaWindow(fTrackEtaWindow);
           taskjetSpectrum->SetJetEtaWindow(fJetEtaWindow);
+          if(i!=1){
+            taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetRecFull,0);
+            taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetGenFull,0);
+          }
           if(iAODanalysis)SetAODInput(taskjetSpectrum);
 
-          
+          TString tmpCut1(kDefaultJetBranch.Data());
+          TString tmp2Cut1(kDefaultJetBranch.Data());
+          tmp2Cut1.ReplaceAll(Form(kJetSubtractMask2.Data(),1),Form(kJetSubtractMask1.Data(),0));
+          tmp2Cut1.ReplaceAll("00150","01000");
+          tmpCut1.ReplaceAll("00150","01000");
+
+       
+          taskjetSpectrum = AddTaskJetSpectrum2(tmpCut1.Data(),tmp2Cut1.Data(),kDeltaAODJetName.Data(),kHighPtFilterMask,AliVEvent::kMB,0,i);
+          taskjetSpectrum->SetBranchBkgRec(bkgClusters.Data());
+          taskjetSpectrum->SetBranchBkgGen(bkgClusters.Data());
+          taskjetSpectrum->SetTrackEtaWindow(fTrackEtaWindow);
+          taskjetSpectrum->SetJetEtaWindow(fJetEtaWindow);
+          if(i!=1){
+            taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetRecFull,0);
+            taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetGenFull,0);
+          }
+       
+          // 
+
+          // the random jets...
+          taskjetSpectrum = AddTaskJetSpectrum2("clustersAOD_KT04_B0_Filter00256_Cut00150_Skip00RandomConeSkip00",
+                                                "clustersAOD_KT04_B0_Filter00256_Cut00150_Skip00RandomCone_random",
+                                                kDeltaAODJetName.Data(),kHighPtFilterMask,AliVEvent::kMB,0,i);
+
+          taskjetSpectrum->SetBranchBkgRec(bkgClusters.Data());
+          taskjetSpectrum->SetBranchBkgGen(bkgClusters.Data());
+          taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetRecFull,0);
+          taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetGenFull,0);
+          taskjetSpectrum->SetDebugLevel(3);
+
+          taskjetSpectrum = AddTaskJetSpectrum2("clustersAOD_KT04_B0_Filter00256_Cut01000_Skip00RandomConeSkip00",
+                                                "clustersAOD_KT04_B0_Filter00256_Cut01000_Skip00RandomCone_random",
+                                                kDeltaAODJetName.Data(),kHighPtFilterMask,AliVEvent::kMB,0,i);
+          taskjetSpectrum->SetBranchBkgRec(bkgClusters.Data());
+          taskjetSpectrum->SetBranchBkgGen(bkgClusters.Data());
+          taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetRecFull,0);
+          taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetGenFull,0);
+
+
+
           // check the old subtracted vs. the new subtracted
+          /*
           TString tmp2(kDefaultJetBranch.Data());
           tmp2.ReplaceAll(Form(kJetSubtractMask2.Data(),1),Form(kJetSubtractMask2.Data(),2));
           taskjetSpectrum = AddTaskJetSpectrum2(tmp2.Data(),kDefaultJetBranch.Data(),kDeltaAODJetName.Data(),kHighPtFilterMask,AliVEvent::kMB,0,i);
@@ -769,8 +840,9 @@ if(iPWG4JetSpectrum){
           taskjetSpectrum->SetTrackEtaWindow(fTrackEtaWindow);
           taskjetSpectrum->SetJetEtaWindow(fJetEtaWindow);
           if(iAODanalysis)SetAODInput(taskjetSpectrum);
-
+          */
           // check the old subtracted vs. the new subtracted
+          /*
           TString tmp3(kDefaultJetBranch.Data());
           tmp3.ReplaceAll(Form(kJetSubtractMask2.Data(),1),Form(kJetSubtractMask2.Data(),3));
 
@@ -782,7 +854,7 @@ if(iPWG4JetSpectrum){
           taskjetSpectrum->SetTrackEtaWindow(0.8);
           taskjetSpectrum->SetJetEtaWindow(0.4);
           if(iAODanalysis)SetAODInput(taskjetSpectrum);
-
+          */
 
           taskjetSpectrum = AddTaskJetSpectrum2(Form("jetsAOD_UA104_B2_Filter%05d_Cut01000",kHighPtFilterMask),kDefaultJetBranch.Data(),kDeltaAODJetName.Data(),kHighPtFilterMask,AliVEvent::kMB,0,i);
           //      taskjetSpectrum->SetDebugLevel(3);
@@ -791,6 +863,10 @@ if(iPWG4JetSpectrum){
           taskjetSpectrum->SetJetEtaWindow(fJetEtaWindow);
           taskjetSpectrum->SetBranchBkgRec(bkgClusters.Data());
           taskjetSpectrum->SetBranchBkgGen(bkgClusters.Data());
+          if(i!=1){
+            taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetRecFull,0);
+            taskjetSpectrum->SetFlagJetType(AliAnalysisTaskJetSpectrum2::kJetGenFull,0);
+          }
           if(iAODanalysis)SetAODInput(taskjetSpectrum);
 
         }
index 6921158..4e590b3 100644 (file)
@@ -37,6 +37,8 @@ using namespace std;
 Int_t PlotSpectraPbPb();
 Int_t PlotJetBFluctuations();
 Int_t PlotJetQA();
+Int_t PlotPlay();
+
 // helpers
 
 
@@ -49,6 +51,10 @@ void SetStyleH1(TH1 *h,Int_t color = 0,Int_t fillStyle = 0,Int_t markerStyle = 0
 void ScaleH2(TH2* h,char* cX = "",char* cY = "",char* cZ = "",Float_t fScale = 1,Bool_t bWidth = true);
 TH1F* GetPythia8Spectrum(const char *cHist,const char *cPythia,Float_t deta,Int_t iRebin = 1,Int_t iMask = 0);
 TF1* FitLHSgaus(TH1D *hDeltaPt, double &sigma, double &sigma_error, double &mean);
+Double_t Gamma(Double_t *x,Double_t *par);
+Double_t Gamma0(Double_t *x,Double_t *par);
+
+
 const char *cPrintMask = "110116_%s.eps";
 void set_plot_style();
 TCanvas *cTmp = 0;
@@ -107,8 +113,8 @@ TH1F* GetPythia8Spectrum(const char *cHist,const char *cPythia,Float_t deta,Int_
 
 }
 
+
 void PlotNote(){
   PlotJetBFluctuations();
 }
 
@@ -1020,7 +1026,8 @@ Int_t PlotJetBFluctuations(){
   Double_t sigma = 0;
   Double_t sigma_err = 0;
   TF1* tmpGaus = 0;
-
+  TF1 *gamma0 = new TF1("gamma0",Gamma0,-40,40,4);
+  gamma0->SetParameters(1,1,1,1);
 
   for(int ic = 0;ic < nCen;ic++){
     TLegend *leg1 = new TLegend(0.2,0.65,0.3,0.93);
@@ -1036,7 +1043,7 @@ Int_t PlotJetBFluctuations(){
     leg1->AddEntry(hBiaL[ic],Form("BiA anti-k_{T}"),"P");
     */
 
-    
+    /*
     hBiaRC2[ic]->DrawCopy("psame");
     tmpGaus = FitLHSgaus(hBiaRC2[ic],sigma,sigma_err,mean);
     tmpGaus->SetRange(-40,40);
@@ -1045,8 +1052,15 @@ Int_t PlotJetBFluctuations(){
     tmpGaus->Draw("same");
     leg1->AddEntry(hBiaRC2[ic],sBiaRC2.Data(),"P");
     leg1->AddEntry(tmpGaus,Form("LHS Gaus fit: #mu = %2.2f, #sigma = %2.2f",mean,sigma),"L");
-    
+    */    
+
     hBiaRC[ic]->DrawCopy("psame");
+    hBiaRC[ic]->Fit(gamma0);
+    c1->Update();
+    if(getchar()=='q')return 1;
+    continue;
+
+
     tmpGaus = FitLHSgaus(hBiaRC[ic],sigma,sigma_err,mean);
     tmpGaus->SetRange(-40,40);
     tmpGaus->SetLineColor( hBiaRC[ic]->GetLineColor());
@@ -1282,6 +1296,40 @@ void ScaleNevent(TH1* h,TFile *fIn,Int_t iCond,Int_t iMC,Int_t iCen){
 
 }
 
+Double_t Gamma0(Double_t *x,Double_t *par){
+
+  // Fixed multiplicity M
+
+
+  Double_t p = par[0];
+  Double_t b = par[1];
+  Double_t A = par[2];
+  Double_t xval = x[0] + p/b;
+  // xval += par[3];
+  if(xval<0)return 0;
+
+  // take log to avoid zeros and NANs
+  Double_t f1 = TMath::Log(A*b);
+  Double_t f11 = ROOT::Math::lgamma(p); 
+  Double_t f2 = TMath::Log(b * xval)*(p-1);
+  Double_t f3 = -1.*b*xval;
+  Double_t f = f1-f11+f2+f3;
+  f = TMath::Exp(f);
+  return f;
+}
+
+Double_t Gamma(Double_t *x,Double_t *par){
+
+  // Fixed multiplicity M
+  Double_t M = par[0];
+  // Paramter p = <p_T>^2/sigma_p_T^2 adn b = <p_T>/sigma_p_T^2
+  Double_t b = par[1]/par[2]/par[2];
+  Double_t p = par[1] * b;
+
+  Double_t f = M * b / ROOT::Math::tgamma(M*p) * TMath::Power(M * b * x[0],M*(p-1)) * TMath::Exp(M*b*x[0]); 
+
+}
+
 TF1* FitLHSgaus(TH1D *hDeltaPt, double &sigma, double &sigma_error, double &mean)
 {
 
index 3062005..50f89e4 100644 (file)
@@ -52,6 +52,7 @@ Printf("%s \n>>> Setting style",(char*)__FILE__);
  gSystem->Load("libVMC");
  gSystem->Load("libCGAL");
  gSystem->Load("libfastjet");
+
  gSystem->Load("libsiscone");
  gSystem->Load("libSISConePlugin");
  gSystem->Load("libSTEERBase");