#include <EdbShowerAlg.h>

Inheritance diagram for EdbShowerAlgIDSimple:

Collaboration diagram for EdbShowerAlgIDSimple:

Public Member Functions
	ClassDef (EdbShowerAlgIDSimple, 1)

void	CreateANN ()

Int_t	DetermineIDFlag (Int_t IDType)

Int_t	DetermineIDType (Double_t val, Int_t SpecificationType)

void	DoRun ()

void	DoRun (EdbTrackP *)

void	DoRun (TObjArray *trackarray)

	EdbShowerAlgIDSimple ()

	EdbShowerAlgIDSimple (EdbTrackP *track)
	CONSTRUCTOR USED FOR libShowRec with EdbShowerP class available. More...

	EdbShowerAlgIDSimple (EdbTrackP track, EdbVertex vtx)

	EdbShowerAlgIDSimple (TObjArray *RecoShowerArray)

	EdbShowerAlgIDSimple (TObjArray RecoShowerArray, TObjArray VtxArray)

void	Help ()

void	Init ()

void	InitStrings ()

void	Print ()

void	PrintSpecifications ()

void	ReadTables ()

void	ReadTables_ID ()

void	SetANNIDCutValue (Double_t ANNIDCutValue)

void	SetSpecificationType (Int_t SpecificationType, Int_t SpecificationTypeVal)

void	Update ()

virtual	~EdbShowerAlgIDSimple ()

Public Member Functions inherited from EdbShowerAlgESimple
	ClassDef (EdbShowerAlgESimple, 1)

void	CreateANN ()

void	DoRun ()

void	DoRun (EdbTrackP *shower)

void	DoRun (TObjArray *trackarray)

void	DumpNeuralNetworkWeight (TString weight, Int_t ANNType=0)

	EdbShowerAlgESimple ()

	EdbShowerAlgESimple (EdbTrackP *track)
	CONSTRUCTOR USED FOR libShowRec with EdbShowerP class available. More...

	EdbShowerAlgESimple (TObjArray *RecoShowerArray)

Int_t	FindClosestEfficiencyParametrization (Double_t TestAngle, Double_t ReferenceEff)

TF1 *	GetEffFunc_all ()

TF1 *	GetEffFunc_edefault ()

TF1 *	GetEffFunc_elletroni ()

TF1 *	GetEffFunc_LowEff ()

TF1 *	GetEffFunc_MiddleFix ()

TF1 *	GetEffFunc_neuchmicro ()

TF1 *	GetEffFunc_UserEfficiency ()

TSpline3 *	GetEfficiencyParametrisation ()

Float_t	GetEnergy ()

Float_t	GetEnergy (EdbTrackP *track)

TArrayF *	GetEnergyArray ()

TArrayF *	GetEnergyArrayUnCorrected ()

TMultiLayerPerceptron *	GetNeuralNetwork (Int_t ANNType=0)

void	GetNplIndexNr (Int_t sh_npl, Int_t &check_Npl_index, Int_t ePlateNumberType)

void	GetPara (EdbTrackP *track)

Int_t	GetRecoShowerArrayN () const

void	GetSpecifications ()

Int_t	GetSpecType (Int_t SpecificationType)

void	Help ()

void	InitStrings ()

void	LoadSpecificationWeightFile ()

void	Print ()

void	PrintEfficiencyParametrisation ()

void	PrintSpecifications ()

void	ReadCorrectionFactors (TString weigthstring, Float_t &p0, Float_t &p1)

void	ReadTables ()

void	ReadTables_Energy ()

void	SetCalibrationOffset (Float_t CalibrationOffset)

void	SetCalibrationSlope (Float_t CalibrationSlope)

void	SetEfficiencyParametrisation (TSpline3 *EfficiencyParametrisation)

void	SetEfficiencyParametrisationAngles ()

void	SetEfficiencyParametrisationValues (Double_t Angles, Double_t EffValuesAtAngles)

void	SetForceSpecificationReload ()

void	SetPlateNumber (Int_t PlateNumber)

void	SetPlateNumberType (Int_t PlateNumberType)

void	SetRecoShowerArray (TObjArray *RecoShowerArray)

void	SetRecoShowerArrayN (Int_t RecoShowerArrayN)

void	SetSpecifications (Int_t sp0, Int_t sp1, Int_t sp2, Int_t sp3, Int_t sp4, Int_t sp5)

void	SetSpecificationType (Int_t SpecificationType, Int_t SpecificationTypeVal)

void	SetWeightFileString (TString weightstring)

void	TrainNeuralNetwork (TString weight, Int_t ANNType=0)

void	UnSetForceSpecificationReload ()

void	Update ()

void	WriteNewRootFile ()

void	WriteNewRootFile (TString sourcefilename)

void	WriteNewRootFile (TString sourcefilename, TString treename)

virtual	~EdbShowerAlgESimple ()

Protected Attributes
Double_t	eANNIDCutValue

Double_t	eIDCutTypeValue [3]

Double_t	eIDEnergySigmaCorr

Protected Attributes inherited from EdbShowerAlgESimple
TString	ANN_Layout

TMultiLayerPerceptron *	ANN_MLP

TMultiLayerPerceptron *	ANN_MLP_ARRAY [15]

Int_t	ANN_n_InputNeurons

Int_t	ANN_n_InputNeurons_ARRAY [15]

Int_t	ANN_nPlates_ARRAY [15]

TString	ANN_WeightFile_ARRAY [15]

TTree *	ANNTree

Float_t	eANN_MLP_CORR_0 [15]

Float_t	eANN_MLP_CORR_1 [15]

Float_t	eCalibrationOffset

Float_t	eCalibrationSlope

TSpline3 *	eEfficiencyParametrisation

Float_t	eEnergy

TArrayF *	eEnergyArray

Int_t	eEnergyArrayCount

TArrayF *	eEnergyArraySigmaCorrected

TArrayF *	eEnergyArrayUnCorrected

Float_t	eEnergyCorr

Float_t	eEnergySigmaCorr

Float_t	eEnergyUnCorr

TF1 *	EffFunc_all

TF1 *	EffFunc_edefault

TF1 *	EffFunc_elletroni

TF1 *	EffFunc_LowEff

TF1 *	EffFunc_MiddleFix

TF1 *	EffFunc_neuchmicro

TF1 *	EffFunc_UserEfficiency

Bool_t	eForceSpecificationReload

TH1D *	eHisto_deltaR

TH1D *	eHisto_deltaR_mean

TH1D *	eHisto_deltaR_rms

TH1D *	eHisto_deltaT

TH1D *	eHisto_deltaT_mean

TH1D *	eHisto_deltaT_rms

TH1D *	eHisto_longprofile

TH1D *	eHisto_longprofile_av

TH1D *	eHisto_nbtk

TH1D *	eHisto_nbtk_av

TH1D *	eHisto_transprofile

TH1D *	eHisto_transprofile_av

Float_t	eParaBT_deltaR_mean

Float_t	eParaBT_deltaR_rms

Float_t	eParaBT_deltaT_mean

Float_t	eParaBT_deltaT_rms

Int_t	eParalongprofile [57]

Int_t	eParaName

Int_t	eParanseg

Float_t	eParaShowerAxisAngle

Int_t	ePlateBinning [15]

Int_t	ePlateNumber

Int_t	ePlateNumberType

TObjArray *	eRecoShowerArray

Int_t	eRecoShowerArrayN

Bool_t	eSpecificationIsChanged

Int_t	eSpecificationType [7]

TString	eSpecificationTypeString [7]

TString	eSpecificationTypeStringArray [7][7]

TObjArray *	eSplineArray_Energy_Stat_Electron

TObjArray *	eSplineArray_Energy_Stat_Gamma

TObjArray *	eSplineArray_Energy_Sys_Electron

TObjArray *	eSplineArray_Energy_Sys_Gamma

TSpline3 *	eSplineCurrent

TString	eWeightFileString

Double_t	inANN [70]

Double_t	outANN

Additional Inherited Members
Protected Member Functions inherited from EdbShowerAlgESimple
void	Init ()

void	Set0 ()

Constructor & Destructor Documentation

◆ EdbShowerAlgIDSimple() [1/5]

EdbShowerAlgIDSimple::EdbShowerAlgIDSimple ( )

 {
     // Default Constructor
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple()   Default Constructor"<<endl;
     Set0();
     Init();
     cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple()   Default Constructor...done."<<endl;
 }

◆ EdbShowerAlgIDSimple() [2/5]

EdbShowerAlgIDSimple::EdbShowerAlgIDSimple ( EdbTrackP * track )

CONSTRUCTOR USED FOR libShowRec with EdbShowerP class available.

 {
     // Default Constructor...
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(EdbTrackP* track)   Default Constructor"<<endl;
     Set0();
     Init();
     // ... with Shower:
     eRecoShowerArray=new TObjArray();
     eRecoShowerArray->Add(track);
     SetRecoShowerArrayN(eRecoShowerArray->GetEntries());
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(EdbTrackP* track)   Default Constructor...done."<<endl;
 }

◆ EdbShowerAlgIDSimple() [3/5]

EdbShowerAlgIDSimple::EdbShowerAlgIDSimple ( TObjArray * RecoShowerArray )

 {
     // Default Constructor
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(TObjArray* RecoShowerArray)   Default Constructor"<<endl;
     Set0();
     Init();
     // ... with ShowerArray:
     eRecoShowerArray=RecoShowerArray;
     eRecoShowerArrayN=eRecoShowerArray->GetEntries();
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(TObjArray* RecoShowerArray)   Default Constructor...done."<<endl;
 }

◆ EdbShowerAlgIDSimple() [4/5]

EdbShowerAlgIDSimple::EdbShowerAlgIDSimple	(	EdbTrackP *	track,
		EdbVertex *	vtx
	)

 {
     // Default Constructor...
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(EdbTrackP* track, EdbVertex* vtx)   Default Constructor"<<endl;
     Set0();
     Init();
     // ... with Shower:
     eRecoShowerArray=new TObjArray();
     eRecoShowerArray->Add(track);
     SetRecoShowerArrayN(eRecoShowerArray->GetEntries());
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(EdbTrackP* track, EdbVertex* vtx)   Default Constructor...done."<<endl;
 }

◆ EdbShowerAlgIDSimple() [5/5]

EdbShowerAlgIDSimple::EdbShowerAlgIDSimple	(	TObjArray *	RecoShowerArray,
		TObjArray *	VtxArray
	)

 {
     // Default Constructor
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(TObjArray* RecoShowerArray, TObjArray* VtxArray)   Default Constructor"<<endl;
     Set0();
     Init();
     // ... with ShowerArray:
     eRecoShowerArray=RecoShowerArray;
     eRecoShowerArrayN=eRecoShowerArray->GetEntries();
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(TObjArray* RecoShowerArray, TObjArray* VtxArray)   Default Constructor...done."<<endl;
 }

◆ ~EdbShowerAlgIDSimple()

EdbShowerAlgIDSimple::~EdbShowerAlgIDSimple ( )

virtual

 {
     // Default Destructor
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::~EdbShowerAlgIDSimple()"<<endl;
     // Delete Histograms (on heap):
     delete eHisto_nbtk_av;
     delete eHisto_longprofile_av;
     delete eHisto_transprofile_av;
     delete eHisto_deltaR_mean;
     delete eHisto_deltaT_mean;
     delete eHisto_deltaR_rms;
     delete eHisto_deltaT_rms;
     delete eHisto_nbtk;
     delete eHisto_longprofile;
     delete eHisto_transprofile;
     delete eHisto_deltaR;
     delete eHisto_deltaT;
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::~EdbShowerAlgIDSimple()...done."<<endl;
 }

Member Function Documentation

◆ ClassDef()

EdbShowerAlgIDSimple::ClassDef	(	EdbShowerAlgIDSimple	,
		1
	)

◆ CreateANN()

void EdbShowerAlgIDSimple::CreateANN ( )

 {
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::CreateANN()"<<endl;
  
     // Compared to EdbShowerAlgESimple this ANN has two BIG differences:
     // a) It 2 InputNeurons more (vtxIP, vtxDeltaZ)
     // b) :inANN[0]/1000 is changed to :inANN[0]  (int instead of float variable)
  
  
     // Create ANN Tree and MLP:
     eParaName=2; // Standard Labeling taken over from EdbShowerP
     ANNTree = new TTree("ANNTree", "ANNTree");
     ANNTree->SetDirectory(0);
     ANNTree->Branch("inANN", inANN, "inANN[70]/D");
  
     //  Plate Binning: 10,12,14,16,18,20,23,26,29,32,35,40,45,45
     // see ePlateBinning[]
     for (int k=0; k<15; k++) {
 //      for (int k=0; k<1; k++) {
         cout << "EdbShowerAlgIDSimple::CreateANN() Number k= " << k << endl;
         ANN_Layout="";
         ANN_n_InputNeurons_ARRAY[k]=6+ePlateBinning[k]+2; // Two extra for ID (vtxIP, vtxDeltaZ)
         ANN_n_InputNeurons=ANN_n_InputNeurons_ARRAY[k];
         cout << "EdbShowerAlgIDSimple::CreateANN() Number ANN_n_InputNeurons= " << ANN_n_InputNeurons << endl;
  
         for (Int_t i=1; i<ANN_n_InputNeurons; ++i) ANN_Layout += "@inANN["+TString(Form("%d",i))+"],";
         ANN_Layout += "@inANN["+TString(Form("%d",ANN_n_InputNeurons))+"]:"+TString(Form("%d",ANN_n_InputNeurons+1))+":"+TString(Form("%d",ANN_n_InputNeurons));
         ANN_Layout+=":inANN[0]";
         //---------------------------
         ANN_MLP_ARRAY[k]    =   new TMultiLayerPerceptron(ANN_Layout,ANNTree,"","");
         ANN_MLP=ANN_MLP_ARRAY[k];
         ANN_Layout=ANN_MLP_ARRAY[k]->GetStructure();
         //---------------------------
         if (gEDBDEBUGLEVEL>2) {
             ANN_MLP->Print();
             cout << ANN_Layout << endl;
         }
     }
     //---------------------------
  
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::CreateANN()...done."<<endl;
     return;
 }

◆ DetermineIDFlag()

Int_t EdbShowerAlgIDSimple::DetermineIDFlag ( Int_t IDType )

 {
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::DetermineIDFlag()" << endl;
  
     // Depending on the SHID<->SPID type we assign the ID output into
     // 22, -11, -211 (PdgIds of Gamma,Electron,Pion)
  
     Int_t Flag=0;
     if (eSpecificationType[1]==2) {
         if (IDType==0) Flag=-211;
         if (IDType==1) Flag=22;
     }
     else if (eSpecificationType[1]==1) {
         if (IDType==0) Flag=-11;
         if (IDType==1) Flag=-211;
     }
     else {
         if (IDType==0) Flag=22;
         if (IDType==1) Flag=-11;
     }
     return Flag;
 }

◆ DetermineIDType()

Int_t EdbShowerAlgIDSimple::DetermineIDType	(	Double_t	val,
		Int_t	SpecificationType
	)

 {
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::DetermineIDType()" << endl;
  
     // Depending on cutthresholdvalue and -type we assign the ANN output into
     // 0 (i.e. SHID) or 1 (i.e. SPID):
     if (val>eIDCutTypeValue[SpecificationType]) return 1;
     return 1;
 }

◆ DoRun() [1/3]

void EdbShowerAlgIDSimple::DoRun ( )

 {
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::DoRun()" << endl;
  
     for (int i=0; i<eRecoShowerArrayN; i++) {
         if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   Doing i= " << i << endl;
         continue;
         eEnergyArrayCount=i;
  
         //EdbShowerP* shower=(EdbShowerP*) eRecoShowerArray->At(i);
         EdbTrackP* shower=(EdbTrackP*) eRecoShowerArray->At(i);
  
         DoRun(shower);
  
     } // (int i=0; i<eRecoShowerArrayN; i++)
  
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::DoRun()...done." << endl;
     return;
 }

◆ DoRun() [2/3]

void EdbShowerAlgIDSimple::DoRun ( EdbTrackP * shower )

Last two variables are reserved for vtxIP and vtxDeltaZ values!!

