diff --git a/Simulation/G4Extensions/Sleptons/Sleptons/G4SElectronMinus.hh b/Simulation/G4Extensions/Sleptons/Sleptons/G4SElectronMinus.hh
new file mode 100755
index 0000000000000000000000000000000000000000..1aefaa50fce083f9904cd9ae29afddd2fe5b5a3c
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/Sleptons/G4SElectronMinus.hh
@@ -0,0 +1,28 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef G4SElectronMinus_H
+#define G4SElectronMinus_H
+
+#include "globals.hh"
+#include "G4ios.hh"
+#include "G4ParticleDefinition.hh"
+
+class G4SElectronMinus: public G4ParticleDefinition
+{
+ // singleton implementation
+
+private:
+
+ static G4SElectronMinus* theInstance;
+ G4SElectronMinus(){}
+ ~G4SElectronMinus(){}
+
+public:
+
+ static G4SElectronMinus* Definition();
+
+};
+
+#endif
diff --git a/Simulation/G4Extensions/Sleptons/Sleptons/G4SElectronPlus.hh b/Simulation/G4Extensions/Sleptons/Sleptons/G4SElectronPlus.hh
new file mode 100755
index 0000000000000000000000000000000000000000..1c1969dda687712adc82599f0d81a61448c9c67f
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/Sleptons/G4SElectronPlus.hh
@@ -0,0 +1,28 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef G4SElectronPlus_H
+#define G4SElectronPlus_H
+
+#include "globals.hh"
+#include "G4ios.hh"
+#include "G4ParticleDefinition.hh"
+
+class G4SElectronPlus: public G4ParticleDefinition
+{
+ // singleton implementation
+
+private:
+
+ static G4SElectronPlus* theInstance;
+ G4SElectronPlus(){}
+ ~G4SElectronPlus(){}
+
+public:
+
+ static G4SElectronPlus* Definition();
+
+};
+
+#endif
diff --git a/Simulation/G4Extensions/Sleptons/Sleptons/G4SMuonMinus.hh b/Simulation/G4Extensions/Sleptons/Sleptons/G4SMuonMinus.hh
new file mode 100755
index 0000000000000000000000000000000000000000..e121f0491cc6d057e42e2639a1032ca5c745bc08
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/Sleptons/G4SMuonMinus.hh
@@ -0,0 +1,28 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef G4SMuonMinus_H
+#define G4SMuonMinus_H
+
+#include "globals.hh"
+#include "G4ios.hh"
+#include "G4ParticleDefinition.hh"
+
+class G4SMuonMinus: public G4ParticleDefinition
+{
+ // singleton implementation
+
+private:
+
+ static G4SMuonMinus* theInstance;
+ G4SMuonMinus(){}
+ ~G4SMuonMinus(){}
+
+public:
+
+ static G4SMuonMinus* Definition();
+
+};
+
+#endif
diff --git a/Simulation/G4Extensions/Sleptons/Sleptons/G4SMuonPlus.hh b/Simulation/G4Extensions/Sleptons/Sleptons/G4SMuonPlus.hh
new file mode 100755
index 0000000000000000000000000000000000000000..a0517a590f0687872057e63971ddd762dea6c616
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/Sleptons/G4SMuonPlus.hh
@@ -0,0 +1,28 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef G4SMuonPlus_H
+#define G4SMuonPlus_H
+
+#include "globals.hh"
+#include "G4ios.hh"
+#include "G4ParticleDefinition.hh"
+
+class G4SMuonPlus: public G4ParticleDefinition
+{
+ // singleton implementation
+
+private:
+
+ static G4SMuonPlus* theInstance;
+ G4SMuonPlus(){}
+ ~G4SMuonPlus(){}
+
+public:
+
+ static G4SMuonPlus* Definition();
+
+};
+
+#endif
diff --git a/Simulation/G4Extensions/Sleptons/Sleptons/G4STauMinus.hh b/Simulation/G4Extensions/Sleptons/Sleptons/G4STauMinus.hh
new file mode 100755
index 0000000000000000000000000000000000000000..82216d4dd7c2c819d261d70966cbf9c64803cbc5
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/Sleptons/G4STauMinus.hh
@@ -0,0 +1,28 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef G4STauMinus_H
+#define G4STauMinus_H
+
+#include "globals.hh"
+#include "G4ios.hh"
+#include "G4ParticleDefinition.hh"
+
+class G4STauMinus: public G4ParticleDefinition
+{
+ // singleton implementation
+
+private:
+
+ static G4STauMinus* theInstance;
+ G4STauMinus(){}
+ ~G4STauMinus(){}
+
+public:
+
+ static G4STauMinus* Definition();
+
+};
+
+#endif
diff --git a/Simulation/G4Extensions/Sleptons/Sleptons/G4STauPlus.hh b/Simulation/G4Extensions/Sleptons/Sleptons/G4STauPlus.hh
new file mode 100755
index 0000000000000000000000000000000000000000..ecbb8a0a92e5fd9ee78791ebbd1ccc8d93296349
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/Sleptons/G4STauPlus.hh
@@ -0,0 +1,28 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef G4STauPlus_H
+#define G4STauPlus_H
+
+#include "globals.hh"
+#include "G4ios.hh"
+#include "G4ParticleDefinition.hh"
+
+class G4STauPlus: public G4ParticleDefinition
+{
+ // singleton implementation
+
+private:
+
+ static G4STauPlus* theInstance;
+ G4STauPlus(){}
+ ~G4STauPlus(){}
+
+public:
+
+ static G4STauPlus* Definition();
+
+};
+
+#endif
diff --git a/Simulation/G4Extensions/Sleptons/Sleptons/SleptonsProcessDefinition.h b/Simulation/G4Extensions/Sleptons/Sleptons/SleptonsProcessDefinition.h
new file mode 100755
index 0000000000000000000000000000000000000000..c7d6063c2a03dd69d33bd5e4c247776a801ca456
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/Sleptons/SleptonsProcessDefinition.h
@@ -0,0 +1,27 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef SleptonsProcessDefinition_H
+#define SleptonsProcessDefinition_H
+
+#include "FadsPhysics/UserPhysicsDefinition.h"
+#include "G4ParticleDefinition.hh"
+
+class SleptonsProcessDefinition : public UserPhysicsDefinition {
+
+public:
+ SleptonsProcessDefinition();
+ void ConstructParticle();
+ void ConstructProcess();
+
+private:
+ G4ParticleDefinition *theSElectronMinus;
+ G4ParticleDefinition *theSElectronPlus;
+ G4ParticleDefinition *theSMuonMinus;
+ G4ParticleDefinition *theSMuonPlus;
+ G4ParticleDefinition *theSTauMinus;
+ G4ParticleDefinition *theSTauPlus;
+};
+
+#endif
diff --git a/Simulation/G4Extensions/Sleptons/cmt/requirements b/Simulation/G4Extensions/Sleptons/cmt/requirements
new file mode 100755
index 0000000000000000000000000000000000000000..87b377a3a694bd780c5edb9a6555d8e10ff9b8de
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/cmt/requirements
@@ -0,0 +1,15 @@
+package Sleptons
+
+author ADA <dellacqu@mail.cern.ch>
+
+public
+use AtlasPolicy AtlasPolicy-*
+use FadsPhysics FadsPhysics-* Simulation/G4Sim/FADS
+use Geant4 Geant4-* External
+
+include_dirs "$(Sleptons_root)"
+
+library Sleptons ../src/*.cxx
+
+apply_pattern linked_library
+
diff --git a/Simulation/G4Extensions/Sleptons/doc/mainpage.h b/Simulation/G4Extensions/Sleptons/doc/mainpage.h
new file mode 100644
index 0000000000000000000000000000000000000000..463a21e32187f8fa2f67901953580caf5bddeec5
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/doc/mainpage.h
@@ -0,0 +1,34 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+/**
+@mainpage Sleptons
+@author Andrea Dell'Acqua (dellacqu@mail.cern.ch)
+
+@section Sleptons Introduction
+
+This package was created in order to introduce sleptons (supersymmetric partners of leptons which could be stable) into the Geant4 simulation. The properties of these sleptons are not set by this package, since they may vary considerably depending on the supersymmetric theory being considered. They should be set up in user job options.
+
+@section Sleptons Class Overview
+
+The package contains only one class:
+
+ - SleptonsProcessDefinition : This class adds the appropriate physics processes to the slepton (for example transportation).
+
+@ref used_Sleptons
+
+@ref requirements_Sleptons
+
+*/
+
+/**
+@page used_Sleptons Used Packages
+@htmlinclude used_packages.html
+*/
+
+/**
+@page requirements_Sleptons Requirements
+@include requirements
+*/
+
diff --git a/Simulation/G4Extensions/Sleptons/src/G4SElectronMinus.cxx b/Simulation/G4Extensions/Sleptons/src/G4SElectronMinus.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..0ac73816d6c539d1586e5ae86f2cd62cabf4c766
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/src/G4SElectronMinus.cxx
@@ -0,0 +1,35 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+
+#include "Sleptons/G4SElectronMinus.hh"
+// ######################################################################
+// ### SElectronMinus ###
+// ######################################################################
+
+G4SElectronMinus* G4SElectronMinus::theInstance = Definition();
+
+G4SElectronMinus* G4SElectronMinus::Definition()
+{
+ if (theInstance !=0) return theInstance;
+
+// Arguments for constructor are as follows
+// name mass width charge
+// 2*spin parity C-conjugation
+// 2*Isospin 2*Isospin3 G-parity
+// type lepton number baryon number PDG encoding
+// stable lifetime decay table
+// shortlived subType anti_encoding
+
+ G4ParticleDefinition* anInstance = new G4ParticleDefinition(
+ "s_e_minus_R", 100.00*CLHEP::GeV, 0.0*CLHEP::MeV, -1.*CLHEP::eplus,
+ 0, 0, 0,
+ 0, 0, 0,
+ "slepton", 1, 0, 2000011,
+ true, -1.0, NULL,
+ false,"SElectronMinus");
+
+ theInstance = reinterpret_cast<G4SElectronMinus*>(anInstance);
+ return theInstance;
+}
diff --git a/Simulation/G4Extensions/Sleptons/src/G4SElectronPlus.cxx b/Simulation/G4Extensions/Sleptons/src/G4SElectronPlus.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..91a634e3d32d9f268f743f06b94daaba17139d1a
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/src/G4SElectronPlus.cxx
@@ -0,0 +1,35 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+
+#include "Sleptons/G4SElectronPlus.hh"
+// ######################################################################
+// ### SElectronPlus ###
+// ######################################################################
+
+G4SElectronPlus* G4SElectronPlus::theInstance = Definition();
+
+G4SElectronPlus* G4SElectronPlus::Definition()
+{
+ if (theInstance !=0) return theInstance;
+
+// Arguments for constructor are as follows
+// name mass width charge
+// 2*spin parity C-conjugation
+// 2*Isospin 2*Isospin3 G-parity
+// type lepton number baryon number PDG encoding
+// stable lifetime decay table
+// shortlived subType anti_encoding
+
+ G4ParticleDefinition* anInstance = new G4ParticleDefinition(
+ "s_e_plus_R", 100.00*CLHEP::GeV, 0.0*CLHEP::MeV, +1.*CLHEP::eplus,
+ 0, 0, 0,
+ 0, 0, 0,
+ "slepton", 1, 0, -2000011,
+ true, -1.0, NULL,
+ false, "SElectronPlus");
+
+ theInstance = reinterpret_cast<G4SElectronPlus*>(anInstance);
+ return theInstance;
+}
diff --git a/Simulation/G4Extensions/Sleptons/src/G4SMuonMinus.cxx b/Simulation/G4Extensions/Sleptons/src/G4SMuonMinus.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..c565c259d8e480d5a81eb40682266cdb0fc7da96
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/src/G4SMuonMinus.cxx
@@ -0,0 +1,35 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+
+#include "Sleptons/G4SMuonMinus.hh"
+// ######################################################################
+// ### SMuonMinus ###
+// ######################################################################
+
+G4SMuonMinus* G4SMuonMinus::theInstance = Definition();
+
+G4SMuonMinus* G4SMuonMinus::Definition()
+{
+ if (theInstance !=0) return theInstance;
+
+// Arguments for constructor are as follows
+// name mass width charge
+// 2*spin parity C-conjugation
+// 2*Isospin 2*Isospin3 G-parity
+// type lepton number baryon number PDG encoding
+// stable lifetime decay table
+// shortlived subType anti_encoding
+
+ G4ParticleDefinition* anInstance = new G4ParticleDefinition(
+ "s_mu_minus_R", 100.00*CLHEP::GeV, 0.0*CLHEP::MeV, -1.*CLHEP::eplus,
+ 0, 0, 0,
+ 0, 0, 0,
+ "slepton", 1, 0, 2000013,
+ true, -1.0, NULL,
+ false, "SMuonMinus");
+
+ theInstance = reinterpret_cast<G4SMuonMinus*>(anInstance);
+ return theInstance;
+}
diff --git a/Simulation/G4Extensions/Sleptons/src/G4SMuonPlus.cxx b/Simulation/G4Extensions/Sleptons/src/G4SMuonPlus.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..861dddfeaccb00e5ad1586e6b6079665edff9827
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/src/G4SMuonPlus.cxx
@@ -0,0 +1,35 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+
+#include "Sleptons/G4SMuonPlus.hh"
+// ######################################################################
+// ### SMuonPlus ###
+// ######################################################################
+
+G4SMuonPlus* G4SMuonPlus::theInstance = Definition();
+
+G4SMuonPlus* G4SMuonPlus::Definition()
+{
+ if (theInstance !=0) return theInstance;
+
+// Arguments for constructor are as follows
+// name mass width charge
+// 2*spin parity C-conjugation
+// 2*Isospin 2*Isospin3 G-parity
+// type lepton number baryon number PDG encoding
+// stable lifetime decay table
+// shortlived subType anti_encoding
+
+ G4ParticleDefinition* anInstance = new G4ParticleDefinition(
+ "s_mu_plus_R", 100.00*CLHEP::GeV, 0.0*CLHEP::MeV, +1.*CLHEP::eplus,
+ 0, 0, 0,
+ 0, 0, 0,
+ "slepton", 1, 0, -2000013,
+ true, -1.0, NULL,
+ false, "SMuonPlus");
+
+ theInstance = reinterpret_cast<G4SMuonPlus*>(anInstance);
+ return theInstance;
+}
diff --git a/Simulation/G4Extensions/Sleptons/src/G4STauMinus.cxx b/Simulation/G4Extensions/Sleptons/src/G4STauMinus.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..e2869dbf47c6e1b9690a020334a512084299ee86
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/src/G4STauMinus.cxx
@@ -0,0 +1,35 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+
+#include "Sleptons/G4STauMinus.hh"
+// ######################################################################
+// ### STauMinus ###
+// ######################################################################
+
+G4STauMinus* G4STauMinus::theInstance = Definition();
+
+G4STauMinus* G4STauMinus::Definition()
+{
+ if (theInstance !=0) return theInstance;
+
+// Arguments for constructor are as follows
+// name mass width charge
+// 2*spin parity C-conjugation
+// 2*Isospin 2*Isospin3 G-parity
+// type lepton number baryon number PDG encoding
+// stable lifetime decay table
+// shortlived subType anti_encoding
+
+ G4ParticleDefinition* anInstance = new G4ParticleDefinition(
+ "s_tau_minus_1", 100.*CLHEP::GeV, 0.0*CLHEP::MeV, -1.*CLHEP::eplus,
+ 0, 0, 0,
+ 0, 0, 0,
+ "slepton", 1, 0, 1000015,
+ true, -1.0, NULL,
+ false, "STauMinus");
+
+ theInstance = reinterpret_cast<G4STauMinus*>(anInstance);
+ return theInstance;
+}
diff --git a/Simulation/G4Extensions/Sleptons/src/G4STauPlus.cxx b/Simulation/G4Extensions/Sleptons/src/G4STauPlus.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..e15db0b51821de86be1c3c9793e4b39815d6b4f1
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/src/G4STauPlus.cxx
@@ -0,0 +1,35 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+
+#include "Sleptons/G4STauPlus.hh"
+// ######################################################################
+// ### STauPlus ###
+// ######################################################################
+
+G4STauPlus* G4STauPlus::theInstance = Definition();
+
+G4STauPlus* G4STauPlus::Definition()
+{
+ if (theInstance !=0) return theInstance;
+
+// Arguments for constructor are as follows
+// name mass width charge
+// 2*spin parity C-conjugation
+// 2*Isospin 2*Isospin3 G-parity
+// type lepton number baryon number PDG encoding
+// stable lifetime decay table
+// shortlived subType anti_encoding
+
+ G4ParticleDefinition* anInstance = new G4ParticleDefinition(
+ "s_tau_plus_1", 100.*CLHEP::GeV, 0.0*CLHEP::MeV, +1.*CLHEP::eplus,
+ 0, 0, 0,
+ 0, 0, 0,
+ "slepton", 1, 0, -1000015,
+ true, -1.0, NULL,
+ false, "STauPlus");
+
+ theInstance = reinterpret_cast<G4STauPlus*>(anInstance);
+ return theInstance;
+}
diff --git a/Simulation/G4Extensions/Sleptons/src/SleptonsProcessDefinition.cxx b/Simulation/G4Extensions/Sleptons/src/SleptonsProcessDefinition.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..7341b8259d34552b4c19191aa9fad06a414b1a49
--- /dev/null
+++ b/Simulation/G4Extensions/Sleptons/src/SleptonsProcessDefinition.cxx
@@ -0,0 +1,62 @@
+/*
+ Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
+*/
+
+#include "Sleptons/SleptonsProcessDefinition.h"
+
+#include "Sleptons/G4SElectronMinus.hh"
+#include "Sleptons/G4SElectronPlus.hh"
+#include "Sleptons/G4SMuonMinus.hh"
+#include "Sleptons/G4SMuonPlus.hh"
+#include "Sleptons/G4STauMinus.hh"
+#include "Sleptons/G4STauPlus.hh"
+
+#include "FadsPhysics/UserPhysicsDefinitionProxy.h"
+
+#include "G4hIonisation.hh"
+#include "G4hMultipleScattering.hh"
+
+#include "G4Transportation.hh"
+#include "G4MuIonisation.hh"
+#include "G4ProcessManager.hh"
+
+static UserPhysicsDefinitionProxy<SleptonsProcessDefinition> plugin;
+
+SleptonsProcessDefinition::SleptonsProcessDefinition()
+ : UserPhysicsDefinition(),
+ theSElectronMinus(0), theSElectronPlus(0), theSMuonMinus(0), theSMuonPlus(0), theSTauMinus(0), theSTauPlus(0)
+{
+ std::cout <<" physics list for the Sleptons being built "<<std::endl;
+}
+
+void SleptonsProcessDefinition::ConstructProcess()
+{
+ std::cout <<" ConstructProcess for Sleptons being run"<<std::endl;
+ theParticleIterator->reset();
+ while( (*theParticleIterator)() ){
+ G4ParticleDefinition* particle = theParticleIterator->value();
+ int pdgcode=std::abs(particle->GetPDGEncoding());
+ if (pdgcode==2000011 || pdgcode==2000013 || pdgcode==1000015)
+ {
+ G4ProcessManager* proc=particle->GetProcessManager();
+ proc->AddProcess(new G4hMultipleScattering,-1,1,1);
+ proc->AddProcess(new G4hIonisation,-1,2,2);
+ }
+ }
+
+
+}
+
+void SleptonsProcessDefinition::ConstructParticle()
+{
+ std::cout << " ConstructParticle for the Sleptons being run" << std::endl;
+
+ theSElectronMinus = G4SElectronMinus::Definition();
+ theSElectronPlus = G4SElectronPlus::Definition();
+
+ theSMuonMinus = G4SMuonMinus::Definition();
+ theSMuonPlus = G4SMuonPlus::Definition();
+
+ theSTauMinus=G4STauMinus::Definition();
+ theSTauPlus=G4STauPlus::Definition();
+}