FEDRA emulsion software from the OPERA Collaboration
Public Member Functions | Private Member Functions | Private Attributes | List of all members
EdbShowAlgE Class Reference

#include <EdbShowAlgE.h>

Inheritance diagram for EdbShowAlgE:
Inheritance graph
[legend]
Collaboration diagram for EdbShowAlgE:
Collaboration graph
[legend]

Public Member Functions

void AddWeightFile (TString tmpString, Int_t Plate, TString tmpString2)
 
Int_t CheckPlateAndPlateBinning (Int_t Plate)
 
 ClassDef (EdbShowAlgE, 1)
 
void Create_MLP ()
 DEBUG ---------------------------— TEST DEVOPLEMENT. More...
 
void Create_MLPTree ()
 
void CreateInputNeuronsArray ()
 
void CreateInputNeuronsArray_XX ()
 
void CreateInputNeuronsArray_YC ()
 
void CreateInputNeuronsArray_YY ()
 
void CreateMLPArray ()
 
void CreateMLPLayout_XX (Int_t nr)
 
void CreateMLPLayout_YC (Int_t nr)
 
void CreateMLPLayout_YY (Int_t nr)
 
void CreateMLPLayouts ()
 
void CreateMLPLayouts_FJ ()
 
void CreateMLPLayouts_JC ()
 
void CreateMLPLayouts_LT ()
 
void CreateMLPLayouts_XX ()
 
void CreateMLPLayouts_YC ()
 
void CreateMLPLayouts_YY ()
 
void CreateMLPs ()
 
void CreateStandardMLPArrayN ()
 
void CreateStandardWeightFileStrings ()
 
void DumpWeights ()
 
 EdbShowAlgE ()
 
 EdbShowAlgE (TObjArray *RecoShowerArray)
 
 EdbShowAlgE (TObjArray *RecoShowerArray, TString ShowAlgEParaName)
 
 EdbShowAlgE (TString ShowAlgEParaName)
 
virtual void Execute ()
 
void Fill_MLPTree ()
 
void FillInputNeuronsArray_XX ()
 
void FillInputNeuronsArray_YC ()
 
void FillInputNeuronsArray_YY ()
 
void FillMLP_SimpleTree ()
 
void FillMLP_SimpleTree_XX ()
 
void FillMLP_SimpleTree_XX (EdbShowerP *show, Int_t ShowAlgEParaNr)
 
void FillMLP_SimpleTree_YC ()
 
void FillMLP_SimpleTree_YC (EdbShowerP *show, Int_t ShowAlgEParaNr)
 
void FillMLP_SimpleTree_YY ()
 
void FillMLP_SimpleTree_YY (EdbShowerP *show, Int_t ShowAlgEParaNr)
 
virtual void Finalize ()
 
TH1F * GetHistInANN (Int_t nr) const
 
TCanvas * GetHistoCanvas ()
 
TMultiLayerPerceptron * GetMLP () const
 
int GetMLPArrayNr (Int_t NPlToCheck)
 
TString GetMLPInputFileName ()
 
TString GetMLPLayout ()
 
TString GetMLPWeightFileName ()
 
virtual void Initialize ()
 
void LoadMLPWeights ()
 
void LoadMLPWeights (Int_t nr)
 
void Print ()
 
void PrintWeightFileStrings ()
 
void ReadMLPInputFile ()
 
int ReadShowAlgEPar (const char *file)
 DEBUG ---------------------------— TEST DEVOPLEMENT. More...
 
void ResetPlateBinning ()
 
void Run ()
 
void SetBinsInputHistos ()
 
void SetBinsInputHistos_FJ ()
 
void SetBinsInputHistos_XX ()
 
void SetBinsInputHistos_YC ()
 
void SetBinsInputHistos_YY ()
 
void SetDoRun (Int_t type)
 
void SetDoTrain (Int_t type)
 
void SetInputHistos (TH1F *h, Int_t nbin, Double_t lower, Double_t upper)
 
void SetInputHistosPara (Int_t ParaType)
 
void SetInputHistosPara_FJ ()
 
void SetInputHistosPara_XX ()
 
void SetInputHistosPara_YC ()
 
void SetInputHistosPara_YY ()
 
void SetMLPArrayN (Int_t MLPArrayN)
 
void SetMLPInputFileName (TString MLP_InputFileName)
 
void SetMLPLayout (TString layout)
 
Int_t SetMLPTrainEpochs ()
 
void SetMLPTrainEpochs (Int_t nep)
 
void SetMLPWeightFileName (TString weightfilename)
 
void SetPlatesN (Int_t ShowAlgEPlatesN)
 
void SetReadMLPInputFile (Int_t type)
 
void SetRecoShowerArray (TObjArray *RecoShowerArray)
 
void SetWriteMLPInputFile (Int_t type)
 
void SetWriteMLPWeightFile (Int_t type)
 
void TESTRUN ()
 
void Train ()
 
void WriteMLP_SimpleTree (TString filename)
 
void WriteMLPInputFile ()
 
virtual ~EdbShowAlgE ()
 

Private Member Functions

void AddWeightFilePlate (TString tmpString, Int_t dospec)
 
Bool_t CheckStringToShowAlgEParaName (TString tmpString2)
 
void ClearPlateBinning ()
 
void CreateInputHistos ()
 
TString CreateMLPLayout (Int_t NrOfANNInputNeurons)
 
void FillInputHistosPara_XX ()
 
int GetLastBinHistoFilles (TH1 *h)
 
void Init ()
 
void Init (TString ParaName)
 
void ResetWeightFileNameArray ()
 
void Set0 ()
 

Private Attributes

Bool_t eDoANNRun
 
Bool_t eDoANNTrain
 
Float_t * eInANN = new Float_t[70]
 
TH1F * eInANN_Hist [70]
 
Float_t * eInfoANN = new Float_t[70]
 
TString eMLP_InputFileName
 
Int_t eMLP_InputNeurons
 
TArrayI * eMLP_InputNeuronsArray
 
Int_t eMLP_InputNeuronsArrayN
 
TString eMLP_Layout
 
TObjArray * eMLP_LayoutArray
 
Bool_t eMLP_LayoutArrayIsDone
 
Bool_t eMLP_ReadMLPInputFile
 
TMultiLayerPerceptron * eMLP_Simple
 
Bool_t eMLP_SimpleIsDone
 
TTree * eMLP_SimpleTree
 
Int_t eMLP_TrainEpochs
 
TString eMLP_WeightFileName
 
TString eMLP_WeightFileNameArray [10]
 
Bool_t eMLP_WriteMLPInputFile
 
Bool_t eMLP_WriteWeightFile
 
TObjArray * eMLPArray
 
Int_t eMLPArrayN
 
Float_t eOutANN
 
TObjArray * eRecoShowerArray
 
Int_t eRecoShowerArrayN
 
TString eShowAlgEParaName
 
Int_t eShowAlgEParaNr
 
TArrayI * eShowAlgEPlateBinning
 
Int_t eShowAlgEPlatesN
 

Constructor & Destructor Documentation

◆ EdbShowAlgE() [1/4]

EdbShowAlgE::EdbShowAlgE ( )

◆ EdbShowAlgE() [2/4]

EdbShowAlgE::EdbShowAlgE ( TString  ShowAlgEParaName = "XX")
31{
32 // Default Constructor
33 cout << "EdbShowAlgE::EdbShowAlgE() Default Constructor with Paraname ="<< paraname<< endl;
34
35 // Reset all:
36 Set0();
37
38 // Init now with given parametrization
39 Init(paraname);
40}
EdbShowAlgE::Set0
void Set0()
Definition: EdbShowAlgE.cxx:100
EdbShowAlgE::Init
void Init()
Definition: EdbShowAlgE.cxx:163

◆ EdbShowAlgE() [3/4]

EdbShowAlgE::EdbShowAlgE ( TObjArray *  RecoShowerArray)
49{
50 // Default Constructor
51 cout << "EdbShowAlgE::EdbShowAlgE(TObjArray* RecoShowerArray) Default Constructor"<<endl;
52
53 // Reset all:
54 Set0();
55
56 // Init now with default parametrization (XX)
57 Init();
58
59 eRecoShowerArray=RecoShowerArray;
60 eRecoShowerArrayN=eRecoShowerArray->GetEntries();
61
62 cout << "EdbShowAlgE::EdbShowAlgE(TObjArray* RecoShowerArray) Default Constructor...done."<<endl;
63}
RecoShowerArray
TObjArray * RecoShowerArray
Definition: Shower_E_FromShowerRoot.C:12
EdbShowAlgE::eRecoShowerArray
TObjArray * eRecoShowerArray
Definition: EdbShowAlgE.h:58
EdbShowAlgE::eRecoShowerArrayN
Int_t eRecoShowerArrayN
Definition: EdbShowAlgE.h:59

◆ EdbShowAlgE() [4/4]

EdbShowAlgE::EdbShowAlgE ( TObjArray *  RecoShowerArray,
TString  ShowAlgEParaName = "XX" 
)
71{
72 // Default Constructor
73 cout << "EdbShowAlgE::EdbShowAlgE(TObjArray* RecoShowerArray,TString paraname) Default Constructor with Paraname ="<< paraname<< endl;
74
75 // Reset all:
76 Set0();
77
78 // Init now with given parametrization
79 Init(paraname);
80
81 eRecoShowerArray=RecoShowerArray;
82 eRecoShowerArrayN=eRecoShowerArray->GetEntries();
83
84 cout << "EdbShowAlgE::EdbShowAlgE(TObjArray* RecoShowerArray) Default Constructor...done."<<endl;
85}

◆ ~EdbShowAlgE()

EdbShowAlgE::~EdbShowAlgE ( )
virtual
92{
93 // Default Destructor
94 cout << "EdbShowAlgE::~EdbShowAlgE()"<<endl;
95}

Member Function Documentation

◆ AddWeightFile()

void EdbShowAlgE::AddWeightFile ( TString  tmpString,
Int_t  Plate,
TString  tmpString2 
)
1499{
1500 cout << "void EdbShowAlgE::AddWeightFile()." << endl;
1501
1502 // Check if the weightfile type and the ShowAlgeE paratye match!
1503 // (otherwise a weightfile for a wrong parametrisation is maybe taken)
1504 if (tmpString!=eShowAlgEParaName) {
1505 cout << tmpString << " != "<< eShowAlgEParaName << endl;
1506 cout << "void EdbShowAlgE::AddWeightFilePlate() WARNING: WRONG PARANAME. WeightFilePlate NOT ADDED!" << endl;
1507 return;
1508 }
1509
1510 // Things to do when called the first time....
1511 if (eMLPArrayN<1) ClearPlateBinning();
1512
1513
1514
1515
1516 // Not more than 10 values of weightfiles will be allowed!!!
1517 if (eMLPArrayN>=10) {
1518 cout << "void EdbShowAlgE::AddWeightFile() WARNING: Not more than 10 values of weightfiles will be allowed!!. WeightFilePlate NOT ADDED!" << endl;
1519 return;
1520 }
1521
1522 // Check if this weightfile to be added (for the plates) is not added already:
1523 Int_t tmpint;
1524 tmpint=CheckPlateAndPlateBinning(Plate);
1525
1526 if (tmpint==1) {
1527 cout << "void EdbShowAlgE::AddWeightFile() WARNING: A weightfile for this plate Number was alread added. NOT ADDED!" << endl;
1528 return;
1529 }
1530
1531
1532 // Set weightfilename string and increase number of weightfile options
1533 eShowAlgEPlateBinning->AddAt(Plate,eMLPArrayN);
1534 tmpint=eShowAlgEPlateBinning->At(eMLPArrayN);
1535 eMLP_WeightFileNameArray[eMLPArrayN]=tmpString2;
1536
1537 cout << "eMLP_WeightFileNameArray[eMLPArrayN] = " << eMLP_WeightFileNameArray[eMLPArrayN] << endl;
1538
1539
1540
1541 // Increase now...
1542 ++eMLPArrayN;
1543
1544 cout << "------"<< endl;
1545 cout << "eMLPArrayN = " << eMLPArrayN << endl;
1546
1547 // These two variables are always same!
1548 eMLP_InputNeuronsArrayN=eMLPArrayN;
1549
1550
1551 cout << "void EdbShowAlgE::AddWeightFile() done." << endl;
1552
1553 return;
1554}
EdbShowAlgE::eShowAlgEParaName
TString eShowAlgEParaName
Definition: EdbShowAlgE.h:47
EdbShowAlgE::eMLP_InputNeuronsArrayN
Int_t eMLP_InputNeuronsArrayN
Definition: EdbShowAlgE.h:115
EdbShowAlgE::eMLPArrayN
Int_t eMLPArrayN
Definition: EdbShowAlgE.h:64
EdbShowAlgE::CheckPlateAndPlateBinning
Int_t CheckPlateAndPlateBinning(Int_t Plate)
Definition: EdbShowAlgE.cxx:1581
EdbShowAlgE::eMLP_WeightFileNameArray
TString eMLP_WeightFileNameArray[10]
Definition: EdbShowAlgE.h:109
EdbShowAlgE::ClearPlateBinning
void ClearPlateBinning()
Definition: EdbShowAlgE.cxx:1338
EdbShowAlgE::eShowAlgEPlateBinning
TArrayI * eShowAlgEPlateBinning
Definition: EdbShowAlgE.h:53

◆ AddWeightFilePlate()

void EdbShowAlgE::AddWeightFilePlate ( TString  tmpString,
Int_t  dospec 
)
private
1438{
1439 cout << "void EdbShowAlgE::AddWeightFilePlate() " << endl;
1440
1441 cout << "eMLPArrayN = " << eMLPArrayN << endl;
1442 cout << "tmpString = " << tmpString << endl;
1443 cout << "Plate = " << Plate << endl;
1444
1445
1446 // Check if the weightfile type and the ShowAlgeE paratye match!
1447 // (otherwise a weightfile for a wrong parametrisation is maybe taken)
1448 if (tmpString!=eShowAlgEParaName) {
1449 cout << tmpString << " != "<< eShowAlgEParaName << endl;
1450 cout << "void EdbShowAlgE::AddWeightFilePlate() WARNING: WRONG PARANAME. WeightFilePlate NOT ADDED!" << endl;
1451 return;
1452 }
1453
1454
1455 // Things to do when called the first time....
1456 if (eMLPArrayN<1) ClearPlateBinning();
1457
1458 // Not more than 10 values of weightfiles will be allowed!!!
1459 if (eMLPArrayN>=10) {
1460 cout << "void EdbShowAlgE::AddWeightFilePlate() WARNING: Not more than 10 values of weightfiles will be allowed!!. WeightFilePlate NOT ADDED!" << endl;
1461 return;
1462 }
1463
1464
1465 // Check if this weightfile to be added (for the plates) is not added already:
1466 Int_t tmpint;
1467 tmpint=CheckPlateAndPlateBinning(Plate);
1468
1469 if (tmpint==1) {
1470 cout << "void EdbShowAlgE::AddWeightFilePlate() WARNING: A weightfile for this plate Number was alread added. NOT ADDED!" << endl;
1471 return;
1472 }
1473
1474
1475 // Set weightfilename string and increase number of weightfile options
1476 eShowAlgEPlateBinning->AddAt(Plate,eMLPArrayN);
1477 tmpint=eShowAlgEPlateBinning->At(eMLPArrayN);
1478 eMLP_WeightFileNameArray[eMLPArrayN]="ANN_Weights_Energy_"+tmpString+TString(Form("_%d.txt",tmpint));
1479
1480 cout << "eMLP_WeightFileNameArray[eMLPArrayN] = " << eMLP_WeightFileNameArray[eMLPArrayN] << endl;
1481
1482 // Increase now...
1483 ++eMLPArrayN;
1484
1485 cout << "------"<< endl;
1486 cout << "eMLPArrayN = " << eMLPArrayN << endl;
1487
1488 // These two variables are always same!
1489 eMLP_InputNeuronsArrayN=eMLPArrayN;
1490 cout << "void EdbShowAlgE::AddWeightFilePlate() done." << endl;
1491 return;
1492}

◆ CheckPlateAndPlateBinning()

Int_t EdbShowAlgE::CheckPlateAndPlateBinning ( Int_t  Plate)
1582{
1583 if (eMLPArrayN<1) return 0;
1584
1585 Int_t tmpint;
1586 for (Int_t i=0; i<eMLPArrayN; ++i) {
1587 tmpint=eShowAlgEPlateBinning->At(i);
1588 //cout <<"CheckPlateAndPlateBinning , Plate, i =eShowAlgEPlateBinning->At(i) " << Plate << " " << i << " " << tmpint<< endl;
1589 if (Plate==tmpint) return 1;
1590 }
1591 return 0;
1592}

◆ CheckStringToShowAlgEParaName()

Bool_t EdbShowAlgE::CheckStringToShowAlgEParaName ( TString  tmpString2)
private
2481{
2482 if ( tmpString2 == eShowAlgEParaName ) return kTRUE;
2483 cout <<"Bool_t EdbShowAlgE::CheckStringToShowAlgEParaName() " <<tmpString2 << " does not matcht to " << eShowAlgEParaName;
2484 cout << " Return False."<< endl;
2485 return kFALSE;
2486}

◆ ClassDef()

EdbShowAlgE::ClassDef ( EdbShowAlgE  ,
1   
)

◆ ClearPlateBinning()

void EdbShowAlgE::ClearPlateBinning ( )
private
1339{
1340 cout << "void EdbShowAlgE::ClearPlateBinning() " << endl;
1341 // Restore Standard Plate Binning
1342 if (eShowAlgEPlateBinning) {
1343 delete eShowAlgEPlateBinning;
1344 eShowAlgEPlateBinning=0;
1345
1346 }
1347 else {
1348 // do nothing;
1349 }
1350
1351 // Not more than 10 values of weightfiles will be allowed!!!
1352 eShowAlgEPlateBinning = new TArrayI(10);
1353 for (Int_t i=0; i<10; ++i) eShowAlgEPlateBinning->AddAt(0,i);
1354
1355 cout << "void EdbShowAlgE::ClearPlateBinning()...done. " << endl;
1356 return;
1357}

◆ Create_MLP()

void EdbShowAlgE::Create_MLP ( )

DEBUG ---------------------------— TEST DEVOPLEMENT.

775{
776 cout << "EdbShowAlgE::Create_MLP()" << endl;
777
778 TString testlayout="HALLO";
779 // This class implementation of TMultiLayerPerceptron has strange behaviors...
780 // The usual eMLP_Simple=0 resetings dont take effect and in the header file this has to be declared with
781 // Brackets.....I do not know why and try to keep it working in the standard way as much as possible!
782 TMultiLayerPerceptron* eMLP_Simple = new TMultiLayerPerceptron(testlayout);
783
784 eMLP_Simple->Print();
785
786 eMLP_SimpleIsDone=kTRUE;
787
788 cout << "EdbShowAlgE::Create_MLP()...done." << endl;
789 return;
790}
EdbShowAlgE::eMLP_SimpleIsDone
Bool_t eMLP_SimpleIsDone
Definition: EdbShowAlgE.h:99
EdbShowAlgE::eMLP_Simple
TMultiLayerPerceptron * eMLP_Simple
Definition: EdbShowAlgE.h:97

◆ Create_MLPTree()

void EdbShowAlgE::Create_MLPTree ( )
1048{
1049
1050 Log(2,"EdbShowAlgE::Create_MLPTree","EdbShowAlgE::Create_MLPTree.");
1051
1052 if (eMLP_SimpleTree) {
1053 delete eMLP_SimpleTree;
1054 eMLP_SimpleTree=0;
1055 cout << "deletion done"<<endl;
1056 }
1057
1058 if (!eMLP_SimpleTree) {
1059 eMLP_SimpleTree = new TTree("eMLP_SimpleTree","eMLP_SimpleTree");
1060 eMLP_SimpleTree -> Branch("eInANN", eInANN, "eInANN[70]/F");
1061 eMLP_SimpleTree -> Branch("eInfoANN", eInfoANN, "eInfoANN[70]/F");
1062 }
1063
1064 Log(2,"EdbShowAlgE::Create_MLPTree","EdbShowAlgE::Create_MLPTree....done");
1065 return;
1066}
Branch
ParaSet Branch("CUT_BACK_DMIN",&cut_back_dmin,"CUT_BACK_DMIN/D")
Log
bool Log(int level, const char *location, const char *fmt,...)
Definition: EdbLog.cxx:75
EdbShowAlgE::eInANN
Float_t * eInANN
Definition: EdbShowAlgE.h:85
EdbShowAlgE::eInfoANN
Float_t * eInfoANN
Definition: EdbShowAlgE.h:89
EdbShowAlgE::eMLP_SimpleTree
TTree * eMLP_SimpleTree
Definition: EdbShowAlgE.h:100

◆ CreateInputHistos()

void EdbShowAlgE::CreateInputHistos ( )
private
1146{
1147 // Create Input Histos:
1148 // Bin Width and Bounds have to be overwritten later when parametrization is specified
1149 for (Int_t j=0; j<70; j++) {
1150 char hname[40];
1151 sprintf(hname,"Hist_eInANN_Hist_%d",j);
1152 if (!eInANN_Hist[j]) eInANN_Hist[j] = new TH1F(hname,hname,100,0, 1000.0);
1153 }
1154 return;
1155}
EdbShowAlgE::eInANN_Hist
TH1F * eInANN_Hist[70]
Definition: EdbShowAlgE.h:92

◆ CreateInputNeuronsArray()

void EdbShowAlgE::CreateInputNeuronsArray ( )
1916{
1917 cout << "void EdbShowAlgE::CreateInputNeuronsArray() " << endl;
1918
1919 switch (eShowAlgEParaNr) {
1920 case 0 : //CreateInputNeuronsArray_FJ();
1921 break;
1922 case 1 : //CreateInputNeuronsArray_LT();
1923 break;
1924 case 2 :
1925 CreateInputNeuronsArray_YC();
1926 break;
1927 case 3 : //CreateInputNeuronsArray_JC();
1928 break;
1929 case 4 :
1930 CreateInputNeuronsArray_XX();
1931 break;
1932
1933 case 5 :
1934 CreateInputNeuronsArray_YY();
1935 break;
1936
1937 default :
1938 CreateInputNeuronsArray_XX();
1939 break;
1940 }
1941
1942 cout << "void EdbShowAlgE::CreateInputNeuronsArray() done..." << endl;
1943 return;
1944}
EdbShowAlgE::CreateInputNeuronsArray_YC
void CreateInputNeuronsArray_YC()
Definition: EdbShowAlgE.cxx:1951
EdbShowAlgE::CreateInputNeuronsArray_YY
void CreateInputNeuronsArray_YY()
Definition: EdbShowAlgE.cxx:2044
EdbShowAlgE::eShowAlgEParaNr
Int_t eShowAlgEParaNr
Definition: EdbShowAlgE.h:48
EdbShowAlgE::CreateInputNeuronsArray_XX
void CreateInputNeuronsArray_XX()
Definition: EdbShowAlgE.cxx:1996

◆ CreateInputNeuronsArray_XX()

void EdbShowAlgE::CreateInputNeuronsArray_XX ( )
1997{
1998 cout << "void EdbShowAlgE::CreateInputNeuronsArray_XX() " << endl;
1999 // Restore Standard Plate Binning
2000 if (eMLP_InputNeuronsArray) {
2001 delete eMLP_InputNeuronsArray;
2002 eMLP_InputNeuronsArray=0;
2003 }
2004 else {
2005 // do nothing;
2006 }
2007
2008 // Not more than 10 values of weightfiles will be allowed!
2009 eMLP_InputNeuronsArray = new TArrayI(10);
2010 for (Int_t i=0; i<10; ++i) eMLP_InputNeuronsArray->AddAt(0,i);
2011
2012 cout << "void EdbShowAlgE::CreateInputNeuronsArray_XX()...done. " << endl;
2013 return;
2014}
EdbShowAlgE::eMLP_InputNeuronsArray
TArrayI * eMLP_InputNeuronsArray
Definition: EdbShowAlgE.h:114

◆ CreateInputNeuronsArray_YC()

void EdbShowAlgE::CreateInputNeuronsArray_YC ( )
1952{
1953 cout << "void EdbShowAlgE::CreateInputNeuronsArray_YC() " << endl;
1954 // Restore Standard Plate Binning
1955 if (eMLP_InputNeuronsArray) {
1956 delete eMLP_InputNeuronsArray;
1957 eMLP_InputNeuronsArray=0;
1958 }
1959 else {
1960 // do nothing;
1961 }
1962
1963 // Not more than 10 values of weightfiles will be allowed!
1964 eMLP_InputNeuronsArray = new TArrayI(10);
1965 for (Int_t i=0; i<10; ++i) eMLP_InputNeuronsArray->AddAt(0,i);
1966
1967 cout << "void EdbShowAlgE::CreateInputNeuronsArray_YC()...done. " << endl;
1968 return;
1969}

◆ CreateInputNeuronsArray_YY()

void EdbShowAlgE::CreateInputNeuronsArray_YY ( )
2045{
2046 cout << "void EdbShowAlgE::CreateInputNeuronsArray_YY() " << endl;
2047 // Restore Standard Plate Binning
2048 if (eMLP_InputNeuronsArray) {
2049 delete eMLP_InputNeuronsArray;
2050 eMLP_InputNeuronsArray=0;
2051 }
2052 else {
2053 // do nothing;
2054 }
2055
2056 // Not more than 10 values of weightfiles will be allowed!
2057 eMLP_InputNeuronsArray = new TArrayI(10);
2058 for (Int_t i=0; i<10; ++i) eMLP_InputNeuronsArray->AddAt(0,i);
2059
2060 cout << "void EdbShowAlgE::CreateInputNeuronsArray_YY()...done. " << endl;
2061 return;
2062}

◆ CreateMLPArray()

void EdbShowAlgE::CreateMLPArray ( )
2148{
2149 if (!eMLPArray) eMLPArray = new TObjArray(10);
2150 return;
2151}
EdbShowAlgE::eMLPArray
TObjArray * eMLPArray
Definition: EdbShowAlgE.h:63

◆ CreateMLPLayout()

TString EdbShowAlgE::CreateMLPLayout ( Int_t  NrOfANNInputNeurons)
private
1300{
1301 TString layout;
1302
1303 for (Int_t i=1; i<NrOfANNInputNeurons; ++i) layout += "@eInANN["+TString(Form("%d",i))+"],";
1304 layout +="@eInANN["+TString(Form("%d",NrOfANNInputNeurons))+"]:"+TString(Form("%d",NrOfANNInputNeurons+1))+":"+TString(Form("%d",NrOfANNInputNeurons));
1305
1306 //layout+=":eInANN[0]/1000"; //if (Geant4) ([E]=MeV) //else layout+=":inANN[0]"; (Geant3) ([E]=GeV)
1307 layout+=":eInfoANN[0]/1000"; //if (Geant4) ([E]=MeV) //else layout+=":inANN[0]"; (Geant3) ([E]=GeV)
1308
1309 return layout;
1310}

◆ CreateMLPLayout_XX()

void EdbShowAlgE::CreateMLPLayout_XX ( Int_t  nr)
1800{
1801 // This Implementation if now depending on the parametrization:
1802 // Either it is depending on the numer of reconstructed plates;
1803 // FJ,LT,YY it is.
1804 // YC,JC,XX it is not.
1805
1806 if (gEDBDEBUGLEVEL>2) cout << "void EdbShowAlgE::CreateMLPLayout_XX("<<nr<<") " << endl;
1807
1808// cout << "eShowAlgEPlateBinning->GetSize() " << eShowAlgEPlateBinning->GetSize() <<endl;
1809// cout << "eMLPArrayN " <<eMLPArrayN <<endl;
1810// cout << "eShowAlgEPlateBinning->At(nr) " << eShowAlgEPlateBinning->At(nr) <<endl;
1811// cout << "eMLP_InputNeuronsArray->At(nr) " << eMLP_InputNeuronsArray->At(nr) <<endl;
1812
1813
1814 // Here the important change it acoording to the needs of the parametrization:
1815 Int_t NrOfANNInputNeurons=eMLP_InputNeuronsArray->At(nr);
1816 //XX:
1817 NrOfANNInputNeurons=7;
1818 //FJ:
1819 //NrOfANNInputNeurons=eMLP_InputNeuronsArray->At(nr)+1+2+2;
1820 // Simple Test
1821 //NrOfANNInputNeurons=eMLP_InputNeuronsArray->At(nr);
1822
1823
1824 TString layout;
1825 for (Int_t i=1; i<NrOfANNInputNeurons; ++i) layout += "@eInANN["+TString(Form("%d",i))+"],";
1826 layout +="@eInANN["+TString(Form("%d",NrOfANNInputNeurons))+"]:"+TString(Form("%d",NrOfANNInputNeurons+1))+":"+TString(Form("%d",NrOfANNInputNeurons));
1827 //layout+=":eInANN[0]/1000"; //if (Geant4) ([E]=MeV) //else layout+=":inANN[0]"; (Geant3) ([E]=GeV)
1828 layout+=":eInfoANN[0]/1000"; //if (Geant4) ([E]=MeV) //else layout+=":inANN[0]"; (Geant3) ([E]=GeV)
1829
1830 if (!eMLP_LayoutArray) eMLP_LayoutArray = new TObjArray(eMLP_InputNeuronsArrayN);
1831
1832
1833 // Add layout to the arrays of layouts:
1834 // -------------
1835 // TString
1836 // Basic string class.
1837 // Cannot be stored in a TCollection... use TObjString instead.
1838 // -------------
1839 eMLP_Layout=layout;
1840 TObjString* objstr = new TObjString(eMLP_Layout.Data());
1841 eMLP_LayoutArray->Add(objstr);
1842
1843
1844
1845 if (gEDBDEBUGLEVEL>3) eMLP_LayoutArray->Print();
1846
1847 if (gEDBDEBUGLEVEL>2) cout << "void EdbShowAlgE::CreateMLPLayout_XX("<<nr<<") ...done." << endl;
1848 return;
1849}
EdbShowAlgE::eMLP_Layout
TString eMLP_Layout
Definition: EdbShowAlgE.h:102
EdbShowAlgE::eMLP_LayoutArray
TObjArray * eMLP_LayoutArray
Definition: EdbShowAlgE.h:103
gEDBDEBUGLEVEL
gEDBDEBUGLEVEL
Definition: energy.C:7

◆ CreateMLPLayout_YC()

