Commit 9f2ad8f2 authored by Emilio Depero's avatar Emilio Depero
Browse files

updated script

parent c6777935
......@@ -4,7 +4,7 @@ INCLUDE := include/
DEP := include/Micromega.cc include/MicromegasPlane.cc include/physicalStrip.cc include/stripsCluster.cc
LIBS := `root-config --cflags --glibs`
all: map_multiplex ConvertOutputMM analyze_PositionData reconstructAndPlot
all: map_multiplex ConvertOutputMM analyze_PositionData reconstructAndPlot PlotMacro
%: %.cc
$(CXX) -o $@ $< -I$(INCLUDE) $(DEP) $(FLAGS) $(LIBS)
......
......@@ -52,6 +52,9 @@ do
#run the reco
echo "$multiplexedfile not yet reconstructed, running reco"
$ANALEXE $recofileout $analfileout
$ANALEXE $recofileout $analfileout 45
$ANALEXE $recofileout $analfileout 90
$ANALEXE $recofileout $analfileout 135
fi
echo "FILE $file SCANNED"
......
......@@ -62,14 +62,15 @@ double Xrot(const double X, const double Y, const double& thisangle)
{
const double Kx = cos(thisangle);
const double Ky = sin(thisangle);
return -Kx * X - Ky * Y;
return Kx * X + Ky * Y;
}
double Yrot(const double X, const double Y, const double& thisangle)
{
const double Kx = cos(thisangle);
const double Ky = sin(thisangle);
return Ky * X - Kx * Y;
//return -Ky * X + Kx * Y;
return Ky * X - Kx * Y; //to match Jonas geometry
}
// the format to store the sources
......@@ -220,7 +221,10 @@ void Analyze(int argc, char *argv[])
TH1D* decayPositions_z = new TH1D("decayPositions_z","reconstructed decay positions in z direction; z [mm] ",100,-32,32);
TH1D* decayPositions_z_target = new TH1D("decayPositions_z_target","reconstructed decay positions in z direction for -10<x1,x2,x3<10; z [mm] ",100,-32,32);
TH1D* decayPositions_z_20deg = new TH1D("decayPositions_z_20deg","reconstructed decay positions in z direction for -10<x<10, angle<20deg; z [mm] ",120,-32,32);
TH1D* decayPositions_z_20deg = new TH1D("decayPositions_z_20deg","reconstructed decay positions in z direction, angle<20deg; z [mm] ",120,-32,32);
TH1D* decayPositions_z_17deg = new TH1D("decayPositions_z_17deg","reconstructed decay positions in z direction, angle<17deg; z [mm] ",120,-32,32);
TH1D* decayPositions_z_13deg = new TH1D("decayPositions_z_13deg","reconstructed decay positions in z direction, angle<13deg; z [mm] ",120,-32,32);
TH1D* decayPositions_z_10deg = new TH1D("decayPositions_z_10deg","reconstructed decay positions in z direction, angle<10deg; z [mm] ",120,-32,32);
TH1D* qMax_All = new TH1D("qMax_All","distribution of (average) qmax; qmax ",200,0,1.4e6);
TH1D* qMax_double = new TH1D("qMax_double","distribution of average qmax (only double hits); qmax ",200,0,1.4e6);
......@@ -374,9 +378,21 @@ void Analyze(int argc, char *argv[])
}
if (TMath::Abs(Delta_y)/dist_MM12 < TMath::Tan(maxAngle_rad[0]))
{
if (std::abs(reconstructed_x)<10) decayPositions_z_20deg->Fill(reconstructed_y);
}
{
decayPositions_z_20deg->Fill(reconstructed_y);
}
if (TMath::Abs(Delta_y)/dist_MM12 < TMath::Tan(maxAngle_rad[1]))
{
decayPositions_z_17deg->Fill(reconstructed_y);
}
if (TMath::Abs(Delta_y)/dist_MM12 < TMath::Tan(maxAngle_rad[2]))
{
decayPositions_z_13deg->Fill(reconstructed_y);
}
if (TMath::Abs(Delta_y)/dist_MM12 < TMath::Tan(maxAngle_rad[3]))
{
decayPositions_z_10deg->Fill(reconstructed_y);
}
}
else
{
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment