EdbShowerAlgIDSimple Class Reference
#include <EdbShowerAlg.h>
Inheritance diagram for EdbShowerAlgIDSimple:
Collaboration diagram for EdbShowerAlgIDSimple:
Additional Inherited Members | |
 Protected Member Functions inherited from EdbShowerAlgESimple | |
| void | Init () |
| void | Set0 () |
Constructor & Destructor Documentation
◆ EdbShowerAlgIDSimple() [1/5]
| EdbShowerAlgIDSimple::EdbShowerAlgIDSimple | ( | ) |
Default Constructor
4122{
4124 if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple() Default Constructor"<<endl;
4125 Set0();
4126 Init();
4127 cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple() Default Constructor...done."<<endl;
4128}
◆ EdbShowerAlgIDSimple() [2/5]
| EdbShowerAlgIDSimple::EdbShowerAlgIDSimple | ( | EdbTrackP * | track | ) |
CONSTRUCTOR USED FOR libShowRec with EdbShowerP class available.
Default Constructor...
4150{
4152 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(EdbTrackP* track) Default Constructor"<<endl;
4153 Set0();
4154 Init();
4155 // ... with Shower:
4159 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(EdbTrackP* track) Default Constructor...done."<<endl;
4160}
void SetRecoShowerArrayN(Int_t RecoShowerArrayN)
Definition: EdbShowerAlg.h:452
Definition: bitview.h:14
◆ EdbShowerAlgIDSimple() [3/5]
| EdbShowerAlgIDSimple::EdbShowerAlgIDSimple | ( | TObjArray * | RecoShowerArray | ) |
Default Constructor
4166{
4168 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(TObjArray* RecoShowerArray) Default Constructor"<<endl;
4169 Set0();
4170 Init();
4171 // ... with ShowerArray:
4174 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(TObjArray* RecoShowerArray) Default Constructor...done."<<endl;
4175}
◆ EdbShowerAlgIDSimple() [4/5]
Default Constructor...
4180{
4182 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(EdbTrackP* track, EdbVertex* vtx) Default Constructor"<<endl;
4183 Set0();
4184 Init();
4185 // ... with Shower:
4189 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(EdbTrackP* track, EdbVertex* vtx) Default Constructor...done."<<endl;
4190}
◆ EdbShowerAlgIDSimple() [5/5]
| EdbShowerAlgIDSimple::EdbShowerAlgIDSimple | ( | TObjArray * | RecoShowerArray, |
| TObjArray * | VtxArray | ||
| ) |
Default Constructor
4196{
4198 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(TObjArray* RecoShowerArray, TObjArray* VtxArray) Default Constructor"<<endl;
4199 Set0();
4200 Init();
4201 // ... with ShowerArray:
4204 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::EdbShowerAlgIDSimple(TObjArray* RecoShowerArray, TObjArray* VtxArray) Default Constructor...done."<<endl;
4205}
◆ ~EdbShowerAlgIDSimple()
|
virtual |
virtual constructor due to inherited class
Default Destructor
4210{
4213 // Delete Histograms (on heap):
4227}
TH1D * eHisto_transprofile_av
Definition: EdbShowerAlg.h:389
TH1D * eHisto_longprofile_av
Definition: EdbShowerAlg.h:388
Member Function Documentation
◆ ClassDef()
| EdbShowerAlgIDSimple::ClassDef | ( | EdbShowerAlgIDSimple | , |
| 1 | |||
| ) |
◆ CreateANN()
| void EdbShowerAlgIDSimple::CreateANN | ( | ) |
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)
4592{
4594
4598
4599
4600 // Create ANN Tree and MLP:
4603 ANNTree->SetDirectory(0);
4605
4606 // Plate Binning: 10,12,14,16,18,20,23,26,29,32,35,40,45,45
4607 // see ePlateBinning[]
4608 for (int k=0; k<15; k++) {
4609// for (int k=0; k<1; k++) {
4610 cout << "EdbShowerAlgIDSimple::CreateANN() Number k= " << k << endl;
4614 cout << "EdbShowerAlgIDSimple::CreateANN() Number ANN_n_InputNeurons= " << ANN_n_InputNeurons << endl;
4615
4617 ANN_Layout += "@inANN["+TString(Form("%d",ANN_n_InputNeurons))+"]:"+TString(Form("%d",ANN_n_InputNeurons+1))+":"+TString(Form("%d",ANN_n_InputNeurons));
4619 //---------------------------
4623 //---------------------------
4625 ANN_MLP->Print();
4626 cout << ANN_Layout << endl;
4627 }
4628 }
4629 //---------------------------
4630
4632 return;
4633}
Int_t ANN_n_InputNeurons_ARRAY[15]
Definition: EdbShowerAlg.h:367
TMultiLayerPerceptron * ANN_MLP_ARRAY[15]
Definition: EdbShowerAlg.h:362
◆ DetermineIDFlag()
| Int_t EdbShowerAlgIDSimple::DetermineIDFlag | ( | Int_t | IDType | ) |
4783{
4785
4786 // Depending on the SHID<->SPID type we assign the ID output into
4787 // 22, -11, -211 (PdgIds of Gamma,Electron,Pion)
4788
4789 Int_t Flag=0;
4791 if (IDType==0) Flag=-211;
4792 if (IDType==1) Flag=22;
4793 }
4795 if (IDType==0) Flag=-11;
4796 if (IDType==1) Flag=-211;
4797 }
4798 else {
4799 if (IDType==0) Flag=22;
4800 if (IDType==1) Flag=-11;
4801 }
4802 return Flag;
4803}
◆ DetermineIDType()
| Int_t EdbShowerAlgIDSimple::DetermineIDType | ( | Double_t | val, |
| Int_t | SpecificationType | ||
| ) |
4771{
4773
4774 // Depending on cutthresholdvalue and -type we assign the ANN output into
4775 // 0 (i.e. SHID) or 1 (i.e. SPID):
4777 return 1;
4778}
◆ DoRun() [1/3]
| void EdbShowerAlgIDSimple::DoRun | ( | ) |
4808{
4810
4813 continue;
4814 eEnergyArrayCount=i;
4815
4816 //EdbShowerP* shower=(EdbShowerP*) eRecoShowerArray->At(i);
4818
4819 DoRun(shower);
4820
4821 } // (int i=0; i<eRecoShowerArrayN; i++)
4822
4824 return;
4825}
Definition: EdbPattern.h:113
◆ DoRun() [2/3]
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...
4642{
4643
4645
4646 // Get internal variables from the parametrisation filled:
4647 GetPara(shower);
4648
4649 /*
4650 int EffNr=FindClosestEfficiencyParametrization(eParaShowerAxisAngle,eEfficiencyParametrisation->Eval(eParaShowerAxisAngle));
4651
4652 // Check If Efficiency is the one we have or if we have to change/reload
4653 // the ANN weightfiles:
4654 if (GetSpecType(2)!=EffNr&&eForceSpecificationReload==kFALSE) {
4655 cout << "EdbShowerAlgIDSimple::DoRun() INFO! Calulated Efficiency is more compatible with another one: Change Specifiaction! Call SetSpecificationType(2,"<< EffNr <<")." << endl;
4656 SetSpecificationType(2,EffNr);
4657 }
4658 if (GetSpecType(2)!=EffNr&&eForceSpecificationReload==kTRUE) {
4659 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;
4660 }
4661 */
4662
4663
4664
4665 // This is to select the suited ANN to the shower, i.e. the one that matches closest
4666 // the number of plates:
4667 int check_Npl_index =0;
4669
4670
4676
4677 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun() ANN_n_InputNeurons_ARRAY[check_Npl_index]="<< ANN_n_InputNeurons_ARRAY[check_Npl_index] << endl;
4678 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun() ANN_MLP_ARRAY[check_Npl_index]="<< ANN_MLP_ARRAY[check_Npl_index] << endl;
4679 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun() Using the following layout: " << endl;
4681 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun() And the following weightfile: " << endl;
4683
4684 //ANN_MLP->LoadWeights(eWeightFileString);
4685
4686 // Reset InputVariables:
4687 for (int k=0; k<70; k++) {
4688 inANN[k]=0;
4689 }
4690
4691
4692 // inANN[0] is ALWAYS Reseverd for the quantity value to be estimated
4693 // (E,Id,...)
4694 // Test with private variables:
4701 for (int ii=0; ii<57; ii++) {
4703 }
4705
4709
4710
4711 // Fill Tree:
4712 ANNTree->Fill();
4713 ANNTree->Show(0);
4714
4715 Double_t params[70];
4716 Double_t val;
4717
4718
4719 cout << "EdbShowerAlgIDSimple::DoRun() ------- TODO HERE: INTERFACE WITH VERTEX TO CALCULATE VTXIP AND VTXDELTAZ TO SHOWER ----------"<<endl;
4722 cout << "EdbShowerAlgIDSimple::DoRun() ------- TODO HERE: INTERFACE WITH VERTEX TO CALCULATE VTXIP AND VTXDELTAZ TO SHOWER ----------"<<endl;
4723
4724
4725
4728 }
4729
4730 for (Int_t j=1; j<=ANN_n_InputNeurons; ++j) cout << "EdbShowerAlgIDSimple::DoRun() : j params[j-1]=inANN[j] " << j<< " " << params[j-1] << endl;
4731
4732
4733 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun() Print Inputvariables...done." << endl;
4734 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun() Evaluate Neural Network Output:" << endl;
4735
4736 // ---------------------------------
4737 // Evaluation of the ANN:
4739 // ---------------------------------
4740
4741 if ( ::isnan(val) ) {
4742 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun() ERROR! ANN_MLP->Evaluate(0,params) is NAN! Setting value to -1. " << endl;
4743 val=-1;
4744 }
4745
4746 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::DoRun() After Correction: ANN_MLP->Evaluate(0,params)= " << val << " (Inputvar=" << inANN[0] << ")." << endl;
4747
4748
4749 // Determining the type: if its SHID like or if its SPID like...
4752
4753// if (val>eANNIDCutValue) cout << "val>eANNIDCutValue" << endl;
4754//
4755// cout << "UUUUUUUUUUUUUUUUUU val= " << val << endl;
4756//
4757// shower->SetP(val);
4758//
4759// cout << "EdbShowerAlgIDSimple::DoRun() Finally, set values..." << endl;
4760// cout << "EdbShowerAlgIDSimple::DoRun() /// VERY PRELIMINARY...." << endl;
4761// cout << "EdbShowerAlgIDSimple::DoRun() Finally, set values..." << endl;
4762
4763
4765 return;
4766}
void GetNplIndexNr(Int_t sh_npl, Int_t &check_Npl_index, Int_t ePlateNumberType)
Definition: EdbShowerAlg.cxx:2993
TString ANN_WeightFile_ARRAY[15]
Definition: EdbShowerAlg.h:370
Int_t DetermineIDFlag(Int_t IDType)
Definition: EdbShowerAlg.cxx:4782
Int_t DetermineIDType(Double_t val, Int_t SpecificationType)
Definition: EdbShowerAlg.cxx:4770
◆ DoRun() [3/3]
| void EdbShowerAlgIDSimple::DoRun | ( | TObjArray * | trackarray | ) |
4831{
4833
4834 eRecoShowerArrayN=trackarray->GetEntries();
4835 eRecoShowerArray=trackarray;
4836
4839
4840 eEnergyArrayCount=i;
4841
4842 //EdbShowerP* shower=(EdbShowerP*) eRecoShowerArray->At(i);
4844
4845 DoRun(shower);
4846
4847 } // (int i=0; i<eRecoShowerArrayN; i++)
4848
4849 if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple::DoRun(TObjArray* trackarray)...done." << endl;
4850 return;
4851}
◆ Help()
| void EdbShowerAlgIDSimple::Help | ( | ) |
4270{
4271 cout << "EdbShowerAlgIDSimple::Help"<<endl;
4272 cout << "EdbShowerAlgIDSimple::Help...done"<<endl;
4273 return;
4274}
◆ Init()
This Plate Binnning is the _final_binning and will not be refined furthermore!
4459{
4461 ePlateBinning[0]=10;
4462 ePlateBinning[1]=12;
4463 ePlateBinning[2]=14;
4464 ePlateBinning[3]=16;
4465 ePlateBinning[4]=18;
4466 ePlateBinning[5]=20;
4467 ePlateBinning[6]=23;
4468 ePlateBinning[7]=26;
4469 ePlateBinning[8]=29;
4470 ePlateBinning[9]=32;
4471 ePlateBinning[10]=35;
4472 ePlateBinning[11]=40;
4473 ePlateBinning[12]=45;
4474 ePlateBinning[13]=45;
4475 ePlateBinning[14]=45; // 12,13,14 are the same, cause I have seen that it makes no sense to produce these,
4476 // since E,and sigmaE for 45 plates doesnt change anymore at all...
4477 // and how often !can! we scan more than 45 plates for a shower?
4478
4487
4495
4496 eForceSpecificationReload=kFALSE;
4497
4498 for (int i=0; i<15; i++) {
4500 }
4501
4504 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);
4506 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);
4507 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);
4508 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();
4509
4510
4511 // Standard supposed efficiency of the showers/tracks that are to be evaluated:
4512 // The correct EffFunc_* is then choosen by taking the closes eff func to that one.
4513 // Initial parameters are the efficiency parameters for the "Neuch" parametrisation
4514 Double_t xarr[7]= {0,0.1,0.2,0.3,0.4,0.5,0.6};
4515 Double_t yarr[7]= {0.95,0.9,0.80,0.75,0.6,0.5,0.4};
4518
4520 eHisto_longprofile_av = new TH1D("eHisto_longprofile_av","Basetracks per emulsion number",57,0.0,57.0);
4521 eHisto_transprofile_av = new TH1D("eHisto_transprofile_av","Basetracks in trans distance",8,0.0,800.0);
4522 eHisto_deltaR_mean = new TH1D("eHisto_deltaR_mean","Mean #deltar of all BTs in one shower",100,0.0,100.0);
4523 eHisto_deltaT_mean = new TH1D("eHisto_deltaT_mean","Mean #delta#theta of all BTs in one shower",100,0.0,0.1);
4524 eHisto_deltaR_rms = new TH1D("eHisto_deltaR_rms","RMS #deltar of all BTs in one shower",100,0.0,100.0);
4525 eHisto_deltaT_rms = new TH1D("eHisto_deltaT_rms","RMS #delta#theta of all BTs in one shower",100,0.0,0.1);
4527 eHisto_longprofile = new TH1D("eHisto_longprofile","Basetracks per emulsion number",57,0.0,57.0);
4529 eHisto_deltaT = new TH1D("eHisto_deltaT","Single #delta#theta of all BTs in Shower",150,0.0,0.15);
4530 eHisto_transprofile = new TH1D("eHisto_transprofile","Basetracks in trans distance",8,0.0,800.0);
4531
4532 // CreateANN
4533 CreateANN();
4534
4535 //First Time, we also Update()
4536 Update();
4537
4538 // Set Strings:
4539 InitStrings();
4540
4541 // Create EnergyArray
4545
4546 // Read ID Lookup tables:
4547 ReadTables_ID();
4548
4549 cout << "EdbShowerAlgIDSimple::Init()...done."<< endl;
4550 return;
4551}
Bool_t eForceSpecificationReload
Definition: EdbShowerAlg.h:354
TSpline3 * eEfficiencyParametrisation
Definition: EdbShowerAlg.h:412
TArrayF * eEnergyArraySigmaCorrected
Definition: EdbShowerAlg.h:417
TArrayF * eEnergyArrayUnCorrected
Definition: EdbShowerAlg.h:416
TF1 * EffFunc_UserEfficiency
Definition: EdbShowerAlg.h:407
TString eSpecificationTypeString[7]
Definition: EdbShowerAlg.h:351
◆ InitStrings()
| void EdbShowerAlgIDSimple::InitStrings | ( | ) |
4557{
4560
4564
4569
4575
4578
4580
4584
4585 cout << "EdbShowerAlgIDSimple::InitStrings()...done."<< endl;
4586 return;
4587}
TString eSpecificationTypeStringArray[7][7]
Definition: EdbShowerAlg.h:352
◆ Print()
| void EdbShowerAlgIDSimple::Print | ( | ) |
4260{
4261 cout << "EdbShowerAlgIDSimple::Print"<<endl;
4262 cout << "EdbShowerAlgIDSimple::Print...done"<<endl;
4263 return;
4264}
◆ PrintSpecifications()
| void EdbShowerAlgIDSimple::PrintSpecifications | ( | ) |
4857 {
4858 cout << "EdbShowerAlgIDSimple::PrintSpecifications" << endl;
4859 cout << "EdbShowerAlgIDSimple:: eSpecificationType[0] (Dataset: (linked tracks/full cp) LT/CP) = " << GetSpecType(0) << endl;
4860 cout << "EdbShowerAlgIDSimple:: eSpecificationType[1] (ShowerID: gamma_vs_electron/electron_vs_pion/pion_vs_gamma) = " << GetSpecType(1) << endl;
4861 cout << "EdbShowerAlgIDSimple:: eSpecificationType[2] (ScanEff: Neuch/All/MiddleFix/LowEff) = " << GetSpecType(2) << endl;
4862 cout << "EdbShowerAlgIDSimple:: eSpecificationType[3] (E range: 0..20/0..40) = " << GetSpecType(3) << endl;
4863 cout << "EdbShowerAlgIDSimple:: eSpecificationType[4] (Npl weight: next before/after) = " <<GetSpecType(4) << endl;
4864 cout << "EdbShowerAlgIDSimple:: eSpecificationType[5] (Npl weight: 10,12,...,45) = " << GetSpecType(5) << endl;
4865 cout << "EdbShowerAlgIDSimple:: eSpecificationType[6] (CutType: MinDist/HighEff/HighPur) = " << GetSpecType(6) << endl;
4866 cout << "EdbShowerAlgIDSimple:: " << endl;
4867 cout << "EdbShowerAlgIDSimple:: eSpecificationTypeString[0] (take CP or linked_tracks) = " << eSpecificationTypeString[0] << endl;
4868 cout << "EdbShowerAlgIDSimple:: eSpecificationTypeString[1] (gamma_vs_electron/electron_vs_pion/pion_vs_gamma) = " << eSpecificationTypeString[1] << endl;
4869 cout << "EdbShowerAlgIDSimple:: eSpecificationTypeString[2] (ScanEff: Neuch/All/MiddleFix/LowEff) = " << eSpecificationTypeString[2] << endl;
4870 cout << "EdbShowerAlgIDSimple:: eSpecificationTypeString[3] (E range: 0..20/0..40) = " << eSpecificationTypeString[3] << endl;
4871 cout << "EdbShowerAlgIDSimple:: eSpecificationTypeString[4] (Npl weight: next before/after) = " << eSpecificationTypeString[4] << endl;
4872 cout << "EdbShowerAlgIDSimple:: eSpecificationTypeString[5] (Npl weight: 10,12,...,45) = " << eSpecificationTypeString[5] << endl;
4873 cout << "EdbShowerAlgIDSimple:: eSpecificationTypeString[6] (CutType: MinDist/HighEff/HighPur) = " << eSpecificationTypeString[6] << endl;
4874 cout << "EdbShowerAlgIDSimple:: " << endl;
4875 cout << "EdbShowerAlgIDSimple:: In case you want to change a specification then do for example:" << endl;
4876 cout << "EdbShowerAlgIDSimple:: EdbShowerAlgIDSimple->SetSpecification(2,3) " << endl;
4877 cout << "EdbShowerAlgIDSimple:: It will change the ScanEff from the actual value to _Low_." << endl;
4878 cout << "EdbShowerAlgIDSimple::" << endl;
4879 cout << "EdbShowerAlgIDSimple::PrintSpecifications...done." << endl;
4880 return;
4881}
Int_t GetSpecType(Int_t SpecificationType)
Definition: EdbShowerAlg.h:527
◆ ReadTables()
| void EdbShowerAlgIDSimple::ReadTables | ( | ) |
4233{
4234 cout << "EdbShowerAlgIDSimple::ReadTables"<<endl;
4235 ReadTables_ID();
4236 cout << "EdbShowerAlgIDSimple::ReadTables...done."<<endl;
4237 return;
4238}
◆ ReadTables_ID()
| void EdbShowerAlgIDSimple::ReadTables_ID | ( | ) |
First guess until tables are fully read in.
4244{
4245 cout << "EdbShowerAlgIDSimple::ReadTables_ID"<<endl;
4246
4248 cout << "First guess until tables are fully read in." << endl;
4249 eIDCutTypeValue[0]=0.5;
4250 eIDCutTypeValue[1]=0.75;
4251 eIDCutTypeValue[2]=0.25;
4252 cout << "EdbShowerAlgIDSimple::ReadTables_ID...done"<<endl;
4253 return;
4254}
◆ SetANNIDCutValue()
|
inline |
◆ SetSpecificationType()
| void EdbShowerAlgIDSimple::SetSpecificationType | ( | Int_t | SpecificationType, |
| Int_t | SpecificationTypeVal | ||
| ) |
4280{
4281
4282 if (gEDBDEBUGLEVEL >1) cout << "EdbShowerAlgIDSimple:: Change Specification (" << SpecificationType << ") from " << GetSpecType(SpecificationType) << " -> " << SpecificationTypeVal << " . Reprint the changed Specification String: " << endl;
4283
4284 eSpecificationType[SpecificationType]=SpecificationTypeVal;
4285 eSpecificationTypeString[SpecificationType] = eSpecificationTypeStringArray[SpecificationType][SpecificationTypeVal];
4286
4288
4289 eSpecificationIsChanged=kTRUE;
4290 eForceSpecificationReload=kTRUE;
4291
4292 cout << "EdbShowerAlgIDSimple::SetSpecificationType eSpecificationIsChanged: " << eSpecificationIsChanged << endl;
4293 cout << "EdbShowerAlgIDSimple::SetSpecificationType eForceSpecificationReload: " << eForceSpecificationReload << endl;
4294
4295 cout << "Reprint Specifications: " << endl;
4297
4298 Update();
4299 return;
4300}
Bool_t eSpecificationIsChanged
Definition: EdbShowerAlg.h:353
void PrintSpecifications()
Definition: EdbShowerAlg.cxx:4857
◆ Update()
| void EdbShowerAlgIDSimple::Update | ( | ) |
4304 {
4305 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update Does the following things in the order:" << endl;
4306 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update * According to the switch: set the right ANN of the Array as generic one." << endl;
4307 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update * According to the switch: load the right weightfile as generic one." << endl;
4308 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update * According to the switch: set the right correction parameters...." << endl;
4309
4310 TString basicstring = TString(gSystem->ExpandPathName("$FEDRA_ROOT"));
4311 TString addstring = ("/src/libShower/weights/ID/");
4312
4313 // (*) SpecType 0 :
4315 addstring+="volumeSpec_CP/";
4317 }
4318 else {
4319 addstring+="volumeSpec_LT/";
4321 cout << "EdbShowerAlgIDSimple::Update WARNING::eSpecificationTypeString[0]=LT NOT YET SUPPORTED!" << endl;
4322 }
4323
4324 // (*) SpecType 1 :
4325 // SHID_vs_SPID
4327 addstring+="separationSpec_pion_vs_gamma/";
4329 }
4331 addstring+="separationSpec_electron_vs_pion/";
4333 }
4334 else {
4335 addstring+="separationSpec_gamma_vs_electron/";
4337 }
4338
4339 // (*) SpecType 2 :
4341 addstring+="efficiencySpec_MiddleVar/";
4343 }
4345 addstring+="efficiencySpec_LowEff/";
4347 }
4349 addstring+="efficiencySpec_MiddleFix/";
4351 }
4353 addstring+="efficiencySpec_All/";
4355 }
4356 else {
4357 addstring+="efficiencySpec_Neuch/";
4359 }
4360
4361 // (*) SpecType 3 :
4362 // Which addiditonal variables to use....
4363// 0: showershape only
4364// 1: and vtx IP
4365// 2: and vtx delta Z
4366// 3: shoowershape only, No first plate
4367// 4: shoowershape only, No first & second plate
4369 addstring+="showerInfoSpec_SHAPEWOPL12/";
4371 }
4373 addstring+="showerInfoSpec_SHAPEWOPL1/";
4375 }
4377 addstring+="showerInfoSpec_ANDVTXDELTAZ/";
4379 }
4381 addstring+="showerInfoSpec_ANDVTXIP/";
4383 }
4384 else {
4385 addstring+="showerInfoSpec_SHAPEONLY/";
4387 }
4388
4390
4392 CreateANN();
4393 }
4394
4395 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update Now load the different ANN weightstrings:" << endl;
4396 // This was the loading part for the ANN and
4397 // this is the loading part for the Correction Factors, which we take from the
4398 // generic file EnergyCorrections_Npl_%d.txt
4399 // For the ID class we do NOT need the correction factors, so we skip these files anyway....
4400 for (int ll=0; ll<15; ll++) {
4401 TString mkdirstring=basicstring+addstring;
4402 TString weigthstring=basicstring+addstring+TString(Form("weights_Npl_%d.txt",ANN_nPlates_ARRAY[ll]));
4404 ANN_MLP_ARRAY[ll]->LoadWeights(weigthstring);
4405 ANN_WeightFile_ARRAY[ll]=weigthstring;
4406 Float_t p0,p1;
4407 p0=0.0;
4408 p1=1.0;
4409 eANN_MLP_CORR_0[ll]=p0;
4410 eANN_MLP_CORR_1[ll]=p1;
4411 }
4412 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update * According to the switch: set the right ANN of the Array as generic one." << endl;
4413
4414 // (*) SpecType 4:
4415 // Get NplIndexNr:
4416 int check_Npl_index=0;
4418
4422
4423 // (*) SpecType 5:
4424 // HigEfficiency Cut, HighPurityCut, MinDistCut
4426 addstring+="cutType/";
4428 }
4430 addstring+="cutType/";
4432 }
4433 else {
4434 addstring+="cutType/";
4436 }
4437
4438 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple:: The path for the ID weightfile is. " << endl;
4440
4441
4442 // Set Generic ANN_Layout
4444
4445
4446 // Re-Read ID lookup tables, because of the possibility to
4447 // change the cutType threshold.
4448 ReadTables_ID();
4449
4450 if (gEDBDEBUGLEVEL >2) cout << "EdbShowerAlgIDSimple::Update WARNING:: Weightfiles for _LT_ not produced yet!...." << endl;
4451 cout << "EdbShowerAlgIDSimple::Update...done." << endl;
4452 return;
4453}
Member Data Documentation
◆ eANNIDCutValue
|
protected |
◆ eIDCutTypeValue
|
protected |
◆ eIDEnergySigmaCorr
|
protected |
The documentation for this class was generated from the following files:
- /home/antonio/fedra_doxygen/src/libShower/EdbShowerAlg.h
- /home/antonio/fedra_doxygen/src/libShower/EdbShowerAlg.cxx