ALSO, IF first (two) plate information shall be used... DONE IN THE WEIGHTFILE (weight of pl1&2 is then zero), so dont have to do it here...

 {
  
     if (gEDBDEBUGLEVEL>2) cout << "EdbShowerAlgIDSimple::DoRun(EdbTrackP* shower)" << endl;
  
     // Get internal variables from the parametrisation filled:
     GetPara(shower);
  
     /*
          int EffNr=FindClosestEfficiencyParametrization(eParaShowerAxisAngle,eEfficiencyParametrisation->Eval(eParaShowerAxisAngle));
  
         // Check If Efficiency is the one we have or if we have to change/reload
         // the ANN weightfiles:
          if (GetSpecType(2)!=EffNr&&eForceSpecificationReload==kFALSE) {
             cout << "EdbShowerAlgIDSimple::DoRun()   INFO! Calulated Efficiency is more compatible with another one: Change Specifiaction! Call SetSpecificationType(2,"<< EffNr <<")." << endl;
             SetSpecificationType(2,EffNr);
         }
         if (GetSpecType(2)!=EffNr&&eForceSpecificationReload==kTRUE) {
             cout << "EdbShowerAlgIDSimple::DoRun()   INFO! Calulated Efficiency would be more compatible with another one: But eForceSpecificationReload is set to kFALSE. So we keep current Efficiency used." << endl;
         }
     */
  
  
  
     // This is to select the suited ANN to the shower, i.e. the one that matches closest
     // the number of plates:
     int check_Npl_index =0;
     GetNplIndexNr(shower->Npl(),check_Npl_index,ePlateNumberType);
  
  
     ANN_n_InputNeurons=ANN_n_InputNeurons_ARRAY[check_Npl_index];
     ANN_MLP=ANN_MLP_ARRAY[check_Npl_index];
     eWeightFileString=ANN_WeightFile_ARRAY[check_Npl_index];
     eCalibrationOffset = eANN_MLP_CORR_0[check_Npl_index];
     eCalibrationSlope = eANN_MLP_CORR_1[check_Npl_index];
  
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   ANN_n_InputNeurons_ARRAY[check_Npl_index]="<< ANN_n_InputNeurons_ARRAY[check_Npl_index] << endl;
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   ANN_MLP_ARRAY[check_Npl_index]="<< ANN_MLP_ARRAY[check_Npl_index] << endl;
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   Using the following layout: " << endl;
     if (gEDBDEBUGLEVEL >2) cout << ANN_MLP->GetStructure() << endl;
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   And the following weightfile: " << endl;
     if (gEDBDEBUGLEVEL >2) cout << eWeightFileString << endl;
  
     //ANN_MLP->LoadWeights(eWeightFileString);
  
     //  Reset InputVariables:
     for (int k=0; k<70; k++) {
         inANN[k]=0;
     }
  
  
     // inANN[0] is ALWAYS Reseverd for the quantity value to be estimated
     // (E,Id,...)
     // Test with private variables:
     inANN[1]=eParaShowerAxisAngle;
     inANN[2]=eParanseg;
     inANN[3]=eParaBT_deltaR_mean;
     inANN[4]=eParaBT_deltaR_rms;
     inANN[5]=eParaBT_deltaT_mean;
     inANN[6]=eParaBT_deltaT_rms;
     for (int ii=0; ii<57; ii++) {
         inANN[7+ii]= eParalongprofile[ii];
     }
  
  
  
     // Fill Tree:
     ANNTree->Fill();
     ANNTree->Show(0);
  
     Double_t params[70];
     Double_t val;
  
  
     cout << "EdbShowerAlgIDSimple::DoRun()   ------- TODO HERE:  INTERFACE WITH VERTEX TO CALCULATE VTXIP AND VTXDELTAZ TO SHOWER ----------"<<endl;
     params[ANN_n_InputNeurons-3]=-1;
     params[ANN_n_InputNeurons-2]=-1;
     cout << "EdbShowerAlgIDSimple::DoRun()   ------- TODO HERE:  INTERFACE WITH VERTEX TO CALCULATE VTXIP AND VTXDELTAZ TO SHOWER ----------"<<endl;
  
  
  
     for (Int_t j=1; j<=ANN_n_InputNeurons; ++j) {
         params[j-1]=inANN[j];
     }
  
     for (Int_t j=1; j<=ANN_n_InputNeurons; ++j) cout << "EdbShowerAlgIDSimple::DoRun()   :  j  params[j-1]=inANN[j]   " << j<< "  " << params[j-1] << endl;
  
  
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   Print Inputvariables...done." << endl;
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   Evaluate Neural Network Output:" << endl;
  
     // ---------------------------------
     // Evaluation of the ANN:
     val=(ANN_MLP->Evaluate(0,params));
     // ---------------------------------
  
     if ( ::isnan(val) ) {
         if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   ERROR! ANN_MLP->Evaluate(0,params) is NAN! Setting value to -1. " << endl;
         val=-1;
     }
  
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   After  Correction: ANN_MLP->Evaluate(0,params)= " <<  val << " (Inputvar=" << inANN[0]  << ")." << endl;
  
  
     // Determining the type: if its SHID like or if its SPID like...
     Int_t type=DetermineIDType(val,eSpecificationType[6]);
     DetermineIDFlag(type);
  
 //      if (val>eANNIDCutValue) cout << "val>eANNIDCutValue" << endl;
 //
 //      cout << "UUUUUUUUUUUUUUUUUU   val= " <<  val << endl;
 //
 //      shower->SetP(val);
 //
 //     cout << "EdbShowerAlgIDSimple::DoRun() Finally, set values..." <<  endl;
 //      cout << "EdbShowerAlgIDSimple::DoRun() /// VERY PRELIMINARY...." <<  endl;
 //      cout << "EdbShowerAlgIDSimple::DoRun() Finally, set values..." <<  endl;
  
  
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   ...Done." << endl;
     return;
 }