void EdbShowAlgE::CreateMLPLayout_YC ( Int_t  nr)
1741{
1742 // This Implementation if now depending on the parametrization:
1743 // Either it is depending on the numer of reconstructed plates;
1744 // FJ,LT,YY it is.
1745 // YC,JC,XX it is not.
1746
1747 if (gEDBDEBUGLEVEL>2) cout << "void EdbShowAlgE::CreateMLPLayout_YC("<<nr<<") " << endl;
1748
1749// cout << "eShowAlgEPlateBinning->GetSize() " << eShowAlgEPlateBinning->GetSize() <<endl;
1750// cout << "eMLPArrayN " <<eMLPArrayN <<endl;
1751// cout << "eShowAlgEPlateBinning->At(nr) " << eShowAlgEPlateBinning->At(nr) <<endl;
1752// cout << "eMLP_InputNeuronsArray->At(nr) " << eMLP_InputNeuronsArray->At(nr) <<endl;
1753
1754
1755 // Here the important change it acoording to the needs of the parametrization:
1756 Int_t NrOfANNInputNeurons=0;
1757 //XX:
1758 //NrOfANNInputNeurons=7;
1759 //FJ:
1760 //NrOfANNInputNeurons=eMLP_InputNeuronsArray->At(nr)+1+2+2;
1761 // Simple Test
1762 //NrOfANNInputNeurons=eMLP_InputNeuronsArray->At(nr);
1763 //YY:
1764 //NrOfANNInputNeurons=1+5+eMLP_InputNeuronsArray->At(nr);
1765 //YC:
1766 NrOfANNInputNeurons=1+2+2;
1767 //---
1768
1769 TString layout;
1770 for (Int_t i=1; i<NrOfANNInputNeurons; ++i) layout += "@eInANN["+TString(Form("%d",i))+"],";
1771 layout +="@eInANN["+TString(Form("%d",NrOfANNInputNeurons))+"]:"+TString(Form("%d",NrOfANNInputNeurons+1))+":"+TString(Form("%d",NrOfANNInputNeurons));
1772 //layout+=":eInANN[0]/1000"; //if (Geant4) ([E]=MeV) //else layout+=":inANN[0]"; (Geant3) ([E]=GeV)
1773 layout+=":eInfoANN[0]/1000"; //if (Geant4) ([E]=MeV) //else layout+=":inANN[0]"; (Geant3) ([E]=GeV)
1774
1775 if (!eMLP_LayoutArray) eMLP_LayoutArray = new TObjArray(eMLP_InputNeuronsArrayN);
1776
1777
1778 // Add layout to the arrays of layouts:
1779 // -------------
1780 // TString
1781 // Basic string class.
1782 // Cannot be stored in a TCollection... use TObjString instead.
1783 // -------------
1784 eMLP_Layout=layout;
1785 TObjString* objstr = new TObjString(eMLP_Layout.Data());
1786 eMLP_LayoutArray->Add(objstr);
1787
1788
1789 if (gEDBDEBUGLEVEL>3) eMLP_LayoutArray->Print();
1790 if (gEDBDEBUGLEVEL>3) cout << "void EdbShowAlgE::CreateMLPLayout_YC layout = " << layout << endl;
1791
1792 if (gEDBDEBUGLEVEL>2) cout << "void EdbShowAlgE::CreateMLPLayout_YC("<<nr<<") ...done." << endl;
1793 return;
1794}

◆ CreateMLPLayout_YY()

void EdbShowAlgE::CreateMLPLayout_YY ( Int_t  nr)
1856{
1857 // This Implementation if now depending on the parametrization:
1858 // Either it is depending on the numer of reconstructed plates;
1859 // FJ,LT,YY it is.
1860 // YC,JC,XX it is not.
1861
1862 if (gEDBDEBUGLEVEL>2) cout << "void EdbShowAlgE::CreateMLPLayout_YY("<<nr<<") " << endl;
1863
1864// cout << "eShowAlgEPlateBinning->GetSize() " << eShowAlgEPlateBinning->GetSize() <<endl;
1865// cout << "eMLPArrayN " <<eMLPArrayN <<endl;
1866// cout << "eShowAlgEPlateBinning->At(nr) " << eShowAlgEPlateBinning->At(nr) <<endl;
1867// cout << "eMLP_InputNeuronsArray->At(nr) " << eMLP_InputNeuronsArray->At(nr) <<endl;
1868
1869
1870 // Here the important change it acoording to the needs of the parametrization:
1871 Int_t NrOfANNInputNeurons=0;
1872 //XX:
1873 //NrOfANNInputNeurons=7;
1874 //FJ:
1875 //NrOfANNInputNeurons=eMLP_InputNeuronsArray->At(nr)+1+2+2;
1876 // Simple Test
1877 //NrOfANNInputNeurons=eMLP_InputNeuronsArray->At(nr);
1878 //YY:
1879 NrOfANNInputNeurons=1+5+eMLP_InputNeuronsArray->At(nr);
1880 //---
1881
1882
1883 TString layout;
1884 for (Int_t i=1; i<NrOfANNInputNeurons; ++i) layout += "@eInANN["+TString(Form("%d",i))+"],";
1885 layout +="@eInANN["+TString(Form("%d",NrOfANNInputNeurons))+"]:"+TString(Form("%d",NrOfANNInputNeurons+1))+":"+TString(Form("%d",NrOfANNInputNeurons));
1886 //layout+=":eInANN[0]/1000"; //if (Geant4) ([E]=MeV) //else layout+=":inANN[0]"; (Geant3) ([E]=GeV)
1887 layout+=":eInfoANN[0]/1000"; //if (Geant4) ([E]=MeV) //else layout+=":inANN[0]"; (Geant3) ([E]=GeV)
1888
1889 if (!eMLP_LayoutArray) eMLP_LayoutArray = new TObjArray(eMLP_InputNeuronsArrayN);
1890
1891
1892 // Add layout to the arrays of layouts:
1893 // -------------
1894 // TString
1895 // Basic string class.
1896 // Cannot be stored in a TCollection... use TObjString instead.
1897 // -------------
1898 eMLP_Layout=layout;
1899 TObjString* objstr = new TObjString(eMLP_Layout.Data());
1900 eMLP_LayoutArray->Add(objstr);
1901
1902
1903 if (gEDBDEBUGLEVEL>3) eMLP_LayoutArray->Print();
1904 if (gEDBDEBUGLEVEL>3) cout << "void EdbShowAlgE::CreateMLPLayout_YY layout = " << layout << endl;
1905
1906 if (gEDBDEBUGLEVEL>2) cout << "void EdbShowAlgE::CreateMLPLayout_YY("<<nr<<") ...done." << endl;
1907 return;
1908}

◆ CreateMLPLayouts()

void EdbShowAlgE::CreateMLPLayouts ( )
1600{
1601 cout << "void EdbShowAlgE::CreateMLPLayouts() " << endl;
1602
1603 cout << "void EdbShowAlgE::CreateMLPLayouts() Calling fuction suited to the eShowAlgEParaName" << endl;
1604 cout << "void EdbShowAlgE::CreateMLPLayouts() Respectively to the eShowAlgEParaNr" << endl;
1605
1606 switch (eShowAlgEParaNr) {
1607 case 0 :
1608 CreateMLPLayouts_FJ();
1609 break;
1610 case 1 :
1611 CreateMLPLayouts_LT();
1612 break;
1613 case 2 :
1614 CreateMLPLayouts_YC();
1615 break;
1616 case 3 :
1617 CreateMLPLayouts_JC();
1618 break;
1619 case 4 :
1620 CreateMLPLayouts_XX();
1621 break;
1622 case 5 :
1623 CreateMLPLayouts_YY();
1624 break;
1625 default :
1626 CreateMLPLayouts_XX();
1627 break;
1628 }
1629
1630 return;
1631}
EdbShowAlgE::CreateMLPLayouts_YC
void CreateMLPLayouts_YC()
Definition: EdbShowAlgE.cxx:1660
EdbShowAlgE::CreateMLPLayouts_JC
void CreateMLPLayouts_JC()
Definition: EdbShowAlgE.cxx:1684
EdbShowAlgE::CreateMLPLayouts_LT
void CreateMLPLayouts_LT()
Definition: EdbShowAlgE.cxx:1650
EdbShowAlgE::CreateMLPLayouts_YY
void CreateMLPLayouts_YY()
Definition: EdbShowAlgE.cxx:1717
EdbShowAlgE::CreateMLPLayouts_XX
void CreateMLPLayouts_XX()
Definition: EdbShowAlgE.cxx:1694
EdbShowAlgE::CreateMLPLayouts_FJ
void CreateMLPLayouts_FJ()
Definition: EdbShowAlgE.cxx:1640

◆ CreateMLPLayouts_FJ()

void EdbShowAlgE::CreateMLPLayouts_FJ ( )
1641{
1642 cout << "void EdbShowAlgE::CreateMLPLayouts_FJ() " << endl;
1643 return;
1644}

◆ CreateMLPLayouts_JC()

void EdbShowAlgE::CreateMLPLayouts_JC ( )
1685{
1686 cout << "void EdbShowAlgE::CreateMLPLayouts_JC() " << endl;
1687 return;
1688}

◆ CreateMLPLayouts_LT()

void EdbShowAlgE::CreateMLPLayouts_LT ( )
1651{
1652 cout << "void EdbShowAlgE::CreateMLPLayouts_LT() " << endl;
1653 return;
1654}

◆ CreateMLPLayouts_XX()

void EdbShowAlgE::CreateMLPLayouts_XX ( )
1695{
1696 cout << "void EdbShowAlgE::CreateMLPLayouts_XX() " << endl;
1697
1698 if (!eMLP_InputNeuronsArray) CreateInputNeuronsArray();
1699
1700 FillInputNeuronsArray_XX();
1701
1702 // This is the loop for the number of eMLPs to be availible. For each plate number weightfile
1703 // a new layout is created.
1704 for (Int_t i=0; i<eMLPArrayN; ++i) CreateMLPLayout_XX(i);
1705
1706 // Now we know that it was done.
1707 // But no check if layouts are correctly, only the function steps have been done!
1708 eMLP_LayoutArrayIsDone=kTRUE;
1709
1710 cout << "void EdbShowAlgE::CreateMLPLayouts_XX()...done." << endl;
1711 return;
1712}
EdbShowAlgE::CreateInputNeuronsArray
void CreateInputNeuronsArray()
Definition: EdbShowAlgE.cxx:1915
EdbShowAlgE::CreateMLPLayout_XX
void CreateMLPLayout_XX(Int_t nr)
Definition: EdbShowAlgE.cxx:1799
EdbShowAlgE::eMLP_LayoutArrayIsDone
Bool_t eMLP_LayoutArrayIsDone
Definition: EdbShowAlgE.h:104
EdbShowAlgE::FillInputNeuronsArray_XX
void FillInputNeuronsArray_XX()
Definition: EdbShowAlgE.cxx:2020

◆ CreateMLPLayouts_YC()

void EdbShowAlgE::CreateMLPLayouts_YC ( )
1661{
1662 cout << "void EdbShowAlgE::CreateMLPLayouts_YY() " << endl;
1663
1664 if (!eMLP_InputNeuronsArray) CreateInputNeuronsArray();
1665
1666 FillInputNeuronsArray_YC();
1667
1668 // This is the loop for the number of eMLPs to be availible. For each plate number weightfile
1669 // a new layout is created.
1670 for (Int_t i=0; i<eMLPArrayN; ++i) CreateMLPLayout_YC(i);
1671
1672 // Now we know that it was done.
1673 // But no check if layouts are correctly, only the function steps have been done!
1674 eMLP_LayoutArrayIsDone=kTRUE;
1675
1676 cout << "void EdbShowAlgE::CreateMLPLayouts_YC()...done." << endl;
1677 return;
1678}
EdbShowAlgE::FillInputNeuronsArray_YC
void FillInputNeuronsArray_YC()
Definition: EdbShowAlgE.cxx:1975
EdbShowAlgE::CreateMLPLayout_YC
void CreateMLPLayout_YC(Int_t nr)
Definition: EdbShowAlgE.cxx:1740

◆ CreateMLPLayouts_YY()

void EdbShowAlgE::CreateMLPLayouts_YY ( )
1718{
1719 cout << "void EdbShowAlgE::CreateMLPLayouts_YY() " << endl;
1720
1721 if (!eMLP_InputNeuronsArray) CreateInputNeuronsArray();
1722
1723 FillInputNeuronsArray_YY();
1724
1725 // This is the loop for the number of eMLPs to be availible. For each plate number weightfile
1726 // a new layout is created.
1727 for (Int_t i=0; i<eMLPArrayN; ++i) CreateMLPLayout_YY(i);
1728
1729 // Now we know that it was done.
1730 // But no check if layouts are correctly, only the function steps have been done!
1731 eMLP_LayoutArrayIsDone=kTRUE;
1732
1733 return;
1734}
EdbShowAlgE::CreateMLPLayout_YY
void CreateMLPLayout_YY(Int_t nr)
Definition: EdbShowAlgE.cxx:1855
EdbShowAlgE::FillInputNeuronsArray_YY
void FillInputNeuronsArray_YY()
Definition: EdbShowAlgE.cxx:2068

◆ CreateMLPs()

void EdbShowAlgE::CreateMLPs ( )
2094{
2095 // Since the layout if sixed by the parametrisation, this function
2096 // can be equall for all types of parametrisations:
2097
2098 // Check if layouts have been done:
2099 if (!eMLP_LayoutArrayIsDone) {
2100 cout << "void EdbShowAlgE::CreateMLPs() ERROR: eMLP_LayoutArrayIsDone HAS NOT BEEN DONE. RETURN!" << endl;
2101 return;
2102 }
2103
2104
2105 Log(3,"EdbShowAlgE::CreateMLPs()","Since the layout if sixed by the parametrisation, this function ");
2106 Log(3,"EdbShowAlgE::CreateMLPs()","can be equall for all types of parametrisations.");
2107 Log(3,"EdbShowAlgE::CreateMLPs()","It is now only different for some parametrisations by the number of plates. ");
2108 Log(2,"EdbShowAlgE::CreateMLPs()","If you see now a message like:");
2109 Log(2,"EdbShowAlgE::CreateMLPs()","Data not set. Cannot define datasets.");
2110 Log(2,"EdbShowAlgE::CreateMLPs()","Then this means that we have not yet defined the eMLP_SimpleTree for the ANN. This has to be done later. In any case the ANN can be created without it.");
2111
2112
2113 Int_t tmpint;
2114
2115 if (!eMLPArray) CreateMLPArray();
2116
2117
2118 for (Int_t i=0; i<eMLP_InputNeuronsArrayN; ++i) {
2119
2120 tmpint= eMLP_InputNeuronsArray->At(i);
2121
2122 TString layout = ((TObjString*)eMLP_LayoutArray->At(i))->GetString();
2123
2124 eMLP_Simple = new TMultiLayerPerceptron(layout, eMLP_SimpleTree,"Entry$%2","(Entry$+1)%2");
2125
2126 //cout << " eMLP_Simple = " << eMLP_Simple << endl;
2127
2128 eMLPArray->Add(eMLP_Simple);
2129 }
2130
2131 cout << " eMLP_Simple = " << eMLP_Simple << endl;
2132 cout << " eMLPArray->GetEntries() = " << eMLPArray->GetEntries() << endl;
2133 cout << " eMLPArrayN = " << eMLPArrayN << endl;
2134
2135 // This should have no effect since eMLPArrayN is already filled with correct value;
2136 SetMLPArrayN(eMLPArray->GetEntries());
2137
2138 cout << " eMLPArrayN = " << eMLPArrayN << endl;
2139
2140 return;
2141}
EdbShowAlgE::SetMLPArrayN
void SetMLPArrayN(Int_t MLPArrayN)
Definition: EdbShowAlgE.h:187
EdbShowAlgE::CreateMLPArray
void CreateMLPArray()
Definition: EdbShowAlgE.cxx:2147

◆ CreateStandardMLPArrayN()

void EdbShowAlgE::CreateStandardMLPArrayN ( )
1375{
1376 cout << "void EdbShowAlgE::CreateStandardMLPArrayN() " << endl;
1377 cout << "void EdbShowAlgE::CreateStandardMLPArrayN() According to the standard weightfile strings" << endl;
1378 cout << "void EdbShowAlgE::CreateStandardMLPArrayN() (weights_energy_11.txt,...,weights_energy_46.txt)." << endl;
1379
1380 if (!eShowAlgEPlateBinning) ResetPlateBinning();
1381
1382 // Set Number that tells how many MLPs are instantiated.
1383 // 11,16,21,26,31,36,41,46 (51,56)
1384 // the array itsself contains 10 entries but for 51 and 56 plate the weightfiles do not exist yet.
1385 eMLPArrayN=10;
1386
1387 // These two variables are always same!
1388 eMLP_InputNeuronsArrayN=eMLPArrayN;
1389
1390// cout << "------"<< endl;
1391// cout << "eMLPArrayN = " << eMLPArrayN << endl;
1392// cout << "eMLP_InputNeuronsArrayN = " << eMLP_InputNeuronsArrayN << endl;
1393
1394 cout << "void EdbShowAlgE::CreateStandardMLPArrayN() done." << endl;
1395 return;
1396}
EdbShowAlgE::ResetPlateBinning
void ResetPlateBinning()
Definition: EdbShowAlgE.cxx:1316

◆ CreateStandardWeightFileStrings()

