Added class for map configuration

Author: contuandrea

include/analysis.h

@@ -45,51 +45,7 @@
 
 using namespace hydra::placeholders;
 namespace analysis{
-    template <typename Data>
-    TH3D *getHist3DDraw(Data mapx, Data mapy, Data mapz, std::set<double> _vx, std::set<double> _vy, std::set<double> _vz){
-        
-        std::vector<double> vx(_vx.begin(),_vx.end());
-        std::vector<double> vy(_vy.begin(),_vy.end());
-        std::vector<double> vz(_vz.begin(),_vz.end());
-        std::vector<double> vex(vx.size()+1);
-        std::vector<double> vey(vy.size()+1);
-        std::vector<double> vez(vz.size()+1);
-        
-        vex[0]=vx[0]-fabs(vx[0]-vx[1])/10.;
-        vey[0]=vy[0]-fabs(vy[0]-vy[1])/10.;
-        vez[0]=vz[0]-fabs(vz[0]-vz[1])/10.;
-        
-        for(size_t ix=0;ix<vx.size()-1;ix++){
-            vex[ix+1]=(vx[ix]+vx[ix+1])/2.;
-            
-        }
-        for(size_t iy=0;iy<vy.size()-1;iy++){
-            vey[iy+1]=(vy[iy]+vy[iy+1])/2.;
-            
-        }
-        for(size_t iz=0;iz<vz.size()-1;iz++){
-            vez[iz+1]=(vz[iz]+vz[iz+1])/2.;
-        }
-        
-        vex[vx.size()]=vx[vx.size()-1]+fabs(vx[vx.size()-2]-vx[vx.size()-1])/10.;
-        vey[vy.size()]=vy[vy.size()-1]+fabs(vy[vy.size()-2]-vy[vy.size()-1])/10.;
-        vez[vz.size()]=vz[vz.size()-1]+fabs(vz[vz.size()-2]-vz[vz.size()-1])/10.;
-        
-        TH3D *h3=new TH3D("physmap","physmap;x [#mum];y [#mum];z [#mum]",vx.size(),vex.data(),vy.size(),vey.data(),vz.size(),vez.data());
-        
-
-        size_t count=0;
-        for(size_t iz=0;iz<vz.size();iz++){
-            for(size_t ix=0;ix<vx.size();ix++){
-                for(size_t iy=0;iy<vy.size();iy++){
-                    h3->SetBinContent(ix+1,iy+1,iz+1, sqrt(mapx[count]*mapx[count]+mapy[count]*mapy[count]+mapz[count]*mapz[count]));
-                    count++;
-                }
-            }
-        }
-        return h3;
-    }
-
+    
     TGraph *setGraph(size_t points, std::string nm, int color, int mstyle=8, int msize=1, int lwidth=2){
         TGraph *g = new TGraph(points);
         g->SetName(nm.c_str());
@@ -736,7 +692,7 @@ namespace analysis{
             }
 
 
-            
+            extradir->WriteTObject(hdraw,"h3bkg","Overwrite");
 
             can->Print(Form("%s/%s.gif++",settings["outputdir"].c_str(),settings["name"].c_str()));
             delete can;

include/details/utils.h

@@ -81,6 +81,24 @@ if(flag){\
 #include <hydra/Zip.h>
 // #include <hydra/detail/external/Eigen/Dense>
 
+
+//ROOT
+#ifdef _ROOT_AVAILABLE_
+#include "TCanvas.h"
+#include "TROOT.h"
+#include "TGraph.h"
+#include "TGraph2D.h"
+#include "TH3D.h"
+#include "TH2D.h"
+#include "TApplication.h"
+#include "TGraphTime.h"
+#include "TMarker.h"
+#include "TFile.h"
+#include "TPad.h"
+#include "TLegend.h"
+#include "TTree.h"
+#endif //_ROOT_AVAILABLE_
+
 //printouts
 #define TOPROW    "===================================================================================================="
 #define MIDDLEROW "----------------------------------------------------------------------------------------------------"
@@ -425,4 +443,16 @@ int sgn(T val) {
     return (T(0) < val) - (val < T(0));
 }
 
+struct vectorLength
+{
+    
+    template<typename Particle>
+    __hydra_dual__
+    double operator()(Particle p){
+        return sqrt(hydra::get<0>(p)*hydra::get<0>(p) + hydra::get<1>(p)*hydra::get<1>(p) +hydra::get<2>(p)*hydra::get<2>(p));
+                                                                                
+    }
+    
+};
+
 #endif /* UTILS_H_ */