◆ DoRun() [3/3]

void EdbShowerAlgIDSimple::DoRun ( TObjArray * trackarray )

 {
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::DoRun(TObjArray* trackarray)" << endl;
  
     eRecoShowerArrayN=trackarray->GetEntries();
     eRecoShowerArray=trackarray;
  
     for (int i=0; i<eRecoShowerArrayN; i++) {
         if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun()   Doing i= " << i << endl;
  
         eEnergyArrayCount=i;
  
         //EdbShowerP* shower=(EdbShowerP*) eRecoShowerArray->At(i);
         EdbTrackP* shower=(EdbTrackP*) eRecoShowerArray->At(i);
  
         DoRun(shower);
  
     } // (int i=0; i<eRecoShowerArrayN; i++)
  
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::DoRun(TObjArray* trackarray)...done." << endl;
     return;
 }

◆ Help()

void EdbShowerAlgIDSimple::Help ( )

 {
     cout << "EdbShowerAlgIDSimple::Help"<<endl;
     cout << "EdbShowerAlgIDSimple::Help...done"<<endl;
     return;
 }

◆ Init()

void EdbShowerAlgIDSimple::Init ( void )

 {
     // This Plate Binnning is the _final_binning and will not be refined furthermore!
     ePlateBinning[0]=10;
     ePlateBinning[1]=12;
     ePlateBinning[2]=14;
     ePlateBinning[3]=16;
     ePlateBinning[4]=18;
     ePlateBinning[5]=20;
     ePlateBinning[6]=23;
     ePlateBinning[7]=26;
     ePlateBinning[8]=29;
     ePlateBinning[9]=32;
     ePlateBinning[10]=35;
     ePlateBinning[11]=40;
     ePlateBinning[12]=45;
     ePlateBinning[13]=45;
     ePlateBinning[14]=45; // 12,13,14 are the same, cause I have seen that it makes no sense to produce these,
     // since E,and sigmaE for 45 plates doesnt change anymore at all...
     // and how often !can! we scan more than 45 plates for a shower?
  
     ePlateNumber=20; // 20 plates by default...
     eSpecificationType[0] = 1;  // cp files by default...
     eSpecificationType[1] = 0;  // gamma_vs_electron by default...
     eSpecificationType[2] = 0;  // neuch eff by default...
     eSpecificationType[3] = 0;  // SHAPEONLY by default...
     eSpecificationType[4] = 0;  // next before weight by default...
     eSpecificationType[5] = 5;  // 20 plates by default...
     eSpecificationType[6] = 0;  // cutType_minDist by default...
  
     eSpecificationTypeString[0]="CP";
     eSpecificationTypeString[1]="gamma_vs_electron";
     eSpecificationTypeString[2]="Neuch";
     eSpecificationTypeString[3]="SHAPEONLY";
     eSpecificationTypeString[4]="Next";
     eSpecificationTypeString[5]="20";
     eSpecificationTypeString[6]="cutType_minDist";
  
     eForceSpecificationReload=kFALSE;
  
     for (int i=0; i<15; i++) {
         ANN_nPlates_ARRAY[i]=ePlateBinning[i];
     }
  
     EffFunc_all = new TF1("all","1.0-0.00000001*x",0,0.95);
     EffFunc_edefault = new TF1("default","0.94-0.216*x-0.767*x*x+1.865*x*x*x",0,0.95);
     EffFunc_elletroni = new TF1("elletroni","0.79+0.38*x-7.63*x*x+25.13*x*x*x-24.6*x*x*x*x",0,0.95);
     EffFunc_neuchmicro = new TF1("neuchmicro","0.94-0.955*x+1.80*x*x-0.95*x*x*x",0,0.95);
     EffFunc_MiddleFix = new TF1("MiddleFix","0.5*(0.888361-1.299*x+1.49198*x*x+1.64668*x*x*x-2.63758*x*x*x*x+0.79+0.38*x-7.63*x*x+25.13*x*x*x-24.6*x*x*x*x)",0,0.95);
     EffFunc_LowEff = new TF1("LowEff","0.85*0.5*(0.888361-1.299*x+1.49198*x*x+1.64668*x*x*x-2.63758*x*x*x*x+0.79+0.38*x-7.63*x*x+25.13*x*x*x-24.6*x*x*x*x)",0,0.95);
     EffFunc_UserEfficiency = new TF1("EffFunc_UserEfficiency","0.94-0.955*x+1.80*x*x-0.95*x*x*x",0,0.95); //new TSpline3();
  
  
     // Standard supposed efficiency of the showers/tracks that are to be evaluated:
     // The correct EffFunc_* is then choosen by taking the closes eff func to that one.
     // Initial parameters are the efficiency parameters for the "Neuch" parametrisation
     Double_t xarr[7]= {0,0.1,0.2,0.3,0.4,0.5,0.6};
     Double_t yarr[7]= {0.95,0.9,0.80,0.75,0.6,0.5,0.4};
     if (NULL!=eEfficiencyParametrisation) delete eEfficiencyParametrisation;
     eEfficiencyParametrisation = new TSpline3("",xarr,yarr,5,0,0,0.65);
  
     eHisto_nbtk_av          = new TH1D("eHisto_nbtk_av","Average basetrack numbers",21,0.0,10.0);
     eHisto_longprofile_av   = new TH1D("eHisto_longprofile_av","Basetracks per emulsion number",57,0.0,57.0);
     eHisto_transprofile_av  = new TH1D("eHisto_transprofile_av","Basetracks in trans distance",8,0.0,800.0);
     eHisto_deltaR_mean      = new TH1D("eHisto_deltaR_mean","Mean #deltar of all BTs in one shower",100,0.0,100.0);
     eHisto_deltaT_mean      = new TH1D("eHisto_deltaT_mean","Mean #delta#theta of all BTs in one shower",100,0.0,0.1);
     eHisto_deltaR_rms       = new TH1D("eHisto_deltaR_rms","RMS #deltar of all BTs in one shower",100,0.0,100.0);
     eHisto_deltaT_rms       = new TH1D("eHisto_deltaT_rms","RMS #delta#theta of all BTs in one shower",100,0.0,0.1);
     eHisto_nbtk             = new TH1D("eHisto_nbtk","Basetracks in the shower",50,0.0,100.0);
     eHisto_longprofile      = new TH1D("eHisto_longprofile","Basetracks per emulsion number",57,0.0,57.0);
     eHisto_deltaR           = new TH1D("eHisto_deltaR","Single #deltar of all BTs in Shower",100,0.0,150.0);
     eHisto_deltaT           = new TH1D("eHisto_deltaT","Single #delta#theta of all BTs in Shower",150,0.0,0.15);
     eHisto_transprofile     = new TH1D("eHisto_transprofile","Basetracks in trans distance",8,0.0,800.0);
  
     // CreateANN
     CreateANN();
  
     //First Time, we also Update()
     Update();
  
     // Set Strings:
     InitStrings();
  
     // Create EnergyArray
     eEnergyArray=new TArrayF(99999); // no way to adapt tarrayF on the fly...
     eEnergyArrayUnCorrected=new TArrayF(99999);
     eEnergyArraySigmaCorrected=new TArrayF(99999);
  
     // Read ID Lookup tables:
     ReadTables_ID();
  
     cout << "EdbShowerAlgIDSimple::Init()...done."<< endl;
     return;
 }

