From 27e9e1162e567aadf47e5a540752faf67e11b4cb Mon Sep 17 00:00:00 2001
From: Morfouace <pierre.morfouace@gmail.com>
Date: Thu, 23 Feb 2023 13:50:52 +0100
Subject: [PATCH] Updating Example4 analysis SetNuclei3 -> Setparticle3
---
Examples/Example4/Analysis.cxx | 436 +++++++++++-----------
Examples/Example4/configs/ConfigActar.dat | 2 +-
2 files changed, 219 insertions(+), 219 deletions(-)
diff --git a/Examples/Example4/Analysis.cxx b/Examples/Example4/Analysis.cxx
index a9a9a1ab0..0e96776b9 100644
--- a/Examples/Example4/Analysis.cxx
+++ b/Examples/Example4/Analysis.cxx
@@ -38,187 +38,187 @@ Analysis::~Analysis(){
////////////////////////////////////////////////////////////////////////////////
void Analysis::Init(){
- Actar= (TActarPhysics*) m_DetectorManager->GetDetector("Actar");
- ReactionConditions = new TReactionConditions();
-
- Actar->ReadAnalysisConfig();
- if(Actar->GetRansacStatus()){
- Actar->SetRansacParameter("./configs/RansacConfig.dat");
- }
- else if(Actar->GetClusterStatus()){
- Actar->SetClusterParameter("./configs/ClusterConfig.dat");
- }
-
- DriftVelocity = Actar->GetDriftVelocity();
- PadSizeX = Actar->GetPadSizeX();
- PadSizeY = Actar->GetPadSizeY();
- NumberOfPadsX = Actar->GetNumberOfPadsX();
- NumberOfPadsY = Actar->GetNumberOfPadsY();
-
- EnergyLoss_1H = NPL::EnergyLoss("./EnergyLossTable/proton_iC4H10_6.24151e+07_295.G4table","G4Table",100);
- EnergyLoss_18O = NPL::EnergyLoss("./EnergyLossTable/O18_iC4H10_6.24151e+07_295.G4table","G4Table",100);
- TheReaction = new NPL::Reaction("18O(p,p)18O@59");
-
- InitInputBranch();
- InitOutputBranch();
+ Actar= (TActarPhysics*) m_DetectorManager->GetDetector("Actar");
+ ReactionConditions = new TReactionConditions();
+
+ Actar->ReadAnalysisConfig();
+ if(Actar->GetRansacStatus()){
+ Actar->SetRansacParameter("./configs/RansacConfig.dat");
+ }
+ else if(Actar->GetClusterStatus()){
+ Actar->SetClusterParameter("./configs/ClusterConfig.dat");
+ }
+
+ DriftVelocity = Actar->GetDriftVelocity();
+ PadSizeX = Actar->GetPadSizeX();
+ PadSizeY = Actar->GetPadSizeY();
+ NumberOfPadsX = Actar->GetNumberOfPadsX();
+ NumberOfPadsY = Actar->GetNumberOfPadsY();
+
+ EnergyLoss_1H = NPL::EnergyLoss("./EnergyLossTable/proton_iC4H10_6.24151e+07_295.G4table","G4Table",100);
+ EnergyLoss_18O = NPL::EnergyLoss("./EnergyLossTable/O18_iC4H10_6.24151e+07_295.G4table","G4Table",100);
+ TheReaction = new NPL::Reaction("18O(p,p)18O@59");
+
+ InitInputBranch();
+ InitOutputBranch();
}
////////////////////////////////////////////////////////////////////////////////
void Analysis::TreatEvent(){
- ReInitValue();
-
- LightName="";
- if(ReactionConditions->GetParticleMultiplicity()>0){
- InitXVertex = ReactionConditions->GetVertexPositionZ();
- InitTheta3 = ReactionConditions->GetTheta(0);
- InitE3 = ReactionConditions->GetKineticEnergy(0);
- LightName = ReactionConditions->GetParticleName(0);
+ ReInitValue();
+
+ LightName="";
+ if(ReactionConditions->GetParticleMultiplicity()>0){
+ InitXVertex = ReactionConditions->GetVertexPositionZ();
+ InitTheta3 = ReactionConditions->GetTheta(0);
+ InitE3 = ReactionConditions->GetKineticEnergy(0);
+ LightName = ReactionConditions->GetParticleName(0);
+ }
+ int TrackMult = Actar->GetTrackMult();
+
+ TVector3 vX = TVector3(1,0,0);
+ TVector3 aTrack, vB;
+ if(TrackMult>1){
+ vTrack = Actar->GetTracks();
+ double scalarproduct=0;
+ int BeamTrack=0;
+ for(unsigned int i=0; i<TrackMult; i++){
+ TVector3 vtest = TVector3(vTrack[i].GetDirectionVector().X(),vTrack[i].GetDirectionVector().Y(),vTrack[i].GetDirectionVector().Z());
+ TVector3 vunit = vtest.Unit();
+ double scalar = abs(vunit.Dot(vX));
+ vScalar.push_back(scalar);
+ //cout << scalar << endl;
+ //cout << scalarproduct << endl;
+ if(scalar>scalarproduct){
+ BeamTrack=i;
+ scalarproduct=scalar;
+ }
}
- int TrackMult = Actar->GetTrackMult();
-
- TVector3 vX = TVector3(1,0,0);
- TVector3 aTrack, vB;
- if(TrackMult>1){
- vTrack = Actar->GetTracks();
- double scalarproduct=0;
- int BeamTrack=0;
- for(unsigned int i=0; i<TrackMult; i++){
- TVector3 vtest = TVector3(vTrack[i].GetDirectionVector().X(),vTrack[i].GetDirectionVector().Y(),vTrack[i].GetDirectionVector().Z());
- TVector3 vunit = vtest.Unit();
- double scalar = abs(vunit.Dot(vX));
- vScalar.push_back(scalar);
- //cout << scalar << endl;
- //cout << scalarproduct << endl;
- if(scalar>scalarproduct){
- BeamTrack=i;
- scalarproduct=scalar;
- }
- }
-
- double XBeam = vTrack[BeamTrack].GetDirectionVector().X();
- double YBeam = vTrack[BeamTrack].GetDirectionVector().Y();
- double ZBeam = vTrack[BeamTrack].GetDirectionVector().Z();
- TVector3 vBeam = TVector3(XBeam,YBeam,ZBeam);
-
- double XBeamPoint = vTrack[BeamTrack].GetXh();
- double YBeamPoint = vTrack[BeamTrack].GetYh();
- double ZBeamPoint = vTrack[BeamTrack].GetZh();
- TVector3 vBeamPos = TVector3(XBeamPoint,YBeamPoint,ZBeamPoint);
-
- //vB = TVector3(XBeam*PadSizeX, YBeam*PadSizeY,ZBeam*DriftVelocity);
- vB = TVector3(XBeam*PadSizeX, YBeam*PadSizeY,ZBeam);
- BeamAngle = (vX.Angle(vB))*180/TMath::Pi();
-
- for(unsigned int i=0; i<TrackMult; i++){
- if(i!=BeamTrack){
- double Xdir = vTrack[i].GetDirectionVector().X();
- double Ydir = vTrack[i].GetDirectionVector().Y();
- double Zdir = vTrack[i].GetDirectionVector().Z();
-
- double vertex_x = vTrack[i].GetVertexPostion(vBeam,vBeamPos).X()*PadSizeX;
- double vertex_y = vTrack[i].GetVertexPostion(vBeam,vBeamPos).Y()*PadSizeY;
- //double vertex_z = vTrack[i].GetVertexPostion(vBeam,vBeamPos).Z()*DriftVelocity;
- double vertex_z = vTrack[i].GetVertexPostion(vBeam,vBeamPos).Z();
-
- //aTrack = TVector3(Xdir*PadSizeX, Ydir*PadSizeY, Zdir*DriftVelocity);
- aTrack = TVector3(Xdir*PadSizeX, Ydir*PadSizeY, Zdir);
- double angle = vX.Angle(aTrack)*180/TMath::Pi();
- //double angle = vB.Angle(aTrack)*180/TMath::Pi();
- if(angle>90) angle = 180-angle;
-
- double x1 = vTrack[i].GetXm()*PadSizeX;
- double x2 = vTrack[i].GetXh()*PadSizeX;
- double y1 = vTrack[i].GetYm()*PadSizeY-0.5*NumberOfPadsY*PadSizeY;
- double y2 = vTrack[i].GetYh()*PadSizeY-0.5*NumberOfPadsY*PadSizeY;
- //double z1 = -(vTrack[i].GetZm()-256)*DriftVelocity;
- //double z2 = -(vTrack[i].GetZh()-256)*DriftVelocity;
- //double z1 = vTrack[i].GetZm()*DriftVelocity;
- double z1 = vTrack[i].GetZm();
- //double z2 = vTrack[i].GetZh()*DriftVelocity;
- double z2 = vTrack[i].GetZh();
-
-
- if(vertex_x>0 && vertex_x<256){
- double LengthInGas = fSiDistanceX - vertex_x;
- for(unsigned int k=0; k<Actar->Si_E.size(); k++){
- XVertex.push_back(vertex_x);
- YVertex.push_back(vertex_y);
- ZVertex.push_back(vertex_z);
-
- GetMayaSiHitPosition(x1,x2,y1,y2,z1,z2);
- ESi.push_back(Actar->Si_E[k]);
- SiNumber.push_back(Actar->Si_Number[k]);
-
- DE.push_back(vTrack[i].GetPartialCharge(108,128)/(20./cos(angle*TMath::Pi()/180)));
- double E3;
-
- if(LightName=="proton")E3 = EnergyLoss_1H.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
- if(LightName=="deuteron")E3 = EnergyLoss_2H.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
- if(LightName=="triton")E3 = EnergyLoss_3H.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
- if(LightName=="He3")E3 = EnergyLoss_3He.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
- if(LightName=="alpha")E3 = EnergyLoss_4He.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
- //double BeamEnergy = EnergyLoss_18O.Slow(59.4*MeV,(vertex_x[i]+60)*mm, BeamAngle*TMath::Pi()/180);
- double Energy_beam = EnergyLoss_18O.Slow(59.4*MeV,(vertex_x+60)*mm, 0);
- BeamEnergy.push_back(Energy_beam);
- TheReaction->SetBeamEnergy(Energy_beam);
- ELab.push_back(E3);
- ThetaLab.push_back(angle);
- TheReaction->SetNuclei3(E3,angle*TMath::Pi()/180);
- Ex.push_back(TheReaction->GetExcitation4());
- ThetaCM.push_back(TheReaction->GetThetaCM()*180./TMath::Pi());
- }
- }
- }
+
+ double XBeam = vTrack[BeamTrack].GetDirectionVector().X();
+ double YBeam = vTrack[BeamTrack].GetDirectionVector().Y();
+ double ZBeam = vTrack[BeamTrack].GetDirectionVector().Z();
+ TVector3 vBeam = TVector3(XBeam,YBeam,ZBeam);
+
+ double XBeamPoint = vTrack[BeamTrack].GetXh();
+ double YBeamPoint = vTrack[BeamTrack].GetYh();
+ double ZBeamPoint = vTrack[BeamTrack].GetZh();
+ TVector3 vBeamPos = TVector3(XBeamPoint,YBeamPoint,ZBeamPoint);
+
+ //vB = TVector3(XBeam*PadSizeX, YBeam*PadSizeY,ZBeam*DriftVelocity);
+ vB = TVector3(XBeam*PadSizeX, YBeam*PadSizeY,ZBeam);
+ BeamAngle = (vX.Angle(vB))*180/TMath::Pi();
+
+ for(unsigned int i=0; i<TrackMult; i++){
+ if(i!=BeamTrack){
+ double Xdir = vTrack[i].GetDirectionVector().X();
+ double Ydir = vTrack[i].GetDirectionVector().Y();
+ double Zdir = vTrack[i].GetDirectionVector().Z();
+
+ double vertex_x = vTrack[i].GetVertexPostion(vBeam,vBeamPos).X()*PadSizeX;
+ double vertex_y = vTrack[i].GetVertexPostion(vBeam,vBeamPos).Y()*PadSizeY;
+ //double vertex_z = vTrack[i].GetVertexPostion(vBeam,vBeamPos).Z()*DriftVelocity;
+ double vertex_z = vTrack[i].GetVertexPostion(vBeam,vBeamPos).Z();
+
+ //aTrack = TVector3(Xdir*PadSizeX, Ydir*PadSizeY, Zdir*DriftVelocity);
+ aTrack = TVector3(Xdir*PadSizeX, Ydir*PadSizeY, Zdir);
+ double angle = vX.Angle(aTrack)*180/TMath::Pi();
+ //double angle = vB.Angle(aTrack)*180/TMath::Pi();
+ if(angle>90) angle = 180-angle;
+
+ double x1 = vTrack[i].GetXm()*PadSizeX;
+ double x2 = vTrack[i].GetXh()*PadSizeX;
+ double y1 = vTrack[i].GetYm()*PadSizeY-0.5*NumberOfPadsY*PadSizeY;
+ double y2 = vTrack[i].GetYh()*PadSizeY-0.5*NumberOfPadsY*PadSizeY;
+ //double z1 = -(vTrack[i].GetZm()-256)*DriftVelocity;
+ //double z2 = -(vTrack[i].GetZh()-256)*DriftVelocity;
+ //double z1 = vTrack[i].GetZm()*DriftVelocity;
+ double z1 = vTrack[i].GetZm();
+ //double z2 = vTrack[i].GetZh()*DriftVelocity;
+ double z2 = vTrack[i].GetZh();
+
+
+ if(vertex_x>0 && vertex_x<256){
+ double LengthInGas = fSiDistanceX - vertex_x;
+ for(unsigned int k=0; k<Actar->Si_E.size(); k++){
+ XVertex.push_back(vertex_x);
+ YVertex.push_back(vertex_y);
+ ZVertex.push_back(vertex_z);
+
+ GetMayaSiHitPosition(x1,x2,y1,y2,z1,z2);
+ ESi.push_back(Actar->Si_E[k]);
+ SiNumber.push_back(Actar->Si_Number[k]);
+
+ DE.push_back(vTrack[i].GetPartialCharge(108,128)/(20./cos(angle*TMath::Pi()/180)));
+ double E3;
+
+ if(LightName=="proton")E3 = EnergyLoss_1H.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
+ if(LightName=="deuteron")E3 = EnergyLoss_2H.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
+ if(LightName=="triton")E3 = EnergyLoss_3H.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
+ if(LightName=="He3")E3 = EnergyLoss_3He.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
+ if(LightName=="alpha")E3 = EnergyLoss_4He.EvaluateInitialEnergy(Actar->Si_E[k]*MeV,LengthInGas*mm,angle*TMath::Pi()/180);
+ //double BeamEnergy = EnergyLoss_18O.Slow(59.4*MeV,(vertex_x[i]+60)*mm, BeamAngle*TMath::Pi()/180);
+ double Energy_beam = EnergyLoss_18O.Slow(59.4*MeV,(vertex_x+60)*mm, 0);
+ BeamEnergy.push_back(Energy_beam);
+ TheReaction->SetBeamEnergy(Energy_beam);
+ ELab.push_back(E3);
+ ThetaLab.push_back(angle);
+ TheReaction->SetParticle3(E3,angle*TMath::Pi()/180);
+ Ex.push_back(TheReaction->GetExcitation4());
+ ThetaCM.push_back(TheReaction->GetThetaCM()*180./TMath::Pi());
+ }
}
+ }
}
+ }
}
////////////////////////////////////////////////////////////////////////////////
void Analysis::GetMayaSiHitPosition(double xm, double xh, double ym, double yh, double zm, double zh)
{
- double X1, X2, Y1, Y2, Z1, Z2;
-
- if(xm>xh){
- X1 = xh;
- Y1 = yh;
- Z1 = zh;
-
- X2 = xm;
- Y2 = ym;
- Z2 = zm;
- }
- else if(xh>xm){
- X1 = xm;
- Y1 = ym;
- Z1 = zm;
-
- X2 = xh;
- Y2 = yh;
- Z2 = zh;
- }
-
- double l, L, t;
- double zf, yf;
-
- if(fSiDistanceX>X2){
- L = fSiDistanceX-X2;
- l = X2 - X1;
- t = (l+L)/l;
-
- zf = Z1 + (Z2-Z1)*t;
- yf = Y1 + (Y2-Y1)*t;
- }
- else if(fSiDistanceX<X2){
- L = X2 - fSiDistanceX;
- l = fSiDistanceX - X1;
- t = (l+L)/l;
-
- zf = Z1 + (Z2-Z1)/t;
- yf = Y1 + (Y2-Y1)/t;
- }
-
- SiPosY.push_back(yf);
- SiPosZ.push_back(zf);
+ double X1, X2, Y1, Y2, Z1, Z2;
+
+ if(xm>xh){
+ X1 = xh;
+ Y1 = yh;
+ Z1 = zh;
+
+ X2 = xm;
+ Y2 = ym;
+ Z2 = zm;
+ }
+ else if(xh>xm){
+ X1 = xm;
+ Y1 = ym;
+ Z1 = zm;
+
+ X2 = xh;
+ Y2 = yh;
+ Z2 = zh;
+ }
+
+ double l, L, t;
+ double zf, yf;
+
+ if(fSiDistanceX>X2){
+ L = fSiDistanceX-X2;
+ l = X2 - X1;
+ t = (l+L)/l;
+
+ zf = Z1 + (Z2-Z1)*t;
+ yf = Y1 + (Y2-Y1)*t;
+ }
+ else if(fSiDistanceX<X2){
+ L = X2 - fSiDistanceX;
+ l = fSiDistanceX - X1;
+ t = (l+L)/l;
+
+ zf = Z1 + (Z2-Z1)/t;
+ yf = Y1 + (Y2-Y1)/t;
+ }
+
+ SiPosY.push_back(yf);
+ SiPosZ.push_back(zf);
}
////////////////////////////////////////////////////////////////////////////////
@@ -227,73 +227,73 @@ void Analysis::End(){
////////////////////////////////////////////////////////////////////////////////
void Analysis::InitInputBranch() {
- RootInput::getInstance()->GetChain()->SetBranchStatus("ReactionConditions",true);
- RootInput::getInstance()->GetChain()->SetBranchStatus("fRC_*",true);
- RootInput::getInstance()->GetChain()->SetBranchAddress("ReactionConditions",&ReactionConditions);
+ RootInput::getInstance()->GetChain()->SetBranchStatus("ReactionConditions",true);
+ RootInput::getInstance()->GetChain()->SetBranchStatus("fRC_*",true);
+ RootInput::getInstance()->GetChain()->SetBranchAddress("ReactionConditions",&ReactionConditions);
}
////////////////////////////////////////////////////////////////////////////////
void Analysis::InitOutputBranch() {
- RootOutput::getInstance()->GetTree()->Branch("DE",&DE);
- RootOutput::getInstance()->GetTree()->Branch("SiNumber",&SiNumber);
- RootOutput::getInstance()->GetTree()->Branch("ESi",&ESi);
- RootOutput::getInstance()->GetTree()->Branch("ELab",&ELab);
- RootOutput::getInstance()->GetTree()->Branch("ThetaLab",&ThetaLab);
- RootOutput::getInstance()->GetTree()->Branch("Ex",&Ex);
- RootOutput::getInstance()->GetTree()->Branch("ThetaCM",&ThetaCM);
- RootOutput::getInstance()->GetTree()->Branch("vScalar",&vScalar);
- RootOutput::getInstance()->GetTree()->Branch("XVertex",&XVertex);
- RootOutput::getInstance()->GetTree()->Branch("YVertex",&YVertex);
- RootOutput::getInstance()->GetTree()->Branch("ZVertex",&ZVertex);
- RootOutput::getInstance()->GetTree()->Branch("BeamAngle",&BeamAngle,"BeamAngle/D");
- RootOutput::getInstance()->GetTree()->Branch("BeamEnergy",&BeamEnergy);
- RootOutput::getInstance()->GetTree()->Branch("SiPosY",&SiPosY);
- RootOutput::getInstance()->GetTree()->Branch("SiPosZ",&SiPosZ);
- RootOutput::getInstance()->GetTree()->Branch("InitXVertex",&InitXVertex,"InitXVertex/D");
- RootOutput::getInstance()->GetTree()->Branch("InitE3",&InitE3,"InitE3/D");
- RootOutput::getInstance()->GetTree()->Branch("InitTheta3",&InitTheta3,"InitTheta3/D");
-
+ RootOutput::getInstance()->GetTree()->Branch("DE",&DE);
+ RootOutput::getInstance()->GetTree()->Branch("SiNumber",&SiNumber);
+ RootOutput::getInstance()->GetTree()->Branch("ESi",&ESi);
+ RootOutput::getInstance()->GetTree()->Branch("ELab",&ELab);
+ RootOutput::getInstance()->GetTree()->Branch("ThetaLab",&ThetaLab);
+ RootOutput::getInstance()->GetTree()->Branch("Ex",&Ex);
+ RootOutput::getInstance()->GetTree()->Branch("ThetaCM",&ThetaCM);
+ RootOutput::getInstance()->GetTree()->Branch("vScalar",&vScalar);
+ RootOutput::getInstance()->GetTree()->Branch("XVertex",&XVertex);
+ RootOutput::getInstance()->GetTree()->Branch("YVertex",&YVertex);
+ RootOutput::getInstance()->GetTree()->Branch("ZVertex",&ZVertex);
+ RootOutput::getInstance()->GetTree()->Branch("BeamAngle",&BeamAngle,"BeamAngle/D");
+ RootOutput::getInstance()->GetTree()->Branch("BeamEnergy",&BeamEnergy);
+ RootOutput::getInstance()->GetTree()->Branch("SiPosY",&SiPosY);
+ RootOutput::getInstance()->GetTree()->Branch("SiPosZ",&SiPosZ);
+ RootOutput::getInstance()->GetTree()->Branch("InitXVertex",&InitXVertex,"InitXVertex/D");
+ RootOutput::getInstance()->GetTree()->Branch("InitE3",&InitE3,"InitE3/D");
+ RootOutput::getInstance()->GetTree()->Branch("InitTheta3",&InitTheta3,"InitTheta3/D");
+
}
////////////////////////////////////////////////////////////////////////////////
void Analysis::ReInitValue(){
- DE.clear();
- SiNumber.clear();
- ESi.clear();
- ELab.clear();
- ThetaLab.clear();
- Ex.clear();
- ThetaCM.clear();
- vScalar.clear();
- XVertex.clear();
- YVertex.clear();
- ZVertex.clear();
- SiPosY.clear();
- SiPosZ.clear();
- BeamEnergy.clear();
-
- BeamAngle=-1000;
- InitE3=-1000;
- InitTheta3=-1000;
- InitXVertex=-1000;
+ DE.clear();
+ SiNumber.clear();
+ ESi.clear();
+ ELab.clear();
+ ThetaLab.clear();
+ Ex.clear();
+ ThetaCM.clear();
+ vScalar.clear();
+ XVertex.clear();
+ YVertex.clear();
+ ZVertex.clear();
+ SiPosY.clear();
+ SiPosZ.clear();
+ BeamEnergy.clear();
+
+ BeamAngle=-1000;
+ InitE3=-1000;
+ InitTheta3=-1000;
+ InitXVertex=-1000;
}
////////////////////////////////////////////////////////////////////////////////
// Construct Method to be pass to the DetectorFactory //
////////////////////////////////////////////////////////////////////////////////
NPL::VAnalysis* Analysis::Construct(){
- return (NPL::VAnalysis*) new Analysis();
+ return (NPL::VAnalysis*) new Analysis();
}
////////////////////////////////////////////////////////////////////////////////
// Registering the construct method to the factory //
////////////////////////////////////////////////////////////////////////////////
extern "C"{
- class proxy{
+ class proxy{
public:
- proxy(){
- NPL::AnalysisFactory::getInstance()->SetConstructor(Analysis::Construct);
- }
- };
-
- proxy p;
+ proxy(){
+ NPL::AnalysisFactory::getInstance()->SetConstructor(Analysis::Construct);
+ }
+ };
+
+ proxy p;
}
diff --git a/Examples/Example4/configs/ConfigActar.dat b/Examples/Example4/configs/ConfigActar.dat
index b1cbe69c7..c6061f287 100644
--- a/Examples/Example4/configs/ConfigActar.dat
+++ b/Examples/Example4/configs/ConfigActar.dat
@@ -1,7 +1,7 @@
ConfigActar
RecoRansac= 1
RecoCluster= 0
-RecoVisu= 1
+RecoVisu= 0
HIT_THRESHOLD= 2
Q_THRESHOLD= 0
T_THRESHOLD= 0
--
GitLab