void EdbShowAlgE::CreateStandardWeightFileStrings ( )
1403{
1404 cout << "void EdbShowAlgE::CreateStandardWeightFileStrings() " << endl;
1405
1406
1407 if (!eShowAlgEPlateBinning) ResetPlateBinning();
1408
1409 // There can be some cases when default par is read in and eShowAlgEPlateBinning was creadted, but due to wrong input settings
1410 // of default.par. In this case, eShowAlgEPlateBinning is there, but Cleared (filled with zeros).
1411 // So if this, we have to ResetPlateBinning again:
1412 if (eShowAlgEPlateBinning && eShowAlgEPlateBinning->At(0)==0 && eShowAlgEPlateBinning->At(9)==0 ) {
1413 cout << "void EdbShowAlgE::CreateStandardWeightFileStrings() ResetPlateBinning again, cause eShowAlgEPlateBinning is filled with zeros." << endl;
1414 ResetPlateBinning();
1415 }
1416
1417 cout << "void EdbShowAlgE::CreateStandardWeightFileStrings() ResetPlateBinning done" << endl;
1418
1419 // Create Standard weight file strings:
1420 // weight_energy_YY_"i".txt, i=11,16,21,26,31,36,41,46,51,56
1421 // .Data() to convert to const char
1422 Int_t tmpint;
1423 for (Int_t i=0; i<10; ++i) {
1424 tmpint=eShowAlgEPlateBinning->At(i);
1425 eMLP_WeightFileNameArray[i]=TString(Form("weights_energy_%s_%d.txt",eShowAlgEParaName.Data(),tmpint));
1426 //cout << eMLP_WeightFileNameArray[i] << endl;
1427 }
1428
1429 cout << "void EdbShowAlgE::CreateStandardWeightFileStrings() ...done." << endl;
1430 return;
1431}

◆ DumpWeights()

void EdbShowAlgE::DumpWeights ( )
669{
670 if (!eMLP_WriteWeightFile) {
671 cout << "void EdbShowAlgE::DumpWeights() eMLP_WriteWeightFile="<< eMLP_WriteWeightFile << " . So write nothing. RETURN"<<endl;
672 return;
673 }
674 if (GetMLPWeightFileName()==""||GetMLPWeightFileName()=="EMPTY") SetMLPWeightFileName("weights.txt");
675
676 cout << "void EdbShowAlgE::DumpWeights() "<< endl;
677 cout << eMLP_WeightFileName << endl;
678 eMLP_Simple->DumpWeights(eMLP_WeightFileName);
679 cout << "void EdbShowAlgE::DumpWeights() ... done."<< endl;
680 return;
681}
EdbShowAlgE::GetMLPWeightFileName
TString GetMLPWeightFileName()
Definition: EdbShowAlgE.h:225
EdbShowAlgE::SetMLPWeightFileName
void SetMLPWeightFileName(TString weightfilename)
Definition: EdbShowAlgE.h:178
EdbShowAlgE::eMLP_WeightFileName
TString eMLP_WeightFileName
Definition: EdbShowAlgE.h:108
EdbShowAlgE::eMLP_WriteWeightFile
Bool_t eMLP_WriteWeightFile
Definition: EdbShowAlgE.h:107

◆ Execute()

void EdbShowAlgE::Execute ( )
virtual
249{
250 return;
251}

◆ Fill_MLPTree()

void EdbShowAlgE::Fill_MLPTree ( )

DRPECIATED !!!!!!!!!!!!

Get Suited Parametrization: TO BE ADAPTED FOR ALL PARAMETRIZATIONS !!!!!

1073{
1074
1076
1077 return;
1078
1079
1080 cout << "EdbShowAlgE::Fill_MLPTree()" << endl;
1081
1082 cout << "EdbShowAlgE::Fill_MLPTree() // Get correct Parametrisation. ONLY XX IS SUPPORTED AT THE RIGHT MOMENT !!!" << endl;
1083
1084 // Get correct Parametrisation. ONLY XX IS SUPPORTED AT THE RIGHT MOMENT !!!
1085
1086
1087 // Create and Clear Tree:
1088 if (!eMLP_SimpleTree) Create_MLPTree();
1089 eMLP_SimpleTree->Clear();
1090
1091
1092
1093 EdbShowerP* show=0;
1094 for (int i=0; i<eRecoShowerArrayN; i++ ) {
1095 cout << "EdbShowAlgE::TESTTRAIN i= " << i << endl;
1096 show=(EdbShowerP*)eRecoShowerArray->At(i);
1097 show->PrintNice();
1098 if (! show->GetParametrisationIsDone(4)) show->BuildParametrisation(4);
1099
1102 EdbShowerP::Para_XX ePara_XX;
1103 ePara_XX=show->GetPara_XX();
1104
1105 // Get each Variables: InputType
1106 // Energy
1107 eInANN[0]=show->GetSegment(0)->P();
1108 // This can be different in later modifications (for example gamma Showers, where the Energy is not defined by single BT P(), but
1109 // by both starting BT close to each other (e+,e-)pair.
1110
1111 // Variables: TestingType
1112 eInANN[1]=ePara_XX.nseg;
1113 eInANN[2]=ePara_XX.Mean_ProfileLongitudinal;
1114 eInANN[3]=ePara_XX.RMS_ProfileLongitudinal;
1115 eInANN[4]=ePara_XX.Max_ProfileLongitudinal;
1116 eInANN[5]=ePara_XX.Mean_ProfileTransversal;
1117 eInANN[6]=ePara_XX.RMS_ProfileTransversal;
1118 eInANN[7]=ePara_XX.Max_ProfileTransversal;
1119
1120 if (gEDBDEBUGLEVEL>3) {
1121 for (int i=0; i<eRecoShowerArrayN; i++ ) cout << "EdbShowAlgE::_XX() eInANN["<<i<<"]= " << eInANN[i] << endl;
1122 cout << "EdbShowAlgE:: ...."<<endl;
1123 cout << "EdbShowAlgE:: ...."<<endl;
1124 cout << "EdbShowAlgE:: eInANN[69]= " << eInANN[69] << endl;
1125 }
1126
1127 eMLP_SimpleTree->Fill();
1128 }
1129
1130
1131 if (gEDBDEBUGLEVEL>2) cout << "EdbShowAlgE:: Fill_MLPTree done." << endl;
1132 if (gEDBDEBUGLEVEL>2) eMLP_SimpleTree->Print();
1133
1134
1135 cout << "EdbShowAlgE::Fill_MLPTree()...done." << endl;
1136 return;
1137}
EdbSegP::P
Float_t P() const
Definition: EdbSegP.h:152
EdbShowAlgE::Create_MLPTree
void Create_MLPTree()
Definition: EdbShowAlgE.cxx:1047
EdbShowerP
Definition: EdbShowerP.h:28
EdbShowerP::GetSegment
EdbSegP * GetSegment(int i) const
Definition: EdbShowerP.h:435
EdbShowerP::PrintNice
void PrintNice()
Definition: EdbShowerP.cxx:2339
EdbShowerP::GetParametrisationIsDone
Bool_t GetParametrisationIsDone(Int_t nr) const
Definition: EdbShowerP.h:546
EdbShowerP::BuildParametrisation
void BuildParametrisation(Int_t ParaNr)
Definition: EdbShowerP.cxx:1190
EdbShowerP::GetPara_XX
Para_XX GetPara_XX() const
Definition: EdbShowerP.h:563
EdbShowerP::Para_XX
Definition: EdbShowerP.h:200
EdbShowerP::Para_XX::Max_ProfileLongitudinal
Float_t Max_ProfileLongitudinal
Definition: EdbShowerP.h:206
EdbShowerP::Para_XX::nseg
Int_t nseg
Definition: EdbShowerP.h:203
EdbShowerP::Para_XX::Max_ProfileTransversal
Float_t Max_ProfileTransversal
Definition: EdbShowerP.h:209
EdbShowerP::Para_XX::RMS_ProfileLongitudinal
Float_t RMS_ProfileLongitudinal
Definition: EdbShowerP.h:205
EdbShowerP::Para_XX::Mean_ProfileLongitudinal
Float_t Mean_ProfileLongitudinal
Definition: EdbShowerP.h:204
EdbShowerP::Para_XX::RMS_ProfileTransversal
Float_t RMS_ProfileTransversal
Definition: EdbShowerP.h:208
EdbShowerP::Para_XX::Mean_ProfileTransversal
Float_t Mean_ProfileTransversal
Definition: EdbShowerP.h:207

◆ FillInputHistosPara_XX()

void EdbShowAlgE::FillInputHistosPara_XX ( )
private
1270{
1271 // Reset
1272 for (Int_t j=0; j<70; j++) eInANN_Hist[j]->Reset();
1273 //Set bins:
1274 SetInputHistosPara_XX();
1275
1276 //Fill with entries from the tree:
1277
1278
1279 for (int i=0; i<eMLP_SimpleTree->GetEntries(); ++i) {
1280 eMLP_SimpleTree->GetEntry(i);
1281 // Fill InputHistograms now...
1282 for (Int_t j=0; j<70; j++) {
1283 eInANN_Hist[j]->Fill(eInANN[j]);
1284 }
1285 }
1286
1287
1288 return;
1289}
EdbShowAlgE::SetInputHistosPara_XX
void SetInputHistosPara_XX()
Definition: EdbShowAlgE.cxx:1223

◆ FillInputNeuronsArray_XX()

void EdbShowAlgE::FillInputNeuronsArray_XX ( )
2021{
2022 cout << "void EdbShowAlgE::FillInputNeuronsArray_XX() " << endl;
2023 // Restore Standard Plate Binning
2024 if (!eMLP_InputNeuronsArray) return;
2025
2026 // Not more than 10 values of weightfiles will be allowed!
2027 Int_t tmpint;
2028
2029 for (Int_t i=0; i<10; ++i) {
2030 tmpint=eShowAlgEPlateBinning->At(i);
2031 eMLP_InputNeuronsArray->AddAt(tmpint,i);
2032 //cout << i << " " << tmpint << endl;
2033 }
2034
2035 cout << "void EdbShowAlgE::FillInputNeuronsArray_XX()...done. " << endl;
2036 return;
2037}

◆ FillInputNeuronsArray_YC()

void EdbShowAlgE::FillInputNeuronsArray_YC ( )
1976{
1977 cout << "void EdbShowAlgE::FillInputNeuronsArray_YC() " << endl;
1978 // Restore Standard Plate Binning
1979 if (!eMLP_InputNeuronsArray) return;
1980
1981 // Not more than 10 values of weightfiles will be allowed!
1982 Int_t tmpint;
1983
1984 for (Int_t i=0; i<10; ++i) {
1985 tmpint=eShowAlgEPlateBinning->At(i);
1986 eMLP_InputNeuronsArray->AddAt(tmpint,i);
1987 //cout << i << " " << tmpint << endl;
1988 }
1989
1990 cout << "void EdbShowAlgE::FillInputNeuronsArray_YC()...done. " << endl;
1991 return;
1992}

◆ FillInputNeuronsArray_YY()

void EdbShowAlgE::FillInputNeuronsArray_YY ( )
2069{
2070 cout << "void EdbShowAlgE::FillInputNeuronsArray_YY() " << endl;
2071 // Restore Standard Plate Binning
2072 if (!eMLP_InputNeuronsArray) return;
2073
2074 // Not more than 10 values of weightfiles will be allowed!
2075 Int_t tmpint;
2076
2077 for (Int_t i=0; i<10; ++i) {
2078 tmpint=eShowAlgEPlateBinning->At(i);
2079 eMLP_InputNeuronsArray->AddAt(tmpint,i);
2080 //cout << i << " " << tmpint << endl;
2081 }
2082
2083 cout << "void EdbShowAlgE::FillInputNeuronsArray_YY()...done. " << endl;
2084 return;
2085}

◆ FillMLP_SimpleTree()

void EdbShowAlgE::FillMLP_SimpleTree ( )
541{
542 cout << "EdbShowAlgE::FillMLP_SimpleTree()" << endl;
543
544
545 cout <<"===============" << GetMLPInputFileName() << endl;
546
547 // ---------
548 if (!eMLP_SimpleTree) {
549 cout << "EdbShowAlgE::FillMLP_SimpleTree() NO eMLP_SimpleTree. Create now. " << endl;
550 Create_MLPTree();
551 }
552
553
554 cout <<"===============" << GetMLPInputFileName() << endl;
555
556 if (!eMLPArray) {
557 cout << "EdbShowAlgE::FillMLP_SimpleTree() NO eMLPArray " << endl;
558 return;
559 }
560
561
562 cout <<"===============" << GetMLPInputFileName() << endl;
563
564 // Do this in case that eMLP_SimpleTree is read by ReadMLPInputFile
565 if (eMLP_ReadMLPInputFile) {
566 cout << "EdbShowAlgE::FillMLP_SimpleTree() Read eMLP_SimpleTree from file:" << endl;
567 ReadMLPInputFile();
568 return;
569 }
570
571 cout << "EdbShowAlgE::FillMLP_SimpleTree() Call FillMLP_SimpleTree___() now:"<< endl;
572 if (eShowAlgEParaNr==2) FillMLP_SimpleTree_YC();
573 if (eShowAlgEParaNr==4) FillMLP_SimpleTree_XX();
574 if (eShowAlgEParaNr==5) FillMLP_SimpleTree_YY();
575
576
577 cout <<"===============" << GetMLPInputFileName() << endl;
578
579 cout << "EdbShowAlgE::FillMLP_SimpleTree() done." << endl;
580 return;
581}
EdbShowAlgE::FillMLP_SimpleTree_YY
void FillMLP_SimpleTree_YY()
Definition: EdbShowAlgE.cxx:629
EdbShowAlgE::GetMLPInputFileName
TString GetMLPInputFileName()
Definition: EdbShowAlgE.h:228
EdbShowAlgE::ReadMLPInputFile
void ReadMLPInputFile()
Definition: EdbShowAlgE.cxx:527
EdbShowAlgE::FillMLP_SimpleTree_XX
void FillMLP_SimpleTree_XX()
Definition: EdbShowAlgE.cxx:607
EdbShowAlgE::eMLP_ReadMLPInputFile
Bool_t eMLP_ReadMLPInputFile
Definition: EdbShowAlgE.h:119
EdbShowAlgE::FillMLP_SimpleTree_YC
void FillMLP_SimpleTree_YC()
Definition: EdbShowAlgE.cxx:584

◆ FillMLP_SimpleTree_XX() [1/2]

void EdbShowAlgE::FillMLP_SimpleTree_XX ( )
608{
609 cout << "EdbShowAlgE::FillMLP_SimpleTree_XX()" << endl;
610
611
612
613 // Do this in case that eMLP_SimpleTree is filled by eRecoShowerArrayN
614 EdbShowerP* show=0;
615 for (int i=0; i<eRecoShowerArrayN; i++ ) {
616 //cout << "EdbShowAlgE::FILL i= " << i << endl;
617 show=(EdbShowerP*)eRecoShowerArray->At(i);
618 FillMLP_SimpleTree_XX(show,eShowAlgEParaNr);
619 }
620 eMLP_SimpleTree->Print();
621 cout << "eMLP_SimpleTree>GetEntries():"<< eMLP_SimpleTree->GetEntries() << endl;
622
623 eMLP_SimpleTree->Show(eMLP_SimpleTree->GetEntries()-1);
624 cout << "EdbShowAlgE::FillMLP_SimpleTree_XX()...done." << endl;
625 return;
626}

◆ FillMLP_SimpleTree_XX() [2/2]