◆ InitStrings()

void EdbShowerAlgIDSimple::InitStrings ( )

 {
     eSpecificationTypeStringArray[0][0]="LT";
     eSpecificationTypeStringArray[0][1]="CP";
  
     eSpecificationTypeStringArray[1][0]="gamma_vs_electron";
     eSpecificationTypeStringArray[1][1]="electron_vs_pion";
     eSpecificationTypeStringArray[1][2]="pion_vs_gamma";
  
     eSpecificationTypeStringArray[2][0]="Neuch";
     eSpecificationTypeStringArray[2][1]="All";
     eSpecificationTypeStringArray[2][2]="MiddleFix";
     eSpecificationTypeStringArray[2][3]="Low";
  
     eSpecificationTypeStringArray[3][0]="SHAPEONLY";
     eSpecificationTypeStringArray[3][1]="ANDVTXIP";
     eSpecificationTypeStringArray[3][2]="ANDVTXDELTAZ";
     eSpecificationTypeStringArray[3][3]="SHAPEWOPL1";
     eSpecificationTypeStringArray[3][4]="SHAPEWOPL12";
  
     eSpecificationTypeStringArray[4][0]="Next";
     eSpecificationTypeStringArray[4][1]="Before";
  
     eSpecificationTypeStringArray[5][0]="10";
  
     eSpecificationTypeStringArray[6][0]="cutType_minDist";
     eSpecificationTypeStringArray[6][1]="cutType_HighEff";
     eSpecificationTypeStringArray[6][2]="cutType_HighPur";
  
     cout << "EdbShowerAlgIDSimple::InitStrings()...done."<< endl;
     return;
 }

◆ Print()

void EdbShowerAlgIDSimple::Print ( )

 {
     cout << "EdbShowerAlgIDSimple::Print"<<endl;
     cout << "EdbShowerAlgIDSimple::Print...done"<<endl;
     return;
 }

◆ PrintSpecifications()

void EdbShowerAlgIDSimple::PrintSpecifications ( )

                                                {
     cout << "EdbShowerAlgIDSimple::PrintSpecifications" << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationType[0]  (Dataset:    (linked tracks/full cp) LT/CP)   = " << GetSpecType(0) << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationType[1]  (ShowerID:   gamma_vs_electron/electron_vs_pion/pion_vs_gamma) = " << GetSpecType(1) << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationType[2]  (ScanEff:    Neuch/All/MiddleFix/LowEff) = " << GetSpecType(2) << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationType[3]  (E range:    0..20/0..40) = " << GetSpecType(3) << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationType[4]  (Npl weight: next before/after) = " <<GetSpecType(4) << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationType[5]  (Npl weight: 10,12,...,45) = " << GetSpecType(5) << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationType[6]  (CutType: MinDist/HighEff/HighPur) = " << GetSpecType(6) << endl;
     cout << "EdbShowerAlgIDSimple::   " << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationTypeString[0]  (take CP or linked_tracks)   = " << eSpecificationTypeString[0] << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationTypeString[1]  (gamma_vs_electron/electron_vs_pion/pion_vs_gamma) = " << eSpecificationTypeString[1] << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationTypeString[2]  (ScanEff: Neuch/All/MiddleFix/LowEff) = " << eSpecificationTypeString[2] << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationTypeString[3]  (E range: 0..20/0..40) = " << eSpecificationTypeString[3] << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationTypeString[4]  (Npl weight: next before/after) = " << eSpecificationTypeString[4] << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationTypeString[5]  (Npl weight: 10,12,...,45) = " << eSpecificationTypeString[5] << endl;
     cout << "EdbShowerAlgIDSimple::   eSpecificationTypeString[6]  (CutType: MinDist/HighEff/HighPur) = " << eSpecificationTypeString[6] << endl;
     cout << "EdbShowerAlgIDSimple::   " << endl;
     cout << "EdbShowerAlgIDSimple::   In case you want to change a specification then do for example:" << endl;
     cout << "EdbShowerAlgIDSimple::   EdbShowerAlgIDSimple->SetSpecification(2,3) " << endl;
     cout << "EdbShowerAlgIDSimple::   It will change the ScanEff from the actual value to _Low_." << endl;
     cout << "EdbShowerAlgIDSimple::" << endl;
     cout << "EdbShowerAlgIDSimple::PrintSpecifications...done." << endl;
     return;
 }

