3246 Log(2,
"ShowRec.cpp",
"--- void ReconstructShowers_NN() ---");
3264 Float_t dT_InBT_To_TestBT=0;
3265 Float_t dR_InBT_To_TestBT=0;
3266 Float_t dR_TestBT_To_InBT=0;
3267 Float_t zDist_TestBT_To_InBT=0;
3268 Float_t SpatialDist_TestBT_To_InBT=0;
3269 Float_t zDiff_TestBT_To_InBT=0;
3270 Float_t dT_NextBT_To_TestBT=0;
3271 Float_t dR_NextBT_To_TestBT=0;
3272 Float_t mean_dT_2before=0;
3273 Float_t mean_dR_2before=0;
3274 Float_t mean_dT_before=0;
3275 Float_t mean_dR_before=0;
3276 Float_t mean_dT_same=0;
3277 Float_t mean_dR_same=0;
3278 Float_t mean_dT_after=0;
3279 Float_t mean_dR_after=0;
3280 Float_t mean_dT_2after=0;
3281 Float_t mean_dR_2after=0;
3283 Float_t min_dT_2before=0;
3284 Float_t min_dR_2before=0;
3285 Float_t min_dT_before=0;
3286 Float_t min_dR_before=0;
3287 Float_t min_dT_same=0;
3288 Float_t min_dR_same=0;
3289 Float_t min_dT_after=0;
3290 Float_t min_dR_after=0;
3291 Float_t min_dT_2after=0;
3292 Float_t min_dR_2after=0;
3293 Int_t nseg_1before=0;
3294 Int_t nseg_2before=0;
3295 Int_t nseg_3before=0;
3296 Int_t nseg_1after=0;
3297 Int_t nseg_2after=0;
3298 Int_t nseg_3after=0;
3303 TTree *
simu =
new TTree(
"MonteCarlo",
"Filtered Monte Carlo Events");
3304 simu->Branch(
"dT_InBT_To_TestBT", &dT_InBT_To_TestBT,
"dT_InBT_To_TestBT/F");
3305 simu->Branch(
"dR_InBT_To_TestBT", &dR_InBT_To_TestBT,
"dR_InBT_To_TestBT/F");
3306 simu->Branch(
"dR_TestBT_To_InBT", &dR_TestBT_To_InBT,
"dR_TestBT_To_InBT/F");
3307 simu->Branch(
"zDiff_TestBT_To_InBT", &zDiff_TestBT_To_InBT,
"zDiff_TestBT_To_InBT/F");
3308 simu->Branch(
"SpatialDist_TestBT_To_InBT", &SpatialDist_TestBT_To_InBT,
"SpatialDist_TestBT_To_InBT/F");
3309 simu->Branch(
"nseg_1before", &nseg_1before,
"nseg_1before/I");
3310 simu->Branch(
"nseg_2before", &nseg_2before,
"nseg_2before/I");
3311 simu->Branch(
"nseg_same", &nseg_same,
"nseg_same/I");
3312 simu->Branch(
"nseg_1after", &nseg_1after,
"nseg_1after/I");
3313 simu->Branch(
"nseg_2after", &nseg_2after,
"nseg_2after/I");
3315 simu->Branch(
"mean_dT_2before", &mean_dT_2before,
"mean_dT_2before/F");
3316 simu->Branch(
"mean_dT_before", &mean_dT_before,
"mean_dT_before/F");
3317 simu->Branch(
"mean_dT_same", &mean_dT_same,
"mean_dT_same/F");
3318 simu->Branch(
"mean_dT_after", &mean_dT_after,
"mean_dT_after/F");
3319 simu->Branch(
"mean_dT_2after", &mean_dT_2after,
"mean_dT_2after/F");
3321 simu->Branch(
"mean_dR_2before", &mean_dR_2before,
"mean_dR_2before/F");
3322 simu->Branch(
"mean_dR_before", &mean_dR_before,
"mean_dR_before/F");
3323 simu->Branch(
"mean_dR_same", &mean_dR_same,
"mean_dR_same/F");
3324 simu->Branch(
"mean_dR_after", &mean_dR_after,
"mean_dR_after/F");
3325 simu->Branch(
"mean_dR_2after", &mean_dR_2after,
"mean_dR_2after/F");
3327 simu->Branch(
"min_dT_2before", &min_dT_2before,
"min_dT_2before/F");
3328 simu->Branch(
"min_dT_before", &min_dT_before,
"min_dT_before/F");
3329 simu->Branch(
"min_dT_same", &min_dT_same,
"min_dT_same/F");
3330 simu->Branch(
"min_dT_after", &min_dT_after,
"min_dT_after/F");
3331 simu->Branch(
"min_dT_2after", &min_dT_2after,
"min_dT_2after/F");
3333 simu->Branch(
"min_dR_2before", &min_dR_2before,
"min_dR_2before/F");
3334 simu->Branch(
"min_dR_before", &min_dR_before,
"min_dR_before/F");
3335 simu->Branch(
"min_dR_same", &min_dR_same,
"min_dR_same/F");
3336 simu->Branch(
"min_dR_after", &min_dR_after,
"min_dR_after/F");
3337 simu->Branch(
"min_dR_2after", &min_dR_2after,
"min_dR_2after/F");
3339 simu->Branch(
"type", &
type,
"type/I");
3342 TTree *simu_SG = (TTree *)
simu->Clone();
3343 simu_SG->SetName(
"TreeSignalBT");
3344 simu_SG->SetTitle(
"TreeSignalBT");
3345 TTree *simu_BG = (TTree *)
simu->Clone();
3346 simu_BG->SetName(
"TreeBackgroundBT");
3347 simu_BG->SetTitle(
"TreeBackgroundBT");
3355 cout <<
"histos created"<<endl;
3361 cout <<
TMlpANN->GetStructure() << endl;
3368 Float_t local_gAli_pat_interim_halfsize=0;
3393 Int_t local_NBTMC=0;
3394 Int_t local_NBTallMC=0;
3395 Int_t local_NBTeMC=0;
3396 float_t local_pure=-1;
3397 float_t local_purall=-1;
3404 Int_t NBTMCe_Neff=0;
3408 cout << endl << endl <<
"--- Starting Shower for Number " << i <<
" now: "<<endl;
3425 dT_InBT_To_TestBT=0;
3426 dR_InBT_To_TestBT=0;
3427 dR_TestBT_To_InBT=0;
3428 zDiff_TestBT_To_InBT=0;
3461 Int_t local_gAli_npat=local_gAli->
Npatterns();
3462 if (
gEDBDEBUGLEVEL>2) cout <<
"--- local_gAli_npat= " << local_gAli_npat << endl;
3466 for (Int_t patterloop_cnt=local_gAli_npat-1; patterloop_cnt>=0; --patterloop_cnt) {
3467 if (
gEDBDEBUGLEVEL>3) cout <<
"--- --- Doing patterloop_cnt= " << patterloop_cnt << endl;
3473 for (Int_t btloop_cnt=0; btloop_cnt<local_gAli->
GetPattern(patterloop_cnt)->GetN(); ++btloop_cnt) {
3480 zDiff_TestBT_To_InBT=seg->
Z()-InBT->
Z();
3491 if (ann_inputneurons==10) {
3500 if (ann_inputneurons==20) {
3509 if (ann_inputneurons==30) {
3523 if (seg->
MCEvt()>0) {
3559 if (seg->
MCEvt()<0) {
3606 params[0]=dT_InBT_To_TestBT;
3607 params[1]=dR_InBT_To_TestBT;
3608 params[2]=dR_TestBT_To_InBT;
3609 params[3]=zDiff_TestBT_To_InBT;
3610 params[4]=SpatialDist_TestBT_To_InBT;
3619 params[10]=mean_dT_2before;
3620 params[11]=mean_dT_before;
3622 params[13]=mean_dT_after;
3623 params[14]=mean_dT_2after;
3625 params[15]=mean_dR_2before;
3626 params[16]=mean_dR_before;
3628 params[18]=mean_dR_after;
3629 params[19]=mean_dR_2after;
3631 params[20]=min_dT_2before;
3632 params[21]=min_dT_before;
3635 params[24]=min_dT_2after;
3637 params[25]=min_dR_2before;
3638 params[26]=min_dR_before;
3641 params[29]=min_dR_2after;
3660 if (seg->
X()==InBT->
X()&&seg->
Y()==InBT->
Y()) {
3672 CalcTrackDensity(pat_interim,local_gAli_pat_interim_halfsize,npat_int,npat_total,npatN);
3691 Double_t eff, purall, pure;
3742 cout <<
" Now write NN SG/BG Trees with the Inputvariables: --------------------"<<endl;
3744 TFile* fil =
new TFile(
"NN_ALG_ANN_TrainingsFile_DefaultName.root",
"RECREATE");
3746 canv_SG->Divide(5,4);
3907 TCanvas* canv_SG_2 =
new TCanvas(
"canv_SG_2",
"canv_SG_2",800,800);
3908 canv_SG_2->Divide(4,4);
void Get_NN_ALG_MLP(TTree *simu, Int_t parasetnr)
Definition: ShowRec.cpp:5990
void Load_NN_ALG_MLP_weights(TMultiLayerPerceptron *mlp, Int_t parasetnr)
Definition: ShowRec.cpp:6088
void Create_NN_Alg_Histograms()
Definition: ShowRec.cpp:8425
Double_t GetSpatialDist(EdbSegP *s1, EdbSegP *s2)
Definition: ShowRec.cpp:7276
Int_t GetMinsBeforeAndAfter(Float_t &min_dT, Float_t &min_dR, EdbPVRec *local_gAli, Int_t patterloop_cnt, EdbSegP *seg, Int_t n_patterns, Int_t BeforeOrAfter)
Definition: ShowRec.cpp:7295
Int_t GetMeansBeforeAndAfter(Float_t &mean_dT, Float_t &mean_dR, EdbPVRec *local_gAli, Int_t patterloop_cnt, EdbSegP *seg, Int_t n_patterns, Int_t BeforeOrAfter)
Definition: ShowRec.cpp:7383
Int_t GetNSegBeforeAndAfter(EdbPVRec *local_gAli, Int_t patterloop_cnt, EdbSegP *seg, Int_t n_patterns, Int_t BeforeOrAfter)
Definition: ShowRec.cpp:7462
TH1F * var_NN__BG__mean_dT_TestBT_ToSamePlate
Definition: ShowRec.h:275
TH1F * var_NN__SG__min_dR_TestBT_To2AfterPlate
Definition: ShowRec.h:326
TH1F * var_NN__BG__mean_dR_TestBT_To2BeforePlate
Definition: ShowRec.h:259
TH1F * var_NN__BG__nseg_TestBT_ToSamePlate
Definition: ShowRec.h:241
TH1F * var_NN__SG__min_dR_TestBT_ToSamePlate
Definition: ShowRec.h:310
TH1F * var_NN__BG__mean_dT_TestBT_To2BeforePlate
Definition: ShowRec.h:267
TH1F * var_NN__BG__mean_dT_TestBT_ToAfterPlate
Definition: ShowRec.h:283
TH1F * var_NN__SG__dT_InBT_To_TestBT
Definition: ShowRec.h:223
TH1F * var_NN__BG__mean_dR_TestBT_To2AfterPlate
Definition: ShowRec.h:287
TH1F * var_NN__SG__min_dT_TestBT_ToBeforePlate
Definition: ShowRec.h:302
TH1F * var_NN__BG__dR_InBT_To_TestBT
Definition: ShowRec.h:226
TH1F * var_NN__SG__mean_dR_TestBT_ToSamePlate
Definition: ShowRec.h:270
TMultiLayerPerceptron * TMlpANN
Definition: ShowRec.h:340
TH1F * var_NN__SG__min_dR_TestBT_ToAfterPlate
Definition: ShowRec.h:318
TH1F * var_NN__SG__mean_dR_TestBT_To2AfterPlate
Definition: ShowRec.h:286
TH1F * var_NN__BG__min_dR_TestBT_ToSamePlate
Definition: ShowRec.h:311
TH1F * var_NN__BG__nseg_TestBT_ToBeforePlate
Definition: ShowRec.h:234
TH1F * var_NN__SG__mean_dT_TestBT_To2AfterPlate
Definition: ShowRec.h:290
TH1F * var_NN__SG__mean_dR_TestBT_ToBeforePlate
Definition: ShowRec.h:254
TH1F * var_NN__SG__mean_dT_TestBT_ToBeforePlate
Definition: ShowRec.h:262
TH1F * var_NN__SG__mean_dT_TestBT_To2BeforePlate
Definition: ShowRec.h:266
TH1F * var_NN__SG__nseg_TestBT_To2AfterPlate
Definition: ShowRec.h:249
TH1F * var_NN__BG__mean_dR_TestBT_ToSamePlate
Definition: ShowRec.h:271
TH1F * var_NN__SG__min_dT_TestBT_ToSamePlate
Definition: ShowRec.h:314
TH1F * var_NN__SG__mean_dT_TestBT_ToAfterPlate
Definition: ShowRec.h:282
TH1F * var_NN__SG__SpatialDist_TestBT_To_InBT
Definition: ShowRec.h:214
TH1F * var_NN__SG__zDiff_TestBT_To_InBT
Definition: ShowRec.h:217
TH1F * var_NN__BG__min_dT_TestBT_To2BeforePlate
Definition: ShowRec.h:307
TH1F * var_NN__BG__mean_dR_TestBT_ToBeforePlate
Definition: ShowRec.h:255
TH1F * var_NN__BG__mean_dT_TestBT_ToBeforePlate
Definition: ShowRec.h:263
TH1F * var_NN__BG__dR_TestBT_To_InBT
Definition: ShowRec.h:229
TH1F * var_NN__BG__SpatialDist_TestBT_To_InBT
Definition: ShowRec.h:213
TH1F * var_NN__BG__nseg_TestBT_To2AfterPlate
Definition: ShowRec.h:250
TH1F * var_NN__BG__dT_InBT_To_TestBT
Definition: ShowRec.h:222
TH1F * var_NN__SG__dR_TestBT_To_InBT
Definition: ShowRec.h:228
TH1F * var_NN__BG__min_dR_TestBT_To2AfterPlate
Definition: ShowRec.h:327
TH1F * var_NN__BG__min_dR_TestBT_ToAfterPlate
Definition: ShowRec.h:319
TH1F * var_NN__SG__nseg_TestBT_ToBeforePlate
Definition: ShowRec.h:233
TH1F * var_NN__SG__mean_dR_TestBT_ToAfterPlate
Definition: ShowRec.h:278
TH1F * var_NN__BG__min_dT_TestBT_ToSamePlate
Definition: ShowRec.h:315
TH1F * var_NN__BG__zDiff_TestBT_To_InBT
Definition: ShowRec.h:216
Bool_t var_NN_DoTrain
Definition: ShowRec.h:211
TH1F * var_NN__SG__mean_dR_TestBT_To2BeforePlate
Definition: ShowRec.h:258
TH1F * var_NN__SG__min_dT_TestBT_To2BeforePlate
Definition: ShowRec.h:306
TH1F * var_NN__BG__min_dT_TestBT_ToAfterPlate
Definition: ShowRec.h:323
TH1F * var_NN__BG__mean_dR_TestBT_ToAfterPlate
Definition: ShowRec.h:279
TH1F * var_NN__BG__min_dR_TestBT_To2BeforePlate
Definition: ShowRec.h:299
TH1F * var_NN__SG__nseg_TestBT_ToAfterPlate
Definition: ShowRec.h:245
TH1F * var_NN__BG__mean_dT_TestBT_To2AfterPlate
Definition: ShowRec.h:291
TH1F * var_NN__BG__min_dR_TestBT_ToBeforePlate
Definition: ShowRec.h:295
TH1F * var_NN__SG__dR_InBT_To_TestBT
Definition: ShowRec.h:225
TH1F * var_NN__SG__nseg_TestBT_ToSamePlate
Definition: ShowRec.h:242
TH1F * var_NN__SG__mean_dT_TestBT_ToSamePlate
Definition: ShowRec.h:274
TH1F * var_NN__SG__min_dR_TestBT_To2BeforePlate
Definition: ShowRec.h:298
TH1F * var_NN__SG__nseg_TestBT_To2BeforePlate
Definition: ShowRec.h:237
TH1F * var_NN__BG__min_dT_TestBT_ToBeforePlate
Definition: ShowRec.h:303
TH1F * var_NN__BG__nseg_TestBT_ToAfterPlate
Definition: ShowRec.h:246
TH1F * var_NN__BG__nseg_TestBT_To2BeforePlate
Definition: ShowRec.h:238
TH1F * var_NN__SG__min_dR_TestBT_ToBeforePlate
Definition: ShowRec.h:294
TH1F * var_NN__SG__min_dT_TestBT_ToAfterPlate
Definition: ShowRec.h:322
TH1F * var_NN__BG__min_dT_TestBT_To2AfterPlate
Definition: ShowRec.h:331
TH1F * var_NN__SG__min_dT_TestBT_To2AfterPlate
Definition: ShowRec.h:330
Int_t type
Definition: testBGReduction_By_ANN.C:15
Double_t params[3]
Definition: testBGReduction_By_ANN.C:84