void EdbShowAlgE::FillMLP_SimpleTree_XX ( EdbShowerP *  show,
Int_t  ShowAlgEParaNr 
)
2340 {
2341
2342 if (! show->GetParametrisationIsDone(ShowAlgEParaNr)) show->BuildParametrisation(ShowAlgEParaNr);
2343
2344 // Get Suited Parametrization:
2345 EdbShowerP::Para_XX ePara_XX;
2346 ePara_XX=show->GetPara_XX();
2347
2348 // Get each Variables: InputType
2349 // Energy
2350 eInANN[0]=show->GetSegment(0)->P();
2351 // This can be different in later modifications (for example gamma Showers, where the Energy is not defined by single BT P(), but
2352 // by both starting BT close to each other (e+,e-)pair.
2353
2354 // Variables: TestingType
2355 eInANN[1]=ePara_XX.nseg;
2356 eInANN[2]=ePara_XX.Mean_ProfileLongitudinal;
2357 eInANN[3]=ePara_XX.RMS_ProfileLongitudinal;
2358 eInANN[4]=ePara_XX.Max_ProfileLongitudinal;
2359 eInANN[5]=ePara_XX.Mean_ProfileTransversal;
2360 eInANN[6]=ePara_XX.RMS_ProfileTransversal;
2361 eInANN[7]=ePara_XX.Max_ProfileTransversal;
2362
2363 // Fill also related eInfoANN variables:
2364 eInfoANN[0]=show->GetSegment(0)->P();
2365 eInfoANN[1]=show->Npl();
2366
2367
2368 if (gEDBDEBUGLEVEL>3) {
2369 for (int i=0; i<13; i++ ) cout << "EdbShowAlgE::_XX() eInANN["<<i<<"]= " << eInANN[i] << endl;
2370 cout << "EdbShowAlgE:: ...."<<endl;
2371 for (int i=70-13; i<70; i++ ) cout << "EdbShowAlgE::_XX() eInANN["<<i<<"]= " << eInANN[i] << endl;
2372
2373 cout << "EdbShowAlgE::_XX() eInfoANN[0]= " << eInfoANN[0] << endl;
2374 cout << "EdbShowAlgE::_XX() eInfoANN[1]= " << eInfoANN[1] << endl;
2375 }
2376
2377 // Fill InputHistograms now...
2378 for (Int_t j=0; j<70; j++) {
2379 eInANN_Hist[j]->Fill(eInANN[j]);
2380 }
2381
2382
2383 if (gEDBDEBUGLEVEL>3) {
2384 cout << "EdbShowAlgE:: eMLP_SimpleTree " << eMLP_SimpleTree << endl;
2385 cout << "EdbShowAlgE:: eMLP_SimpleTree->ClassName() " << eMLP_SimpleTree->ClassName() << endl;
2386 cout << "EdbShowAlgE:: eMLP_SimpleTree->GetEntries():" << endl;
2387 cout << eMLP_SimpleTree->GetEntries() <<endl;
2388 }
2389
2390
2391 if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlgE:: eMLP_SimpleTree->Fill() do:" << endl;
2392// gDebug=2;
2393 eMLP_SimpleTree->Fill();
2394// gDebug=1;
2395 if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlgE:: eMLP_SimpleTree->Fill() done." << endl;
2396
2397 return;
2398}
EdbShowerP::Npl
Int_t Npl() const
Definition: EdbShowerP.h:409

◆ FillMLP_SimpleTree_YC() [1/2]

void EdbShowAlgE::FillMLP_SimpleTree_YC ( )
585{
586 cout << "EdbShowAlgE::FillMLP_SimpleTree_YC()" << endl;
587
588
589
590 // Do this in case that eMLP_SimpleTree is filled by eRecoShowerArrayN
591 EdbShowerP* show=0;
592 for (int i=0; i<eRecoShowerArrayN; i++ ) {
593 //cout << "EdbShowAlgE::FILL i= " << i << endl;
594 show=(EdbShowerP*)eRecoShowerArray->At(i);
595 FillMLP_SimpleTree_YC(show,eShowAlgEParaNr);
596 }
597 eMLP_SimpleTree->Print();
598 cout << "eMLP_SimpleTree>GetEntries():"<< eMLP_SimpleTree->GetEntries() << endl;
599
600 eMLP_SimpleTree->Show(eMLP_SimpleTree->GetEntries()-1);
601 cout << "EdbShowAlgE::FillMLP_SimpleTree_YC()...done." << endl;
602 return;
603}

◆ FillMLP_SimpleTree_YC() [2/2]

void EdbShowAlgE::FillMLP_SimpleTree_YC ( EdbShowerP *  show,
Int_t  ShowAlgEParaNr 
)
2270 {
2271
2272 if (! show->GetParametrisationIsDone(ShowAlgEParaNr)) show->BuildParametrisation(ShowAlgEParaNr);
2273
2274 // Get Suited Parametrization:
2275 EdbShowerP::Para_YC ePara_YC;
2276 ePara_YC=show->GetPara_YC();
2277
2278
2279 cout << "EdbShowerP::FillMLP_SimpleTree_YC() ePara_YC.nseg= " << ePara_YC.nseg << endl;
2280 cout << "EdbShowerP::FillMLP_SimpleTree_YC() ePara_YC.C1= " << ePara_YC.C1 << endl;
2281 cout << "EdbShowerP::FillMLP_SimpleTree_YC() ePara_YC.a1= " << ePara_YC.a1 << endl;
2282 cout << "EdbShowerP::FillMLP_SimpleTree_YC() ePara_YC.alpha= " << ePara_YC.alpha << endl;
2283 cout << "EdbShowerP::FillMLP_SimpleTree_YC() ePara_YC.nmax= " << ePara_YC.nmax << endl;
2284
2285 // Get each Variables: InputType
2286 // Energy
2287 eInANN[0]=show->GetSegment(0)->P();
2288 // This can be different in later modifications (for example gamma Showers, where the Energy is not defined by single BT P(), but
2289 // by both starting BT close to each other (e+,e-)pair.
2290
2291 // Variables: TestingType
2292 eInANN[1]=ePara_YC.nseg;
2293 eInANN[2]=ePara_YC.C1;
2294 eInANN[3]=ePara_YC.a1;
2295 eInANN[4]=ePara_YC.alpha;
2296 eInANN[5]=ePara_YC.nmax;
2297
2298 // Fill also related eInfoANN variables:
2299 eInfoANN[0]=show->GetSegment(0)->P();
2300 eInfoANN[1]=show->Npl();
2301
2302
2303 if (gEDBDEBUGLEVEL>3) {
2304 for (int i=0; i<13; i++ ) cout << "EdbShowAlgE::_YC() eInANN["<<i<<"]= " << eInANN[i] << endl;
2305 cout << "EdbShowAlgE:: ...."<<endl;
2306 for (int i=70-13; i<70; i++ ) cout << "EdbShowAlgE::_YC() eInANN["<<i<<"]= " << eInANN[i] << endl;
2307
2308 cout << "EdbShowAlgE::_YC() eInfoANN[0]= " << eInfoANN[0] << endl;
2309 cout << "EdbShowAlgE::_YC() eInfoANN[1]= " << eInfoANN[1] << endl;
2310 }
2311
2312 // Fill InputHistograms now...
2313 for (Int_t j=0; j<70; j++) {
2314 eInANN_Hist[j]->Fill(eInANN[j]);
2315 }
2316
2317
2318 if (gEDBDEBUGLEVEL>3) {
2319 cout << "EdbShowAlgE:: eMLP_SimpleTree " << eMLP_SimpleTree << endl;
2320 cout << "EdbShowAlgE:: eMLP_SimpleTree->ClassName() " << eMLP_SimpleTree->ClassName() << endl;
2321 cout << "EdbShowAlgE:: eMLP_SimpleTree->GetEntries():" << endl;
2322 cout << eMLP_SimpleTree->GetEntries() <<endl;
2323 }
2324
2325
2326 if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlgE:: eMLP_SimpleTree->Fill() do:" << endl;
2327// gDebug=2;
2328 eMLP_SimpleTree->Fill();
2329// gDebug=1;
2330 if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlgE:: eMLP_SimpleTree->Fill() done." << endl;
2331
2332 return;
2333}
EdbShowerP::GetPara_YC
Para_YC GetPara_YC() const
Definition: EdbShowerP.h:557
EdbShowerP::Para_YC
Definition: EdbShowerP.h:169
EdbShowerP::Para_YC::C1
Float_t C1
Definition: EdbShowerP.h:173
EdbShowerP::Para_YC::alpha
Float_t alpha
Definition: EdbShowerP.h:175
EdbShowerP::Para_YC::nseg
Int_t nseg
Definition: EdbShowerP.h:172
EdbShowerP::Para_YC::a1
Float_t a1
Definition: EdbShowerP.h:174
EdbShowerP::Para_YC::nmax
Float_t nmax
Definition: EdbShowerP.h:176

◆ FillMLP_SimpleTree_YY() [1/2]

void EdbShowAlgE::FillMLP_SimpleTree_YY ( )
630{
631 cout << "EdbShowAlgE::FillMLP_SimpleTree_YY()" << endl;
632
633 // Do this in case that eMLP_SimpleTree is filled by eRecoShowerArrayN
634 EdbShowerP* show=0;
635 for (int i=0; i<eRecoShowerArrayN; i++ ) {
636 //cout << "EdbShowAlgE::FILL i= " << i << endl;
637 show=(EdbShowerP*)eRecoShowerArray->At(i);
638 FillMLP_SimpleTree_YY(show,eShowAlgEParaNr);
639 }
640 eMLP_SimpleTree->Print();
641 cout << "eMLP_SimpleTree>GetEntries():"<< eMLP_SimpleTree->GetEntries() << endl;
642
643 eMLP_SimpleTree->Show(eMLP_SimpleTree->GetEntries()-1);
644 cout << "EdbShowAlgE::FillMLP_SimpleTree_YY()...done." << endl;
645 return;
646}

◆ FillMLP_SimpleTree_YY() [2/2]

void EdbShowAlgE::FillMLP_SimpleTree_YY ( EdbShowerP *  show,
Int_t  ShowAlgEParaNr 
)
2408{
2409
2410 cout << "DEBUGDEBUGDEBUGDEBUGDEBUGDEBUGDEBUGDEBUGDEBUGDEBUGDEBUGDEBUGDEBUGDEBUG"<<endl;
2411 if (! show->GetParametrisationIsDone(ShowAlgEParaNr)) {
2412 show->PrintNice();
2413 cout << "NO PARAMETRISATION DONE FOR show= " << show << "with " << show->N() << "Baestracks"<< endl;
2414 }
2415
2416 if (! show->GetParametrisationIsDone(ShowAlgEParaNr)) show->BuildParametrisation(ShowAlgEParaNr);
2417
2418 // Get Suited Parametrization:
2419 EdbShowerP::Para_YY ePara_YY;
2420 ePara_YY=show->GetPara_YY();
2421
2422 // Get each Variables: InputType
2423 // Energy
2424 eInANN[0]=show->GetSegment(0)->P();
2425 // This can be different in later modifications (for example gamma Showers, where the Energy is not defined by single BT P(), but
2426 // by both starting BT close to each other (e+,e-)pair.
2427
2428 // Variables: TestingType
2429 eInANN[1]=ePara_YY.nseg;
2430 for (Int_t j=0; j<5; j++) eInANN[2+j]=ePara_YY.ProfileTransversalBincontent[1+j]; // 0 is Underflow bin.
2431 for (Int_t j=0; j<57; j++) eInANN[7+j]=ePara_YY.ProfileLongitudinalBincontent[1+j]; // 0 is Underflow bin.
2432
2433 // Fill also related eInfoANN variables:
2434 eInfoANN[0]=show->GetSegment(0)->P();
2435 eInfoANN[1]=show->Npl();
2436
2437
2438 if (gEDBDEBUGLEVEL>3) {
2439 for (int i=0; i<13; i++ ) cout << "EdbShowAlgE::_YY() eInANN["<<i<<"]= " << eInANN[i] << endl;
2440 cout << "EdbShowAlgE:: ...."<<endl;
2441 for (int i=70-13; i<70; i++ ) cout << "EdbShowAlgE::_YY() eInANN["<<i<<"]= " << eInANN[i] << endl;
2442
2443 cout << "EdbShowAlgE::_YY() eInfoANN[0]= " << eInfoANN[0] << endl;
2444 cout << "EdbShowAlgE::_YY() eInfoANN[1]= " << eInfoANN[1] << endl;
2445 }
2446
2447
2448 if (gEDBDEBUGLEVEL>3) show->PrintBasics();
2449 if (gEDBDEBUGLEVEL>3) show->PrintNice();
2450 if (gEDBDEBUGLEVEL>3) ;
2451 show->PrintParametrisation_YY();
2452
2453
2454
2455 // Fill InputHistograms now...
2456 for (Int_t j=0; j<70; j++) {
2457 eInANN_Hist[j]->Fill(eInANN[j]);
2458 }
2459
2460
2461 if (gEDBDEBUGLEVEL>3) {
2462 cout << "EdbShowAlgE:: eMLP_SimpleTree " << eMLP_SimpleTree << endl;
2463 cout << "EdbShowAlgE:: eMLP_SimpleTree->ClassName() " << eMLP_SimpleTree->ClassName() << endl;
2464 cout << "EdbShowAlgE:: eMLP_SimpleTree->GetEntries():" << endl;
2465 cout << eMLP_SimpleTree->GetEntries() <<endl;
2466 }
2467
2468
2469 if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlgE:: eMLP_SimpleTree->Fill() do:" << endl;
2470// gDebug=2;
2471 eMLP_SimpleTree->Fill();
2472// gDebug=1;
2473 if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlgE:: eMLP_SimpleTree->Fill() done." << endl;
2474
2475 return;
2476}
EdbShowerP::N
Int_t N() const
Definition: EdbShowerP.h:412
EdbShowerP::GetPara_YY
Para_YY GetPara_YY() const
Definition: EdbShowerP.h:566
EdbShowerP::PrintParametrisation_YY
void PrintParametrisation_YY()
Definition: EdbShowerP.cxx:2533
EdbShowerP::PrintBasics
void PrintBasics()
Definition: EdbShowerP.cxx:2360
for
for(int i=0;i< nentries;i++)
Definition: check_shower.C:42
EdbShowerP::Para_YY
Definition: EdbShowerP.h:215
EdbShowerP::Para_YY::ProfileLongitudinalBincontent
Int_t ProfileLongitudinalBincontent[57]
Definition: EdbShowerP.h:219
EdbShowerP::Para_YY::ProfileTransversalBincontent
Int_t ProfileTransversalBincontent[8]
Definition: EdbShowerP.h:220
EdbShowerP::Para_YY::nseg
Int_t nseg
Definition: EdbShowerP.h:218

◆ Finalize()

void EdbShowAlgE::Finalize ( )
virtual
256{
257 return;
258}

◆ GetHistInANN()

TH1F * EdbShowAlgE::GetHistInANN ( Int_t  nr) const
inline
218 {
219 return eInANN_Hist[nr];
220 }

◆ GetHistoCanvas()

TCanvas * EdbShowAlgE::GetHistoCanvas ( )
1162{
1163 TCanvas* canv= new TCanvas();
1164
1165 eMLP_InputNeurons=7;
1166 Int_t n=TMath::Sqrt(eMLP_InputNeurons)+1;
1167 //cout << "n= " << n << endl;
1168 canv->Divide(n,n);
1169 for (Int_t j=1; j<=eMLP_InputNeurons+1; j++) {
1170 canv->cd(j);
1171 eInANN_Hist[j-1]->Draw();
1172 }
1173
1174 return canv;
1175
1176}
EdbShowAlgE::eMLP_InputNeurons
Int_t eMLP_InputNeurons
Definition: EdbShowAlgE.h:113
TCanvas
new TCanvas()
canv
TCanvas * canv
Definition: testBGReduction_By_ANN.C:43

◆ GetLastBinHistoFilles()

int EdbShowAlgE::GetLastBinHistoFilles ( TH1 *  h)
private
2491{
2492 int nbin=h1->GetNbinsX(); //overflow=nbin
2493 for (int i=nbin-1; i>0; i--) if (h1->GetBinContent(i)>0) return i;
2494 return 0;
2495}
h1
TH1F * h1
Definition: energy.C:16

◆ GetMLP()

TMultiLayerPerceptron * EdbShowAlgE::GetMLP ( ) const
inline
214 {
215 return eMLP_Simple;
216 }

◆ GetMLPArrayNr()