◆ ReadTables()

void EdbShowerAlgIDSimple::ReadTables ( )

 {
     cout << "EdbShowerAlgIDSimple::ReadTables"<<endl;
     ReadTables_ID();
     cout << "EdbShowerAlgIDSimple::ReadTables...done."<<endl;
     return;
 }

◆ ReadTables_ID()

void EdbShowerAlgIDSimple::ReadTables_ID ( )

First guess until tables are fully read in.

 {
     cout << "EdbShowerAlgIDSimple::ReadTables_ID"<<endl;
  
     cout << "First guess until tables are fully read in." << endl;
     eIDCutTypeValue[0]=0.5;
     eIDCutTypeValue[1]=0.75;
     eIDCutTypeValue[2]=0.25;
     cout << "EdbShowerAlgIDSimple::ReadTables_ID...done"<<endl;
     return;
 }

◆ SetANNIDCutValue()

void EdbShowerAlgIDSimple::SetANNIDCutValue ( Double_t ANNIDCutValue )

inline

                                                          {
         eANNIDCutValue=ANNIDCutValue;
     }

◆ SetSpecificationType()

void EdbShowerAlgIDSimple::SetSpecificationType	(	Int_t	SpecificationType,
		Int_t	SpecificationTypeVal
	)

 {
  
     if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple:: Change Specification (" << SpecificationType << ") from " << GetSpecType(SpecificationType) << " -> " << SpecificationTypeVal << " . Reprint the changed Specification String: " << endl;
  
     eSpecificationType[SpecificationType]=SpecificationTypeVal;
     eSpecificationTypeString[SpecificationType] = eSpecificationTypeStringArray[SpecificationType][SpecificationTypeVal];
  
     if (gEDBDEBUGLEVEL >1) cout << eSpecificationTypeString[SpecificationType]  << endl;
  
     eSpecificationIsChanged=kTRUE;
     eForceSpecificationReload=kTRUE;
  
     cout << "EdbShowerAlgIDSimple::SetSpecificationType eSpecificationIsChanged: " << eSpecificationIsChanged << endl;
     cout << "EdbShowerAlgIDSimple::SetSpecificationType eForceSpecificationReload: " << eForceSpecificationReload << endl;
  
     cout << "Reprint Specifications: " << endl;
     if (gEDBDEBUGLEVEL >1) PrintSpecifications();
  
     Update();
     return;
 }

◆ Update()

void EdbShowerAlgIDSimple::Update ( )

                                   {
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update  Does the following things in the order:" << endl;
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update  * According to the switch: set the right ANN of the Array as generic one." << endl;
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update  * According to the switch: load the right weightfile as generic one." << endl;
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update  * According to the switch: set the right correction parameters...." << endl;
  
     TString basicstring = TString(gSystem->ExpandPathName("$FEDRA_ROOT"));
     TString addstring = ("/src/libShower/weights/ID/");
  
     // (*) SpecType 0 :
     if (eSpecificationType[0]==1) {
         addstring+="volumeSpec_CP/";
         eSpecificationTypeString[0]="CP";
     }
     else {
         addstring+="volumeSpec_LT/";
         eSpecificationTypeString[0]="LT";
         cout << "EdbShowerAlgIDSimple::Update  WARNING::eSpecificationTypeString[0]=LT NOT YET SUPPORTED!" << endl;
     }
  
     // (*) SpecType 1 :
     // SHID_vs_SPID
     if (eSpecificationType[1]==2) {
         addstring+="separationSpec_pion_vs_gamma/";
         eSpecificationTypeString[1]="pion_vs_gamma";
     }
     else if (eSpecificationType[1]==1)  {
         addstring+="separationSpec_electron_vs_pion/";
         eSpecificationTypeString[1]="electron_vs_pion";
     }
     else {
         addstring+="separationSpec_gamma_vs_electron/";
         eSpecificationTypeString[1]="gamma_vs_electron";
     }
  
     // (*) SpecType 2 :
     if (eSpecificationType[2]==4) {
         addstring+="efficiencySpec_MiddleVar/";
         eSpecificationTypeString[2]="MiddleVar";
     }
     if (eSpecificationType[2]==3) {
         addstring+="efficiencySpec_LowEff/";
         eSpecificationTypeString[2]="LowEff";
     }
     else if (eSpecificationType[2]==2)  {
         addstring+="efficiencySpec_MiddleFix/";
         eSpecificationTypeString[2]="MiddleFix";
     }
     else if (eSpecificationType[2]==1)  {
         addstring+="efficiencySpec_All/";
         eSpecificationTypeString[2]="All";
     }
     else {
         addstring+="efficiencySpec_Neuch/";
         eSpecificationTypeString[2]="Neuch";
     }
  
     // (*) SpecType 3 :
     // Which addiditonal variables to use....
 //  0: showershape only
 //  1: and vtx IP
 //  2: and vtx delta Z
 //  3: shoowershape only, No first plate
 //  4: shoowershape only, No first & second plate
     if (eSpecificationType[3]==4) {
         addstring+="showerInfoSpec_SHAPEWOPL12/";
         eSpecificationTypeString[3]="SHAPEWOPL12";
     }
     else if (eSpecificationType[3]==3) {
         addstring+="showerInfoSpec_SHAPEWOPL1/";
         eSpecificationTypeString[3]="SHAPEWOPL1";
     }
     else if (eSpecificationType[3]==2) {
         addstring+="showerInfoSpec_ANDVTXDELTAZ/";
         eSpecificationTypeString[3]="ANDVTXDELTAZ";
     }
     else if (eSpecificationType[3]==1) {
         addstring+="showerInfoSpec_ANDVTXIP/";
         eSpecificationTypeString[3]="ANDVTXIP";
     }
     else {
         addstring+="showerInfoSpec_SHAPEONLY/";
         eSpecificationTypeString[3]="SHAPEONLY";
     }
  
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update  " << addstring << endl;
  
     if (NULL==ANN_MLP_ARRAY[0]) {
         CreateANN();
     }
  
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update  Now load the different ANN weightstrings:" << endl;
     // This was the loading part for the ANN and
     // this is the loading part for the Correction Factors, which we take from the
     // generic file EnergyCorrections_Npl_%d.txt
     // For the ID class we do NOT need the correction factors, so we skip these files anyway....
     for (int ll=0; ll<15; ll++) {
         TString mkdirstring=basicstring+addstring;
         TString weigthstring=basicstring+addstring+TString(Form("weights_Npl_%d.txt",ANN_nPlates_ARRAY[ll]));
         if (gEDBDEBUGLEVEL >1) cout << "weigthstring = " << weigthstring << endl;
         ANN_MLP_ARRAY[ll]->LoadWeights(weigthstring);
         ANN_WeightFile_ARRAY[ll]=weigthstring;
         Float_t p0,p1;
         p0=0.0;
         p1=1.0;
         eANN_MLP_CORR_0[ll]=p0;
         eANN_MLP_CORR_1[ll]=p1;
     }
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update  * According to the switch: set the right ANN of the Array as generic one." << endl;
  
     // (*) SpecType 4:
     // Get NplIndexNr:
     int check_Npl_index=0;
     GetNplIndexNr(ePlateNumber,check_Npl_index, eSpecificationType[4]);
  
     if (check_Npl_index==0) eSpecificationTypeString[4]="Next";
     if (check_Npl_index==1) eSpecificationTypeString[4]="Before";
     eSpecificationTypeString[5]=TString(Form("%d",ePlateNumber));
  
     // (*) SpecType 5:
     // HigEfficiency Cut, HighPurityCut, MinDistCut
     if (eSpecificationType[6]==1) {
         addstring+="cutType/";
         eSpecificationTypeString[6]="cutType_HighEff";
     }
     else if (eSpecificationType[6]==2) {
         addstring+="cutType/";
         eSpecificationTypeString[6]="cutType_HighPur";
     }
     else {
         addstring+="cutType/";
         eSpecificationTypeString[6]="cutType_minDist";
     }
  
     if (gEDBDEBUGLEVEL >2)  cout << "EdbShowerAlgIDSimple:: The path for the ID weightfile is. " << endl;
     if (gEDBDEBUGLEVEL >2)  cout << addstring << endl;
  
  
     // Set Generic ANN_Layout
     ANN_MLP=ANN_MLP_ARRAY[check_Npl_index];
  
  
     // Re-Read ID lookup tables, because of the possibility to
     // change the cutType threshold.
     ReadTables_ID();
  
     if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update  WARNING:: Weightfiles for _LT_  not produced yet!...." << endl;
     cout << "EdbShowerAlgIDSimple::Update...done." << endl;
     return;
 }

Member Data Documentation

◆ eANNIDCutValue

Double_t EdbShowerAlgIDSimple::eANNIDCutValue

protected

◆ eIDCutTypeValue

Double_t EdbShowerAlgIDSimple::eIDCutTypeValue[3]

protected

◆ eIDEnergySigmaCorr

Double_t EdbShowerAlgIDSimple::eIDEnergySigmaCorr

protected

The documentation for this class was generated from the following files:

src/libShower/EdbShowerAlg.h
src/libShower/EdbShowerAlg.cxx

Public Member Functions

Protected Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ EdbShowerAlgIDSimple() [1/5]

◆ EdbShowerAlgIDSimple() [2/5]

◆ EdbShowerAlgIDSimple() [3/5]

◆ EdbShowerAlgIDSimple() [4/5]

◆ EdbShowerAlgIDSimple() [5/5]

◆ ~EdbShowerAlgIDSimple()

Member Function Documentation

◆ ClassDef()

◆ CreateANN()

◆ DetermineIDFlag()

◆ DetermineIDType()

◆ DoRun() [1/3]

◆ DoRun() [2/3]

◆ DoRun() [3/3]

◆ Help()

◆ Init()

◆ InitStrings()

◆ Print()

◆ PrintSpecifications()

◆ ReadTables()

◆ ReadTables_ID()

◆ SetANNIDCutValue()

◆ SetSpecificationType()

◆ Update()

Member Data Documentation

◆ eANNIDCutValue

◆ eIDCutTypeValue

◆ eIDEnergySigmaCorr