Int_t EdbShowAlgE::GetMLPArrayNr ( Int_t  NPlToCheck)
398{
399 cout <<"dbShowAlgE::GetMLPArray"<<endl;
400
401 Int_t n=-1;
402 for (Int_t i=0; i<10; ++i) {
403 //cout << "i = eShowAlgEPlateBinning->At(i) = NPlToCheck " << i << " " << eShowAlgEPlateBinning->At(i) << " " << NPlToCheck << endl;
404 n=i;
405 if ( eShowAlgEPlateBinning->At(i) >= NPlToCheck) return n;
406 }
407
408 // In standard case it retuns 0 for the first MLP
409 return 0;
410}

◆ GetMLPInputFileName()

TString EdbShowAlgE::GetMLPInputFileName ( )
inline
228 {
229 return eMLP_InputFileName;
230 }
EdbShowAlgE::eMLP_InputFileName
TString eMLP_InputFileName
Definition: EdbShowAlgE.h:117

◆ GetMLPLayout()

TString EdbShowAlgE::GetMLPLayout ( )
inline
222 {
223 return eMLP_Layout;
224 }

◆ GetMLPWeightFileName()

TString EdbShowAlgE::GetMLPWeightFileName ( )
inline
225 {
226 return eMLP_WeightFileName;
227 }

◆ Init() [1/2]

void EdbShowAlgE::Init ( void  )
private

eInANN = new Float_t[70]; eInfoANN = new Float_t[70];

164{
165 cout << "EdbShowAlgE::Init()"<<endl;
166 // default ShowAlge name is set to "XX" (cause it is implemented now.)
167 // later its gonna be "FJ".. or the one parametrisation, which performs best.
168 eShowAlgEParaName="XX";
169 eShowAlgEParaNr=4;
170
171 // Create the Arrays
174
175
176 // Init to set Paraname and ParaNr
177 Init(eShowAlgEParaName);
178
179 return;
180}

◆ Init() [2/2]

void EdbShowAlgE::Init ( TString  ParaName)
private
185{
186
187 cout << "EdbShowAlgE::Init("<<ParaName<<")"<<endl;
188
189 if (ParaName=="FJ") eShowAlgEParaNr=0;
190 if (ParaName=="LT") eShowAlgEParaNr=1;
191 if (ParaName=="YC") eShowAlgEParaNr=2;
192 if (ParaName=="JC") eShowAlgEParaNr=3;
193 if (ParaName=="XX") eShowAlgEParaNr=4;
194 if (ParaName=="YY") eShowAlgEParaNr=5;
195 if (ParaName=="PP") eShowAlgEParaNr=6;
196 if (ParaName=="AS") eShowAlgEParaNr=7;
197
198 eShowAlgEParaName=ParaName;
199
200 cout << "EdbShowAlgE::Init("<<eShowAlgEParaNr<<")"<<endl;
201
202
203 // Now be careful !!
204 // These functions have to be called in THE RIGHT ORDER to have the necessary variables correctly filled:
205
206 // Read Parameter Values for settings to be taken into account:
207 ReadShowAlgEPar("default.par");
208
209
210 // If by "default.par" still no weightfiles and therefore no eMLPArrays are given,
211 // create them by the standard way:
212 if (eMLPArrayN==0) {
213 CreateStandardWeightFileStrings();
214 CreateStandardMLPArrayN();
215 }
216 if (eMLP_InputFileName=="EMPTY") {
217 // jaja ich weiss, super stil...
218 // Set Input file name (default)
219 TString tmpString="MLPInput_"+TString(Form("%d_",eShowAlgEPlatesN))+ParaName+".root";
220 SetMLPInputFileName(tmpString);
221 }
222
223 // Create (on heap) input histograms:
224 CreateInputHistos();
225 // Specific bins for Parametrization:
226 SetBinsInputHistos();
227
228 // Make Layouts (default);
229 CreateMLPLayouts();
230
231 // Make MLPs (default);
232 CreateMLPs();
233
234 return;
235}
EdbShowAlgE::CreateInputHistos
void CreateInputHistos()
Definition: EdbShowAlgE.cxx:1145
EdbShowAlgE::ReadShowAlgEPar
int ReadShowAlgEPar(const char *file)
DEBUG ---------------------------— TEST DEVOPLEMENT.
Definition: EdbShowAlgE.cxx:802
EdbShowAlgE::CreateMLPLayouts
void CreateMLPLayouts()
Definition: EdbShowAlgE.cxx:1599
EdbShowAlgE::CreateMLPs
void CreateMLPs()
Definition: EdbShowAlgE.cxx:2093
EdbShowAlgE::eShowAlgEPlatesN
Int_t eShowAlgEPlatesN
Definition: EdbShowAlgE.h:49
EdbShowAlgE::SetBinsInputHistos
void SetBinsInputHistos()
Definition: EdbShowAlgE.cxx:2185
EdbShowAlgE::CreateStandardWeightFileStrings
void CreateStandardWeightFileStrings()
Definition: EdbShowAlgE.cxx:1402
EdbShowAlgE::CreateStandardMLPArrayN
void CreateStandardMLPArrayN()
Definition: EdbShowAlgE.cxx:1374
EdbShowAlgE::SetMLPInputFileName
void SetMLPInputFileName(TString MLP_InputFileName)
Definition: EdbShowAlgE.h:181

◆ Initialize()

void EdbShowAlgE::Initialize ( )
virtual
243{
244 return;
245}

◆ LoadMLPWeights() [1/2]

void EdbShowAlgE::LoadMLPWeights ( )
2158{
2159 // Since the layout if sixed by the parametrisation, this function
2160 // can be equall for all types of parametrisations:
2161
2162 cout << "void EdbShowAlgE::LoadMLPWeights() " << endl;
2163 if (!eMLPArray) {
2164 cout << "void EdbShowAlgE::LoadMLPWeights() ERROR: eMLPArray IS EMPTY. RETURN!" << endl;
2165 return;
2166 }
2167 for (Int_t i=0; i<eMLPArrayN; ++i) LoadMLPWeights(i);
2168
2169 return;
2170}
EdbShowAlgE::LoadMLPWeights
void LoadMLPWeights()
Definition: EdbShowAlgE.cxx:2157

◆ LoadMLPWeights() [2/2]

void EdbShowAlgE::LoadMLPWeights ( Int_t  nr)
2175{
2176 TMultiLayerPerceptron* mlp;
2177 mlp = (TMultiLayerPerceptron*)eMLPArray->At(nr);
2178 cout << "void EdbShowAlgE::LoadMLPWeights() Name ("<< nr <<"): "<< eMLP_WeightFileNameArray[nr] << endl;
2179 mlp->LoadWeights( eMLP_WeightFileNameArray[nr]);
2180 return;
2181}
mlp
TMultiLayerPerceptron * mlp
Definition: testBGReduction_By_ANN.C:61

◆ Print()

void EdbShowAlgE::Print ( )
263{
264 cout << "EdbShowAlgE::Print()" << endl;
265 cout << "EdbShowAlgE::Print()...done." << endl;
266 return;
267}

◆ PrintWeightFileStrings()

void EdbShowAlgE::PrintWeightFileStrings ( )
1561{
1562 if (!eShowAlgEPlateBinning) return;
1563 if (eMLPArrayN<1) return;
1564
1565 cout << "void EdbShowAlgE::PrintWeightFileStrings() " << endl;
1566
1567 // Create Standard weight file strings:
1568 // weights_energy_"i".txt, i=11,16,21,26,31,36,41,46
1569 for (Int_t i=0; i<eMLPArrayN; ++i) {
1570 cout << "void EdbShowAlgE::PrintWeightFileStrings() Name ("<<i<<"): "<< eMLP_WeightFileNameArray[i] << endl;
1571 }
1572
1573 return;
1574}

◆ ReadMLPInputFile()

void EdbShowAlgE::ReadMLPInputFile ( )
528{
529 if (!eMLP_ReadMLPInputFile) return;
530 cout << eMLP_InputFileName << endl;
531 TFile* file= new TFile(eMLP_InputFileName,"READ");
532 eMLP_SimpleTree=(TTree*)file->Get("eMLP_SimpleTree");
533 eMLP_SimpleTree->Print();
534 cout <<"void EdbShowAlgE::ReadMLPInputFile()...done."<<endl;
535 return;
536}
file
TFile * file
Definition: write_pvr.C:3

◆ ReadShowAlgEPar()

int EdbShowAlgE::ReadShowAlgEPar ( const char *  file = "default.par")

DEBUG ---------------------------— TEST DEVOPLEMENT.

TO IMPLEMENT TAHTS ITS REALLY WRITTEN !!!! sakjdhskjadhlaskjdhlasdlsakdhaslkd

TO IMPLEMENT TAHTS ITS REALLY WRITTEN !!!! sakjdhskjadhlaskjdhlasdlsakdhaslkd

803{
804 // read parameters from par-file
805 // return: 0 if ok
806 // -1 if file access failed
807 char buf[256];
808 char key[256];
809 char name[256];
810
811 FILE *fp = fopen(file,"r");
812 if (!fp) {
813 Log(1,"ReadShowAlgEPar","ERROR open file: %s", file);
814 return -1;
815 }
816 else Log(2,"ReadShowAlgEPar","Read shower parameters from file: %s", file );
817
818 int dospec;
819 int narg;
820 TString tmpString;
821 char tmpchar[256];
822 TString tmpString2;
823 char tmpchar2[256];
824
825 // Since ShowerParametrisation if Fixed by default in the Constructor,
826 // we can check here already if in a String the correct parametrisation is contained.
827
828 while (fgets(buf, sizeof(buf), fp)) {
829 for (Int_t i = 0; i < (Int_t)strlen(buf); i++)
830 if (buf[i]=='#') {
831 buf[i]='\0'; // cut out comments starting from #
832 break;
833 }
834
835 if ( sscanf(buf,"%s",key)!=1 ) continue;
836
837 if ( !strcmp(key,"INCLUDE") )
838 {
839 sscanf(buf+strlen(key),"%s",name);
840 ReadShowAlgEPar(name);
841 }
842
843 else if ( !strcmp(key,"SHOW_ENERGY_TRAIN") )
844 {
845 narg=sscanf(buf+strlen(key),"%d",&dospec);
846 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_TRAIN = " << dospec << " narg= " << narg << endl;
847 if (dospec==0 || dospec==1) SetDoTrain(dospec);
848 }
849 else if ( !strcmp(key,"SHOW_ENERGY_RUN") )
850 {
851 sscanf(buf+strlen(key),"%d",&dospec);
852 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_RUN = " << dospec << endl;
853 if (dospec==0 || dospec==1) SetDoRun(dospec);
854 }
855
856 else if ( !strcmp(key,"SHOW_ENERGY_PLATES") )
857 {
858 sscanf(buf+strlen(key),"%d",&dospec);
859 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_PLATES = " << dospec << endl;
860 SetPlatesN(dospec);
861 }
862
863
864
865
866
867 else if ( !strcmp(key,"SHOW_ENERGY_FILE_WEIGHTFILE") )
868 {
869 sscanf(buf+strlen(key),"%s %s",tmpchar2, tmpchar);
870 tmpString=TString(tmpchar);
871 tmpString2=TString(tmpchar2);
872 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_FILE_WEIGHTFILE = " << tmpString << endl;
873 if (!CheckStringToShowAlgEParaName(tmpString2)) continue;
874// if (!CheckStringToShowAlgEParaName(tmpString2)) break;
875
876 SetMLPWeightFileName(tmpString);
877 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_FILE_WEIGHTFILE = " << eMLP_WeightFileName << endl;
878 }
879
880
881 else if ( !strcmp(key,"SHOW_ENERGY_WRITE_WEIGHTFILE") )
882 {
883 sscanf(buf+strlen(key),"%s %d",tmpchar2, &dospec);
884 tmpString2=TString(tmpchar2);
885 if (!CheckStringToShowAlgEParaName(tmpString2)) continue;
886 SetWriteMLPWeightFile(dospec);
887 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_WRITE_WEIGHTFILE = " << dospec << endl;
888 }
889
890
893
894
895
896
897 else if ( !strcmp(key,"SHOW_ENERGY_FILE_MLPINPUTFILE") )
898 {
899 sscanf(buf+strlen(key),"%s %s",tmpchar2, tmpchar);
900 tmpString=TString(tmpchar);
901 tmpString2=TString(tmpchar2);
902 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_FILE_MLPINPUTFILE = " << tmpString << endl;
903 cout << CheckStringToShowAlgEParaName(tmpString2) << endl;
904 if (!CheckStringToShowAlgEParaName(tmpString2)) continue;
905// if (!CheckStringToShowAlgEParaName(tmpString2)) break;
906 SetMLPInputFileName(tmpString);
907 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_FILE_MLPINPUTFILE = " << eMLP_InputFileName << endl;
908 }
909
910
911 else if ( !strcmp(key,"SHOW_ENERGY_WRITE_MLPINPUTFILE") )
912 {
913 sscanf(buf+strlen(key),"%s %d",tmpchar2, &dospec);
914 tmpString2=TString(tmpchar2);
915 if (!CheckStringToShowAlgEParaName(tmpString2)) continue;
916 SetWriteMLPInputFile(dospec);
917 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_WRITE_MLPINPUTFILE = " << dospec << endl;
918 }
919
920
921 else if ( !strcmp(key,"SHOW_ENERGY_READ_MLPINPUTFILE") )
922 {
923 sscanf(buf+strlen(key),"%s %d",tmpchar2, &dospec);
924 tmpString2=TString(tmpchar2);
925 if (!CheckStringToShowAlgEParaName(tmpString2)) continue;
926 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_READ_MLPINPUTFILE = " << dospec << endl;
927 SetReadMLPInputFile(dospec);
928 }
929
930
933
934
935
936
937
938
939 else if ( !strcmp(key,"SHOW_ENERGY_TRAINEPOCHS") )
940 {
941 sscanf(buf+strlen(key),"%d",&dospec);
942 SetMLPTrainEpochs(dospec);
943 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_TRAINEPOCHS = " << dospec << endl;
944 }
945
946
947
948 else if ( !strcmp(key,"SHOW_ENERGY_ADD_WEIGHTFILE_PLATE") )
949 {
950 sscanf(buf+strlen(key),"%s %d",tmpchar, &dospec);
951 tmpString=TString(tmpchar);
952 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_ADD_WEIGHTFILE_PLATE = " << tmpString << " " << dospec << endl;
953 AddWeightFilePlate(tmpString,dospec);
954 }
955
956 else if ( !strcmp(key,"SHOW_ENERGY_READ_WEIGHTFILE_PLATES") )
957 {
958 int nspec;
959 int spec[10];
960 nspec=sscanf(buf+strlen(key),"%s %d %d %d %d %d %d %d %d %d %d " ,tmpchar, &spec[0], &spec[1], &spec[2], &spec[3], &spec[4], &spec[5], &spec[6], &spec[7], &spec[8], &spec[9]);
961 cout << "nspec= " << nspec << endl;
962 tmpString=TString(tmpchar);
963 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_ADD_WEIGHTFILE_PLATES = " << tmpString;
964 for (int i=0; i<nspec-1; i++ ) {
965 cout <<" " << spec[i];
966 }
967 for (int i=0; i<nspec-1; i++ ) AddWeightFilePlate(tmpString,spec[i]);
968 }
969
970 else if ( !strcmp(key,"SHOW_ENERGY_ADD_WEIGHTFILE") )
971 {
972 sscanf(buf+strlen(key),"%s %d %s",tmpchar, &dospec, tmpchar2);
973 tmpString=TString(tmpchar);
974 tmpString2=TString(tmpchar2);
975 cout << "EdbShowRec::ReadShowAlgEPar SHOW_ENERGY_ADD_WEIGHTFILE_PLATE = " << tmpString << " " << dospec << " " << tmpString2 <<endl;
976 AddWeightFile(tmpString,dospec,tmpString2);
977 }
978
979
980
981 }
982 fclose(fp);
983
984
985 PrintWeightFileStrings();
986 cout <<"===============" << GetMLPInputFileName() << endl;
987
988 Log(2,"ReadShowAlgEPar","Read shower parameters...done.");
989 return 0;
990}
EdbShowAlgE::SetMLPTrainEpochs
Int_t SetMLPTrainEpochs()
Definition: EdbShowAlgE.h:231
EdbShowAlgE::SetWriteMLPWeightFile
void SetWriteMLPWeightFile(Int_t type)
Definition: EdbShowAlgE.cxx:1015
EdbShowAlgE::PrintWeightFileStrings
void PrintWeightFileStrings()
Definition: EdbShowAlgE.cxx:1560
EdbShowAlgE::SetDoTrain
void SetDoTrain(Int_t type)
Definition: EdbShowAlgE.cxx:995
EdbShowAlgE::SetWriteMLPInputFile
void SetWriteMLPInputFile(Int_t type)
Definition: EdbShowAlgE.cxx:1026
EdbShowAlgE::AddWeightFile
void AddWeightFile(TString tmpString, Int_t Plate, TString tmpString2)
Definition: EdbShowAlgE.cxx:1498
EdbShowAlgE::SetDoRun
void SetDoRun(Int_t type)
Definition: EdbShowAlgE.cxx:1004
EdbShowAlgE::SetReadMLPInputFile
void SetReadMLPInputFile(Int_t type)
Definition: EdbShowAlgE.cxx:1036
EdbShowAlgE::SetPlatesN
void SetPlatesN(Int_t ShowAlgEPlatesN)
Definition: EdbShowAlgE.h:191
EdbShowAlgE::CheckStringToShowAlgEParaName
Bool_t CheckStringToShowAlgEParaName(TString tmpString2)
Definition: EdbShowAlgE.cxx:2480
EdbShowAlgE::AddWeightFilePlate
void AddWeightFilePlate(TString tmpString, Int_t dospec)
Definition: EdbShowAlgE.cxx:1437
fclose
fclose(pFile)
name
const char * name
Definition: merge_Energy_SytematicSources_Electron.C:24

◆ ResetPlateBinning()

void EdbShowAlgE::ResetPlateBinning ( )
1317{
1318 Log(2,"EdbShowAlgE::ResetPlateBinning","EdbShowAlgE::ResetPlateBinning.");
1319
1320
1321 // Restore Standard Plate Binning:
1322 // (11,16,21..51,56)
1323 if (!eShowAlgEPlateBinning) eShowAlgEPlateBinning = new TArrayI(10);
1324
1325 cout << "void EdbShowAlgE::ResetPlateBinning(): Steps for:";
1326 for (Int_t i=0; i<10; ++i) {
1327 eShowAlgEPlateBinning->AddAt(11+5*i,i);
1328 cout << eShowAlgEPlateBinning->At(i);
1329 }
1330 cout << " plates availible."<<endl;
1331 Log(2,"EdbShowAlgE::ResetPlateBinning","EdbShowAlgE::ResetPlateBinning...done.");
1332 return;
1333}

◆ ResetWeightFileNameArray()

void EdbShowAlgE::ResetWeightFileNameArray ( )
private
1364{
1365 for (Int_t i=0; i<10; ++i) eMLP_WeightFileNameArray[i]="";
1366 return;
1367}

◆ Run()

void EdbShowAlgE::Run ( )

Load Weights : Absolutely neccessary! MISSING: A CHECK IF THE WEIGHTFILE IS IN THE CURRENT DIRECTORY AND THERE AT ALL !!! (otherwise link from source dir)

Weights have to be already loaded ??????

286{
287 cout << "void EdbShowAlgE::Run() "<< endl;
288
289 if (!eMLP_SimpleTree) {
290 cout << "void EdbShowAlgE::Run() NO eMLP_SimpleTree Read eMLP_SimpleTree from file:";
291 cout << "void EdbShowAlgE::Run() "<< eMLP_InputFileName << endl;
292 if (eMLP_ReadMLPInputFile) ReadMLPInputFile();
293 }
294
295 if (!eMLP_SimpleTree) {
296 cout << "void EdbShowAlgE::Run() NO eMLP_SimpleTree Some Error has occured! RETURN."<< endl;
297 return;
298 }
299
300 // Set Branch Addresses for the tree:
301 eMLP_SimpleTree->SetBranchAddress("eInANN",eInANN);
302 eMLP_SimpleTree->SetBranchAddress("eInfoANN",eInfoANN);
303
304
305 PrintWeightFileStrings(); // DEBUG()
306
307
308
309 //------------------------------
313 LoadMLPWeights();
314 //------------------------------
315
316
317 // Set TTree as source for the MLP:
318 // Absolutely neccessary otherwise its empty.
319 for (Int_t i=0; i<10; ++i) {
320 TEventList* list;
321 list= (TEventList*)eMLP_SimpleTree->GetEventList();
322 eMLP_Simple= (TMultiLayerPerceptron*)eMLPArray->At(i);
323 cout << eMLP_Simple->GetStructure() << endl;
324 cout << eMLP_SimpleTree->GetEntries() << endl;
325 eMLP_Simple->Print();
326 eMLP_Simple->SetData(eMLP_SimpleTree);
327 eMLP_Simple->SetTrainingDataSet(list);
328 eMLP_Simple->SetTestDataSet(list);
329 eMLP_Simple->SetTrainingDataSet("Entry$%2");
330 eMLP_Simple->SetTestDataSet("(Entry$+1)%2");
331 }
332
333
334 Int_t nMLPArray=-1;
335 Double_t params[70]; // MLP works only with Double_t
336 Float_t val;
337
338
339 TH1F* h1= new TH1F("h1","h1",100,0,10);
340 // Loop over entries:
341// for (Int_t i=1; i<eMLP_SimpleTree->GetEntries(); ++i) {
342 for (Int_t i=0; i<3; ++i) {
343 eMLP_SimpleTree->GetEntry(i);
344
345 eMLP_SimpleTree->Show(i);
346 cout << "i= " << i << " " << eInfoANN[0] << " " << eInfoANN[1] << endl;
347 cout << "i= " << i << " " << eInANN[0] << " " << eInANN[1] << endl;
348
349 // Get Right ANN (depending on reconstructed plates).
351 nMLPArray=GetMLPArrayNr(eInfoANN[1]);
352
353 cout << nMLPArray << endl;
354
355 eMLP_Simple= (TMultiLayerPerceptron*)eMLPArray->At(nMLPArray);
356 cout << eMLP_Simple->GetStructure() << endl;
357
358
359
360
361
362
363
364 for (Int_t j=1; j<=70; ++j) {
365// cout << "DEBUG eInANN[j] " << eInANN[j] << endl;
366 params[j-1]=eInANN[j];
367 }
368 val=(eMLP_Simple->Evaluate(0,params));
369
370
371 h1->Fill(val);
372
373 cout << "DEBUG " << eInANN[1] << " " << eInANN[2] << " " << eInANN[3] << " " << eInANN[4] << " " << eInANN[5] << " " << eInANN[6] << " " << eInANN[7] << " " << endl;
374 cout << "DEBUG: i, val (=mlp1->Evaluate(0,params)) == real energy " << i << " " << val << " == " << eInANN[0] << endl;
375
376
377 }
378
379
380
381 h1->Draw();
382
383
384// eMLP_Simple= (TMultiLayerPerceptron*)eMLPArray->At(0);
385
386// cout << eMLP_Simple->GetStructure() << endl;
387// cout << eMLP_SimpleTree->GetEntries() << endl;
388
389
390 return;
391}
EdbShowAlgE::GetMLPArrayNr
int GetMLPArrayNr(Int_t NPlToCheck)
Definition: EdbShowAlgE.cxx:397
params
Double_t params[3]
Definition: testBGReduction_By_ANN.C:84

◆ Set0()

void EdbShowAlgE::Set0 ( )
private
101{
102 // Set0()
103 eRecoShowerArray=0;
104 eRecoShowerArrayN=0;
105
106 eDoANNTrain=kFALSE;
107 eDoANNRun=kFALSE;
108
109 cout << " Test if Neural Network is loaded: " << endl;
110 if (!gROOT->GetClass("TMultiLayerPerceptron")) cout << "NOT"<<endl;
111 if (gROOT->GetClass("TMultiLayerPerceptron")) cout << "YES"<<endl;
112 cout << " So now we know if the TMultiLayerPerceptron class is loaded... " << endl;
113
114 eMLP_Simple=0;
115 eMLP_Layout="";
116 eMLP_LayoutArray=0;
117 eMLP_LayoutArrayIsDone=kFALSE;
118 eMLP_WeightFileName="";
119 eMLP_SimpleTree=0;
120 eMLP_TrainEpochs=100;
121
122
123 // Reset All eInANN Values to Zero:
124
125 for (Int_t j=0; j<70; j++) {
126 eInANN[j]=0;
127 eInfoANN[j]=0;
128 eInANN_Hist[j]=0;
129 }
130 eOutANN=0;
131
132 // WHEN COMPILING THE LIBRARY, AN ERROR APPEARS HERE:
133// Error: *** Datamember EdbShowAlgE::eInANN: no size indication!
134// Error: *** Datamember EdbShowAlgE::eInANN: pointer to fundamental type (need manual intervention)
135// Error: *** Datamember EdbShowAlgE::eInfoANN: no size indication!
136// Error: *** Datamember EdbShowAlgE::eInfoANN: pointer to fundamental type (need manual intervention)
137
138
139
140 eMLP_InputNeurons=0;
141 eMLP_InputNeuronsArray=0;
142 eMLP_InputNeuronsArrayN=0;
143
144 eMLP_ReadMLPInputFile=kFALSE;
145 eMLP_WriteMLPInputFile=kFALSE;
146 eMLP_InputFileName="EMPTY";
147
148 eMLP_WriteWeightFile=kFALSE;
149 eMLP_WeightFileName="EMPTY";
150
151 eMLPArrayN=0;
152 eMLPArray=0;
153
154 eShowAlgEPlateBinning=0;
155
156 eShowAlgEPlatesN=0;
157
158 return;
159}
EdbShowAlgE::eMLP_TrainEpochs
Int_t eMLP_TrainEpochs
Definition: EdbShowAlgE.h:111
EdbShowAlgE::eOutANN
Float_t eOutANN
Definition: EdbShowAlgE.h:87
EdbShowAlgE::eDoANNRun
Bool_t eDoANNRun
Definition: EdbShowAlgE.h:68
EdbShowAlgE::eMLP_WriteMLPInputFile
Bool_t eMLP_WriteMLPInputFile
Definition: EdbShowAlgE.h:118
EdbShowAlgE::eDoANNTrain
Bool_t eDoANNTrain
Definition: EdbShowAlgE.h:67

◆ SetBinsInputHistos()

void EdbShowAlgE::SetBinsInputHistos ( )
2185 {
2186
2187 switch (eShowAlgEParaNr) {
2188 case 0 :
2189 SetBinsInputHistos_FJ();
2190 break;
2191 case 1 : //CreateInputNeuronsArray_LT();
2192 break;
2193 case 2 :
2194 SetBinsInputHistos_YY();
2195 break;
2196 case 3 : //CreateInputNeuronsArray_JC();
2197 break;
2198 case 4 :
2199 SetBinsInputHistos_XX();
2200 break;
2201 case 5 :
2202 SetBinsInputHistos_YY();
2203 break;
2204 default :
2205 SetBinsInputHistos_XX();
2206 break;
2207 }
2208 return;
2209}
EdbShowAlgE::SetBinsInputHistos_FJ
void SetBinsInputHistos_FJ()
Definition: EdbShowAlgE.cxx:2216
EdbShowAlgE::SetBinsInputHistos_YY
void SetBinsInputHistos_YY()
Definition: EdbShowAlgE.cxx:2254
EdbShowAlgE::SetBinsInputHistos_XX
void SetBinsInputHistos_XX()
Definition: EdbShowAlgE.cxx:2239

◆ SetBinsInputHistos_FJ()

void EdbShowAlgE::SetBinsInputHistos_FJ ( )
2216 {
2217 // Specific bins for FJ Parametrization:
2218 return;
2219}

◆ SetBinsInputHistos_XX()

void EdbShowAlgE::SetBinsInputHistos_XX ( )
2239 {
2240 // Specific bins for XX Parametrization:
2241 SetInputHistos(eInANN_Hist[0],100,0,20000);
2242 SetInputHistos(eInANN_Hist[1],200,0,1000);
2243 SetInputHistos(eInANN_Hist[2],100,0,70000);
2244 SetInputHistos(eInANN_Hist[3],100,0,30000);
2245 SetInputHistos(eInANN_Hist[4],30,0,30);
2246 SetInputHistos(eInANN_Hist[5],100,0,5000);
2247 SetInputHistos(eInANN_Hist[6],100,0,5000);
2248 SetInputHistos(eInANN_Hist[7],30,0,30);
2249 return;
2250}
EdbShowAlgE::SetInputHistos
void SetInputHistos(TH1F *h, Int_t nbin, Double_t lower, Double_t upper)
Definition: EdbShowAlgE.cxx:1180

◆ SetBinsInputHistos_YC()

void EdbShowAlgE::SetBinsInputHistos_YC ( )
2225 {
2226 // Specific bins for YC Parametrization:
2227 SetInputHistos(eInANN_Hist[0],100,0,20000);
2228 SetInputHistos(eInANN_Hist[1],200,0,1000);
2229 SetInputHistos(eInANN_Hist[2],100,0,100);
2230 SetInputHistos(eInANN_Hist[3],100,0,100);
2231 SetInputHistos(eInANN_Hist[4],100,0,100);
2232 SetInputHistos(eInANN_Hist[5],100,0,100);
2233 return;
2234}

◆ SetBinsInputHistos_YY()

void EdbShowAlgE::SetBinsInputHistos_YY ( )
2254 {
2255 // Specific bins for YY Parametrization:
2256 SetInputHistos(eInANN_Hist[0],100,0,20000);
2257 SetInputHistos(eInANN_Hist[1],200,0,1000);
2258 SetInputHistos(eInANN_Hist[2],100,0,100);
2259 SetInputHistos(eInANN_Hist[3],100,0,100);
2260 SetInputHistos(eInANN_Hist[4],100,0,100);
2261 SetInputHistos(eInANN_Hist[5],100,0,100);
2262 SetInputHistos(eInANN_Hist[6],100,0,100);
2263 for (int i=0; i<56; i++) SetInputHistos(eInANN_Hist[6+i],100,0,100);
2264 return;
2265}

◆ SetDoRun()

void EdbShowAlgE::SetDoRun ( Int_t  type)
1005{
1006 if (type==1) eDoANNRun=kTRUE;
1007 else eDoANNRun=kFALSE;
1008 return;
1009}
type
Int_t type
Definition: testBGReduction_By_ANN.C:15

◆ SetDoTrain()

void EdbShowAlgE::SetDoTrain ( Int_t  type)
996{
997 if (type==1) eDoANNTrain=kTRUE;
998 else eDoANNTrain=kFALSE;
999 return;
1000}

◆ SetInputHistos()

void EdbShowAlgE::SetInputHistos ( TH1F *  h,
Int_t  nbin,
Double_t  lower,
Double_t  upper 
)
1181{
1182 h->SetBins(nbin,lower,upper);
1183 return;
1184}

◆ SetInputHistosPara()

void EdbShowAlgE::SetInputHistosPara ( Int_t  ParaType)

TO IMPLEMENT ALSO THE OTHERS....

1191{
1192 if (ParaType==0) SetInputHistosPara_FJ();
1193
1194 if (ParaType==2) SetInputHistosPara_YC();
1195 if (ParaType==4) SetInputHistosPara_XX();
1196 if (ParaType==5) SetInputHistosPara_YY();
1197
1198
1200
1201 return;
1202}
EdbShowAlgE::SetInputHistosPara_YC
void SetInputHistosPara_YC()
Definition: EdbShowAlgE.cxx:1205
EdbShowAlgE::SetInputHistosPara_YY
void SetInputHistosPara_YY()
Definition: EdbShowAlgE.cxx:1247
EdbShowAlgE::SetInputHistosPara_FJ
void SetInputHistosPara_FJ()
Definition: EdbShowAlgE.cxx:1214

◆ SetInputHistosPara_FJ()

void EdbShowAlgE::SetInputHistosPara_FJ ( )
1215{
1216
1217 return;
1218}

◆ SetInputHistosPara_XX()

void EdbShowAlgE::SetInputHistosPara_XX ( )
1224{
1225
1226 // ---------
1227 // Specific bins for XX Parametrization:
1228 SetInputHistos(eInANN_Hist[0],100,0,20000);
1229
1230 SetInputHistos(eInANN_Hist[1],200,0,1000);
1231
1232 SetInputHistos(eInANN_Hist[2],100,0,70000);
1233 SetInputHistos(eInANN_Hist[3],100,0,30000);
1234 SetInputHistos(eInANN_Hist[4],30,0,30);
1235
1236 SetInputHistos(eInANN_Hist[5],100,0,5000);
1237 SetInputHistos(eInANN_Hist[6],100,0,5000);
1238 SetInputHistos(eInANN_Hist[7],100,0,100);
1239 // ---------
1240
1241
1242 return;
1243}

◆ SetInputHistosPara_YC()

void EdbShowAlgE::SetInputHistosPara_YC ( )
1206{
1207 // ---------
1208
1209 return;
1210}

◆ SetInputHistosPara_YY()

void EdbShowAlgE::SetInputHistosPara_YY ( )
1248{
1249 // ---------
1250 // Specific bins for XX Parametrization:
1251 SetInputHistos(eInANN_Hist[0],100,0,20000);
1252
1253 SetInputHistos(eInANN_Hist[1],200,0,1000);
1254
1255 SetInputHistos(eInANN_Hist[2],100,0,100);
1256 SetInputHistos(eInANN_Hist[3],100,0,100);
1257 SetInputHistos(eInANN_Hist[4],100,0,100);
1258 SetInputHistos(eInANN_Hist[5],100,0,100);
1259 SetInputHistos(eInANN_Hist[6],100,0,100);
1260 for (int i=0; i<56; i++) SetInputHistos(eInANN_Hist[6+i],100,0,100);
1261 // ---------
1262 return;
1263}

◆ SetMLPArrayN()

void EdbShowAlgE::SetMLPArrayN ( Int_t  MLPArrayN)
inline
187 {
188 eMLPArrayN = MLPArrayN;
189 }

◆ SetMLPInputFileName()

void EdbShowAlgE::SetMLPInputFileName ( TString  MLP_InputFileName)
inline
181 {
182 eMLP_InputFileName = MLP_InputFileName;
183 }

◆ SetMLPLayout()

void EdbShowAlgE::SetMLPLayout ( TString  layout)
inline
175 {
176 eMLP_Layout = layout;
177 }

◆ SetMLPTrainEpochs() [1/2]

Int_t EdbShowAlgE::SetMLPTrainEpochs ( )
inline
231 {
232 return eMLP_TrainEpochs;
233 }

◆ SetMLPTrainEpochs() [2/2]

void EdbShowAlgE::SetMLPTrainEpochs ( Int_t  nep)
inline
184 {
185 eMLP_TrainEpochs = nep;
186 }

◆ SetMLPWeightFileName()

void EdbShowAlgE::SetMLPWeightFileName ( TString  weightfilename)
inline
178 {
179 eMLP_WeightFileName = weightfilename;
180 }

◆ SetPlatesN()

void EdbShowAlgE::SetPlatesN ( Int_t  ShowAlgEPlatesN)
inline
191 {
192 eShowAlgEPlatesN = ShowAlgEPlatesN;
193 }

◆ SetReadMLPInputFile()

void EdbShowAlgE::SetReadMLPInputFile ( Int_t  type)
1037{
1038 if (type==1) eMLP_ReadMLPInputFile=kTRUE;
1039 else eMLP_ReadMLPInputFile=kFALSE;
1040 return;
1041}

◆ SetRecoShowerArray()

void EdbShowAlgE::SetRecoShowerArray ( TObjArray *  RecoShowerArray)
inline
171 {
172 eRecoShowerArray=RecoShowerArray;
173 }

◆ SetWriteMLPInputFile()

void EdbShowAlgE::SetWriteMLPInputFile ( Int_t  type)
1027{
1028 if (type==1) eMLP_WriteMLPInputFile=kTRUE;
1029 else eMLP_WriteMLPInputFile=kFALSE;
1030 return;
1031}

◆ SetWriteMLPWeightFile()

void EdbShowAlgE::SetWriteMLPWeightFile ( Int_t  type)
1016{
1017 if (type==1) eMLP_WriteWeightFile=kTRUE;
1018 else eMLP_WriteWeightFile=kFALSE;
1019 return;
1020}

◆ TESTRUN()

void EdbShowAlgE::TESTRUN ( )

DEBUG DEBUG DEBUG

DEBUG DEBUG DEBUG

421{
422
423 // Fill Tree with existing entries of para of eRecoShowerArray :
424 Fill_MLPTree();
425
426
427
428 FillInputHistosPara_XX();
429
430
431
432
433
434 TString layout;
435 int NrOfANNInputNeurons=7;
436 layout=CreateMLPLayout(NrOfANNInputNeurons);
437
438
439
440
441 SetMLPLayout(layout);
442
443
444 cout <<"__________"<<endl;
445 cout << layout << endl;
446 cout <<"__________"<<endl;
447
448 Double_t params[70]; // MLP works only with Double_t
449 Float_t val;
450
451
452
454
455 cout << " eMLP_Simple = " << eMLP_Simple << endl;
456
457 eMLP_Simple =0;
458
459 cout << " eMLP_Simple = " << eMLP_Simple << endl;
460
461// TMultiLayerPerceptron*
462 eMLP_Simple = new TMultiLayerPerceptron(layout,eMLP_SimpleTree,"Entry$%2","(Entry$+1)%2");
463
464
465
466 cout << " eMLP_Simple = " << eMLP_Simple << endl;
467
468
470
471
472 eMLP_Simple->LoadWeights("weights.txt");
473
474 eMLP_Simple->Print();
475
476
477 TH1F* h1 = new TH1F("hh","hh",100,0,10);
478
479
480 for (Int_t i=1; i<eMLP_SimpleTree->GetEntries(); ++i) {
481 eMLP_SimpleTree->GetEntry(i);
482
483 for (Int_t j=1; j<=NrOfANNInputNeurons; ++j) {
484 params[j-1]=eInANN[j];
485 }
486 val=(eMLP_Simple->Evaluate(0,params));
487
488
489 h1->Fill(val);
490
491 cout << "DEBUG " << eInANN[1] << " " << eInANN[2] << " " << eInANN[3] << " " << eInANN[4] << " " << eInANN[5] << " " << eInANN[6] << " " << eInANN[7] << " " << endl;
492 cout << "DEBUG: i, val (=mlp1->Evaluate(0,params)) == real energy " << i << " " << val << " == " << eInANN[0] << endl;
493 }
494
495
496 TCanvas* cc= new TCanvas();
497 h1->Draw();
498
499 return;
500}
EdbShowAlgE::Fill_MLPTree
void Fill_MLPTree()
Definition: EdbShowAlgE.cxx:1072
EdbShowAlgE::SetMLPLayout
void SetMLPLayout(TString layout)
Definition: EdbShowAlgE.h:175
EdbShowAlgE::CreateMLPLayout
TString CreateMLPLayout(Int_t NrOfANNInputNeurons)
Definition: EdbShowAlgE.cxx:1299
EdbShowAlgE::FillInputHistosPara_XX
void FillInputHistosPara_XX()
Definition: EdbShowAlgE.cxx:1269

◆ Train()

void EdbShowAlgE::Train ( )
691{
692 cout << "void EdbShowAlgE::Train() "<< endl;
693 cout << "void EdbShowAlgE::Train() We have to set the tree data from the ANN since at "<< endl;
694 cout << "void EdbShowAlgE::Train() creation time the eMLP_SimpleTree was empty."<< endl;
695
696 if (!eMLP_SimpleTree) {
697 cout << "void EdbShowAlgE::Train() NO eMLP_SimpleTree Read eMLP_SimpleTree from file:"<< endl;
698 if (eMLP_ReadMLPInputFile) ReadMLPInputFile();
699 }
700 if (!eMLP_SimpleTree) {
701 cout << "void EdbShowAlgE::Train() NO eMLP_SimpleTree Some ERROR has occured! RETURN."<< endl;
702 return;
703 }
704
705
706 // In order not to tell which plates we have to train the ANN, we grep out the information
707 // out of eINfoAnn[1], which contains plates. The maximumBin is used for the layout structure:
708 eMLP_SimpleTree->Draw("eInfoANN[1]>>h1");
709 TH1F *h1 = (TH1F*)gDirectory->Get("h1");
710 h1->Draw();
711
712 int MaxBinCenter = h1->GetBinCenter(h1->GetMaximumBin()) ;
713 int MaxBinEntryCenter=h1->GetBinCenter(GetLastBinHistoFilles(h1));
714 int MLPArrayNr = GetMLPArrayNr(MaxBinEntryCenter);
715
716
717
718 cout << " MaxBinCenter = " << MaxBinCenter << endl;
719 cout << " MLPArrayNr = " << MLPArrayNr << endl;
720 cout << " MaxBinEntryCenter = " << MaxBinEntryCenter << endl;
721
722
723
724 // check MaxBinCenter for Bounds
725 if (MaxBinEntryCenter>57 || MaxBinEntryCenter<1) {
726 cout << "void EdbShowAlgE::Train() Something is wrong with MaxBinEntryCenter. Check your data (especially eInfoANN[1]) RETURN! "<<endl;
727 return;
728 }
729 cout << "void EdbShowAlgE::Train() MLPArrayNr= " << MLPArrayNr << " So use this eMLPArray for Training."<<endl;
730
731
732 TEventList* list;
733 list= (TEventList*)eMLP_SimpleTree->GetEventList();
734 eMLP_Simple= (TMultiLayerPerceptron*)eMLPArray->At(MLPArrayNr);
735
736 cout << eMLP_Simple->GetStructure() << endl;
737 cout << eMLP_SimpleTree->GetEntries() << endl;
738
739// eMLP_SimpleTree->Show(0);
740// eMLP_SimpleTree->Show(eMLP_SimpleTree->GetEntries()-1);
741// cout << eMLP_TrainEpochs << endl;
742
743 eMLP_Simple->Print();
744
745 eMLP_Simple->SetData(eMLP_SimpleTree);
746 eMLP_Simple->SetTrainingDataSet(list);
747 eMLP_Simple->SetTestDataSet(list);
748 eMLP_Simple->SetTrainingDataSet("Entry$%2");
749 eMLP_Simple->SetTestDataSet("(Entry$+1)%2");
750 eMLP_Simple->Train(eMLP_TrainEpochs, "text,graph,update=10");
751
752
753
754 TCanvas* mlpa_canvas = new TCanvas("mlpa_canvas","Network analysis");
755 mlpa_canvas->Divide(1,2);
756 TMLPAnalyzer ana(eMLP_Simple);
757 ana.GatherInformations();
758 ana.CheckNetwork();
759 mlpa_canvas->cd(1);
760 ana.DrawDInputs();
761 mlpa_canvas->cd(2);
762 eMLP_Simple->Draw();
763
764 cout << "void EdbShowAlgE::Train() ...done."<< endl;
765 return;
766}
EdbShowAlgE::GetLastBinHistoFilles
int GetLastBinHistoFilles(TH1 *h)
Definition: EdbShowAlgE.cxx:2490
mlpa_canvas
TCanvas * mlpa_canvas
Definition: testBGReduction_By_ANN.C:66
ana
TMLPAnalyzer ana(mlp)

◆ WriteMLP_SimpleTree()

void EdbShowAlgE::WriteMLP_SimpleTree ( TString  filename)
653{
654 if (!eMLP_WriteMLPInputFile) return;
655// TFile* file= new TFile(filename,"RECREATE");
656 TFile* file= new TFile(eMLP_InputFileName,"RECREATE");
657 eMLP_SimpleTree->Write();
658 file->Close();
659 delete file;
660 return;
661}

◆ WriteMLPInputFile()

void EdbShowAlgE::WriteMLPInputFile ( )
507{
508 cout <<"void EdbShowAlgE::WriteMLPInputFile() " << eMLP_InputFileName << endl;
509
510 if (!eMLP_WriteMLPInputFile) {
511 cout <<"void EdbShowAlgE::WriteMLPInputFile() WARNING::eMLP_WriteMLPInputFile= " << eMLP_WriteMLPInputFile << " RETURN." << endl;
512 return;
513 }
514
515 TFile* file= new TFile(eMLP_InputFileName,"RECREATE");
516 eMLP_SimpleTree->Write();
517 eMLP_SimpleTree->Print();
518 file->Close();
519 delete file;
520 cout <<"void EdbShowAlgE::WriteMLPInputFile()...done."<<endl;
521 return;
522}

Member Data Documentation

◆ eDoANNRun

Bool_t EdbShowAlgE::eDoANNRun
private

◆ eDoANNTrain

Bool_t EdbShowAlgE::eDoANNTrain
private

◆ eInANN

Float_t* EdbShowAlgE::eInANN = new Float_t[70]
private

◆ eInANN_Hist

TH1F* EdbShowAlgE::eInANN_Hist[70]
private

◆ eInfoANN

Float_t* EdbShowAlgE::eInfoANN = new Float_t[70]
private

◆ eMLP_InputFileName

TString EdbShowAlgE::eMLP_InputFileName
private

◆ eMLP_InputNeurons

Int_t EdbShowAlgE::eMLP_InputNeurons
private

◆ eMLP_InputNeuronsArray

TArrayI* EdbShowAlgE::eMLP_InputNeuronsArray
private

◆ eMLP_InputNeuronsArrayN

Int_t EdbShowAlgE::eMLP_InputNeuronsArrayN
private

◆ eMLP_Layout

TString EdbShowAlgE::eMLP_Layout
private

◆ eMLP_LayoutArray

TObjArray* EdbShowAlgE::eMLP_LayoutArray
private

◆ eMLP_LayoutArrayIsDone

Bool_t EdbShowAlgE::eMLP_LayoutArrayIsDone
private

◆ eMLP_ReadMLPInputFile

Bool_t EdbShowAlgE::eMLP_ReadMLPInputFile
private

◆ eMLP_Simple

TMultiLayerPerceptron* EdbShowAlgE::eMLP_Simple
private

◆ eMLP_SimpleIsDone

Bool_t EdbShowAlgE::eMLP_SimpleIsDone
private

◆ eMLP_SimpleTree

TTree* EdbShowAlgE::eMLP_SimpleTree
private

◆ eMLP_TrainEpochs

Int_t EdbShowAlgE::eMLP_TrainEpochs
private

◆ eMLP_WeightFileName

TString EdbShowAlgE::eMLP_WeightFileName
private

◆ eMLP_WeightFileNameArray

TString EdbShowAlgE::eMLP_WeightFileNameArray[10]
private

◆ eMLP_WriteMLPInputFile

Bool_t EdbShowAlgE::eMLP_WriteMLPInputFile
private

◆ eMLP_WriteWeightFile

Bool_t EdbShowAlgE::eMLP_WriteWeightFile
private

◆ eMLPArray

TObjArray* EdbShowAlgE::eMLPArray
private

◆ eMLPArrayN

Int_t EdbShowAlgE::eMLPArrayN
private

◆ eOutANN

Float_t EdbShowAlgE::eOutANN
private

◆ eRecoShowerArray

TObjArray* EdbShowAlgE::eRecoShowerArray
private

◆ eRecoShowerArrayN

Int_t EdbShowAlgE::eRecoShowerArrayN
private

◆ eShowAlgEParaName

TString EdbShowAlgE::eShowAlgEParaName
private

◆ eShowAlgEParaNr

Int_t EdbShowAlgE::eShowAlgEParaNr
private

◆ eShowAlgEPlateBinning

TArrayI* EdbShowAlgE::eShowAlgEPlateBinning
private

◆ eShowAlgEPlatesN

Int_t EdbShowAlgE::eShowAlgEPlatesN
private
The documentation for this class was generated from the following files: